• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于貝葉斯理論的孔隙流體模量疊前AVA反演

    2015-06-27 05:54:48印興耀張廣智張世鑫
    石油物探 2015年4期
    關(guān)鍵詞:模量反演巖石

    李 超,印興耀,張廣智,劉 倩,張世鑫

    (1.中國石油大學(華東)地球科學與技術(shù)學院,山東青島266580;2.中海油研究總院,北京100028)

    基于貝葉斯理論的孔隙流體模量疊前AVA反演

    李 超1,2,印興耀1,張廣智1,劉 倩1,張世鑫2

    (1.中國石油大學(華東)地球科學與技術(shù)學院,山東青島266580;2.中海油研究總院,北京100028)

    孔隙流體模量是一種極其敏感的流體指示因子,對該參數(shù)進行地震尺度的估算在含油氣儲層識別中具有重要的意義。首先在孔隙彈性介質(zhì)理論的指導下,推導了包含孔隙流體模量的Zoeppritz方程的線性近似方程;然后在貝葉斯理論框架下,假設(shè)先驗分布服從四元柯西分布、似然函數(shù)服從高斯分布,建立了包含正則化約束的疊前AVA反演方法。模型試算結(jié)果表明,利用該方法估算孔隙流體模量,能夠提高流體識別的精度;實際資料的試應(yīng)用取得了較好效果,證明了該方法的有效性和實用性。

    流體識別;孔隙流體模量;疊前AVA反演;貝葉斯理論

    基于地震資料的儲層流體識別技術(shù)始于20世紀70年代,利用高于周圍振幅的強反射振幅尋找地下含氣儲層的“亮點”技術(shù)得到廣泛應(yīng)用[1-2]。1987年,Smith等[3]率先提出利用縱、橫波速度相對變化量的加權(quán)疊加來判別儲層含烴異常,并給出了流體因子的概念。Goodway等[4]提出了使用拉梅參數(shù)和密度進行流體指示的LMR方法。Gray等[5]消除了密度項影響,直接以拉梅參數(shù)為流體因子進行儲層含流體識別。Quakenbush等[6]提出了泊松阻抗概念,通過組合縱、橫波阻抗來區(qū)分巖性和流體。基于Biot-Gassmann理論,Russell等[7-8]指出流體項f可以直接作為一項流體因子參與流體檢測。近些年,疊前AVO/AVA反演技術(shù)的進步極大地促進了地震流體識別技術(shù)的發(fā)展[9-11]。印興耀等[12-13]提出利用疊前反演方法直接反演Gassmann流體項,進行精確的流體識別。宗兆云等[14]建立了利用縱、橫波模量表示的孔隙彈性理論和反射系數(shù)近似方程,并發(fā)展了基于縱、橫波模量反演的流體識別方法。

    目前,常用的流體因子在流體識別的過程中會存在一定的誤差和假象,主要是受巖石固體部分的影響[15-16]。因此要提高流體識別的精度,需要一個不受巖石固體部分影響的流體因子??紫读黧w模量Kf是表征孔隙中充填流體的抗壓縮能力的物理量,它獨立于巖石骨架,不受固體部分的影響,因而在指示流體性質(zhì)方面有著獨特的優(yōu)勢[17]。本文在孔隙彈性介質(zhì)理論的基礎(chǔ)上,提出了包含孔隙流體模量的Zoeppritz方程的線性近似方程;然后利用該近似方程,基于貝葉斯理論框架構(gòu)建了疊前反演方法,形成了地震尺度的孔隙流體模量估算體系。

    1 流體因子的提出

    1.1 常規(guī)流體因子敏感性分析

    基于地震資料的流體識別通常是利用流體因子對儲層孔隙內(nèi)的流體性質(zhì)進行檢測。流體因子的選擇是流體識別成敗的關(guān)鍵,因此需要一種能夠充分反映孔隙流體特征的流體因子。在此,我們建立一個簡單的巖石模型,在其孔隙內(nèi)充填油(氣)和水,并計算孔隙度變化在0~0.35,含水飽和度變化在0~1時的幾種常用的流體指示因子的值,計算結(jié)果如圖1和圖2所示。

    圖1是孔隙流體為油和水混合物的情況下的計算結(jié)果。圖1a至圖1d中的流體因子分別是縱波阻抗、縱橫波速度比、拉梅參數(shù)λ和Gassmann流體項。從圖1中可以看出,如果含水飽和度一定,這幾種流體因子的值都會隨孔隙度的變化而劇烈變化,也就是說,這些流體指示因子的值受孔隙度的影響,而不能與含水飽和度的值一一對應(yīng),表明這幾種常用的流體因子很容易受巖石固體部分的影響,不能很好地指示孔隙流體的性質(zhì)。

    圖2是孔隙流體為氣和水混合物的情況下的計算結(jié)果。圖2中的情況與圖1類似,幾種流體因子受孔隙度影響都較大,其中圖2d中含水飽和度較低時,Gassmann流體項的變化曲線近乎于直線,這是由于尺度差異造成的,若放大顯示低含水飽和度部分(圖2e),其變化趨勢與圖1d類似,這表明Gassmann流體項本身的流體敏感性高于其它幾種流體因子,但是仍然會受固體部分的影響。

    1.2 孔隙流體模量

    這些常用的流體因子顯然不適用于高精度的流體識別,因此要尋找受巖石固體部分影響小甚至不受巖石固體部分影響的流體因子??紫读黧w模量Kf是表征孔隙中充填流體的抗壓縮能力的彈性參數(shù),這是流體本身固有的性質(zhì),不受孔隙及巖石骨架的影響。圖3顯示的是孔隙流體模量在不同孔隙度和含水飽和度情況下的值。圖3a中的孔隙流體是油和水混合物,圖3b中的孔隙流體是氣和水混合物。比較圖3a和圖3b可以看出,在含水飽和度一定時,孔隙流體模量不隨孔隙度的變化而變化;而當含水飽和度變化時,孔隙流體模量相應(yīng)地發(fā)生變化,因而孔隙流體模量與含水飽和度一一對應(yīng),這是理想的流體指示因子的性質(zhì)。因此,若可以進行地震尺度的孔隙流體模量估計,則可以更加精確地進行儲層流體識別。

    圖1 流體因子隨孔隙度和含水飽和度變化(流體為油和水)曲線

    圖2 流體因子隨孔隙度和含水飽和度變化(流體為氣和水)曲線

    圖3 孔隙流體模量隨孔隙度和含水飽和度的變化曲線

    2 孔隙流體模量直接反演

    2.1 包含孔隙流體模量的近似方程的建立

    地下巖石實際上是由固體礦物和孔隙流體組成的多孔雙相介質(zhì)。孔隙流體的存在必然會影響多孔巖石的力學特征,從而使介質(zhì)巖石表現(xiàn)出一定的彈性異常。Gassmann在一系列假設(shè)前提下,研究了不排水情況下的巖石力學性質(zhì),推導了表征雙相介質(zhì)彈性模量參數(shù)的Gassmann方程[18]。直到現(xiàn)在,在基于地震信息估算巖石彈性模量方面,特別是針對碎屑巖儲層來說,Biot-Gassmann理論的表現(xiàn)仍然最穩(wěn)健,在實際應(yīng)用中的使用頻率也最高。Gassmann方程的表達式為:

    (1)

    式中:Ksat,μsat分別為飽和巖石的體積模量和剪切模量;Kf為孔隙流體模量;Kdry,μdry分別為干燥巖石體積模量和剪切模量;Km為巖石基質(zhì)體積模量;φ為孔隙度。Gassmann方程的體積模量方程等式右邊的第2項是流體項f,即:

    (2)

    Han等[19]研究了Biot-Gassmann理論中的孔隙流體與巖石骨架的固體效應(yīng)對巖石模量的影響,提出由于Km?Kf,所以有:

    (3)

    因此,可以得到流體項的近似表達式為:

    (4)

    我們將G(φ)看做一個整體,稱為骨架因子G(無量綱)。假設(shè)流體項f中包含的巖石骨架礦物與孔隙度綜合作用均由骨架因子G來表示(G與孔隙度φ的具體關(guān)系在此不作討論),由(2)式近似得到(3)式的過程中產(chǎn)生的誤差也包含在G中,因此公式(4)可以變?yōu)椋?/p>

    f=GKf

    (5)

    公式(4)在孔隙度較大時比較精確,而在孔隙度較小時會有一定的誤差,并且隨著孔隙度降低誤差逐漸增大;公式(5)由于不是近似式,因此不會受孔隙度大小的影響。

    Russell等基于Biot-Gassmann理論研究了多孔飽和流體多孔介質(zhì)的AVO理論,推導了新的反

    射系數(shù)近似方程:

    (6)

    式中:θ為平均角度;γdry,γsat分別為干巖、飽和巖石縱橫波速度比;Δf/f,Δμ/μ,Δρ/ρ分別為流體項、剪切模量和密度的相對變化率。

    將公式(5)代入公式(6),可得:

    (7)

    令η=Gμ,則:

    (8)

    經(jīng)過簡單的運算即可得到:

    (9)

    公式(9)就是包含孔隙流體模量的反射系數(shù)近似方程。

    2.2 近似方程精度驗證

    為了檢驗推導出的AVO近似方程的精確程度,建立4個兩層地層模型,模型參數(shù)如表1所示。4個模型的地層界面分別屬于第1類至第4類AVO類型。針對這幾個模型分別用精確的Zoeppritz方程、Russell近似方程以及公式(9)計算兩個地層界面處的縱波反射系數(shù),結(jié)果如圖4所示。

    圖4a至圖4d中的橫坐標是平均角度,縱坐標是縱波反射系數(shù)。紅色曲線是利用精確Zoeppritz方程計算的結(jié)果,粉色點線是利用Russell近似方程計算的結(jié)果,藍色點線是利用公式(9)計算的結(jié)果。從圖4中可以看出,采用公式(9)計算的精度與采用Russell近似方程計算的精度接近。對于第1,2,3類AVO模型,采用公式(9)計算的精度接近于采用精確Zoeppritz方程計算的精度;計算第4類AVO模型時,在角度較小時(小于30°),采用公式(9)計算的精度較高,隨著角度的增加,與采用精確Zoeppritz方程計算的誤差逐漸增大。總的來說,公式(9)的精度能夠滿足地震反演的需求。

    表1 模型參數(shù)

    圖4 采用不同模型計算的縱波反射系數(shù)

    2.3 測井尺度孔隙流體模量估計

    在實際應(yīng)用中,無法利用常規(guī)的測井方法直接測量孔隙流體模量Kf和骨架因子G,因此需要根據(jù)巖石物理理論進行間接計算。有多種可以用來計算孔隙流體模量的方法,本文選取基于巖石物理建模的計算方法。在Gassmann方程中,若可以準確的計算出飽和巖石體積模量、干燥巖石體積模量以及巖石基質(zhì)的體積模量等參數(shù),就可以準確的計算孔隙流體模量Kf。

    計算巖石基質(zhì)的模量最常用的方法是VRH理論。Voigt模型提出的N種礦物組成的復合介質(zhì)的有效體積模量KV的計算公式為:

    (10)

    式中:KV是Voigt有效體積模量;fi是第i種組分的體積分量;Ki是第i種組分的體積模量。Voigt有效體積模量代表上邊界,利用Reuss模型給出有效體積模量的下邊界KR:

    (11)

    利用Voigt-Reuss-Hill平均得到巖石基質(zhì)的體積模量:

    (12)

    干巖石體積模量Kdry和巖石基質(zhì)體積模量Km的關(guān)系可以通過建立巖石物理模型得到。巖石物理模型需要根據(jù)實際儲層特點建立,這里假設(shè)所研究的儲層為典型的砂泥巖儲層,以最為常用的Xu-White模型[20]為例進行計算。

    Keys和Xu提出了砂巖的干巖石體積模量的表達式[21]:

    (13)

    由測井資料可以獲取縱、橫波速度以及密度信息,在此基礎(chǔ)上根據(jù)方程(14)可以得到飽和巖石體積模量:

    (14)

    利用公式(12)至公式(14)得到飽和巖石體積模量Ksat,干燥巖石體積模量Kdry和巖石基質(zhì)體積模量Km后,結(jié)合公式(1)即可求得測井曲線對應(yīng)的孔隙流體模量曲線;然后,將孔隙流體模量帶入公式(5) 即可得到骨架因子G。

    (15)

    (16)

    2.4 基于貝葉斯理論的AVA反演

    疊前AVA反演需要綜合地質(zhì)、測井和地震資料,它以包含豐富地下信息的疊前地震反射資料為主要資料,以地質(zhì)和測井資料作為約束,來揭示地下儲層的屬性及其含流體特征。為了更好地利用待反演參數(shù)的先驗信息和提高反演結(jié)果的合理性,本文采用基于貝葉斯理論的AVA反演方法。

    由于地震數(shù)據(jù)是4個不同角度的部分角度疊加道集,因此公式(9)按照入射角不同,以矩陣形式表示為:

    (17)

    考慮n個界面的情況,將矩陣進行塊化處理并加入子波矩陣,得到:

    (18)

    (18)式可簡記為:

    Gm=d

    (19)

    基于貝葉斯理論得到的反射系數(shù)的后驗概率密度函數(shù)可表示為:

    (20)

    式中:d表示部分角度疊加地震數(shù)據(jù);m表示待反演角度反射系數(shù);P(m|d)表示后驗概率密度函數(shù);P(d|m)表示描述數(shù)據(jù)d和反射系數(shù)m之間關(guān)系的似然函數(shù);P(m)是先驗概率密度函數(shù);P(d)表示邊緣分布,在只考慮后驗分布的形狀時,其數(shù)值可以取常數(shù)。

    假設(shè)地震資料背景噪聲服從高斯分布,地震數(shù)據(jù)和反射系數(shù)之間的關(guān)系可以用服從高斯分布的似然函數(shù)來描述,即

    (21)

    由于柯西分布能夠突出小反射系數(shù),能提高反演結(jié)果的分辨率,因此,假設(shè)先驗信息服從柯西分布。考慮到各待反演參數(shù)之間的相關(guān)性,這里采用四元柯西分布:

    (22)

    式中:N為地震數(shù)據(jù)樣點數(shù);Φi=(Di)TΨ-1Di;Ψ是一個4×4的系數(shù)矩陣,Ψ需要使用最大期望(Expectation Maximization,EM)算法計算;Di是4×4N的矩陣,定義為:

    (23)

    由貝葉斯公式可以得到后驗概率密度分布為:

    (24)

    將(24)式代入邊緣化公式,取對數(shù)后得到最大化后驗概率分布目標函數(shù)為:

    J(m)=JG(m)+JCauchy(m)

    (25)

    (25)式通過柯西先驗分布加入稀疏約束,在壓制噪聲影響的同時,提高了反演結(jié)果的稀疏性。但是,由于地震資料是帶限的,所以反演結(jié)果也是帶限的。為了改善反演剖面的橫向連續(xù)性,提高反演結(jié)果的預(yù)測精度,通過加入平滑正則約束項和點約束項,對反演結(jié)果低頻分量進行補償。

    為了對反演結(jié)果低頻分量進行補償,在目標函數(shù)中加入平滑正則約束項和點約束項,得到

    (26)

    式中:平滑正則約束項JSmooth(m)=α‖C1m-ξ1‖2;點約束項JPoint(m)=β‖C2m-ξ2‖2。其中,α,β為加權(quán)系數(shù),控制約束信息的相對使用量。可以看出,平滑正則約束項和點約束項的數(shù)學表達式相同,但物理意義卻不同,主要體現(xiàn)在積分矩陣C的構(gòu)建上。在平滑正則約束項中,C1和ξ1是由反射系數(shù)序列時間采樣點的數(shù)目和低頻模型或非常平滑單井模型計算得到;而在點約束項中,C2和ξ2是由約束點的位置和對應(yīng)位置處測井數(shù)據(jù)或井插值數(shù)據(jù)計算得到。

    因此,(26)式即為最終的目標函數(shù),對目標函數(shù)求梯度并令其等于0,得到反演方程:

    (27)

    式中:λ為加權(quán)系數(shù)。

    從(27)式所示的反演方程可以看出,該方程有一定的弱非線性性,因此采用重加權(quán)迭代最小二乘算法進行求解即可得到各項的相對變化量,再進行積分即可得到最終結(jié)果。

    反演流程如圖5所示。

    圖5 基于貝葉斯理論的AVA反演流程

    3 模型試算

    為驗證本文方法的可行性和穩(wěn)定性,用某實際工區(qū)的井數(shù)據(jù)建立模型進行試算。利用實際測井資料計算得到的模型曲線如圖6所示。采用主頻為35Hz的雷克子波和精確Zoeppritz方程進行正演模擬得未加噪聲的合成地震記錄(圖7)。

    圖6 用于模型試算的井數(shù)據(jù)

    圖7為未加噪聲的合成地震記錄。通過合理的角度分析進行疊加得到部分角度疊加道集,然后利用本文提出的反演方法實現(xiàn)孔隙流體模量反演。圖8為無噪聲合成數(shù)據(jù)反演結(jié)果,圖中實線表示模型的值,虛線表示反演結(jié)果。反演結(jié)果表明,無噪聲情況下,該方法獲取的孔隙流體模量與真實值吻合程度較高。

    圖7 未加噪聲的合成地震記錄

    圖8 未加噪聲合成地震記錄的反演結(jié)果

    為了進一步檢驗該方法的穩(wěn)定性,對合成地震記錄中加入信噪比為2∶1的高斯隨機噪聲,圖9是加入噪聲后的合成地震記錄。同樣,經(jīng)過角度分析進行疊加得到部分角度疊加道集并進行孔隙流體模量反演。圖10是在地震數(shù)據(jù)中加入高斯隨機噪聲(信噪比為2∶1)后的反演結(jié)果。從圖10中可以看出,孔隙流體模量反演結(jié)果仍與模型值基本吻合,但密度項的誤差相比其它3項略大,這是因為密度項的系數(shù)權(quán)值低,因而對反射系數(shù)的影響相對較小;其它兩個參數(shù)的反演結(jié)果也較為準確,從而證明該方法具有一定的抗噪性。

    圖9 加入信噪比為2∶1的高斯隨機噪聲后的合成地震記錄

    圖10 加入信噪比為2∶1的高斯隨機噪聲后的反演結(jié)果

    4 實例應(yīng)用

    將本文所述方法應(yīng)用于中國東部某油田的實際資料。在進行疊前地震反演之前,需要對地震數(shù)據(jù)進行保幅處理,包括精細的波前擴散補償、震源組合與檢波器組合效應(yīng)的校正、反Q濾波、地表一致性處理、疊前去噪處理、去除多次波等,并假設(shè)處理后的層間多次波、各向異性的影響可以忽略不計。由于輸入地震數(shù)據(jù)是部分角度疊加道集,所以要根據(jù)已知的速度信息,將疊前CMP或CRP道集轉(zhuǎn)換到角度域,再根據(jù)實際入射角范圍,進行分角度疊加。圖11是由小到大的4個角度部分疊加剖面。圖中,在CDP55處有一口井(A井),紫色直線是A井的井軌跡,井軌跡處的柱狀體表示測井解釋成果,其中,藍色部分表示非儲層,綠色部分表示含水砂巖,紅色部分表示含油砂巖;深度為1.56s處有一套含水儲層,1.58s處有一套含油儲層。首先利用常規(guī)反演方法反演得到Gassmann流體項f(圖12a),從圖12a中可以看出,在兩個儲層位置均有低值異常顯示,但是含油儲層和含水儲層的異常值一樣,無法區(qū)分。然后利用本文提出的反演方法實現(xiàn)孔隙流體模量反演Kf(圖12b),從圖12b中可以看出,在兩個儲層位置均有低值異常顯示,而且含油儲層的值比含水儲層的值更低,可以進行區(qū)分。該反演結(jié)果與巖石物理分析結(jié)果以及實際鉆井結(jié)果吻合,從而證明了方法的實用性和優(yōu)越性。

    圖11 不同中心角度的疊加道集剖面

    圖12 疊前AVA反演結(jié)果剖面

    5 結(jié)論與討論

    儲層流體識別是地震勘探的重要目標之一,孔隙流體模量作為獨立的流體指示因子,不受巖石固體部分的影響,流體敏感性高。本文在孔隙彈性介質(zhì)理論的基礎(chǔ)上,提出利用孔隙流體模量作為流體因子,并推導了包含孔隙流體模量的Zoeppritz方程的線性近似方程。在貝葉斯理論的框架下,以柯西分布作為先驗分布,高斯分布作為似然函數(shù),建立了相應(yīng)的疊前AVA反演方法。利用該方法進行疊前反演可以直接估算孔隙流體模量,提高儲層流體識別的精度。模型試算可以看出該方法精度高,而且有一定的抗噪性。將其應(yīng)用于實際工區(qū)資料反演,流體識別效果很好,說明了該方法應(yīng)用于實際生產(chǎn)的可行性。

    [1] Backus M M,Chen R L.Flat spot exploration[J].Geophysical Prospecting,1975,23(3):533-577

    [2] Ostrander W J.Plane-wave reflection coefficients for gas sands at nonnormal angles of incidence[J].Geophysics,1984,49(10):1637-1648

    [3] Smith G C,Gidlow P M.Weighted stacking for rock property estimation and detection of gas[J].Geophysical Prospecting,1987,35(9):993-1014

    [4] Goodway B,Chen T,Downton J.Improved AVO fluid detection and lithology discrimination using Lamé petrophysical parameters: “λρ”,“μρ”,“λμfluid stack”,from P and S inversions[J].Expanded Abstracts of 67thAnnual Internat SEG Mtg,1997,183-186

    [5] Gray D,Goodway B,Chen T.Bridging the gap:using AVO to detect changes in fundamental elastic constants[J].Expanded Abstracts of 69thAnnual Internat SEG Mtg,1999,852-855

    [6] Quakenbush M,Shang B,Tuttle C.Poisson impedance[J].The Leading Edge,2006,25(2):128-138

    [7] Russell B H,Hedlin K,Hilterman F J,et al.Fluid-property discrimination with AVO:a Biot-Gassmann perspective[J].Geophysics,2003,68(1):29-39

    [8] Russell B H,Gray D,Hampson D P.Linearized AVO and poroelasticity[J].Geophysics,2011,76(3):C19-C29

    [9] Buland A,Omre H.Bayesian linearized AVO inversion[J].Geophysics,2003,68(1):185-198

    [10] Downton J E.Seismic parameter estimation from AVO inversion[D].Calgary:University of Calgary,2005

    [11] Yin X Y,Yang P J,Zhang G Z.A novel prestack AVO inversion and its application[J].Expanded Abstracts of 78thAnnual Internat SEG Mtg,2008,2041-2045

    [12] 印興耀,張世鑫,張繁昌,等.利用基于Russell近似的彈性波阻抗反演進行儲層描述和流體識別[J].石油地球物理勘探,2010,45(3):373-380 Yin X Y,Zhang S X,Zhang F C,et al.Utilizing Russell approximation based elastic wave impedance inversion to conduct reservoir description and fluid identification[J].Oil Geophysical Prospecting,2010,45(3):373-380

    [13] 印興耀,宗兆云,吳國忱.非均勻介質(zhì)孔隙流體參數(shù)

    地震散射波反演[J].中國科學:地球科學,2013,43(12):1934-1942 Yin X Y,Zong Z Y,Wu G C.Seismic wave scattering inversion for fluid factor of heterogeneous media[J].Science in China:Earth Sciences,2013,43(12):1934-1942

    [14] 宗兆云,印興耀,吳國忱.基于疊前地震縱橫波模量直接反演的流體檢測方法[J].地球物理學報,2012,55(1):284-292 Zong Z Y,Yin X Y,Wu G C.Fluid identification method based on compressional and shear modulus direct inversion[J].Chinese Journal of Geophysics,2012,55(1):284-292

    [15] 張世鑫.基于地震信息的流體識別方法研究與應(yīng)用[D].青島:中國石油大學(華東),2012 Zhang S X.Methodology and application of fluid identification with seismic information[D].Qingdao:China University of Petroleum,2012

    [16] 印興耀,曹丹平,王保麗,等.基于疊前地震反演的流體識別方法研究進展[J].石油地球物理勘探,2014,49(1):22-34 Yin X Y,Cao D P,Wang B L,et al.Research progress of fluid discrimination with pre-stack seismic inversion[J].Oil Geophysical Prospecting,2014,49(1):22-34

    [17] 印興耀,李超,張世鑫.基于雙相介質(zhì)的地震流體識別[J].中國石油大學學報(自然科學版),2013,37(5):38-43 Yin X Y,Li C,Zhang S X.Seismic fluid discrimination based on two-phase theory[J].Journal of China University of Petroleum (Edition of Natural Sciences),2013,37(5):38-43

    [18] Gassmann F.Uber die elastizitat,poroser medien[J].Vierteljahrsschrift der Natruforschenden Gesellschaft in Zurich,1951,96:1-23

    [19] Han D H,Batzle M.Gassmann’s equation and fluid-saturation effects on seismic velocities[J].Geophysics,2004,69(2):398-405

    [20] Xu S,White R E.A new velocity model for clay-sand mixtures[J].Geophysical Prospecting,1995,43(1):91-118

    [21] Keys R G,Xu S Y.An approximation for the Xu-White velocity model[J].Geophysics,2002,67(5):1406-1414

    [22] Eshelby J D.The determination of the elastic field of an ellipsoidal inclusion and related problems[J].Proceedings of the Royal Society of London,1957,241(1226):376-396

    (編輯:顧石慶)

    Prestack AVA inversion for pore fluid modulus based on the Bayesian theory

    Li Chao1,2,Yin Xingyao1,Zhang Guangzhi1,Liu Qian1,Zhang Shixin2

    (1.CollegeofGeosciences,ChinaUniversityofPetroleum,Qingdao266580,China; 2.CNOOCResearchInstitute,Beijing100028,China)

    Pore fluid modulus is an extremely sensitive fluid indicator.It is significant to estimate pore fluid modulus for reservoir fluid discrimination in seismic scale.Based on the poroelastic theory,a Zoeppritz approximation equation containing the pore fluid modulus is derived.Then in the theoretical framework of Bayesian theory,we assume that the prior distribution can be modelled via the four-variable Cauchy distribution and the likelihood function follows the Gaussian distribution.The objective function of prestack AVA inversion is constructed containing the regularized constraint.A synthetic data test shows that the estimation of pore fluid modulus with this method improves the accuracy of fluid identification,and the application results of actual seismic data demonstrates its effectiveness and practicability.

    fluid discrimination,pore fluid modulus,prestack AVA inversion,Bayesian theory

    2014-11-12;改回日期:2015-01-29。

    李超(1987—),男,博士,主要從事儲層預(yù)測和流體識別方面的研究工作。

    國家重點基礎(chǔ)研究發(fā)展計劃(973計劃)項目(2013CB228604,2014CB239201)、國家油氣重大專項(2011ZX05014-001-010HZ)、中國石油科技創(chuàng)新基金項目(2011D-5006-0301)、中央高?;究蒲袠I(yè)務(wù)費專項資金項目(11CX06003A)共同資助。

    P631

    A

    1000-1441(2015)04-0467-10

    10.3969/j.issn.1000-1441.2015.04.014

    猜你喜歡
    模量反演巖石
    反演對稱變換在解決平面幾何問題中的應(yīng)用
    第五章 巖石小專家
    3深源巖石
    一種叫做煤炭的巖石
    高勁度模量瀝青混合料在京臺高速車轍維修段的應(yīng)用
    室內(nèi)回彈模量和回彈再壓縮模量試驗參數(shù)探討
    山西建筑(2020年11期)2020-06-04 00:09:48
    海藻與巖石之間
    關(guān)于現(xiàn)行規(guī)范路基頂面回彈模量的理解和應(yīng)用
    上海公路(2018年4期)2018-03-21 05:57:24
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    欧美精品高潮呻吟av久久| 美女高潮到喷水免费观看| 亚洲精品中文字幕在线视频| 久久国产精品影院| 国产精品电影一区二区三区 | 操出白浆在线播放| 国产色视频综合| 色精品久久人妻99蜜桃| 高清黄色对白视频在线免费看| 欧美激情高清一区二区三区| 90打野战视频偷拍视频| 亚洲av日韩精品久久久久久密| 国产免费现黄频在线看| 美女高潮到喷水免费观看| 欧美激情极品国产一区二区三区| 免费观看人在逋| 色94色欧美一区二区| 中文字幕另类日韩欧美亚洲嫩草| 国产男女超爽视频在线观看| 一夜夜www| 久久中文字幕一级| 日韩视频一区二区在线观看| 国产亚洲精品久久久久5区| 亚洲欧美色中文字幕在线| 日本五十路高清| 亚洲精华国产精华精| 久久久精品94久久精品| 精品一区二区三卡| 999久久久国产精品视频| 一区二区三区激情视频| 女人被躁到高潮嗷嗷叫费观| 国产真人三级小视频在线观看| 一区二区日韩欧美中文字幕| 热99久久久久精品小说推荐| av视频免费观看在线观看| 国产免费av片在线观看野外av| 欧美黄色片欧美黄色片| 深夜精品福利| 成在线人永久免费视频| 肉色欧美久久久久久久蜜桃| 亚洲精品中文字幕一二三四区 | 精品欧美一区二区三区在线| √禁漫天堂资源中文www| 久久久精品94久久精品| 男男h啪啪无遮挡| 欧美性长视频在线观看| 久久久精品区二区三区| 亚洲全国av大片| 日韩视频在线欧美| 在线永久观看黄色视频| 精品少妇黑人巨大在线播放| 99国产精品免费福利视频| 女性生殖器流出的白浆| 满18在线观看网站| 一边摸一边抽搐一进一小说 | 日韩有码中文字幕| 99久久精品国产亚洲精品| 欧美一级毛片孕妇| 成人18禁高潮啪啪吃奶动态图| 精品国产乱子伦一区二区三区| 一级毛片精品| 国产亚洲精品一区二区www | 日韩免费av在线播放| 亚洲欧美激情在线| 在线十欧美十亚洲十日本专区| 国产视频一区二区在线看| 久久天躁狠狠躁夜夜2o2o| 免费女性裸体啪啪无遮挡网站| 五月开心婷婷网| 国产伦人伦偷精品视频| 最新的欧美精品一区二区| 国产av又大| 女人精品久久久久毛片| www.熟女人妻精品国产| 一级毛片女人18水好多| 国产精品秋霞免费鲁丝片| 精品欧美一区二区三区在线| 熟女少妇亚洲综合色aaa.| 国产一区二区 视频在线| 亚洲中文日韩欧美视频| 激情在线观看视频在线高清 | 欧美精品亚洲一区二区| 国产精品1区2区在线观看. | av一本久久久久| 怎么达到女性高潮| 久久久久网色| 另类精品久久| tube8黄色片| 他把我摸到了高潮在线观看 | 久久久久久久久免费视频了| 日日爽夜夜爽网站| 久久精品亚洲精品国产色婷小说| 无限看片的www在线观看| 欧美日韩精品网址| av免费在线观看网站| 亚洲精品中文字幕在线视频| 嫩草影视91久久| 国产精品麻豆人妻色哟哟久久| kizo精华| 亚洲精品粉嫩美女一区| 亚洲中文av在线| 一夜夜www| 91国产中文字幕| 精品久久蜜臀av无| 久久国产精品人妻蜜桃| 午夜日韩欧美国产| 一本大道久久a久久精品| 日韩一区二区三区影片| 18禁观看日本| 在线观看一区二区三区激情| 热99久久久久精品小说推荐| 黄色视频不卡| 欧美在线一区亚洲| 欧美日韩一级在线毛片| 精品午夜福利视频在线观看一区 | 亚洲欧美日韩另类电影网站| 天堂俺去俺来也www色官网| 99久久人妻综合| 在线观看人妻少妇| 一个人免费看片子| 久久精品aⅴ一区二区三区四区| av视频免费观看在线观看| 人人妻人人澡人人看| 天天操日日干夜夜撸| 亚洲一卡2卡3卡4卡5卡精品中文| 成人国产av品久久久| 在线十欧美十亚洲十日本专区| 美女国产高潮福利片在线看| 久久人人97超碰香蕉20202| 80岁老熟妇乱子伦牲交| 日本精品一区二区三区蜜桃| 99精国产麻豆久久婷婷| 男女无遮挡免费网站观看| 亚洲久久久国产精品| 欧美精品人与动牲交sv欧美| 日本wwww免费看| 欧美久久黑人一区二区| 免费日韩欧美在线观看| 午夜两性在线视频| 日本五十路高清| 美女高潮到喷水免费观看| 一级片免费观看大全| aaaaa片日本免费| 看免费av毛片| 精品少妇一区二区三区视频日本电影| 成人av一区二区三区在线看| 国产精品久久久av美女十八| 五月开心婷婷网| 在线观看免费日韩欧美大片| 高清av免费在线| 国产av国产精品国产| 一级a爱视频在线免费观看| 亚洲国产毛片av蜜桃av| videosex国产| 国产精品 国内视频| 久久精品国产99精品国产亚洲性色 | 欧美激情久久久久久爽电影 | 欧美激情 高清一区二区三区| 狠狠婷婷综合久久久久久88av| 69精品国产乱码久久久| 午夜成年电影在线免费观看| 亚洲色图综合在线观看| 大片电影免费在线观看免费| 欧美久久黑人一区二区| 久久国产亚洲av麻豆专区| 女人爽到高潮嗷嗷叫在线视频| 蜜桃国产av成人99| 精品一品国产午夜福利视频| 人人妻人人澡人人爽人人夜夜| 夜夜爽天天搞| 99久久精品国产亚洲精品| 一级毛片精品| 亚洲精品成人av观看孕妇| 女人高潮潮喷娇喘18禁视频| 亚洲国产看品久久| 亚洲精品国产精品久久久不卡| 老鸭窝网址在线观看| 久久久久久亚洲精品国产蜜桃av| 精品卡一卡二卡四卡免费| 大型黄色视频在线免费观看| 国产高清激情床上av| 黄色视频不卡| 在线观看免费日韩欧美大片| 久久久久久亚洲精品国产蜜桃av| 精品高清国产在线一区| 中文字幕高清在线视频| 久久人人97超碰香蕉20202| 久久毛片免费看一区二区三区| 最近最新中文字幕大全免费视频| 日韩有码中文字幕| 50天的宝宝边吃奶边哭怎么回事| 一边摸一边抽搐一进一出视频| 天天躁夜夜躁狠狠躁躁| 亚洲成人国产一区在线观看| a级毛片在线看网站| 欧美成人免费av一区二区三区 | www.熟女人妻精品国产| 中文字幕制服av| 国产伦人伦偷精品视频| 国产成人免费无遮挡视频| 亚洲第一av免费看| 精品国产一区二区久久| 热re99久久精品国产66热6| 香蕉丝袜av| 制服人妻中文乱码| 久久精品熟女亚洲av麻豆精品| 午夜精品久久久久久毛片777| 757午夜福利合集在线观看| 午夜福利,免费看| 国产区一区二久久| 成人av一区二区三区在线看| 天天操日日干夜夜撸| 日日摸夜夜添夜夜添小说| 99re在线观看精品视频| 老司机午夜十八禁免费视频| 天堂俺去俺来也www色官网| 精品国产乱码久久久久久男人| 精品第一国产精品| 亚洲精品久久成人aⅴ小说| 18禁黄网站禁片午夜丰满| 国产精品99久久99久久久不卡| 欧美日本中文国产一区发布| 国产激情久久老熟女| 久久久久久人人人人人| 又大又爽又粗| a级毛片黄视频| 成人免费观看视频高清| 人人妻人人添人人爽欧美一区卜| 精品人妻在线不人妻| 国产极品粉嫩免费观看在线| 妹子高潮喷水视频| 国产成人精品久久二区二区91| av国产精品久久久久影院| av欧美777| 亚洲国产av影院在线观看| 亚洲欧美激情在线| 一级片免费观看大全| 99久久99久久久精品蜜桃| 宅男免费午夜| 免费av中文字幕在线| 国产在线一区二区三区精| 这个男人来自地球电影免费观看| 岛国在线观看网站| 老司机影院毛片| 国产日韩一区二区三区精品不卡| 在线观看免费视频网站a站| 久久久精品国产亚洲av高清涩受| 丰满饥渴人妻一区二区三| 亚洲七黄色美女视频| 免费一级毛片在线播放高清视频 | 窝窝影院91人妻| 欧美黑人精品巨大| 日韩中文字幕视频在线看片| 动漫黄色视频在线观看| 电影成人av| 精品卡一卡二卡四卡免费| 蜜桃在线观看..| 中文字幕色久视频| videos熟女内射| 亚洲天堂av无毛| 一二三四在线观看免费中文在| 国产一区二区三区在线臀色熟女 | 久久久国产一区二区| 亚洲国产成人一精品久久久| 最新的欧美精品一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久亚洲精品国产蜜桃av| 丝瓜视频免费看黄片| 久久亚洲真实| 中文字幕精品免费在线观看视频| 免费高清在线观看日韩| 精品高清国产在线一区| 国产三级黄色录像| 国精品久久久久久国模美| 亚洲视频免费观看视频| 国产精品欧美亚洲77777| 12—13女人毛片做爰片一| 亚洲一区二区三区欧美精品| 欧美精品亚洲一区二区| 久久精品国产亚洲av高清一级| 少妇猛男粗大的猛烈进出视频| 免费在线观看完整版高清| 国产男女内射视频| 少妇精品久久久久久久| 99久久精品国产亚洲精品| 丝袜在线中文字幕| 波多野结衣一区麻豆| 水蜜桃什么品种好| 母亲3免费完整高清在线观看| 黑人操中国人逼视频| 亚洲熟妇熟女久久| 国产人伦9x9x在线观看| 国产免费av片在线观看野外av| 成年女人毛片免费观看观看9 | 搡老岳熟女国产| 精品国产乱码久久久久久小说| 成人特级黄色片久久久久久久 | 一本综合久久免费| 日韩熟女老妇一区二区性免费视频| 两个人看的免费小视频| 久久亚洲精品不卡| 免费久久久久久久精品成人欧美视频| 婷婷成人精品国产| 亚洲天堂av无毛| 久久精品亚洲熟妇少妇任你| 精品久久久久久久毛片微露脸| xxxhd国产人妻xxx| 亚洲精品美女久久av网站| 最近最新中文字幕大全电影3 | 1024香蕉在线观看| 国产高清国产精品国产三级| 国产一卡二卡三卡精品| 黄色丝袜av网址大全| 男女无遮挡免费网站观看| 亚洲成a人片在线一区二区| 999久久久国产精品视频| 后天国语完整版免费观看| 97人妻天天添夜夜摸| 肉色欧美久久久久久久蜜桃| 咕卡用的链子| 久久中文字幕一级| 亚洲av成人一区二区三| 亚洲午夜理论影院| 免费在线观看完整版高清| 在线播放国产精品三级| 精品国产一区二区三区四区第35| 欧美日韩亚洲国产一区二区在线观看 | 在线观看免费视频日本深夜| 俄罗斯特黄特色一大片| 成人精品一区二区免费| 女人久久www免费人成看片| 久久99热这里只频精品6学生| 久热这里只有精品99| 可以免费在线观看a视频的电影网站| 午夜免费成人在线视频| 午夜福利在线观看吧| 亚洲色图av天堂| 日韩三级视频一区二区三区| xxxhd国产人妻xxx| 午夜福利,免费看| 丰满人妻熟妇乱又伦精品不卡| 国产精品香港三级国产av潘金莲| 高清av免费在线| 巨乳人妻的诱惑在线观看| 国产精品亚洲av一区麻豆| 日本五十路高清| 国产一区二区在线观看av| 自线自在国产av| 女警被强在线播放| 9191精品国产免费久久| 黄色视频,在线免费观看| √禁漫天堂资源中文www| 精品国内亚洲2022精品成人 | 每晚都被弄得嗷嗷叫到高潮| 99香蕉大伊视频| 精品国产超薄肉色丝袜足j| 老鸭窝网址在线观看| 黄色视频在线播放观看不卡| 亚洲av国产av综合av卡| 妹子高潮喷水视频| 在线观看免费午夜福利视频| 黄色 视频免费看| 免费女性裸体啪啪无遮挡网站| 久久国产亚洲av麻豆专区| 国产主播在线观看一区二区| 国产精品欧美亚洲77777| 人妻 亚洲 视频| 一级毛片电影观看| 黄片小视频在线播放| 欧美国产精品一级二级三级| 亚洲第一青青草原| 日本av手机在线免费观看| 国产在线视频一区二区| 国产成人精品久久二区二区免费| 丝袜美腿诱惑在线| 在线 av 中文字幕| 久久精品91无色码中文字幕| 国产亚洲一区二区精品| 亚洲人成电影免费在线| 亚洲精品在线观看二区| 老司机午夜福利在线观看视频 | 国产国语露脸激情在线看| 欧美+亚洲+日韩+国产| 精品久久久久久电影网| 色在线成人网| 国产精品免费大片| 黄片小视频在线播放| 国产在线精品亚洲第一网站| 国产一区二区三区综合在线观看| 精品一区二区三区av网在线观看 | 丰满人妻熟妇乱又伦精品不卡| 热99久久久久精品小说推荐| 午夜免费鲁丝| 露出奶头的视频| 咕卡用的链子| 久久久久久久精品吃奶| 日韩一区二区三区影片| 亚洲一卡2卡3卡4卡5卡精品中文| 精品人妻1区二区| aaaaa片日本免费| 波多野结衣av一区二区av| 18禁裸乳无遮挡动漫免费视频| 国产一区有黄有色的免费视频| 老司机福利观看| 亚洲国产看品久久| 91大片在线观看| 精品乱码久久久久久99久播| 日日爽夜夜爽网站| 香蕉丝袜av| 亚洲一区中文字幕在线| 欧美日韩成人在线一区二区| 日韩免费av在线播放| 精品久久久久久电影网| 国产真人三级小视频在线观看| 欧美中文综合在线视频| 欧美黄色片欧美黄色片| 真人做人爱边吃奶动态| 叶爱在线成人免费视频播放| 757午夜福利合集在线观看| 久久精品国产亚洲av高清一级| 老司机在亚洲福利影院| 中文字幕人妻丝袜制服| 久热爱精品视频在线9| 国产精品偷伦视频观看了| 久久久久久久大尺度免费视频| 一区二区av电影网| 成人国产一区最新在线观看| 777米奇影视久久| 午夜福利乱码中文字幕| 亚洲精品在线观看二区| 成人国产av品久久久| 一本—道久久a久久精品蜜桃钙片| 搡老岳熟女国产| 免费av中文字幕在线| 深夜精品福利| 精品国产国语对白av| 中文字幕人妻丝袜制服| 又紧又爽又黄一区二区| 老鸭窝网址在线观看| 国产日韩欧美视频二区| 99热国产这里只有精品6| 涩涩av久久男人的天堂| 免费av中文字幕在线| 亚洲久久久国产精品| 国产成人精品久久二区二区免费| 久久99一区二区三区| 亚洲五月色婷婷综合| 欧美成狂野欧美在线观看| 日韩制服丝袜自拍偷拍| 欧美日韩亚洲国产一区二区在线观看 | 日韩欧美国产一区二区入口| 大片电影免费在线观看免费| 日本精品一区二区三区蜜桃| 99久久99久久久精品蜜桃| 狠狠狠狠99中文字幕| 国产视频一区二区在线看| av福利片在线| 国产成人系列免费观看| 欧美精品av麻豆av| 视频区欧美日本亚洲| 高清欧美精品videossex| 亚洲人成77777在线视频| 欧美国产精品va在线观看不卡| 国产黄色免费在线视频| 欧美乱妇无乱码| 成人av一区二区三区在线看| 国产单亲对白刺激| 另类亚洲欧美激情| 91精品国产国语对白视频| 18禁美女被吸乳视频| 国产精品亚洲一级av第二区| 三级毛片av免费| 美女午夜性视频免费| 国产日韩一区二区三区精品不卡| 国产精品一区二区在线观看99| 欧美国产精品va在线观看不卡| 99国产综合亚洲精品| 中文字幕精品免费在线观看视频| 肉色欧美久久久久久久蜜桃| 在线观看免费高清a一片| 日韩欧美免费精品| 亚洲专区国产一区二区| 人人妻人人澡人人看| 精品国产一区二区三区四区第35| 国产成人av激情在线播放| 久久国产精品大桥未久av| 97人妻天天添夜夜摸| 久久香蕉激情| 国产又爽黄色视频| 日本av手机在线免费观看| 日韩欧美免费精品| 成人特级黄色片久久久久久久 | 欧美 日韩 精品 国产| 国产精品亚洲av一区麻豆| 国产不卡av网站在线观看| 啦啦啦视频在线资源免费观看| 狠狠狠狠99中文字幕| 国产精品1区2区在线观看. | 国产精品秋霞免费鲁丝片| 日韩精品免费视频一区二区三区| 亚洲色图综合在线观看| 男女下面插进去视频免费观看| 中文亚洲av片在线观看爽 | 99久久99久久久精品蜜桃| 老鸭窝网址在线观看| 另类精品久久| 成年女人毛片免费观看观看9 | 老司机影院毛片| 九色亚洲精品在线播放| av福利片在线| av欧美777| 丁香六月欧美| 久久久久久久久久久久大奶| 欧美日韩亚洲高清精品| 熟女少妇亚洲综合色aaa.| 亚洲熟妇熟女久久| 一边摸一边抽搐一进一出视频| 亚洲天堂av无毛| 天天影视国产精品| 黄色怎么调成土黄色| 亚洲av第一区精品v没综合| 精品国产亚洲在线| 久久精品国产99精品国产亚洲性色 | 两性夫妻黄色片| 新久久久久国产一级毛片| av天堂久久9| 极品教师在线免费播放| 亚洲成人免费电影在线观看| 十分钟在线观看高清视频www| 亚洲,欧美精品.| 成人永久免费在线观看视频 | 男女午夜视频在线观看| 搡老乐熟女国产| 99精国产麻豆久久婷婷| 天堂俺去俺来也www色官网| 成人国语在线视频| 99国产精品免费福利视频| 国产精品熟女久久久久浪| 免费观看a级毛片全部| 国产xxxxx性猛交| 日日夜夜操网爽| 母亲3免费完整高清在线观看| 亚洲 国产 在线| www.精华液| 日韩熟女老妇一区二区性免费视频| 777久久人妻少妇嫩草av网站| 国产极品粉嫩免费观看在线| 欧美变态另类bdsm刘玥| 咕卡用的链子| 1024视频免费在线观看| 国产高清视频在线播放一区| 亚洲欧洲日产国产| 黑丝袜美女国产一区| 在线观看免费视频日本深夜| 一本—道久久a久久精品蜜桃钙片| 日本精品一区二区三区蜜桃| 国产老妇伦熟女老妇高清| 成人手机av| 啦啦啦免费观看视频1| 啦啦啦视频在线资源免费观看| a级毛片在线看网站| 纯流量卡能插随身wifi吗| 免费av中文字幕在线| 大陆偷拍与自拍| 色老头精品视频在线观看| 99热国产这里只有精品6| 少妇被粗大的猛进出69影院| 色综合婷婷激情| 亚洲av国产av综合av卡| 丝袜喷水一区| 亚洲色图综合在线观看| 色精品久久人妻99蜜桃| 国产男靠女视频免费网站| 久久久久久久久免费视频了| 大型av网站在线播放| 丝瓜视频免费看黄片| 最近最新中文字幕大全电影3 | 99香蕉大伊视频| 天天影视国产精品| 少妇猛男粗大的猛烈进出视频| 大型av网站在线播放| 国产无遮挡羞羞视频在线观看| av网站在线播放免费| 亚洲熟妇熟女久久| 手机成人av网站| 热99re8久久精品国产| 精品乱码久久久久久99久播| 在线 av 中文字幕| 国产麻豆69| 欧美黑人欧美精品刺激| 美女国产高潮福利片在线看| 精品久久久久久电影网| 中文字幕人妻丝袜一区二区| 手机成人av网站| e午夜精品久久久久久久| 国产一区二区三区综合在线观看| 黄色成人免费大全| 首页视频小说图片口味搜索| 欧美激情极品国产一区二区三区| 侵犯人妻中文字幕一二三四区| 国产单亲对白刺激| 女性生殖器流出的白浆| 天堂俺去俺来也www色官网| 高清av免费在线| 国产av精品麻豆| 国产精品久久久久久精品古装| 免费久久久久久久精品成人欧美视频| 国产男女超爽视频在线观看| 欧美亚洲日本最大视频资源| 一级毛片电影观看| 18禁美女被吸乳视频| 日韩制服丝袜自拍偷拍| 91老司机精品| 久久天躁狠狠躁夜夜2o2o| 欧美av亚洲av综合av国产av|