丁 芳 段冬平 宋剛祥 劉英輝 陳 晨

(中海石油(中國)有限公司上海分公司 上海 200335)

東海盆地西湖凹陷花港組發(fā)育河流相和河流—三角洲相，其中河流相儲層是油氣聚集最多的儲層，因此開展河流相儲層精細(xì)研究，準(zhǔn)確預(yù)測河流形態(tài)及其在三維空間的分布，建立精細(xì)的儲層地質(zhì)模型，具有重要的現(xiàn)實(shí)意義[1-2]?；谙笤突谀繕?biāo)建模方法是現(xiàn)在比較流行的儲層表征方法。基于象元的序貫指示建模方法主要是利用兩點(diǎn)間變差函數(shù)來反映儲層結(jié)構(gòu)的變異性，由于海上油田開發(fā)井?dāng)?shù)少，井距大，存在變差函數(shù)求不準(zhǔn)問題，特別是對于幾何形態(tài)特征復(fù)雜的河道砂體，基于象元的隨機(jī)建模方法無法很好地表征出砂體的接觸關(guān)系。而基于目標(biāo)建模方法可以回避統(tǒng)計(jì)學(xué)空間關(guān)系特征，利用一些先驗(yàn)地質(zhì)認(rèn)識作為條件信息加入到模型中去，這樣可最大限度地綜合地質(zhì)家的認(rèn)識[3-4]，因此，基于目標(biāo)建模方法得到了廣泛應(yīng)用。

目前關(guān)于河道物性的建模方法較少，主要有局部變化均值方法、多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)方法，這些方法可實(shí)現(xiàn)河道物性中心比邊部好的地質(zhì)效果，在一定程度上提高了建模結(jié)果的合理性，但是實(shí)現(xiàn)過程較繁雜且僅針對單個(gè)方向進(jìn)行了約束，在井距小的陸上油田實(shí)現(xiàn)效果可能較好，但對于井?dāng)?shù)少的海上油田，這些方法效果不理想[5-6]。在前人研究的基礎(chǔ)上，本文首先通過單井河道期次精細(xì)劃分，采用分級和基于目標(biāo)建模方法建立河道相，然后在相控建模的基礎(chǔ)上，考慮多方向約束，基于平面分布趨勢、垂向沉積韻律趨勢對物性進(jìn)行模擬，同時(shí)發(fā)揮地震資料平面分辨率高的特點(diǎn)，融合了泊松阻尼因子屬性的趨勢，實(shí)現(xiàn)了河道相特殊的地質(zhì)效果。

1 研究區(qū)概況

D油田位于東海陸架盆地浙東坳陷西湖凹陷中央背斜帶南部，經(jīng)歷過多次大的構(gòu)造運(yùn)動，為斷陷—斷坳結(jié)構(gòu)基礎(chǔ)上的反轉(zhuǎn)背斜構(gòu)造，D油田漸新統(tǒng)花港組上段(E3h1)3層分為3a、3b和3c共3小層，其中3b和3c以辮狀河道沉積為主，砂體的連續(xù)性和物性總體較好，孔隙度為15%～30%，滲透率為50～1 000 mD，屬中孔中滲儲層。2014年9月該油田進(jìn)入投產(chǎn)階段3b層是重點(diǎn)開發(fā)層位，目前該層有4口水平井同時(shí)生產(chǎn)，2015年6月水平井D2H井快速見水，2016年2月達(dá)到水淹狀態(tài)，2015年10月水平井D1H井也開始產(chǎn)水。目前該油田已鉆井11口(含4口水平井)，隨著開發(fā)的不斷深入，生產(chǎn)矛盾越來越突出，儲層局部非均質(zhì)性差異對開發(fā)的影響也日益明顯。

2 多期河道疊置的厚砂體刻畫

2.1 井點(diǎn)河道期次劃分

通過上述4種方法，把D油田花上段3層3b和3c小層的河道劃分出7期河道(圖1)，復(fù)合砂體內(nèi)部夾層發(fā)育不均勻，連井剖面反映出研究區(qū)井順河道分布的特點(diǎn)，單河道砂體厚度均勻，無法反映河道邊界信息。由于3c層儲層為油水同層或水層，本文不作重點(diǎn)闡述。

2.2 河道砂體表征

由于研究區(qū)鉆井少，且順河道方向，河道邊界刻畫存在很大的不確定性，加之受地震資料分辨率的影響，單砂體無法在地震資料上很好地刻畫出來，因此本次研究中河道參數(shù)參考了相似條件下的露頭資料——山西柳林地區(qū)辮狀河露頭[16]，同時(shí)參考巖心和地震資料等多種數(shù)據(jù)，確定河道寬度范圍為200～1 000 m。根據(jù)現(xiàn)代沉積和經(jīng)驗(yàn)值給定振幅和波長，振幅范圍為500～800 m，波長范圍為1 000～2 500 m。

本次相模型采用分級和基于目標(biāo)建模方法，首先根據(jù)地震解釋的砂泥分布，建立巖相模型，再在砂巖中根據(jù)井上單期河道的劃分結(jié)果，在模型中建立3b小層砂體的連通關(guān)系和垂向上河道的概率曲線，設(shè)置河道的參數(shù)，通過反復(fù)調(diào)整河道參數(shù)與人機(jī)交互，最終得到與研究區(qū)相符的沉積相模型，如圖2所示。從模擬結(jié)果來看，基于目標(biāo)建模方法能完全忠實(shí)于沉積來源，并且保證了模擬的單河道沿河道帶的延伸方向展布。

表1東海D油田單一期次河道識別的4個(gè)標(biāo)志
Table1FourmarksofidentificationofasinglechannelinDoilfield,theEastChinaSea

圖1 東海D油田小層河道期次劃分結(jié)果Fig.1 Classification results of channel identification of layers of D oilfield,the East China Sea

圖2 東海D油田小層河道模擬結(jié)果Fig.2 Simulation results of channel of layer of D oilfield，the East China Sea

3 多源趨勢融合建模方法

河道砂體地下情況復(fù)雜，采用單一趨勢約束建模不能完全體現(xiàn)復(fù)雜的非均質(zhì)性，本文應(yīng)用多源趨勢融合的方法建模。

1) 平面趨勢變換。平面趨勢是覆蓋建模工區(qū)范圍的趨勢面，趨勢面的每個(gè)網(wǎng)格值即為當(dāng)前位置處的趨勢值。根據(jù)河道特殊的地質(zhì)現(xiàn)象，河道中心比兩側(cè)物性好，趨勢面產(chǎn)生的方式主要是根據(jù)河道中心軸線向兩側(cè)，隨距離變大，趨勢面值變小。圖3是建立平面趨勢的變化過程，首先建立地質(zhì)體的中心線，該處值為0，越遠(yuǎn)離該處，距離越大，趨勢值越小，通過數(shù)學(xué)變換，得到了河道平面的物性分布趨勢圖。

2) 垂向趨勢變換。考慮到河流特殊的地質(zhì)現(xiàn)象，物性受沉積作用影響大，物性的變化是沿垂向等時(shí)地層單位，而非垂直深度[17]。因此，在趨勢分析過程中從地層頂面起始網(wǎng)格開始，往下分析深度與物性的關(guān)系。以第6期河道孔隙度模擬為例說明，將7口井的孔隙度數(shù)據(jù)與垂向網(wǎng)格層序號散點(diǎn)作圖，如圖4所示，可以看出隨著垂向深度變大，孔隙度有變大的趨勢，這與河道底部物性好于頂部物性的規(guī)律相符；利用圖4中擬合公式Y(jié)=1.859 54X+3.089 19，得到了孔隙度垂向趨勢(圖4)。

圖3 平面趨勢變換過程Fig.3 Plane trend transformation process

3) 多趨勢融合。為避免建模軟件算法的不足，在相控物性建模過程中融合多種趨勢的控制作用，在沒有井點(diǎn)約束的位置也能很好地體現(xiàn)地質(zhì)上對工區(qū)的認(rèn)識，如整合基于平面位置分布趨勢和垂向沉積韻律趨勢獲取河道物性的特殊地質(zhì)意義。泊松阻尼因子，它是基于縱波阻抗和橫波阻抗的數(shù)學(xué)運(yùn)算得到的，本區(qū)泊松阻尼因子與孔隙度相關(guān)性較好，能較好地識別相對高孔砂巖區(qū)，因此在模擬過程中，也將泊松阻尼因子反演趨勢加入到模型中。首先將得到的平面趨勢和垂向趨勢做歸一化處理；然后2個(gè)趨勢體相乘得到了平面和垂向趨勢的合體；最后通過建模軟件Petrel的trend modeling功能，將平面和垂向趨勢的合體與地震孔隙度反演融合，設(shè)置2個(gè)屬性體的權(quán)重即可得到整合3個(gè)趨勢的約束體。這樣的多趨勢融合體使模型滿足地質(zhì)意義的前提下，與實(shí)際情況相符，地震屬性中孔隙度高值的區(qū)域在模型中也表現(xiàn)為高值(圖5)。圖6為過D5、D6、D1和D4井河道和孔隙度模擬圖，在單河道內(nèi)部，每一期河道孔隙度的分布基本都符合河道底部比頂部好，中心比邊部好的特殊物性分布規(guī)律。

圖4 垂向趨勢變換過程Fig.4 Vertical trend transformation process

圖5 東海D油田小層第6期河道孔隙度模擬過程圖Fig.5 Porosity simulation process of the sixth channel of layer of D oilfield，the East China Sea

圖6 東海D油田小層河道期次與孔隙度模擬過井剖面圖Fig.6 Channel periods and porosity simulation diagram through-well profile of layer of D oilfield，the East China Sea

4 模型驗(yàn)證

驗(yàn)證模型的可靠性非常重要[18]。 從以下2個(gè)方面來驗(yàn)證。

二是從油田和油水井生產(chǎn)動態(tài)擬合精度方面。利用模型對單井的生產(chǎn)效果和氣油比進(jìn)行了擬合(圖7)，從擬合效果來看，準(zhǔn)確率較高，符合生產(chǎn)要求，模型預(yù)測結(jié)果可以為現(xiàn)今生產(chǎn)提供參考。

圖7 東海D油田開發(fā)效果擬合圖Fig.7 History matching of development effect of D oilfield， the East China Sea

5 結(jié)論

1) 在東海盆地西湖凹陷D油田花港組河流相儲層物性建模過程中，應(yīng)用基于平面位置分布趨勢、垂向沉積韻律趨勢和孔隙度敏感參數(shù)泊松阻尼因子反演屬性趨勢等多個(gè)級次融合，有效地刻畫了河道物性底部比頂部好，中部比邊部好的特殊地質(zhì)效果，降低了少井條件下隨機(jī)模擬方法的不確定性。

2) 模型驗(yàn)證結(jié)果表明，本文方法儲量擬合誤差在10%以內(nèi)，D油田4口水平井生產(chǎn)動態(tài)模擬結(jié)果和實(shí)際較為吻合，本文方法可在海上油田井少、井距大區(qū)域進(jìn)行推廣應(yīng)用。

[1] 徐安娜，穆龍新，裘懌楠，等.我國不同沉積類型儲集層中的儲量和可動剩余油分布規(guī)律[J].石油勘探與開發(fā)，1998，25(5):41-44.

XU Anna,MU Longxin,QIU Yinan,et al.Distribution pattern of OOIP and remaining mobile oil in different types of sedimentary reservoir of China[J].Petroleum Exploration and Development,1998,25(5):41-44.

[2] 安桂榮，許家峰，周文勝，等.海上復(fù)雜河流相水驅(qū)稠油油田井網(wǎng)優(yōu)化：以B-2 油田為例[J].中國海上油氣，2013，25(3)：28-31.

AN Guirong,XU Jiafeng,ZHOU Wensheng,et al.Well pattern optimization in offshore water-drive heavy oilfields with complicated fluvial facies[J].China Offshore Oil and Gas，2013,25(3)：28-31.

[3] 王家華，劉衛(wèi)麗，白軍衛(wèi)，等.基于目標(biāo)的隨機(jī)建模方法[J].重慶科技學(xué)院院報(bào)(自然科學(xué)版)，2012,14(1):162-163.

WANG Jiahua,LIU Weili,BAI Junwei,et al.On the object-based stochastic modeling method[J].Journal of Chongqing University of Science and Technology(Natural Sciences Edition)，2012,14(1):162-163.

[4] 吳勝和，金振奎，黃滄鈿，等.儲層建模[M].北京：石油工業(yè)出版社，1999.

[5] 高博禹，孫立春，胡光義，等.基于單砂體的河流相儲層地質(zhì)建模方法探討[J].中國海上油氣,2008，20(1):34-37.

GAO Boyu,SUN Lichun,HU Guangyi,et al.A discussion on a fluvial reservoir modeling method based on a single sandbody[J].China Offshore Oil and Gas,2008，20(1):34-37.

[6] 林承焰，陳仕臻，張憲國，等.多趨勢融合的概率體約束方法及其在儲層建模中的應(yīng)用[J].石油學(xué)報(bào)，2015,36(6):730-739.

LIN Chengyan,CHEN Shizhen,ZHANG Xianguo,et al.Probability constraint method based on multiple trend integration and its application in reservoir medeling[J].Acta Petrolei Sinica，2015,36(6):730-739.

[7] 葛云龍,逯徑鐵,廖保方,等.辮狀河相儲集層地質(zhì)模型:“泛連通體”[J].石油勘探與開發(fā),1998,25(5):77-79.

GE Yunlong,LU Jingtie,LIAO Baofang,et al.A braided river reservoir geological model：“pan-communicated sandbody”[J].Petroleum Exploration and Development,1998,25(5):77-79.

[8] 李順明，宋新民，蔣有偉，等.高尚堡油田砂質(zhì)辮狀河儲集層構(gòu)型與剩余油分布[J].石油勘探與開發(fā)，2011，38(4):474-482.

LI Shunming,SONG Xinmin,JIANG Youwei,et al.Architecture and remaining oil distribution of the sandy braided river reservoir in the Gaoshangpu Oilfield[J].Petroleum Exploration and Development,2011，38(4):474-482.

[9] 金振奎，時(shí)曉章，何苗.單河道砂體的識別方法[J].新疆石油地質(zhì)，2010,31(6):572-575.

JIN Zhenkui,SHI Xiaozhang,HE miao.Identification methods for single-channel sand body[J].Xinjiang Petroleum Geology,2010,31(6):572-575.

[10] 胡光義，王加瑞，武士堯.利用地震分頻處理技術(shù)預(yù)測河流相儲層:基于精細(xì)儲層預(yù)測調(diào)整海上高含水油田開發(fā)方案實(shí)例[J].中國海上油氣，2005，17(4)：237-241.

HU Guangyi,WANG Jiarui,WU Shiyao.Prediction of fluvial-facies reservoirs using a technique of seismic frequency demultiplication:a case of adjusting development plan for an offshore oilfield with high water cut based on fine reservoir prediction[J].China Offshore Oil and Gas，2005，17(4)：237-241.

[11] ZHANG Jianning,HAN Wengong,YAN Zhaomin,et al.The geologic model and seimic forward response characteristic of the boundaries of single river channel[J].Petroleum Geophysics，2006，4(4)：32-35．

[12] LEEDER M R.Fluviatile fining upwards cycles and the magnitude of paleochannels[J].Geological Magazine,1973,110(3):265-276.

[13] 封從軍，鮑志東，單啟銅,等.三角洲平原復(fù)合分流河道內(nèi)部單砂體劃分:以扶余油田中區(qū)南部泉頭組四段為例[J].石油與天然氣地質(zhì),2012,33(1):77-83.

FENG Congjun,BAO Zhidong,SHAN Qitong,et al.Single sand body identification in compound distributary channel ofdelta plain:a case study from the fourth member of Quantou Formation in the southern part of central Fuyu oilfield[J].Oil & Gas Geology,2012,33(1):77-83.

[14] 何宇航,于開春.分流平原相復(fù)合砂體單一河道識別及效果分析[J].大慶石油地質(zhì)與開發(fā),2005,24(2):17-19.

HE Yuhang,YU Kaichun.Recognition and its effect analysis of single river channel in composite sand body with distributary plain facies[J].Petroleum Geology & Oilfield Development in Daqing,2005,24(2):17-19.

[15] 孫寶京.孤島油田特高含水期復(fù)雜非均質(zhì)油藏剩余油描述方法研究[D].東營：中國石油大學(xué)(華東)，2007.

SUN Baojing.Study on remaining oil description in complex heterogeneous reservoir at super-high water cut period of Gudao oilfield[D].Dongying:China University of Petroleum,2007.

[16] 金振奎，楊有星，尚建林，等.辮狀河砂體構(gòu)型及定量參數(shù)研究[J].天然氣地球科學(xué)，2014,25(3):311-317.

JIN Zhenkui,YANG Youxing,SHANG Jianlin,et al.Sandbody architecture and quantitative parameters of single channel sandbodies of braided river[J].Natural Gas Geoscience，2014,25(3):311-317.

[17] 吳勝和.儲層表征與建模[M].北京:石油工業(yè)出版社，2010:408-413.

[18] 丁芳，張金亮.復(fù)雜斷塊油藏隨機(jī)建模技術(shù)的應(yīng)用:以文留油田文25斷塊為例[J].北京師范大學(xué)學(xué)報(bào)(自然科學(xué)版)，2014，50(2):205-211.

DING Fang,ZHANG Jinliang.Application of stochastic modeling in complex fault block：case of Wen25 Block in Wenliu oilfield[J].Journal of Beijing Normal University(Natural Science)，2014，50(2):205-211.