AVO地震反演中文培訓(xùn)

1、1會(huì)計(jì)學(xué)AVO地震反演中文培訓(xùn)地震反演中文培訓(xùn)2345678112211220VVVVR地震路徑地震路徑深度為深度為d的地質(zhì)的地質(zhì)界面界面 1 V1 2 V2t在時(shí)間為在時(shí)間為t = 2d/V1 地震反射地震反射地表地震子波地震子波泥巖泥巖氣砂巖氣砂巖9淺層泥巖和氣砂巖的速度差別.101120sin)(BRR123地表反射層 1 VP1 VS1 2 VP2 VS2 B 是梯度,是它產(chǎn)生了 AVO 響應(yīng). 它與密度, P波速度VP, 以及橫波速度VS 有關(guān).Zoeppritz1213P-waveDensityS-wavePoissons ratioSyntheticOffset Stack14巖

2、石物理研究開發(fā)了一些方法使得我們可以對(duì)巖石物理研究開發(fā)了一些方法使得我們可以對(duì)不同流體飽和度正演出測(cè)井曲線不同流體飽和度正演出測(cè)井曲線。在我們。在我們的軟件包中這個(gè)過程稱為的軟件包中這個(gè)過程稱為 FRM： 流體替換建模。（無橫波怎么辦？）流體替換建模。（無橫波怎么辦？） 然后對(duì)不同的流體情況創(chuàng)建合成道集并判斷該油氣藏然后對(duì)不同的流體情況創(chuàng)建合成道集并判斷該油氣藏是否適宜于做是否適宜于做AVO分析分析。BRINEGASOILINSITU正演油氣藏頂正演油氣藏頂?shù)渍穹渍穹?15Intercept: AGradient: B16Gradient (B)Intercept (A)17Near Inv

3、ersionFar Inversion將將AVO反演和反演和AVO分析結(jié)合分析結(jié)合起來進(jìn)行分析是一個(gè)非常重起來進(jìn)行分析是一個(gè)非常重要的識(shí)別油氣藏的手段要的識(shí)別油氣藏的手段18 輸入NMO動(dòng)校正后的道集處理分析反演 正 演截距、梯度屬性部分疊加Zoeppritz波動(dòng)方程法交繪圖LMR彈性阻抗聯(lián)合反演19SAMPLE IMAGE21孔隙流體孔隙流體Pores / Fluid基質(zhì)基質(zhì)Rock MatrixAVO響應(yīng)受孔隙型儲(chǔ)層的響應(yīng)受孔隙型儲(chǔ)層的縱波速度、橫波速度縱波速度、橫波速度及及密度密度的影響，如下所示，這涉及的影響，如下所示，這涉及到基質(zhì)、孔隙、及孔隙中的流體的性質(zhì)。到基質(zhì)、孔隙、及孔隙中的

4、流體的性質(zhì)。22)S1(S)1(whcwwmsat.,subscriptswaternhydrocarbomatrix,saturated,wsat,m,hc,saturationwaterwSporosity,density,where:下圖進(jìn)行解釋下圖進(jìn)行解釋.對(duì)密度的影響因素可用以下方程來描述對(duì)密度的影響因素可用以下方程來描述:232425P-wavesS-waves這就導(dǎo)致有兩種不同類型的速度這就導(dǎo)致有兩種不同類型的速度:P波（壓縮波）粒子運(yùn)動(dòng)方向和波傳輸方向一樣波（壓縮波）粒子運(yùn)動(dòng)方向和波傳輸方向一樣S波（剪切波）粒子運(yùn)動(dòng)方向和波傳輸方向成波（剪切波）粒子運(yùn)動(dòng)方向和波傳輸方向成90度

5、角度角262VPsV = 第一拉梅常量, = 第二拉梅常量, = 密度.P波和波和S波速度最簡(jiǎn)單的方程是從波速度最簡(jiǎn)單的方程是從非孔狀、各向同性的巖石中求得非孔狀、各向同性的巖石中求得。這個(gè)就是。這個(gè)就是用拉梅系數(shù)表示的方程用拉梅系數(shù)表示的方程:2734KVPsV另外一種常用的方程是用體模量或剪切模量來描述另外一種常用的方程是用體模量或剪切模量來描述:282222SPVV此處以上方程的反推，可以通過泊松比反求出VP與 VS 的比值:1222觀察觀察VP與與 VS 的比值通常是用泊松比的方法的比值通常是用泊松比的方法, 泊松比定義為泊松比定義為:29Vp/Vs vs Poissons Ratio

6、-0.2-012345678910Vp/VsPoissons RatioGas CaseWet Case30If VP/VS = 2, 則則 = 0If VP/VS = 1.5, 則則 = 0.1 (氣藏氣藏)If VP/VS = 2, 則則 = 1/3 (含水含水)If VP/VS = , 則則 = 0.5 (VS = 0)31)S(tSt)(ttwhcwwmsat11但以上方程，需要求解新的速度方程。Vtwhere/132Velocity vs Sw with Volume Avg. Eq.Por = 33%, Voil = 1300 m/s, Vga

7、s = 300 m/s 10001500200025003000350000.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.91Water SaturationVelocity (m/sec)OilGas33satsatsatsatPKV34_satsatsatSV_注意 sat是由前面討論的體平均方程式計(jì)算的34巖石基質(zhì)孔狀和流體干巖框架或骨架(孔隙間是空的) 飽和狀態(tài)巖石(孔隙中充滿流體) 35drysat干巖石的剪切模量飽和巖石的剪切模量drysat 在在Hampson-Russell AVO 軟件包中，軟件包中，Biot-Gassmann 分析是通過分析是通過 FRM

8、 (流體替流體替換建模換建模) 來完成的。首先看一些原理，然后考慮實(shí)際使用來完成的。首先看一些原理，然后考慮實(shí)際使用FRM選項(xiàng)時(shí)需要考慮的一選項(xiàng)時(shí)需要考慮的一些問題。些問題。在在Biot-Gassmann 方程中，如果方程中，如果孔隙度不變孔隙度不變，那么飽和度變化時(shí)剪切模量不改變，那么飽和度變化時(shí)剪切模量不改變:362211mdrymflmdrydrysatKKKKKKKK)(flmfldrymdrysatmsatKKKKKKKKKsat = 飽和巖石；飽和巖石；dry = 干巖框架；干巖框架；m = 巖石基質(zhì)；巖石基質(zhì)；fl = 流體；流體； =孔隙度孔隙度.(1)(2) Biot-Gas

9、smann 體模量方程如下體模量方程如下:37Ksandstone = 40 GPa, Klimestone = 60 GPa.現(xiàn)在分析在Biot-Gassmann方程中如何估算各種體模量，先計(jì)算固體巖石基質(zhì)的體模量。值是以 (GPa)為單位，等于1010 dynes/cm2。固體巖石基質(zhì)的體模量固體巖石基質(zhì)的體模量Km 通常都是從巖心鉆探樣本通常都是從巖心鉆探樣本(試驗(yàn)室試驗(yàn)室)測(cè)量所公布的數(shù)據(jù)中測(cè)量所公布的數(shù)據(jù)中得來，標(biāo)準(zhǔn)值有得來，標(biāo)準(zhǔn)值有:38hcwwwflKSKSK11Batzle 和 Wang，1992給出了估算水、氣、油體模量的方程。論文名：Seismic Properties o

10、f Pore Fluids, Geophysics, 57, 1396-1408. 標(biāo)準(zhǔn)值有：Kgas = 0.021 GPa； Koil = 0.79 GPa； Kw = 2.38 GPaflwhcwhere the bulk modulus of the fluid, the bulk modulus of the water,and the bulk modulus of the hydrocarbon.KKK流體體模量流體體模量 Kfl 可用以下方程計(jì)算可用以下方程計(jì)算:流體的體模量（要求的）流體的體模量（要求的）水的體模量水的體模量油氣的體模量油氣的體模量39)(flmfldrymd

11、rysatmsatKKKKKKKKKMavko 方程40后面的幾張片子后面的幾張片子,我們將看一下利用我們將看一下利用Biot-Gassmann 方程方程,計(jì)算出的含油和含氣的計(jì)算出的含油和含氣的不同結(jié)果不同結(jié)果.我們還將分析飽和度對(duì)速度我們還將分析飽和度對(duì)速度 (VP 和和VS)以及泊松比的影響以及泊松比的影響.這個(gè)模型是假定氣在流體中是均勻分布的，不調(diào)和飽和狀態(tài)將提供不同的函數(shù)這個(gè)模型是假定氣在流體中是均勻分布的，不調(diào)和飽和狀態(tài)將提供不同的函數(shù). (參見參見 Mavko et al: The Rock Physics Handbook.)41Velocity vs Sw - Gas Cas

12、e, Por = 33%Ks = 40, Kgas = .021, Kdry = 3.25, Mu = 3.3 GPa10001200140016001800200022002400260001Sw Velocity (m/s)VpVs42Poissons Ratio vs Water Saturation Gas Case00.40.501Sw Poissons RatioPoissons Ratio43Velocity vs Sw - Oil CasePorosit

13、y = 33%, Koil = 1.0 GPa1000150020002500300001Sw Velocity (m/s)VsVp44Poissons Ratio vs Water Saturation Oil Case00.40.501Sw Poissons RatioPoissons Ratio4546當(dāng)出現(xiàn)多種孔狀流體當(dāng)出現(xiàn)多種孔狀流體, Kfl 通常是用通常是用 Reuss 平均方程進(jìn)行計(jì)算平均方程進(jìn)行計(jì)算:這個(gè)平均方程這個(gè)平均方程假定流體假定流體平均分布平均

14、分布!氣和流體在每個(gè)孔隙中氣和流體在每個(gè)孔隙中均勻分布。均勻分布。這種方法認(rèn)為組合后流體的可壓縮性很大程度上取決于最易壓縮的那部分流體項(xiàng)：這種方法認(rèn)為組合后流體的可壓縮性很大程度上取決于最易壓縮的那部分流體項(xiàng)：ggoowwflKSKSKSK147不調(diào)和的狀態(tài)，Voigt 平均方程可以最好地估算流體模量 Kfl (Domenico, 1976).當(dāng)不調(diào)和狀態(tài)出現(xiàn)當(dāng)不調(diào)和狀態(tài)出現(xiàn)介于兩者的中等情形介于兩者的中等情形, 可對(duì)每個(gè)不調(diào)和區(qū)域進(jìn)行可對(duì)每個(gè)不調(diào)和區(qū)域進(jìn)行Gassmann 流體流體替換然后求一個(gè)體平均替換然后求一個(gè)體平均 (Dvorkin et al, 1999)，這做法類似于使用，這做法類

15、似于使用power-law 平均技平均技術(shù)估算的結(jié)果術(shù)估算的結(jié)果.ggoowwflKSKSKSK48Gassmann 預(yù)測(cè)的速度非固化砂巖基質(zhì)Porosity = 30%100% Gas to 100% Brine saturation49SPV1222V后面進(jìn)行解釋后面進(jìn)行解釋對(duì)于常泊松比對(duì)于常泊松比, 截距為零截距為零:smVVSP/136016. 1泥巖線是泥巖線是VP 和和VS之間的線性關(guān)系之間的線性關(guān)系,是由是由Castagna et al (1985)提出的提出的:50ARCOs original mudrock derivation(Castagna et al, Geophys

16、ics, 1985)51020002000400060001000300040000100030005000VP (m/s)VS(m/s)Mudrock LineGas Sand52020002000400060001000300040000100030005000VP (m/s)VS(m/s)Mudrock LineGas Sand= 1/3 （VP/VS = 2）含水含水53VP (m/s)020002000400060001000300040000100030005000VS(m/s)泥巖線泥巖線氣砂巖氣砂巖 = 0.1 或或VP/VS = 1.5= 1/3 （VP/VS = 2）含水含

17、水54VPVSSW = 50%VPVSSW = 100%輸入測(cè)井曲線輸入測(cè)井曲線輸出測(cè)井曲線輸出測(cè)井曲線(1) 流體替換流體替換Biot-Gassmann 方程的基本應(yīng)用是在一組目的層內(nèi)將某種流體替換成另一種流體方程的基本應(yīng)用是在一組目的層內(nèi)將某種流體替換成另一種流體. 這種情形這種情形VP VS 都已知，流體的飽和度也已知都已知，流體的飽和度也已知 (SW)。通常替換后三種曲線在目標(biāo)處都有變。通常替換后三種曲線在目標(biāo)處都有變化?；?5VPSW = 50%VPVS輸入曲線輸入曲線輸出曲線輸出曲線(2) 計(jì)算Vs應(yīng)用 Biot-Gassmann 方程式的第二個(gè)目的是計(jì)算 VS 曲線( 測(cè)井中沒

18、有測(cè)量)，干巖 Kdry 已知或假定泥巖線方程適合于水砂巖。 SW = 50%56VPVPVS-Cast事實(shí)上，計(jì)算事實(shí)上，計(jì)算 Vs需要兩個(gè)不同的方程：需要兩個(gè)不同的方程：Castagnas 方程用于計(jì)算水砂巖方程用于計(jì)算水砂巖/泥巖泥巖(目標(biāo)區(qū)域外），目標(biāo)區(qū)域外），Biot-Gassmann 方程用于目標(biāo)區(qū)域內(nèi)。計(jì)算分為兩步方程用于目標(biāo)區(qū)域內(nèi)。計(jì)算分為兩步:(1)首先使用首先使用Castagnas 方方程計(jì)算程計(jì)算 VS 曲線曲線, 結(jié)果在結(jié)果在目標(biāo)區(qū)域外準(zhǔn)確目標(biāo)區(qū)域外準(zhǔn)確.(2)利用利用 Biot-Gassmann 計(jì)算目標(biāo)區(qū)域內(nèi)準(zhǔn)確的計(jì)算目標(biāo)區(qū)域內(nèi)準(zhǔn)確的 VS 曲線曲線.VS-Cast

19、VPVSVPSAMPLE IMAGE58第一組練習(xí)數(shù)據(jù)來自于加拿大西部地區(qū)第一組練習(xí)數(shù)據(jù)來自于加拿大西部地區(qū)colony白堊系砂巖地層。白堊系砂巖地層。目標(biāo)是一個(gè)目標(biāo)是一個(gè)8米厚的薄層氣砂米厚的薄層氣砂巖。巖。 我們只有一口井，這口井有測(cè)我們只有一口井，這口井有測(cè)量的縱波和密度曲線；量的縱波和密度曲線；砂巖已知含水飽和度砂巖已知含水飽和度50% ，含，含氣飽和度氣飽和度50%.在這個(gè)練習(xí)中，我們要通過在這個(gè)練習(xí)中，我們要通過Castagnas 方程和方程和Biot-Gassmann分析生成一個(gè)橫波分析生成一個(gè)橫波。59啟動(dòng)Geoview軟件接下來創(chuàng)建一個(gè)新的數(shù)據(jù)庫，點(diǎn)擊New生成一個(gè)新的數(shù)據(jù)庫

20、。下一個(gè)菜單，輸入數(shù)據(jù)庫的名字下一個(gè)菜單，輸入數(shù)據(jù)庫的名字avo_class，然后點(diǎn)擊，然后點(diǎn)擊 OK.60數(shù)據(jù)庫是一個(gè)井集，每口井可能有多條各種類型的曲線，它們需要轉(zhuǎn)換成深度域的值，深度數(shù)據(jù)庫是一個(gè)井集，每口井可能有多條各種類型的曲線，它們需要轉(zhuǎn)換成深度域的值，深度域不一定等間隔采樣，域不一定等間隔采樣，GEOVIEW窗口允許我們查看這個(gè)數(shù)據(jù)庫，此時(shí)是空的：窗口允許我們查看這個(gè)數(shù)據(jù)庫，此時(shí)是空的：61與與Well Explorer window一起出現(xiàn)的還有軟件主菜單如下：一起出現(xiàn)的還有軟件主菜單如下：62要導(dǎo)入井曲線，點(diǎn)擊 63在導(dǎo)入菜單，選擇含有在導(dǎo)入菜單，選擇含有AVO工作數(shù)據(jù)的目錄，

21、選擇工作數(shù)據(jù)的目錄，選擇 avo_well.las 這個(gè)文件：這個(gè)文件：確定曲線格式設(shè)置成了確定曲線格式設(shè)置成了LAS，然后點(diǎn)擊，然后點(diǎn)擊Next . 64下一頁要求你定義這口井的參數(shù)，點(diǎn)擊下一頁要求你定義這口井的參數(shù)，點(diǎn)擊Next 接受默認(rèn)選項(xiàng)。接受默認(rèn)選項(xiàng)。65666768如果你點(diǎn)擊曲線上方的顯示窗口上的如果你點(diǎn)擊曲線上方的顯示窗口上的 View / Display Options，會(huì)看到一個(gè)菜單允許你改變顯，會(huì)看到一個(gè)菜單允許你改變顯示：示：點(diǎn)擊上面的標(biāo)簽，可以點(diǎn)擊上面的標(biāo)簽，可以變換修改項(xiàng)，注意此時(shí)變換修改項(xiàng)，注意此時(shí)可用的選項(xiàng)是可用的選項(xiàng)是Layout, Scale & De

22、tails, and Synthetics.69Layout菜單允許你重疊顯示菜單允許你重疊顯示曲線，如右圖選擇三條曲線曲線，如右圖選擇三條曲線 Computed Impedance 、Computed Reflectivity 與與P-wave 顯示在同一列，點(diǎn)顯示在同一列，點(diǎn)擊擊Apply應(yīng)用應(yīng)用7071現(xiàn)在點(diǎn)擊現(xiàn)在點(diǎn)擊Curves 標(biāo)簽設(shè)置值域范圍標(biāo)簽設(shè)置值域范圍 Density 曲線設(shè)為曲線設(shè)為 1.8 到到2.7 ， SP 曲線曲線 -150 到到 -50 ：設(shè)置完后，點(diǎn)擊設(shè)置完后，點(diǎn)擊Save settings as project template. 這意這意味著將來的顯示味著將

23、來的顯示中默認(rèn)的這些曲中默認(rèn)的這些曲線將一直沿用這線將一直沿用這個(gè)值域范圍。個(gè)值域范圍。72點(diǎn)擊點(diǎn)擊Ok ，顯示窗口出，顯示窗口出現(xiàn)：現(xiàn)：73 Well Explorer 提供了另外一種看數(shù)據(jù)庫的方式：提供了另外一種看數(shù)據(jù)庫的方式：數(shù)據(jù)庫中每口井對(duì)應(yīng)一行，可以修改相應(yīng)的參數(shù)。數(shù)據(jù)庫中每口井對(duì)應(yīng)一行，可以修改相應(yīng)的參數(shù)。74點(diǎn)擊井名左側(cè)箭頭：點(diǎn)擊井名左側(cè)箭頭：現(xiàn)在可以看到這口井的一系列曲線，這些參數(shù)也可以進(jìn)行修改現(xiàn)在可以看到這口井的一系列曲線，這些參數(shù)也可以進(jìn)行修改.75現(xiàn)在對(duì)井現(xiàn)在對(duì)井AVO_WELL應(yīng)用應(yīng)用Check shot校正。首先創(chuàng)建一個(gè)初始校正。首先創(chuàng)建一個(gè)初始 check shot

24、 曲線。點(diǎn)擊曲線。點(diǎn)擊Log Options / Create a new log in table / Check Shot:76會(huì)看到該井會(huì)出現(xiàn)新一行井曲線，它有相應(yīng)的單位及類型會(huì)看到該井會(huì)出現(xiàn)新一行井曲線，它有相應(yīng)的單位及類型:To begin entering values, click on the arrow to the left of the Log Name field:77填入值填入值:這個(gè)這個(gè) check shot 曲線包含兩列，深度及時(shí)間。點(diǎn)擊曲線包含兩列，深度及時(shí)間。點(diǎn)擊 Update 更新，然后點(diǎn)擊更新，然后點(diǎn)擊 Yes 就可以把就可以把曲線值更改進(jìn)去了。曲線值更改

25、進(jìn)去了。78即使添加了即使添加了check shot 值，值，聲波曲線也不會(huì)被應(yīng)用，除聲波曲線也不會(huì)被應(yīng)用，除非我們應(yīng)用時(shí)深校正非我們應(yīng)用時(shí)深校正 。79要應(yīng)用時(shí)深校正，點(diǎn)擊要應(yīng)用時(shí)深校正，點(diǎn)擊Option / Check Shot Correction: 時(shí)深校正參數(shù)菜單和校正分析窗口如圖：時(shí)深校正參數(shù)菜單和校正分析窗口如圖：80有許多地方要進(jìn)行設(shè)置才能讓它自動(dòng)運(yùn)算。一處要求刪除有許多地方要進(jìn)行設(shè)置才能讓它自動(dòng)運(yùn)算。一處要求刪除Check shot數(shù)據(jù)點(diǎn)重新應(yīng)用相關(guān)數(shù)據(jù)點(diǎn)重新應(yīng)用相關(guān)，另一處要求在曲線校正中應(yīng)用點(diǎn)與點(diǎn)之間的多項(xiàng)式插值來預(yù)防假的不連續(xù)現(xiàn)象。現(xiàn)在我們，另一處要求在曲線校正中應(yīng)用點(diǎn)

26、與點(diǎn)之間的多項(xiàng)式插值來預(yù)防假的不連續(xù)現(xiàn)象?，F(xiàn)在我們采用默認(rèn)的參數(shù)，點(diǎn)擊采用默認(rèn)的參數(shù)，點(diǎn)擊Ok ，如圖：，如圖：點(diǎn)完點(diǎn)完 Ok ，將會(huì)創(chuàng)建一個(gè)新的校正后的聲波曲線，名稱，將會(huì)創(chuàng)建一個(gè)新的校正后的聲波曲線，名稱 P-wave_chk81現(xiàn)在我們開始用剛才導(dǎo)入的井曲線創(chuàng)建一個(gè)現(xiàn)在我們開始用剛才導(dǎo)入的井曲線創(chuàng)建一個(gè)AVO模模型，點(diǎn)擊型，點(diǎn)擊 AVO/AVO 按鈕啟動(dòng)按鈕啟動(dòng)AVO選擇選項(xiàng)選擇選項(xiàng) Start New Project:新名稱新名稱 colony:82出現(xiàn)一個(gè)初始的空白出現(xiàn)一個(gè)初始的空白AVO程序窗口程序窗口83選擇井選擇井 (AVO_WELL) ，點(diǎn)擊，點(diǎn)擊 Open:當(dāng)當(dāng)AVO窗口出

27、現(xiàn)時(shí)，點(diǎn)擊窗口出現(xiàn)時(shí)，點(diǎn)擊 Modeling / Single Well:84這個(gè)菜單，我們希望激活默認(rèn)的輸入曲線這個(gè)菜單，我們希望激活默認(rèn)的輸入曲線:下一個(gè)信息窗口告訴你需要一個(gè)橫波曲下一個(gè)信息窗口告訴你需要一個(gè)橫波曲線，它將要被創(chuàng)建，點(diǎn)擊線，它將要被創(chuàng)建，點(diǎn)擊 OK.點(diǎn)擊點(diǎn)擊OK.85在下一菜單，接受默認(rèn)選項(xiàng)，選擇應(yīng)用在下一菜單，接受默認(rèn)選項(xiàng)，選擇應(yīng)用Castagnas equation:Click Next 86我們應(yīng)用經(jīng)過校正后的縱波曲線進(jìn)行換算并且給輸出的橫波曲線命名我們應(yīng)用經(jīng)過校正后的縱波曲線進(jìn)行換算并且給輸出的橫波曲線命名 S-wave_Castagna. 改變名稱如下：改變名稱

28、如下：Click Next 87應(yīng)用默認(rèn)的應(yīng)用默認(rèn)的Castagna coefficients參數(shù)參數(shù):Finally, click OK.88 AVO模型窗口出現(xiàn)，顯示了生成的橫波曲線及運(yùn)算得到的泊松比曲線。模型窗口出現(xiàn)，顯示了生成的橫波曲線及運(yùn)算得到的泊松比曲線。89橫波曲線是用橫波曲線是用Castagna泥巖方程運(yùn)算得到的，它只適用于水砂巖和泥巖地層。這意味著對(duì)泥巖方程運(yùn)算得到的，它只適用于水砂巖和泥巖地層。這意味著對(duì)于目標(biāo)砂巖層這個(gè)結(jié)果是錯(cuò)誤的。于目標(biāo)砂巖層這個(gè)結(jié)果是錯(cuò)誤的。我們將要應(yīng)用我們將要應(yīng)用Biot-Gassmann方程來替換此段生成的適合含氣砂巖的橫波速度。方程來替換此段生成

29、的適合含氣砂巖的橫波速度。 對(duì)于這個(gè)運(yùn)算，我們需要關(guān)于實(shí)際砂巖流體部分進(jìn)一步的信息，如含水飽各度對(duì)于這個(gè)運(yùn)算，我們需要關(guān)于實(shí)際砂巖流體部分進(jìn)一步的信息，如含水飽各度50%.要進(jìn)行流體替換，點(diǎn)擊要進(jìn)行流體替換，點(diǎn)擊FRM 按鈕按鈕90919293949596979899100101SAMPLE IMAGE103以上示意圖顯示了以上示意圖顯示了 (a) P波波, 壓縮波壓縮波 (b) SH波波, 水平橫波水平橫波(c) SV波波, 垂向橫波垂向橫波 （用橫波震源激發(fā)接收到的）（用橫波震源激發(fā)接收到的）. (Ensley, 1984)(a)(b)(c)104105反射P-波 = RP反射SV-波透射

30、 P-波入射 P-波透射 SV-波VP1 , VS1 , 1VP2 , VS2 , 2irrtt如果入射角如果入射角 0, 入射入射 P波將產(chǎn)生波將產(chǎn)生P 波和波和 SV 反射以及透射反射以及透射. 這就稱之為波形轉(zhuǎn)這就稱之為波形轉(zhuǎn)換換.106107112211222221221121122211121112222112211111sincossincoscossincossinsin2cos2sin2cos2cos2sin2cos2sin2PSSPSPPPSSPSSSSPPPPRRV VV VVTVV VVTVVVVVV11111sincossin2cos2Zoeppritz 應(yīng)用應(yīng)力應(yīng)變方

31、程導(dǎo)出了不同地層界面對(duì)應(yīng)的反射波和透射波振幅變化方程，它給出了4個(gè)未知數(shù)的4個(gè)方程，Zoeppritz方程給我們隨角度變化的準(zhǔn)確的振幅。 which gives four equations with four unknowns. Inverting the matrix form of theZoeppritz equations gives us the exact amplitudes as a function of angle:108SSPPVVcbVVaR)(此處此處:,sin4,sin25 . 0,cos2122222PSPSVVcVVba.2,2,2,2121212121212

32、tiSSSSSSPPPPPPandVVVVVVVVVVVV109222sintanCsinBA)(Rwhere:pP2PSSS2PSpPpPVV21CVV2VVVV4VV21BVV21AWiggins導(dǎo)出了一種更直觀，基本上相同的一種形式，他把方程劃分為三個(gè)反射項(xiàng)，后面的比前面的影響要?。?10第一項(xiàng)A是零偏移距反射系數(shù)，與縱波速度Vp與密度 乘積有關(guān)第二項(xiàng)B是梯度，與 sin2相乘，在隨偏移距變化的方程中它起最大作用。它受VP、VS與 變化的影響。第三項(xiàng)C稱作曲率curvature 只與縱波速度有關(guān)。它與tan2*sin2相乘對(duì)振幅影響在30度角度以內(nèi)幾乎可以忽略。 (注意: 其實(shí)tan2*

33、sin2 應(yīng)該是 tan2 - sin2, 兩者近似相等，方程通常寫成第一種形式.)111Ostrander (1984) 是最早提出含氣砂巖有AVO效應(yīng)的人之一，并且提出了一個(gè)簡(jiǎn)單的雙層模型，兩層高阻抗高泊松比泥巖中間夾一層低阻抗、低泊松比砂巖.下一頁會(huì)顯示這個(gè)模型。Ostrander的模型在 Sacramento 河谷氣砂巖處應(yīng)用非常好。然而，它只對(duì)一種AVO （3類）異常起作用，其它類型將在接下來的章節(jié)介紹。112注意：模型包含一套低波阻抗低泊松比氣砂巖及上下兩套泥巖注意：模型包含一套低波阻抗低泊松比氣砂巖及上下兩套泥巖113(a) Well log responses for the

34、model.(b) Synthetic seismic.114115(a) Wet model(b) Gas model116接下來兩頁，我們將要從模型接下來兩頁，我們將要從模型A和和B（應(yīng)用下面的值）生成頂?shù)淄噍S，水和氣層（應(yīng)用下面的值）生成頂?shù)淄噍S，水和氣層段要用段要用Biot-Gassmann 方程來計(jì)算。方程來計(jì)算。Wet: VP= 2500 m/s, VS= 1250 m/s, = 2.11 g/cc, = 0.33 Gas: VP= 2000 m/s, VS= 1310 m/s, = 1.95 g/cc, = 0.12 Shale: VP= 2250 m/s, VS= 1125

35、 m/s, = 2.0 g/cc, = 0.33 我們將要考慮我們將要考慮Aki-Richards方程中有無第三項(xiàng)的情況下的方程中有無第三項(xiàng)的情況下的AVO效應(yīng)。效應(yīng)。117這個(gè)表顯示了這個(gè)表顯示了水砂巖水砂巖模型應(yīng)模型應(yīng)用用Aki-Richards三項(xiàng)式三項(xiàng)式方程方程（A，B，C）在模型頂?shù)咨┰谀Ｐ晚數(shù)咨傻某傻腁VO曲線，然后是曲線，然后是2項(xiàng)式項(xiàng)式。注意在。注意在25度度以后兩者才分以后兩者才分離。離。118這個(gè)表顯示了這個(gè)表顯示了含氣砂巖含氣砂巖模型模型應(yīng)用應(yīng)用Aki-Richards三項(xiàng)式三項(xiàng)式方方程程（A，B，C）在模型頂?shù)祝┰谀Ｐ晚數(shù)咨傻纳傻腁VO曲線，然后是曲線，然后是2

36、項(xiàng)項(xiàng)式。注意在式。注意在25度度以后兩者才以后兩者才分離。分離。1192)1(121)D1(2DAB1212PPPP2,/V/VV/VD:where（A+B）/2=22222SPSPVVVV120Gas Sand ModelAki-Richards vs Shuey-0.250-0.200-0.150-0.100-0.0500.0000.0500.1000.1500.2000.250051015202530354045Angle (degrees)AmplitudeA-R TopShuey TopA-R BaseShuey Base這個(gè)表顯示了對(duì)于氣砂巖兩種這個(gè)表顯示了對(duì)于氣砂巖兩種方程（方程

37、（Aki與與Shuey）的對(duì)比）的對(duì)比。121SAMPLE IMAGE123124125126Well to Seismic 菜單顯示菜單顯示:這個(gè)菜單告訴程序這個(gè)菜單告訴程序Geoview數(shù)據(jù)庫中如何定位井在地震上的位置，改變參數(shù)如上，輸入數(shù)據(jù)庫中如何定位井在地震上的位置，改變參數(shù)如上，輸入CDP330。點(diǎn)擊。點(diǎn)擊Ok127移動(dòng)下部的滾動(dòng)條到移動(dòng)下部的滾動(dòng)條到 CDP 325-335.道集顯示道集顯示:128釋放鼠標(biāo)后顯示菜單如下，在釋放鼠標(biāo)后顯示菜單如下，在X330處顯示的井曲線是縱波速度：處顯示的井曲線是縱波速度：你會(huì)注意到同相軸在你會(huì)注意到同相軸在630ms時(shí)間處線的中間比左側(cè)要強(qiáng)時(shí)間

38、處線的中間比左側(cè)要強(qiáng)，遠(yuǎn)偏移距比近偏移距振幅，遠(yuǎn)偏移距比近偏移距振幅有明顯的增強(qiáng)。這是有明顯的增強(qiáng)。這是AVO異異常，在接下來的練習(xí)中會(huì)研常，在接下來的練習(xí)中會(huì)研究它。究它。129接下來，我們將在目標(biāo)區(qū)附近創(chuàng)建一個(gè)超道集，或者共偏接下來，我們將在目標(biāo)區(qū)附近創(chuàng)建一個(gè)超道集，或者共偏移距道集。選擇移距道集。選擇 Process / Super Gather 如右所示：如右所示：調(diào)出超道集創(chuàng)建窗口，保持默認(rèn)參調(diào)出超道集創(chuàng)建窗口，保持默認(rèn)參數(shù)數(shù)Next .130131132133134135136137138139點(diǎn)擊兩次點(diǎn)擊兩次 Next , 然后改變子波長(zhǎng)度為然后改變子波長(zhǎng)度為100然后點(diǎn)擊然后點(diǎn)

39、擊OK.140141現(xiàn)在我們將要?jiǎng)?chuàng)建一個(gè)偏移距合成道現(xiàn)在我們將要?jiǎng)?chuàng)建一個(gè)偏移距合成道點(diǎn)擊點(diǎn)擊 Synthetic / Zoeppritz:142點(diǎn)擊點(diǎn)擊 Next 然后然后OK 在在AVO模型窗口中將會(huì)顯示合成道集模型窗口中將會(huì)顯示合成道集注意合成記錄與實(shí)際地震道相比有一點(diǎn)過比例化了，要修正這個(gè)，點(diǎn)眼睛修改參數(shù)菜單。注意合成記錄與實(shí)際地震道相比有一點(diǎn)過比例化了，要修正這個(gè)，點(diǎn)眼睛修改參數(shù)菜單。143指到指到 Seismic Views 標(biāo)簽標(biāo)簽,如右圖所示如右圖所示.點(diǎn)亮點(diǎn)亮 avo_well_syn 并且改并且改變變 Excursion 為為 0.70.然后點(diǎn)擊然后點(diǎn)擊 Ok 應(yīng)用改變應(yīng)用改

40、變144比例化后新的顯示如下：比例化后新的顯示如下：SAMPLE IMAGE1461472sinBA)(R舍棄掉舍棄掉C項(xiàng)，定義項(xiàng)，定義A和和B如下：如下：,VV2VVVV4VV21B2PSSS2PSpP./V/VV/VD,)1(121)D1(2DABPPPP2,VV21ApP或者應(yīng)用或者應(yīng)用Shuey近似方程：近似方程：148應(yīng)用應(yīng)用2項(xiàng)式比較普遍，因?yàn)椋喉?xiàng)式比較普遍，因?yàn)椋?1) 分析相對(duì)簡(jiǎn)單分析相對(duì)簡(jiǎn)單It simplifies the analysis considerably(2) 對(duì)于角度小于對(duì)于角度小于40度，第三項(xiàng)沒什么意義，如前所述：度，第三項(xiàng)沒什么意義，如前所述：1492s

41、inBA)(R150Offset (m)4506000Angle (degrees)090151從偏移距道集轉(zhuǎn)換成角度道集可用以下方法之一從偏移距道集轉(zhuǎn)換成角度道集可用以下方法之一:(1) 直線法直線法 (常速度常速度)(2) 射線追蹤射線追蹤 (可變速度可變速度)(3) 射線參數(shù)逼近射線參數(shù)逼近(可變速度逼近可變速度逼近).velocityRMSV,timeway2t,2tVdepthd,offsetXwhere,VtXd2Xtan:RayStraight)1(RMS00RMS0.traveltimetotalt,velocityIntervalVwhere,tVXVsin:Paramete

42、rRay)2(INT2RMSINT152分析地震數(shù)據(jù)時(shí)通常是應(yīng)用射線追蹤分析地震數(shù)據(jù)時(shí)通常是應(yīng)用射線追蹤法（法（Ray parameter）.因?yàn)樗纫驗(yàn)樗萊ay Tracing法運(yùn)算快法運(yùn)算快.兩者相似，但當(dāng)角度大于兩者相似，但當(dāng)角度大于50度時(shí)射線度時(shí)射線追蹤法的結(jié)果衰退相對(duì)微弱。追蹤法的結(jié)果衰退相對(duì)微弱。右邊的實(shí)例顯示了用直線法和射線追右邊的實(shí)例顯示了用直線法和射線追蹤法計(jì)算的角度等值線。蹤法計(jì)算的角度等值線。Ray TracingRay Parameter5050353543436565153154155CDP疊加道在含氣井處顯示如下，井曲線是聲波曲線。注意“亮點(diǎn)”可能對(duì)應(yīng)著氣層（也

43、可能是特殊巖性引起的）“Bright Spot”156(a) 在含氣區(qū)域在含氣區(qū)域?qū)?yīng)的一系列校對(duì)應(yīng)的一系列校正后的正后的CDP道集道集(b) 同樣的道集，同樣的道集，但是疊加顯示了但是疊加顯示了顏色振幅包絡(luò)。顏色振幅包絡(luò)。157(a) 共偏移距疊加共偏移距疊加(b) 沿波谷提取的振幅(c) 沿波峰提取的振幅158Offset+A-A+B- Bsin2TimeAki-Richards 方程預(yù)測(cè)了振幅隨方程預(yù)測(cè)了振幅隨sin2變變化的線性關(guān)系。化的線性關(guān)系。對(duì)同一個(gè)時(shí)間采樣，會(huì)生成回歸曲線，并對(duì)同一個(gè)時(shí)間采樣，會(huì)生成回歸曲線，并以此來預(yù)測(cè)以此來預(yù)測(cè)A和和B值。值。可以在任意時(shí)間拾取振幅，其中有兩

44、個(gè)如下可以在任意時(shí)間拾取振幅，其中有兩個(gè)如下所示：所示：159運(yùn)算的結(jié)果是生成兩種基本屬性體Intercept: AGradient: B160應(yīng)用最廣的應(yīng)用最廣的AVO屬性有：屬性有：(1) AVO 乘積屬性：乘積屬性：A*B(2) 比例化的泊松比變化屬性：比例化的泊松比變化屬性：A+B(3) 橫波阻抗反射系數(shù)：橫波阻抗反射系數(shù)： A-B(4) 流體因子流體因子原始的原始的A和和B屬性體用的相對(duì)較少，相對(duì)而言由它們生成的其它的屬性體用的相對(duì)較少，相對(duì)而言由它們生成的其它的AVO屬性經(jīng)常會(huì)屬性經(jīng)常會(huì)被使用。被使用。161許多許多AVO異常通常有如右圖所示形式異常通常有如右圖所示形式.這種情況下

45、，截距和梯度值都比較大，即所這種情況下，截距和梯度值都比較大，即所謂的謂的“亮點(diǎn)亮點(diǎn)”。同時(shí)，它們的極性相同。同時(shí)，它們的極性相同。這是一個(gè)這是一個(gè)3類類AVO異常的例子。異常的例子。生成生成A*B: 砂巖頂：砂巖頂：(-A)*(-B) = +AB 砂巖底砂巖底: (+A)*(+B) = +AB它給出了一個(gè)在頂?shù)锥加械恼o出了一個(gè)在頂?shù)锥加械恼傲咙c(diǎn)亮點(diǎn)” 響應(yīng)響應(yīng)。+A-A+B- Bsin2162AVO乘積在儲(chǔ)層頂?shù)罪@示了一個(gè)正響應(yīng)乘積在儲(chǔ)層頂?shù)罪@示了一個(gè)正響應(yīng)TopBase1632( )sinRAB12PpVAV2)1(121)D1(2DAB2192(1)2(2/3)4BA DDA從從S

46、huey方程得到的第二個(gè)組合方程得到的第二個(gè)組合如果我們假設(shè)背景泊松比趨勢(shì)如果我們假設(shè)背景泊松比趨勢(shì) = 1/3, 那么那么:或或: 94AB求和求和 A+B 是比例化后的泊松比變化是比例化后的泊松比變化164AVO求和屬性求和屬性 (A+B)在儲(chǔ)層頂顯示負(fù)響應(yīng)在儲(chǔ)層頂顯示負(fù)響應(yīng) (泊松比泊松比 降低降低) ，在儲(chǔ)層底顯示正響應(yīng)，在儲(chǔ)層底顯示正響應(yīng) (泊松比泊松比 增加增加)：TopBase1652( )sinRAB第三個(gè)組合是從第三個(gè)組合是從 Aki_Richards 方程求得的：方程求得的：如果假定背景趨勢(shì)如果假定背景趨勢(shì)VP/VS = 2, 那么那么:22121422PpPSSSpPSP

47、VAVVVVVBVVVV22142212,:,2PSSSpPSPSPSPSSVVVVBVVVVVRRwhere RAand RVor: A-B 是比例化后的橫波阻抗反射系數(shù)是比例化后的橫波阻抗反射系數(shù)BA21RS166TopBaseAVO差屬性差屬性 (A-B)在儲(chǔ)層頂對(duì)應(yīng)橫波阻抗增加，這個(gè)運(yùn)算經(jīng)常用到更準(zhǔn)確的在儲(chǔ)層頂對(duì)應(yīng)橫波阻抗增加，這個(gè)運(yùn)算經(jīng)常用到更準(zhǔn)確的Fatti方程，接下來會(huì)看到。方程，接下來會(huì)看到。1671222212( ):1tan,8 sin,11,.22PSSPPSPSPSRc Rc RVwhereccVVVRRVV 允許我們（對(duì)地震數(shù)據(jù)）用生成A和B時(shí)相同的方式來生成RP 與

48、 RS 體168RPRS169RP/RS 屬性經(jīng)?？梢赞D(zhuǎn)變成以下屬性經(jīng)常可以轉(zhuǎn)變成以下2種新屬性種新屬性(1) 流體因子流體因子, 接下來會(huì)討論接下來會(huì)討論.(2) AVO反演或反演或Lambda-Mu-Rho，下面的章節(jié)會(huì)討論，下面的章節(jié)會(huì)討論170注意因子注意因子 經(jīng)常被標(biāo)準(zhǔn)化以符合本區(qū)數(shù)據(jù)經(jīng)常被標(biāo)準(zhǔn)化以符合本區(qū)數(shù)據(jù)smVVSP/136016. 1流體因子屬性流體因子屬性 (Smith and Gidlow, 1987, Fatti et al., 1994) 是基于是基于Castagnas 泥巖方程，它假設(shè)非油氣區(qū)地層是正確的。泥巖方程，它假設(shè)非油氣區(qū)地層是正確的。1.16SPSPVRR

49、V從這個(gè)方程式，我們可以導(dǎo)出以下同等的方程從這個(gè)方程式，我們可以導(dǎo)出以下同等的方程, SPSPRVV16.1RF流體因子的定義是為了突出顯示流體因子的定義是為了突出顯示Castagna方程不符合的地層，如潛在的油氣區(qū)域：方程不符合的地層，如潛在的油氣區(qū)域：1.16SPVV171AVO流體因子顯示了在儲(chǔ)層頂與儲(chǔ)層底流體因子顯示了在儲(chǔ)層頂與儲(chǔ)層底 相對(duì)相對(duì) Castagna方程的嚴(yán)重背離程度。方程的嚴(yán)重背離程度。TopBase172角度計(jì)算中運(yùn)用的角度計(jì)算中運(yùn)用的速度速度角度范圍角度范圍最小相關(guān)系數(shù)最小相關(guān)系數(shù)for handling noisy picks.分析類型分析類型SAMPLE IMAG

50、E174首先，我們要用一套真實(shí)數(shù)據(jù)來計(jì)算角度道集并且可以對(duì)輸入的速度在結(jié)果上進(jìn)行檢查，超首先，我們要用一套真實(shí)數(shù)據(jù)來計(jì)算角度道集并且可以對(duì)輸入的速度在結(jié)果上進(jìn)行檢查，超道集在道集在AVO地震數(shù)據(jù)窗口的顯示如下：地震數(shù)據(jù)窗口的顯示如下：這個(gè)數(shù)據(jù)是有一個(gè)強(qiáng)這個(gè)數(shù)據(jù)是有一個(gè)強(qiáng)AVO異常的異常的CDP道集數(shù)據(jù)，為了做道集數(shù)據(jù)，為了做AVO分析，我們必須確定在目的層的分析，我們必須確定在目的層的入射角范圍要合理。要看這個(gè)，接下來我們要將數(shù)據(jù)轉(zhuǎn)換成角度域數(shù)據(jù)。入射角范圍要合理。要看這個(gè)，接下來我們要將數(shù)據(jù)轉(zhuǎn)換成角度域數(shù)據(jù)。.Anomaly175點(diǎn)擊點(diǎn)擊 Process / Angle Gather:在第一

51、頁菜單，保持默認(rèn)選項(xiàng)，點(diǎn)擊在第一頁菜單，保持默認(rèn)選項(xiàng)，點(diǎn)擊Next : 176在這一頁，我們必須定在這一頁，我們必須定義速度數(shù)據(jù)并用它來運(yùn)義速度數(shù)據(jù)并用它來運(yùn)算將偏移距數(shù)據(jù)轉(zhuǎn)成角算將偏移距數(shù)據(jù)轉(zhuǎn)成角度域數(shù)據(jù)，我們用聲波度域數(shù)據(jù)，我們用聲波曲線曲線 點(diǎn)擊點(diǎn)擊 Open Well Log 按鈕：按鈕：177選擇在角度轉(zhuǎn)換中用到的速度（井速度曲線）選擇在角度轉(zhuǎn)換中用到的速度（井速度曲線）如上菜單填寫，然后點(diǎn)擊如上菜單填寫，然后點(diǎn)擊Ok.178之前的菜單如下圖所示：之前的菜單如下圖所示：點(diǎn)擊點(diǎn)擊 Next and OK 生成結(jié)果生成結(jié)果179任意時(shí)刻的角度范圍是輸入速度域的函數(shù)，這個(gè)運(yùn)算結(jié)果顯示目的層

52、任意時(shí)刻的角度范圍是輸入速度域的函數(shù)，這個(gè)運(yùn)算結(jié)果顯示目的層630ms處角度范圍大約到處角度范圍大約到30度。這對(duì)接下來的分析足夠了，我們將要用這個(gè)角度道集體進(jìn)行接下來的聯(lián)合反演。現(xiàn)在可度。這對(duì)接下來的分析足夠了，我們將要用這個(gè)角度道集體進(jìn)行接下來的聯(lián)合反演。現(xiàn)在可以關(guān)掉它了，點(diǎn)擊以關(guān)掉它了，點(diǎn)擊 File/Exit Window 。180現(xiàn)在，我們將要生成現(xiàn)在，我們將要生成AVO屬性并且觀察改變參數(shù)后的效果，這個(gè)練習(xí)用的數(shù)據(jù)是超道集：屬性并且觀察改變參數(shù)后的效果，這個(gè)練習(xí)用的數(shù)據(jù)是超道集：181首先，首先， 我們要拾取異常的同相軸。點(diǎn)擊我們要拾取異常的同相軸。點(diǎn)擊 Horizon / Pic

53、k Horizons:在第一個(gè)菜單，點(diǎn)擊在第一個(gè)菜單，點(diǎn)擊Ok 接受默認(rèn)的名字接受默認(rèn)的名字修改如下兩個(gè)參數(shù)，它允許你拾取修改如下兩個(gè)參數(shù)，它允許你拾取道集上的波谷。道集上的波谷。182183184點(diǎn)擊點(diǎn)擊 OK 會(huì)看到沿會(huì)看到沿層振幅變化會(huì)顯示層振幅變化會(huì)顯示在窗口下方在窗口下方要?jiǎng)h除振幅顯示，點(diǎn)擊要?jiǎng)h除振幅顯示，點(diǎn)擊AVO Analysis / View Pick Analysis，把勾去掉就可以了。，把勾去掉就可以了。185要對(duì)這個(gè)數(shù)據(jù)生成截距和梯度體，在包要對(duì)這個(gè)數(shù)據(jù)生成截距和梯度體，在包含超道集的窗口選擇含超道集的窗口選擇 AVO Attribute / AVO Attribute

54、Volume接受第一頁的默認(rèn)選項(xiàng)，點(diǎn)擊接受第一頁的默認(rèn)選項(xiàng)，點(diǎn)擊 Next :現(xiàn)在我們要生成現(xiàn)在我們要生成AVO屬性體屬性體186看一下默認(rèn)菜單選項(xiàng)看一下默認(rèn)菜單選項(xiàng)我們只用我們只用AKI方程方程的前兩項(xiàng)。的前兩項(xiàng)。用井速度曲線用井速度曲線點(diǎn)擊點(diǎn)擊 Next 后后OK 開始運(yùn)算開始運(yùn)算187注意在剖面中間有一個(gè)很強(qiáng)的注意在剖面中間有一個(gè)很強(qiáng)的“紅色紅色”異常，按色標(biāo)顯示，紅色對(duì)應(yīng)強(qiáng)正振幅?？瓷戏郊t框異常，按色標(biāo)顯示，紅色對(duì)應(yīng)強(qiáng)正振幅?？瓷戏郊t框可知顯示的是截距梯度的乘積屬性可知顯示的是截距梯度的乘積屬性A*B，它的強(qiáng)正值對(duì)應(yīng)異常。，它的強(qiáng)正值對(duì)應(yīng)異常。AVO屬性運(yùn)算的結(jié)果如下屬性運(yùn)算的結(jié)果如下

55、:188實(shí)際上，顯示的層位指示氣砂巖層的頂，臨時(shí)不顯示這個(gè)層位。實(shí)際上，顯示的層位指示氣砂巖層的頂，臨時(shí)不顯示這個(gè)層位。 點(diǎn)擊眼睛：點(diǎn)擊眼睛：在參數(shù)菜單，選擇層位在參數(shù)菜單，選擇層位Horizons頁面：頁面：然后，選擇選項(xiàng)然后，選擇選項(xiàng)Do Not Show, 點(diǎn)擊點(diǎn)擊 Ok:189現(xiàn)在，我們可以清楚的看到在氣砂巖層頂和氣砂巖層底的紅色（正）強(qiáng)反射，這是我們所期望現(xiàn)在，我們可以清楚的看到在氣砂巖層頂和氣砂巖層底的紅色（正）強(qiáng)反射，這是我們所期望的三類的三類AVO異常的特點(diǎn)：異常的特點(diǎn)：190另外，要看非常有用的比例化泊松比屬性另外，要看非常有用的比例化泊松比屬性A+B，點(diǎn)擊眼睛：，點(diǎn)擊眼睛：

56、在在General菜單頁，改變菜單頁，改變 Color Data Volume 項(xiàng)如圖，點(diǎn)擊項(xiàng)如圖，點(diǎn)擊Ok:191結(jié)果如圖：結(jié)果如圖：再一次看到了氣濤巖的期望響應(yīng)再一次看到了氣濤巖的期望響應(yīng)在頂部對(duì)應(yīng)負(fù)（橘黃）變化，底部對(duì)應(yīng)正（黃色）變化。在頂部對(duì)應(yīng)負(fù)（橘黃）變化，底部對(duì)應(yīng)正（黃色）變化。SAMPLE IMAGE193Rutherford/Williams 分類 194Shale ShaleSandAcoustic Impedance = VPClass 1Class 2Class 3分類只與目的層砂巖與圍巖（泥巖）的聲波阻抗差異有關(guān)：195氣砂巖頂對(duì)應(yīng)的基本的AVO曲線196Rutherf

57、ord and Williams (1989)(a) 數(shù)據(jù)實(shí)例數(shù)據(jù)實(shí)例.(b) 模型實(shí)例模型實(shí)例197 (a) Class 2 sand.(b) Class 3 sand.Rutherford and Williams (1989)198Ross 與與 Kinman (1995) 建議創(chuàng)建一個(gè)近道疊加建議創(chuàng)建一個(gè)近道疊加(NTS)與遠(yuǎn)道疊加與遠(yuǎn)道疊加 (FTS). Class 2p: 最終疊加最終疊加 = FTS - NTS Class 2: 最終疊加最終疊加 = FTS 對(duì)已知的對(duì)已知的2類與類與2P類異常，有些處理選項(xiàng)可能提高解釋：類異常，有些處理選項(xiàng)可能提高解釋：2P199Ross an

58、d Kinman (1995)(a) 2類砂巖的全道集疊加類砂巖的全道集疊加(b) 2類砂巖的遠(yuǎn)道集疊加類砂巖的遠(yuǎn)道集疊加200Ross and Kinman (1995)(a) 2p類砂巖的全道集疊加類砂巖的全道集疊加(b) 2p 類砂巖的類砂巖的FTS - NTS2019,4(1) If 0.3 and A0.1, then B -0.575 (Class 3)(2) If 0.1 and A0.3, then B 0.075 (Class 4)BA 202這個(gè)圖闡述了這個(gè)圖闡述了4類類AVO異常的更多異常的更多詳細(xì)信息。詳細(xì)信息。203VP =1.16 VS+1360 m/sec204P

59、PVV41:Gardner,VV2VVVV4VV21B2PSSS2PSpPpPVV21A如果我們假定 VP = cVS, 我們可以推得下式：2c91A54B205如果 c = 3，梯度=0，在截距與梯度的交匯圖上就是一條水平線不同的c值會(huì)生成不同的直線（“水砂”趨勢(shì)）接下來會(huì)將其顯示在截距梯度交匯圖上A491A54B0991A54B206不同的VP/VS 比值對(duì)應(yīng)的泥巖線 (Castagna and Swan, 1998).2c91A54B207208GradientIntercept“Wet” Trend 2spVVBase IIIBase IIBase II PTop IVTop IIIT

60、op IITop II PBase IVTop IBase I交匯圖顯示異常交匯圖顯示異常209Foster et al (1993)(a) Cross-plot of well log derived A and B.(b) Cross-plot of seismically derived A and B.地震數(shù)據(jù)求得的點(diǎn)是地震數(shù)據(jù)求得的點(diǎn)是 “分散分散”的原因：的原因：(1) VP/VS 比值變化比值變化(2) 子波干擾子波干擾210(b) 解釋含氣區(qū)域(a) 未解釋含氣區(qū)域211(a) 解釋前(b) 解釋后212213至今為止，我們已經(jīng)應(yīng)用相至今為止，我們已經(jīng)應(yīng)用相應(yīng)大分析時(shí)窗生成了應(yīng)大分析時(shí)