付 超，謝玉洪，宋來明，王 暉，徐 偉，苑志旺，陳國寧

1 中海油研究總院有限責(zé)任公司；2 中國海洋石油集團(tuán)有限公司

0 前 言

隨著油氣勘探向深水區(qū)進(jìn)軍，深水區(qū)沉積儲(chǔ)層的刻畫與表征成為了目前研究的熱點(diǎn)。濁積砂體為深水區(qū)的主要儲(chǔ)層類型，占目前深水區(qū)已發(fā)現(xiàn)油氣田探明儲(chǔ)量的75%以上［1-2］。由于較高的勘探成本，深水區(qū)勘探前期的鉆井?dāng)?shù)量少，因此儲(chǔ)層研究通常需要基于稀井網(wǎng)展開［3-4］；同時(shí)，深水區(qū)的地震工區(qū)范圍有限，無資料覆蓋區(qū)的范圍較大［5］。在陸相湖盆的早期研究時(shí)遇到相似的問題，前人提出了相控逐級(jí)建模的“兩步法”建模思路，即首先建立沉積相、儲(chǔ)層結(jié)構(gòu)模型，然后根據(jù)不同沉積相的儲(chǔ)層參數(shù)定量分布規(guī)律，按相帶的不同進(jìn)行井間插值或隨機(jī)模擬，以達(dá)到精細(xì)描述儲(chǔ)層空間分布格局的目的［6-9］。付超等［10］和張文彪等［11］將該方法應(yīng)用于深水峽谷建模過程中，實(shí)現(xiàn)了融合地震資料建模。趙迎月等［12］提出融合地震正演模擬的建模思路，實(shí)現(xiàn)了深水區(qū)的無井條件下的儲(chǔ)層建模。近些年隨著多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)方法的成熟［13］，在利用地震驅(qū)動(dòng)建模的過程中，分別引入Snesim 算法和Alluvsim 算法對(duì)高彎曲度峽谷進(jìn)行三維圖像訓(xùn)練，結(jié)合基于目標(biāo)建模的方法，提高了深水淺層峽谷體系的模型精度［14-15］。然而，瓊東南盆地中央峽谷陵水段埋深較大（超3 000 m），峽谷順直且古地形變化較快，造成濁積水道的單期次砂體難以識(shí)別［16-17］，對(duì)這類深水濁積砂體亟需引入新的建模方法。

深水濁積砂體建模存在2 個(gè)方面特有的難題：①濁積砂體通常表現(xiàn)為幕式發(fā)育的特點(diǎn)，在垂向上表現(xiàn)為厚層砂巖和厚層泥巖互層的沉積規(guī)律，由此造成不同分層內(nèi)數(shù)據(jù)點(diǎn)密度的分布范圍差異較大。②濁流在形成過程中，流動(dòng)性受限于峽谷壁的發(fā)育程度和穩(wěn)定程度，濁積砂體在平面上表現(xiàn)為短變程-高塊金值的特點(diǎn)［18-19］。卜范青等［20-21］提出了利用分區(qū)建立變差函數(shù)的方法，克服該類型沉積體在平面上相變快的難題。付超等［10］和張佳佳等［22］提出利用分層相控約束、多重屬性協(xié)同的方法構(gòu)建深水沉積模型，展現(xiàn)了模型的垂向演化。Ebong 等［23］綜合指出分區(qū)分層位使用多種建模方法可以較大地提高模型的可靠程度。

本文在前人提出的“兩步法”分層分區(qū)建模的基礎(chǔ)上，針對(duì)深水濁積砂體搬運(yùn)過程的特點(diǎn)，引入“沉積物輸送體系理論”，這不僅可以在宏觀尺度計(jì)算沉積物通量和輸送效率，預(yù)測濁積砂體的儲(chǔ)層結(jié)構(gòu)，而且可以在微觀尺度模擬顆粒搬運(yùn)過程和輸砂機(jī)制，厘清微觀儲(chǔ)層結(jié)構(gòu)特征、遞變規(guī)律及影響因素。利用該理論進(jìn)行地質(zhì)建模的方法如下：首先，使用沉積正演模擬早中新世沉積充填過程；而后，沿峽谷建立虛擬井，通過神經(jīng)網(wǎng)絡(luò)方法分析其與地震數(shù)據(jù)、測井?dāng)?shù)據(jù)等資料的相關(guān)性，擬合得到虛擬井的沉積相與物性參數(shù)等數(shù)據(jù)；最后，按照分區(qū)分層的建模思路建立復(fù)合儲(chǔ)層模型。該方法解決了兩點(diǎn)變差函數(shù)建模和基于目標(biāo)建模對(duì)密井網(wǎng)的需求問題，同時(shí)也彌補(bǔ)了多點(diǎn)建模中訓(xùn)練圖像過于主觀的缺陷，為深水稀井網(wǎng)密度下儲(chǔ)層建模提供了解決方案。

1 中央峽谷沉積背景與地層格架

瓊東南盆地處于海南島以南、西沙群島以北的海域（圖1a），是一個(gè)大型的新生代裂谷型大陸邊緣盆地，呈北東向展布。研究區(qū)所在的中央峽谷主要發(fā)育在盆地中部，根據(jù)其所在的凹陷可以劃分為樂東段、陵水段、松濤寶島段和長昌段。中央峽谷發(fā)育于中新世，主要經(jīng)歷了早中新世局部抬升、晚中新世沉降2 個(gè)大的構(gòu)造演化階段，整體構(gòu)造活動(dòng)相對(duì)穩(wěn)定，工區(qū)無大型斷裂發(fā)育［24-27］。峽谷內(nèi)部濁積砂體的快速侵蝕，使峽谷具有相對(duì)穩(wěn)定的峽谷壁和全盆地可對(duì)比追蹤的沉積邊界。峽谷內(nèi)部沉積地層包括新近系中新統(tǒng)三亞組、梅山組、黃流組，上新統(tǒng)鶯歌海組，第四系更新統(tǒng)樂東組［27］，其中黃流組時(shí)期以峽谷充填沉積為主（圖1b）。

圖1 瓊東南盆地中央峽谷發(fā)育位置與充填期次劃分Fig.1 Location of Central Canyon in Qiongdongnan Basin and its filling stage division

瓊東南盆地中央峽谷的充填演化可劃分為3個(gè)演化階段：初始形成階段、穩(wěn)定發(fā)育階段和后期改造階段。根據(jù)巖性變化將本次研究的目的層段黃流組由上至下劃分為6 個(gè)氣組（圖1b），即I、II1、II2、III、IV 和V 氣組，其中IV—V 氣組對(duì)應(yīng)峽谷初始形成階段，II2—III 氣組對(duì)應(yīng)峽谷穩(wěn)定發(fā)育階段，I—II1氣組對(duì)應(yīng)峽谷后期改造階段，三者之間不存在明顯的沉積間斷。

2 深水峽谷樣式與沉積參數(shù)

2.1 深水峽谷刻畫與定量參數(shù)統(tǒng)計(jì)

選擇中央峽谷陵水段由上游至下游的6條地震剖面（圖2），依次刻畫峽谷形態(tài)。

圖2 瓊東南盆地不同階段深水峽谷及其砂巖厚度平面圖Fig.2 Seismic profiles showing deep water canyon at different stages and its sand thickness contours in Qiongdongnan Basin

峽谷上游（圖2a，2b） IV—V 氣組地層相對(duì)較薄，無法識(shí)別出明顯的峽谷邊界。非限制性峽谷內(nèi)部可以識(shí)別少量切疊的濁積水道帶，水道帶寬度較大，平面上同樣無明確邊界，水道帶厚度為20～30 m。I—III 氣組厚度由峽谷內(nèi)部向兩翼快速減?。? 個(gè)氣組厚度分別為60 m、40 m 和20 m）。限制性峽谷內(nèi)部可以識(shí)別出3 條明顯的水道。其中，III 氣組水道呈對(duì)稱寬“U”形，水道深度較大（約為60 m）；II氣組水道變?yōu)椤癢”形和不對(duì)稱“U”形，水道深度減?。s為10 m）。

峽谷中游（圖2c,2d） IV—V 氣組表現(xiàn)出與上游峽谷和水道相近的特點(diǎn)，但是水道帶寬度增大。I—III 氣組峽谷充填厚度較上游有所增加，同樣可以識(shí)別3～4 條水道，表現(xiàn)出對(duì)稱的寬“U”形，水道厚度減?。s為10 m）。參考龔承林等［17］對(duì)研究區(qū)沉積展布的研究，認(rèn)為該時(shí)期峽谷樣式較為順直，濁積砂體多在峽谷底部充填。

峽谷下游（圖2e,2f） IV—V 氣組地層厚度增大，其中除水道帶外，可見少量條帶狀透鏡體，厚度范圍為5～10 m。參考程聰?shù)龋?8］對(duì)研究區(qū)的沉積認(rèn)識(shí)，將透鏡體解釋為MTDs。I—III 氣組峽谷充填厚度較上游有明顯減薄，可以識(shí)別2～3條水道，水道內(nèi)部連續(xù)性增強(qiáng)，表現(xiàn)出寬“U”形和碟形。

通過單井砂體厚度標(biāo)定，結(jié)合古地貌趨勢繪制研究區(qū)砂巖厚度平面圖（圖2g—2i）。砂巖厚度反映在IV—V 氣組沉積時(shí)期，研究區(qū)表現(xiàn)出北高南低的特點(diǎn)，峽谷未發(fā)育，僅存在少量沖溝。在II2—III 氣組沉積時(shí)期，濁積砂體快速充填進(jìn)入半限制性峽谷內(nèi)部，砂厚高值區(qū)顯示沉積物在峽谷上游和中上游堆積；該時(shí)期峽谷進(jìn)入半限制性背景，峽谷面積約為450 km2。在I—II1氣組沉積時(shí)期，隨著濁積砂體的持續(xù)充填，沉積物向下游搬運(yùn)；該時(shí)期中游區(qū)域砂體厚度較大，同時(shí)隨著峽谷的進(jìn)一步充填，峽谷范圍減小，其面積約為360 km2。

基于上述對(duì)于地震剖面的解釋，對(duì)研究區(qū)深水峽谷沉積體參數(shù)（圖3）進(jìn)行統(tǒng)計(jì)，統(tǒng)計(jì)參數(shù)包括形態(tài)參數(shù)（峽谷方向、相對(duì)彎曲度）和剖面參數(shù)（水道的擺動(dòng)幅度、擺動(dòng)距離、寬度和深度等）［17～18］（表1）。其中形態(tài)參數(shù)基于砂巖厚度圖中砂巖厚度高值區(qū)走向測量得到；剖面參數(shù)基于連續(xù)地震剖面（時(shí)間域）統(tǒng)計(jì)得到，結(jié)合時(shí)深轉(zhuǎn)換公式，得到深度域參數(shù)。表1 中的Zone 為本次建模研究的垂向單元，在地層劃分中對(duì)應(yīng)氣層組，由下至上劃分為Zone1（對(duì)應(yīng)IV—V 氣組）、Zone2（對(duì)應(yīng)II2—III 氣組）、Zone3（對(duì)應(yīng)I—II1氣組）三個(gè)單元。

表1 瓊東南盆地深水峽谷沉積體參數(shù)統(tǒng)計(jì)表Table 1 Parameters of the sedimentary bodies of deep water canyon in Qiongdongnan Basin

圖3 深水峽谷剖面參數(shù)示意圖Fig.3 A sketch illustrating deep water canyon profile parameters

2.2 深水峽谷充填過程模擬與虛擬井建立

根據(jù)上述建立的中央峽谷陵水段充填規(guī)律和沉積體參數(shù)統(tǒng)計(jì)，利用DionisosFlow?軟件進(jìn)行沉積正演模擬（圖4）。參考胡勇等［29］和劉建良等［30］的沉積正演模擬方法，步驟可分為沉積參數(shù)輸入、古地貌輸入、沉降速率輸入、模擬步長設(shè)定和模擬結(jié)果分析。其中：①沉積參數(shù)輸入（包括砂泥占比等）可依據(jù)目前鉆井資料統(tǒng)計(jì)得到。②古地貌數(shù)據(jù)可基于三維地震資料解釋得到（圖4a），其具體做法為通過井震標(biāo)定和等時(shí)界面追蹤，對(duì)于泥質(zhì)含量較高的區(qū)域需進(jìn)行去壓實(shí)校正。③沉降速率輸入?yún)⒖糑onsoer等［31］建立的砂體幾何形態(tài)和沉積通量關(guān)系，通過地震資料解釋（與反演）數(shù)據(jù)識(shí)別并測量砂體幾何參數(shù)，計(jì)算出研究區(qū)沉積物供給速率。此外，基于沉積通道附近的巖（壁）心粒度和沉積構(gòu)造（交錯(cuò)層理、遞變層理等）規(guī)模，改進(jìn)了Komar［32］利用沉積構(gòu)造和粒度遞變計(jì)算濁流流速的方法，恢復(fù)中央峽谷濁流流動(dòng)過程。計(jì)算結(jié)果與袁玉松等［33］的計(jì)算結(jié)果（179 m/Ma）相近。④模擬步長設(shè)置為為0.5 Ma，滑塌過程設(shè)定的臨界坡度為10°，模擬得到黃流組沉積時(shí)期的充填過程（圖4b）。然后，沿沉積物流動(dòng)方向設(shè)定虛擬井，虛擬井井距為1 000 m，在峽谷內(nèi)均勻分布。

圖4 瓊東南盆地深水峽谷基于沉積正演模擬建立的虛擬井與虛擬井?dāng)?shù)據(jù)校對(duì)Fig.4 Pseudo-wells and pseudo-well data verification based on sedimentary forward modeling of the deep water canyon in Qiongdongnan Basin

分三步驗(yàn)證沉積模擬結(jié)果的可靠性（圖4c）：①擬合井點(diǎn)地震（RMS）屬性與井點(diǎn)沉積模擬結(jié)果（砂地比）（擬合度高于70%）；②疊合RMS 屬性高值區(qū)與沉積模擬結(jié)果富砂區(qū)域（疊合占比大于70%）；③通過單井對(duì)比沉積模擬結(jié)果驗(yàn)證關(guān)鍵界面發(fā)育位置（相對(duì)可靠）。

最終選取可靠度較高的模擬結(jié)果用于虛擬井?dāng)?shù)據(jù)學(xué)習(xí)中。按照Banchs 等［34］提出的神經(jīng)網(wǎng)絡(luò)方法：首先對(duì)虛擬井和實(shí)際井的參數(shù)進(jìn)行相關(guān)性分析，選擇其中相關(guān)性較低的參數(shù)（均方根振幅、波阻抗和密度）進(jìn)行擬合，確定儲(chǔ)層巖石類型（圖5a—5c），明確濁積砂體的泥質(zhì)含量范圍為0.2～0.5（平均值0.3）（圖5d）、有效孔隙度范圍為0.2～0.8（平均值0.25）（圖5e）、含氣飽和度范圍0.4～0.9（平均值0.85）（圖5f）。最終根據(jù)測井曲線與地震屬性值交會(huì)圖版，確定對(duì)巖性與物性較為敏感的曲線類型，并擬合得到虛擬井的物性分布。

圖5 瓊東南盆地深水峽谷沉積體測井參數(shù)-測井參數(shù)與測井參數(shù)-地震屬性交會(huì)圖Fig.5 Crossplot charts of logging parameter-logging parameter and logging parameter-seismic attribute of deep water canyon sedimentary body in Qiongdongnan Basin

3 深水峽谷分區(qū)分層建模策略

上述對(duì)峽谷充填過程的恢復(fù)分析表明，深水峽谷表現(xiàn)出從早期非限制性峽谷到后期半限制性峽谷的演化特征。因此，需要考慮差異沉積背景下不同類型峽谷的建模過程。對(duì)于限制性峽谷，目前通常采用改進(jìn)算法的多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)建模方法［11,14］，其基本過程為：①根據(jù)峽谷流動(dòng)方向確定主要發(fā)育的水道中線；②根據(jù)幾何參數(shù)給出水道的彎曲程度、下切程度等參數(shù)進(jìn)行圖像訓(xùn)練；③在輔助參數(shù)（例如砂地比）約束的背景下建立對(duì)應(yīng)的模型。而對(duì)于非限制性峽谷，主要采用相帶邊界約束下的序貫高斯（如邊界清晰可采用截?cái)喔咚梗┙７椒ǎ?0,35］。

本文采用針對(duì)不同層位、不同區(qū)域進(jìn)行分區(qū)分層建模的思路（圖6）。三維構(gòu)造模型平面上采用50 m×50 m 的網(wǎng)格；垂向上，根據(jù)界面接觸關(guān)系和氣組分布特點(diǎn)，劃分了3 個(gè)Zone：Zone1 包含IV—V 氣組，對(duì)應(yīng)峽谷初始形成期，厚度范圍為30～40 m，劃分為40 個(gè)Layers；Zone2 包含II2—III氣組，對(duì)應(yīng)峽谷最大充填期，厚度范圍為20～100 m，劃分為80 個(gè)Layers；Zone3 包含I—II1氣組，對(duì)應(yīng)峽谷后期改造期，厚度范圍為0～40 m，劃分為40 個(gè)Layers。由于研究區(qū)內(nèi)峽谷存在快速尖滅，因此采用角點(diǎn)網(wǎng)格進(jìn)行建模。在此基礎(chǔ)上建立沉積相模型、巖性模型和物性（孔隙度、滲透率和含氣飽和度）模型。對(duì)于不同分層根據(jù)其沉積充填特征選擇對(duì)應(yīng)的建模方法。

圖6 瓊東南盆地深水峽谷沉積體建模方法流程圖Fig.6 Flow chart of modeling method of deep water canyon sedimentary body in Qiongdongnan Basin

Zone1（IV—V 氣組）峽谷初始形成階段 該時(shí)期峽谷陸坡走向平行發(fā)育，表現(xiàn)為非限制性沉積特點(diǎn)，砂體多表現(xiàn)出分區(qū)分布的特點(diǎn)［16］，主要集中在峽谷中游的W4井附近（圖2g）。該時(shí)期的沉積體無法識(shí)別出其準(zhǔn)確的邊界，因此主要利用井震數(shù)據(jù)擬合虛擬井?dāng)?shù)據(jù)，進(jìn)行序貫高斯模擬。

Zone2（II2—III 氣組）峽谷穩(wěn)定發(fā)育階段 該時(shí)期峽谷內(nèi)部快速充填，多表現(xiàn)出厚層的濁積砂體（圖2），砂體較為連片分布，此外根據(jù)地震剖面解釋，該時(shí)期復(fù)合水道帶寬度接近峽谷寬度，故為了簡化計(jì)算，采用截?cái)喔咚鼓M進(jìn)行儲(chǔ)層建模。

Zone3（I—II1氣組）峽谷后期改造階段 上游地震剖面（圖2a—2b）顯示具有明確的孤立濁積水道砂體，并且水道發(fā)育較為連續(xù)，因此可采用基于訓(xùn)練圖像的多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)建模方法；而中下游（圖2c—2e）則表現(xiàn)出快速充填的特征，故采用基于趨勢面的截?cái)喔咚鼓M方法進(jìn)行建模。

4 相控儲(chǔ)層建模

4.1 沉積相模型隨機(jī)模擬

根據(jù)前述沉積相帶認(rèn)識(shí)，劃分出3 種類型的相帶（濁積水道、天然堤與峽谷壁、遠(yuǎn)洋細(xì)粒泥巖）按照上述建模策略進(jìn)行逐層建模（圖7）。

圖7 瓊東南盆地深水峽谷沉積體RMS屬性模型與沉積相隨機(jī)模型Fig.7 RMS attribute model and sedimentary facies stochastic model of deep water canyon sedimentary body in Qiongdongnan Basin

Zone1（IV—V 氣組） 單井（包含實(shí)際井與虛擬井）粗化數(shù)據(jù)顯示濁積水道占比為24.8%，天然堤與峽谷占比為15.0%，遠(yuǎn)洋細(xì)粒泥巖占比為60.2%。分別對(duì)3 個(gè)相帶的變差函數(shù)進(jìn)行分析，利用球形曲線，設(shè)定主變程范圍為165 m，次級(jí)變程范圍為25 m（主變程范圍參考水道長度，次級(jí)變程范圍參考峽谷寬度［36-38］）。在沉積相建模中選取序貫指示模擬，利用Cell visitation 引入均方根振幅屬性體（圖7c）進(jìn)行聯(lián)合建模，得到IV—V 氣組的沉積相模型。統(tǒng)計(jì)模型中濁積水道、天然堤與峽谷、遠(yuǎn)洋細(xì)粒泥巖的比例為2:1:5。

Zone2（II2—III 氣組） 選取與Zone1 相同的方法進(jìn)行數(shù)據(jù)分析，同樣利用Cell visitation 引入均方根振幅屬性體進(jìn)行截?cái)喔咚菇?，截?cái)噙吔绺鶕?jù)圖2 中砂巖邊界（砂厚約為20 m）進(jìn)行確定，其中截?cái)噙吔鐑?nèi)部濁積水道、天然堤與峽谷、遠(yuǎn)洋細(xì)粒泥巖的比例為5:2:3，截?cái)噙吔缤馊叩谋壤秊?:4:5。

Zone3（I—II1氣組） 劃分為上游部分和中下游部分進(jìn)行分區(qū)建模。上游部分基于Petrel 軟件中Snesim 算法，模擬過程主要包括數(shù)據(jù)準(zhǔn)備、掃描訓(xùn)練圖像并構(gòu)建穩(wěn)定的搜索樹、調(diào)整匹配參數(shù)、選擇隨機(jī)路徑，利用基于條件概率分布函數(shù)的多點(diǎn)模擬方法，最終建立相模型。其結(jié)果驗(yàn)證可采用對(duì)比地震數(shù)據(jù)與模型剖面中的相帶分布，進(jìn)而確定其模型的合理性。對(duì)于中下游部分，采用截?cái)喔咚惯M(jìn)行建模。

4.2 物性模型隨機(jī)模擬

在沉積相模型的基礎(chǔ)上，對(duì)研究區(qū)所有井的泥質(zhì)含量、孔隙度、滲透率和含氣飽和度4個(gè)連續(xù)數(shù)據(jù)進(jìn)行粗化，在分區(qū)分層分析變差函數(shù)（表2）的基礎(chǔ)上建立對(duì)應(yīng)的模型（圖8）。

表2 瓊東南盆地深水峽谷沉積體分區(qū)分層位建模參數(shù)統(tǒng)計(jì)Table 2 Zoned and layered modeling parameters of deep water canyon sedimentary body in Qiongdongnan Basin

圖8 瓊東南盆地深水峽谷沉積體泥質(zhì)含量、孔隙度、滲透率、含氣飽和度模型建立Fig.8 Establishment of mud content,porosity,permeability and gas saturation model of deep water canyon sedimentary body in Qiongdongnan Basin

4.2.1 泥質(zhì)含量模型

針對(duì)Zone1（IV—V 氣組）時(shí)期峽谷處于非限制性背景這一特征，在變差函數(shù)分析方面，與上述沉積相模型一樣考慮長變程高塊金值的特點(diǎn)進(jìn)行建模［38］。結(jié)合RMS屬性進(jìn)行隨機(jī)建模，生成泥質(zhì)含量模型的100 個(gè)隨機(jī)實(shí)現(xiàn)，最終通過算術(shù)平均處理得到泥質(zhì)含量模型（圖8c）。針對(duì)Zone2（II2—III 氣組）水道快速充填且較為順直的特點(diǎn)［17］，借鑒Rongier等［40］對(duì)于疊置水道的建模策略，使用上述相模型作為約束邊界，采用球狀模型確立主副變程范圍（基臺(tái)值為5，塊金值為2.5，主變程為2 000 m，次變程為250 m），得到該層對(duì)應(yīng)的泥質(zhì)含量模型（圖8b）。針對(duì)Zone3（I—II1氣組）同樣采用分區(qū)建模（圖8a），上游區(qū)在沉積相模型的約束下，采用球狀模型確定變差函數(shù)范圍（基臺(tái)值為2.5，塊金值為1），變程范圍內(nèi)數(shù)據(jù)量吻合程度較好，中下游濁積砂體沉積相變較快，其主變程為1 800 m，基臺(tái)值為2.5，塊金值為1。

建立的泥質(zhì)含量模型顯示，Zone1（IV—V 氣組）峽谷充填早期高砂地比范圍主要集中在峽谷上游（圖8c），表現(xiàn)出多個(gè)孤立的區(qū)域，指示了濁積砂體和MTDs 在峽谷中的滯留沉積。該認(rèn)識(shí)與李超等［16］，通過地震剖面解釋的結(jié)論一致。Zone2（II2—III 氣組）的砂地比高值區(qū)集中于峽谷內(nèi)部和峽谷壁一側(cè)（圖8b），結(jié)合上述地震剖面解釋結(jié)論［16］認(rèn)為這是由于該時(shí)期較為發(fā)育的濁積砂體和MTDs 快速充填峽谷所致。Zone3（I—II1氣組）砂體比高值區(qū)集中在峽谷內(nèi)部（圖8a），并且上游區(qū)域砂地比高值區(qū)范圍小于中下游區(qū)域，認(rèn)為其沉積成因?yàn)闈岱e砂體經(jīng)搬運(yùn)后在中下游局部平坦區(qū)域堆積。

4.2.2 孔隙度、滲透率與含氣飽和度模型

孔隙度、滲透率與含氣飽和度的建模方法與泥質(zhì)含量模型的建立方法基本一致，并且分布相似。不同的是：泥質(zhì)含量的原始數(shù)據(jù)是單峰式分布，轉(zhuǎn)化成正態(tài)數(shù)據(jù)橫坐標(biāo)值變化較?。豢紫抖?、滲透率和含氣飽和度數(shù)據(jù)為多峰式，因此轉(zhuǎn)換成正態(tài)數(shù)據(jù)橫坐標(biāo)值變化較大。同樣采用沉積相進(jìn)行約束，得到對(duì)應(yīng)的模型（圖8）。

4.3 模型結(jié)果質(zhì)控

本次研究從2 個(gè)方面對(duì)模型建立質(zhì)控：①通過物性對(duì)比圖進(jìn)行質(zhì)檢（圖9a）。對(duì)比測井解釋的孔隙度、滲透率與含水飽和度的取值范圍和模型結(jié)果，偏差小于10%，認(rèn)為模型較為可靠。②通過物性分布圖進(jìn)行質(zhì)檢（圖9b）。對(duì)比模型參數(shù)、粗化參數(shù)與實(shí)際井?dāng)?shù)據(jù)之間趨勢分布的相關(guān)性，粗化結(jié)果與單井實(shí)測結(jié)果差異小于10%。

圖9 基于物性對(duì)比和物性分布的模型結(jié)果質(zhì)控Fig.9 Quality control of model results based on the property comparison and property distribution

該方法的模型結(jié)果與井點(diǎn)存在差異的主要原因是上游區(qū)域沉積地層相對(duì)下游較薄，但是各層采用統(tǒng)一的網(wǎng)格密度，造成上游網(wǎng)格數(shù)據(jù)粗化結(jié)果與實(shí)際井存在誤差。模型參數(shù)與粗化參數(shù)之間的差異是由于虛擬井主要沿峽谷設(shè)置，對(duì)于峽谷壁兩側(cè)的物性高值區(qū)難以進(jìn)行控制所致?？傊c傳統(tǒng)建模方法得到的峽谷模型相比，該方法采用融合沉積正演模型設(shè)定虛擬井，并利用分區(qū)分層建模的策略，最大程度地保證了對(duì)峽谷內(nèi)物性展布規(guī)律的控制，后續(xù)建模方法可以在物性敏感性分析、多點(diǎn)建模策略選擇等方面進(jìn)行改進(jìn)［41-43］。

5 結(jié) 論

（1）以瓊東南盆地中央峽谷陵水段為研究對(duì)象，引入基于峽谷充填參數(shù)的沉積正演模擬，模擬結(jié)果顯示峽谷初始形成時(shí)期濁積物在上游非限制性峽谷內(nèi)發(fā)育，后期在限制性峽谷內(nèi)部充填。在此基礎(chǔ)上，利用正演模擬結(jié)合地震屬性和測井?dāng)?shù)據(jù)建立虛擬井，通過交會(huì)圖分析地震屬性與測井?dāng)?shù)據(jù)之間的關(guān)系，建立虛擬井物性分布規(guī)律，明確密度、波阻抗、自然伽馬在不同巖性間的截?cái)喾秶?/p>

（2）提出分層分區(qū)“兩步法”的建模策略，明確了對(duì)于非限制性峽谷采用序貫高斯模擬方法，早期快速充填順直水道采用截?cái)喔咚鼓M方法，晚期擺動(dòng)水道上游采用多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)方法、中下游采用截?cái)喔咚鼓M方法；實(shí)現(xiàn)了“兩步法”在瓊東南盆地中央峽谷的應(yīng)用，刻畫了峽谷內(nèi)部物性差異，并分析其沉積成因。