海域油氣蓋層疊前地震反演預(yù)測方法及應(yīng)用

摘要：

蓋層對有效儲層刻畫、儲層質(zhì)量評估和油氣藏開發(fā)具有重要意義。現(xiàn)有的海域地震預(yù)測方法主要針對儲層，缺乏對蓋層空間展布的研究與分析，亟需建立海域油氣蓋層疊前地震反演預(yù)測方法，為海域油氣蓋層的地震預(yù)測提供技術(shù)支撐。因此，本文提出一種基于巖石物理約束的海域蓋層巖相以及物性參數(shù)疊前地震預(yù)測方法，對海域蓋層巖相以及物性參數(shù)的地震描述進行深入研究。首先，通過海域地層的巖石物理建模分析，優(yōu)選海域油氣蓋層巖相以及物性敏感彈性參數(shù)。其次，基于巖石物理模型，分別推導(dǎo)巖相敏感參數(shù)以及物性參數(shù)直接表征的地震反射系數(shù)方程和彈性阻抗方程，為海域蓋層巖相的空間展布地震刻畫和物性參數(shù)的地震預(yù)測提供技術(shù)支撐。為了解決模型參數(shù)之間的相關(guān)性導(dǎo)致反演結(jié)果不穩(wěn)定的問題，在地震預(yù)測過程中通過奇異值分解方法實現(xiàn)模型參數(shù)的去相關(guān)性，并在貝葉斯反演理論的框架下分別推導(dǎo)海域蓋層敏感彈性參數(shù)、物性參數(shù)的反演目標函數(shù)。將該方法應(yīng)用于實際資料，結(jié)果表明泊松阻抗預(yù)測精度可以達到95%，孔隙度預(yù)測精度可以達到85%，驗證了海域油氣蓋層疊前地震反演預(yù)測方法的合理性與有效性。

關(guān)鍵詞：海域油氣；蓋層；孔隙度；泊松阻抗；疊前反演；巖石物理模型

doi：10.13278/j.cnki.jjuese.20240255

中圖分類號：P631.4

文獻標志碼：A

鄭軒，宗兆云，付亞群，等. 海域油氣蓋層疊前地震反演預(yù)測方法及應(yīng)用. 吉林大學(xué)學(xué)報（地球科學(xué)版），2024，54（6）：21282141. doi：10.13278/j.cnki.jjuese.20240255.

Zheng Xuan， Zong Zhaoyun， Fu Yaqun， et al. Pre-Stack Seismic Inversion Prediction Method and Application for Marine Oil and Gas Cap Rocks. Journal of Jilin University （Earth Science Edition）， 2024， 54 （6）： 21282141. doi：10.13278/j.cnki.jjuese.20240255.

收稿日期：20240924

作者簡介：鄭軒（1998—），男，博士研究生，主要從事油氣地球物理理論方法與應(yīng)用等方面的研究，E-mail： b23010048@s.upc.edu.cn

基金項目：國家自然科學(xué)基金項目（42174139，41974119，42030103）；嶗山實驗室科技創(chuàng)新計劃（LSKJ202203406）

Supported by the National Natural Science Foundation of China （42174139， 41974119， 42030103） and the Laoshan Laboratory Science and Technology Innovation Program （LSKJ202203406）

Pre-Stack Seismic Inversion Prediction Method and Application for Marine Oil and Gas Cap Rocks

Zheng Xuan1， Zong Zhaoyun1， Fu Yaqun2， Luo Kun1

1. School of Geosciences， China University of Petroleum （East China）， Qingdao 266580， Shandong， China

2. Sinopec Exploration Company， Chengdu 610041， China

Abstract：

Cap rocks are

of great significance for effective reservoir characterization， reservoir quality assessment， and oil and gas reservoir development. The existing marine seismic prediction methods mainly focus on reservoirs， lacking research and analysis on the spatial distribution of cap rocks. It is urgent to establish a pre-stack seismic inversion prediction method for marine oil and gas cap rocks to provide technical support for the seismic prediction of marine oil and gas cap rocks. Therefore， this paper proposes a pre-stack seismic prediction method for marine cap rock facies and physical property parameters based on rock physics constraints， and conducts an in-depth study on the seismic description of" marine cap rock facies and physical property parameters. First， through the rock physics analysis of marine strata， the marine oil and gas cap rock facies and physical property sensitive elastic parameters are optimized. Secondly， based on the rock physics model， the seismic reflection coefficient equation and elastic impedance equation directly characterized by the lithofacies sensitive parameters and physical property parameters are derived respectively， providing technical support for the seismic characterization of the spatial distribution of the marine cap rock facies and the seismic prediction of the physical property parameters. In order to solve the problem of unstable inversion results caused by the correlation between model parameters， the singular value decomposition method is used to realize the decorrelation of model parameters in the seismic prediction process， and the inversion objective functions of sensitive elastic parameters and physical parameters of"" marine cap rocks" are derived respectively under the framework of Bayesian inversion theory. This method is applied to actual data to verify the rationality and effectiveness of the pre-stack seismic inversion prediction method for marine oil and gas cap rocks. The results show that the Poisson impedance prediction accuracy can reach 95%， and the porosity prediction accuracy can reach 85%.

Key words：

marine oil and gas; caprocks; porosity; Poisson impedance; pre-stack inversion; rock physics model

0" 引言

蓋層作為海域有利油氣儲層形成的重要環(huán)節(jié)，對海域油氣運移、聚集以及成藏都有顯著的影響［12］。針對蓋層的地震預(yù)測方法研究，能夠為海域有效儲層的地震勘探開發(fā)提供重要的技術(shù)基礎(chǔ)。蓋層的非滲透性特征對易揮發(fā)的油氣具有重要的保護以及封蓋作用，但目前對海域蓋層的地震識別與預(yù)測方面涉及較少［3］，亟待開展海域蓋層的地震反演方法研究，為海域有利儲層地震勘探提供技術(shù)支撐。

蓋層在海域油氣勘探開發(fā)中的研究地位越發(fā)突出，同時也面臨著一些亟待解決的難題。我國海域大型油氣田成藏模式復(fù)雜［4］。近年來，湖相、海相、海陸過渡相儲蓋組合方式以及有利儲蓋控富集、構(gòu)造脊蓋層控運等多種成藏模式被提出并成功指導(dǎo)了海域油氣的勘探與開發(fā)［5］。在珠江口盆地珠三凹陷，由于太平洋板塊俯沖減弱、海底擴張逐漸停止以及構(gòu)造活動較弱等因素，沉積了穩(wěn)定的淺海泥巖蓋層。近期在珠江口盆地的新近系儲層中，發(fā)現(xiàn)了可能為沉積成因的鈣質(zhì)砂巖夾層［6］，使得泥巖和鈣質(zhì)砂巖成為該海域蓋層的主要成分。在渤海灣盆地，大氣田的勘探開發(fā)一直是研究的重點目標，其區(qū)域蓋層為古近系下部厚層泥巖，其成巖程度高、封堵能力強的特征，有助于渤海灣盆地大氣田的形成［7］。南海北部大陸邊緣瓊東南盆地有利油氣聚集帶也擁有較好的儲蓋組合條件［8］，近期研究表明，斷層封堵失效、蓋層欠發(fā)育等原因成為了探井過程中的重要阻礙［9］。因此，形成有效的海域蓋層疊前地震反演預(yù)測方法，對我國海域油氣增儲上產(chǎn)具有重要意義。

20世紀70年代以前，人們很少關(guān)注蓋層。80年代以來，隨著復(fù)雜油氣藏勘探開發(fā)需求的增多，蓋層的研究成果日益豐富。1984年，Downey［10］首先指出復(fù)雜油氣勘探開發(fā)過程不應(yīng)忽視蓋層，應(yīng)從微觀和宏觀兩個角度評估油氣蓋層，以降低勘探風(fēng)險。1981年，Grunau［11］對蓋層進行了廣泛的地質(zhì)調(diào)查，指出擴散速率和裂縫的形成與演化是蓋層評價的重要參數(shù)。1993年，張樹林等［12］系統(tǒng)地介紹了蓋層的封閉能力和評價蓋層的方法，確定了包括蓋層厚度、平面展布、塑性、孔隙流體壓力、排驅(qū)壓力、滲透率等12個主要蓋層評價參數(shù)，并給出了蓋層的分類方法。而后，付廣等［1314］豐富了蓋層封閉性理論，并利用鉆井、測井和分析測試資料對蓋層類型、識別標志、分布及封閉能力進行了系統(tǒng)研究。Wu等［15］基于密封強度指數(shù)，實現(xiàn)了對井上油氣蓋層封閉性的定量評估。以上研究受限于一維油井提供的蓋層封閉能力等信息，無法準確估計蓋層平面展布。近年來，利用地震資料評估蓋層密封性逐漸成為蓋層研究的重點。2009年，Davis等［16］利用時移多分量地震資料探測了科羅拉多州魯利森氣田盆地中心氣藏內(nèi)的封閉性。2010年，Nourollah等［17］運用地震疊后屬性指示了封閉性。后來，Nourollah等［18］開發(fā)了一種利用地震各向異性研究頁巖層序封閉潛力的方法。但是，目前利用地震資料評估蓋層的方法具有局限性：一方面，依賴于疊后地震數(shù)據(jù)處理的質(zhì)量及準確性，疊后地震屬性精度較低，可能不適用于復(fù)雜地質(zhì)結(jié)構(gòu)和速度變化較大的區(qū)域；另一方面，地震各向異性方法適用于具有各向異性特征的頁巖，而不適用于其他蓋層類型，例如具有非常小孔喉的致密硬石膏、鈣質(zhì)夾層等。目前，地質(zhì)測井實驗室分析測試相結(jié)合的蓋層綜合評價方法已經(jīng)逐漸成熟，但仍缺少地球物理方法如地震反演技術(shù)在蓋層識別方面的運用，特別針對海域油氣蓋層空間展布預(yù)測的成功案例比較匱乏。

地震反演是預(yù)測地下巖性及其空間展布的重要手段，為蓋層的預(yù)測提供了有利的工具。相較機器學(xué)習(xí)等方法對測井?dāng)?shù)據(jù)質(zhì)量和數(shù)量依賴較大的特點［19］，地震反演方法更適合海域油氣勘探初期少井的情況。基于彈性阻抗反演的參數(shù)預(yù)測是目前應(yīng)用較多的地震反演方法［2021］。1999年，Connolly［22］首先提出了彈性阻抗的概念。Whitcombe等［23］提出了擴展彈性阻抗的概念，擴大了彈性阻抗的定義域。在此之后，基于彈性阻抗反演的參數(shù)預(yù)測方法被廣泛提出和運用。王保麗等［24］推導(dǎo)了拉梅參數(shù)表征的彈性阻抗方程，實現(xiàn)了基于彈性阻抗的拉梅參數(shù)直接反演；宗兆云等［25］在此基礎(chǔ)上將拉梅參數(shù)表征的彈性阻抗方程應(yīng)用于碳酸鹽巖縫洞型儲層流體檢測中，并取得了良好的效果；李超等［26］提出了基于入射角的彈性阻抗方程并建立了相應(yīng)的彈性阻抗反演方法；張廣智等［27］通過縱橫波速度彈性阻抗聯(lián)合反演，得到了較為穩(wěn)定的密度；馬妮等［28］基于方位彈性阻抗反演得到了各向異性參數(shù)并預(yù)測了地應(yīng)力。此外，各類反演算法為線性、非線性方程反演奠定了有利基礎(chǔ)。李博南等［29］基于頻變AVO（amplitude variation with offset）理論，采用遺傳算法實現(xiàn)了裂縫密度和流體時間尺度因子的同步非線性反演；鄧煒等［30］基于逆算子估計算法實現(xiàn)了Zoeppritz方程高階近似式的非線性反演；印興耀等［3132］基于貝葉斯理論實現(xiàn)了時頻域聯(lián)合反演，并將高斯混合模型與序貫指示模擬相結(jié)合，準確地將離散變量進行歸類。

針對海域油氣蓋層巖相刻畫和蓋層物性參數(shù)地震識別難度大以及地震預(yù)測穩(wěn)定性差等問題，本文提出一種基于巖石物理約束的海域蓋層巖相以及物性參數(shù)疊前地震預(yù)測方法，分別從巖相預(yù)測以及物性預(yù)測兩個角度，對海域蓋層的地震描述進行深入研究。首先，在巖石物理的基礎(chǔ)上，優(yōu)選巖相以及物性敏感參數(shù)，推導(dǎo)由蓋層敏感參數(shù)直接表征的彈性阻抗方程，為海域蓋層敏感參數(shù)的地震直接反演奠定理論基礎(chǔ)；另外，通過建立針對海域油氣蓋層預(yù)測的巖石物理模型，得到孔隙度與彈性阻抗的定量關(guān)系，為直接預(yù)測蓋層的物性參數(shù)展布規(guī)律提供解決方法。然后，在貝葉斯反演理論的框架下分別推導(dǎo)海域蓋層巖相敏感彈性參數(shù)、物性參數(shù)的反演目標函數(shù)，預(yù)測海域油氣蓋層的橫向展布。最后，通過實際資料的應(yīng)用驗證海域油氣蓋層疊前地震反演預(yù)測方法的合理性與有效性。

1" 技術(shù)方法

1.1" 海域蓋層巖相疊前地震反演

巖相的類型和分布可以幫助預(yù)測海域油氣蓋層的橫向展布及連續(xù)性等特征。通過了解巖石的物理屬性和分布，可以更準確地定位油氣儲層的位置，從而提高勘探的成功率，降低鉆探的風(fēng)險和成本。蓋層的巖相特征直接影響油氣藏的封閉性，進而影響油氣的保存。通過地震預(yù)測可以評估蓋層巖相的連續(xù)性、厚度和巖性，從而判斷其作為有效封閉層的能力，避免油氣泄漏。構(gòu)建蓋層巖相的反演目標函數(shù)是地震反演過程中的關(guān)鍵，主要涉及將地震數(shù)據(jù)與巖石物理模型相結(jié)合，并通過優(yōu)化算法來求解地下巖相分布的最優(yōu)解。首先基于巖石物理分析，優(yōu)選巖相敏感參數(shù)；其次，建立巖相敏感參數(shù)表征的反射系數(shù)及彈性阻抗方程；最后，建立反演目標函數(shù)，通過概率化算法實現(xiàn)蓋層巖相的地震預(yù)測。特別地，由于與傳統(tǒng)AVO反演相比，彈性阻抗反演具有更好的穩(wěn)定性［22］，因此通過稀疏脈沖方法和道積分反演彈性阻抗數(shù)據(jù)體。并且，為了保證地震反演的穩(wěn)定性，采取奇異值分解方法實現(xiàn)模型參數(shù)的去相關(guān)性。具體技術(shù)流程如圖1所示。

蓋層巖相預(yù)測的首要任務(wù)是巖石物理敏感參數(shù)分析，這為后續(xù)的反演提供依據(jù)。使用我國南海某海域油田的數(shù)據(jù)對蓋層敏感巖石物理參數(shù)進行識別和排序。該工區(qū)具有有利儲蓋控富集的油氣成藏模式，發(fā)育了多期、多種類型的儲蓋組合，蓋層主要包括淺湖相泥巖、潮坪泥巖、淺海相泥巖以及主力儲層中普遍發(fā)育的鈣質(zhì)砂巖夾層。研究該工區(qū)對海域油氣蓋層敏感參數(shù)的確定具有重要指導(dǎo)意義。相關(guān)數(shù)據(jù)參考“2.1巖石物理分析”的實際測井資料。我們通過確定蓋層巖相與其他巖相在各種巖石物性參數(shù)中的差異，直觀地顯示每個巖石物性參數(shù)對蓋層的敏感性。巖性差異越高，代表該參數(shù)對蓋層越敏感，能夠更好地區(qū)分蓋層。圖2展示了各種彈性參數(shù)對蓋層的敏感性分析，其中泊松阻抗的巖性差異最高（圖2），表明該參數(shù)區(qū)分鈣質(zhì)砂巖的能力較好，泊松阻抗對蓋層的敏感性最高；這表明，與常規(guī)縱橫波速度、密度、彈性阻抗、縱波阻抗等相比，泊松阻抗在蓋層預(yù)測中更具優(yōu)勢。

Zong等［33］推導(dǎo)了基于泊松阻抗的縱波（P波）近似反射系數(shù)方程：

RPP（θ）=14sec2θ－2ksin2θΔIYIY+14sec2θ·

（2k－3）（2k－1）2k（4k－3）+2ksin2θ1－2k3－4kΔIPIP+12－14sec2θ4k2k2－3k+2－34k2－3k+

4ksin2θ2－3k3－4kΔρρ。（1）

式中：RPP為P波反射系數(shù)；θ為入射角；IY為楊氏阻抗；k為橫波速度與縱波速度之比的平方；Δ表示變化量。當(dāng)波阻抗變化較小時，以下近似是合理的：

Rpp≈12ΔIeIe≈12Δln Ie。（2）

由式（1）（2）得

Ie（θ）=Ia（θ）YIb（θ）Pρc（θ）。（3）

其中：

a（θ）=12sec2θ－4ksin2θ；b（θ）=12sec2θ（2k－3）（2k－1）2k（4k－3）+4ksin2θ1－2k3－4k；c（θ）=1－12sec2θ4k（2k2－3k+2）－34k2－3k+""""""" 8ksin2θ2－3k3－4k。（4）

式（3）表示了楊氏阻抗、泊松阻抗、密度與彈性

阻抗之間的非線性指數(shù)關(guān)系。為了方便求解，式（3）

λ. 拉梅常數(shù)；vP. 縱波速度；vS. 橫波速度；ρ. 密度；IP. 泊松阻抗；Ie.彈性阻抗。

被標準化并簡化為線性形式：

ln IeIe0=12sec2θ－4ksin2θIYIY0+" 12sec2θ（2k－3）（2k－1）2k（4k－3）+

4ksin2θ1－2k3－4kIPIP0+" 1－12sec2θ4k2k2－3k+2－34k2－3k+

8ksin2θ2－3k3－4kρρ0。（5）

式中，Ie0、IY0、IP0和ρ0分別為彈性阻抗、楊氏阻抗、泊松阻抗和密度的平均值。通過三個不同入射角度的彈性阻抗可構(gòu)建彈性阻抗與楊氏阻抗、泊松阻抗和密度的關(guān)系：

lnIeθ1Ie0=a（θ1）IYIY0+b（θ1）IPIP0+c（θ1）ρρ0；lnIe（θ2）Ie0=a（θ2）IYIY0+b（θ2）IPIP0+c（θ2）ρρ0；lnIe（θ3）Ie0=a（θ3）IYIY0+b（θ3）IPIP0+c（θ3）ρρ0。（6）

假設(shè)地層的反射系數(shù)序列是稀疏分布的，則估計彈性阻抗的目標函數(shù)可以寫成

J=minL11（rEI）+α1L22［d－G（m）］+

α2L22（t－Ie）+α3L22（SymbolQC@Ie）。（7）

式中：Lji（i=1，2，j=1，2）代表i范數(shù)的j次方；rEI為彈性阻抗反射系數(shù)；αj（j=1， 2， 3）為權(quán)重系數(shù)；d為實際地震記錄；G（m）為合成地震記錄；m為速度、密度等各種與地震傳播相關(guān)的參數(shù)；t為彈性阻抗低頻模型；Ie為彈性阻抗；SymbolQC@Ie為測井計算得到的彈性阻抗梯度。通過求解式（7）獲得彈性阻抗反射系數(shù)，再利用道積分得到最終的彈性阻抗：

Ie=Ie0exp∫tt0rτdτ。（8）

式中：Ie0為初始彈性阻抗；t為采樣時間；t0為初始時間。

為了解決模型參數(shù)之間相關(guān)性導(dǎo)致反演結(jié)果不穩(wěn)定的問題，采用奇異值分解方法實現(xiàn)模型參數(shù)的去相關(guān)［34］。假設(shè)模型參數(shù)的協(xié)方差矩陣為

Cr=σ2m1σm1m2σm1m3σm2m1σ2m2σm2m3σm1m3σm2m3σ2m3。（9）

式中：σ2m1為方差；σm2m3為協(xié)方差。通過奇異值分解可以寫成

Cr=uΣuT=uσ21σ22σ23uT。（10）

式中：u為由特征向量組成的單位正交矩陣；Σ為包含奇異值平方的對角矩陣；σ21、σ22、σ23為奇異值的平方。稀疏協(xié)方差矩陣可以寫為

e－1=e110…e120…e130…0e110…e120…e130e210…e220…e230…0e210…e220…e230e310…e320…e330…0e310…e320…e3302t×2t。（11）

因此，地球物理問題的正演算子可以寫成

S=A′R′。（12）

其中：

A′=Ae；R′=e-1R。（13）

式中：A′為去相關(guān)待反演參數(shù)系數(shù)矩陣；R′為去相關(guān)待反演參數(shù)矩陣。

貝葉斯反演方法是一種基于概率與統(tǒng)計的不確定性參數(shù)估計方法。該方法利用待反演參數(shù)的先驗分布信息，通過求解模型參數(shù)的后驗概率實現(xiàn)模型參數(shù)的預(yù)測，在處理線性和非線性反演問題中應(yīng)用較廣。貝葉斯公式可以表示為［35］

P（R′|S）=" P（R′）P（S|R′）∫P（R′）P（S|R′）dR′∝p（R′）P（S|R′）。（14）

式中：P（R′|S）為后驗概率；P（R′）為待估計參數(shù)的先驗分布；P（S|R′）為似然函數(shù)；∫P（R′）P（S|R′）dR′為觀測樣本的全概率，可近似為常數(shù)。為提高反演的穩(wěn)定性，可假設(shè)待估計參數(shù)滿足柯西分布［36］：

PCauchy（R′）=1（πσm）N∏Ni=111+R′i2/σm2。（15）

式中：σm2為模型參數(shù)的方差；N為待估計參數(shù)的數(shù)量；R′i為R′的第i個分量。假設(shè)似然函數(shù)符合高斯分布，σn2是隨機噪聲的方差，則似然函數(shù)可表示為

P（S|R′）∝exp－（A′R′－S）T（A′R′－S）2σn2。（16）

綜合式（14）（15）（16），待反演參數(shù)的后驗概率可表示為

P（Ρ′|S）∝exp［－∑Ni=1ln（1+R′i2/σm2）］·

exp－（A′R′－S）T（A′R′－S）2σn2。（17）

將模型約束加入式（17），并求取其最大值可得到估計楊氏阻抗、泊松阻抗、密度的最終目標函數(shù)：

F=（A′R′－S）T（A′R′－S）+2σn2∑Ni=1ln（1+R′i2/σm2）+

B12lnR′（t）R′（t0）－∫tt0R′（τ）dτ。（18）

式中，B為模型參數(shù)的約束系數(shù)。通過迭代重加權(quán)最小二乘方法［37］求解式（18），即可得到最終解。

1.2" 基于巖石物理的蓋層物性參數(shù)疊前地震反演

現(xiàn)有研究缺少針對海域油氣蓋層的理論巖石物理模型，難以充分反映蓋層彈性阻抗與物性參數(shù)的關(guān)系。因此基于巖石物理建模，充分利用工區(qū)的實際地質(zhì)特征，最大程度地得到彈性阻抗與物性參數(shù)之間的關(guān)系。以下簡述巖石物理建模過程：

第一步，計算各向同性巖石基質(zhì)等效彈性模量。經(jīng)調(diào)研，研究工區(qū)的礦物以黏土、石英、鉀長石、方解石為主，含少量白云石。采用VRH平均模型［38］計算巖石基質(zhì)等效彈性模量：

MVRH=12∑ni=1Mifi+∑ni=1fiMi－1。（19）

式中：MVRH為巖石基質(zhì)等效彈性模量；Mi為第i種固體礦物的彈性模量；fi為第i種礦物的體積分數(shù)；n為固體礦物的數(shù)目。

第二步，耦合孔隙，形成背景骨架。采用KT模型［39］進行孔隙耦合，通過形狀因子控制孔隙的縱橫比：

（KKT－Km）Km+43μmKKT+43μm=∑i=H/Sxi（Ki－KKT）Pmi；（20）

（μKT－μm）μm+ζmμKT+ζm=∑i=H/Sxi（μi－μKT）Qmi。（21）

其中：

ζm=μm69Km+8μmKm+2μm。（22）

式中：K*KT和μ*KT分別為有效介質(zhì)的體積模量和剪切模量；Km和μm分別為基質(zhì)的體積模量和剪切模量；S和H分別為在孔隙壓力下形變的軟孔隙和不宜形變的硬孔隙；xi為包含物材料i的體積分數(shù)；Ki和μi分別為包含物的體積模量和剪切模量；Pmi和Qmi分別為與包含物模量和形狀有關(guān)的極化因子，描述在基質(zhì)中加入包含物材料i后的效果。

第三步，利用含流體孔隙介質(zhì)理論添加孔隙流體，形成含流體孔隙巖石。各向異性情況下的流體替換公式［40］是Gassmann方程的各向異性形式：

Sdryijkl－Ssatijkl=Sdryijαα－SmijααSdrykααα－SmklααSdryααβββ－Smααββ+βfl－βmφ。（23）

式中：Sijkl為巖石的柔度；β為可壓縮量；角標sat、dry、m和fl分別代表飽和巖石、干骨架、基質(zhì)和流體；φ為孔隙度。當(dāng)巖石骨架滿足各向同性假設(shè)時，式（23）可以寫成：

K－1sat=K－1dry－K－1dry－K－1m2φK－1i－K－1φ+K－1dry－K－1m；μ－1sat=μ－1dry－μ－1dry－μ－1m2φμ－1i－μ－1φ+μ－1dry－μ－1m。（24）

式中：Ksat、μsat分別為飽和巖石的體積模量和剪切模量；Kdry、μdry分別為干骨架的體積模量和剪切模量；Ki、μi分別為第i種孔隙填充材料的體積模量和剪切模量；Kφ、μφ分別為孔隙空間的體積模量和剪切模量。

因此，飽和流體巖石的縱、橫波速度可表示為：

vP=Ksat+43μsatρ；（25）

vS=μsatρ。（26）

基于上述地震巖石物理模型，將孔隙度等物性參數(shù)表征的飽和流體巖石彈性參數(shù)代入歸一化彈性阻抗方程，可得到孔隙度等物性參數(shù)與彈性阻抗的定量關(guān)系。其中，物性參數(shù)中孔隙度對彈性阻抗和AVO反射系數(shù)的貢獻度最大［41］，因此可以利用地震數(shù)據(jù)對蓋層孔隙度進行定量地震預(yù)測，采取貝葉斯反演方法，算法與前文類似。具體步驟如圖3所示。

2" 應(yīng)用實例

2.1" 巖石物理分析

巖石物理分析是優(yōu)選蓋層敏感參數(shù)的重要手段。以中國南海某海域工區(qū)的實際測井資料為例，圖4顯示了工區(qū)內(nèi)井上巖性解釋結(jié)果及物性參數(shù)真實值。從井上選取具有代表性的一組儲蓋組合，即“a蓋層 b儲層 c蓋層”，可以看出：儲層以砂巖為主，孔隙度較高；蓋層以泥巖、鈣質(zhì)砂巖為主，孔隙度較低。

圖5顯示了工區(qū)內(nèi)井上彈性參數(shù)的真實值，可以看到：可能由于地層壓實等作用，一些彈性參數(shù)

的值隨著深度的增加而逐漸增加，這會降低反演結(jié)果的可解釋性；而隨著深度的增加，泊松阻抗未受壓實左右而顯著變小，有助于提高巖性判別的可靠性。另外，可以看出：在a蓋層、b儲層和c蓋層位置，縱波速度、密度、縱波阻抗、彈性阻抗規(guī)律不清晰，蓋層同時對應(yīng)著各參數(shù)的高值與低值，難以有效區(qū)分儲層與蓋層；而泊松阻抗對蓋層有著較好的敏感性，主要表現(xiàn)為儲層泊松阻抗較低、蓋層泊松阻抗較高。

為了驗證泊松阻抗識別蓋層的可靠性，圖6a、b分別定量展示了泊松阻抗與方解石體積分數(shù)和泥質(zhì)體積分數(shù)的巖性交會圖，可以看出，泊松阻抗可以較好地區(qū)分蓋層（藍色、黑色）和儲層（黃色）。圖6c展示了泊松阻抗與彈性阻抗的巖性交會圖，可以看出，泊松阻抗識別蓋層的能力明顯高于彈性阻抗。圖6d展示了泊松阻抗與孔隙度的巖性交會圖，可以看出，蓋層具有低孔隙度和高泊松阻抗的特征。因此，可以采取通過泊松阻抗巖相刻畫、孔隙度定量預(yù)測相結(jié)合的方法對海域油氣蓋層進行識別。

2.2" 海域蓋層疊前地震反演預(yù)測方法的應(yīng)用

常規(guī)泊松阻抗反演方法利用泊松阻抗與彈性阻抗之間的統(tǒng)計關(guān)系進行反演，或通過縱橫波速度與密度的提取計算泊松阻抗，存在物理意義不明確、統(tǒng)計性關(guān)系適用范圍較小及誤差累計等問題。針對這些問題，依據(jù)前文建立的泊松阻抗；表征的彈性阻抗方程，直接反演得到泊松阻抗。該方法物理意義明確、試用范圍大、累計誤差小。在海域蓋層物性參數(shù)預(yù)測方面，依據(jù)前文建立的針對海域油氣蓋層預(yù)測的巖石物理模型，得到彈性阻抗與物性參數(shù)之間的定量關(guān)系，并基于彈性阻抗實現(xiàn)孔隙度反演。

為了驗證本文提出的海域油氣蓋層疊前地震反演預(yù)測方法的適用性，將該方法應(yīng)用于中國南海海域某工區(qū)。圖7為過A井和B井的三個角度部分疊加地震剖面。其中：A為已知井，根據(jù)測井巖性解釋，已知存在三個較好蓋層；B井是驗證井，是蓋層

a. 8°；b. 18°；c. 28°。

預(yù)測方法應(yīng)用后新鉆的井，未參與反演，也存在三個較好蓋層。

圖8分別展示了泊松阻抗和孔隙度反演結(jié)果剖面。在圖8a中，剖面的藍色代表泊松阻抗較高，指示蓋層，剖面的黃色代表泊松阻抗較低，指示儲層。A井展示的是泊松阻抗井曲線，可以看出反演結(jié)果與井曲線匹配較好；B井是驗證井，井上藍色代表測井解釋為蓋層，黃色代表測井解釋為儲層，可以看出泊松阻抗反演剖面與驗證井匹配良好，說明基于彈性阻抗反演的泊松阻抗預(yù)測方法在蓋層定性刻畫蓋層空間展布方面是有效的。在圖8b中，孔隙度反演結(jié)果與A井孔隙度曲線對應(yīng)較好，與驗證井B井的測井巖性解釋結(jié)果也有很好的吻合，說明基于彈性阻抗反演的孔隙度預(yù)測方法在定量刻畫蓋層空間展布方面是有效的。

圖9展示了A井位置泊松阻抗、孔隙度井旁道對比及反演誤差分析，可以看出反演結(jié)果與真實值吻合良好，誤差在合理范圍。

平均絕對百分比誤差是指將預(yù)測值與實際值之間的絕對誤差轉(zhuǎn)化為百分比的平均值，該參數(shù)可以較好地評估預(yù)測值的準確性和誤差程度。平均絕對百分比誤差越小，說明預(yù)測準確性越高。其計算公式為

EMAP=1n∑ni=1yi－y^iyi。（27）

式中：n為值的數(shù)量；yi為實際值；y^i為預(yù)測值。為了進一步定量分析預(yù)測精度，我們計算得到了井旁道泊松阻抗預(yù)測值和孔隙度預(yù)測值的平均絕對百分比誤差，分別為4.91%和14.88%。由此，我們可以認為，泊松阻抗的預(yù)測精度可以達到95%，孔隙度的預(yù)測精度可以達到85%。

圖10顯示了海域蓋層泊松阻抗和孔隙度的反演結(jié)果切片，工區(qū)內(nèi)另一口驗證井C井的巖性綜合

解釋結(jié)果在切片左側(cè)展示。通過圖10左側(cè)的測井綜合解釋可以看出，C井存在明顯的油氣蓋層，蓋層處鈣質(zhì)砂巖體積分數(shù)較高，蓋層下含油體積分數(shù)較高，是良好的儲層。通過反演切片可以看到，在C井位置處，泊松阻抗為高值、孔隙度為低值，反演結(jié)果可以有效識別驗證井C井的油氣蓋層，這與巖石物理交會分析結(jié)果一致，進一步驗證了海域油氣

蓋層疊前地震反演預(yù)測方法的有效性和適用性。另外，反演結(jié)果可以較好地指示蓋層所在位置及其橫向展布，能夠為后續(xù)海域油氣的勘探開發(fā)提供幫助。

3" 結(jié)論

針對海域油氣蓋層地震預(yù)測方法較少、預(yù)測不準確、橫向展布刻畫較難等問題，本文提出了一種基于巖石物理約束的海域油氣蓋層巖相以及物性參數(shù)疊前地震預(yù)測方法。基于本文的研究得到以下結(jié)論：

1）通過海域地層的巖石物理建模分析及敏感參數(shù)優(yōu)選，發(fā)現(xiàn)泊松阻抗對海域油氣蓋層巖相更加敏感，孔隙度可以定量表征蓋層的物性。

2）在蓋層巖相刻畫方面，利用泊松阻抗表征的彈性阻抗方程，在貝葉斯理論的框架下可以實現(xiàn)泊松阻抗直接反演。在蓋層物性定量預(yù)測方面，通過建立針對海域油氣蓋層預(yù)測的巖石物理近似模型和物性參數(shù)解析表征，能夠?qū)崿F(xiàn)孔隙度預(yù)測。

3）實際資料驗證了所提方法在海域油氣蓋層識別方面的可行性。泊松阻抗預(yù)測精度可以達到95%，孔隙度預(yù)測精度可以達到85%。

4）本方法的提出為海域油氣蓋層預(yù)測提供了新方法和可靠的地球物理依據(jù)，能夠為海域油氣勘探與開發(fā)提供有利指導(dǎo)。

參考文獻（References）：

［1］" 周雁，金之鈞，朱東亞，等. 油氣蓋層研究現(xiàn)狀與認識進展［J］. 石油實驗地質(zhì)，2012，34（3）：234245.

Zhou Yan， Jin Zhijun， Zhu Dongya， et al. Current Status and Progress in Research of Hydrocarbon Cap Rocks［J］. Petroleum Geology amp; Experiment， 2012， 34（3）： 234245.

［2］" 付廣，劉江濤. 中國高效大中型氣田形成的封蓋保存條件［J］. 石油勘探與開發(fā)，2006，33（6）：662666.

Fu Guang， Liu Jiangtao. Sealing and Preservation Conditions for Large and Medium Gas Fields of China［J］. Petroleum Exploration and Development， 2006， 33（6）： 662666.

［3］" Dott R H， Reynolds M J. Source Book for Petroleum Geology：" Semicentennial Commemorative Volume［M］. Tulsa： American Association of Petroleum Geologists， 1969： 1471.

［4］" 劉宗賓，程奇，呂坐彬，等. 渤海海域錦州南太古宇變質(zhì)巖潛山儲層特征及其發(fā)育模式［J］. 吉林大學(xué)學(xué)報（地球科學(xué)版）， 2022， 52（3）：829839.

Liu Zongbin，Cheng Qi，Lü Zuobin，et al. Reservoir Characteristics and Development Model of Archean Metamorphic Buried Hill in Jinzhou South Oilfield， Bohai Sea［J］. Journal of Jilin University （Earth Science Edition）， 2022， 52（3）：829839.

［5］" 楊計海，楊希冰，游君君，等. 珠江口盆地珠三坳陷油氣成藏規(guī)律及勘探方向［J］. 石油學(xué)報，2019，40（增刊1）：1125.

Yang Jihai， Yang Xibing， You Junjun， et al. Hydrocarbon Accumulation Law and Exploration Direction in Zhu Ⅲ Depression， Pearl River Mouth Basin［J］. Acta Petrolei Sinica， 2019， 40（Sup.1）： 1125.

［6］" 戴建文，柴愈坤，王冠民，等. 珠江口盆地恩平凹陷A油田新近系儲層鈣質(zhì)砂巖夾層成因［J］. 海相油氣地質(zhì)，2023，28（4）：391400.

Dai Jianwen， Chai Yukun， Wang Guanmin， et al.Genesis of Calcareous Sandstone Interlayer in Neogene Reservoir of A Oilfield in Enping Sag， Pearl River Mouth Basin［J］. Marine Origin Petroleum Geology， 2023， 28（4）： 391400.

［7］" 謝玉洪，張功成，沈樸，等. 渤海灣盆地渤中凹陷大氣田形成條件與勘探方向［J］. 石油學(xué)報，2018，39（11）：11991210.

Xie Yuhong， Zhang Gongcheng， Shen Pu， et al. Formation Conditions and Exploration Direction of Large Gas Field in Bozhong Sag of Bohai Bay Basin［J］. Acta Petrolei Sinica， 2018， 39（11）： 11991210.

［8］" 胡楊，張偉，張景茹，等. 南海北部大陸邊緣瓊東南盆地有利油氣聚集帶及勘探方向［J］. 海洋地質(zhì)與第四紀地質(zhì)，2015，35（4）：7386.

Hu Yang， Zhang Wei， Zhang Jingru， et al.The Favorable Petroleum Accumulation Belts and Exploration Targets in Qiongdongnan Basin on the Margin of Northern South China Sea［J］. Marine Geology amp; Quaternary Geology， 2015， 35（4）： 7386.

［9］" 張為彪，易浩，鐘輝，等. 珠江口盆地東部地區(qū)預(yù)探井失利原因分析及啟示［J］. 中國石油勘探，2018，23（3）：1827.

Zhang Weibiao， Yi Hao， Zhong Hui， et al. Analysis on the Causes of Failed Wildcat Wells in the Eastern Pearl River Mouth Basin and Its Enlightenment［J］. China Petroleum Exploration， 2018， 23（3）： 1827.

［10］" Downey M W. Evaluating Seals for Hydrocarbon Accumulation［J］. AAPG Bulletin， 1984， 68（11）： 17521763.

［11］" Grunau H R. Worldwide Review of Seals for Major Accumulations of Natural Gas： ABSTRACT［J］. AAPG Bulletin， 1981， 65（5）： 933．

［12］" 張樹林，田世澄. 蓋層的研究方法及其在油氣勘探中的意義［J］. 地質(zhì)科技情報，1993，12（1）：7378.

Zhang Shulin， Tian Shicheng. Methods for Reserch of Caprocks and Their Significance in Hydrocarbon Exploration［J］. Geological Science and Technology Information， 1993， 12（1）： 7378.

［13］" 付廣，康德江，段海鳳. 泥質(zhì)巖蓋層對水溶相天然氣封閉能力的綜合評價方法及其應(yīng)用［J］. 吉林大學(xué)學(xué)報（地球科學(xué)版），2006，36（1）：6065.

Fu Guang， Kang Dejiang， Duan Haifeng. A Comprehensive Evaluation Method on Seal Ability of Mudstone Caprock to Gas in Water-Dissolving Phase and Its Application［J］. Journal of Jilin University （Earth Science Edition）， 2006， 36（1）： 6065.

［14］" 付廣，胡明，于丹. 火山巖蓋層類型及封氣能力：以松遼盆地徐家圍子斷陷為例［J］. 吉林大學(xué)學(xué)報（地球科學(xué)版），2010，40（2）：237244.

Fu Guang， Hu Ming， Yu Dan.Volcanic Cap Rock Type and Evaluation of Sealing Gas Ability： An Example of Xujiaweizi Depression［J］. Journal of Jilin University （Earth Science Edition）， 2010， 40（2）： 237244.

［15］" Wu S Y， Sun S Q， Pollitt D A. Quantitative Evaluation of Hydrocarbon Containment Using Seal Strength Index （SSI） and Benchmarking Against Global Analogues［C］// Aria A， Aiman M B. Third International Meeting for Applied Geoscience amp; Energy. Houston： Society of Exploration Geophysicists and American Association of Petroleum Geologists， 2023： 526530.

［16］" Davis T L， Benson R D， Stroker T. Seal Identification from Time-Lapse Multicomponent Seismic Data， Rulison Field， Piceance Basin， Colorado［C］// John C. SEG Technical Program Expanded Abstracts. Houston： Society of Exploration Geophysicists， 2009： 18751879.

［17］" Nourollah H， Keetley J， O’Brien G. Gas Chimney Identification Through Seismic Attribute Analysis in the Gippsland Basin， Australia［J］. The Leading Edge， 2010， 29（8）： 896901.

［18］" Nourollah H， Urosevic M， Keetley J. Seal Potential of Shale Sequences Through Seismic Anisotropy： Case Study from Exmouth Sub-Basin， Australia［J］. Interpretation， 2015， 3（4）： T257T267.

［19］" 侯賢沐，王付勇，宰蕓，等. 基于機器學(xué)習(xí)和測井?dāng)?shù)據(jù)的碳酸鹽巖孔隙度與滲透率預(yù)測［J］. 吉林大學(xué)學(xué)報（地球科學(xué)版），2022，52（2）：644653.

Hou Xianmu， Wang Fuyong， Zai Yun， et al.Prediction of Carbonate Porosity and Permeability Based on Machine Learning and Logging Data［J］. Journal of Jilin University （Earth Science Edition）， 2022， 52 （2）： 644653.

［20］" 欒穎，馮晅，劉財，等. 波阻抗反演技術(shù)的研究現(xiàn)狀及發(fā)展［J］ . 吉林大學(xué)學(xué)報（地球科學(xué)版），2008，38（增刊1）：9498.

Luan Ying， Feng Xuan， Liu Cai， et al. The Research Present and Future of Wave Impedance Inversion Technique［J］. Journal of Jilin University （Earth Science Edition）， 2008， 38（Sup. 1）： 9498.

［21］" Ji L X， Zong Z Y. Lithology Discrimination Based on Direct Inversion of Poisson Impedance for Deep Tight-Sandstone Reservoirs［J］. Journal of Geophysics and Engineering， 2023， 20（1）： 3848.

［22］" Connolly P. Elastic Impedance［J］. The Ieading Edge， 1999， 18（4）： 438452.

［23］" Whitcombe D N， Connolly P A， Reagan R L， et al. Extended Elastic Impedance for Fluid and Lithology Prediction［J］. Geophysics， 2002， 67（1）： 6367.

［24］" 王保麗，印興耀，張繁昌. 基于Gray近似的彈性波阻抗方程及反演［J］. 石油地球物理勘探，2007，42（4）：435439.

Wang Baoli， Yin Xingyao， Zhang Fanchang. Gray Approximation： Based Elastic Wave Impedance Equation and Inversion ［J］. Oil Geophysical Prospecting， 2007， 42（4）： 435439.

［25］" 宗兆云，印興耀，吳國忱. 拉梅參數(shù)直接反演技術(shù)在碳酸鹽巖縫洞型儲層流體檢測中的應(yīng)用［J］. 石油物探，2011，50（3）：241246.

Zong Zhaoyun， Yin Xingyao， Wu Guochen. Application of Lame Parameter Direct Inversion in Fluid Identification for Carbonate Reservoir［J］. Geophysical Prospecting for Petroleum， 2011， 50（3）： 241246.

［26］" 李超，印興耀，張廣智，等. 基于入射角的AVO近似彈性波阻抗方程及反演［J］. 吉林大學(xué)學(xué)報（地球科學(xué)版），2013，43（5）：16631671.

Li Chao， Yin Xingyao， Zhang Guangzhi， et al. Elastic Impedance Equation Based on the Incident-Angle AVO Approximation and Inversion［J］. Journal of Jilin University （Earth Science Edition）， 2013， 43（5）： 16631671.

［27］" 張廣智，杜炳毅，陳懷震，等. 縱橫波彈性阻抗聯(lián)合反演方法［J］. 吉林大學(xué)學(xué)報（地球科學(xué)版），2014，44（5）：16951704.

Zhang Guangzhi， Du Bingyi， Chen Huaizhen， et al. Joint Inversion of PP and PS Waves Elastic Impedance［J］. Journal of Jilin University （Earth Science Edition）， 2014， 44（5）： 16951704.

［28］" 馬妮，印興耀，孫成禹，等. 基于方位地震數(shù)據(jù)的地應(yīng)力反演方法［J］. 地球物理學(xué)報，2018，61（2）：697706.

Ma Ni， Yin Xingyao， Sun Chengyu， et al. Inversion for Crustal Stress Based on Azimuthal Seismic Data［J］. Chinese Journal of Geophysics， 2018， 61（2）： 697706.

［29］" 李博南，劉財，郭智奇. 基于等效介質(zhì)模型和頻變AVO反演的裂縫儲層參數(shù)估算方法［J］. 吉林大學(xué)學(xué)報（地球科學(xué)版），2017，47（1）：234244.

Li Bonan， Liu Cai， Guo Zhiqi. Estimation of Fractured Reservoir Parameters Based on Equivalent Media Model and Frequncy-Dependent AVO Inversion［J］. Journal of Jilin University （Earth Science Edition）， 2017， 47（1）： 234244.

［30］" 鄧煒，印興耀，宗兆云，等. 基于逆算子估計的高階AVO非線性反演［J］. 石油地球物理勘探，2016，51（5）：955964.

Deng Wei， Yin Xingyao， Zong Zhaoyun， et al.High-Order Nonlinear AVO Inversion Based on Estimated Inverse Operator［J］. Oil Geophysical Prospecting， 2016， 51（5）： 955964.

［31］" 印興耀，李坤，宗兆云，等. 時頻聯(lián)合域貝葉斯地震反演方法［J］. 石油物探，2017，56（2）：250260.

Yin Xingyao， Li Kun， Zong Zhaoyun， et al. Seismic Inversion in Joint Time-Frequency Domain Based on Bayesian Scheme［J］. Geophysical Prospecting for Petroleum， 2017， 56（2）： 250260.

［32］" 印興耀，賀東陽，宗兆云，等. 結(jié)合序貫指示模擬的貝葉斯高斯混合線性反演［J］. 地球物理學(xué)報，2019，62（12）：48054816.

Yin Xingyao， He Dongyang， Zong Zhaoyun， et al. Bayesian Gaussian Mixture Linear Inversion Combined with Sequential Indicator Simulation［J］. Chinese Journal of Geophysics， 2019， 62（12）： 48054816.

［33］" Zong Z Y， Yin X Y. Direct Inversion of Young’s and Poisson Impedances for Fluid Discrimination［J］. Geofluids， 2016， 16（5）： 10061016.

［34］" Ji L X， Zong Z Y. Lithology Discrimination Based on Direct Inversion of Poisson Impendance for Tight-Stone Reservoirs［J］. Journal of Geophysics and Engineering， 2023， 20： 3848.

［35］" Bayes T. XLIII： A Letter from the Late Reverend Mr Thomas Bayes， F R S to John Canton， M A and F R S［J］. Philosophical Transactions of the Royal Society of London， 1763， 53： 269271.

［36］" Zong Z Y， Wang Y R， Li K， et al. Broadband Seismic Inversion for Low-Frequency Component of the Model Parameter［J］. IEEE Transactions on Geoscience and Remote Sensing， 2018， 56（9）： 51775184.

［37］" Daubechies I， DeVore R， Fornasier M， et al. Iteratively Reweighted Least Squares Minimization for Sparse Recovery［J］. Communications on Pure and Applied Mathematics： A Journal Issued by the Courant Institute of Mathematical Sciences， 2010， 63（1）： 138.

［38］" Hill R. The Elastic Behaviour of a Crystalline Aggregate［J］. Proceedings of the Physical Society： Section A， 1952， 65（5）： 349.

［39］" Kuster G T， Toksz M N. Velocity and Attenuation of Seismic Waves in Two-Phase Media： Part I： Theoretical Formulations［J］. Geophysics， 1974， 39（5）： 587606.

［40］" Brown R J S， Korringa J. On the Dependence of the Elastic Properties of a Porous Rock on the Compressibility of the Pore Fluid［J］. Geophysics， 1975， 40（4）： 608616.

［41］" 李紅兵，張佳佳，潘豪杰，等. 基于彈性阻抗的孔隙結(jié)構(gòu)與物性參數(shù)非線性同步反演［J］. 中國科學(xué)： 地球科學(xué)，2021，51（7）：11661180.

Li Hongbing， Zhang Jiajia， Pan Haojie， et al. Nonlinear Simultaneous Inversion of Pore Structure and Physical Parameters Based on Elastic Impedance［J］. Science China： Earth Sciences， 2021， 51（7）： 11661180.