海上復(fù)雜碎屑巖儲(chǔ)層油氣藏地質(zhì)建模關(guān)鍵技術(shù)*

葉小明 王鵬飛 霍春亮 高振南 徐 靜

(中海石油(中國(guó))有限公司天津分公司渤海石油研究院 天津 300459)

渤海油田在區(qū)域構(gòu)造上屬于渤海灣盆地的海域部分,在已探明和已開(kāi)發(fā)的油氣田中，各類復(fù)雜碎屑巖儲(chǔ)層油氣藏占有相當(dāng)大的比重，以河流相及淺水三角洲相油田為例，這兩類油氣藏的探明地質(zhì)儲(chǔ)量約占渤海油田總探明地質(zhì)儲(chǔ)量的60%[1]，因此實(shí)現(xiàn)復(fù)雜碎屑巖儲(chǔ)層油氣藏的高效開(kāi)發(fā)，對(duì)于渤海油田上產(chǎn)、穩(wěn)產(chǎn)意義重大。近年來(lái)的開(kāi)發(fā)實(shí)踐表明，在不同開(kāi)發(fā)階段，針對(duì)渤海海域的復(fù)雜碎屑巖儲(chǔ)層開(kāi)展精細(xì)地質(zhì)研究、充分把握影響油水運(yùn)動(dòng)規(guī)律的重要地質(zhì)特征，并建立能夠反映地質(zhì)認(rèn)識(shí)與滲流特征的儲(chǔ)層三維模型，是高效開(kāi)發(fā)此類油氣田的必要條件和重要基礎(chǔ)[2-4]。本文從不同開(kāi)發(fā)階段、不同儲(chǔ)層類型儲(chǔ)層表征及建模過(guò)程中的難點(diǎn)出發(fā)，提出了海上復(fù)雜碎屑巖儲(chǔ)層油氣藏地質(zhì)建模關(guān)鍵技術(shù)。

1 海上復(fù)雜碎屑巖儲(chǔ)層地質(zhì)建模面臨的挑戰(zhàn)

1.1 開(kāi)發(fā)早期地質(zhì)風(fēng)險(xiǎn)量化難度大

海上油田受高開(kāi)發(fā)成本的制約，井網(wǎng)密度低、資料相對(duì)少[4]，同時(shí)由于碎屑巖儲(chǔ)層本身較強(qiáng)的非均質(zhì)性及研究人員認(rèn)識(shí)的局限性，使得海上復(fù)雜碎屑巖油氣藏的開(kāi)發(fā)存在較大不確定性。這種不確定性在油田開(kāi)發(fā)早期表現(xiàn)得更為明顯，與鉆前相比，較多新投產(chǎn)油田鉆后認(rèn)識(shí)在油藏模式、構(gòu)造形態(tài)及儲(chǔ)量規(guī)模等方面都發(fā)生了較大的變化，因此只建立一套地質(zhì)模型的常規(guī)方法難以實(shí)現(xiàn)對(duì)地質(zhì)風(fēng)險(xiǎn)的量化表征。

1.2 中深層儲(chǔ)層精細(xì)表征難度大

隨著油田勘探開(kāi)發(fā)的深入，在渤海發(fā)現(xiàn)了大量中深層油田。這類油田儲(chǔ)層大多以三角洲相沉積為主，儲(chǔ)層橫向變化大，非均質(zhì)性強(qiáng)；同時(shí)，中深層儲(chǔ)層埋藏深，地震資料分辨率低，資料品質(zhì)難以滿足儲(chǔ)層預(yù)測(cè)的需求，給鉆前儲(chǔ)層定量化描述帶來(lái)了較大困難，增加了油田開(kāi)發(fā)的風(fēng)險(xiǎn)，使得開(kāi)發(fā)井失利增多。

1.3 儲(chǔ)層構(gòu)型建模難度大

隨著油田開(kāi)發(fā)的進(jìn)行，目前渤海各大主力油田相繼進(jìn)入中高含水階段，為進(jìn)一步挖潛剩余油，做好精細(xì)儲(chǔ)層表征很有必要。油田開(kāi)發(fā)實(shí)踐表明，諸如曲流河點(diǎn)壩內(nèi)部側(cè)積夾層等儲(chǔ)層內(nèi)部的小尺度構(gòu)型界面，對(duì)于油藏內(nèi)部的流體運(yùn)動(dòng)具有重要的影響，但如何將其定量表征到模型中還存在較大挑戰(zhàn)。國(guó)內(nèi)外學(xué)者在構(gòu)型建模方面也開(kāi)展了許多工作，但主要是基于加密網(wǎng)格來(lái)進(jìn)行表征，由于模型網(wǎng)格規(guī)模過(guò)大，在實(shí)際應(yīng)用中存在數(shù)值模擬運(yùn)算耗時(shí)巨大等問(wèn)題[5-13]。

2 海上復(fù)雜碎屑巖儲(chǔ)層地質(zhì)建模關(guān)鍵技術(shù)

針對(duì)目前面臨的挑戰(zhàn)，以多專業(yè)聯(lián)合驅(qū)動(dòng)地質(zhì)建模為手段，以最大程度提高儲(chǔ)層預(yù)測(cè)精度為目標(biāo)，形成了海上復(fù)雜碎屑巖儲(chǔ)層建模關(guān)鍵技術(shù)。

2.1 基于模型的地質(zhì)風(fēng)險(xiǎn)定量評(píng)價(jià)技術(shù)

從儲(chǔ)層表征與建模的角度來(lái)看，油田開(kāi)發(fā)的地質(zhì)風(fēng)險(xiǎn)主要表現(xiàn)在儲(chǔ)量規(guī)模及儲(chǔ)層連通性兩個(gè)方面。地質(zhì)儲(chǔ)量的多少直接決定了油田開(kāi)發(fā)的規(guī)模與效益，可驅(qū)替體積的多少則反映了在同一井網(wǎng)下儲(chǔ)層連通性的好與壞[14]。影響儲(chǔ)量和可驅(qū)替體積的不確定性變量有地層體積、儲(chǔ)層體積、孔隙體積、有效體積和含油體積等，在地質(zhì)建模時(shí)通常用構(gòu)造層面位置、砂地比、變差函數(shù)、NTG下限值和流體界面來(lái)表征。

從各不確定性變量分析入手，系統(tǒng)開(kāi)展了模型不確定性分析(建模參數(shù)選取)、實(shí)驗(yàn)設(shè)計(jì)及模型實(shí)現(xiàn)、儲(chǔ)量和可驅(qū)替體積概率分布統(tǒng)計(jì)及模型排隊(duì)優(yōu)選等工作。

首先應(yīng)用概率統(tǒng)計(jì)的方法全面考察所有影響儲(chǔ)層地質(zhì)模型的因素，對(duì)不確定變量進(jìn)行顯著性分析，優(yōu)選出主要不確定性變量并量化其分布范圍，為建模參數(shù)選取提供依據(jù)。渤海BZ34-X油田主要含油層位為東營(yíng)組和沙河街組，主要為辮狀河三角洲沉積，其間發(fā)育大套火山巖，導(dǎo)致地震分辨率較低，構(gòu)造解釋存在較大不確定性；按照鉆井的進(jìn)度統(tǒng)計(jì)鉆前地震解釋預(yù)測(cè)誤差，發(fā)現(xiàn)深度預(yù)測(cè)誤差在-180～107 m。為了量化構(gòu)造層面位置不確定性范圍，根據(jù)井上統(tǒng)計(jì)的誤差值，模擬得到井間誤差范圍；然后將井上的誤差值校正為0，得到每個(gè)解釋層面誤差模擬圖；最后將基于所有鉆井分層標(biāo)定進(jìn)行解釋的構(gòu)造層面作為中方案，將誤差模擬圖加或減這套構(gòu)造層面深度得到的構(gòu)造層面作為高低方案(圖1)。砂地比、變差函數(shù)、NTG下限值和流體界面等參數(shù)不確定性分布范圍的確定，需要充分利用地震、測(cè)井、巖心等多種資料來(lái)進(jìn)行研究，并最大程度縮小不確定性變化范圍。

然后通過(guò)實(shí)驗(yàn)設(shè)計(jì)，選出代表性強(qiáng)的少數(shù)實(shí)驗(yàn)方案，最大限度地反映不確定分布范圍。表1為BZ34-X油田實(shí)驗(yàn)設(shè)計(jì)方案， 表中-1、0和1分別代表各不確定性變量的最小值、期望值和最大值。按照實(shí)驗(yàn)設(shè)計(jì)方案中對(duì)各不確定性變量的取值，建立了9套地質(zhì)模型，并開(kāi)展流線模擬，即可計(jì)算得到各套模型的儲(chǔ)量和可驅(qū)替體積(表1)。

圖1 BZ34-X油田構(gòu)造方案示意圖Fig.1 Sketch map of structural scheme of BZ 34-X oilfield

表1 BZ 34-X油田實(shí)驗(yàn)設(shè)計(jì)方案Table 1 Experimental design scheme of BZ 34-X oilfield

依據(jù)各個(gè)模型計(jì)算的儲(chǔ)量與可驅(qū)替體積，應(yīng)用統(tǒng)計(jì)學(xué)響應(yīng)曲面法建立地質(zhì)變量與儲(chǔ)量及儲(chǔ)層連通性的響應(yīng)方程，為概率儲(chǔ)量及概率連通性計(jì)算提供數(shù)學(xué)模型。BZ34-X油田儲(chǔ)量響應(yīng)曲面方程為

STOIIP=3 638+335V1+128.75V2+53V3+

429.75V4+200V5-384V1V2-2.25V2V3-135.25V2V5

(1)

式(1)中：STOIIP為地質(zhì)儲(chǔ)量；V1為構(gòu)造層面位置；V2為砂地比；V3為變差函數(shù)；V4為有效厚度下限值；V5為油水界面。

利用響應(yīng)曲面公式，即可通過(guò)蒙特卡洛方法計(jì)算概率儲(chǔ)量及概率可驅(qū)替體積。用概率儲(chǔ)量和概率可驅(qū)替體積繪制交會(huì)圖，將實(shí)驗(yàn)設(shè)計(jì)中每個(gè)模型對(duì)應(yīng)的儲(chǔ)量及可驅(qū)替體積投到交會(huì)圖上，即可選取出具有一定概率分布的最大、期望和最小值的模型作為樂(lè)觀、期望和悲觀模型。

對(duì)BZ34-X油田地質(zhì)風(fēng)險(xiǎn)進(jìn)行了定量評(píng)價(jià)，利用選取的3套地質(zhì)模型開(kāi)展油藏?cái)?shù)值模擬，評(píng)價(jià)開(kāi)發(fā)指標(biāo)的差異。通過(guò)評(píng)價(jià)得出該油田最大累產(chǎn)為1 250×104m3，最小累產(chǎn)為782×104m3?；诖耍俳Y(jié)合后期經(jīng)濟(jì)評(píng)價(jià)等工作，便可初步判別油田開(kāi)發(fā)是否有經(jīng)濟(jì)效益，減小了油田開(kāi)發(fā)過(guò)程中的風(fēng)險(xiǎn)。

2.2 基于沉積演化模擬的中深層儲(chǔ)層定量表征技術(shù)

沉積演化數(shù)值模擬起源于20世紀(jì)60年代[15]，近年來(lái)主要在勘探階段應(yīng)用較多[16-17]，本文將其與儲(chǔ)層隨機(jī)建模相結(jié)合，實(shí)現(xiàn)對(duì)中深層儲(chǔ)層的定量表征。該方法主要包含以下幾步：①通過(guò)沉積物沉積過(guò)程參數(shù)定量化(包括可容納空間、物源供應(yīng)和沉積搬運(yùn)參數(shù)等)[17]，還原沉積物的形成過(guò)程，建立初始沉積演化模型；②用實(shí)鉆井的巖性分布特征來(lái)進(jìn)行驗(yàn)證，在地質(zhì)認(rèn)識(shí)合理范圍內(nèi)調(diào)整沉積演化模擬輸入?yún)?shù)，通過(guò)多次迭代使得模擬結(jié)果與井上特征基本吻合，得到最終的沉積演化模型；③通過(guò)網(wǎng)格等比例剖分、井點(diǎn)數(shù)據(jù)標(biāo)定，將沉積演化模型轉(zhuǎn)換為三維地質(zhì)建模約束趨勢(shì)體；④在趨勢(shì)體約束下建立儲(chǔ)層巖相及屬性模型，實(shí)現(xiàn)中深層儲(chǔ)層的定量化表征。

系統(tǒng)開(kāi)展了BZ34-X油田單井巖石學(xué)及古生物學(xué)特征分析、沉積環(huán)境研究、可容空間變化規(guī)律研究、物源模式及供應(yīng)速率研究、顆粒搬運(yùn)方式研究，量化了沉積物沉積過(guò)程參數(shù)。其中，物源供應(yīng)速率確定方法是：根據(jù)地層厚度、沉積相展布、沉積時(shí)間估算物源供給總速率，依據(jù)單井的厚度和巖性進(jìn)行劈分，得到各個(gè)物源流入點(diǎn)的供給速率和巖性比例。

在基底初始水深、可容空間變化、物源供給和搬運(yùn)方式研究的基礎(chǔ)上，進(jìn)行儲(chǔ)層沉積演化模擬研究，將模擬結(jié)果與單井的厚度和巖性進(jìn)行對(duì)比，調(diào)整模擬參數(shù)直至兩者近似一致，最終獲得經(jīng)過(guò)校正的沉積演化模型(圖2)。經(jīng)過(guò)校正的沉積演化模型較好地模擬了砂泥巖在三維空間的分布特征，但由于其在模擬過(guò)程中并不直接使用井?dāng)?shù)據(jù)，模擬結(jié)果在井點(diǎn)處難以達(dá)到完全匹配，因此需要將沉積演化模擬與傳統(tǒng)的建模方法相結(jié)合，以沉積模擬的砂泥巖分布結(jié)果作為約束條件，在井點(diǎn)數(shù)據(jù)的控制下進(jìn)行常規(guī)的地質(zhì)統(tǒng)計(jì)學(xué)建模。實(shí)際操作中，需要將沉積演化模型輸入到地質(zhì)建模軟件中，將其轉(zhuǎn)化為約束數(shù)據(jù),通過(guò)協(xié)同模擬的方法來(lái)建立砂泥巖巖相模型，然后通過(guò)相控逐步建立孔隙度、滲透率等物性參數(shù)模型。圖3為基于本文方法建立的模型預(yù)測(cè)的砂泥巖分布與Well7井鉆后實(shí)際結(jié)果的對(duì)比，可以看出，模型預(yù)測(cè)結(jié)果在砂體垂向發(fā)育部位及厚度都與實(shí)鉆結(jié)果具有較好的一致性。

圖2 BZ34-X油田砂巖分布模型三維圖Fig.2 Sandstone distribution model of BZ34-X oilfield

圖3 BZ34-X油田砂巖分布模型預(yù)測(cè)與實(shí)鉆結(jié)果對(duì)比圖Fig.3 Comparison of sandstone distribution between model prediction and actual drilling results of BZ 34-X oilfield

2.3 基于油藏?cái)?shù)值模擬動(dòng)態(tài)響應(yīng)的構(gòu)型界面定量表征技術(shù)

針對(duì)油田開(kāi)發(fā)中后期儲(chǔ)層構(gòu)型建模難度大的問(wèn)題，提出了基于油藏?cái)?shù)值模擬動(dòng)態(tài)響應(yīng)的構(gòu)型界面定量表征技術(shù),其方法步驟為：①基于數(shù)模動(dòng)態(tài)響應(yīng)確定需要進(jìn)行構(gòu)型界面定量表征的井區(qū)(范圍)；②小尺度構(gòu)型界面定量表征；③基于數(shù)模動(dòng)態(tài)響應(yīng)確定構(gòu)型界面對(duì)流體滲流影響的大小。

研究過(guò)程中，首先在前期地質(zhì)認(rèn)識(shí)基礎(chǔ)上，利用常規(guī)建模方法建立微相級(jí)別地質(zhì)模型，基于該模型進(jìn)行油藏動(dòng)態(tài)歷史擬合。將完成歷史擬合后修改的物性參數(shù)與初始模型參數(shù)進(jìn)行求差，開(kāi)展油藏模型動(dòng)態(tài)響應(yīng)的地質(zhì)因素判別與分析，識(shí)別由于小尺度構(gòu)型界面的存在影響流體運(yùn)動(dòng)的井區(qū)及層段，將其確定為開(kāi)展構(gòu)型建模目標(biāo)區(qū)域。這樣便避免了依據(jù)常規(guī)構(gòu)型分布模式來(lái)建立全區(qū)分布的構(gòu)型模型，而只是對(duì)流體滲流有影響的小尺度構(gòu)型界面進(jìn)行定量表征，極大地減小了研究工作量。由于地質(zhì)研究的不確定性，研究過(guò)程中還需要結(jié)合油藏動(dòng)態(tài)響應(yīng)來(lái)判別是否需要修正前期儲(chǔ)層構(gòu)型研究的成果，以確保表征到油藏模型中的構(gòu)型界面的合理性。

明確需要表征的區(qū)域后，便可以開(kāi)展小尺度構(gòu)型界面的定量表征工作。提出了一種等效表征方法，即將小尺度構(gòu)型界面對(duì)流體滲流的影響通過(guò)網(wǎng)格界面?zhèn)鲗?dǎo)率表征到油藏?cái)?shù)值模擬模型中去，而不是將其幾何參數(shù)反映到油藏?cái)?shù)值模擬模型中[2]。通過(guò)該方法，最終僅需提供給油藏?cái)?shù)值模擬人員一套符合數(shù)模數(shù)據(jù)規(guī)范的傳導(dǎo)率乘數(shù)數(shù)據(jù)卡，利用該數(shù)據(jù)卡實(shí)現(xiàn)在油藏模型中小尺度構(gòu)型界面的定量表征。

傳導(dǎo)率乘數(shù)數(shù)據(jù)卡給油藏人員提供了構(gòu)型界面在模型中的空間分布位置，其對(duì)流體滲流影響的大小(傳導(dǎo)率乘數(shù)數(shù)值大小)還需要基于數(shù)模動(dòng)態(tài)響應(yīng)來(lái)確定。研究中，需要根據(jù)注采井組的注采動(dòng)態(tài)通過(guò)井組之間的歷史擬合，來(lái)確定傳導(dǎo)率乘數(shù)的初始值，然后再開(kāi)展歷史擬合，若無(wú)法很好地?cái)M合，則重新對(duì)傳導(dǎo)率乘數(shù)進(jìn)行調(diào)整，直到最終能很好擬合生產(chǎn)動(dòng)態(tài)為止。該方法實(shí)現(xiàn)了構(gòu)型界面在三維模型中的定量表征，使得油藏?cái)?shù)值模擬人員能更好地依據(jù)地質(zhì)認(rèn)識(shí)來(lái)開(kāi)展歷史擬合。

基于以上步驟，編制了一套油藏模型及地質(zhì)模型雙模耦合及等效表征軟件，在提高工作效率的同時(shí)也方便了技術(shù)的推廣使用。對(duì)QHD32-X油田點(diǎn)壩砂體內(nèi)部側(cè)積夾層進(jìn)行了三維模型定量表征，大幅提高了動(dòng)態(tài)歷史擬合符合率及井間剩余油分布預(yù)測(cè)精度。圖4a為該油田A03井的初始模型含水率與實(shí)際生產(chǎn)含水率曲線，圖4b為采用常規(guī)修改物性參數(shù)的方法來(lái)進(jìn)行歷史擬合的結(jié)果?？梢钥闯?，通過(guò)常規(guī)歷史擬合，A03井的含水率與實(shí)際生產(chǎn)含水率的吻合程度有了一定提高，但在開(kāi)發(fā)后期，含水率的趨勢(shì)不一致。圖4c為利用本文方法完成側(cè)積夾層等效表征后進(jìn)行歷史擬合的結(jié)果，可見(jiàn)A03井不僅保證了前中期含水率的高吻合度，后期含水率曲線下降趨勢(shì)也擬合得很好。

圖4 QHD 32-X油田A03井不同方法歷史擬合結(jié)果對(duì)比Fig.4 History matching results contrast of Well A03 in QHD 32-X oilfield

3 應(yīng)用效果

本文建模關(guān)鍵技術(shù)在渤海8個(gè)油田進(jìn)行了應(yīng)用，在油田開(kāi)發(fā)方案編制及實(shí)施、綜合調(diào)整等方面成效顯著，3年間已見(jiàn)效增油量467.70×104m3。

基于模型的地質(zhì)風(fēng)險(xiǎn)定量評(píng)價(jià)技術(shù)從儲(chǔ)量規(guī)模與儲(chǔ)層連通性進(jìn)行了儲(chǔ)層定量地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)，在BZ25-XS、KL10-X及BZ34-X等油田開(kāi)發(fā)方案編制及實(shí)施中為開(kāi)發(fā)決策及方案優(yōu)化提供了可靠依據(jù)，累積動(dòng)用原油探明儲(chǔ)量1.8×108m3。

應(yīng)用基于沉積演化模擬的中深層儲(chǔ)層定量表征技術(shù)，預(yù)測(cè)儲(chǔ)層鉆遇率得到了顯著提高，在勘探開(kāi)發(fā)一體化背景下，成功指導(dǎo)了BZ34-X油田勘探評(píng)價(jià)井位部署、儲(chǔ)量評(píng)價(jià)及KL10-X、CFD6-X油田等多個(gè)油田開(kāi)發(fā)方案研究。

基于油藏?cái)?shù)值模擬動(dòng)態(tài)響應(yīng)的構(gòu)型界面定量表征技術(shù)，實(shí)現(xiàn)了油田開(kāi)發(fā)中后期儲(chǔ)層構(gòu)型建模，提高了儲(chǔ)層建模精度及準(zhǔn)確度，模型預(yù)測(cè)水淹程度與鉆后測(cè)井解釋符合率達(dá)到90%以上，成功指導(dǎo)了SZ36-X、QHD32-X及JZ9-X等油田236口調(diào)整井的部署實(shí)施，實(shí)現(xiàn)了油田高產(chǎn)高效開(kāi)發(fā)。

4 結(jié)論

提出了海上復(fù)雜碎屑巖儲(chǔ)層油氣藏地質(zhì)建模技術(shù)：基于模型的地質(zhì)風(fēng)險(xiǎn)定量評(píng)價(jià)技術(shù)，為油田開(kāi)發(fā)風(fēng)險(xiǎn)決策及優(yōu)化開(kāi)發(fā)方案提供了量化依據(jù)；基于沉積演化模擬的中深層儲(chǔ)層定量表征技術(shù)，大幅提高了中深層儲(chǔ)層預(yù)測(cè)精度；基于油藏?cái)?shù)值模擬動(dòng)態(tài)響應(yīng)的構(gòu)型界面定量表征技術(shù)，解決了小尺度構(gòu)型界面在油藏模型中的表征難題，提高了模型準(zhǔn)確度。渤海8個(gè)油田應(yīng)用表明，本文建模技術(shù)實(shí)用有效，3年間已見(jiàn)效增油量467.70×104m3。今后將逐步推廣應(yīng)用于渤海多個(gè)復(fù)雜碎屑巖油氣田，為“十三五”期間渤海油田年均超過(guò)億噸新增石油天然氣探明儲(chǔ)量的方案編制提供技術(shù)支撐，指導(dǎo)渤海近50個(gè)在生產(chǎn)油田及多個(gè)在建油田的高效開(kāi)發(fā)。

[1] 郭太現(xiàn)，楊慶紅，黃凱，等．海上河流相油田高效開(kāi)發(fā)技術(shù)[J]．石油勘探與開(kāi)發(fā)，2013， 40(6):708-714．

GUO Taixian,YANG Qinghong,HUANG Kai，et al．Techniques for high-efficient development of offshore fluvial oilfields[J].Petroleum Exploration and Development，2013，40(6):708-714．

[2] 霍春亮，葉小明，高振南，等．儲(chǔ)層內(nèi)部小尺度構(gòu)型單元界面等效表征方法[J]．中國(guó)海上油氣，2016，28(1)：54-59.DOI：10.11935/j.issn.1673-1506.2016.01.008.

HUO Chunliang，YE Xiaoming，GAO Zhennan，et al.Equivalent characterization method of small scale reservoir configuration unit interface[J]．China Offshore Oil and Gas，2016，28(1)：54-59.DOI：10.11935/j.issn.1673-1506.2016.01.008.

[3] 趙永軍,舒曉,胡勇，等．一種復(fù)雜曲流帶儲(chǔ)層三維構(gòu)型建模新方法[J]．中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版)，2015，39(1)：1-7．

ZHAO Yongjun,SHU Xiao,HU Yong,et al.A new 3D architecture modeling method of complex meander belt reservoir[J].Journal of China University of Petroleum(Edition of Natural Science)，2015，39(1)：1-7．

[4] 劉超，趙春明，廖新武，等．海上油田大井距條件下曲流河儲(chǔ)層內(nèi)部構(gòu)型精細(xì)解剖及應(yīng)用分析[J]．中國(guó)海上油氣，2014，26(1)：58-64．

LIU Chao，ZHAO Chunming，LIAO Xinwu，et al．A refined anatomy of the internal structure of meandering river reservoirs under large well spacing in offshore oilfields and its application[J]．China Offshore Oil and Gas，2014，26(1)：58-64．

[5] 周銀邦,吳勝和,計(jì)秉玉,等.曲流河儲(chǔ)層構(gòu)型表征研究進(jìn)展[J].地球科學(xué)進(jìn)展,2011,26(7):695-702.

ZHOU Yinbang,WU Shenghe,JI Bingyu,et al.Research program on the characterization of fluvial reservoir architecture[J].Advances in Earth Science,2011,26(7):695-702.

[6] 閆百泉，馬世忠，王龍，等．曲流點(diǎn)壩內(nèi)部剩余油形成與分布規(guī)律物理模擬[J]．地學(xué)前緣，2008，15(1)：65-70．

YAN Baiquan，MA Shizhong，WANG Long，et al．The formation and distribution of residual oil in meander point bar by physical modeling[J]．Earth Science Frontiers，2008，15(1)：65-70．

[7] 岳大力，吳勝和，劉建民．曲流河點(diǎn)壩地下儲(chǔ)層構(gòu)型精細(xì)解剖方法[J]．石油學(xué)報(bào)，2007，28(4)：99-103．

YUE Dali，WU Shenghe，LIU Jianmin．An accurate method for anatomizing architecture of subsurface reservoir in point bar of meandering river[J]．Acta Petrolei Sinica，2007，28(4)：99-103．

[8] 白振強(qiáng)，王清華，杜慶龍，等．曲流河砂體三維構(gòu)型地質(zhì)建模及數(shù)值模擬研究[J]．石油學(xué)報(bào)，2009，30(6)：898-902．

BAI Zhenqiang，WANG Qinghua，DU Qinglong，et al．Study on 3D architecture geology modeling and digital simulation in meandering reservoir[J]．Acta Petrolei Sinica，2009，30(6)：898-902．

[9] 吳勝和，岳大力，劉建民，等．地下古河道儲(chǔ)層構(gòu)型的層次建模研究[J]．中國(guó)科學(xué)D 輯：地球科學(xué)，2008，38(增刊1 )：111-121．

WU Shenghe， YUE Dali， LIU Jianmin，et al.Hierarchy modeling of subsurface palaeochannel reservoir architecture[J].Science in China Series D：Earth Sciences，2008，51(S2)：126-137．

[10] 尹艷樹(shù)，張昌民，李少華，等．一種新的曲流河點(diǎn)壩側(cè)積層建模方法[J]．石油學(xué)報(bào)，2011，32(2)：315-319．

YIN Yanshu，ZHANG Changmin，LI Shaohua，et al．A new stochastic modeling method for 3-D forecasting lateral accretion bedding of point bars in meandering rivers[J]．Acta Petrolei Sinica，2011，32(2)：315-319．

[11] 高博禹，孫立春，胡光義，等．基于單砂體的河流相儲(chǔ)層地質(zhì)建模方法探討[J]．中國(guó)海上油氣，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．

[12] 范崢，吳勝和，岳大力,等．曲流河點(diǎn)壩內(nèi)部構(gòu)型的嵌入式建模方法研究[J]．中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版)，2012，36(3)：1-6．

FAN Zheng ,WU Shenghe，YUE Dali，et al．Embedding modeling method for internal architecture of point bar sand body in meandering river reservoir[J].Journal of China University of Petroleum(Edition of Natural Science)，2012，36(3)：1-6．

[13] 李宇鵬,吳勝和,耿麗慧,等.基于空間矢量的點(diǎn)壩砂體儲(chǔ)層構(gòu)型建模[J].石油學(xué)報(bào),2013,34(1):133-139.

LI Yupeng,WU Shenghe,GENG Lihui,et al.Spatial-vector-based reservoir architecture modeling of point-bar sand[J].Acta Petrolei Sinica,2013,34(1):133-139.

[14] 霍春亮,劉松,古莉,等.一種定量評(píng)價(jià)儲(chǔ)集層地質(zhì)模型不確定性的方法[J]．石油勘探與開(kāi)發(fā)，2007， 34(5):574-579．

HUO Chunliang，LIU Song，GU Li，et al．A quantitative method for appraising reservoir geological model uncertainty[J].Petroleum Exploration and Development，2007，34(5):574-579．

[15] DUAN Taizhong.Similarity measure of sedimentary successions and its application in inverse stratigraphic modeling[J].Petroleum Science,2017，14(3):484-492．

[16] 徐長(zhǎng)貴，杜曉峰．陸相斷陷盆地源匯理論工業(yè)化應(yīng)用初探:以渤海海域?yàn)槔齕J]．中國(guó)海上油氣，2017，29(4)：9-18．DOI：10.11935/j.issn.1673-1506.2017.04.002.

XU Changgui，DU Xiaofeng.Industrial application of source-to-sink theory in continental rift basin:a case study of Bohai sea area[J]．China Offshore Oil and Gas，2017，29(4)：9-18．DOI：10.11935/j.issn.1673-1506.2017.04.002.

[17] 魏洪濤.遼中凹陷北部東二下亞段湖底扇沉積數(shù)值模擬及應(yīng)用[J].巖性油氣藏,2015,27(5):183-188.

WEI Hongtao.Numerical simulation of sublacustrine fan deposition of lower Ed2Formation and its application in northern Liaozhong Depression[J].Lithologic Reservoirs,2015,27(5):183-188.