陳更新，趙凡，王建功，鄭紅軍，嚴(yán)耀祖，王愛萍，李積永，胡云鵬

（1.中國地質(zhì)大學(xué)（北京）能源學(xué)院；2.中國石油勘探開發(fā)研究院西北分院；3.中國石油青海油田勘探開發(fā)研究院；4.中國石油勘探開發(fā)研究院）

分區(qū)域多點(diǎn)統(tǒng)計(jì)隨機(jī)地質(zhì)建模方法
——以柴達(dá)木盆地辮狀河三角洲沉積儲集層為例

陳更新1,2，趙凡2，王建功2，鄭紅軍2，嚴(yán)耀祖2，王愛萍2，李積永3，胡云鵬4

（1.中國地質(zhì)大學(xué)（北京）能源學(xué)院；2.中國石油勘探開發(fā)研究院西北分院；3.中國石油青海油田勘探開發(fā)研究院；4.中國石油勘探開發(fā)研究院）

通過對密井網(wǎng)區(qū)精細(xì)地質(zhì)研究，結(jié)合沉積物理模擬實(shí)驗(yàn)，對柴達(dá)木盆地昆北油田辮狀河三角洲沉積模式開展研究，獲得沉積微相單元的統(tǒng)計(jì)特征參數(shù)和平面配置關(guān)系，建立辮狀河三角洲不同沉積區(qū)域的訓(xùn)練圖像，最終建立辮狀河三角洲儲集層地質(zhì)模型。采用角度旋轉(zhuǎn)、比例縮放、鄰區(qū)序貫條件模擬方法，分區(qū)分訓(xùn)練圖像建立精細(xì)的辮狀河三角洲沉積微相模型。分區(qū)域多點(diǎn)統(tǒng)計(jì)隨機(jī)建模間接放寬多點(diǎn)統(tǒng)計(jì)隨機(jī)建模平穩(wěn)性要求，并改善了河道連續(xù)性，提高了訓(xùn)練圖像的可移植性。在柴達(dá)木盆地昆北油田油藏描述應(yīng)用中，模型與生產(chǎn)實(shí)際吻合較好，取得較好應(yīng)用效果，證明該方法的可靠性和可行性，探索了采用隨機(jī)模擬方法準(zhǔn)確刻畫辮狀河三角洲沉積儲集層的新途徑。圖8參29

分區(qū)域多點(diǎn)統(tǒng)計(jì)；隨機(jī)地質(zhì)建模；訓(xùn)練圖像；沉積微相；辮狀河三角洲；柴達(dá)木盆地

0 引言

隨著石油勘探與開發(fā)程度不斷提高，研究領(lǐng)域不斷擴(kuò)大，勘探目標(biāo)地質(zhì)體的空間結(jié)構(gòu)和幾何形態(tài)日益復(fù)雜，油藏和儲集層類型更加多樣，因此區(qū)域化地質(zhì)變量不僅與矢量距離有關(guān)，更與空間位置相關(guān)。傳統(tǒng)的兩點(diǎn)地質(zhì)統(tǒng)計(jì)模擬方法（如截?cái)喔咚鼓M、序貫指示模擬）基于變差函數(shù)表達(dá)兩點(diǎn)間相關(guān)性，其優(yōu)點(diǎn)是易于忠實(shí)條件數(shù)據(jù),不足之處是只能在二階平穩(wěn)或本征假設(shè)前提下，把握空間任意兩點(diǎn)之間的相關(guān)性，因而難以表征復(fù)雜空間結(jié)構(gòu)和幾何形態(tài)（如彎曲河道）。因此Guardiano等提出一種新的基于象元的多點(diǎn)地質(zhì)統(tǒng)計(jì)隨機(jī)模擬方法[1-2]。

相對于兩點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)，多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)應(yīng)用“訓(xùn)練圖像”代替變差函數(shù)表達(dá)多點(diǎn)間相關(guān)性，即地質(zhì)變量的空間結(jié)構(gòu)性。訓(xùn)練圖像是能夠表達(dá)實(shí)際儲集層結(jié)構(gòu)、幾何形態(tài)和分布模式的數(shù)字化圖像，它能夠反映先驗(yàn)地質(zhì)概念和沉積模式。多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)方法仍然采用基于象元的序貫算法，因而易于忠實(shí)硬數(shù)據(jù)，與兩點(diǎn)統(tǒng)計(jì)本質(zhì)區(qū)別在于它通過給定數(shù)據(jù)事件對訓(xùn)練圖像進(jìn)行

掃描來獲取條件概率分布函數(shù)。Strebelle等在前人研究的基礎(chǔ)上，提出Snesim算法[3-4]，將多點(diǎn)統(tǒng)計(jì)概率保存于“搜索樹”中，大大提高計(jì)算效率，解決了需要反復(fù)掃描訓(xùn)練圖像的問題，自此多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)真正成為一種實(shí)用的隨機(jī)建模方法。此后，國外研究者又提出許多新的多點(diǎn)算法或改進(jìn)算法，如Simpat、Filtersim、Growth-sim算法等[5-6]。以吳勝和、王家華為代表的學(xué)者將多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)的概念、內(nèi)涵、算法引入國內(nèi)，并取得一些研究進(jìn)展[7-10]。目前，主流儲集層建模軟件Petrel和RMS中均有多點(diǎn)統(tǒng)計(jì)隨機(jī)建模模塊。

多點(diǎn)地質(zhì)統(tǒng)計(jì)建模在再現(xiàn)儲集層構(gòu)型及數(shù)據(jù)條件化方面，相對于傳統(tǒng)的兩點(diǎn)統(tǒng)計(jì)建模有較大提高，但也存在一些不足。該方法要求訓(xùn)練圖像平穩(wěn)，即訓(xùn)練圖像內(nèi)目標(biāo)體的幾何構(gòu)型及目標(biāo)形態(tài)在全區(qū)基本不變，沒有明顯趨勢或明顯局部變異性。而實(shí)際儲集層受多種因素影響，往往具有復(fù)雜的非均質(zhì)性，即非平穩(wěn)性。非平穩(wěn)訓(xùn)練圖像會導(dǎo)致雜亂無章的模擬實(shí)現(xiàn)[11-12]。Zhang提出通過旋轉(zhuǎn)和比例壓縮將非平穩(wěn)訓(xùn)練圖像變?yōu)槠椒€(wěn)訓(xùn)練圖像，并建立多個訓(xùn)練圖像以獲取未取樣點(diǎn)條件概率分布函數(shù)，是一種簡單化的解決途徑，但具有局限性[13]。另外受遍歷性假設(shè)和隨機(jī)性影響，Snesim算法再現(xiàn)目標(biāo)體連續(xù)性的能力有限，尤其是穿越整個工區(qū)的河道[14]。尹艷樹等針對河流相儲集層提出了基于儲集層骨架的多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)方法，比Snesim和Simpat方法在河道模擬方面更加有效、穩(wěn)定[15]；另外如何建立適宜的訓(xùn)練圖像，還沒有統(tǒng)一的成熟方法。

在辮狀河三角洲的不同沉積區(qū)域，微相類型多，平面配置關(guān)系復(fù)雜，微相單元形態(tài)迥異，利用多點(diǎn)統(tǒng)計(jì)隨機(jī)模擬建立此類儲集層沉積微相模型是研究難點(diǎn)[16]。本次研究針對辮狀河三角洲沉積特點(diǎn)，以及多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)存在的問題，采用分區(qū)域建立訓(xùn)練圖像的方法，對辮狀河三角洲儲集層進(jìn)行多點(diǎn)統(tǒng)計(jì)地質(zhì)建模，精細(xì)表征研究區(qū)沉積微相分布。

圖1 研究區(qū)位置圖

1 多點(diǎn)統(tǒng)計(jì)地質(zhì)建模訓(xùn)練圖像的建立

訓(xùn)練圖像是多點(diǎn)統(tǒng)計(jì)隨機(jī)模擬的基本工具和輸入條件。建立準(zhǔn)確的、包含更多普適地質(zhì)信息的訓(xùn)練圖像是多點(diǎn)建模成功的關(guān)鍵。一般根據(jù)先驗(yàn)認(rèn)識、專業(yè)知識或類比地質(zhì)知識庫，通過露頭建?；颥F(xiàn)代沉積考察、密井網(wǎng)區(qū)精細(xì)地質(zhì)研究、巖心觀察分析、地震信息挖掘等，建立儲集層微相單元結(jié)構(gòu)或空間組合模式[17-21]。研究區(qū)昆北油田位于柴達(dá)木盆地西緣（柴西緣）南段（見圖1），被阿爾金山和昆侖山所圍，兩大山系為柴西緣古近系提供大量的陸源碎屑物質(zhì)，是沉積物的主

要來源。研究區(qū)主要儲集層為古近系路樂河組，沉積物主要來自昆侖山祁漫塔格物源，發(fā)育辮狀河三角洲—湖泊沉積體系，沉積相帶寬，巖性以中—粗砂巖為主，沉積層理和構(gòu)造發(fā)育[22-23]。

本次研究首先根據(jù)沉積物理模擬實(shí)驗(yàn)、現(xiàn)代沉積研究成果以及辮狀河三角洲認(rèn)識，結(jié)合巖心測井資料確定研究區(qū)粗略沉積模式，主要包括沉積相類型、形態(tài)特征，平面、垂向配置關(guān)系等。采用多井對比、井組解剖、生產(chǎn)動態(tài)驗(yàn)證等方法對該模式進(jìn)行擬合，確定研究區(qū)各沉積微相單元統(tǒng)計(jì)特征參數(shù)（寬度、厚度、寬厚比等）。

對研究區(qū)進(jìn)行沉積物理模擬實(shí)驗(yàn)，有助于分析路樂河組（E1+2）沉積期辮狀河三角洲沉積特征（見圖2）。模擬實(shí)驗(yàn)考慮了古地形、古鹽度、古水深、湖盆邊界、構(gòu)造運(yùn)動強(qiáng)度、古氣候、物源供應(yīng)、水流速度和攜砂量等因素對沉積過程的控制作用，其中考慮古鹽度影響是本實(shí)驗(yàn)的特殊之處。采用活動底板分多次進(jìn)行沉降過程模擬，實(shí)驗(yàn)初期地形坡度相對較陡，以充填作用為主，沉積速率快，河道遷移擺動，各類砂壩相繼形成。隨著沉積過程不斷進(jìn)行，地形坡度逐漸變緩，限制了辮狀河三角洲前緣的發(fā)育空間。推測柴西緣古湖水深約25～30 m，實(shí)驗(yàn)裝置內(nèi)最大水深可達(dá)60 cm，設(shè)計(jì)水深比尺為1∶50。柴西緣第三紀(jì)湖泊以微咸—咸水湖為主。下干柴溝組下段（E31）平均氯離子濃度為3 056×10?6，下油砂山組（N21）平均氯離子濃度為2 723×10?6，據(jù)水體鹽度計(jì)算得到柴西緣E31、N21時(shí)期湖盆水體鹽度分別為15‰～16‰、21‰～30‰，因此，實(shí)驗(yàn)中第三紀(jì)湖盆水體鹽度設(shè)計(jì)為15‰～30‰。路樂河組沉積過程模擬設(shè)定邊界條件：西部為阿爾金山脈，南部邊界為現(xiàn)今盆緣以南20～30 km處，北部為賽什騰山—綠梁山，按照一定比例（縱向、橫向和深度上比例分別為1∶28 000，1∶20 000，1∶2 000）設(shè)計(jì)實(shí)驗(yàn)?zāi)Ｐ?。陸上沉積泥巖顏色、蒸發(fā)巖和古生物證據(jù)都表明，柴西緣第三紀(jì)處于干旱環(huán)境[24-25]。路樂河組至下干柴溝組下部沉積時(shí)期氣候干熱，實(shí)驗(yàn)過程按此控制洪水期、中水期和枯水期時(shí)間比例為2∶4∶6。沉積物理模擬實(shí)驗(yàn)結(jié)果為沉積微相的形態(tài)特征和平面配置關(guān)系提供了定性認(rèn)識[26-28]。

圖2 辮狀河三角洲沉積物理模擬實(shí)驗(yàn)結(jié)果

物理實(shí)驗(yàn)結(jié)果表明辮狀河三角洲平原河道外部構(gòu)型呈窄條狀，三角洲前緣河道呈寬帶狀，壩體呈狹長型或透鏡狀。通過對不同微相砂體規(guī)模的定量測量和統(tǒng)計(jì)，結(jié)合前人研究成果，建立咸化湖盆水深、砂體厚度、寬度、流速與湖水密度（鹽度）等參數(shù)之間的經(jīng)驗(yàn)公式，為咸化湖盆砂體預(yù)測提供參考。

式中 d——顆粒直徑，m；g——重力加速度，m/s2； L——砂體厚度，m；M——粉砂及黏土含量，%，根據(jù)數(shù)據(jù)計(jì)算定義為7%；V——河流流速，m/s；W——單河道寬度，m；ρs——泥沙密度，kg/m3，泥沙的濕密度一般為（1.6～1.7）×10?3kg/m3；ρ——咸水密度，kg/m3；t——時(shí)間，s；μ0——清水粘結(jié)系數(shù)，N·s/m2；μ——咸水粘結(jié)系數(shù)，N·s/m2。

本次建模對象為昆北油田E1+2油藏。通過對密井網(wǎng)區(qū)的精細(xì)地質(zhì)研究和構(gòu)型分析，主力含油層段上部為三角洲前緣亞相（Ⅱ-5—Ⅱ-8小層），下部為三角洲平原亞相（Ⅱ-9—Ⅱ-11小層），整體為水進(jìn)退積過程，自下而上，砂體從水上分流河道逐漸過渡到水下分流河道，

到濱淺湖的河口壩與席狀砂。Ⅱ-5—Ⅱ-11小層的微相接觸關(guān)系包括垂向重疊、側(cè)向拼接和側(cè)向分割（見圖3）。根據(jù)單井相以及生產(chǎn)動態(tài)資料，判斷砂體連通關(guān)系，勾畫砂體邊界，進(jìn)而獲取沉積微相單元規(guī)模統(tǒng)計(jì)（見圖4）。

圖3 研究區(qū)垂直物源方向沉積微相剖面圖（圖中數(shù)據(jù)為深度）

圖4 辮狀河三角洲沉積主力砂體規(guī)模概率統(tǒng)計(jì)圖

精細(xì)地質(zhì)研究認(rèn)為研究區(qū)E1+2油藏辮狀河三角洲平原水上分流河道砂體寬度為200～400 m，多點(diǎn)統(tǒng)計(jì)隨機(jī)建模結(jié)果見圖5，寬厚比為30～55，平均41.7；三角洲前緣水下分流河道寬度為300～600 m，河道寬深比為50～80，平均59.2；河口壩長度為300～800 m，寬度為200～400 m，長寬比1.3～2.4，平均1.9；辮狀河三角洲為近物源沉積，坡度較陡，河道形態(tài)以順直或低彎度為主（見圖6）。

通過設(shè)定前期研究獲取的地質(zhì)體形態(tài)、規(guī)模等參數(shù)，縱向以相比例曲線為約束，采用目標(biāo)體隨機(jī)模擬、多重實(shí)現(xiàn)方式初步建立多個訓(xùn)練圖像，然后人工參與進(jìn)行優(yōu)選和后期修改，以求最大限度地反映沉積微相單元的形態(tài)特征和空間組合關(guān)系（見圖5、圖6）。

圖5 Ⅱ-10小層辮狀河三角洲平原多點(diǎn)統(tǒng)計(jì)隨機(jī)建模結(jié)果

圖6 Ⅱ-6小層辮狀河三角洲前緣分區(qū)域多點(diǎn)統(tǒng)計(jì)隨機(jī)建模

2 分區(qū)域多點(diǎn)統(tǒng)計(jì)隨機(jī)建模的實(shí)現(xiàn)

辮狀河三角洲平原主要發(fā)育泛濫平原和水上分流河道兩種微相，平面配置關(guān)系簡單，交織疊置的水上分流河道為主要儲集層，泛濫平原則發(fā)育隔夾層，整個工區(qū)可采用一種訓(xùn)練圖像進(jìn)行模擬。將訓(xùn)練圖像旋轉(zhuǎn)25°，與祁漫塔格物源方向一致；也可通過設(shè)定縱向、橫向、深度3個方向因子在建模過程中對訓(xùn)練圖像進(jìn)行比例縮放。通過角度旋轉(zhuǎn)和比例縮放，可以使訓(xùn)練圖像適應(yīng)相似的沉積模式、不同沉積規(guī)模和不同物源方向的儲集層建模，因此訓(xùn)練圖像具有較強(qiáng)的可移植性。

多點(diǎn)地質(zhì)統(tǒng)計(jì)建模易于實(shí)現(xiàn)數(shù)據(jù)條件化（建模結(jié)果與井?dāng)?shù)據(jù)吻合較好），通過學(xué)習(xí)訓(xùn)練圖像包含的相模式完成隨機(jī)建模。該方法避免人工繪制沉積相的主觀差異性，而且速度比手工勾畫要快（見圖5），當(dāng)然也需要訓(xùn)練圖像包含更多普適性地質(zhì)信息。

辮狀河三角洲前緣微相類型多，平面配置關(guān)系復(fù)雜，形態(tài)迥異，一種訓(xùn)練圖像難以表達(dá)如此復(fù)雜的信息，必須分區(qū)域采用不同訓(xùn)練圖像進(jìn)行隨機(jī)建模，并在分區(qū)邊界上進(jìn)行拼接，可取得較好效果。分區(qū)域建模大大提高了多點(diǎn)統(tǒng)計(jì)建模適用性，可用于建立復(fù)雜多尺度的沉積微相模型。

首先將單井相在平面上投影，同時(shí)參考砂體厚度等值線圖（類似于手工勾畫沉積亞相平面圖），勾畫水下分流河道、河口壩及席狀砂微相的分布區(qū)域。3個區(qū)域分別采用不同訓(xùn)練圖像進(jìn)行模擬：三角洲前緣上游，儲集層砂體主要是水下分流河道交織疊置為主，透鏡狀或狹長的河口壩和水下分流河道在河口處交替分布，河口壩和分流河道末端發(fā)育連片分布的席狀砂體。理論上，各區(qū)域可同時(shí)或按任意順序進(jìn)行模擬，但這樣會造成分區(qū)邊界上明顯的相不連續(xù)。采用鄰區(qū)序貫條件模擬：從已經(jīng)模擬完的區(qū)域的鄰區(qū)中（除第1個區(qū)）選擇一個區(qū)域進(jìn)行模擬，并以模擬完的區(qū)域值為條件序貫進(jìn)行模擬，至所有區(qū)域完成模擬。這樣每個模擬區(qū)內(nèi)，模擬條件不僅包含本區(qū)內(nèi)的屬性值，而且包含模擬完成區(qū)域的值（主要是邊界附近的），從而保證區(qū)域邊界的相連續(xù)性[29]。前提是兩個區(qū)域的訓(xùn)練圖像包含相同的相類型。

多點(diǎn)統(tǒng)計(jì)隨機(jī)建模結(jié)果既反映了訓(xùn)練圖像所包含的沉積模式，又與單井?dāng)?shù)據(jù)完全吻合，更加精細(xì)和定量地表征微相單元，刻畫辮狀河三角洲前緣二級分岔河道，再現(xiàn)河口壩受湖水改造的沉積樣式（見圖6）。

對單井相進(jìn)行數(shù)據(jù)分析，繪制垂向各微相體積百分比曲線（見圖7a），反映相比例隨深度變化規(guī)律。后續(xù)模擬中，各微相體積百分比曲線作為輸入條件約束模型垂向變化，從而保證相模型和井點(diǎn)數(shù)據(jù)統(tǒng)計(jì)的一致性。Ⅱ-6小層上部主要發(fā)育分流間灣，其次為分流河道和席狀砂（見圖7b）；下部主要發(fā)育分流間灣（見圖7c），中部大量發(fā)育分流河道、席狀砂和河口壩（見圖7d）。多點(diǎn)建模結(jié)果不僅平面上符合沉積樣式，也較好地表征了垂向非均質(zhì)性。

需要說明的是，訓(xùn)練圖像和井點(diǎn)統(tǒng)計(jì)的相比例往往不一致，這是建立的模式和實(shí)際儲集層間存在誤差的表現(xiàn)。模擬時(shí)，允許誤差存在，但誤差不能太大，最好小于10%。因此建立訓(xùn)練圖像需考慮井點(diǎn)數(shù)據(jù)的分析結(jié)果。目前還不能通過三維訓(xùn)練圖像體現(xiàn)模型的

縱向非均質(zhì)性，因?yàn)橄囝愋碗S深度變化幾乎不滿足平穩(wěn)性要求，因此借鑒兩點(diǎn)建模采用垂向比例曲線作為約束條件體現(xiàn)模型的縱向非均質(zhì)性：當(dāng)全局相比例和垂向比例同時(shí)設(shè)定時(shí)，全局相比例會被忽略。

圖7 Ⅱ-6小層分區(qū)域多點(diǎn)統(tǒng)計(jì)隨機(jī)建模實(shí)現(xiàn)的縱向沉積微相變化

根據(jù)地質(zhì)研究結(jié)果，Ⅱ-7小層和Ⅱ-8小層的沉積環(huán)境、樣式、規(guī)模非常相似，直接把Ⅱ-7小層的河道訓(xùn)練圖像用來模擬Ⅱ-8小層水下分流河道，也達(dá)到理想建模效果。因此充分證明其可移植性，這也是訓(xùn)練圖像優(yōu)于手工繪制相圖之處。

與辮狀河三角洲相比，以曲流河為供源的正常三角洲，在三角洲平原部分，河道數(shù)量少，形態(tài)相對穩(wěn)定，根據(jù)厚度較容易劃分出河道發(fā)育區(qū)域和泛濫平原區(qū)域。河道發(fā)育區(qū)域基本滿足多點(diǎn)統(tǒng)計(jì)建模的穩(wěn)定性要求；三角洲前緣部分，由于水體與河流動力強(qiáng)度變化，砂體形態(tài)為鳥足狀或朵狀，因此需根據(jù)實(shí)際情況劃分模擬區(qū)域，從而滿足穩(wěn)定性要求。

圖8 切6-325井組示蹤劑監(jiān)測結(jié)果

3 效果分析

分區(qū)域多點(diǎn)統(tǒng)計(jì)隨機(jī)建模在柴達(dá)木盆地昆北油田取得較好應(yīng)用效果。從切6-325井組的示蹤劑監(jiān)測結(jié)果（見圖8）可看出，砂體連通區(qū)域的監(jiān)測井可見示蹤劑，被分流間灣隔擋區(qū)域的切6-327等井則未見示蹤劑顯

示。水驅(qū)數(shù)據(jù)顯示水驅(qū)前緣推進(jìn)具有明顯的方向性：順物源方向水驅(qū)速度快，切6-345井和切6-315井分別為4.76 m/d、3.35 m/d；垂直物源方向明顯較慢，切6-323井和切6-343井分別為1.13 m/d、1.80 m/d。示蹤劑監(jiān)測結(jié)果和水驅(qū)速度證明建模結(jié)果的可靠性，也證明分區(qū)多點(diǎn)統(tǒng)計(jì)建模方法的可行性。

砂體在不同沉積區(qū)域展布特征不同：在河道區(qū)，辮狀分流河道砂體順物源方向交織疊置；席狀砂區(qū)，席狀砂體垂直物源方向連片疊置，呈迷宮狀；河口壩區(qū)，河口壩呈橢圓狀向湖方向延伸，長軸與物源方向一致；分流間灣主要發(fā)育隔夾層，且具有明顯方向性。分流間灣在河道區(qū)域順物源方向呈條帶狀分布，在席狀砂區(qū)與湖岸線方向平行（見圖7）。

多點(diǎn)地質(zhì)統(tǒng)計(jì)建模充分展現(xiàn)了辮狀河三角洲沉積模式，同時(shí)滿足與單井?dāng)?shù)據(jù)吻合較好。模型受井點(diǎn)統(tǒng)計(jì)相概率和訓(xùn)練圖像約束。精細(xì)刻畫發(fā)育儲集層砂體和隔夾層的沉積微相，對注水高效開發(fā)具有較強(qiáng)的指導(dǎo)性，可為注采井網(wǎng)調(diào)整提供依據(jù)，指導(dǎo)井網(wǎng)優(yōu)化部署。

4 結(jié)論

針對辮狀河三角洲多點(diǎn)地質(zhì)統(tǒng)計(jì)建模中建立訓(xùn)練圖像和滿足其平穩(wěn)性的問題，采用分區(qū)域多點(diǎn)隨機(jī)統(tǒng)計(jì)對咸化湖盆辮狀河三角洲進(jìn)行地質(zhì)建模。通過密井網(wǎng)區(qū)精細(xì)地質(zhì)研究和構(gòu)型分析，結(jié)合室內(nèi)辮狀河三角洲沉積物理模擬實(shí)驗(yàn)，獲得研究區(qū)沉積微相單元平面配置關(guān)系、形態(tài)特征及規(guī)模概率分布，采用目標(biāo)體隨機(jī)模擬方法，整合多種信息建立辮狀河三角洲不同沉積區(qū)域的訓(xùn)練圖像。采用鄰區(qū)序貫條件模擬方法實(shí)現(xiàn)分區(qū)分訓(xùn)練圖像多點(diǎn)建模，保證分區(qū)邊界上相的連續(xù)性；通過角度旋轉(zhuǎn)及比例縮放，使得訓(xùn)練圖像適應(yīng)不同物源方向、不同沉積規(guī)模的儲集層建模，具有較好的可移植性；間接地放寬了多點(diǎn)統(tǒng)計(jì)隨機(jī)建模對平穩(wěn)性的要求，更好地體現(xiàn)了河道的連續(xù)性。

分區(qū)域多點(diǎn)統(tǒng)計(jì)隨機(jī)建模既滿足了數(shù)據(jù)條件化，又實(shí)現(xiàn)了訓(xùn)練圖像包含的沉積模式，避免人工繪制相圖的主觀差異性，且模擬速度快，對沉積微相的刻畫更加精細(xì)和定量化。沉積微相模型與生產(chǎn)實(shí)際吻合性較好，證明該方法的可靠性、可行性和實(shí)用性。該方法不僅適用于辮狀河三角洲沉積儲集層建模，對其他沉積類型的微相模擬也具有借鑒意義。

致謝：感謝中國石油勘探開發(fā)研究院西北分院油氣地質(zhì)研究所杜斌山、劉應(yīng)如、張平、張婷靜提供的大量幫助，感謝中國石油勘探開發(fā)研究院西北分院袁劍英教授和中國石油青海油田研究院屈信忠教授的指導(dǎo)。

[1] Cuardiano F,Srivastava R M.Multivariate geostatistics:Beyond bivariate moments[M]//Snares A.Geostatistics Troia’92.Dordrecht:Kluwer Academic Publishers,1993:133-144.

[2] Caers J.Petroleum Geostatistics[M].Texas:Society of Petroleum Engineers,2005:1-12.

[3] Strebelle S,Journel A.Reservoir modeling using multiple-point statistics[R].SPE 71324,2001.

[4] Strebelle S.Conditional simulation of complex geological structures using multiple-point statistics[J].Mathematical Geology,2002,34(1):1-21.

[5] Arpat G B,Cares J.A multiple-scale,pattern-based approach to sequential simulation[J].Quantitative Geology &Geostatistics,2005,12:255-264.

[6] Arpat G B.Sequential simulation with patterns[D].California:Stanford University,2005:1-85.

[7] 王家華.迎接油氣儲層建模理論、應(yīng)用的大發(fā)展:從2007年國際石油地質(zhì)統(tǒng)計(jì)學(xué)大會談起[J].地學(xué)前緣,2008,15(1):16-25.Wang Jiahua.Greeting the great development of theories and applications of reservoir modeling[J].Earth Science Frontiers,2008,15(1):16-25.

[8] 吳勝和,李宇鵬.儲層地質(zhì)建模的現(xiàn)狀與展望[J].海相油氣地質(zhì),2007,12(3):53- 60.Wu Shenghe,Li Yupeng.Reservoir modeling:Current situation and development prospect[J].Marine Origin Petroleum Geology,2007,12(3):53- 60.

[9] 吳勝和,李文克.多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué):理論、應(yīng)用與展望[J].古地理學(xué)報(bào),2005,7(1):137-144.Wu Shenghe,Li Wenke.Multiple-point geo-statistics:Theory,application and perspective[J].Journal of Palaeogeography,2005,7(1):137-144.

[10] 尹艷樹,吳勝和.儲層建模方法研究進(jìn)展[J].天然氣地球科學(xué),2006,17(2):200-206.Yin Yanshu,Wu Shenghe.The progress of reservoir stochastic modeling[J].Natural Gas Geoscience,2006,17(2):200-206.

[11] Manuel L,Carrera J.Application of multiple point geostatistics to non-stationary images[J].Mathematical Geosciences,2009,36(8):29-42.

[12] 李宇鵬,吳勝和.儲集層構(gòu)型分級套合模擬方法[J].石油勘探與開發(fā),2013,40(5):630-635.Li Yupeng,Wu Shenghe.Hierarchical nested simulation approach in reservoir architecture modeling[J].Petroleum Exploration and Development,2013,40(5):630-635.

[13] Zhang T.Rotation and affinity invariance in multiple-point geostatistics[R].California:Stanford Center of Reservoir Forecast Annual Meeting,2002.

[14] Arpat G B,Caers J.Conditional simulation with patterns[J].Mathematical Geology,2007,39(2):177-203.

[15] 尹艷樹,吳勝和,張昌民,等.基于儲層骨架的多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)方法[J].中國科學(xué):地球科學(xué),2008,38(ZKⅡ):157-164.Yin Yanshu,Wu Shenghe,Zhang Changmin,et al.A reservoir skeleton-based multiple point geostatistics method[J].SCIENCE CHINA Earth Sciences,2009,52(zkⅠ):171-178.

[16] 張團(tuán)峰.在儲層建模中利用多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)整合地質(zhì)概念模型及其解釋[J].地學(xué)前緣,2008,15(1):26-34.Zhang Tuanfeng.Incorporating geological conceptual modes and interpretations into reservoir modeling using multiple-point geo-statistics[J].Earth Science Frontiers,2008,15(1):26-34.

[17] Zhang T,Switzer P,Journel A.Filter-based classification of training image patterns for spatial simulation[J].Mathematical Geology,2006,38(1):63-80.

[18] 李少華,張昌民.儲層建模中幾種原型模型的建立[J].沉積與特提斯地質(zhì),2004,24(3):102-106.Li Shaohua,Zhang Changmin.The construction of prototype models in reservoir modeling[J].Sedimentary Geology and Tethyan Geology,2004,24(3):102-106.

[19] 王家華,趙巍.基于地震約束的地質(zhì)統(tǒng)計(jì)學(xué)建模方法研究[J].海洋石油,2010,30(4):46-49.Wang Jiahua,Zhao Wei.Research on geo-statistical modeling method constrained by seismic data[J].Offshore Oil,2010,30(4):46-49.

[20] GonzáLez E F,Mukerji T,Mavko G.Seismic inversion combining rock physics and multiplepoint geostatisitics[J].Geophysics,2008,73(1):11-21.

[21] Strebelle S.Simulation of petro-physical property trends within facies geobodies[J].Petroleum Geostatistics,2007,9:10-14.

[22] 羅曉容,孫盈,汪立群,等.柴達(dá)木盆地北緣西段油氣成藏動力學(xué)研究[J].石油勘探與開發(fā),2013,40(2):159-170.Luo Xiaorong,Sun Ying,Wang Liqun,et al.Dynamics of hydrocarbon accumulation in the west section of the northern margin of the Qaidam Basin,NW China[J].Petroleum Exploration and Development,2013,40(2):159-170.

[23] 孫平,郭澤清,劉衛(wèi)紅,等.柴達(dá)木盆地英東一號油氣田成藏機(jī)理[J].石油勘探與開發(fā),2013,40(4):429-435.Sun Ping,Guo Zeqing,Liu Weihong,et al.Accumulation mechanism of the Yingdong I field in the Qaidam Basin,NW China[J].Petroleum Exploration and Development,2013,40(4):429-435.

[24] Zhao F.Lacustrine algal limestone reservoir in western Qaidam Basin,China[J].Carbonates and Evaporates,2013,28:182-191.

[25] Zhao F.Identification of Neogene mixed lacustrine carbonate in western Qaidam basin[J].Carbonates and Evaporates,2014,29:210-217.

[26] 張春生.碎屑巖沉積模擬技術(shù)[M].北京:石油工業(yè)出版社,2003:3-12.Zhang Chunsheng.Sedimentary physical simulation technology of clastic rock[M].Beijing:Petroleum Industry Press,2003:3-12.

[27] 張春生,劉忠保,施冬,等.利用數(shù)值模擬方法預(yù)測辮狀河砂體幾何形態(tài)[J].石油勘探與開發(fā),2004,31(增刊):85-88.Zhang Chunsheng,Liu Zhongbao,Shi Dong,et al.Braid river sand body shape forecast by numerical simulation method[J].Petroleum Exploration and Development,2004,31(Supp.):85-88.

[28] 王俊輝,姜在興.三角洲沉積的物理模擬[J].石油與天然氣地質(zhì),2013,34(6):758-763.Wang Junhui,Jiang Zaixing.Physical simulation of deltaic deposits[J].Oil &Gas Geology,2013,34(6):758-763.

[29] Remy N,Boucher A,Wu Jianbing.Applied geostatistics with SGeMS:A user’s guide[M].Cambridge:Cambridge University Press,2009:132-214.

（編輯 林敏捷）

Regionalized multiple-point stochastic geological modeling:A case from braided delta sedimentary reservoirs in Qaidam Basin,NW China

Chen Gengxin1,2,Zhao Fan2,Wang Jiangong2,Zheng Hongjun2,Yan Yaozu2,Wang Aiping2,Li Jiyong3,Hu Yunpeng4
(1.School of Energy Resources,China University of Geosciences,Beijing 100083,China;2.PetroChina Research Institute of Petroleum Exploration &Development-Northwest,Lanzhou 730020,China;3.PetroChina Qinghai Oilfield,Dunhuang 736200,China;4.PetroChina Research Institute of Petroleum Exploration &Development,Beijing 100083,China)

Research on the depositional model of a braided delta in Kunbei oilfield,Qaidam Basin,was conducted with a fine geologic study in a dense well network area by combining sedimentary physical simulation to obtain statistical characteristic parameters and allocation relations of plane sedimentary facies units.On this basis,training images in different sedimentary regions of the braided delta were developed for reservoir modeling.In each sub-region,sequential conditional simulation in adjacent regions used different training images with angle revolving and scaling.This simulation succeeded in finely modeling the sedimentary microfacies of the braided delta.This method relaxed the requirement for a stationary training image in multipoint simulation and improved channel continuity as well as portability of the training image.In the fine reservoir description of the Kunbei oilfield,the model results were consistent with practical production and have the potential to guide effective water-flooding development.Effective application demonstrated the feasibility and practicability of the method.This paper presents a new method for accurately characterizing the braided delta using stochastic simulation.

regionalized multi-point statistics;stochastic geological modeling;training image;sedimentary microfacies;braided delta;Qaidam Basin

國家科技重大專項(xiàng)（2011ZX05007-006）；中國石油重大科技專項(xiàng)（2011E-0308）

TE122

A

1000-0747(2015)05-0638-08

10.11698/PED.2015.05.11

陳更新（1984-），男，山東菏澤人，中國地質(zhì)大學(xué)（北京）在讀博士，主要從事開發(fā)地質(zhì)、儲集層建模和油藏?cái)?shù)值模擬研究。地址：甘肅省蘭州市雁兒灣路535號，中國石油勘探開發(fā)研究院西北分院油氣地質(zhì)研究所，郵政編碼：730020。E-mail：chen_gx@petrochina.com.cn

2014-04-11

2015-07-26