本申請要求2014年5月15日提交的第61/994,084號(hào)美國臨時(shí)專利申請的優(yōu)先權(quán)，所述美國臨時(shí)專利申請以全文引用的方式并入本文。

技術(shù)領(lǐng)域

本公開涉及地質(zhì)力學(xué)建模，且特定來說涉及使用來自時(shí)移地震數(shù)據(jù)的動(dòng)態(tài)邊界條件來更新或精細(xì)化地質(zhì)力學(xué)模型。

背景技術(shù)：

地下的地質(zhì)力學(xué)建模可用于研究由注入和生產(chǎn)引起的應(yīng)力和應(yīng)變。為了避免鉆井危險(xiǎn)性、最大化回收率以及確保儲(chǔ)層完整性，對于應(yīng)力和應(yīng)變的正確理解是必要的。

時(shí)移地震數(shù)據(jù)可以提供關(guān)于在兩次地震勘測之間地下的動(dòng)態(tài)行為的信息，包含地震事件的密度和速度變化以及位移?？梢允褂梅囱輥硎箷r(shí)移變化與巖石性質(zhì)、壓力、溫度、飽和度和巖石位移的變化相關(guān)聯(lián)。

地下的變化暗示了所述變化發(fā)生的位置中及周圍的改動(dòng)的應(yīng)力和應(yīng)變。在地質(zhì)力學(xué)模擬模型中，改動(dòng)的應(yīng)力和應(yīng)變狀態(tài)通常受到應(yīng)用于模型的應(yīng)力增量或位移增量的控制。應(yīng)力增量可以從壓力、溫度和飽和度的變化導(dǎo)出，而從時(shí)移地震數(shù)據(jù)導(dǎo)出的位移增量尚未得到廣泛研究。

從時(shí)移地震位移估計(jì)實(shí)際巖石位移依賴于巖石的速度的估計(jì)以及巖石位移如何改動(dòng)速度。相對速度變化與實(shí)際巖石應(yīng)變之間的關(guān)系的常用近似是如下定義的R因子(Hatchell和Bourne，2005a，2005b)：

此處δv/v是相對速度變化，而ε_zz是垂直應(yīng)變。假定變化是小的，則在雙向行進(jìn)時(shí)間上的相對變化(稱為時(shí)間應(yīng)變)可表達(dá)為：

從時(shí)移地震數(shù)據(jù)測得的雙向行進(jìn)時(shí)間偏移是累積時(shí)間應(yīng)變的作用。需要R因子的估計(jì)來將雙向行進(jìn)時(shí)間偏移轉(zhuǎn)換為實(shí)際巖石位移的估計(jì)。

地質(zhì)力學(xué)建模的重要目標(biāo)是使用模擬結(jié)果與時(shí)移觀測值之間的失配來更新材料性質(zhì)、斷層/斷裂模型，和/或巖石應(yīng)變-速度變化關(guān)系的模型。經(jīng)過適當(dāng)校準(zhǔn)的地質(zhì)力學(xué)模型可用于預(yù)測。另外，來自地質(zhì)力學(xué)建模的結(jié)果可以幫助解譯時(shí)移數(shù)據(jù)。

技術(shù)實(shí)現(xiàn)要素：

提供此發(fā)明內(nèi)容以介紹一系列概念，所述概念在以下詳細(xì)描述中進(jìn)一步描述。此發(fā)明內(nèi)容既不意在識(shí)別所要求主題的關(guān)鍵特征或本質(zhì)特征，也不意在用以輔助限制所要求主題的范圍。

根據(jù)本公開的一個(gè)實(shí)施方案，提供對地下的地質(zhì)力學(xué)效應(yīng)進(jìn)行建模以及針對時(shí)移觀測值調(diào)節(jié)地質(zhì)力學(xué)模型參數(shù)的新工作流程。通過從所觀測的時(shí)移行進(jìn)時(shí)間偏移和時(shí)間應(yīng)變導(dǎo)出的位移邊界條件來驅(qū)動(dòng)所述模型。從時(shí)移數(shù)據(jù)提取所述模型的邊界處的位移，從行進(jìn)時(shí)間偏移轉(zhuǎn)換為深度偏移和橫向偏移，且作為位移增量應(yīng)用于初始地質(zhì)力學(xué)模型。在一些方面中，由地質(zhì)力學(xué)模擬器計(jì)算應(yīng)力和應(yīng)變的增量，且可以將所述模型內(nèi)部的時(shí)移相關(guān)參數(shù)與所述時(shí)移觀測值進(jìn)行比較。本公開的實(shí)施方案實(shí)現(xiàn)對模擬與觀測之間的失配的全面研究，這可用于更新材料性質(zhì)、斷層、斷裂以及巖石應(yīng)變-速度變化關(guān)系(R因子)。根據(jù)本公開的一些實(shí)施方案，邊界條件方法使得包圍儲(chǔ)層的巖石的地質(zhì)力學(xué)模型可以從儲(chǔ)層模型的歷史匹配去耦。這使得根據(jù)這些實(shí)施方案的地質(zhì)力學(xué)模型的分析和歷史匹配更簡單且更準(zhǔn)確。

根據(jù)本發(fā)明的實(shí)施方案的一種方法包含：建立目標(biāo)體積的地質(zhì)力學(xué)模型；獲得所述目標(biāo)體積的一對時(shí)移地震數(shù)據(jù)體且計(jì)算偏移以使所述地震數(shù)據(jù)體中的事件對準(zhǔn)；估計(jì)邊界處的實(shí)際垂直位移；以從上一步驟導(dǎo)出的所述位移作為位移增量而運(yùn)行地質(zhì)力學(xué)模擬器；將模擬位移和應(yīng)變與時(shí)移觀測值進(jìn)行比較；更新R因子或模型；以及重復(fù)一個(gè)或多個(gè)步驟直到收斂。所述地質(zhì)力學(xué)模型可以在3維(3D)、2D或1D中。

附圖說明

參考附圖描述本公開的實(shí)施方案。在全部圖中使用相同標(biāo)號(hào)來提及相同的特征和組件。當(dāng)結(jié)合附圖考慮以下對若干實(shí)施方案的詳細(xì)說明時(shí)可得到對所述方法或設(shè)備的更好理解，附圖中：

圖1圖示了具有層位、斷層、斷裂和材料性質(zhì)的3D地質(zhì)力學(xué)模型。

圖2圖示了時(shí)間偏移：在左邊的地震圖像中由點(diǎn)線識(shí)別的事件在右邊的時(shí)移地震圖像中可見地偏移。使用非剛性匹配算法來計(jì)算3D偏移(底部)。

圖3圖示了從模擬結(jié)果確定因子α(位移的拉普拉斯算子對比垂直應(yīng)變)的圖表。

圖4圖示了從時(shí)間偏移估計(jì)的實(shí)際位移(左)對比從時(shí)間應(yīng)變估計(jì)的實(shí)際位移(右)的比較。

圖5圖示了位移增量作為具有屬性的點(diǎn)。

圖6圖示了在邊界條件更新之前(左)和在更新之后(右)的海床沉降失配。

圖7圖示了1D時(shí)間偏移比較：沿著軌跡/柵格列的模擬(藍(lán))對比所觀測時(shí)間偏移(黑)和平均化所觀測時(shí)間偏移(紅)。

圖8圖示了2D時(shí)間偏移比較：模擬時(shí)間偏移(左)對比所觀測時(shí)間偏移(右)。

圖9圖示了3D時(shí)間偏移比較：閾值模擬時(shí)間偏移(右)對比閾值所觀測時(shí)間偏移(右)。

圖10圖示了模擬時(shí)間應(yīng)變對比所觀測時(shí)間應(yīng)變的比較。

圖11圖示了針對一層(左)、相交視圖(右)的R因子估計(jì)的實(shí)施例。

圖12圖示了模擬應(yīng)變(左)和位移的拉普拉斯算子(右)。

圖13圖示了與模擬應(yīng)變共同可視化的井塌陷位置(藍(lán)點(diǎn))的圖表-紅色指示壓實(shí)，藍(lán)色指示擴(kuò)張，且白色指示其中可能發(fā)生強(qiáng)剪力的過渡區(qū)域。

圖14圖示了根據(jù)本發(fā)明的實(shí)施方案的用于地質(zhì)力學(xué)建模的工作流程圖。

圖15圖示了可以實(shí)現(xiàn)本申請的一個(gè)方法的計(jì)算機(jī)系統(tǒng)。

具體實(shí)施方式

現(xiàn)在將詳細(xì)參考實(shí)施方案，在附圖和圖式中說明實(shí)施方案的實(shí)施例。在以下詳細(xì)描述中，陳述許多具體細(xì)節(jié)以便提供對本文的主題的詳盡理解。然而本領(lǐng)域的技術(shù)人員將了解，可以在沒有這些具體細(xì)節(jié)的情況下實(shí)踐所述主題。在其它實(shí)例中，未詳細(xì)描述眾所周知的方法、程序、組件和系統(tǒng)，以免不必要地混淆實(shí)施方案的方面。

還將了解，雖然本文可使用術(shù)語“第一”、“第二”等等來描述各種元件，但這些元件不應(yīng)當(dāng)受這些術(shù)語限制。這些術(shù)語僅用以區(qū)分一個(gè)元件與另一個(gè)元件。舉例來說，第一對象或步驟可稱為第二對象或步驟，且類似地，第二對象或步驟可稱為第一對象或步驟。第一對象或步驟以及第二對象或步驟分別都是對象或步驟，但它們不應(yīng)視為相同對象或步驟。

在本文的公開的描述中使用的術(shù)語是僅用于描述特定實(shí)施方案的目的，且不希望限制主題。除非上下文另外清楚指示，否則如本說明書和所附權(quán)利要求書中所使用，單數(shù)形式“一”、“一個(gè)”和“所述”希望也包含復(fù)數(shù)形式。還將了解，如本文使用的術(shù)語“和/或”指代且涵蓋相關(guān)聯(lián)列出項(xiàng)目中的一個(gè)或多個(gè)的任何及所有可能的組合。進(jìn)一步將了解，術(shù)語“包含”和/或“包括”當(dāng)在本說明書中使用時(shí)指定所陳述特征、整數(shù)、步驟、操作、元件和/或組件的存在，但不排除一個(gè)或多個(gè)其它特征、整數(shù)、步驟、操作、元件、組件和/或其群組的存在或添加。

如本文使用，取決于上下文，術(shù)語“如果”可以解釋為意味著“當(dāng)...時(shí)”或“在...時(shí)”或“響應(yīng)于確定”或“響應(yīng)于檢測到”。類似地，取決于上下文，短語“如果確定”或“如果檢測到[所陳述條件或事件]”可以解釋為意味著“在確定...時(shí)”或“響應(yīng)于確定”或“在檢測到[所陳述條件或事件]時(shí)”或“響應(yīng)于檢測到[所陳述條件或事件]”。

圖14說明工作流程1400，其中：

(1)在步驟1(1410)中，建立目標(biāo)體積的地質(zhì)力學(xué)模型(參見圖1)。所述地質(zhì)力學(xué)模型可以包括1D、2D或3D模型。模型幾何形狀被配置為與來自地震數(shù)據(jù)的層位與斷層一致，適當(dāng)?shù)剞D(zhuǎn)換為深度。使用可用的信息(例如，鉆井日志、巖心測量、地震反演等等)將材料性質(zhì)指派于模型中的單元/元素(彈性性質(zhì)和控制塑性行為的性質(zhì))。使用地質(zhì)力學(xué)情形來確定模型中的初始應(yīng)力、位移和應(yīng)變；

(2)對于多個(gè)時(shí)移地震數(shù)據(jù)體，在步驟2(1420)中，計(jì)算使地震數(shù)據(jù)體中的事件對準(zhǔn)所必要的偏移(參見圖2)。舉例來說，在3D模型中，所述偏移是3D偏移。根據(jù)本公開的用于計(jì)算3D偏移的一種算法是非剛性匹配(參見Nickel等人的美國專利No.6,574,563)。對于3D情況，觀測到的3D偏移是與多次勘測之間的速度變化相關(guān)聯(lián)的3D位移場以及與實(shí)際巖石位移相關(guān)聯(lián)的3D位移場的總和。從作為3D偏移的垂直分量的行進(jìn)時(shí)間偏移，可以使用垂直導(dǎo)數(shù)來計(jì)算時(shí)間應(yīng)變；

(3)在步驟3(1430)中，估計(jì)在模型的邊界處的時(shí)移之間發(fā)生的實(shí)際位移。實(shí)際位移可以通過應(yīng)用以下關(guān)系中的一個(gè)或多個(gè)而確定：

其中：

δz是實(shí)際垂直偏移，δt是時(shí)間偏移，v是速度且R是R因子；

下標(biāo)K指示估計(jì)位移處的深度水平；

下標(biāo)k指示海床與目標(biāo)層級(jí)之間的層；

下表SB指示海床處的測量值；

橫杠表示海床與目標(biāo)層級(jí)之間的代表性值；

因子α是位移的差分斜率與垂直應(yīng)變之間的比率，其在簡單的板塊理論成立的情況下對于每一層是常數(shù)(因子α是依賴于材料性質(zhì)的，且可從地質(zhì)力學(xué)模擬結(jié)果導(dǎo)出，參見圖3)；

算子是拉普拉斯算子(橫向二階導(dǎo)數(shù)之和)，且時(shí)間偏移的時(shí)間導(dǎo)數(shù)是時(shí)間應(yīng)變；且

前四個(gè)關(guān)系提供了從時(shí)間偏移對實(shí)際位移的估計(jì)，而最后的關(guān)系給出了基于時(shí)間應(yīng)變的估計(jì)(參見圖4)。

應(yīng)注意，將在地質(zhì)力學(xué)模型的邊界處應(yīng)用的垂直(且也可能為水平)位移可以通過讀取邊界的位置處的估計(jì)(從時(shí)移地址數(shù)據(jù)體)偏移且針對“速度變化效應(yīng)”校正這些值(除以(1+R_average))而產(chǎn)生。其中此R_average可能必須通過例如鉆井日志和/或可能傾斜儀從沉降測量值來猜測、估計(jì)和/或外推。根據(jù)本公開的方面的方法接著繼續(xù)，且估計(jì)局部R因子(逐個(gè)單元、逐個(gè)巖性層)。這些局部R因子是依賴于材料的，即材料性質(zhì)。此外，它們也可以是依賴于狀態(tài)的，且是工作點(diǎn)(狀態(tài))周圍的非線性曲線的線性化的比例因子。

在估算局部R因子之后，可以計(jì)算且更新R平均因子。在正確的條件下，R平均值無變化，且過程收斂于其中R因子是最終結(jié)果的值。這些R因子可隨后用來預(yù)測地下的未來應(yīng)力狀態(tài)；

(4)任選地，在步驟(4)(1440)中，從步驟1中計(jì)算的位移場的水平分量估計(jì)實(shí)際的水平位移。替代地，從垂直位移與水平位移之間的某種關(guān)系可以估計(jì)水平位移；

(5)在步驟5(1450)中，運(yùn)行地質(zhì)力學(xué)模擬器，其中在步驟3和/或步驟4中導(dǎo)出的位移作為在步驟1中描述的初始狀態(tài)上的位移增量(參見圖5)；

(6)在步驟6(1460)中，將模擬的位移和應(yīng)變與時(shí)移觀測值進(jìn)行比較。更新材料模型，即斷層、斷裂和/或R因子模型，以減少失配。失配分析可以包含以下步驟中的一個(gè)或多個(gè)：

a.留下邊界的某個(gè)部分無位移約束。如果在邊界的此部分處存在位移的額外測量值，例如來自測深數(shù)據(jù)的表面沉降測量值或海床沉降測量值，那么可以分析失配。圖6圖示了海床沉降失配分析。

b.在1D、2D和3D中(參見圖7、圖8和圖9)將觀測到的時(shí)間偏移與模擬時(shí)間偏移進(jìn)行比較。對于時(shí)間應(yīng)變進(jìn)行相同的比較(參見圖10)。識(shí)別失配是否與材料性質(zhì)、斷層/斷裂或應(yīng)變-速度變化關(guān)系(R因子)相關(guān)。評估從時(shí)移數(shù)據(jù)測得的時(shí)間偏移是否充分平滑而與模擬結(jié)果進(jìn)行比較，或者是否平均化是必要的。此步驟還可以揭露處理地震數(shù)據(jù)中的問題，例如是否已經(jīng)引入非物理時(shí)間應(yīng)變。

c.使用以下關(guān)系計(jì)算使模擬應(yīng)變和觀測時(shí)間應(yīng)變相等所需的R因子。

圖11圖示了R因子計(jì)算的結(jié)果。

d.將時(shí)間偏移/位移的拉普拉斯算子與模擬應(yīng)變/時(shí)間應(yīng)變進(jìn)行比較(參見圖12)。如果簡單板塊行為中的假設(shè)是有效的，那么拉普拉斯算子和應(yīng)變應(yīng)當(dāng)是成比例的。識(shí)別其中假設(shè)似乎有效的區(qū)以及其中假設(shè)似乎無效的區(qū)。

e.針對所報(bào)告的井問題比較壓實(shí)與擴(kuò)張之間的過渡區(qū)(參見圖13)。評估剪力是否會(huì)造成井問題，以及是否應(yīng)當(dāng)更新模型以匹配井問題的區(qū)；以及

(7)在步驟7(1470)中，如果分析導(dǎo)致R因子的修改，那么方法回到步驟3且使用被更新的R因子來計(jì)算新位移邊界條件。如果材料性質(zhì)、斷層或斷裂被更新，那么重新運(yùn)行模型且方法回到步驟6以進(jìn)行進(jìn)一步分析。

應(yīng)注意，可以不同次序執(zhí)行步驟，且個(gè)別步驟可為任選的。在前五個(gè)步驟期間，可以使用地震時(shí)移數(shù)據(jù)來找到用以限定地質(zhì)力學(xué)模型的邊界的時(shí)間偏移。步驟6(1460)可用于將來自內(nèi)部體積(不是邊界)中的模型的模擬性質(zhì)與來自同一體積的時(shí)移觀測的那些性質(zhì)進(jìn)行比較。在一些方面中，可以分析失配以找到它們的原因。

一旦從以上工作流程得到精細(xì)模型，便可使用所述模型用于許多不同應(yīng)用，包含例如：

通過將模擬運(yùn)行的結(jié)果與觀測到的時(shí)間偏移進(jìn)行比較且引入從模擬考慮預(yù)期趨勢的對所觀測時(shí)間偏移的自適應(yīng)平滑/平均化，而使得時(shí)間偏移估計(jì)更物理上合理；

識(shí)別處理假象，其中如果在方法中找到的失配是物理上合理的，那么可將它們識(shí)別為處理假象且可移除，這可改善處理方法和所得數(shù)據(jù)質(zhì)量；

通過從時(shí)間應(yīng)變與模擬應(yīng)變的比較局部地估計(jì)R因子來估計(jì)應(yīng)變與速度變化(R因子)之間的關(guān)系。這是針對多個(gè)時(shí)移而完成，使得當(dāng)找到一致的R因子時(shí)，可以將時(shí)間應(yīng)變轉(zhuǎn)換為實(shí)際應(yīng)變，這允許方法用于鉆井危險(xiǎn)性、井完整性和儲(chǔ)層完整性的定量分析。R因子的現(xiàn)實(shí)且可信估計(jì)是來自以上工作流程的極為有用的結(jié)果；以及

使用來自方法的R因子使用它們通過三階彈性常數(shù)與材料參數(shù)(密度和速度)的關(guān)系來更新材料性質(zhì)/三階彈性參數(shù)。在此分析中，可以使用R的估計(jì)來估計(jì)其它材料性質(zhì)，例如與密度和速度數(shù)據(jù)一致的三階彈性常數(shù)。如果所得估計(jì)是不現(xiàn)實(shí)/不合理的，那么更新材料性質(zhì)且重新運(yùn)行模型。

圖1至13圖示了使用上文論述的方法的若干實(shí)施例。舉例來說：

圖1圖示了具有層位、斷層、斷裂和材料性質(zhì)的3D地質(zhì)力學(xué)模型；

圖2圖示了時(shí)間偏移，其中在地震圖像中由點(diǎn)線識(shí)別的事件在時(shí)移地震圖像中可見地偏移。使用根據(jù)本申請的非剛性匹配算法來計(jì)算底部圖示中的3D偏移。

圖3圖示了從模擬結(jié)果確定因子α(位移的拉普拉斯算子對比垂直應(yīng)變)；

圖4圖示了從時(shí)間偏移估計(jì)的實(shí)際位移(左圖像)對比從時(shí)間應(yīng)變估計(jì)的實(shí)際位移(右圖像)；

圖5圖示了位移增量作為具有屬性的點(diǎn)。在此實(shí)施例中，僅非零位移在模型的基底處(頂部儲(chǔ)層)。在側(cè)面處的位移增量設(shè)定為零，而在頂部(海床)處的位移不受約束，使得可以對照測深數(shù)據(jù)比較模擬位移；

圖6圖示了在邊界條件更新之前(左圖像)和在更新之后(右圖像)的海床沉降失配；

圖7圖示了1D時(shí)間偏移比較：沿著軌跡/柵格列的模擬(藍(lán))710對比所觀測時(shí)間偏移(黑)720和平均化所觀測時(shí)間偏移(紅)730；

圖8圖示了2D時(shí)間偏移比較：模擬時(shí)間偏移(左圖像)對比所觀測時(shí)間偏移(右圖像)；

圖9圖示了3D時(shí)間偏移比較：閾值模擬時(shí)間偏移(右圖像)對比閾值所觀測時(shí)間偏移(左圖像)；

圖10圖示了模擬時(shí)間應(yīng)變對比所觀測時(shí)間應(yīng)變；

圖11圖示了針對一層(左圖像)和相交視圖(右圖像)的R因子估計(jì)的實(shí)施例；

圖12圖示了模擬應(yīng)變(左圖像)和位移的拉普拉斯算子(右圖像)；以及

圖13圖示了與模擬應(yīng)變共同可視化的井塌陷位置(藍(lán)點(diǎn))1310-紅色區(qū)1320指示壓實(shí)，藍(lán)色區(qū)1330指示擴(kuò)張，且白色區(qū)1340指示其中可能發(fā)生強(qiáng)剪力的過渡區(qū)域。

如本領(lǐng)域的技術(shù)人員將了解，可以組合上文論述的方法的步驟中的一個(gè)或多個(gè)，步驟可以是任選的，和/或可以改變一些操作的次序。此外，方法中的一些操作可以與本文公開的其它實(shí)施例實(shí)施方案的方面組合，和/或可以改變一些操作的次序。測量過程、其解釋以及操作者采取的動(dòng)作可以重復(fù)方式完成；此概念適用于本文論述的方法。

圖15圖示了用于執(zhí)行本文描述的方法的計(jì)算系統(tǒng)1500。系統(tǒng)計(jì)算機(jī)1530可以與可為外部硬盤存儲(chǔ)裝置的磁盤存儲(chǔ)裝置1529、1531、1533和1535通信。預(yù)期磁盤存儲(chǔ)裝置1529、1531、1533和1535是常規(guī)硬盤驅(qū)動(dòng)器，且因此將借助于局域網(wǎng)或通過遠(yuǎn)程接入來實(shí)施。當(dāng)然，雖然將磁盤存儲(chǔ)裝置說明為單獨(dú)的裝置，但可以按需要使用單個(gè)磁盤存儲(chǔ)裝置來存儲(chǔ)任何及所有程序指令、測量數(shù)據(jù)和結(jié)果。

在一種實(shí)現(xiàn)中，數(shù)據(jù)可以存儲(chǔ)在磁盤存儲(chǔ)裝置1531中。來自不同源的各種數(shù)據(jù)可以存儲(chǔ)在磁盤存儲(chǔ)裝置1533中。系統(tǒng)計(jì)算機(jī)1530可以從磁盤存儲(chǔ)裝置1531或1533檢索適當(dāng)數(shù)據(jù)以根據(jù)對應(yīng)于本文描述的各種技術(shù)的實(shí)現(xiàn)的程序指令來處理數(shù)據(jù)。程序指令可以用例如C++、Java和類似語言等計(jì)算機(jī)編程語言編寫。程序指令可以存儲(chǔ)在例如程序磁盤存儲(chǔ)裝置1535等計(jì)算機(jī)可讀媒介中。此計(jì)算機(jī)可讀媒介可以包含計(jì)算機(jī)存儲(chǔ)媒介。計(jì)算機(jī)存儲(chǔ)媒介可以包含在用于信息存儲(chǔ)的任何方法或技術(shù)中實(shí)現(xiàn)的易失性和非易失性以及可裝卸式和不可裝卸式媒介，例如計(jì)算機(jī)可讀指令、數(shù)據(jù)結(jié)構(gòu)、程序模塊或其它數(shù)據(jù)。計(jì)算機(jī)存儲(chǔ)媒介可以進(jìn)一步包含RAM、ROM、可擦除可編程只讀存儲(chǔ)器(EPROM)、電可擦除可編程只讀存儲(chǔ)器(EEPROM)、快閃存儲(chǔ)器或其它固態(tài)存儲(chǔ)器技術(shù)、CD-ROM、數(shù)字多功能光盤(DVD)、或其它光學(xué)存儲(chǔ)裝置、磁性卡帶、磁帶、磁盤存儲(chǔ)裝置或其它磁性存儲(chǔ)裝置，或者可用于存儲(chǔ)所需信息且可由系統(tǒng)計(jì)算機(jī)1530存取的任何其它媒介。以上各項(xiàng)中的任一者的組合也可以包含在計(jì)算機(jī)可讀媒介的范圍內(nèi)。

在一種實(shí)現(xiàn)中，系統(tǒng)計(jì)算機(jī)1530可以將輸出主要呈現(xiàn)到圖形顯示器1527上，或替代地經(jīng)由打印機(jī)1528呈現(xiàn)(未圖示)。系統(tǒng)計(jì)算機(jī)1530可以將上文描述的方法的結(jié)果存儲(chǔ)在磁盤存儲(chǔ)裝置1529上以用于稍后使用和進(jìn)一步分析。鍵盤1526和指向裝置(例如，鼠標(biāo)、軌跡球或類似物)1525可以隨系統(tǒng)計(jì)算機(jī)1530提供以實(shí)現(xiàn)交互式操作。

系統(tǒng)計(jì)算機(jī)1530可以位于遠(yuǎn)離開采場的數(shù)據(jù)中心。系統(tǒng)計(jì)算機(jī)1530可以與現(xiàn)場的設(shè)備通信以接收各種測量的數(shù)據(jù)。系統(tǒng)計(jì)算機(jī)1530還可以位于場地的現(xiàn)場以提供較快的反饋和現(xiàn)場操作的指導(dǎo)。此數(shù)據(jù)在常規(guī)格式化和其它初始處理之后可以由系統(tǒng)計(jì)算機(jī)1530存儲(chǔ)為磁盤存儲(chǔ)裝置1531或1533中的數(shù)字?jǐn)?shù)據(jù)，用于以上文描述的方式的后續(xù)檢索和處理。雖然圖15將磁盤存儲(chǔ)裝置(例如，1531)圖示為直接連接到系統(tǒng)計(jì)算機(jī)1530，但也預(yù)期磁盤存儲(chǔ)裝置可以通過局域網(wǎng)或遠(yuǎn)程存取而存取。此外，雖然將磁盤存儲(chǔ)裝置1529、1531圖示為用于存儲(chǔ)輸入數(shù)據(jù)和分析結(jié)果的單獨(dú)裝置，但磁盤存儲(chǔ)裝置1529、1531可以在單個(gè)磁盤驅(qū)動(dòng)器(與程序磁盤存儲(chǔ)裝置1533一起或分開)內(nèi)實(shí)現(xiàn)，或者以參考本說明書的本領(lǐng)域的技術(shù)人員將完全理解的任何其它常規(guī)方式來實(shí)現(xiàn)。

以上公開的特定實(shí)施方案僅是說明性的，因?yàn)榭梢杂玫靡嬗诒疚慕淌镜谋绢I(lǐng)域的技術(shù)人員了解的不同但等效方式來修改和實(shí)踐本發(fā)明。此外，除了所附權(quán)利要求書中所描述的之外，不希望對本文所示的構(gòu)造或設(shè)計(jì)的細(xì)節(jié)存在限制。因此顯然可以更改或修改上文公開的特定實(shí)施方案，且所有這些變化視為在本發(fā)明的范圍內(nèi)。因此，本文尋求的保護(hù)范圍如所附權(quán)利要求書中陳述。

雖然上文僅詳細(xì)描述了幾個(gè)實(shí)施例實(shí)施方案，但本領(lǐng)域的技術(shù)人員將容易了解，在不實(shí)質(zhì)上脫離本發(fā)明的情況下在實(shí)施例實(shí)施方案中許多修改是可能的。因此，希望所有此類修改包含于如所附權(quán)利要求書中界定的本公開的范圍內(nèi)。在權(quán)利要求書中，裝置加功能子句希望涵蓋本文描述為執(zhí)行所陳述功能的結(jié)構(gòu)，且不僅涵蓋結(jié)構(gòu)等效物也包含等效結(jié)構(gòu)。因此，雖然因?yàn)獒斪硬捎脠A柱形表面來將木制零件固定在一起，而螺釘采用螺旋形表面，所以釘子和螺釘不會(huì)是結(jié)構(gòu)等效物，但是在將木制零件緊固的環(huán)境中，釘子和螺釘可以是等效結(jié)構(gòu)。申請人的表達(dá)意圖不是針對本文的任何權(quán)利要求的任何限制調(diào)用35U.S.C.§112第6段，在權(quán)利要求明確使用詞語‘用于...的裝置’以及相關(guān)聯(lián)功能的那些情況除外。

參考文獻(xiàn)：

Hatchell，P.J.，and Bourne，S.J.，2005a.Measuring reservoir compaction using time-lapse time shifts.75^th Annual International MeetingSociety of Exploration Geophysicists，2500-2503.

Hatchell，P.J.，and Bourne，S.J.，2005b.Rocks under strain：Strain-induced time-lapse time shifts are observed for depleting reservoirs.The Leading Edge，24(12)，1244-1248.

Nickel，M.，and Sonneland，L.1999.Non-rigid matching of migrated time-lapse seismic.European Association of Geoscientists and Engineers，Extended Abstracts.

US Patent No.6,574,563