李 苗，耿 懿，李 璇

(1.長(zhǎng)慶油田分公司 第五采油廠，陜西 西安 710200； 2.中國(guó)科學(xué)院 蘭州文獻(xiàn)情報(bào)中心，甘肅 蘭州 730030； 3.長(zhǎng)慶油田分公司 第二采油廠，甘肅 慶陽(yáng) 735101)

?

地震約束降低儲(chǔ)層建模中的不確定性

——以委內(nèi)瑞拉MPE3區(qū)塊超重油油藏O-11a為目標(biāo)層研究

李 苗1，耿 懿2，李 璇3

(1.長(zhǎng)慶油田分公司 第五采油廠，陜西 西安 710200； 2.中國(guó)科學(xué)院 蘭州文獻(xiàn)情報(bào)中心，甘肅 蘭州 730030； 3.長(zhǎng)慶油田分公司 第二采油廠，甘肅 慶陽(yáng) 735101)

委內(nèi)瑞拉奧里諾科重油帶，地下砂體展布復(fù)雜，非均質(zhì)性強(qiáng)，傳統(tǒng)兩步建模法難以取得可靠結(jié)果。結(jié)合地震數(shù)據(jù)通過(guò)刻畫(huà)訓(xùn)練圖像及生成砂體概率曲線，對(duì)研究區(qū)的沉積微相模擬結(jié)果進(jìn)行了不確定性分析；利用多點(diǎn)地質(zhì)統(tǒng)計(jì)法分別建立無(wú)地震數(shù)據(jù)約束及有地震數(shù)據(jù)約束的2組沉積微相模型。對(duì)比結(jié)果表明：有地震約束的模型，砂體展布與地震數(shù)據(jù)基本一致，模型不確定性明顯降低。該地震數(shù)據(jù)約束儲(chǔ)層地質(zhì)建模的正確和可靠，對(duì)研究區(qū)后續(xù)水平井的分布及鄰區(qū)的勘探開(kāi)發(fā)有很好的借鑒作用。

地質(zhì)建模；地震數(shù)據(jù)約束；訓(xùn)練圖像；不確定性

李 苗，耿 懿，李 璇.地震約束降低儲(chǔ)層建模中的不確定性:以委內(nèi)瑞拉MPE3區(qū)塊超重油油藏O-11a為目標(biāo)層研究[J].西安石油大學(xué)學(xué)報(bào)(自然科學(xué)版)，2016，31(3):39-47,56.

LI Miao,GENG Yi,LI Xuan.Reducing uncertainty of reservoir modeling by seismic constraints:taking O-11a super-heavy oil reservoir in MPE3 block of Venezuela as an example[J].Journal of Xi'an Shiyou University(Natural Science Edition),2016,31(3):39-47,56.

引　言

近年來(lái)，儲(chǔ)層建模技術(shù)從起初的確定性建模、隨機(jī)建模，到測(cè)井資料與地震資料共同約束的建模[1]，再到地震反演波阻抗等約束的建模[2]，有了較大的發(fā)展。

傳統(tǒng)建模方法即為兩步法：首先由測(cè)井?dāng)?shù)據(jù)建立起研究區(qū)的沉積相模型，再由沉積相模型建立出研究區(qū)儲(chǔ)層的物性模型。這種方法產(chǎn)生的結(jié)果不確定性很大。目前國(guó)內(nèi)對(duì)于不確定性的研究主要是針對(duì)不確定性的評(píng)價(jià)，對(duì)于不確定性建模方法的研究較少[3]。相對(duì)于傳統(tǒng)建模方法不能反映復(fù)雜地質(zhì)體等缺點(diǎn)，Bortoli等[4]在20世紀(jì)90年代提出的基于多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)的建模方法在模擬過(guò)程中的優(yōu)勢(shì)明顯[5-7]。2001年由于搜索樹(shù)的提出，多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)進(jìn)入實(shí)質(zhì)應(yīng)用階段[7-8]。這種方法可以利用訓(xùn)練圖像代替目標(biāo)體，也可以克服必須要對(duì)目標(biāo)體有精準(zhǔn)認(rèn)識(shí)的不足[9]。訓(xùn)練圖像是結(jié)合了多種數(shù)據(jù)的綜合體，作為地質(zhì)體的一種數(shù)字化表示[10]。三維訓(xùn)練圖像可以描述物性在空間上的變化，可以反映出儲(chǔ)層在橫向及縱向上的變化。在已有基礎(chǔ)資料的前提下，可以通過(guò)多種方法刻畫(huà)用于建模的訓(xùn)練圖像。這使得復(fù)雜的地質(zhì)形態(tài)可以得到較為正確的描述。根據(jù)砂體厚度和地震信息獲得訓(xùn)練圖像就是其中的2種方法。

委內(nèi)瑞拉奧里諾科重油帶MPE3區(qū)塊砂體看似均勻，但內(nèi)部具有非常復(fù)雜的非均質(zhì)性，采用傳統(tǒng)兩步建模方法難以取得較為可信的結(jié)果，基于此，本文通過(guò)嘗試針對(duì)MPE3區(qū)塊刻畫(huà)訓(xùn)練圖像，對(duì)比不同訓(xùn)練圖像下有無(wú)地震約束建立模型的精度，確定有地震約束的儲(chǔ)層建模中不確定性大為降低，為研究區(qū)后期整體開(kāi)發(fā)方案的確定提供較為準(zhǔn)確的依據(jù)。

1　訓(xùn)練圖像的刻畫(huà)

1.1研究區(qū)概況

MPE3區(qū)塊位于委內(nèi)瑞拉奧里諾科重油帶東部，南北約16 km，東西約10 km，總計(jì)面積約150 km2。區(qū)內(nèi)共有24口探井。Morichal段是本區(qū)最主要的含油層段，是以辮狀河河道為主要沉積體的沉積[11]。本文主要目的層為位于Morichal的頂層的O-11a。

前期研究表明，24口直井探井中共有19口探井具有完整的巖相數(shù)據(jù)，能夠用來(lái)參與建模。由于這19口直井在工區(qū)內(nèi)分布較為均勻，井距比較接近，因此建模結(jié)果對(duì)該區(qū)具有參考意義。圖1為研究區(qū)三維井位分布。

圖1　研究區(qū)三維井位分布Fig.1　3D well distribution of the study area

1.2地震約束建模原理

Journel等[12]提出了一種多數(shù)據(jù)結(jié)合的建模方法。原理如下：假設(shè)有3種數(shù)據(jù)A、B、C，A表示將要被模擬的數(shù)據(jù)，B表示精確的數(shù)據(jù)(比如測(cè)井?dāng)?shù)據(jù))，C作為A的同位協(xié)同數(shù)據(jù)(本文中為地震數(shù)據(jù))。利用數(shù)據(jù)C約束，通過(guò)數(shù)據(jù)B建立模型。

根據(jù)貝葉斯公式，數(shù)據(jù)A、B、C存在數(shù)理關(guān)系：

(1)

(2)

(3)

其中

利用參數(shù)τ1和τ2，可以定義：

(4)

多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)的建模流程是通過(guò)基礎(chǔ)地質(zhì)研究(包括地質(zhì)認(rèn)識(shí)、測(cè)井資料和地震數(shù)據(jù))綜合建立訓(xùn)練圖像，由地震反演數(shù)據(jù)處理得到砂體概率數(shù)據(jù)，再由測(cè)井?dāng)?shù)據(jù)得到垂向砂體概率。多點(diǎn)統(tǒng)計(jì)建模中的地震數(shù)據(jù)約束可以通過(guò)生成砂體概率曲線的方法實(shí)現(xiàn)[13]。砂體概率曲線可以將波阻抗數(shù)據(jù)轉(zhuǎn)化為3種微相的砂體空間分布概率。在多點(diǎn)統(tǒng)計(jì)建模中，利用軟數(shù)據(jù)約束。砂體空間分布概率是由波阻抗得出，所以建模實(shí)質(zhì)依然是利用波阻抗數(shù)據(jù)進(jìn)行的。

1.3垂向比例曲線分析

為了更好的研究目標(biāo)層位，可以將目標(biāo)層位進(jìn)行細(xì)分。本文將目的層細(xì)分為10個(gè)小層。垂向比例曲線可以在垂向上確定每個(gè)小層中各微相的含量。垂向比例曲線可以直接通過(guò)測(cè)井?dāng)?shù)據(jù)計(jì)算得到，圖2即為O-11a層通過(guò)測(cè)井?dāng)?shù)據(jù)得到的垂向比比例曲線。

圖2　O-11a層垂向比例疊加曲線Fig.2　Vertical proportion curve of O-11a layer

圖2中，泛濫平原微相代碼0，用shale表示，灰色。河道微相代碼1，用Sand表示，黃色。心灘微相代碼2，用Sand2表示，橘黃色。從圖2可以看出，O-11a層泛濫平原微相最多，河道和心灘較少，發(fā)育程度較低。泛濫平原微相在該層頂部和底部較多，在該層中部變化較為平緩。心灘在O-11a層的變化比較平穩(wěn)。

垂向比例曲線可以很好地指導(dǎo)訓(xùn)練圖像的繪制。目標(biāo)層在垂向共分為10片，那么相對(duì)應(yīng)的訓(xùn)練圖像也為10片。訓(xùn)練圖像中每小片上3種沉積微相的含量應(yīng)該同垂向比例曲線中對(duì)應(yīng)的相同。

1.4O-11a層訓(xùn)練圖像

訓(xùn)練圖像中的模式是目標(biāo)儲(chǔ)層中反復(fù)出現(xiàn)的特征。在決定將一個(gè)訓(xùn)練圖像作為先驗(yàn)?zāi)Ｐ陀糜诙帱c(diǎn)統(tǒng)計(jì)時(shí)，要判斷其是否滿足平穩(wěn)性的要求。作為訓(xùn)練圖像的另一個(gè)特點(diǎn)就是必須相對(duì)簡(jiǎn)單，也就是訓(xùn)練圖像的構(gòu)型不能太復(fù)雜。

本文利用砂體厚度(圖3)和地震反演波阻抗數(shù)據(jù)(圖4)做出建模所需的2種訓(xùn)練圖像。由于主觀認(rèn)識(shí)，這2組訓(xùn)練圖像會(huì)有一定的不同。這為研究不同方法造成的不確定性提供了幫助。

圖3　O-11a層砂體厚度Fig.3　Sand-body thickness of O-11a layer

1.4.1第1種訓(xùn)練圖像根據(jù)19口直井在O-11a層的砂體數(shù)據(jù)，做出O-11a層的砂體厚度(圖3)。在辮狀河中，河道厚度一般小于心灘的厚度，這為繪制訓(xùn)練圖像提供了依據(jù)。從圖3可以看出，研究區(qū)內(nèi)有幾處砂體最厚的地方(圖3中黃色區(qū)域)，這幾處對(duì)應(yīng)了辮狀河道中的心灘。區(qū)內(nèi)砂體較厚的地方(圖3中淺綠色區(qū)域)對(duì)應(yīng)了河道微相，同時(shí)河道砂體具有很好的連續(xù)性。在研究區(qū)西北部砂巖較薄(圖3中深藍(lán)色、藍(lán)色區(qū)域)，該處對(duì)應(yīng)了泛濫平原。這符合辮狀河沉積特點(diǎn)，可以根據(jù)砂體厚度繪制訓(xùn)練圖像。因此，在根據(jù)砂體厚度做訓(xùn)練圖像時(shí)，可將砂體厚度最大的地方劃為心灘，厚度小的地方劃為泛濫平原，中間部分劃為河道。

圖4　地震反演波阻抗分布Fig.4　Distribution of seismic inversion wave impedance

圖5　O-11a層訓(xùn)練圖像Fig.5　Training images of O-11a layer

從圖5(a)的10片訓(xùn)練圖像可以看出，區(qū)內(nèi)3種沉積微相的分布較穩(wěn)定，在橫向和垂向僅有很小的變化。3種微相在平面上位置的變化，反映了3種微相在橫向上的遷移，在垂向具有很好的連續(xù)性，反應(yīng)了沉積微相的垂向加積。訓(xùn)練圖像很明顯地反映了辮狀河沉積特征，心灘發(fā)育在河道中，同時(shí)也反映了該區(qū)的地質(zhì)特征：河道、心灘發(fā)育較差，泛濫平原微相較為發(fā)育。因此是一組合理的訓(xùn)練圖像。1.4.2第2種訓(xùn)練圖像地震資料具有橫向分辨率高的特點(diǎn)，在橫向上可以連續(xù)追蹤地層信息，可以在較大尺度上反映地層構(gòu)造。通過(guò)地震資料反演的波阻抗數(shù)據(jù)，可以轉(zhuǎn)換為訓(xùn)練圖像。本文利用地震反演波阻抗分布(圖4)繪制了第2種訓(xùn)練圖像，由于篇幅原因，現(xiàn)選取地震波阻抗分布圖第1、4、7、10片。

從圖5中各圖對(duì)比可以看出，O-11a層波阻抗分布差別很小。

波阻抗數(shù)據(jù)與巖性的關(guān)系具有很強(qiáng)的區(qū)域特性[14]。為了利用地震波阻抗資料繪制訓(xùn)練圖像，需建立波阻抗數(shù)據(jù)與沉積微相的關(guān)系[15]。通常，泥巖表現(xiàn)出的波阻抗值較低，而砂巖則較高。通過(guò)波阻抗值數(shù)據(jù)與沉積微相的對(duì)應(yīng)關(guān)系，可以為訓(xùn)練圖像的繪制提供依據(jù)。

從波阻抗分布圖可以看出，在靠近研究區(qū)邊界的地方波阻抗值較低，在中央出現(xiàn)高值。低值與高值之間的連續(xù)性較好，它們之間也出現(xiàn)了明顯的過(guò)度區(qū)域(圖4中灰色部分)。

觀察圖4可以看出，心灘和河道微相對(duì)應(yīng)較高的波阻抗值(灰色、紅色和黃色指示區(qū)域)，泛濫平原對(duì)應(yīng)較低的波阻抗值(淺藍(lán)色區(qū)域)，且3種微相在波阻抗數(shù)據(jù)上的分布較明顯。因此，在繪制訓(xùn)練圖像可以將波阻抗值低的區(qū)域劃為泛濫平原微相(淺藍(lán)色區(qū)域)，阻抗值最高的區(qū)域(黃色、紅色區(qū)域)劃為心灘微相，介于兩者之間的區(qū)域(灰色區(qū)域)劃為河道微相。從波阻抗分布圖也可以看出河道為北西南東向，同時(shí)考慮河道的連續(xù)性，作出圖5(b)所示的第2種訓(xùn)練圖像。

為了研究波阻抗值與沉積微相的對(duì)應(yīng)關(guān)系，選取連井剖面觀察波阻抗值與井上微相的對(duì)應(yīng)關(guān)系(圖6)。圖6中的測(cè)井曲線為自然伽馬曲線。

圖6　波阻抗數(shù)據(jù)剖面Fig.6　Wave impedance data section

2　地震約束的儲(chǔ)層建模

近年來(lái)，地震檢測(cè)水平有了很大的提高。目前地震資料在儲(chǔ)層橫向展布預(yù)測(cè)及油藏精細(xì)描述中占有非常重要的地位，對(duì)確定剩余油的分布也極其重要[16]。

為獲得砂體概率曲線，需要對(duì)研究區(qū)微相與波阻抗數(shù)據(jù)的關(guān)系進(jìn)行分析[17]。統(tǒng)計(jì)沉積微相在自然伽馬-波阻抗交匯圖上的分布[18-19]，可分析波阻抗與沉積微相之間的關(guān)系，如圖7所示。

圖7　沉積微相在自然伽馬-波阻抗交匯圖上的分布Fig.7　Distribution of sedimentary microfacies in natural gamma-wave impedance intersection diagram

3種微相在波阻抗數(shù)據(jù)上的分布具有很大的不同。泛濫平原—河道—心灘的波阻抗值由低到高分布。

從波阻抗數(shù)據(jù)的分布可以看出，波阻抗值主要分布在(3 300～7 800)×106g/(m2·s)之間，高阻抗值主要集中在(4 000～5 200)×106g/(m2·s)。通過(guò)分析對(duì)比確定了泛濫平原與河道的波阻抗分界值為4 650×106g/(m2·s)，河道與心灘的波阻抗分界值為5 650×106g/(m2·s)。因此，得到了3種微相的砂體概率生成曲線，如圖8所示。

砂體概率生成曲線的含義是指在某一個(gè)指定波阻抗值條件下單一微相出現(xiàn)的概率。砂體概率曲線的意義就是3種微相在空間上的概率分布。需要指出的是，砂體概率生成曲線必須滿足數(shù)學(xué)定義，因此3種微相的砂體概率生成曲線在同一個(gè)波阻抗值的和必須為1。根據(jù)3種微相的概率生成曲線可以得到3種微相的空間概率分布，得到3種微相的空間概率分布即可作為軟數(shù)據(jù)用來(lái)約束建模。圖9即為3種沉積微相空間概率分布圖。

僅用測(cè)井?dāng)?shù)據(jù)建模，使用2種訓(xùn)練圖像得到的模擬結(jié)果差異很大，說(shuō)明僅在測(cè)井?dāng)?shù)據(jù)控制下，模擬結(jié)果具有較大的不確定性[20]。為了減少建模過(guò)程中的不確定性，采用地震數(shù)據(jù)約束建模。

在地震數(shù)據(jù)約束下，利用2種訓(xùn)練圖像建模。選用井震比為1:1(即建模過(guò)程中測(cè)井?dāng)?shù)據(jù)和地震數(shù)據(jù)的權(quán)重比相等)，隨機(jī)種子為211175，建模結(jié)果對(duì)比如圖10、圖11所示。

從圖10和圖11可以看出，在沒(méi)有地震約束的情況下，由于2組訓(xùn)練圖像不同，其模擬結(jié)果中各沉積微相在橫向上的分布具有較大的差異。在地震數(shù)據(jù)約束下，模擬結(jié)果的分布和形態(tài)除小部分有差異外基本一致，同時(shí)泛濫平原、河道和心灘微相的分布形態(tài)比較穩(wěn)定。

圖8　3種沉積微相概率生成曲線Fig.8　Sandbody probability generating curves of three kinds of microfacies

圖9　各沉積微相空間概率分布Fig.9　Spatial probability distribution of every sedimentary microfacies

圖10　有無(wú)地震約束的模擬結(jié)果對(duì)比(第1種訓(xùn)練圖像)Fig.10　Simulation results of seismic constraints and non-seismic constraints(1st training images)

圖11　有無(wú)地震約束的模擬結(jié)果對(duì)比(第2種訓(xùn)練圖像)Fig.11　Simulation results of seismic constraints and non-seismic constraints (2nd training images)

對(duì)比結(jié)果表明，地震數(shù)據(jù)約束下可以降低建模中訓(xùn)練圖像對(duì)結(jié)果的影響，不確定性明顯降低。這說(shuō)明地震數(shù)據(jù)約束建模的重要性。

另外分析了地震約束后，兩組訓(xùn)練圖像模擬結(jié)果中各微相的含量，如圖12所示。

圖12　地震約束下兩組訓(xùn)練圖像微相含量對(duì)比Fig.12　Comparison of the microseries content of two groups of training images under seismic constraints

圖中泛濫平原微相用Shale表示，河道微相用Sand表示，心灘微相用Sand2表示。藍(lán)色表示模擬結(jié)果中3種沉積微相的比例、紅色為原始數(shù)據(jù)中各沉積微相的比例。

從圖12可知，在地震數(shù)據(jù)的約束下，2組訓(xùn)練圖像的模擬結(jié)果與原始數(shù)據(jù)(紅色)非常接近，這說(shuō)明利用地震數(shù)據(jù)約束建模的結(jié)果是合理可信的。

本文選取2組訓(xùn)練圖像分別在地震數(shù)據(jù)約束前后的連井剖面進(jìn)行對(duì)比和分析。這種分析和對(duì)比，是以地震波阻抗剖面圖為依據(jù)的。剖面上3口井從左至右依次是CN-42井—CN-41井—CJS-1-0井。圖13中的測(cè)井曲線為自然伽馬曲線。

從圖13的4組剖面對(duì)比可以看出，在地震數(shù)據(jù)約束前，2組訓(xùn)練圖像模擬結(jié)果的CN-42井和CN-41井井間區(qū)域出現(xiàn)較多的心灘和河道微相，尤其是第2種訓(xùn)練圖像模擬結(jié)果中，河道微相大量出現(xiàn)，同時(shí)連續(xù)性很好。在CN-41井和CJS-1-0井井間區(qū)域出現(xiàn)大量泛濫平原微相，僅有少量河道和心灘微相出現(xiàn)。

地震數(shù)據(jù)約束后，2組訓(xùn)練圖像模擬結(jié)果的CN-42井和CN-41井井間區(qū)域出現(xiàn)大量泛濫平原微相，僅有很少的河道微相。在CN-41井和CJS-1-0井井間區(qū)域出現(xiàn)大量河道和心灘微相，同時(shí)都發(fā)育在該層的底部。

從波阻抗數(shù)據(jù)剖面可以看出，由于地震資料橫向上具有較高的分辨率，在CN-42井和CN-41井井間區(qū)域并沒(méi)有河道和心灘發(fā)育。在CN-41井和CJS-1-0井井間區(qū)域，該層有河道和心灘發(fā)育，并出現(xiàn)在該層底部。在地震數(shù)據(jù)約束后，就出現(xiàn)了上面這種現(xiàn)象。

從圖13中4組剖面也可以看出，井上數(shù)據(jù)與井附近模擬結(jié)果對(duì)應(yīng)較好。但在地震數(shù)據(jù)約束下得到兩組的剖面中，CJS-1-0井上數(shù)據(jù)與井附近對(duì)應(yīng)較差，主要原因是由于地震在垂向上分辨率很低造成的。

圖13　有無(wú)地震約束的模擬剖面對(duì)比Fig.13　Simulation sections of seismic constraints and non-seismic constraints

3　地震約束對(duì)模擬結(jié)果的影響

對(duì)比可知，是否有地震數(shù)據(jù)的約束，會(huì)造成模擬結(jié)果的差異，同時(shí)造成的這種差異也是不同的。不同的隨機(jī)種子導(dǎo)致的模擬結(jié)果也有較大的差異，河道形態(tài)、分布變化較大。各模擬實(shí)現(xiàn)之間缺少相似性。在地震數(shù)據(jù)約束下，不同隨機(jī)種子得到的結(jié)果具有較高的相似度。在不同模擬實(shí)現(xiàn)的同一片上，各微相在橫向上僅有很小的變化，這說(shuō)明地震數(shù)據(jù)在橫向上發(fā)揮了很強(qiáng)的控制作用，極大降低了不確定性。

利用2組訓(xùn)練圖像，在地震數(shù)據(jù)約束和沒(méi)有約束的情況下，計(jì)算得到了網(wǎng)格中微相一致的總網(wǎng)格數(shù)與模型總網(wǎng)格數(shù)比值分別為0.542 955、0.529 230、0.807 925、0.806 715。圖14為4個(gè)平均值的直方圖。

圖14　平均值直方圖Fig.14　Mean value histogram

圖14中，第1種代表僅用測(cè)井?dāng)?shù)據(jù)，利用第1種訓(xùn)練圖像得到的40組模擬結(jié)果，兩兩比較后，網(wǎng)格中微相一致的總網(wǎng)格數(shù)與模型總網(wǎng)格數(shù)比值的平均值。第2種代表僅用測(cè)井?dāng)?shù)據(jù)，利用第2種訓(xùn)練圖像得到的40組模擬結(jié)果，兩兩比較后，網(wǎng)格中微相一致的總網(wǎng)格數(shù)與模型總網(wǎng)格數(shù)比值的平均值。第1種(約束)代表在地震約束下，利用第1種訓(xùn)練圖像得到的40組模擬結(jié)果，兩兩比較后，網(wǎng)格中微相一致的總網(wǎng)格數(shù)與模型總網(wǎng)格數(shù)比值的平均值。第2種(約束)代表在地震約束下，利用第2種訓(xùn)練圖像得到的40組模擬結(jié)果，兩兩比較后，網(wǎng)格中微相一致的總網(wǎng)格數(shù)與模型總網(wǎng)格數(shù)比值的平均值。

從圖14可以看出，僅用測(cè)井?dāng)?shù)據(jù)建模，模擬結(jié)果之間具有較大的差異，說(shuō)明模擬實(shí)現(xiàn)具有較大的不確定性。在地震約束下，模擬結(jié)果之間具有較小的差異，說(shuō)明模擬實(shí)現(xiàn)的不確定性較小。

4　結(jié)　論

(1)僅依靠測(cè)井曲線約束的傳統(tǒng)建模方法所得到的結(jié)果不準(zhǔn)確，不確定性很大。

(2)訓(xùn)練圖像主要依賴(lài)人的主觀意識(shí)。由于主觀認(rèn)識(shí)的差異，針對(duì)同一區(qū)域，不同的人會(huì)做出不同的訓(xùn)練圖像，這就會(huì)導(dǎo)致模擬結(jié)果具有較大差別，引起較大的不確定性。

(3)在地震數(shù)據(jù)的約束下，不論是依靠砂體厚度建立的訓(xùn)練圖像或者是依靠波阻抗數(shù)據(jù)建立的訓(xùn)練圖像，模擬結(jié)果與實(shí)際情況接近，說(shuō)明地震數(shù)據(jù)能有效降低儲(chǔ)層建模中的不確定性。這為建模方法的選擇及儲(chǔ)層的后續(xù)開(kāi)發(fā)提供了很好的借鑒。

[1]ARAKTINGI U G,BASHORE W M .Effects of properties in seismic data on reservoir characterization and consequent fluid-flow predictions when integrated with well logs[C].SPE 24752,1992:913-926 .

[2]YANG C T,CHOPRA A K,CHU J,et al.Integrated geostatistical reservoir description using petrophysical,geological,and seismic data for Yacheng 13-1 Gasfield[C].SPE 30566,1995:357-372 .

[3]戴危艷，李少華，譙嘉翼，等.儲(chǔ)層不確定性建模研究進(jìn)展[J].巖性油氣藏，2015,27(4)：127-133.

DAI Weiyan,LI Shaohua,QIAO Jiayi,et al.Progress of reservoir uncertainty modeling[J].Lithologic Reservoir,2015,27(4):127-133.

[4]BORTOLI L J，ALABERT F，HAAS A，et al.Constraining Stochastic Images to Seismic Data[M].Troia:International Geostatistics Congress，1992:1-50.

[5]JOURNEL A G.Geostatistics:roadblocks and challenges[C]//Soares A .Geostatistics-Troia.Dordrecht:Kluwer Academic，1992:475-489.

[6]GUARDIANO F，SRIVASTAVA R M.Multivariate geostatistics:beyond bivariate moments[C]// Soares A.Geostatistics-Troia.Dordrecht:Kluwer Academic，1993:133-144.

[7]STREBELLE S.Conditional simulation of complex geological structures using multiple-point statistics[C].Mathmatical Geology，2001，34(1):1-21

[8]STREBELLE S，JOURNEL A G.Reservoir modeling using multiple-point statistics[C].SPE 71324，2001.

[9]劉超，謝傳禮，Yannick Tepinhi，等.多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)在儲(chǔ)層相建模中的應(yīng)用[J].斷塊油氣田，2015,22(6)：760-764.

LIU Chao,XIE Chuanli,Yannick Tepinhi,et al.Application of multiple-point geostatistics in reservoir facies modeling[J].Fault-Block Oil & Gas Field,2015,22(6):760-764.

[10] 吳小軍，李曉梅，謝丹，等.多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)方法在沖積扇構(gòu)型建模中的應(yīng)用[J].巖性油氣藏，2015,27(5)：87-91.

WU Xiaojun,LI Xiaomei,XIE Dan,et al.Application of multiple-point geostatistics method to structure modeling of alluvial fan[J].Lithologic Reservoirs,2015,27(5):87-91.

[11] 侯君，戴國(guó)汗，危杰，等.委內(nèi)瑞拉奧利諾科重油帶油藏特征及開(kāi)發(fā)潛力[J].石油實(shí)驗(yàn)地質(zhì)，2014,36(6)：725-730.

HOU Jun,DAI Guohan,WEI Jie,et al.Reservoir characterization and exploitation potential in Orinoco heavy oil belt in Venezuela[J].Petroleum Geology Experiment,2014,36(6):725-730.

[12] JOURNEL A G.Combining knowledge from diverse sources：An alternative to traditional data independence hypotheses[J].Mathematical Geology，2002,34(5)：573-596.

[13] 劉玲.大慶宋芳屯油田芳2區(qū)塊地震與地質(zhì)資料綜合儲(chǔ)層地質(zhì)建模研究[D].北京：中國(guó)地質(zhì)大學(xué)，2006.

LIU ling.A study on reservoir geological modeling with seismic and well-log data in Fang 2 area of Daqing Songfangtun Oilfield[D].Beijing:China University of Geosciences,2006.

[14] 周麗清，趙麗敏，趙國(guó)梁，等.高分辨率地震約束相建模[J].石油勘探與開(kāi)發(fā)，2002，29(3)：56-58.

ZHOU Liqing,ZHAO Limin,ZHAO Guoliang,et al.Stochastic facies modeling constrained by high-resolution seismic data[J].Petroleum Exploration and Development,2015,29(3):56-58.

[15] TJOLSEN C B，JOHNSEN Gyrid,HALVORSEN Gunnar，et al.Seismic data can improve the stochastic facies modeling[J].SPE Formation Evaluation,1996,11(3):141-146.

[16] 李桂亮.基于地質(zhì)統(tǒng)計(jì)學(xué)的地震約束建模方法研究[D].西安：西安石油大學(xué)，2010.

LI Guiliang.Reasearch of Method on Seismic Contained Modeling Based on Geostatistics[D].Xi'an:Xi'an Shiyou University,2010.

[17] 劉文嶺.地震約束儲(chǔ)層地質(zhì)建模技術(shù)[J].石油學(xué)報(bào)，2008，29(1)：64-68,74.

LIU Wenling.Geological modeling technique for reservoir constrained by seismic data[J].Acta Petrolei Sinica,2008,29(1):64-68,74.

[18] AVSETH P，TAPAN Mukerji.Seismic lithofacies classification from well logs using statistical rock physics[J].Petrophysics，2002，21：70-81.

[19] SRIVASTARA R M.Reservoir characterization with probability field simulation[C].SPE 24753,1992.

[20] 陳更新，趙凡，曹正林，等.地震反演、地質(zhì)協(xié)同約束儲(chǔ)層精細(xì)建模研究[J].天然氣地球科學(xué)，2014，25(11)：1839-1846.

CHEN Gengxin,ZHAO Fan,CAO Zhenglin,et al.Integrated reservoir modeling based on seismic inverdion and geological data[J].Natural Gas Geoscience,2014,25(11):1839-1846.

責(zé)任編輯：張新寶

Reducing Uncertainty of Reservoir Modeling by Seismic Constraints:Taking O-11a Super-heavy Oil Reservoir in MPE3 Block of Venezuela as an Example

LI Miao1,GENG Yi2,LI Xuan3

(1.The 5th Oil Production Plant,Changqing Oilfield Company,Xi'an 710200,Shaanxi,China;2.Lanzhou Documentation and Information Centre,Chinese Academy of Sciences,Lanzhou 730030,Gansu,China;3.The 2nd Oil Production Plant,Changqing Oilfield Company,Qingyang 735101,Gansu,China)

The sand-body distribution in Venezuela's Orinoco heavy oil belt is complex and the reservoir heterogeneity in this area is strong,and therefore the traditional two-step modeling method is difficult to obtain reliable results.The uncertainty of the simulation results of the sedimentary microfacies in the study area is analyzed by using the characterization of training images and the probability curves of the sand-bodies based on the seismic data;two groups of sedimentary microfacies models with and without the constraint of seismic data were established using multi-point geological statistics,and the comparison of their results shows that the sand-body distribution obtained by the sedimentary microfacies model with the constraint of seismic data is consistent with seismic data,and the uncertainty of the sedimentary microfacies model is obviously reduced.The correct and reliable result of reservoir geological modeling with the constraint of seismic data has a good reference to the distribution of the follow-up horizontal wells and the exploration and development of the adjacent areas in the study area.

geological modeling;seismic data constraint;training image;uncertainty

2015-12-29

國(guó)家科技重大專(zhuān)項(xiàng)(編號(hào)：2011ZX05032-006)

李苗(1989-)，女，助理工程師，碩士，主要從事石油地質(zhì)及油氣田開(kāi)發(fā)研究.E-mail:miaolinda@foxmail.com

10.3969/j.issn.1673-064X.2016.03.006

P618.130.2

1673-064X(2016)03-0039-09

A