劉寶國

(中國石油化工股份有限公司華北分公司，河南鄭州450006)

大牛地氣田位于鄂爾多斯盆地北部，氣層主要集中分布在上古生界下石盒子組、山西組及太原組，儲層具有非常強(qiáng)烈的非均質(zhì)性[1-2]。對于低滲特低滲氣田來說，水平井完鉆井技術(shù)加分段壓裂技術(shù)是動用儲量的最有效手段[3-4]。因此，2012年后大牛地氣田整體轉(zhuǎn)為水平井規(guī)?；_發(fā)[5]。在水平井部署過程中，能否準(zhǔn)確地對單砂體進(jìn)行預(yù)測與雕刻，并得到儲層在地下的準(zhǔn)確埋藏深度是水平井軌跡設(shè)計的根本依據(jù)。由于地震數(shù)據(jù)分辨率較低，測井與地震數(shù)據(jù)之間的時深標(biāo)定不能夠達(dá)到絕對準(zhǔn)確，人為誤差是不可消除的。當(dāng)研究區(qū)內(nèi)一定數(shù)量的井因井徑擴(kuò)張導(dǎo)致合成記錄不能反映真實(shí)地層信息時，井間會存在不同程度的標(biāo)定誤差，反演的預(yù)測能力及深度預(yù)測的準(zhǔn)確性受到很大的制約[6-7]。除此之外，由于水平井鉆井過程中，鉆頭方位角會發(fā)生近90°轉(zhuǎn)變，也就是說水平井A靶點(diǎn)數(shù)據(jù)僅能測到目的層著陸點(diǎn)處，其下伏地層并未鉆穿，測井?dāng)?shù)據(jù)的不完整，致使目的層段無法形成有效的合成記錄，從而難以進(jìn)行時深標(biāo)定，水平井導(dǎo)眼或A靶點(diǎn)測井?dāng)?shù)據(jù)得不到有效利用。

針對上述問題，我們通過攻關(guān)研發(fā)與應(yīng)用試驗(yàn)，形成了適用于大牛地氣田水平井軌跡設(shè)計的高精度深度域隨機(jī)模擬反演技術(shù)，使得水平井導(dǎo)眼及A靶點(diǎn)測井?dāng)?shù)據(jù)無需進(jìn)行合成記錄標(biāo)定，直接在深度域開展隨機(jī)模擬反演，大幅度提高了單砂體描述精度。

1 深度域隨機(jī)模擬反演技術(shù)

時間域隨機(jī)模擬反演基于地質(zhì)統(tǒng)計學(xué)理論的序貫高斯配置協(xié)模擬算法。這種算法首先以高精度層序地層研究與地震資料解釋成果建立精細(xì)初始模型，在進(jìn)行了確定性反演后，以波阻抗數(shù)據(jù)體為背景趨勢，依據(jù)序貫高斯配置協(xié)模擬方法沿著任一隨機(jī)路徑進(jìn)行。隨機(jī)路徑的計算過程是依賴于空間變差函數(shù)的統(tǒng)計學(xué)分布的，不同的隨機(jī)路徑得到不同的結(jié)果和實(shí)現(xiàn)，不同實(shí)現(xiàn)的差異反映了地下地質(zhì)的非均質(zhì)性和隨機(jī)性[8-9]。其關(guān)鍵環(huán)節(jié)可以概述為：應(yīng)用測井資料與地震資料建立高精度等時地層格架模型；在地層格架模型的基礎(chǔ)上開展確定性稀疏脈沖反演；擬合各等時地層空間變差的高斯分布；依據(jù)高斯配置協(xié)模擬算法開展變差函數(shù)約束下的隨機(jī)計算[10]。

從基礎(chǔ)原理及關(guān)鍵環(huán)節(jié)中可以看出，時間域隨機(jī)模擬反演方法能夠包含足夠的地質(zhì)信息，更好地反映儲層空間結(jié)構(gòu)變化特性，同時也綜合了測井?dāng)?shù)據(jù)的垂向分辨率和地震數(shù)據(jù)的橫向分辨率優(yōu)勢。除大牛地氣田外，該技術(shù)在國內(nèi)其他油氣田開發(fā)過程中已發(fā)揮出顯著作用。

在實(shí)際應(yīng)用過程中發(fā)現(xiàn)，時間域地質(zhì)統(tǒng)計學(xué)隨機(jī)模擬方法需要建立在較為準(zhǔn)確的時深標(biāo)定關(guān)系上，對于未鉆穿目的層或鉆穿目的層即停鉆的水平井A靶點(diǎn)、導(dǎo)眼井來說，其測井?dāng)?shù)據(jù)無法得到有效的利用，時深標(biāo)定的微小誤差都將導(dǎo)致單砂體頂、底埋深及水平井軌跡設(shè)計的差異。為此，我們針對大牛地氣田水平井軌跡設(shè)計的技術(shù)需求，研究試驗(yàn)高精度深度域隨機(jī)模擬反演技術(shù)。其主要步驟如圖1所示。

圖1 深度域隨機(jī)模擬反演流程

1) 在擬進(jìn)行反演運(yùn)算的三維地震研究區(qū)內(nèi)選一定比例測井質(zhì)量較好的典型井進(jìn)行合成記錄標(biāo)定；典型井是指測井質(zhì)量高、測井?dāng)?shù)據(jù)完整、未發(fā)生井徑擴(kuò)張的優(yōu)質(zhì)井。

2) 僅利用典型井，通過內(nèi)插計算的方式建立初始模型。讀取典型井的聯(lián)井初始模型剖面，并選擇一個或多個具有標(biāo)志性的巖石層。通過對比標(biāo)志巖石層的橫向連續(xù)性微調(diào)典型井的時間-深度標(biāo)定關(guān)系。如此反復(fù)，直到所建立初始模型中的標(biāo)志巖石層橫向連續(xù)，不存在由于標(biāo)定過程中的人為誤差導(dǎo)致的畸變?yōu)橹?。需要注意的是該步驟僅利用測井質(zhì)量較好的典型井，目的是盡可能減小合成記錄標(biāo)定過程中形成的人為誤差，非典型井(如擴(kuò)徑井、水平井導(dǎo)眼段或A靶點(diǎn)測井)不參與初始模型建立。

3) 僅利用典型井作為約束，結(jié)合三維地震數(shù)據(jù)，進(jìn)行時間域稀疏脈沖波阻抗反演[11-12]計算，非典型井不參與運(yùn)算。

4) 選擇合理的速度場，如需要生成井控速度場[13]，則僅利用典型井，非典型井不參與運(yùn)算。

5) 利用步驟4)建立的井控速度場對步驟3)得到的時間域稀疏脈沖波阻抗數(shù)據(jù)體進(jìn)行時深轉(zhuǎn)換，得到深度域稀疏脈沖波阻抗數(shù)據(jù)體，由于稀疏脈沖反演和井控速度場僅利用典型井生成，典型井標(biāo)定中的誤差在本步驟中是相同的，非典型井未參與運(yùn)算，因此誤差得到有效限定。

6) 深度域稀疏脈沖波阻抗數(shù)據(jù)由典型井及三維地震數(shù)據(jù)通過步驟3)至步驟5)共同作用形成，以此作為約束，將所有井全部參與隨機(jī)模擬反演[14-15]，可以得到深度域巖性預(yù)測數(shù)據(jù)體及儲層的埋藏深度信息。由于用于約束的稀疏脈沖波阻抗數(shù)據(jù)是深度域的，因此對于難以得到理想標(biāo)定結(jié)果的井(如擴(kuò)徑井、水平井導(dǎo)眼段或A靶點(diǎn)測井)無需進(jìn)行時深標(biāo)定即可直接參與反演，深度信息不存在誤差。即對于擴(kuò)徑或者無法制作合成記錄的非典型井無需進(jìn)行標(biāo)定，而是直接利用其絕對的深度信息進(jìn)行反演，從而完全避免了非典型井時間-深度標(biāo)定的誤差。

2 深度域隨機(jī)模擬特點(diǎn)分析

如前所述，不同隨機(jī)模擬方法的區(qū)別之處主要在于作為軟約束而存在的稀疏脈沖波阻抗反演上，即在時間域還是在深度域[16]進(jìn)行相應(yīng)的約束。由于稀疏脈沖反演對反射系數(shù)序列進(jìn)行了稀疏處理，也就是說，對于具體的測井資料約束方面只是做到了較強(qiáng)反射系數(shù)界面的抽稀，在細(xì)節(jié)方面并未得到精細(xì)的約束，這就為測井資料的合成記錄標(biāo)定造成了相對的軟約束處理，即僅僅約束較強(qiáng)反射稀疏序列。而隨機(jī)模擬過程則是對測井資料進(jìn)行了嚴(yán)格的約束，測井資料合成記錄標(biāo)定中所存在的誤差將在稀疏脈沖反演軟約束的基礎(chǔ)上進(jìn)一步累積，在時深轉(zhuǎn)換的過程中，只能對已有的約束成果進(jìn)行相應(yīng)的地質(zhì)框架對比及轉(zhuǎn)換，并不能對此產(chǎn)生任何影響和作用。

相比之下，深度域的隨機(jī)模擬反演方案在合成記錄問題上有了很大的改進(jìn)，在將稀疏脈沖反演結(jié)果轉(zhuǎn)換到深度域后，由于計算過程只是對稀疏處理了的反射系數(shù)進(jìn)行反演處理，所得到的結(jié)果只是體現(xiàn)了主導(dǎo)反射系數(shù)的作用，對主導(dǎo)反射系數(shù)進(jìn)行了嚴(yán)格的約束限制[6]。除此之外，測井資料在深度域是絕對的，在與地震資料的結(jié)合與匹配過程中，可以在深度約束上起到主導(dǎo)作用，而在時間域反演中則由地震資料主導(dǎo)時間約束。因此很大程度上可以改進(jìn)由于合成記錄標(biāo)定所導(dǎo)致的細(xì)微誤差，不僅可以進(jìn)一步提高反演的外推能力，外推的精確度也可以得到相應(yīng)的改進(jìn)。

深度域隨機(jī)模擬主要具備幾個方面特點(diǎn)。

1) 強(qiáng)調(diào)以測井質(zhì)量高、數(shù)據(jù)完整、未發(fā)生井徑擴(kuò)張的優(yōu)質(zhì)測井資料作為反演主體；盡管井震合成記錄標(biāo)定具有很多的不確定性，達(dá)不到絕對準(zhǔn)確，但是經(jīng)過優(yōu)選的典型井井間合成記錄標(biāo)定誤差相對較小。

2) 利用典型井測井曲線內(nèi)插形成的初始模型具有很高的分辨率，但是初始模型中地層對比關(guān)系完全由井間合成記錄標(biāo)定結(jié)果決定，通過對橫向穩(wěn)定性較好的標(biāo)志層進(jìn)行地層對比，進(jìn)而微調(diào)合成記錄標(biāo)定時深關(guān)系，可以將典型井井間合成記錄誤差轉(zhuǎn)化為系統(tǒng)差。

3) 經(jīng)過測井曲線與初始模型精細(xì)標(biāo)定后的典型井具有測井質(zhì)量高、橫向地層對比關(guān)系穩(wěn)定的特點(diǎn)，以此為基礎(chǔ)開展稀疏脈沖反演及井控速度場的建立，可以盡可能將標(biāo)定誤差控制在系統(tǒng)差內(nèi)，不至于進(jìn)一步放大。

4) 轉(zhuǎn)換至深度域的稀疏脈沖反演結(jié)果具有忠實(shí)于原始地震數(shù)據(jù)及構(gòu)造特征的優(yōu)點(diǎn)，同時其相對較低的分辨能力也是優(yōu)選為趨勢約束背景的主要原因。由于合成記錄標(biāo)定無法達(dá)到絕對準(zhǔn)確，典型井時間與深度對比關(guān)系仍然具有一定的不確定性。當(dāng)利用稀疏脈沖波阻抗數(shù)據(jù)體進(jìn)行約束時，由于同樣在深度域進(jìn)行計算，測井資料的絕對深度信息將發(fā)揮二次約束作用，從而達(dá)到進(jìn)一步消除標(biāo)定誤差的效果。在此過程中如果選擇高分辨率波阻抗體作為背景約束，地層對比關(guān)系將發(fā)生明顯扭動。

5) 對于測井質(zhì)量較差的直井或測井?dāng)?shù)據(jù)不完整的A靶點(diǎn)測井、導(dǎo)眼井而言，由于直接參與深度域反演計算，從而完全避免了合成記錄標(biāo)定過程，也就是說完全消除了合成記錄標(biāo)定誤差，井間對比關(guān)系及反演外推能力得到加強(qiáng)。

對于水平井而言，目標(biāo)砂體埋深預(yù)測精度及軌跡設(shè)計是部署及實(shí)施的關(guān)鍵因素，而從以上說明可以看出，深度域隨機(jī)模擬的優(yōu)點(diǎn)在于減小低質(zhì)量測井資料時深標(biāo)定的誤差，過多地依靠優(yōu)質(zhì)典型井開展前期數(shù)據(jù)準(zhǔn)備工作，進(jìn)而達(dá)到提高預(yù)測精度的目的(對于直井部署而言意義相對較小)。以此同樣可以推斷，在合成記錄標(biāo)定能夠得到絕對精確的假設(shè)條件下，深度域隨機(jī)模擬與時間域隨機(jī)模擬結(jié)果是一致的。

3 大牛地氣田應(yīng)用效果分析

為進(jìn)一步說明深度域隨機(jī)模擬反演的技術(shù)優(yōu)勢，針對大牛地氣田某井區(qū)開展了時間域隨機(jī)模擬然后作時深轉(zhuǎn)換與深度域隨機(jī)模擬的對比試驗(yàn)。這里需要說明的是，在二者對比中采用的模擬參數(shù)及顯示參數(shù)完全相同，唯一的區(qū)別僅在于時間域約束還是深度域約束。

圖2給出了時間域隨機(jī)模擬轉(zhuǎn)深度域剖面與深度域隨機(jī)模擬反演剖面的對比結(jié)果。從圖2可以看出，不同反演方法得到的剖面上目的層段砂體之間的空間展布特征完全不同，深度域隨機(jī)模擬反演結(jié)果展現(xiàn)出了砂體的橫向連續(xù)性，砂體之間的薄夾層同樣得到了很好的外推和擴(kuò)展；而時間域隨機(jī)模擬進(jìn)而時深轉(zhuǎn)換的剖面盡管同樣可以得到砂體的走向形態(tài)，但是砂體更多是以獨(dú)立的單砂體形式存在，外推能力受限較大。除了砂巖預(yù)測效果外，砂體之間的泥巖隔夾層同樣可以得到較好的刻畫，在深度域反演剖面上泥巖夾層連續(xù)性相對較好，而在時間域隨機(jī)模擬轉(zhuǎn)深度域剖面上則僅體現(xiàn)為井口處的局部現(xiàn)象。

圖2 大牛地氣田某井區(qū)時間域隨機(jī)模擬轉(zhuǎn)深度域剖面(a)與深度域隨機(jī)模擬剖面(b)

由此可見，深度域隨機(jī)模擬反演預(yù)測的砂體連通性增強(qiáng)，單井外推能力得到有效提高，對于薄夾層也給出了精細(xì)的刻畫。在2012年大牛地氣田水平井建產(chǎn)過程中，以泥巖隔夾層精細(xì)描述為基礎(chǔ)的分段壓裂技術(shù)得到了很好的應(yīng)用，形成了“一井一案”的壓裂措施，達(dá)到了增產(chǎn)降本的目的。

圖3為該井區(qū)某層段兩種反演方法得到的砂體預(yù)測厚度圖(反演參數(shù)及顯示參數(shù)相同)，圖中氣井符號標(biāo)注的均為水平井導(dǎo)眼井或A靶點(diǎn)井，細(xì)黑點(diǎn)為直井，所有井均全部參與了反演。從圖3中可以看出，由于大牛地氣田井控程度相對較高，時間域隨機(jī)模擬轉(zhuǎn)深度域預(yù)測結(jié)果(圖3a)與深度域隨機(jī)模擬預(yù)測結(jié)果(圖3b)總體是一致的，但是對于井控程度相對較低的區(qū)域，深度域隨機(jī)模擬的單井外推能力則得到了明顯的體現(xiàn)。

圖3 大牛地氣田某井區(qū)時間域隨機(jī)模擬轉(zhuǎn)深度域(a)與深度域隨機(jī)模擬(b)反演預(yù)測的砂體厚度

單砂體精細(xì)預(yù)測的深度域隨機(jī)模擬反演技術(shù)在大牛地氣田水平井規(guī)模化部署中發(fā)揮了至關(guān)重要的作用，根據(jù)深度域反演結(jié)果進(jìn)行水平軌跡設(shè)計，成功部署水平井200余口，平均砂巖鉆遇率超過94%，氣層鉆遇率較往年提高兩倍以上。結(jié)合水平井分段壓裂技術(shù)有效動用三類地質(zhì)儲量，2012年及2013年間共新建天然氣產(chǎn)能1.9×109m3(約190×104t油當(dāng)量)，成為國內(nèi)水平井整體規(guī)?；渴鸬某晒ο壤?/p>

4 結(jié)束語

1) 深度域隨機(jī)模擬反演計算可以有效削弱井間合成記錄誤差，效果明顯；從已有的地質(zhì)認(rèn)識及實(shí)鉆井跟蹤結(jié)果來看，深度域反演剖面更為合理，外推能力顯著加強(qiáng)。

2) 利用測井資料絕對深度信息開展二次約束，可以有效提高對目標(biāo)砂體埋深的預(yù)測精度，針對水平井實(shí)施而言具有非常重要的意義。

3) 面對勘探開發(fā)對儲層預(yù)測精度要求的不斷提高，結(jié)合工區(qū)沉積特征和油氣藏類型有針對性地配置反演方法是有效提高儲層預(yù)測能力的有效途徑。

參 考 文 獻(xiàn)

[1] 郝蜀民,惠寬洋,李良.鄂爾多斯盆地大牛地大型低滲氣田成藏特征及其勘探開發(fā)技術(shù)[J].石油與天然氣地質(zhì),2006,27(6):762-768

Hao S M,Xi K Y,Li L.Reservoiring features of Daniudi low-permeability gas field in Ordos basin and its exploration and development technologies[J].Oil & Gas Geology,2006,27(6):762-768

[2] 李良,袁志祥,惠寬洋,等.鄂爾多斯盆地北部上古生界天然氣聚集規(guī)律[J].石油與天然氣地質(zhì),2000,21(3):268-282

Li L,Yuan Z X,Xi K Y,et al.Accumulation regularity of upper Paleozoic gas in north Ordos Basin[J].Oil & Gas Geology,2000,21(3):268-282

[3] 盧濤,張吉,李躍剛.蘇里格氣田致密砂巖氣藏水平井開發(fā)技術(shù)及展望[J].天然氣工業(yè),2013,33(8):38-43

Lu T,Zhang J,Li Y G.Hroizontal well development technology for tight sandstone gas reservoirs in the Sulige Gas Field,Ordos Basin[J].Natural Gas Industry,2013,33(8):38-43

[4] 朱新佳,董建輝,任茵.蘇里格氣田蘇53-4井區(qū)水平井整體部署研究[J].天然氣勘探與開發(fā),2012,35(2):44-49

Zhu X J,Dong J H,Ren Y.Integral deployment of horizontal well in Su53-4 well area of Sulige Gas Field[J].Natural Gas Exploration and Development,2012,35(2):44-49

[5] 侯瑞云,劉忠群。鄂爾多斯盆地大牛地氣田致密低滲儲層評價與開發(fā)對策[J].石油與天然氣地質(zhì),2012,33(1):118-128

Hou R Y,Liu Z Q.Reservoir evaluation and development strategies of Daniudi tight sand gas field in the Ordos Basin[J].Oil & Gas Geology,2012,33(1):118-128

[6] 苑書金.鄂爾多斯盆地大牛地氣田下石盒子組地震儲層預(yù)測技術(shù)的研究和應(yīng)用[J].石油地球物理勘探,2008,43(1):48-52

Yuan S J.Study on technique of Lower Shihezi Formation reservoir prediction in Daniudi gasfield of Ordos Basin[J].Oil Geophysical Prospecting,2008,43(1):48-52

[7] 張超英,周小鷹,董寧.測井約束的地震反演在鄂爾多斯盆地大牛地氣田中的應(yīng)用[J].地球物理學(xué)進(jìn)展,2004,19(4):909-917

Zhang C Y,Zhou X Y,Dong N.Seismic inversion constrained by wells in the Daniudi gas-field in Erduosi Basin[J].Progress in Geophysics,2004,19(4):909-917

[8] 孫思敏,彭仕宓.地質(zhì)統(tǒng)計學(xué)反演方法及其在薄層砂體預(yù)測中的應(yīng)用[J].西安石油大學(xué)學(xué)報(自然科學(xué)版),2007,22(1):41-48

Sun S M,Peng S B.Geostatistical inversion method and its application in the prediction of thin reservoirs[J].Journal of Xi’an Shiyou University(Natural Science Edition),2007,22(1):41-48

[9] 王香文,劉紅,滕彬彬.地質(zhì)統(tǒng)計學(xué)反演技術(shù)在薄儲層中的應(yīng)用[J].石油與天然氣地質(zhì),2012,33(5):730-735

Wang X W,Liu H,Teng B B.Application of geostatistical inversion to thin reservoir prediction[J].Oil&Gas Geology,2012,33(5):730-735

[10] 徐立恒,陳顯森,姜巖,等.不同變差函數(shù)變程下隨機(jī)反演與隨機(jī)模擬對比分析[J].物探與化探,2012,36(2):224-227

Xu L H,Chen X S,Jiang Y,et al.A comparative analysis of stochastic inversion and stochastic simula-tion on different ranges of variograms[J].Geophysical & Geochemical Exploration,2012,36(3):224-227

[11] 郭朝斌,楊小波,陳紅岳，等.約束稀疏脈沖反演在儲層預(yù)測中的應(yīng)用[J].石油物探,2006,45(4):397-400

Guo C B,Yang X B,Chen H Y，et al.Constrained sparse pulse inversion research in north of Haitongji Depression[J].Geophysical Prospecting for Petroleum,2006,45(4):397-400

[12] Wang Q,Guo K.Constrained sparse spike inversion applied in the reservoir prediction[J].Well Logging Technology,2008,32(1):35-37

[13] 凌云,郭建明,郭向宇,等.油藏描述中的井震時深轉(zhuǎn)換技術(shù)研究[J].石油物探,2011,50(1):1-13

Ling Yun,Guo Jianming,Guo Xiangyu,et al.Research on time-depth conversion by well to seismic in reservoir characterization[J].Geophysical Prospecting for Petroleum,2011,50(1):1-13

[14] 劉俊州,周小鷹,王箭波.地質(zhì)統(tǒng)計學(xué)反演在非均質(zhì)薄砂巖氣藏開發(fā)中的實(shí)踐[J].內(nèi)蒙古石油化工,2009(7):16-19

Liu J Z,Zhou X Y,Wang J B.Parctices of opening heterogeneous thin sand rock gases deposit in geometry statistic inversion[J].Inner Mongolia Petrochemical Industry,2009(7):16-19

[15] 常炳章,王雪敏,李永清.隨機(jī)模擬反演技術(shù)在泌陽凹陷中南部儲層預(yù)測中的應(yīng)用[J].石油天然氣學(xué)報,2009,31(1):216-219

Chang B Z,Wang X M,Li Y Q.Application of stochastic simulation and inversion in reservoir prediction in the south of Biyang depression[J].Journal of Oil and Gas Technology,2009,31(1):216-219

[16] 張靜,楊勤林,王天琦.深度域地震反演方法探索[J].石油地球物理勘探,2010,45(S1):114-116

Zhang J,Yang Q L,Wang T Q.Probing for seismic inversion in depth domain[J].Oil Geophysical Prospecting,2010,45(S1):114-116