大慶長垣油田開發地震沉積學研究

李 操 姜 巖 樊曉東 劉宗堡 王元波 齊金成

(①東北石油大學地球科學學院，黑龍江大慶 163318；②大慶油田勘探開發研究院，黑龍江大慶163712)

1 開發地震沉積學概述

三維地震資料較鉆井資料的平面采樣密度大、橫向分辨率高，應用地震資料的橫向分辨率可以提高儲層預測能力。地震沉積學的形成和發展得益于20 世紀 70年代以來三維地震技術的快速發展和工業化應用。Dahm等[1]首次在地震時間切片上識別了曲流河道高分辨率振幅影像，體現了利用地震成像技術直接顯示古地形和古地貌的優點。隨著三維地震數據采集、處理技術的不斷進步，地震資料的面元更小、縱向分辨率更高，地震資料品質明顯提高，提高了對較小尺度地質體的解釋精度。在這種背景下，結合層序地層學基本原理和方法，Zeng等于1998年首次提出“地震沉積學”一詞，認為地震沉積學是綜合地震巖性學、地震地貌學研究地層巖性、沉積成因、沉積體系和盆地充填歷史的學科[2-5]。

2000 年之后，中國石油地質研究人員將地震沉積學理論方法引入中國陸相含油氣沉積盆地沉積體系研究、薄層砂體預測以及油氣精細勘探與開發等方面，并發揮了重要的作用[6-11]。朱筱敏等[12]在現代沉積學和地震沉積學理論指導下，確立松遼盆地齊家地區青山口組主要發育湖盆淺水三角洲沉積體系，利用地層切片刻畫了淺水三角洲沉積砂體的枝狀形態和大面積分布特征，指出近油源和鄰近斷層的三角洲前緣砂體為有利勘探、開發目標。20余年的地震沉積學理論發展和國內、外油氣精細勘探經驗表明，地震沉積學方法具有較高的橫向分辨率，可識別薄砂體(尺度小于1/4波長(厚度為3～5m)的砂體)的空間分布范圍，刻畫薄砂體的形態和分布特征，重建不同類型儲層分布和沉積體系的沉積過程。

利用地震沉積學研究油氣開發地質是近十年出現的全新領域，主要根據地層切片分析井間儲層分布，識別曲流河點壩、分流河道砂壩、廢棄河道等儲層[13-17]。吳勝和等[13]應用地震沉積學方法，通過三維地層切片等方式，利用地震成像特征研究沉積單元分布，研究儲層構型地震響應及解釋方法。張濤等[17]利用地層切片解釋曲流河沉積結構，建立了開發尺度的曲流河內部結構地震沉積學解釋技術。朱筱敏等[10]系統總結了中國地震沉積學研究現狀，提出“開發地震沉積學”一詞，但并未給出開發地震沉積學的具體定義。

大慶長垣油田位于松遼盆地中央坳陷區大慶長垣二級構造帶內，是一個大型背斜構造砂巖油氣田，研究區位于大慶長垣北部薩爾圖油田(圖1)，研究層位為嫩江組下部、姚家組和青山口組上部，發育薩爾圖、葡萄花、高臺子三套主力油層。2008年大慶長垣油田首次實現了覆蓋長垣北部薩爾圖油田的三維地震采集，獲得了高品質的高密度(面元尺寸為10m×10m)地震資料，為密井網區儲層描述提供了新的數據支撐。在密井網區預測厚度小于調諧厚度的砂泥巖薄互層儲層對地震技術提出了嚴峻挑戰，為此大慶油田組建了專業的開發地震團隊開展技術攻關，揭開了大慶長垣油田開發地震沉積學研究的序幕，并取得了較好應用效果。

圖1 大慶長垣油田區域位置圖

王元波等[18]研究了基于地質目標的保幅地震處理技術，根據密井網解剖地質研究成果，利用地震屬性相似性分析優化地震處理流程，定量質控與評價地震數據處理結果，提高了處理成果的保真度，在整體上獲得了與儲層沉積特征吻合、橫向變化信息豐富、窄條帶河流相砂體邊界清晰的地震資料處理成果。齊金成等[19]、姜巖等[20]、李操[21]在長垣油田建立了等時地層格架，在密井網開發區應用地層切片揭示了薩二組15+16(SⅡ15+16)等沉積單元的井間儲層變化，并取得了較好效果。郝蘭英等[22]針對不同砂體分布特點，利用地層切片采用“砂中找泥”和“泥中找砂”的單一河道識別方法，形成了“地震趨勢引導、井點微相控制”的井—震結合儲層精細刻畫方法。隨著地震沉積學研究和應用的不斷深入，出現地層切片優選效率低、井—震信息匹配分析難度大、地層切片與儲層厚度缺乏量化的對應關系等問題。為此，筆者進一步研究了地層切片自動優選、地震屬性平面可信度分析、井—震結合厚度定量預測等方法，結合以往研究成果形成了基于地震沉積學的儲層定量預測技術系列，并研發了井—震結合儲層預測軟件(iRPS)，提高了長垣油田開發地震沉積學儲層預測的精度和效率。

大慶長垣油田相關的研究與應用實踐表明，開發地震沉積學是地震沉積學在油田開發領域的深入研究與應用的產物，在開發區較密井網條件下針對厚度小于1/4波長的砂巖組、小層或沉積單元(單砂層)，注重于高效化、定量化研究井間儲層分布、連通關系、沉積特征等問題[18-24]。

2 開發地震沉積學方法

2.1 開發地震沉積學基礎研究

地震波形是振幅、頻率、相位等的綜合反映。從地質含義上說，波形橫向變化是地下地質體的巖性、物性、流體性質、沉積結構、構造等橫向變化的綜合反映，即地震波形是地下地質體相關地質和地球物理特征的綜合表現，其與地質體的波阻抗特征密切相關，只有了解了地層的波阻抗特征，才能建立有效的地震波形(或地震屬性)與儲層參數間的對應關系，降低儲層預測的不確定性。

2.1.1 砂泥巖波阻抗特征分析

大慶長垣油田薩爾圖、葡萄花、高臺子油層主要由不同速度、密度的細砂巖、粉砂巖、泥質粉砂巖、灰質粉砂巖、粉砂質泥巖和泥巖等組成，不同巖性對應不同波阻抗值。從取心井地層速度、聲波阻抗統計結果看，由于砂巖物性較好，整體上SⅡ及以下細砂巖速度最小，粉砂巖、泥質粉砂巖次之，泥巖、粉砂質泥巖速度較大，灰質粉砂巖速度最大，波阻抗具類似特征。在分析不同油組的聲波阻抗特征時，將巖性簡化為砂、泥巖兩元巖性介質。統計結果表明，SⅡ到高一組(GⅠ)砂巖具有低速度、低密度、低阻抗特征，泥巖具有高速度、高密度、高阻特征[25]，屬于“亮點型”儲層，砂、泥巖阻抗值平均相差約10%。依據波阻抗特征可較好地區分砂、泥巖，其中砂巖以地震波谷反射為主。

2.1.2 基于地質目標的保幅地震數據處理

在密井網區必須對地震資料進行相對保幅處理才能發揮地震技術的優勢。開發地震資料處理過程與常規的地震資料處理過程并無本質區別，但前者面向的地質目標是油田開發動用的儲層，因此對地震資料的保真度和分辨能力要求更高，質量控制過程更精細。

密井網條件的地質研究成果客觀反映了儲層砂體的整體分布特征，可以驗證地震儲層預測結果，進而檢驗地震處理成果的保真度，指導處理方法和參數優選。因此利用油田開發區井網密集、地質模式認知程度高的有利條件，針對地震資料處理過程中的關鍵環節，基于現有的地震資料處理的內部質控流程，采用井—震聯合質控方式，以地震反射特征與合成記錄吻合程度、地層切片的地震屬性整體趨勢與地質研究成果主要特征的相似程度為檢驗標準，進一步優化地震處理方法和參數。

油田區井網密集，可準確描述井點位置的儲層砂體地質特征，因此可以利用典型層位地質研究成果檢驗地震資料的砂體預測效果，進而有效評價地震處理成果保真度，優化地震處理參數。這種質量控制方法以地震數據處理、解釋一體化工作模式為支撐，通過地震資料處理試驗，結合地震解釋成果和效果評價，優化地震資料處理方法和參數。基于地質目標的保幅地震數據處理的關鍵在于處理參數的優選與質控，包括三個方面：一是質控技術系列。在常規地震處理各個環節，應用地震數據處理的內部質控技術縮小地震處理方法和參數的選擇范圍，初步優化方法和參數(圖2a)。二是合成記錄質控。對比、分析地震處理成果與合成記錄(圖2c)，評價地震處理參數效果。三是地質成果質控。對比地質研究成果(如砂巖等厚圖或沉積微相圖)(圖2b上)與地震資料處理某一環節的沿層振幅屬性(圖2b下)，通過相似程度分析、判斷處理成果的優劣。結合上述環節，可進一步判斷地震處理成果的合理性，在很大程度上縮小地震處理參數的選擇范圍，增強處理成果的可靠性。

圖2 基于地質目標的保幅地震處理質控示意圖

“六分法”去噪指分類、分域、分時、分頻、分區、分步的去噪方法；分域“兩步法”反褶積指將炮點域和檢波點域分開進行反褶積處理，可針對性地選擇不同域反褶積參數，從而提高反褶積效果

2.2 開發地震沉積學定性預測方法

2.2.1 等時地層格架建立

松遼盆地薩爾圖、葡萄花油層屬于大型陸相湖盆河流—三角洲沉積體系。上部坳陷超層序SSQ5由姚家組和嫩江組組成，底界面為SSB5(88Ma)，頂界面SSB6(72Ma)是嫩江組與四方臺組之間不整合面。湖平面于嫩一段沉積末期和嫩二段沉積初期再次上升到最高，其間沉積的油頁巖成為超層序的密集段，并形成松遼盆地的一套烴源巖層，這一時期沉積范圍達到最大，遠遠超出現今盆地范圍。SSQ5劃分為6個層序，從下到上依次為SQ19、SQ20、SQ21、SQ22、SQ23、SQ24，相當于姚家組至嫩江組[26](圖3)。

圖3 長垣地層格架對應關系圖

地層切片的核心思想是保證切片的地質等時性，地震同相軸等時地質意義取決于地震反射主體頻率。Zeng等[27]對美國二疊盆地三角洲前積體的正演模擬表明，低頻(15～30Hz)地震反射同相軸可能部分反映了厚層巖性界面信息，地震反射同相軸較易穿時，而高頻(50～80Hz)地震反射同相軸準確反映了較薄層的等時沉積界面信息，地震反射同相軸更趨于地質等時。長垣油田高密度地震資料在薩爾圖、葡萄花、高臺子油層的地震資料主頻約為45Hz，接近正演分析的等時條件，可識別以下層序界面[26]：SB19為地震反射次級標準層T1-1，為晚燕山運動主要幕次形成的區域性不整合面，對應上白堊統青山口組頂面，也對應GⅠ頂面，表現為中—低振幅、中—低連續反射特征；SB20為地震反射次級標準層T11，該界面為分隔姚家組一段和姚二、姚三段的反射界面，對應葡一油層組(PⅠ)頂面，在全區為較穩定、連續性較好的強反射；SB21為地震反射次級標準層T1，該界面為分隔姚家組和嫩江組的反射界面，對應SⅡ頂面，表現為橫向上分布非常穩定、連續性好的強反射。

大慶長垣油田地震標準層與地質標志層對應關系明確，地震標準層的同相軸等時性良好，因此長垣油田T1(SⅡ頂面)、T11(SⅢ底面即PⅠ頂面)、T1-1(GⅠ頂面)可以作為等時標志層，以此為基礎采用等比例線性內插方法建立高精度層序地層格架。根據采樣定理，地震數據采集時間確定了采樣間隔，制作地層切片相當于對地震數據重采樣，為避免制作地層切片重采樣時地震原有采樣點信息的丟失，對兩個標志層內插時，地層切片數不應小于兩標志層時間采樣個數。以SⅡ、薩Ⅲ油層組(SⅢ)為例，井震標定表明兩個油層組的地層總厚度約為90m，對應的兩個地震反射標準層T1與T11之間時間厚度約為60ms，考慮地層的橫向厚度變化，以T1、T11為等時界面應制作60張以上地層切片。

2.2.2 地層切片儲層預測

通過對比、分析地層切片發現，儲層預測效果不僅與目的層砂巖厚度有關，還與該層頂、底泥巖隔層的厚度有關，地層切片反映了儲層分布。如研究區X沉積單元的河道砂體厚度和砂體頂、底泥巖隔層厚度均較大，基于井數據得到的河道砂體邊界較清晰、可靠。應用該單元測井解釋砂巖厚度對地層切片進行標定(圖4)。可見，砂巖厚度大于3m的井多位于紅色區域，砂巖厚度為2～3m的井多落在紅色和白色區域，砂巖厚度小于2m的井多落在藍色區域，由地層切片刻畫的砂體邊界清晰，反映了砂體厚度的明顯變化。165口井數據的統計結果表明，除斷層部位井、震符合程度相對較低外，厚度大于2m的砂巖預測符合率達到75%。

圖4 長垣油田X沉積單元地層切片

2.3 開發地震沉積學定量預測方法

地層切片儲層預測方法在確定沉積單元砂體分布、儲層描述方面發揮了重要作用，但應用過程中也出現了一些問題：一是開發區井網密度大，由人工逐井制作合成記錄標定地震—地質層位的工作量大；二是長垣油田主力油層為砂泥巖薄互層，沉積單元數量多，需要人工優選沉積單元對應的地層切片，不同解釋者的優選結果有差異，同時工作效率較低；三是由于地震分辨率的限制，優選的地層切片除了包含目標沉積單元的有效信息外，還受上、下沉積單元的影響，需要利用測井數據由人工逐井判別地震信息的可靠性，井、震匹配難度大且工作量巨大；四是地層切片中的地震振幅數據雖然能夠反映砂體變化趨勢，但與砂體厚度缺乏量化的對應關系，且地層切片的色標調整會影響砂體的分布范圍。為此，筆者研究了開發地震沉積學定量預測方法。

2.3.1 地震—地質層位快速匹配

針對開發區井數據多、人工標定工作量大的問題，充分利用大慶長垣油田測井曲線齊全、油層組地質分層準確、區域構造解釋成果較可信等有利條件，研究了基于標志層的快速時深轉換方法。該方法的基本流程為：①以區域地震構造解釋標準層(如T1)和對應的地質分層(如SⅡ)建立基準時深對。②以聲波測井曲線計算測井數據段的地層速度，從基準時深對出發應用

(1)

計算全井初始的時深關系。式中：T(N)為第N個采樣點的雙程旅行時；Δd(1)為第N與第N-1個采樣點的深度差；VINT為第N和第N-1個采樣點間的層速度。③利用聲波和密度測井曲線計算波阻抗曲線和反射系數。④利用優選的理論子波或地震數據提取的子波與反射系數褶積得到合成記錄。⑤匹配、調整合成地震記錄與井旁地震道得到最優時深關系。依據以上思路編寫了地震—地質層位批量快速標定軟件模塊，集成于iRPS中。應用傳統的人工標定方法標定20km2的密井網開發區3000口井數據大約需要30個工作日，應用地震—地質層位批量快速標定軟件模塊僅需約1個工作日，大幅提高了井、震標定的工作效率，為密井網區大量深度域的井數據(如分層、測井曲線、巖性解釋結果)在時間域地震剖面快速投影提供了工具。經地震—地質層位快速匹配，各頂面地質分層位置與人工合成記錄標定位置一致，地震—地質層位快速匹配效果較好(圖5)。

圖5 地震—地質層位快速匹配效果剖面

2.3.2 地層切片自動優選方法

地震沉積學研究通常由有經驗的物探或地質解釋人員優選與目的層沉積規律吻合的地層切片，據此開展地震沉積學或地震地貌學研究。上述優選過程通常需要研究人員根據區域地質認識或井資料進行人工對比、分析，通常稱為“專家優選法”。大慶長垣油田儲層研究沉積單元眾多，應用專家優選法存在優選效率低、對專家經驗依賴強等問題。為此，筆者引入皮爾遜相關分析方法。定義兩個變量之間的皮爾遜相關系數為兩個變量之間的協方差和標準差的商，反映了兩變量之間的相關程度。應用皮爾遜相關分析公式

(2)

皮爾遜相關分析方法的應用基于以下兩個假設：一是多個連續的沉積單元在深度域按照深度序列是順序增加的，其地震響應在時間域也是按照反射時間序列順序增加的，利用深度域的位置比例關系可以確定時間域的位置比例關系，進而提供合理的地層切片優選時窗；二是在優選時窗內必能找到一張地層切片與沉積單元儲層信息(如砂巖厚度)達到最優匹配。基于以上假設，應用皮爾遜相關分析方法建立地層切片與砂巖厚度的相關系數曲線(圖6)自動優選地層切片。

圖6 地層切片相關系數優選曲線

2.3.3 地震屬性平面可信度分析方法

應用地層切片自動優選方法確定的地層切片在縱向上反映了目標沉積單元的地質特征，但在平面上不同位置的地震信息的有效性及是否能夠反映該單元的砂體平面展布規律仍需進一步判斷。為此，提出了地震屬性平面可信度分析方法，該方法同樣以皮爾遜相關分析理論為基礎，通過不同位置的地震屬性值與沉積單元儲層參數的相關性判斷地震屬性的可信度。具體的實施方法為：首先提取井點處地層切片地震屬性值和沉積單元儲層參數數據(通常為砂巖厚度數據)，構成平面可信度判別數據集；之后應用數據集分析平面不同位置井、震數據的相關性，進而得到地層切片地震屬性可信度；最后利用井點地層切片可信度繪制研究范圍內可信度平面圖。圖7為基于地層切片可信度的儲層定量預測結果。由圖可見：①基于井數據的沉積相圖(圖7a)反映了密井網條件的儲層認識。②優選的地層切片(圖7b)與圖7a的儲層特征一致。③地層切片可信度平面圖(圖7c)將可信度分為高(紅色)、中(綠色)及低(藍色)三級。④根據圖7c，針對不同的可信度(CL)區域采用不同方法刻畫儲層(圖7d，后文詳述)，即：將低可信度區域自動剔除，利用井數據和地質模式刻畫儲層；對中可信度區域利用井、震綜合信息刻畫儲層；對高可信度區域以地震信息為主，按照“井定位置、震控形態”的思路刻畫儲層。

圖7 基于地層切片可信度的儲層定量預測結果

2.3.4 井—震結合儲層厚度定量預測方法

根據Widess準則，當儲層厚度小于1/4地震波長時，地震振幅與儲層厚度呈正相關關系[28]。該準則為利用振幅數據以及各種地震屬性數據定量預測儲層厚度提供了理論基礎。大慶長垣油田中、淺層速度約為2800m/s，地震主頻約為40Hz，開發地質的研究對象為沉積單元，砂體厚度一般小于10m(小于1/4地震波長(17.5m))，滿足Widess準則的應用條件。地震正演和實際情況表明，受薄互層組合關系、砂巖厚度、隔層厚度等諸多因素的影響，地震屬性與儲層厚度的關系較復雜，利用簡單的線性回歸擬合方法難以得到地震振幅與儲層厚度的關系。為此，將地層切片可信度CL引入儲層定量預測，根據地層切片CL建立儲層厚度與地震振幅或地震屬性的函數關系，實現儲層厚度定量預測。該方法通過調整不同CL區域的井、震數據預測權重(在高CL區域地震數據參與的權重大，在低CL區域井數據參與的權重大)，建立基于地層切片可信度的井—震結合厚度預測方法

RT=CL×RTSEIS+(1-CL)×RTWELL

(3)

式中：RT為基于CL的儲層厚度預測值；RTSEIS為基于地震數據的儲層厚度預測值；RTWELL為基于井數據插值的儲層厚度預測值。

由多個地層切片或多個屬性數據應用多元回歸方法計算RTSEIS，基于井數據應用克里金插值方法得到RTWELL，應用式(3)得到RT(圖7d)。井—震結合儲層厚度預測結果可以得到每個地震面元的砂巖儲層厚度，實現井點吻合、井間地震約束。

2.4 開發地震沉積學的軟件實現

依據開發地震沉積學定性預測和定量預測的技術思路自主研發了井—震結合儲層預測軟件iRPS(圖8)。iRPS不僅具有常規地震、測井、地質數據管理功能及平面、剖面、三維等顯示功能，還具有單井分析功能，可完成常規地震、地質分析工作，且具有三項特色功能模塊：一是切片自動提取與優選，能夠以單層或雙層控制，批量提取振幅切片，并快速計算與儲層密切相關的時間切片、沿層切片和地層切片，提高切片優選的工作效率；二是基于可信度分析的儲層厚度定量預測，能夠自動建立多種地震屬性與儲層厚度之間的關系，通過結合井點儲層厚度與地震屬性數據預測儲層厚度；三是沉積微相預測，通過建立地震屬性、儲層厚度與沉積微相之間的關系，依據井—震結合預測沉積微相平面分布。

圖8 iRPS主要模塊

3 應用實例

應用開發地震沉積學方法和iRPS預測研究區儲層。研究區Y沉積單元屬于三角洲內前緣亞相，儲層橫向變化快，為利用地震橫向分辨率預測儲層提供了地質基礎。研究區的Y沉積單元地層厚度為0～12.3m，砂巖厚度為0～8.6m(小于1/4地震波長(約為17.5m))。對所有井數據的井—震聯合地震—地質層位快速匹配結果表明，Y沉積單元砂巖主要對應波谷反射。以T11、T1-1建立等時地層格架，利用地層切片自動優選和地震屬性平面可信度分析方法自動優選與優化地震屬性，應用基于可信度分析的井—震結合儲層厚度定量預測方法預測Y沉積單元的儲層厚度。后驗井證實，井—震結合砂巖厚度預測符合率達到85.6%。井—震預測結果發現了多條呈帶狀分布的砂體，形態與窄河道砂體一致。圖9為研究區0.5km2局部區域Y沉積單元砂巖厚度預測結果。由圖可見：①基于63口井數據插值得到的砂巖厚度變化趨勢(圖9a)顯示砂體邊界較平滑，東部和西部砂體較發育，中部砂體不發育，東北部有一條自東向西展布的條帶狀砂體，寬度約為150m。②井—震結合預測的砂巖厚度變化趨勢(圖9b)與圖9a基本一致，但存在很多細節變化，砂體邊界變化較快。如：工區東北部的條帶狀砂體寬度變窄，延伸長度變長；工區南部也出現條帶狀砂體。按照Y沉積單元的地質認識，儲層厚度為1.2m的位置可視為河道邊界，因此圖9a的河道邊界位于AA′的1井和2井中間，圖9b的河道邊界位于AA′靠近1井一側，即后者較前者向西移動了40m。③AA′地震剖面(圖9c)顯示，圖9b的河道邊界呈波谷反射，且振幅和相位發生明顯變化。后驗井證實，圖9b的河道邊界位置更準確，有利于快速追蹤河道邊界，分析井間砂體走向、連續性及連通性。④在圖9a中，BB′處的砂巖厚度為0～1m；在圖9b中，BB′處存在一個厚度大于1m的帶狀厚砂體，為井間發育的窄小河道，與東部河道相連；BB′地震剖面(圖9c)在3個位置出現振幅和相位異常，指示寬度為30m的河道。因此井—震結合儲層預測結果突破了井間均勻插值成圖的限制，充分體現了地震信息的約束作用。

圖9 研究區0.5km2局部區域Y沉積單元砂巖厚度預測結果

綜上所述，在密井網開發區，主河道帶內部仍然存在井數據無法確定的砂體變化，在井間可能發育寬度小于半個井間距的窄河道。以此為基礎，通過井—震結合解剖這些地質現象對于油田開發具有重要意義。

4 結論與認識

(1)開發地震沉積學是地震沉積學理論和方法在油田開發領域的深入研究與應用的產物，在開發區較密井網條件下針對厚度小于1/4波長的砂巖組、小層或沉積單元(單砂層)等地質目標，可研究井間儲層厚度變化、儲層邊界分布、儲層連通關系、沉積特征等問題。

(2)目的層波阻抗特征分析、基于地質目標的保幅地震數據處理是開發地震沉積學研究的基礎；開發地震沉積學定性預測方法是經典地震沉積學在開發領域的有效應用；開發地震沉積學定量預測方法注重于高效化、定量化儲層預測，其中地震—地質層位快速匹配、地層切片自動優選是針對開發區井密、層多問題的有效研究方法，地震屬性平面可信度分析、井—震結合儲層厚度定量預測是高精度預測開發區井間砂體的定量研究方法。這些方法初步構成了開發地震沉積學技術的基本框架，隨著研究的深入、地質目標的豐富，開發地震沉積學會逐步發展與完善。

(3)在密井網開發區仍然存在井數據無法確定的砂體變化，在井間可能發育寬度小于半個井間距的窄河道。通過開發地震沉積學有效預測與分析這些地質現象對于老油田開發具有重要意義。