多波聯(lián)合反演技術(shù)在SC盆地PLN地區(qū)油氣勘探開發(fā)中的應(yīng)用研究

王 棟,賀振華,王 瓏,楊海濤,趙 堯,符志國

(1.中國石油天然氣集團(tuán)公司川慶鉆探工程有限公司地球物理勘探公司,四川成都610213;2.中國石油天然氣集團(tuán)公司山地地震技術(shù)試驗(yàn)基地,四川成都610213;3.成都理工大學(xué)地球物理學(xué)院,四川成都610059)

王棟,賀振華,王瓏,等.多波聯(lián)合反演技術(shù)在SC盆地PLN地區(qū)油氣勘探開發(fā)中的應(yīng)用研究[J].石油物探,2017,56(6):-873

WANG Dong,HE Zhenhua,WANG Long,et al.Application of multi-wave joint inversion in exploration and development of oil and gas in the PLN area of SC Basin[J].Geophysical Prospecting for Petroleum,2017,56(6):-873

多波聯(lián)合反演技術(shù)在SC盆地PLN地區(qū)油氣勘探開發(fā)中的應(yīng)用研究

王 棟1,2,3,賀振華3,王 瓏1,2,楊海濤1,2,趙 堯1,2,符志國1,2

(1.中國石油天然氣集團(tuán)公司川慶鉆探工程有限公司地球物理勘探公司,四川成都610213;2.中國石油天然氣集團(tuán)公司山地地震技術(shù)試驗(yàn)基地,四川成都610213;3.成都理工大學(xué)地球物理學(xué)院,四川成都610059)

多波聯(lián)合反演是多波地震勘探的重要環(huán)節(jié)之一。從多波聯(lián)合反演的基本原理出發(fā),研究了該方法的特點(diǎn)及技術(shù)流程,并將該技術(shù)應(yīng)用于SC盆地PLN地區(qū)某工區(qū)多波資料處理。首先對參與反演的子波提取方法進(jìn)行了對比與選擇;然后利用近、中、遠(yuǎn)3個角度域的子波進(jìn)行連井線和全工區(qū)數(shù)據(jù)反演;最后利用該地區(qū)速度與孔隙度、速度與含氣飽和度擬合公式得到孔隙度和含氣飽和度預(yù)測結(jié)果。對多波聯(lián)合反演與單一縱波反演在有效儲層預(yù)測方面的實(shí)際應(yīng)用效果進(jìn)行了對比,結(jié)果表明,多波聯(lián)合反演可以增加對反演結(jié)果的約束條件,降低單一縱波反演的多解性,提高反演精度,對儲層的識別更加客觀可靠,可以有效提高油氣儲層預(yù)測的準(zhǔn)確性。

多波地震勘探;子波提取;多波聯(lián)合反演;速度;孔隙度;含氣飽和度

隨著地震采集儀器精度的提高和地震數(shù)據(jù)處理技術(shù)的進(jìn)步,多波多分量地震勘探技術(shù)得到了迅速發(fā)展和廣泛應(yīng)用。多波聯(lián)合反演是多波地震勘探的重要環(huán)節(jié)之一。STEWART[1]將SMITH等[2]提出的PP波加權(quán)疊加方法推廣到聯(lián)合反演,首次給出了實(shí)用的縱橫波聯(lián)合反演方法,并得到了縱、橫波速度比等彈性參數(shù)。LARSEN等[3]在STEWART[1]研究的基礎(chǔ)上討論了同時反演縱、橫波阻抗的方法。楊紹國等[4]、李錄明等[5]、雍楊等[6]、王明春等[7-8]、黃中玉等[9]對聯(lián)合反演公式進(jìn)行了研究,進(jìn)一步利用多波信息聯(lián)合反演出巖性參數(shù)及更多的彈性參數(shù),據(jù)此對巖性進(jìn)行有效識別。但是,由于在常規(guī)反演過程中增加了橫波參數(shù),反演算法更加復(fù)雜,對橫波資料品質(zhì)的要求更高,因此,學(xué)者們不斷研究推出新的思路、優(yōu)化算法及高質(zhì)量橫波資料處理方法。VEIRE等[10]討論了利用奇異值分解技術(shù)(SVD)進(jìn)行縱橫波聯(lián)合三參數(shù)反演的方法,給出了該方法在實(shí)際數(shù)據(jù)處理中的應(yīng)用效果。陳天勝等[11]提出了一種基于方向加速度最優(yōu)化和縱橫波速度比值掃描的縱橫波聯(lián)合反演方法。張春濤等[12]對縱橫波聯(lián)合反演方法的研究及應(yīng)用進(jìn)行了總結(jié)。HU等[13]利用縱波和轉(zhuǎn)換波聯(lián)合反演來估算密度比和速度比,并在反演過程中根據(jù)貝葉斯理論加入先驗(yàn)信息,以改善反演的不適定性。侯棟甲等[14]進(jìn)一步研究了基于貝葉斯理論的VTI介質(zhì)疊前多波聯(lián)合反演技術(shù)。杜啟振等[15]提出了基于橫縱波速度比值迭代的縱橫波聯(lián)合反演方法。張廣智等[16]提出了縱橫波聯(lián)合疊前自適應(yīng)MCMC反演方法。黃中玉等[17]提出正交各向異性介質(zhì)中的多方位三維轉(zhuǎn)換波疊前時間偏移方法。尋超等[18]提出并利用多方向矢量中值濾波方法進(jìn)行多分量地震數(shù)據(jù)的去噪處理。隨著研究的不斷深入,多波聯(lián)合反演技術(shù)被廣泛應(yīng)用于各大勘探工區(qū),取得了良好的效果。如ZHANG等[19]對Pikes Peak油田3C-2D地震數(shù)據(jù)體進(jìn)行了縱橫波聯(lián)合反演;王興建等[20]利用縱橫波聯(lián)合反演技術(shù)對蘇里格含氣性進(jìn)行了檢測;程冰潔等[21]、黃躍等[22]對川西地區(qū)3C-3D地震數(shù)據(jù)體進(jìn)行了聯(lián)合反演;付雷等[23]利用縱波和轉(zhuǎn)換波聯(lián)合反演技術(shù)對松遼盆地北部地區(qū)油氣藏分布進(jìn)行了預(yù)測;李昂等[24]將縱、橫波聯(lián)合反演方法應(yīng)用于海拉爾盆地貝39井區(qū)儲層描述。

SC盆地PLN地區(qū)須家河組主要為陸相地層,沉積環(huán)境復(fù)雜,埋藏時間較長,砂巖儲層致密化程度高,在測井曲線上呈“三低”特征(低孔隙度、低飽和度、低電阻率)。目前常用的地震反演預(yù)測技術(shù)雖然可以部分解決該區(qū)有效儲層的厚度及橫向分布預(yù)測問題,但因反演技術(shù)的局限性以及該區(qū)儲層的非均質(zhì)性較強(qiáng),預(yù)測成果存在多解性。同時,利用伽馬反演技術(shù)進(jìn)行巖性預(yù)測的結(jié)果存在較多爭議。

常規(guī)反演過程中,縱波參數(shù)來源于實(shí)際縱波資料,而橫波參數(shù)是由縱波參數(shù)和該地區(qū)的經(jīng)驗(yàn)關(guān)系公式計(jì)算得到的,并不是實(shí)測的結(jié)果,勢必與縱波存在一定的相關(guān)性。因此,常規(guī)反演并沒有真正同時利用PP波與PS波的實(shí)測信息,從而忽略了縱橫波阻抗、縱橫波速度和地質(zhì)體密度之間內(nèi)在的聯(lián)系。事實(shí)上,縱波速度與橫波速度之間存在著線性關(guān)系(如Castagna泥巖公式[25]),縱波速度與地質(zhì)體密度之間也存在著密切聯(lián)系(如Gardner公式[26])。忽略縱波阻抗、橫波阻抗與密度之間的內(nèi)在聯(lián)系,將直接降低反演結(jié)果的精度和可靠性,最終影響儲層預(yù)測、巖性判斷、流體識別的準(zhǔn)確性。本文通過多波聯(lián)合反演技術(shù)的應(yīng)用,獲得SC盆地PLN地區(qū)縱波速度、橫波速度、密度、縱橫波速度比、泊松比以及拉梅常數(shù)等巖石物理參數(shù),為該區(qū)儲層和流體預(yù)測提供更為可靠的依據(jù)。

1 方法原理

多波聯(lián)合反演仍然沿用縱波反演框架,只是在反演目標(biāo)函數(shù)上增加了轉(zhuǎn)換橫波誤差項(xiàng),可以用(1)式表示:

(1)

式中:m表示待求解的模型目標(biāo)參數(shù),J(m)表示目標(biāo)函數(shù);SPP表示PP波合成記錄,DPP表示PP波地震記錄;E(SPP-DPP)表示PP波誤差能量函數(shù)項(xiàng),E(SPS-DPS)為PS波誤差能量函數(shù)項(xiàng)。極小化目標(biāo)函數(shù)即求解一組模型的參數(shù)使得合成記錄與地震記錄的誤差為最小。

Zoeppritz方程描述了精確的反射系數(shù)與地層彈性參數(shù)的關(guān)系[27],但該方程為隱式方程,不利于求解。FATTI等[28]對縱波與轉(zhuǎn)換波反射系數(shù)公式進(jìn)行了簡化,使縱波、轉(zhuǎn)換波反射系數(shù)成為與速度、密度、入射角有關(guān)的函數(shù)?；隈薹e模型,利用反射系數(shù)與地震子波合成地震記錄,則反射振幅成為入射角的函數(shù),即:

為增強(qiáng)反演解的穩(wěn)定性,利用以下對數(shù)線性關(guān)系進(jìn)行約束:

對于多個入射角的道集,可將PP波、PS波地震記錄表示為:

(6)

將公式(6)進(jìn)一步表示成:

(7)

2 技術(shù)流程

多波聯(lián)合反演的技術(shù)流程如圖1所示?？梢钥?/p>

出,多波聯(lián)合反演需要輸入6類數(shù)據(jù),分別為:縱波角度道集數(shù)據(jù)、轉(zhuǎn)換波角度道集數(shù)據(jù)、縱波子波、轉(zhuǎn)換波子波、初始縱橫波速度及密度模型、測井速度及密度曲線。其中,初始縱橫波速度模型可以利用縱、橫波疊前時間偏移速度建立,初始密度模型利用密度測井曲線建立。測井曲線在反演過程中用于井控處理。

多波聯(lián)合反演的關(guān)鍵是子波的求取。子波形態(tài)的不同,會直接影響反演的結(jié)果。目前常用的子波求取方法有理論計(jì)算子波、地震資料提取子波和井震聯(lián)合提取子波,需要將不同子波的反演結(jié)果與測井結(jié)果進(jìn)行對比,選擇適合工區(qū)實(shí)際的子波提取方法。其次是利用目的層段的測井曲線計(jì)算擬合對數(shù)線性關(guān)系的系數(shù)k,m。最后是進(jìn)行連井?dāng)?shù)據(jù)試驗(yàn)反演,判斷反演結(jié)果的準(zhǔn)確性,完成整個工區(qū)數(shù)據(jù)體的反演計(jì)算。

圖1 多波聯(lián)合反演技術(shù)流程

3 實(shí)際應(yīng)用效果分析

3.1 多波聯(lián)合反演技術(shù)在PLN地區(qū)的應(yīng)用

在進(jìn)行多波聯(lián)合反演之前,通常需要對縱波、轉(zhuǎn)換波數(shù)據(jù)進(jìn)行層位標(biāo)定、層位對比解釋、層位匹配以及相應(yīng)的巖石物理參數(shù)分析。PLN地區(qū)影響儲層預(yù)測精度的主要因素是難以區(qū)分低速泥巖和低速含氣砂巖。由于轉(zhuǎn)換橫波主要反映巖性變化,流體對其影響很小,因此,利用轉(zhuǎn)換波的特點(diǎn)并結(jié)合縱波資料,可以有效地進(jìn)行巖性識別,去除泥巖的影響。

由于PLN地區(qū)須家河組儲層物性橫向變化較大,理論子波不適合該地區(qū)數(shù)據(jù)反演,因此需要對子波提取的方法進(jìn)行比較和分析。具體做法是將提取的子波和利用測井曲線計(jì)算得到的波阻抗褶積,再計(jì)算與實(shí)際地震道的相關(guān)系數(shù),如果相關(guān)系數(shù)高,說明子波合理,反之則不合理。

圖2對比了井震結(jié)合方法與地震道統(tǒng)計(jì)方法提取子波的效果。其中圖2a是利用井震結(jié)合方法提取的子波制作的過41井、206井合成地震記錄(藍(lán)色)與實(shí)際地震記錄(紅色)對比結(jié)果,41井的相關(guān)系數(shù)為0.82,206井的相關(guān)系數(shù)為0.73;圖2b是采用地震道統(tǒng)計(jì)方法得到的子波制作的過41井、206井合成地震記錄(藍(lán)色)與實(shí)際地震記錄(紅色)對比結(jié)果,41井的相關(guān)系數(shù)為0.64,206井的相關(guān)系數(shù)為0.56?？梢钥闯?井震結(jié)合方法提取子波整體相關(guān)系數(shù)較高,故該方法更適合本工區(qū)數(shù)據(jù)反演。

在多波聯(lián)合反演過程中,需要考慮反演的時效性。通常的做法是將偏移后道集分成近、中、遠(yuǎn)3個角度域,每個角度域提取一個子波,并通過目的層AVO正演分析確定最大角度的取值(圖3)。由圖3可以看出,在偏移距4300m左右,轉(zhuǎn)換波反射系數(shù)接近于0,在5000m以后,縱波反射系數(shù)出現(xiàn)畸變,因此最大偏移距選為4200m,此時的入射角為42°。

圖2 不同子波提取方法比較(藍(lán)色為合成地震記錄,紅色為實(shí)際地震記錄) a 井震結(jié)合方法; b 地震道統(tǒng)計(jì)方法

圖3 縱橫波AVO反射系數(shù)分析

將偏移后的縱波、轉(zhuǎn)換波道集抽取為角度道集,按照14°的間隔分為3個角度域,采用井震結(jié)合全局尋優(yōu)方法提取每個角度域的子波,如圖4和圖5所示?？梢钥闯?縱波子波主頻為45Hz左右,轉(zhuǎn)換波子波主頻為25Hz左右。利用提取的子波對PLN地區(qū)三維三分量資料連井線進(jìn)行多波聯(lián)合反演試驗(yàn),再利用得到的阻抗和密度計(jì)算出相應(yīng)的縱波速度和橫波速度,如圖6和圖7所示?？v波和橫波速度反演剖面均層次清晰,形態(tài)自然,與研究區(qū)河流相沉積砂包泥的地質(zhì)特征和測井資料吻合。

將縱、橫波速度反演結(jié)果與實(shí)際測井曲線進(jìn)行對比,如圖8所示。可以看出,多波聯(lián)合反演得到的橫波速度曲線較單一縱波反演與實(shí)際測井曲線的趨勢一致性更好,從而更能反映地層在縱向上的物性變化,揭示巖性和儲層的分布。

由于縱波速度受流體影響較大,因此在有利巖性和儲層分布區(qū)域內(nèi),可以利用反演得到的縱波速度擬合儲層孔隙度和含水(氣)飽和度等參數(shù)。將研究區(qū)測井孔隙度和含水飽和度與砂巖縱波速度進(jìn)行交會分析,如圖9所示,色標(biāo)為泥質(zhì)含量。通過數(shù)據(jù)擬合得到孔隙度和含水(氣)飽和度與速度的關(guān)系:

圖4 縱波(a)、轉(zhuǎn)換波(b)角度道集

圖5 縱橫波不同角度域子波及頻譜 a 縱波子波; b 轉(zhuǎn)換波子波; c 縱波子波頻譜; d 轉(zhuǎn)換波子波頻譜

圖6 PLN地區(qū)連井線縱波速度聯(lián)合反演

(6)

(7)

圖7 PLN地區(qū)連井線橫波速度聯(lián)合反演

圖8 縱、橫波速度反演結(jié)果與測井曲線對比 a 111井; b 107井

圖9 測井孔隙度(a)和含水飽和度(b)與砂巖縱波速度交會分析

圖10為利用為多波聯(lián)合反演得到的縱波速度由公式(6)計(jì)算得到的連井線孔隙度預(yù)測結(jié)果,可見優(yōu)質(zhì)儲層分布與井對應(yīng)良好。圖11為由公式(7)所示流體飽和度與縱波速度的擬合關(guān)系得到的含氣(水)飽和度預(yù)測剖面,高含氣部位與氣井對應(yīng)關(guān)系較好。儲層和流體預(yù)測結(jié)果均表明,多波聯(lián)合反演能夠較好地刻畫目標(biāo)層地質(zhì)情況。

為了提高多波聯(lián)合反演的儲層厚度預(yù)測精度,通常在地震數(shù)據(jù)基礎(chǔ)上加入測井資料對其進(jìn)行約束,即增加了高頻信息。研究區(qū)目的層主要為薄互層砂體,將多波聯(lián)合反演的結(jié)果與實(shí)際測井資料進(jìn)行對比,如表1所示?？梢钥闯?在該區(qū)已知5口井中,多波聯(lián)合反演得到的目標(biāo)層段總儲層厚度與實(shí)際鉆井情況更加接近,整體效果優(yōu)于單一縱波預(yù)測結(jié)果,預(yù)測符合率達(dá)到100%,厚度預(yù)測誤差小于5m,孔隙度預(yù)測誤差小于1.0%,含氣飽和度預(yù)測誤差小于0.5%。

圖10 連井線孔隙度預(yù)測

圖11 連井線含氣(水)飽和度預(yù)測

井名油氣產(chǎn)量/(104t)測井解釋厚度/m縱波預(yù)測厚度/m多波預(yù)測厚度/m測井解釋孔隙度，%縱波預(yù)測孔隙度，%多波預(yù)測孔隙度，%測井解釋含氣飽和度，%縱波預(yù)測含氣飽和度，%多波預(yù)測含氣飽和度，%9井氣：0．318水：36．00010．87．09．07．17．47．2455548107井氣：52．00017．020．018．08．18．58．2585560111井干井0．54．02．06．16．36．2404040128井干井10．17．511．06．56．56．5495050113井氣：19．50010．68．011．07．67．87．7615056

在上述研究的基礎(chǔ)上,對全區(qū)數(shù)據(jù)進(jìn)行反演計(jì)算,得到了縱波阻抗、橫波阻抗和密度的三維數(shù)據(jù)體,如圖12所示。

利用反演得到的阻抗和密度,計(jì)算出縱、橫波速度,再利用研究區(qū)速度與含水飽和度擬合公式(7)計(jì)算出含水飽和度。圖13為目的層含氣(水)飽和度預(yù)測平面圖,圖中含氣飽和度大于44%的紅色-黃色為含氣飽和度高的含氣有利區(qū),綠色為含氣飽和度中等的含氣較有利區(qū),藍(lán)色為含氣飽和度低的含氣不利區(qū)?？梢钥闯?大部分工業(yè)氣井都處于含氣飽和度大于44%的紅色-黃色含氣有利區(qū)。利用試驗(yàn)區(qū)塊內(nèi)5口測試井資料(未參與反演),對含氣性預(yù)測結(jié)果的準(zhǔn)確性進(jìn)行了評價,如表2所示。5口驗(yàn)證井中有4口井與鉆井吻合,符合率80%,唯有002-2-X1井測試為干井,在圖上卻處于含氣性較高的有利區(qū)。綜合前文所述中參與反演的5口井,共計(jì)10口井中有9口井與鉆井吻合,符合率90%。

圖12 縱橫波阻抗、密度多波聯(lián)合反演數(shù)據(jù)體 a 縱波阻抗; b 橫波阻抗; c 密度

圖13 PLN地區(qū)目的層預(yù)測的含氣飽和度平面分布

3.2 與單一縱波反演方法效果對比

將多波聯(lián)合反演的縱橫波阻抗與單一縱波反演的結(jié)果進(jìn)行了對比,如圖14所示。整體面貌上兩者較為一致,但是聯(lián)合反演結(jié)果的分辨率有所提高(圖14 中紅色圓圈處所示)。兩種方法反演的縱波阻抗差異相對較小,但聯(lián)合反演的橫波阻抗縱、橫向?qū)哟蚊黠@變得更加清晰,表明橫波振幅信息發(fā)揮了較大作用。

圖15對比了多波聯(lián)合反演與單一縱波反演的橫波速度剖面。圖15a中橢圓處單一縱波反演的結(jié)果為橫波速度降低,意味該處存在巖性變化,而此處多波聯(lián)合反演結(jié)果則相反,橫波速度變化不大,說明不存在巖性變化(圖15b)。若基于圖15a所示反演結(jié)果進(jìn)行孔隙度計(jì)算,橢圓處目標(biāo)體將作為非儲層被去除,如圖16a所示。若基于圖15b所示反演結(jié)果進(jìn)行孔隙度計(jì)算,則橢圓處目標(biāo)體將作為高孔儲層被保留,如圖16b所示。與實(shí)際鉆井結(jié)果進(jìn)行對比(該井為日產(chǎn)52×104m3的高產(chǎn)氣井)可知,多波聯(lián)合反演結(jié)果所表現(xiàn)出的高孔儲層特征更加符合實(shí)際情況。同時,單一縱波反演的結(jié)果有大量泥巖沒有被去除;而多波聯(lián)合反演則很好地去除了泥巖等非儲層的影響,有效地提高了儲層預(yù)測的精度。

表2 PLN地區(qū)含氣性預(yù)測結(jié)果準(zhǔn)確性評價

圖14 多波聯(lián)合反演與單一縱波反演縱、橫波阻抗對比 a 單一縱波反演縱波阻抗; b 多波聯(lián)合反演縱波阻抗; c 單一縱波反演橫波阻抗; d 多波聯(lián)合反演橫波阻抗

圖15 多波聯(lián)合反演與單一縱波反演的橫波速度剖面對比 a 單一縱波反演; b 多波聯(lián)合反演

圖16 多波聯(lián)合反演與單一縱波反演的孔隙度剖面對比 a 單一縱波反演; b 多波聯(lián)合反演

4 結(jié)束語

多波聯(lián)合反演由于利用實(shí)際轉(zhuǎn)換橫波資料增加了對反演結(jié)果的約束條件,因此,能獲得較高精度的反演參數(shù),去除泥巖等非儲層的影響,顯著提高了反演結(jié)果的可靠性和精度。多波聯(lián)合反演技術(shù)在SC盆地PLN地區(qū)的應(yīng)用有效地降低了單一縱波反演預(yù)測的多解性,較好地刻畫了目的層的地質(zhì)情況,為須家河組致密砂巖氣藏的預(yù)測提供了可靠的依據(jù)。

[1] STEWART R R.Joint P and P-SV inversion[R].Calgary:The CREWES Project,1990:112-115

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

[3] LARSEN J,MARGRAVE G,LU H X.AVO analysis by simultaneous P-P and P-S weighted stacking applied to 3C-3D seismic data[J].Expanded Abstracts of 69thAnnual Internat SEG Mtg,1999:721-724

[4] 楊紹國,周熙襄.多波多層 AVO 數(shù)據(jù)反演[J].石油地球物理勘探,1994,29(6):695-705

YANG S G,ZHOU X X.Multiwave multilayer AVO inversion[J].Oil Geophysical Prospecting,1994,29(6):695-705

[5] 李錄明,羅省賢,王明春,等.各向異性介質(zhì)三維縱橫波聯(lián)合疊前反演方法及應(yīng)用[J].石油地球物理勘探,2010,45(1):60-65

LI L M,LUO X X,WANG M C,et al.3D PP-PS joint inversion method and application in anisotropic medium[J].Oil Geophysical Prospecting,2010,45(1):60-65

[6] 雍楊,李錄明,羅省賢,等.TI介質(zhì)多波AVA方程及參數(shù)反演[J].石油物探,2004,43(1):11-16

YONG Y,LI L M,LUO X X,et al.Multiwave AVA equation and parameter inversion in TI media[J].Geophysical Prospecting for Petroleum,2004,43(1):11-16

[7] 王明春,李錄明,羅省賢.TI介質(zhì)多波AVA聯(lián)合反演巖性參數(shù)及應(yīng)用[J].天然氣工業(yè),2007,27(4):49-52

WANG M C,LI L M,LUO X X.Method for lithology parameter inversion with joint multiwave AVA in transversely isotropic media and its application[J].Natural Gas Industry,2007,27(4):49-52

[8] 王明春.VTI介質(zhì)多波疊前聯(lián)合反演巖性參數(shù)方法及其應(yīng)用[D].成都:成都理工大學(xué),2007

WANG M C.Method and application of multiwave prestack joint inversion lithologic parameters in VTI media[D].Chengdu :Chengdu University of Technology,2007

[9] 黃中玉,趙金州.縱波和轉(zhuǎn)換波 AVO 聯(lián)合反演技術(shù)[J].石油物探,2004,43(4):319-322

HUANG Z Y,ZHAO J Z.Simultaneous P-and S-wave AVO inversion[J].Geophysical Prospecting for Petroleum,2004,43(4):319-322

[10] VEIRE H H,LANDRO M.Simultaneous inversion of PP and PS seismic data[J].Geophysics,2006,71(3):1-10

[11] 陳天勝,劉洋,魏修成.縱波和轉(zhuǎn)換波聯(lián)合 AVO 反演方法研究[J].中國石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2006,30(1):33-37

CHEN T S,LIU Y,WEI X C.Joint amplitude versus offset inversion of P-P and P-SV seismic data[J].Journal of China University of Petroleum,2006,30(1):33-37

[12] 張春濤,王尚旭,李生杰,等.縱橫波聯(lián)合反演方法研究及應(yīng)用[J].石油地球物理勘探,2010,45(4):520-524

ZHANG C T,WANG S X,LI S J,et al.Compressional and shear wave joint inversion technique research and its application[J].Oil Geophysical Prospecting,2010,45(4):520-524

[13] HU G Q,LIU Y.Joint PP and PS AVO inversion based on Bayes theorem[J].Applied Geophysics,2011,8(4):293-302

[14] 侯棟甲,劉洋,任志明,等.基于貝葉斯理論的VTI介質(zhì)多波疊前聯(lián)合反演[J].石油物探,2014,53(3):294-303

HOU D J,LIU Y,REN Z M,et al.Multi-wave prestack joint inversion in VTI media based on Bayesian theory[J].Geophysical Prospecting for Petroleum,2014,53(3):294-303

[15] 杜啟振,霍國棟,陳剛.基于橫縱波速度比值迭代的縱橫波聯(lián)合反演方法[J].石油地球物理勘探,2016,51(1):97-105

DU Q Z,HUO G D,CHEN G.Joint inversion of PP-and PS-waves based on the iteration of ratio of S-wave velocity and P-wave velocity[J].Oil Geophysical Prospecting,2016,51(1):97-105

[16] 張廣智,潘新朋,孫昌路,等.縱橫波聯(lián)合疊前自適應(yīng)MCMC反演方法[J].石油地球物理勘探,2016,51(5):938-946

ZHANG G Z,PAN X P,SUN C L,et al.PP-& PS-wave prestack nonlinear inversion based on adaptive MCMC algorithm[J].Oil Geophysical Prospecting,2016,51(5):938-946

[17] 黃中玉,余波,王于靜,等.正交各向異性介質(zhì)中多方位三維轉(zhuǎn)換波疊前時間偏移[J].石油物探,2015,54(4):435-442

HUANG Z Y,YU B,WANG Y J,et al.Multi-azimuth 3D converted wave prestack time migration in orthotropic medium[J].Geophysical Prospecting for Petroleum,2015,54(4):435-442

[18] 尋超,汪超,王赟.多方向矢量中值濾波在多分量地震數(shù)據(jù)中的應(yīng)用[J].石油物探,2016,55(5):703-710

XUN C,WANG C,WANG Y.The application of multi-directional vector median filtering in multi-component seismic data[J].Geophysical Prospecting for Petroleum,2016,55(5):703-710

[19] ZHANG H,MARGRAVE G F,BROWN R J.Joint PP-PS inversion at Pikes Peak oilfield,Saskatchewan[R].Calgary:CSEG Annual Convention,2003

[20] 王興建,曹俊興.縱橫波聯(lián)合反演在蘇里格含氣性檢測中的應(yīng)用[J].天然氣工業(yè),2008,28(10):44-45

WANG X J,CAO J X.Application of p-and s-wave inversion to detecting gas potential in the Sulige gas field[J].Natural Gas Industry,2008,28(10):44-45

[21] 程冰潔,徐天吉.轉(zhuǎn)換波資料在川西坳陷儲層預(yù)測中的應(yīng)用[J].石油物探,2009,48(2):181-186

CHENG B J,XU T J.Application of converted wave data in reservoir prediction for Chuanxi depression[J].Geophysical Prospecting for Petroleum,2009,48(2):181-186

[22] 黃躍,徐天吉,程冰潔,等.轉(zhuǎn)換波疊后聯(lián)合反演技術(shù)在川西須家河組氣藏的應(yīng)用[J].天然氣工業(yè),2010,30(4):38-41

HUANG Y,XU T J,CHENG B J,et al.Application of converted-wave poststack joint inversion in the Xujiahe Formation gas reservoir in the western Sichuan Basin[J].Natural Gas Industry,2010,30(4):38-41

[23] 付雷,王建民,許紅,等.縱波和轉(zhuǎn)換波聯(lián)合反演技術(shù)在松遼盆地北部地區(qū)的應(yīng)用[J].石油物探,2009,48(6):563-567

FU L,WANG J M,XU H,et al.Application of P-wave and PS-wave jointed inversion in LMD area of northern Songliao Basin[J].Geophysical Prospecting for Petroleum,2009,48(6):563-567

[24] 李昂,張麗艷,宋宗平,等.縱、橫波聯(lián)合反演方法在海拉爾盆地貝39井區(qū)的應(yīng)用[J].地球物理學(xué)進(jìn)展,2013,28(4):1908-1918

LI A,ZHANG L Y,SONG Z P,et al.The application of the method of joint inversion of PP and PS waves in Bei39 zone of Hailar Basin[J].Progress in Geophysics,2013,28(4):1908-1918

[25] CASTAGNA J P,BATZLE M L,EASTWOOD R L.Relationships between compressional-wave and shear-wave velocities in clastic rocks[J].Geophysics,1985,50(4):571-581

[26] GARDNER G H F,GARDNER L W,GERGORY A R.Formation velocity and density—the diagnostic basics for stratigraphic traps[J].Geophysics,1974,39(6):770-780

[27] ZOEPPRITZ K,ERDBEBENWELLEN Ⅷ B.On the reflection and propagation of seismic waves[J].Gottinger Nachrichten,1919,17(1):66-84

[28] FATTI J L,SMITH G C,VAIL P J,et al.Detection of gas in sandstone reservoirs using AVO analysis:a 3D seismic case history using the geostack technique[J].Geophysics,1994,59(9):1362-1376

(編輯：戴春秋)

Applicationofmulti-wavejointinversioninexplorationanddevelopmentofoilandgasinthePLNareaofSCBasin

WANG Dong1,2,3,HE Zhenhua3,WANG Long1,2,YANG Haitao1,2,ZHAO Yao1,2,FU Zhiguo1,2

(1.GeophysicalProspectingCompany,CNPCChuanqingDrillingEngineeringLimitedCompany,Chengdu610213,China;2.MountainGeophysicalTechnologyTestCenter,CNPC,Chengdu610213,China;3.CollegeofGeophysics,ChengduUniversityofTechnology,Chengdu610059,China)

Multi-wave joint inversion is an important step in multi-wave seismic exploration.In this paper,the basic principle and technical process of multi-wave joint inversion are introduced,and the technique is applied to process multi-wave data from the SC basin.In the process of application,this study initially compares and selects the wavelet extraction method.Subsequently,the inversion of the well profile and the entire work area was carried out using the obtained wavelet in the near,middle,and far angles.Then,the prediction porosity and gas saturation were estimated by the velocity and porosity,and velocity and gas saturation fitting formula,respectively.Finally,the multi-wave joint inversion and the single P-wave inversion were compared with respect to the effective reservoir prediction.The results showed that multi-wave joint inversion could increase the constraints on inversion,reduce the multi-solution of the single P-wave inversion,improve the inversion accuracy,and make the reservoir identification more objective and reliable,thus effectively improving the accuracy of oil and gas reservoir prediction and development.

multi-wave seismic exploration,wavelet extraction,multi-wave joint inversion,velocity,porosity,gas saturation

2017-02-08;改回日期2017-05-08。

王棟(1984—),男,碩士,工程師,現(xiàn)從事多波地震勘探技術(shù)研究工作。

國家科技重大專項(xiàng)(2011ZX05019-008)和國家自然科學(xué)基金項(xiàng)目(41374111)共同資助。

This research is financially supported by the National Science and Technology Major Project of China (Grant No.2011ZX05019-008) and the National Natural Science Foundation of China (Grant No.41374111).

P631

A

1000-1441(2017)06-0863-11

10.3969/j.issn.1000-1441.2017.06.012