尹艷樹，張昌民，石書緣，馮文杰，和景陽(yáng) （長(zhǎng)江大學(xué)地球科學(xué)學(xué)院，湖北荊州434023）

綜合隨機(jī)游走過程與多點(diǎn)統(tǒng)計(jì)的河流相建模新方法

尹艷樹，張昌民，石書緣，馮文杰，和景陽(yáng) （長(zhǎng)江大學(xué)地球科學(xué)學(xué)院，湖北荊州434023）

多點(diǎn)統(tǒng)計(jì)學(xué)由于其算法的限制，不能很好再現(xiàn)連續(xù)河道形態(tài)，其關(guān)鍵是在抽樣過程中隨機(jī)性較強(qiáng)，導(dǎo)致河道連續(xù)性中斷。提出了將隨機(jī)游走過程應(yīng)用于多點(diǎn)統(tǒng)計(jì)建模的新方法。通過隨機(jī)游走產(chǎn)生河道主流線，并利用河道主流線約束多點(diǎn)統(tǒng)計(jì)預(yù)測(cè)，克服多點(diǎn)統(tǒng)計(jì)隨機(jī)抽樣導(dǎo)致的河道不連續(xù)性問題。從實(shí)際建模效果看，基于隨機(jī)游走過程與多點(diǎn)統(tǒng)計(jì)的耦合建模方法較好地再現(xiàn)了河道連續(xù)性；抽稀檢驗(yàn)驗(yàn)證了耦合建模方法的穩(wěn)健性；表明所設(shè)計(jì)的新方法可以應(yīng)用于實(shí)際儲(chǔ)層建模。

隨機(jī)游走過程；多點(diǎn)地質(zhì)統(tǒng)計(jì)；綜合建模；河流相

多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)儲(chǔ)層建模方法已經(jīng)在國(guó)內(nèi)油田應(yīng)用。雖然在儲(chǔ)層形態(tài)再現(xiàn)上較傳統(tǒng)的序貫指示建模方法好，但是與基于目標(biāo)建模方法相比，儲(chǔ)層連續(xù)性仍然存在較大問題［1～5］。最典型的就是河道連續(xù)性中斷，其原因主要有2個(gè)方面：①SNESIM多點(diǎn)統(tǒng)計(jì)方法仍然是單點(diǎn)估計(jì)，不能保證河道連續(xù)性［6］。②SIMPAT多點(diǎn)統(tǒng)計(jì)方法雖然通過數(shù)據(jù)樣板整體替換，一定程度上改善了單點(diǎn)估計(jì)的不足，但其數(shù)據(jù)樣板選擇的隨機(jī)性導(dǎo)致河道連續(xù)性中斷［7］；Filtersim方法是對(duì)SIMPAT計(jì)算效率的完善［8］，對(duì)模擬連續(xù)性改善效果有限。為了解決此問題，將隨機(jī)游走過程耦合到多點(diǎn)地質(zhì)統(tǒng)計(jì)方法中，通過隨機(jī)游走過程產(chǎn)生河道主流線，并作為條件數(shù)據(jù)約束多點(diǎn)統(tǒng)計(jì)預(yù)測(cè)，較好地再現(xiàn)了河道連續(xù)性，提高了儲(chǔ)層建模精度。

1 隨機(jī)游走過程

隨機(jī)游走過程是一種不規(guī)則的變動(dòng)形式。在這種變動(dòng)過程中，每一步的變化都純粹是隨機(jī)性的［9］。其廣泛應(yīng)用于金融、圖像處理與識(shí)別、地表水污染、電子信息等行業(yè)。在油氣建模領(lǐng)域，隨機(jī)游走已經(jīng)應(yīng)用于辮狀河道相的二維隨機(jī)建模［10］。其建模的核心思想是將研究區(qū)域離散成網(wǎng)格系統(tǒng)；然后在邊界區(qū)域內(nèi)尋找河道源頭的位置，并在此基礎(chǔ)上依次獲取河道主流線的位置。一旦主流線形成，就可以沿著主流線加寬處理，形成河道二維平面模型。

在辮狀河道建模中［10］，定義了河道隨機(jī)游走的5個(gè)方向 （圖1（a））。但是從衛(wèi)星圖片觀察現(xiàn)代河流形態(tài)可以發(fā)現(xiàn)，河道流線變化非常復(fù)雜。對(duì)于高彎度河道，河道流向發(fā)生反轉(zhuǎn)情況常有發(fā)生。此時(shí)利用5個(gè)方向隨機(jī)游走模型就不能表征河道流向反轉(zhuǎn)情況。需要增加更多的方向以描述河道流向復(fù)雜變化。在原模型基礎(chǔ)上，增加了2個(gè)反方向游走的特征 （圖1（b）），以描述河道高彎曲甚至反轉(zhuǎn)的情況。新增2個(gè)方向的游走概率為：

圖1 隨機(jī)游走方向示意圖 （實(shí)線箭頭為游走方向，頂部虛線箭頭代表物源方向）

2 多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)

多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)是近10年涌現(xiàn)出的一種新的建模方法，在國(guó)外已經(jīng)得到廣泛研究和應(yīng)用。自2005年引入國(guó)內(nèi)后，迅速引起重視。在方法、原理以及應(yīng)用上都得到了深入研究。其建模過程如下［11］：①掃描訓(xùn)練圖像，建立搜索樹。②條件數(shù)據(jù)分配到最臨近的網(wǎng)格，定義隨機(jī)訪問路徑。③在每一個(gè)未取樣位置，保留包含在最大搜索數(shù)據(jù)樣板內(nèi)的條件數(shù)據(jù)，從搜索樹里提取多點(diǎn)概率分布，計(jì)算局部條件概率，建立未取樣點(diǎn)局部條件概率分布；如果條件數(shù)據(jù)為0，則以全局概率代替局部概率；如果條件數(shù)據(jù)構(gòu)成的數(shù)據(jù)事件在訓(xùn)練圖像中很少出現(xiàn)，則可以通過去掉離中心位置處最遠(yuǎn)的點(diǎn)，從而使得條件數(shù)據(jù)構(gòu)成的數(shù)據(jù)事件在訓(xùn)練圖像中重復(fù)次數(shù)足夠多，進(jìn)而局部概率的估計(jì)較可靠。④蒙特卡羅抽樣決定未取樣位置的值，并將實(shí)現(xiàn)值加入到條件數(shù)據(jù)中。⑤重復(fù)步驟③和步驟④，直到所有未取樣網(wǎng)格都得到訪問；如果需要多個(gè)實(shí)現(xiàn)，只需改變隨機(jī)訪問路徑，重復(fù)步驟③、④、⑤即可。

3 隨機(jī)游走過程與多點(diǎn)統(tǒng)計(jì)耦合

多點(diǎn)統(tǒng)計(jì)學(xué)由于其算法的限制，不能很好再現(xiàn)連續(xù)河道形態(tài)，其關(guān)鍵是在抽樣過程中隨機(jī)性較強(qiáng)，導(dǎo)致河道連續(xù)性中斷。如果將隨機(jī)游走形成的河道主流線作為條件數(shù)據(jù)約束隨機(jī)抽樣，將能保證河道連續(xù)性。同時(shí)，由于河道剖面形態(tài)通過多點(diǎn)統(tǒng)計(jì)來(lái)描述，將能夠更好忠實(shí)于條件數(shù)據(jù)，避免了基于目標(biāo)方法需要對(duì)河道形態(tài)簡(jiǎn)單化和參數(shù)化的困難。此外，在開發(fā)中、后期密井網(wǎng)條件下，其模擬也能很好條件化，避免了反復(fù)迭代收斂的難題。

設(shè)計(jì)了將隨機(jī)游走過程與多點(diǎn)統(tǒng)計(jì)耦合的建模方法。其基本思想是，首先通過條件數(shù)據(jù)進(jìn)行隨機(jī)游走模擬，利用隨機(jī)游走過程確定河道主流線，保留主流線位置；其次，將所有河道主流線穿越網(wǎng)格點(diǎn)條件化，作為河道數(shù)據(jù)，并與初始條件數(shù)據(jù)一起約束多點(diǎn)統(tǒng)計(jì)，從而更好再現(xiàn)河道形態(tài)。以SNESIM方法為例，2種方法耦合步驟如下：①讀取參數(shù)文件，將區(qū)域網(wǎng)格化；②將條件數(shù)據(jù)分配到最近網(wǎng)格點(diǎn)；③定義不同方向遷移系數(shù)，進(jìn)行隨機(jī)游走模擬；④將隨機(jī)游走生成的主流線網(wǎng)格化，并將網(wǎng)格分配為河道代碼；⑤進(jìn)行SNESIM模擬；⑥輸出模擬結(jié)果。如圖2所示。

4 實(shí)例研究

通過Google Earth軟件截取現(xiàn)代河流圖像，并取樣形成條件數(shù)據(jù)點(diǎn)，對(duì)模擬方法進(jìn)行檢驗(yàn)。圖3（a）是用Google Earth截取漢江武漢段高彎曲河道，河道主流線發(fā)生反轉(zhuǎn)。分別應(yīng)用改進(jìn)前后的隨機(jī)游走過程結(jié)合多點(diǎn)統(tǒng)計(jì)對(duì)河道進(jìn)行模擬。從模擬結(jié)果看，改進(jìn)前的隨機(jī)游走過程受方向限制，不能很好反映河道高彎曲反轉(zhuǎn)特征 （圖3（b）），而改進(jìn)后的隨機(jī)游走過程則較好再現(xiàn)河道反轉(zhuǎn)形態(tài)（圖3 （c））。

圖4（a）是長(zhǎng)江南京段局部圖像，反映了多河道特征。圖4（b）是僅用多點(diǎn)統(tǒng)計(jì)SNESIM方法模擬結(jié)果，圖4（c）則是耦合了隨機(jī)游走與多點(diǎn)統(tǒng)計(jì)模擬結(jié)果。從圖4（b）和圖4（c）可以明顯的看出，耦合方法更好反映了河道連續(xù)性特征，表明耦合方法具備優(yōu)勢(shì)。

圖2 建模流程圖

圖3 高彎曲河道模擬結(jié)果比較

圖4 河道模擬效果比較

為了驗(yàn)證算法的穩(wěn)定性，將條件數(shù)據(jù)點(diǎn)減少，檢驗(yàn)新方法的穩(wěn)定性。雖然抽稀符合率都為100%。但是從圖5可以看出，僅多點(diǎn)統(tǒng)計(jì)SNESIM方法模擬河道出現(xiàn)中斷，形態(tài)也有變形；但是耦合后建模方法則仍然較好地反映了河道連續(xù)性形態(tài)，表明耦合方法具備穩(wěn)健性。

圖5 河道條件數(shù)據(jù)點(diǎn)抽稀模擬效果比較

5 結(jié) 論

1）提出了河道隨機(jī)游走的7個(gè)方向，解決了傳統(tǒng)的隨機(jī)游走方法僅設(shè)計(jì)河道游走的5個(gè)方向而沒有考慮高彎曲河道流線反轉(zhuǎn)問題。

2）設(shè)計(jì)了將隨機(jī)游走與多點(diǎn)統(tǒng)計(jì)耦合建模方法。將隨機(jī)游走產(chǎn)生河道主流線作為條件數(shù)據(jù)約束多點(diǎn)統(tǒng)計(jì)預(yù)測(cè)。實(shí)際建模效果表明，耦合建模方法較傳統(tǒng)的多點(diǎn)統(tǒng)計(jì)方法更好再現(xiàn)了河道連續(xù)性。表明利用隨機(jī)游走過程約束確實(shí)提高了多點(diǎn)統(tǒng)計(jì)建模效果。

［1］Liu Y，Journel A G.Improving sequential simulation with a structured path guided by information content［J］.Mathematical Geology，2004，36 （8）：945～964.

［2］Caers J，Avseth P，Mukerji T.Geostatistical integration of rock physics，seismic amplitudes and geological models in North-Sea turbidite systems ［J］.SPE 71321，2001.

［3］Caers J，Zhang T.Multiple-point geostatistics：aquantitative vehicle for integrating geologic analogs into multiple reservoir models［M］.AAPG memoir：Integration of outcrop and modern analog data in reservoir models，2002.

［4］吳勝和，李文克.多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)——理論、應(yīng)用與展望 ［J］.古地理學(xué)報(bào)，2005，7（1）：137～144.

［5］李少華，張昌民，何幼斌，等.河道砂體內(nèi)部物性分布趨勢(shì)的模擬 ［J］.石油天然氣學(xué)報(bào)，2009，31（1）：23～25.

［6］Ezequiel F González.Seismic inversion combining rock physics and multiple-point geostatistics ［J］.Geophysics，2008，73 （1）：11～21.

［7］尹艷樹，吳勝和，張昌民，等.基于儲(chǔ)層骨架的多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)方法 ［J］.中國(guó)科學(xué)，2008，38（增刊Π）：157～164.

［8］Zhang T，Switzer P，Journel A G.Filter-based classification of training image patterns for spatial simulation ［J］.Mathematical Geology，2006，38 （1）：63～80.

［9］雷鈺麗，李陽(yáng)，王崇駿，等.基于權(quán)重的馬爾可夫隨機(jī)游走相似度度量的實(shí)體識(shí)別方法 ［J］.河北師范大學(xué)學(xué)報(bào) （自然科學(xué)版），2010，34 （1）：26～30.

［10］王家華，張團(tuán)峰.油氣儲(chǔ)層隨機(jī)建模 ［M］.北京：石油工業(yè)出版社，2001.90～94.

［11］李少華，尹艷樹，張昌民.儲(chǔ)層隨機(jī)建模系列技術(shù) ［M］.北京：石油工業(yè)出版社，2007.56～57.

Modeling on Integrating Random Walk Process and Multiple-point Geostatistics to Fluvial Reservoirs

YIN Yan-shu，ZHANG Chang-min，SHI Shu-yuan，F(xiàn)ENG Wen-jie，HE Jing－yang（Authors’Address：College of Geosciences，Yangtze University，Jingzhou434023，Hubei，China）

The multiple-point geostatistics could not represent the continuous shape of the channel due to the limitation of its algorithm with random draw of the simulated value.A new method was proposed，which integrated the random walk process with multiple-point geostatistics.First the random walk process produced main streamline of cannel，by which the main streamline was constrained by the multiple-point geostatistics，the method was used to solve the problem of non-continuities of channel.Based on the actual model，the method based on random walk process and the multiple-point geostatistics can well reproduce the continuity of channel.The results show the new method is validate，it proves that the new method is robust and can be used for reservoir modeling.

random walk process；multiple-point geostatistics；integrating modeling；fluvial reservoir

TE122.2

A

1000-9752（2011）08-0044-04

2011-02-21

國(guó)家科技重大專項(xiàng) （2008ZX05011-3）；國(guó)家自然科學(xué)基金項(xiàng)目 （40902043）；湖北省自然科學(xué)基金項(xiàng)目 （2008CDB390）。

尹艷樹 （1978-），男，2000年江漢石油學(xué)院畢業(yè)，博士，副教授，現(xiàn)主要從事儲(chǔ)層建模的方法與教學(xué)研究工作。

［編輯］ 宋換新