1、第五章 油氣綜合檢測技術第一節 疊前流體檢測技術近幾年，隨著地震采集處理技術的進步，尤其疊前偏移技術的發展和推廣應用，使得研究人員可以得到來自地下真實反射點的疊前道集（CRP道集），為疊前烴類檢測技術的發展奠定了資料基礎。目前基于疊前道集的直接烴類檢測方法主要有兩種：一種是在巖石物理建模的基礎上進行疊前道集AVO響應特征分析；一種是利用多個限角疊加數據體進行疊前彈性參數反演，利用縱橫波波阻抗、縱橫波速度比、泊松比、拉梅系數等敏感屬性反映含油氣性。一. AVO分析技術1、 AVO理論簡介AVO（Amplitude Variation with Offset），早先也稱為Amplitude Ver

2、sus Offset，譯為振幅隨炮檢距變化。由此而衍生的有振幅隨入射角變化AVO（Amplitude Variation with Angle），振幅隨方位角變化AVA（Amplitude Variation with Azimuth），振幅隨炮檢距和方位角變化AVOA（Amplitude Variation with Offset and Azimuth）等。AVO作為一種含氣砂巖的異常地球物理現象，最早在20世紀80年代初被Ostrander發現。這一現象表現為：當儲層砂巖含氣后，地震反射振幅隨炮檢距會發生明顯的加大（基于SEG標準極性）。因為AVO現象與含氣砂巖的對應關系，從而引起勘探地

3、球物理界廣泛的重視。后續的研究表明：這種異常現象并非一種特殊的形式，而是遵循Zoepprittz早先所提出的地震反射波動力學方程式，從而對AVO現象的解釋有了完整的理論基礎。針對AVO現象繼而出現的AVO技術是繼亮點之后又一項利用振幅信息研究巖性、檢測油氣的技術手段。AVO技術具有以下特點：A、直接利用CDP道集資料進行分析，這就充分利用的多次覆蓋得到的豐富的原始信息；B、利用振幅隨炮檢距（入射角）的變化的特點，即利用整條曲線的特點。而亮點技術只是利用了這一特殊情況下曲線的一個數值。所以，AVO技術對巖性的分析比亮點技術更為可靠。C、這幾年波動方程對地震剖面的成像有了更大的成果，是對地下構造形

4、態的反演。AVO技術從嚴格意義上說算不上是利用波動方程進行巖性反演分析的方法，但是其理論和思路是對波動方程得到的結果的比較精確的利用。D、AVO技術是一種研究巖性的比較細致的方法，并且需要有測井資料的配合。2、 AVO技術的理論基礎振幅隨炮檢距的變化來自于所謂的“能量分區”。當地震波入射到地層界面時，一部分能量反射，一部分能量透射。如果入射角不等于零度，縱波（P波）能量一部分反射，一部分轉化成透射P波和S波。反射和透射波的振幅能量取決于地層邊界的物理性質差異。縱波速度Vp、橫波速度Vs和密度是非常重要的。同時，需要注意反射振幅也依賴于入射波的入射角（圖5-3-1）。因此，當一個平面縱波非垂直入

5、射到兩種介質的分界面上，就要產生反射縱、橫波和透射縱、橫波。在界面上，根據應力連續性和位移連續性，依據邊界條件并引入反射系數、透射系數，就可以得出四個相應波的位移振幅應當滿足的方程叫做Zoeprritz方程，這個方程是Zoeprritz在1919年解出的。這個方程組比較復雜，不能解出新產生的波的振幅與有關參數明確的函數關系。但是從方程組可以看出，一般反射縱波的反射系數Rpp是入射角界面上部介質的密度1，縱波速度Vp1，橫波速度Vs1以及界面以下的介質密度2，縱波速度Vp2，橫波速度Vs2等七個參數的函數，可以簡單的表示為Rpp（，Vp1，Vs1，1，Vp2，Vs2，2），雖然不能直接從方程中解

6、出Rpp與七個參數的具體關系，但是可以假設以物質的六個物性參數為參變量，以為變量，仔細分析可以得到，六個參數是以兩個參數的比值，例如、等形式出現。這樣就可以把、等分別看作一個參數，再加上在同一種介質中，縱波速度p，橫波速度s，以及泊松比之間又有關系，如：，于是有關系式：Rpp=f（，），這樣來達到減少參數的目的。從理論上說，在實際地震記錄上得到某個界面的反射波的振幅與入射角的變化關系曲線，并且又知道某些參數，就可以利用曲線族作為量板來估算地層參數。圖5-3-1：入射到地層邊界的地震波示意圖通過對地層彈性參數的研究，得到泊松比是一個對巖性和含油氣情況反應比較敏感，參數，所以就要對方程進行適當的化

7、簡，得出以泊松比為參數的以為變量的簡單的近似關系，即Rpp=f（，），這樣通過反算，來達到對儲層的參數測定和檢測。由于Zoeppritz方程過于復雜，因此有許多學者嘗試對其進行簡化或近似，其中比較著名和實用的主要有：1）Koefoed的試算Koefoed在1955年第一個給出了用Zoeppritz方程，以R為參數計算出的Rpp曲線。他用17組縱橫波速度、密度和泊松比參數，較為詳細地研究了泊松比對兩個各向同性介質之間反射和折射面所產生的反射系數的影響，最大的入射角達到30°。他的研究結果被公認為Koefoed五原則。雖然Koefoed的結論說明了利用Rpp曲線是可以反算出泊松比的，但是

8、用未簡化的Zoeppritz方程進行計算太復雜，因此反求彈性參數也是很復雜的。2）AkiKI和RichardsPG的縱波近似公式1980年AkiKI和RichardsPG給出了一個縱波反射振幅的一個近似表達式：R（）=（5-3-1）這個公式是有適用條件的，就是界面兩邊的彈性介質性質的百分比變化小，式中的參數有：Vp=Vp2-Vp1、Vs=Vs2-Vs1、Vp=（Vp2+Vp1）/2、Vs=（Vs2+Vs1）/2、=2-1、=（2+1）/2、=（1+2）/2，1是入射角，2是按斯奈爾定理計算的透射角。這個式子雖然對Zoeppritz方程進行了化簡，但是用了Vs做參數而未用作參數。3)Bortfe

9、ld近似式1961年，Bortfeld提出近似式Rpp=（5-3-2）4）Hilterman近似式 （5-3-3）5）Shuey簡化公式1985年Shuey利用上面的式子進行了簡化，用代替Vs，修改得： R（）= （5-3-4）其中： 、 。垂直入射時的反射振幅，1、2分別為入射介質和透射介質的泊松比。其中，界面兩側泊松比的差是一個至關重要的因素，這就是振幅與炮檢距關系研究的物理基礎。Shuey近似式的特點就是其三項都有明確的物理意義：第一，垂直入射時，1=2=0，R（0）=R0，即R0是垂直入射時的反射振幅。第二，在中等入射的情況( 030°），有近似tansin，于是有： （5-

10、3-5）此時，反射振幅與A有關，前兩項起作用。這時的反射系數與介質的泊松比有密切關系，因此，利用此式更能突出油氣特征。第三，對于大角度入射情況，反射振幅與速度變化有關中的第三項起主要作用。6）Simth 和 Gidlow 加權疊加方法1987年Smith和Gidlow根據Gardner等人給出的水飽和巖石的密度與速度的4次方根成正比的假設，將Aki和Richards（1980）提出的Zoeppritz方程進行了修改化簡： （5-3-6）式中：、Vp、Vs分別為縱橫波速度，B、C為常數。因此，盡管表示反射系數的簡化式有很多，但其最終目的都是為了求得一個最簡表達式來表示反射系數隨入射角的變化。而上

11、、下介質的泊松比又對反射系數隨入射角的變化起重要作用，因此，需要在表達式中包含標志油氣特征的參數。Hilterman的簡化方程與Shuey的簡化方程滿足了以上要求，且適用于各種層狀模型的AVO模擬，因此得到了廣泛的應用。在應用中最常見的一種簡化公式為： （5-3-7）式中：R為反射系數、為入射角；P為近似為零偏移距下的縱波的反射振幅，也稱AVO截距，其大小決定于上下層之間的縱波波阻抗差異；G為縱波反射振幅隨入射角的變化梯度，也稱AVO斜率，其大小取決于泊松比的變化。這個方程不能很好處理入射角較大的情況，但是簡單易懂，至今仍然非常實用。對NMO（正常時差校正）道集中每一個時間采樣點作sin2（為

12、入射角）和振幅交匯圖（圖5-3-2），截距描述了正常入射角時的P波反射率（NIP），同時斜率就是梯度（振幅隨入射角怎么變化）。式（5-3-7）表明，在入射角小于中等角度（一般為30°）時，縱波反射系數近似與入射角正弦值的平方成線性關系。基于該方程的疊前反演可以獲得AVO屬性參數P與G及其各種轉換屬性參數，對這些屬性參數的解釋是：a、由AVO截距組成的P剖面是一個真正的法線入射零炮檢距剖面。b、由AVO斜率組成的G剖面反映的是巖層彈性參數的綜合特征。c、在橫波剖面上，當縱、橫波速度比近似等于2時，P-G可以反映出橫波波阻抗的特征。d、在泊松比剖面上，當縱、橫波速度比近似等于2時，P+G

13、反映的是泊松比的特征。e、截距與梯度的乘積（P×G）剖面，也稱AVO強度剖面，更有利于識別氣層。圖5-3-2：sin2（為入射角）和振幅交匯圖3、 AVO的地質意義AVO應用的基礎是泊松比的變化，而泊松比的變化是不同巖性和不同孔隙流體介質之間存在差異的客觀事實。大量的試驗研究和實踐表明，沉積巖的泊松比值具有如下特點：a、未固結的淺層鹽水飽和沉積巖往往具有非常高的泊松比值（0.4以上）；b、泊松比往往隨孔隙度的減小及沉積固結程度而減少；c、高孔隙度的鹽水飽和砂巖往往具有較高的泊松比值（0.30.4）；d、氣飽和高孔隙度砂巖往往具有極低的泊松比值（如低到0.1）。一般來講，不同巖性按泊松

14、比值從高到低排序依次是：石灰巖、白云巖、泥巖和砂巖，在砂巖中由于孔隙流體性質的差異，依次是：水砂巖、油砂巖和氣砂巖。從波速到泊松比，人們對地層的研究已逐漸深入到其本質。在研究振幅變化與波速的關系時，發現振幅變化與其他彈性參數的關系更有意義，這是因為波速不是一個獨立的彈性參數，而是介質的幾個彈性參數組合的結果。在拉梅常數（）、切變模量（）和泊松比（）等彈性參數中，泊松比在Zoeppritz方程中是以一個獨立的參數出現的。而且，通過試驗研究，發現泊松比是對區分巖性有特殊作用的一個參數。Gregory（1976）分別對各種巖性的泊松比做了大量測量試驗，結果表明不同巖性泊松比的差別比速度的差別大。因此

15、，利用泊松比判別巖性更可靠。所以，AVO技術的地質基礎在于不同巖石以及含有不同流體的同類巖石之間泊松比存在差別。Domenico（1977）研究了含氣、含油、含水砂巖的泊松比隨埋藏深度的變化規律，結果發現含不同流體砂巖的泊松比隨深度的變化特征是不同的：含氣砂巖的泊松比隨著深度的增加而增加，但泊松比的值總是小于含油和含水砂巖的泊松比值；含水砂巖的泊松比隨著深度的增加而減小，但泊松比的值總是大于含油和含氣砂巖的泊松比值；含油砂巖的泊松比也隨著深度的增加而減小，泊松比的值總是介于含水和含氣砂巖泊松比值之間。巖石物性研究發現，當砂巖中含氣時，縱波速度明顯降低，含氣層的泊松比較小，與圍巖的泊松比之差大都

16、在0.20.3之間。因此，一般都可檢測到較明顯的AVO響應。含油層由于泊松比值明顯大于含氣層，與圍巖的泊松比差值較小甚至接近，因此反射系數隨炮檢距的變化程度明顯小于含氣層，AVO響應要比含氣層弱得多，而且包括了所有可能出現的AVO響應，檢測也將更為困難。因此，AVO技術應用在尋找氣藏方面更為有利，更能體現其優越性。一般來講，不同巖性按泊松比值從高到低排序依次是：石灰巖、白云巖、泥巖和砂巖。4、 AVO應用及簡要流程簡單來說，AVO技術是通過研究地下介質的地震反射波振幅隨炮檢距的變化來反映地下介質的巖性和孔隙流體的性質，進而直接預測儲層。在實際應用中，就是利用地震反射的CDP道集資料，分析儲層界

17、面上的反射波振幅隨炮檢距的變化規律，或通過計算反射波振幅隨其入射角的變化參數，估算界面上的AVO屬性參數和泊松比差，進一步推斷儲層巖性和含油氣性質。通常采用泊松比參數來描述反射界面振幅的變化情況：當介質間無明顯泊松比變化時，不論反射系數是正負，振幅都隨入射角的增大而減小；當反射系數為正且泊松比增加（或反射系數為負而泊松比降低時），振幅隨入射角的增大而增加；當反射系數為負且泊松比降低（或如果反射系數為正而泊松比增加時），振幅隨入射角增大先減小，當入射角增大到一定時會出現極性反轉。因此，利用AVO技術中振幅隨入射角變化這一特征可判定巖石物理參數。設各炮檢距采樣的實際振幅為Ai，那么與模型曲線相比，

18、各振幅均方誤差為：（5-3-8）若要使e取最小值，只需e對和偏導數，并令其為0，用曲線擬合實際地震資料就可求得和，即：= （5-3-9）= （5-3-10）式（5-3-9）和（5-3-10）中i=1，2，n為道號、Ai為振幅、Wp和Ws為加權系數。通過這種方法，可以求得每一道上的權系數，權系數的大小變化可反映出振幅隨偏移距變化的情況，加權疊加的結果包含了振幅隨偏移距變化的信息。由曲線擬合可以得到P波速度反射率剖面、S波速度反射率剖面及衍生出的擬泊松比反射率剖面和流體因子剖面。發展出來的AVO油氣檢測技術是利用疊前地震資料中反射波振幅與炮檢距變化關系，研究地下巖性變化并進行油氣檢測和油氣富集帶圈

19、定的一項技術。由于AVO技術是根據振幅隨炮檢距的變化做反映的地下巖性以及孔隙流體的性質來直接檢測油氣和估計巖性參數的一項技術。加上Zoeppritz方程的復雜性和非直觀性，現在很多學者都在研究不同的簡化方法，得到了不同的近似公式，針對不同的公式就提出了不同的分析和反演方法。目前AVO技術較成功的實例大多在碎屑巖中，該技術在碳酸鹽巖中應用的實例不多，這主要是碳酸鹽巖本身的特點所決定的。在碳酸鹽巖中，決定于礦物成分，而非孔隙度；碳酸鹽巖是印膜孔隙或空穴孔隙，不是粒間孔隙，決定Vp的是基質而不是孔隙和充填流體；充氣后，Vp沒有明顯減低，泊松比變化不大，然而Willston（1984）對砂巖分析表明：

20、含氣對碳酸鹽巖巖性不僅有影響，而且影響很大。隨密度減小，P波速度和減小，S波速度輕度增加，另外，隨孔隙度的增加巖石對流體反應更靈敏。通常情況下，碳酸鹽巖振幅隨偏移距變化是速度、孔隙度及含氣性的綜合響應。5、 巖石物理建模塔中地區通過對中古5-8井區已鉆井進行巖石物理建模分析，以及查對相應巖石組分的實驗室試驗數據的綜合分析，得出以下認識：純泥巖，包括（夾泥灰巖）的巖石特征表現為低縱、橫波速度，低密度、低泊松比（在2.83.0）；致密灰巖表現為高縱、橫波速度，高密度；而對于儲層來說，不論是鉆井放空段儲層還是被泥質、螢石或鐵礦充填的儲層，其巖石特征都表現為縱波速度低，橫波速度略低，低密度和低泊松比（

21、2.93.1），可以看出儲層的泊松比與泥巖的泊松比有部分的交叉重合。為了更好的來區分致密灰巖和儲層，進行了彈性參數交匯。通過在交匯圖（圖5-3-3）上圈出特定區域，來影射到密度、縱、橫波等曲線上的不同巖性段來劃分儲層和非儲層。通過圈定劃分認為儲層主要分布于從尾段開始分叉的上部的條帶上，尾段的下部條帶分析認為以泥巖為主；從分叉部位開始至“掃帚”的把柄處為非儲層，圈定結果與上述巖石組分特征分析結果一致，且對應于鉆井主要油層段。同樣，通過縱波阻抗和泊松比的交匯也可以較好地區分儲層。在縱波阻抗與泊松比交匯圖(圖5-3-3)上，所有的數據呈現與縱橫波交匯類似的“掃帚”形狀，儲層同樣發育與“掃帚”的尾端分

22、叉部分。通過交匯圖上圈定的儲層部分投射工區實際鉆井中，結合該井產量證明了圈定范圍的合理性，說明建模方法和參數選取合理。圖5-3-3 圖5-3-3：ZG10、ZG5井區測井曲線交匯圖二. 彈性阻抗及參數反演技術疊前彈性參數反演以描述平面波反射和透射的Zoeppritz方程為基礎，利用反射系數隨入射角變化與地層彈性參數間的關系，采用一定的數學算法，從疊前地震數據中估算巖石的彈性參數（如縱波速度、橫波速度、泊松比、密度等），進而利用這些彈性參數進行巖性分析和含油氣的預測。目前常規方法一般可分為彈性阻抗反演（Elastic Inversion）和疊前同時反演（Simultaneous Inversio

23、n）兩類，分別簡敘如下：1、 彈性阻抗反演彈性阻抗是聲波阻抗的推廣，它是縱波速度、橫波速度、密度以及入射角的函數（式5-3-11）。可以簡單地表示為： （5-3-11）其中，E1、E2分別為地震波能量傳播時反射界面上、下兩層介質的彈性阻抗，是縱波速度Vp、橫波速度Vs、密度以及入射角的函數。R()是以角入射時的反射系數。在應用中，把經過疊前精細保幅處理和偏移的共反射點道集，分成多個入射角疊加，形成不同入射角剖面，將3個及以上不同入射角剖面分別進行子波提取和疊后(彈性)阻抗反演，得到不同角度下的彈性阻抗EI (1 )、EI (2 )、EI (3 )，并形成方程組，就可以聯合求出縱、橫波速度及密度了。其流程如圖5-3-10。可見，彈性阻抗反演能有效地解決地震子波隨炮檢距變化的問題，得到了不同入射角的波阻抗，即彈性阻抗；這是以另一種方式來表示AVO信息的方法。通過聯合求解巖石物理參數后，可以進一步開展多信息交匯進行巖性和含油氣性的綜合解釋。圖5-3-10：疊前彈性阻抗反演流程圖2、 疊前同時反演疊前同時反演就是利用不同角道集的地震數據、層位數據、測井數據進行同時反演，直接得到縱、橫波阻抗和密度。與上述疊前彈性阻抗反演比，理論上并沒有本質區別，僅是