[發明專利]一種定量計算薄儲層厚度的方法

有效

• 申請號: 201310032441.5
• 申請日: 2013-01-22
• 公開（公告）號: CN103412332A
• 公開（公告）日: 2013-11-27
• 發明（設計）人: 孫魯平
• 申請（專利權）人: 中國地質大學(北京)
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 暫無信息
• 代理人: 暫無信息

技術領域

本發明涉及應用地球物理中的地震勘探方法，特別是關于薄儲層勘探技術領域，具體來說是一種利用相對譜峰頻率增量定量計算薄儲層厚度的方法。

背景技術

應用地球物理中的地震勘探方法是指在地表以人工方法激發地震波，在向地下傳播時，遇有介質性質不同的巖層分界面，地震波將發生反射與透射，在地表用檢波器接收地震波，得到地震記錄。地震記錄上的特征與地下巖層的性質和結構有關，通過對地震記錄進行處理和解釋，可以推斷地下巖層的性質和形態。由于在勘查的精度和準確性方面，地震勘探均優于其它地球物理勘探方法，因而在油氣勘探中占有舉足輕重的地位。

但是，近幾十年隨著油氣勘探的不斷深入，勘探難度的增加，規模大的、易于尋找的構造圈閉油氣藏已經被發掘殆盡，我們的勘探目標不得不轉向薄儲層、巖性及復雜構造油氣藏。薄儲層是我國陸相含油氣盆地的重要儲集類型之一，是現階段我國石油勘探的主要研究目標。受制于地震采集、處理各個環節的影響，常規儲層預測手段分辨率不夠，達不到對薄儲層進行精細預測的效果，與實際生產的要求差距較大。因而，研究和發展新的薄(互)層地震定量預測方法具有重要的理論和實際意義。

對于薄儲層識別及厚度定量研究，國內外不少學者作過探討。1973年，Widess指出當薄層厚度小于λ/8時，薄層厚度與反射波形振幅的最大值呈線性關系，并給出了厚度計算公式，但他只討論了頂底反射系數大小相等，極性相反的薄層，且計算公式需要已知反射系數的絕對大小。Koefoed和Voogd(1981)、Robertson和Nogami(1984)等也都做過類似研究，得出了相似的結論。Chung和Lawton(1995)對Widess(1973)的近似公式進行了改進，提高了超薄層的計算精度。孫魯平等(2009，2010)對峰值頻率與薄層時間厚度關系進行了討論，提出了一種利用峰值頻率定量計算薄層時間厚度的思路，但并未給出獲取縱波速度將薄層時間厚度轉換為地層厚度的方法。Partyka和Gridley等(1999)提出譜分解方法，利用薄層頻譜的頻陷周期與薄層時間厚度存在倒數關系進行薄層厚度定量估算(Partyka(1999，2005)；Marfurt及Kirlin(2001)；黃緒德(2002)等)，但該方法受地震頻帶寬度的影響，地震頻寬不夠時頻譜的頻陷周期難以確定。在譜分解方法的基礎上，Castagna(2006)的研究團隊提出了譜反演技術，采用先驗信息和譜分解技術來提高小于調諧厚度的薄層成像效果(Chopra和Castagna等，2007；Puryear和Castagna等，2008)，該方法精度較高，但對子波準確性和資料的信噪比要求很高。此外，曾洪流(2005)提出90°相位轉換技術，更有利于對薄層的解釋；竇易升(1995)提出振幅譜平方比法定量解釋薄層厚度。黃真萍等(1997)，張恩珂等(2006)采用優化的非線性反演算法預測薄層厚度，這類方法需要有已知厚度的薄層樣本參與學習，只能用在有井的地區。綜上所述，目前的薄儲層定量預測方法或事先需要已知薄儲層頂、底的反射系數的絕對大小或對地震資料的采集處理有苛刻的要求，又或者嚴格依賴于鉆井，從而使這些方法在實際應用時受到很大的限制。

發明內容

本發明的目的是，提供了一種定量計算薄儲層厚度的方法。提出了一種新的地震屬性——相對譜峰頻率增量，通過該屬性計算薄儲層地層厚度，該方法不受制于薄儲層頂底面的反射系數絕對大小的限制，適用范圍廣且計算精度更高。

本發明的內容包括：

采集三維地震資料，對所述地震資料進行“三高”(高保真、高信噪比、高分辨率)處理獲得三維疊后地震數據；

利用測井和三維地震數據進行井-震合成記錄標定，確定時深關系，明確目標儲層所在位置，在包含目標儲層段內的地震數據上利用統計方法提取統計地震子波；

在目標層段內進行三維地震層位解釋，獲得時窗提取引導層，在時窗提取引導層的指導下，提取出包含目標儲層段的地震數據子體，對包含目標儲層段的時窗內的地震數據進行相對譜峰頻率增量地震屬性的提取；

利用相對譜峰頻率增量地震屬性定量計算目標薄儲層的雙程旅行時間厚度；

對測井資料進行統計分析，獲得目標層段內縱波阻抗與縱波速度二者的關系；

對三維地震數據進行疊后縱波阻抗反演，獲得三維縱波阻抗數據體，進而在時窗提取引導層的指導下，獲得包含目標儲層段的縱波阻抗數據子體；

提取出目標層段縱波阻抗平面地震屬性，利用目標層段縱波阻抗與縱波速度二者的關系，將目標層段縱波阻抗屬性轉換為縱波速度屬性，獲得目標層段的縱波速度平面屬性；