專利名稱：利用縱波、轉(zhuǎn)換波數(shù)據(jù)聯(lián)合反演儲(chǔ)層縱橫波速度比的方法
技術(shù)領(lǐng)域：
本發(fā)明涉及地球物理勘探技術(shù)，更具體地講，涉及ー種利用縱波、轉(zhuǎn)換波數(shù)據(jù)聯(lián)合反演儲(chǔ)層縱橫波速度比參數(shù)并進(jìn)行儲(chǔ)層含氣性或含水性識(shí)別的方法。
背景技術(shù)：
在油氣檢測(cè)和油藏描述中，介質(zhì)的彈性參數(shù)十分重要，這些參數(shù)與巖性和流體成分有夫。描述彈性介質(zhì)的最主要的3個(gè)參數(shù)是縱波速度、橫波速度和密度。目前地球物理勘探的常規(guī)技術(shù)是利用縱波AVCKAmplitude Versus O ffset,振幅隨炮檢距變化)疊前道集反演上述3個(gè)參數(shù)，反演時(shí)考慮了縱波反射振幅隨入射角變化的因素。但是大量理論研究和實(shí)例應(yīng)用表明，在有限偏移距的縱波資料中，這3個(gè)參數(shù)的準(zhǔn)確求解存在一定的多解性。隨著多波多分量地震勘探技術(shù)的不斷發(fā)展，高質(zhì)量轉(zhuǎn)換波資料的獲取，使得縱波、轉(zhuǎn)換波聯(lián)合反演技術(shù)成為隨之發(fā)展起來的新技木?？v波、轉(zhuǎn)換波聯(lián)合反演方法的基本原理是基于Smith和Gidlow( 1987)提出的加權(quán)疊加方法，通過PP數(shù)據(jù)的加權(quán)疊加來計(jì)算縱波和橫波波阻抗以及其它的相關(guān)參數(shù)。Fatti等(1994)做了進(jìn)ー步的研究，取消了對(duì)于Gardner公式的依賴并且包含了密度項(xiàng)。Stewart等(2001)、Landro等(2006)將加權(quán)疊加的思想用于聯(lián)合反演，給出了實(shí)用的縱橫波聯(lián)合反演方法。但利用上述反演方法，反演結(jié)果得到的是縱波速度、橫波速度和密度的變化率AVp/Vp、AVS/VS、Δ P/P，不能直接獲取儲(chǔ)層的縱橫波速度比值，對(duì)儲(chǔ)層含不同流體成分的彈性參數(shù)特征描述不夠直觀精確，解釋復(fù)雜程度較大。

發(fā)明內(nèi)容
本發(fā)明的示例性實(shí)施例的目的在于克服現(xiàn)有技術(shù)中的問題，例如，在于克服現(xiàn)有縱波、轉(zhuǎn)換波儲(chǔ)層巖石物理弾性參數(shù)聯(lián)合反演方法存在的上述不足。因此，本發(fā)明的示例性實(shí)施例提出了ー種利用縱波、轉(zhuǎn)換波數(shù)據(jù)直接反演儲(chǔ)層縱橫波速度比的方法，從而能夠?qū)?chǔ)層不同含氣性或含水性條件下地球物理響應(yīng)特征進(jìn)行精細(xì)刻畫，并對(duì)縱波、轉(zhuǎn)換波數(shù)據(jù)聯(lián)合反演的計(jì)算方法提出了改進(jìn)方案。此外，本發(fā)明的示例性實(shí)施例采用了基于阻尼最小ニ乘解的廣義線性反演方法，提高了反演算法的穩(wěn)定性，提高了計(jì)算效率和反演精度。本發(fā)明的示例性實(shí)施例提供ー種利用縱波、轉(zhuǎn)換波數(shù)據(jù)聯(lián)合反演儲(chǔ)層縱橫波速度比的方法，所述利用縱波、轉(zhuǎn)換波數(shù)據(jù)聯(lián)合反演儲(chǔ)層縱橫波速度比的方法可以包括下述步驟。(I)輸入地震數(shù)據(jù)，從而基于相同成像點(diǎn)位置的縱波成像道集和轉(zhuǎn)換波成像道集，根據(jù)目的儲(chǔ)層所在時(shí)間范圍來選取縱波反演時(shí)窗[n]和轉(zhuǎn)換波反演時(shí)窗[m(2)根據(jù)選取的時(shí)窗的每ー個(gè)縱波成像道集和轉(zhuǎn)換波成像道集，建立縱波反射系數(shù)表達(dá)式和轉(zhuǎn)換波反射系數(shù)表達(dá)式
r n D ⑷ 1 avP Δρ f I ΔΓ 4 AVs 2 Δρ) . 2ηRpp(O) = - ^Tii + - +Sin26>
2V Vp P j 2 vp r Fs r P
權(quán)利要求
1.ー種利用縱波、轉(zhuǎn)換波數(shù)據(jù)聯(lián)合反演儲(chǔ)層縱橫波速度比的方法，其特征在于，所述方法包括下述步驟 輸入地震數(shù)據(jù)，從而基于相同成像點(diǎn)位置的縱波成像道集和轉(zhuǎn)換波成像道集，根據(jù)目的儲(chǔ)層所在時(shí)間范圍來選取縱波反演時(shí)窗[Ta，Tb]和轉(zhuǎn)換波反演時(shí)窗[Ud]； 根據(jù)選取的時(shí)窗的每ー個(gè)縱波成像道集和轉(zhuǎn)換波成像道集，建立縱波反射系數(shù)表達(dá)式和轉(zhuǎn)換波反射系數(shù)表達(dá)式
2.根據(jù)權(quán)利要求I所述的方法，其特征在于，輸入的地震數(shù)據(jù)包括與ニ維縱波和轉(zhuǎn)換波相關(guān)的數(shù)據(jù)以及與三維縱波和轉(zhuǎn)換波相關(guān)的數(shù)據(jù)中的至少ー種。
3.根據(jù)權(quán)利要求2所述的方法，其特征在干， 所述成像點(diǎn)位置是輸入的地震數(shù)據(jù)中的與縱波和轉(zhuǎn)換波相關(guān)的數(shù)據(jù)體中的任意ー個(gè)相同的共深度點(diǎn)或共轉(zhuǎn)換點(diǎn)位置，或者 所述成像點(diǎn)位置是輸入的地震數(shù)據(jù)中的與縱波和轉(zhuǎn)換波相關(guān)的數(shù)據(jù)體中的一段共深度點(diǎn)或共轉(zhuǎn)換點(diǎn)成像范圍，或者 所述成像點(diǎn)位置是輸入的地震數(shù)據(jù)中的全部數(shù)據(jù)體的所有成像點(diǎn)范圍。
4.根據(jù)權(quán)利要求I所述的方法，其特征在干， 根據(jù)地震エ區(qū)目標(biāo)儲(chǔ)層的地質(zhì)層位來進(jìn)行選取時(shí)窗的步驟。
5.根據(jù)權(quán)利要求I所述的方法，其特征在于，根據(jù)輸入的地震數(shù)據(jù)中的與縱波成像道集和轉(zhuǎn)換波成像道集相關(guān)的振幅隨炮檢距變化信息來計(jì)算縱波成像道集在選定時(shí)窗內(nèi)給定樣點(diǎn)的截距、縱波成像道集在選定時(shí)窗內(nèi)給定樣點(diǎn)的梯度以及轉(zhuǎn)換波成像道集在選定時(shí)窗內(nèi)給定樣點(diǎn)的梯度。
6.根據(jù)權(quán)利要求I至權(quán)利要求5中的任意一項(xiàng)權(quán)利要求所述的方法，其特征在于，所述方法還包括下述步驟 采集包括ニ維多分量和三維多分量中的至少ー種的地震數(shù)據(jù)資料； 通過ニ維或三維縱波、轉(zhuǎn)換波地震資料處理，以得到包括能夠進(jìn)行儲(chǔ)層縱橫波速度比反演的縱波、轉(zhuǎn)換波疊前時(shí)間偏移道集數(shù)據(jù)的將被輸入的地震數(shù)據(jù)。
7.根據(jù)權(quán)利要求I所述的方法，其特征在于，所述方法還包括下述步驟 分別對(duì)于其他的成像點(diǎn)位置進(jìn)行反演以得到相應(yīng)的縱橫波速度比，從而得到整個(gè)地震數(shù)據(jù)的縱橫波速度比反演結(jié)果； 根據(jù)得到的整個(gè)地震數(shù)據(jù)的縱橫波速度比反演結(jié)果，進(jìn)行儲(chǔ)層預(yù)測(cè)和氣水識(shí)別。
全文摘要
本發(fā)明提供一種利用縱波、轉(zhuǎn)換波數(shù)據(jù)聯(lián)合反演儲(chǔ)層縱橫波速度比的方法，所述方法采用道集地震反射信息獲取的γ變化率與旅行時(shí)信息獲取的γ變化率無限逼近這一思想，構(gòu)建反演目標(biāo)函數(shù)，目標(biāo)函數(shù)的關(guān)鍵參數(shù)包括縱波截距、縱波梯度、轉(zhuǎn)換波梯度及旅行時(shí)信息計(jì)算得到的縱橫波速度比，這些參數(shù)均可由實(shí)際輸入的縱波、轉(zhuǎn)換波道集的AVO信息和旅行時(shí)信息獲取，因此目標(biāo)函數(shù)的構(gòu)建更加合理，最終反演的縱橫波速度比彈性參數(shù)也更加精確。
文檔編號(hào)G01V1/28GK102692645SQ20121017791
公開日2012年9月26日 申請(qǐng)日期2012年6月1日 優(yōu)先權(quán)日2012年6月1日
發(fā)明者何光明, 周阿波, 李忠, 楊海濤, 王棟, 謝萬學(xué) 申請(qǐng)人:中國(guó)石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司