劉仕友，陳殿遠(yuǎn)，周家雄

（中海石油（中國(guó)）有限公司湛江分公司，湛江 524057）

調(diào)諧AVO的存在條件及其識(shí)別

劉仕友，陳殿遠(yuǎn)，周家雄

（中海石油（中國(guó)）有限公司湛江分公司，湛江 524057）

自80年代以來(lái)，AVO屬性技術(shù)逐漸成為一種有效的儲(chǔ)層預(yù)測(cè)和烴類(lèi)檢測(cè)技術(shù)，但是在實(shí)際應(yīng)用中，AVO的復(fù)雜性、多解性也越來(lái)越為人們所熟知。疊前CRP道集，各偏移距傳播路徑是不同的，遠(yuǎn)近道頻率也是不同，頻率差異與地層耦合會(huì)產(chǎn)生一種調(diào)諧AVO現(xiàn)象。這里基于波動(dòng)的疊加和干涉原理，推導(dǎo)了薄層的振幅響應(yīng)公式，提出了調(diào)諧頻率的計(jì)算公式。通過(guò)正演模型分析了調(diào)諧AVO的特征：①當(dāng)?shù)貙拥恼{(diào)諧頻率高于地震資料頻率時(shí)，近道調(diào)諧作用大于遠(yuǎn)道，表現(xiàn)為IV類(lèi)異常；②當(dāng)?shù)貙拥恼{(diào)諧頻率低于地震資料頻率時(shí)，近道調(diào)諧作用小于遠(yuǎn)道，表現(xiàn)為III類(lèi)異常；③當(dāng)?shù)貙拥恼{(diào)諧頻率位于地震資料頻率區(qū)間時(shí)，AVO特征表現(xiàn)為“V”特征。調(diào)諧AVO存在的基礎(chǔ)為薄層加變化的子波頻率，最后結(jié)合樂(lè)東X1目標(biāo)實(shí)際例子，探索了調(diào)諧AVO的識(shí)別方法。

調(diào)諧效應(yīng)；AVO分析；薄互層；調(diào)諧頻率；調(diào)諧厚度

0 引言

隨著地震采集與處理技術(shù)的進(jìn)一步發(fā)展，地震資料品質(zhì)越來(lái)越高，進(jìn)而地震勘探成果的精度得到了大幅度提高?；诏B前道集資料的AVO信息，不僅包含了縱波速度、密度信息，還包含了橫波速度信息，可以有效地反應(yīng)了地層的巖性、物性以及孔隙流體性質(zhì)。通過(guò)AVO分析可以有效地降低地震勘探多解性、極大地提高了油藏流體識(shí)別精度，因此在油氣勘探中得到了廣泛重視。

自1982年Ostrander提出AVO技術(shù)以來(lái)，眾多學(xué)者對(duì)地震AVO技術(shù)進(jìn)行了大量的研究［1，2］，包括AVO參數(shù)與地層信息的關(guān)系，Lortzer等［3］討論了根據(jù)疊前CMP道集上的AVO信息進(jìn)行求取地震參數(shù)差異的方法；Sung等［4］研究了油藏壓力和含油飽和度變化對(duì)時(shí)移地震AVO的影響；劉洋［5］和陰可等［6］對(duì)各向異性介質(zhì)的AVO進(jìn)行了研究；趙偉等［7］研究了薄互層調(diào)諧效應(yīng)對(duì)AVO的影響。

AVO研究的理論基礎(chǔ)是Zoeppritz方程，該方程僅對(duì)單一界面地震振幅隨人射角的變化進(jìn)行了精確描述，但沒(méi)有考慮儲(chǔ)層厚度和地震波頻率對(duì)AVO的影響。地震波在傳播過(guò)程中，地震波的頻率成分會(huì)發(fā)生變化，在共反射點(diǎn)道集上，遠(yuǎn)道與近道的頻譜差異明顯，一般表現(xiàn)為規(guī)律性變化，如圖1所示。隨著埋深的增加這種差異也越來(lái)越大，中深層地震資料的近道頻譜與遠(yuǎn)道頻譜差異更為明顯，有時(shí)主頻差異甚至可達(dá)15Hz～20Hz。

當(dāng)遠(yuǎn)近道頻譜差異較小時(shí)，其影響可以忽略不計(jì)，但是對(duì)于遠(yuǎn)近道頻譜差異較大時(shí)，不同頻率的調(diào)諧差異將會(huì)顯現(xiàn)起來(lái)，并在CRP道集上表現(xiàn)出來(lái)，作者將從波動(dòng)的疊加和干涉原理出發(fā)，對(duì)地層厚度和頻率共同作用的調(diào)諧AVO做出論證分析。

圖1 地震資料主頻與偏移距關(guān)系示意圖Fig.1 The relation cartogram between frequency and offset

1 基本原理

基于地震勘探原理中關(guān)于薄層的定義［8］，一般用地震波波長(zhǎng)λ／2作為薄層的界限，即厚度遠(yuǎn)大于λ／2的稱(chēng)為厚層，厚度相當(dāng)于或小于λ／2的屬薄層。當(dāng)?shù)貙铀俣葹? 600m／s，地震主頻為25Hz時(shí)，λ／2將是36m，對(duì)于某些河流相等儲(chǔ)層來(lái)講，存在著大量的所謂薄層。

目前討論調(diào)諧影響采用的最簡(jiǎn)便而實(shí)用的方法［7］是從波動(dòng)的疊加和干涉原理出發(fā)，對(duì)上界面的反射波和下界面的反射波進(jìn)行疊加，利用疊加后波的特點(diǎn)，來(lái)說(shuō)明薄層的存在對(duì)反射波的影響。

圖2 薄層模型Fig.2 The model of thin bed

假設(shè)有一個(gè)薄層夾于均勻?qū)又校▓D2），并認(rèn)為薄層頂、底的反射波相同，振幅一樣，只存在反射時(shí)差。再把子波簡(jiǎn)化為周期T，頻率為f，振幅A的余弦簡(jiǎn)諧波。

假設(shè)設(shè)上界面反射波為R1：

對(duì)應(yīng)的下界面反射波為R2：

其中：V是薄層的縱波速度；b為薄層厚度。上、下界面反射的時(shí)差為，由于式（1）～式（2）中，兩個(gè)波

對(duì)上、下界面反射波進(jìn)行疊加，則合成的薄層反射為：

所以薄層反射振幅的絕對(duì)值是：

經(jīng)過(guò)簡(jiǎn)單的假定推導(dǎo)，我們可以獲得薄層反射振幅Ad與薄層厚度b、頂?shù)讍蝹€(gè)反射振幅A、波長(zhǎng)λ、頻率f之間的簡(jiǎn)單明確的關(guān)系，如式（6）所示。對(duì)薄層而言，振幅是反射系數(shù)、厚度、頻率、地層速度的綜合表征，利用式（6）可以很容易的描述出薄層反射振幅Ad與當(dāng)?shù)貙雍穸然蛘叩卣鹱硬l率中的關(guān)系，其變化特征如圖3所示。

假定地層速度不受頻率影響，可以發(fā)現(xiàn)，當(dāng)頻率一定時(shí)，當(dāng)層厚為λ／4時(shí)，調(diào)諧作用最大，合成振幅為最大，厚度λ／4稱(chēng)之為調(diào)諧厚度；當(dāng)厚度一定時(shí)，子波頻率為V／4b時(shí)，調(diào)諧作用最大，合成振幅為最大，頻率V／4b稱(chēng)之為調(diào)諧頻率。

對(duì)于疊前CRP道集而言，對(duì)于每一個(gè)偏移距可以利用同樣的思路來(lái)進(jìn)行描述，不同的偏移距看成一個(gè)疊后記錄，但其反射系數(shù)有所差異；一般而言，地層厚度不大的情況下，上、下界面的反射角幾乎沒(méi)有什么差異，因此忽略入射角的影響；因此，疊前CRP也基本滿(mǎn)足前文薄層反射的假設(shè)，同樣可以利用公式（6）進(jìn)行調(diào)諧對(duì)AVO影響理論的討論。具體計(jì)算步驟為：

圖3 振幅與頻率、厚度的關(guān)系示意圖Fig.3 The relation cartogram between amplitude，frequency and thickness

1）根據(jù)地層參數(shù)和偏移距，計(jì)算出對(duì)應(yīng)的入射角，進(jìn)而求出上下地層的反射振幅A；

2）每個(gè)偏移距的子波頻率、厚度已知，利用公式（6），可求取對(duì)應(yīng)的調(diào)諧振幅Ad－r；

3）根據(jù)求得的每個(gè)厚度的不同偏移距對(duì)應(yīng)的Ad－r可獲得單點(diǎn)的振幅變化規(guī)律。

為了進(jìn)一步研究調(diào)諧作用對(duì)疊前資料的AVO分析的影響，在不考慮地層衰減吸收的情況下，直接假定遠(yuǎn)近道的地震子波頻率，即子波頻率隨著偏移距逐漸變化，其中近道子波選取主頻為50Hz的，向遠(yuǎn)道逐漸變化為主頻為35Hz的子波。為了測(cè)試調(diào)諧的AVO假象，分別選擇的砂泥二相地層參數(shù)為：①泥巖的參數(shù)：縱波速度3 000m／s，橫波速度2 120m／s，密度2.6g／cc；②砂巖參數(shù)：縱波速度2 000m／s，橫波速度1 415m／s，密度1.9g／cc；楔狀模型如圖4（a）所示。利用AKi公式，可以求出P、G屬性，其中G屬性很小，即此砂泥界面的AVO變化很微弱，近似可以認(rèn)為沒(méi)有AVO異常。采用公式（6），根據(jù)不同偏移距的不同主頻子波，在每一個(gè)厚度點(diǎn)處可以計(jì)算出不同偏移距的合成記錄道，模擬結(jié)果如圖4所示。

圖4 調(diào)諧AVO與地層厚度間的關(guān)系Fig.4 The relationship between tuning AVO and thickness

模擬結(jié)果顯示，當(dāng)?shù)貙幼陨鞟VO微弱時(shí)，調(diào)諧產(chǎn)生的假象讓地層產(chǎn)生類(lèi)似AVO異常，并表現(xiàn)出了隨著厚度變化呈現(xiàn)規(guī)律性變化。

當(dāng)?shù)貙幼畋。ê穸刃∮?0m）時(shí)，近道的調(diào)諧要強(qiáng)于遠(yuǎn)道的調(diào)諧，這時(shí)近道能量要強(qiáng)于遠(yuǎn)道，表現(xiàn)為振幅隨著偏移距的增大而減小的特征，類(lèi)似于四類(lèi)AVO異常。

當(dāng)儲(chǔ)層厚度較大（大于14.3m小于半波長(zhǎng)時(shí)）時(shí)，遠(yuǎn)道的調(diào)諧要強(qiáng)于近道，表現(xiàn)為類(lèi)似Ⅲ類(lèi)AVO異常，振幅表現(xiàn)為振幅隨著偏移距的增大而增大的特點(diǎn)。

當(dāng)儲(chǔ)層厚度位于10m～14.3m間時(shí)，這時(shí)地層的調(diào)諧頻率恰好位于某個(gè)偏移距位置上，這個(gè)偏移距的調(diào)諧振幅最大，而左右兩端都要變?。徽麄€(gè)振幅特征則表現(xiàn)為振幅隨偏移距的增大，先增大后減小的特點(diǎn)。模型正演模擬結(jié)果揭示：隨著儲(chǔ)層厚度變大，則“V”頂點(diǎn)往遠(yuǎn)偏移距位置漂移；反之，則往近偏移距位置漂移（圖5）。

圖5 調(diào)諧AVO的V頂點(diǎn)與地層厚度間的關(guān)系Fig.5 The relation between the peak of tuning－AVO and thickness

2 實(shí)例應(yīng)用

樂(lè)東X－1目標(biāo)在地震上表現(xiàn)為強(qiáng)振幅的亮點(diǎn)特征，高部位表現(xiàn)為一紅一黑的波形組合特征（圖6），低部位振幅明顯增強(qiáng)，表現(xiàn)為兩紅一黑的波組特征。同步彈性阻抗反演顯示，目標(biāo)體表現(xiàn)為低縱波阻抗、低橫波阻抗、低密度的特征、低縱橫波速度比的“四低”特征。尤其是低部位，“四低”現(xiàn)象更為明顯。低部位為III類(lèi)AVO異常，高部位為弱Ⅳ類(lèi)異常。

圖6 疊后純波地震剖面Fig.6 The seismic profile of pure wave

通過(guò)對(duì)地震資料頻譜分析發(fā)現(xiàn)，目標(biāo)層段近道地震資料（offset 400m～1 200m）的主頻約為35 Hz，遠(yuǎn)道地震資料（offset 4 000m～5 200m）的主頻約為25Hz，，遠(yuǎn)近道頻率差異明顯，可能存在調(diào)諧AVO的影響。

逐點(diǎn)AVO屬性分析如圖7所示，分析發(fā)現(xiàn)目標(biāo)體的疊前AVO變化有以下3個(gè)特點(diǎn)：

1）構(gòu)造低部位表現(xiàn)為振幅隨著偏移距增大而增大的III類(lèi)AVO異常，構(gòu)造的高部位表現(xiàn)為振幅隨著偏移距增大而減小的IV類(lèi)AVO異常。

圖7 逐點(diǎn)AVO屬性分析圖Fig.7 The AVO analysis by point－to－point

2）在AVO異常由III類(lèi)逐漸變?yōu)镮V類(lèi)間，目標(biāo)體表現(xiàn)為振幅隨著偏移距的增大，先增大后減?。▓D8）的趨勢(shì)，即成一個(gè)“V”形分布。

3）通過(guò)多點(diǎn)分析發(fā)現(xiàn)，疊前道集的“V”頂點(diǎn)成一種規(guī)律形的漂移，結(jié)合波阻抗反演等認(rèn)識(shí)，目標(biāo)體為一套厚度逐漸發(fā)生變化的砂體（圖9），低部位為一套厚度逐漸減薄的特征，結(jié)合頻率調(diào)諧的分析認(rèn)識(shí)，薄層的厚度變大，它的調(diào)諧頻率則變小，即對(duì)應(yīng)調(diào)諧位置的偏移距更大，通過(guò)分析不同厚度的“V”頂點(diǎn)位置發(fā)現(xiàn)，變化規(guī)律符合理論模型認(rèn)識(shí)。

圖8 單點(diǎn)AVO屬性分析圖Fig.8 The AVO analysis by single point

圖9 “V”形調(diào)諧點(diǎn)與地層厚度對(duì)比圖Fig.9 The comparison diagram between tuning point（V）and thickness

綜合以上分析，目標(biāo)體存在產(chǎn)生調(diào)諧AVO的條件，其空間變化特征也符合調(diào)諧AVO的規(guī)律，認(rèn)為目標(biāo)體低部位的III類(lèi)AVO異常主要為調(diào)諧影響造成的假象，高部位厚層的AVO特征是地層AVO的真實(shí)響應(yīng)，即目標(biāo)體表現(xiàn)振幅隨著偏移距增大而減小的IV類(lèi)AVO異常。

3 結(jié)論

1）當(dāng)頻率一定時(shí)，對(duì)于地層厚度，當(dāng)層厚為λ／4時(shí)，調(diào)諧作用最大，合成振幅為最大；當(dāng)厚度一定時(shí)，子波頻率為V／4b時(shí)，調(diào)諧作用最大，合成振幅為最大。

2）調(diào)諧AVO存在的前提條件是薄層（地層厚度小于半波長(zhǎng)）和變化的地震頻率（遠(yuǎn)近道頻譜有明顯的差異）同時(shí)存在。

3）當(dāng)?shù)貙拥恼{(diào)諧頻率高于地震資料頻率時(shí)，近道調(diào)諧作用大于遠(yuǎn)道，表現(xiàn)為偽IV類(lèi)異常；當(dāng)?shù)貙拥恼{(diào)諧頻率低于地震資料頻率時(shí)，近道調(diào)諧作用小于遠(yuǎn)道，表現(xiàn)為偽III類(lèi)異常；當(dāng)?shù)貙拥恼{(diào)諧頻率位于地震資料頻率區(qū)間時(shí)，AVO特征表現(xiàn)為“V”特征。

4）薄層AVO是調(diào)諧AVO與地層真實(shí)AVO共同作用的結(jié)果，地層真實(shí)AVO的識(shí)別利用需要空間上特征的分析，或者其變化規(guī)律，單點(diǎn)分析的結(jié)果大多都是復(fù)合AVO效應(yīng)。

［1］ OSTRANDER W J.Plane－wave reflection coefficients for gas sands at non－normal angles of incidence ［J］.Geo－physics，1984，49（10）：1637－1648.

［2］ 鄭曉東.AVO理論和方法和一些新進(jìn)展［J］.石油地球物理勘探，1992，27（3）：305－317.

ZHENG X D.Some new development of AVO technique［J］.OGP，1992，27（3）：305－317.（In Chinese）

［2］ LORTZER G J M，HAAS J C，BERKHOUT A J.E－valuation of weighted stacking techniques for AVO inversion［C］.Expanded Abstracts of 58＇"Annual Internet SEG Mtg，1988：1204－1208.

［3］ SUNG H Y，SEONGSIK Y.Effects of pressure and fluid saturation changes on time－lapse AVO response［C］.Expanded Abstracts of 70th Annual Internet SEG Mtg，2000：1481－1484.

［4］ 劉洋，董敏煌.各向異性介質(zhì)中的方位AVO［J］.石油地球物理勘探，1999，34（3）：260－268.

LIU Y，DONG M H.Azimuthal AVO in anisotropic medium［J］.OGP，1999，34（3）：260－268.（In Chinese）

［5］ 陰可，楊慧珠.各向異性介質(zhì)中的AVO［J］.地球物理學(xué)報(bào)，1998，41（3）：382－391.

YIN K，YANG H ZH.AVO IN ANISOTROPIC MEDIA［J］.Chinese Journal Geophysics，1998，41 （03）：382－391.（In Chinese）

［6］ 趙偉，陳小宏，李景葉.薄互層調(diào)諧效應(yīng)對(duì)AVO的影響［J］.石油物探，2006，45（6）：570－573.

ZHAO W，CHEN X H，LI J Y.The effect of thin interbed tuning on AVO［J］.Geophysical Prospecting for Petroleum，2006，45（6）：570－573.（In Chinese）

［7］ 陸基孟，地震勘探原理（下冊(cè)）［M］.青島：石油大學(xué)出版社，1993.

LU J M，The principle of seismic exploration［M］.Qindao：China University of Petroleum Press，1993.（In Chinese）

The existence conditions and identification of tuning AVO

LIU Shi－you，CHEN Dian－yuan，ZHOU Jia－xiong
（CNOOC Ltd.，Zhanjiang，Zhanjiang 524057，China）

Since the 1980s，AVO attribute technology becomes an effective technology for reservoir prediction and hydrocarbon detection，but in practical application，the complexity and multiplicity of the AVO are more and more popular.In Prestack CRP gathers，different offsets have different propagation paths，and gathers of different distances have different frequencies.Frequency difference and stratigraphic coupling can produce a kind of tuning AVO phenomenon.Based on superposition principle and interference principle of fluctuation，formula of the thin layer amplitude response is derived，and the formula of tuning frequency is put forward.The characteristics of tuning AVO is analyzed by forward model：when the tuning frequency of the layer is higher than the frequency of seismic data，near tuning effect had more effects and it shows the AVO of class IV；when the tuning frequency of the layer is lower than the frequency of seismic data，far tuning effect had more effects and it shows the AVO of class III；when the tuning frequency of the layer is in the frequency range of seismic data，AVO feature is the"V"characteristic feature.The basis of tuning AVO is thin layer and changed wavelet frequency.In the end of the paper，the identification method of tuning AVO is explored by the evaluation of Ledong X1target region.

tuning effect；AVO；thin bed；tuning frequency；tuning thickness

P 631.4

：A

10.3969／j.issn.1001－1749.2015.06.16

1001－1749（2015）06－0768－05

2014－09－22改回日期：2014－10－08

國(guó)家科技重大專(zhuān)項(xiàng)課題（2011ZX05023－004）

劉仕友（1982－），男，工程師，主要從事儲(chǔ)層預(yù)測(cè)及烴類(lèi)檢測(cè)工作研究，E－mail：liushiyou＠139.com。