劉秀娟，劉金朋，方中于，梁立鋒 ,唐永杰

(1.中海油能源發(fā)展工程技術(shù)物探技術(shù)研究所，廣東 湛江 524000；2.中國地質(zhì)大學 地球物理與空間信息學院，湖北 武漢 430074)

?

薄砂體AVO曲線計算及特征分析

劉秀娟1，劉金朋1，方中于1，梁立鋒1,唐永杰2

(1.中海油能源發(fā)展工程技術(shù)物探技術(shù)研究所，廣東 湛江 524000；2.中國地質(zhì)大學 地球物理與空間信息學院，湖北 武漢 430074)

常規(guī)薄砂體AVO正演模擬一般采用Zoeppritz方程或者簡化公式進行模擬，由于未考慮薄砂體厚度對AVO的影響，因此方法本身具有局限性。鑒于此，采用了基于薄層的反射率法計算RVA及AVO曲線，并討論了反射率法精度較高的機理。通過對惠州某區(qū)薄砂體理論模型的RVA計算，反射率法可以精細計算不同薄層厚度及不同頻率的RVA曲線，提高了RVA曲線描述薄層的精度，而利用Zoeppritz方程計算的RVA曲線對薄砂體的厚度及入射波頻率不敏感；對比地震道集AVO、基于Zoeppritz方程模擬的AVO及反射率法模擬的AVO結(jié)果表明，反射率法模擬的AVO與地震道集AVO對應更好，說明了本文方法的有效性。

薄砂體；AVO；正演；RVA

1 引 言

目前針對薄儲層的地震模擬方法更多的采用基于Zoeppritz方程進行模擬，從方法本身來講，快捷、方便，易于使用，很多商業(yè)化軟件中也集成了這一模塊或計算功能。通過應用該方法模擬薄儲層模型，可以從理論上解決薄儲層的厚度計算問題，對地震資料的解釋也有一定的指導作用[1]。為了進一步使應用更加方便，鄭曉東[2]等人對該方程進行了簡化，從而使得由解矩陣方程組改為解線性方程的形式，大大簡化了計算過程，且具有明確的物理意義，推動了該方程在地震正演和反演方面的應用。然而，實際上地下薄儲層波場特征非常復雜，使得Zoeppritz方程的前提條件很難得到滿足，即假設(shè)地層分界面上下兩側(cè)為半無限彈性介質(zhì)分界面的假設(shè)條件，顧漢明、蔡希玲、杜勁松、袁子龍等眾多學者都進行了深入的研究[3-6]。在AVO(Amplitude Versus Offset)正反演演應用上，夏媛媛等[7]應用正演模擬技術(shù)指導地震資料波阻抗反演解釋，有效地減少地震解釋的多解性；吳國慶等[8]設(shè)計了賦存不同流體的薄煤層模型，通過地震AVO正演模擬技術(shù)，分析了含不同類別流體儲層的反射系數(shù)曲線隨入射角度變化特征及疊前道集的振幅特性；姜勇等[9]研究了AVO流體反演技術(shù)，并應用在低孔滲砂巖氣藏預測；尹川等[10]對AVO流體的指示因子敏感性進行了研究；梁立鋒等[11,12]研究了突出薄層AVO特性的反射率法，并將其應用在惠州某區(qū)的薄儲層地震響應模擬中。在前人研究的基礎(chǔ)上，本文探討了薄層情況下的RVA(Reflectivity Versus Angle)及AVO計算方法，設(shè)計了抽象薄層地質(zhì)模型，并對模型進行了試算，定量分析了不同頻率的薄砂體RVA特征；對比了惠州某區(qū)地震道集AVO、基于Zoeppritz方程模擬的AVO及反射率法模擬的AVO三者的差異；通過薄層的調(diào)諧效應，驗證了本文方法的正確性。

2 方法原理

2.1 RVA曲線計算方法

針對海上勘探階段少井條件下的薄儲層復雜波場特征，采用逐層遞推的方法獲得薄層系的應力和位移，進而獲得頻率域縱波反射系數(shù)精確解析公式。假設(shè)n-1個薄層組成的彈性層系，它的上面為液體半空間，下面為固體半空間(圖1)。平面聲波從液體半空間入射到層系，入射波的振幅為1，入射角為θn+1。在彈性層系的每個彈性層中，將存在一對縱波和一對橫波，它們分別相對于水平對稱向上和向下傳播。在下半空間中，縱波和橫波向下傳播而離開層系，此模型適合海上地震采集系統(tǒng)。此時，存在關(guān)系式：

ξ=kn+1sinθn+1=kisinθi=Kisinγi;i=1,2,…,n

(1)

圖1 聲波在彈性層系上的反射與透射Fig.1 Acoustic reflection and transmission on the elastic media

其中，ki=ω/ci為縱波波數(shù)，單位為rad·s/m；Ki=ω/bi為橫波波數(shù)，單位是rad·s/m；ci為縱波速度，單位是m/s；bi為橫波速度，單位為m/s；θi為i層縱波入射角，單位為rad；γi為i層橫波入射角，單位為rad。

依據(jù)彈性層系第n+ 1層與第1層的位移與應力的遞推關(guān)系，可以獲得某一頻率下的不同入射角度對應的縱波反射系數(shù)精確解RVA(Reflec-tivity Versus Angle)曲線，縱波反射系數(shù)精確解表達式如方程式(2)所示：

(2)

其中矩陣B中的各元素可以通過文獻[11]求得；式中ω為入射波頻率；bn+1為第n+1層的橫波速度；i為虛數(shù)單位；ρn+1為第n+1層的密度；Kn+1為第n+1層的橫波波數(shù)。求解方程式(2)，即可得到平面諧和波在彈性層系上的縱波反射系數(shù)V，同時可得到橫波反射系數(shù)V′，橫波透射系數(shù)W′和縱波透射系數(shù)W。方程中把多個薄層作為一個系統(tǒng)研究，考慮了單個薄層對反射系數(shù)的貢獻，同時考慮了薄層厚度、入射波頻率對反射系數(shù)的影響，相對于Zoeppritz方程和其他簡化形式更具有優(yōu)越性。

2.2 AVO曲線計算

為了計算薄砂體(薄層系)AVO曲線，首先建立薄層地質(zhì)概念模型，并依據(jù)鉆井及測井資料，對薄砂體賦予彈性參數(shù)，如縱波速度、橫波速度、密度等，依據(jù)方程式(2)計算頻率域反射系數(shù)，可得到頻率域中頂層反射系數(shù)譜R(f,θ)，對R(f,θ)進行反傅里葉變換得到時間域的反射系數(shù)R(τ,θ)，若要求取某一界面的反射系數(shù)，可根據(jù)斯奈爾定律及縱向走時，在R(τ,θ)中尋找。另外，根據(jù)反射系數(shù)與子波褶積得到的理論地震記錄，在該理論地震記錄中，通過走時尋找得到對應的振幅值，為了得到不同偏移距或不同入射角度的振幅，重復該過程，即可得到薄砂體頂界面的AVO曲線。

該方法的優(yōu)點：①該方法模擬的地震記錄中包含反射波、轉(zhuǎn)換波、多次波信息，符合薄儲層地震波傳播規(guī)律，優(yōu)于目前國際上以Zoeppritz方程為基礎(chǔ)的波場正演模擬手段；②可以模擬任意多個薄層的地震波場，適用于多個砂體疊置情況下的AVO曲線分析，具有較強的通用性；③該方法不受頻散及網(wǎng)格大小影響，對薄層及砂體邊界刻劃較好。

3 薄砂體AVO曲線特征分析

3.1 RVA曲線特征計算

惠州某區(qū)薄儲層發(fā)育，儲層厚度從幾米到十幾米不等，AVO特征復雜，經(jīng)常出現(xiàn)地震道集AVO曲線與測井正演的AVO曲線不符合的情況，為地震勘探帶來了很大的困擾。為此，利用本文研究的薄層RVA及薄層AVO計算方法，進行薄層研究。根據(jù)研究區(qū)的地震地質(zhì)特征，設(shè)計了如圖2所示的抽象薄層地質(zhì)模型，其中彈性參數(shù)(縱波速度、橫波速度、密度)的取值參考了HZ21-A井的K08層。

分別利用Zoeppritz方程和本文的反射率法計算頂界面的反射系數(shù)譜RVA曲線(圖3)，其中頻率變化、儲層厚度不變。從圖3看到，兩者計算的RVA曲線在趨勢上是一致的；利用發(fā)射率法計算的RVA曲線，其不同頻率成分對應的結(jié)果是不同的，它隨著頻率的變化而變化。但是，利用Zoeppritz方程計算RVA曲線， 各頻率分量的能量值相同，不能直接分析厚度對各頻率分量的影響效果。

圖2 薄層地質(zhì)模型Fig.2 Model of thin layer geology

圖3 不同單頻率時RVA曲線Fig.3 RVA cure for different single frequencies

當薄層砂體厚度變化時，厚度從波長/4變化為波長/8、波長/16、波長/32，此時入射波頻率為50 Hz，同時用Zoeppritz方程和反射率法計算RVA曲線，圖4為計算結(jié)果， 其中橫坐標為入射角度，縱坐標為反射系數(shù)。從圖4中可以看出，Zoeppritz方程求取的反射系數(shù)僅僅與界面兩側(cè)的物性參數(shù)有關(guān)，與薄砂體的厚度并沒有關(guān)系，說明Zoeppritz方程對薄砂體的厚度不敏感，這也是利用測井資料正演的AVO與地震道集的不符合的一個重要原因。圖4中其余的幾條曲線為利用反射率法計算的RVA曲線，并且用箭頭指出了薄砂體的厚度，可以看到，薄砂體的厚度變化，帶來了RVA曲線的變化，并且隨著厚度減小，中間薄砂體對反射系數(shù)影響變小，底層的影響增加，反射系數(shù)趨近于頂?shù)捉缑鏄?gòu)成的反射系數(shù)。

圖4 不同層厚時的RVA曲線Fig.4 RVA curve for different levels of thickness

3.2 單個薄砂體的AVO曲線特征

惠州某區(qū)地震道集AVO特征與測井曲線正演的AVO曲線對應較差，為了探討這一問題的原因，對HZ-A井K22儲層進行正演研究。圖5(a)為過HZ-A井的偏移后的地震道集，圖中紅線所指的位置為K22儲層的頂界面，對應的雙程旅行時為2 072 ms；圖5(b)為從地震道集上拾取的K22儲層的AVO曲線，屬于第三類AVO反射特征；圖5(c)為利用HZ-A井縱橫波速度及密度曲線，應用Zoeppritz方程，進行AVO曲線模擬的結(jié)果，也屬于第三類AVO反射特征。對比圖5(b)、圖5(c)，可以看出，兩者大致形態(tài)相同，但是當入射角為21°時，正演的AVO曲線(圖5(c))出現(xiàn)拐點，但是實際道集的AVO曲線沒有出現(xiàn)拐點。圖5(d)為基于反射率法模擬的AVO曲線，應用的也是HZ-A井縱橫波速度及密度曲線，對比圖5(b)、圖5(c)、圖5(d)三個結(jié)果，可以看出，反射率法模擬的AVO曲線與實際地震道集的AVO曲線最相似，究其原因，利用反射率法求取的AVO包含了薄層效應，下層介質(zhì)的物性參數(shù)如厚度等都對其有影響，而Zoeppritz方程僅僅受單界面兩側(cè)物性參數(shù)影響。這也驗證了反射率法計算的AVO曲線的正確性及與地震的可對比性。

圖5 地震道集及三條AVO曲線對比Fig.5 Contrast of seismic gathers and the three AVO curves

3.3 方法驗證

為了驗證方法的正確性，對薄層地質(zhì)模型(圖2)進行計算。設(shè)計薄砂體厚度為18 m，計算過程中保持厚度不變，假設(shè)波入射頻率由10 Hz到90 Hz，間隔1 Hz遞變，應用反射率法計算RVA曲線，如圖6所示，圖中兩條曲線分別為0°入射的RVA曲線及30°入射的RVA曲線，從圖中可以看出，對于同一角度，隨著入射波頻率增大，RVA曲線都是先增大后減小。以0°入射角為例，當K08層厚為18 m，層速度按照3 550 m/s計算，則薄砂體厚度為調(diào)諧厚度(λ/4)時，理論上對應的調(diào)諧頻率為49.3 Hz。從圖6中可看到，對于0°入射，在50 Hz附近，振幅達到最大值，驗證了薄層的調(diào)諧效應，從而驗證了方法的正確性。

圖6 某一入射角度的反射系數(shù)隨頻率變化Fig.6 Reflection coefficient variation with frequency of incident wave

4 結(jié) 論

1)通過對惠州某區(qū)薄砂體理論模型的RVA計算，結(jié)果表明，利用Zoeppritz方程計算的RVA曲線對薄砂體的厚度及入射波頻率不敏感，而反射率法可以精細計算不同薄層厚度及不同頻率的RVA曲線，提高了RVA曲線描述薄層的精度，對于后續(xù)的薄層石油勘探及開發(fā)有重要的作用；

2)對比了惠州某區(qū)地震道集AVO、基于Zoeppritz方程模擬的AVO、反射率法模擬的AVO，三者對比結(jié)果表明，反射率法模擬的AVO與地震道集AVO對應更好，說明了本文方法的有效性；

3)分析了反射率法模擬的薄層AVO更精確的原因：從薄層地震反射機理上，利用反射率法求取的AVO包含了薄層效應，包含了層間多次波、轉(zhuǎn)換波對薄砂體振幅的影響，而Zoeppritz方程僅僅受單界面兩側(cè)物性參數(shù)影響。通過薄層的調(diào)諧效應，驗證了本文方法的正確性。

[1]王鵬飛，劉桂林，李嵩陽，等.春光區(qū)塊薄儲層地震響應特征的模型正演分析[J].石油地質(zhì)與工程，2012，26(6)：64-67.

[2]鄭曉東.Zoeppritz方程的近似及其應用[J].石油地球物理勘探,1991,26(2):129-145.

[3]顧漢明,王緯,陳國俊.復雜介質(zhì)中地震多次反射波快速正演模擬[J].地球科學,2001,26(5):541-544.

[4]蔡希玲,閆忠,崔全章，等.砂泥巖薄互層分辨率的理論分析[J].石油物探,2004,43(3):229-233.

[5]杜勁松,韓福民,陳發(fā)亮，等.薄互層正演模擬分析及其應用[J].斷塊油氣田,2004,11(6):4-7.

[6]袁子龍,楊冰,王建國.薄層、薄互層地震反射時間域與頻率域正演模擬研究及應用[J].石油物探,1996,35(3):14-20.

[7]夏媛媛，趙民，藏歌，等.正演模擬技術(shù)在解釋反演中的應用[J].工程地球物理學報，2014,11(6)：842-846.

[8]吳國慶，李剛，董守華，等.薄煤層賦存不同流體的AVO特征分析[J].工程地球物理學報，2015,12(3)：395-400.

[9]姜勇，李鍵，宋傳會.AVO流體反演技術(shù)在低孔滲砂巖氣藏預測中的應用[J].工程地球物理學報，2014,11(5)：673-677.

[10]尹川，顧漢明.AVO流體的指示因子敏感性分析[J]．工程地球物理學報，2008,5(1)：85-88.

[11]梁立鋒，劉秀娟.彈性層系反射系數(shù)正演[J].物探與化探,2007,28(1):70-74.

[12]梁立鋒，張宏兵，唐永杰，等.反射率法波場正演模擬及其在惠州某區(qū)薄儲層地震響應模擬中的應用[J].地質(zhì)科技情報，2016,35(3)：202-206.

The Calculation of Thin Sand Body AVO Curve and Its Feature Analysis

Liu Xiujuan1,Liu Jinpeng1,Fang Zhongyu1,Liang Lifeng1,Tang Yongjie2

(1.ResearchInstituteofCNOOCEnergyDevelopmentEngineeringandGeophysicalTechnology,ZhanjiangGuangdong524000,China;2.InstituteofGeophysics&Geomatics,ChinaUniversityofGeosciences,WuhanHubei430074,China)

Zoeppritz equation or simplified formula is used in conventional AVO forward simulation of thin sand body. It has some limitations in these methods due to not considering the influence of thickness of thin sand to AVO. In view of this,this paper adopts the reflectivity method based on thin layer to compute RVA and AVO curve,and discusses the reason why the reflectivity method has high precision. Through calculating RVA curve of thin sand model of Huizhou,reflectivity method can calculate accurately RVA curve for different thickness of thin layer and different frequencies. So it improves the accuracy of the RVA curve description to thin layer,but the RVA curve is not sensitive to thickness of thin sand and frequency of incident wave by using Zoeppritz equation method. We compare AVO among seismic angle gather and Zoeppritz equation simulation and reflectivity method simulation and the result of comparision shows that AVO by the reflectivity method is better in line with seismic angle gather,so the method presented in this paper is effective.

thin sand body; AVO; forward modeling; RVA

1672—7940(2016)06—0689—05

10.3969/j.issn.1672-7940.2016.06.001

中國海洋石油總公司科技項目(編號：CNOOC-KJ 125 ZDXM 07 LTD NFGC 2014-04)；國家自然科學基金項目(編號：41374116，41674113)

劉秀娟(1980-)，女，工程師,主要從事地震資料處理及正反演工作。E-mail：121436068@qq.com

P631.4

A

2016-06-26