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

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

0 前 言

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

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

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

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

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

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

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

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

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

選擇中央峽谷陵水段由上游至下游的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 氣組地層相對較薄，無法識別出明顯的峽谷邊界。非限制性峽谷內(nèi)部可以識別少量切疊的濁積水道帶，水道帶寬度較大，平面上同樣無明確邊界，水道帶厚度為20～30 m。I—III 氣組厚度由峽谷內(nèi)部向兩翼快速減?。? 個氣組厚度分別為60 m、40 m 和20 m）。限制性峽谷內(nèi)部可以識別出3 條明顯的水道。其中，III 氣組水道呈對稱寬“U”形，水道深度較大（約為60 m）；II氣組水道變?yōu)椤癢”形和不對稱“U”形，水道深度減?。s為10 m）。

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

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

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

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

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

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

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

最終選取可靠度較高的模擬結(jié)果用于虛擬井?dāng)?shù)據(jù)學(xué)習(xí)中。按照Banchs 等［34］提出的神經(jīng)網(wǎng)絡(luò)方法：首先對虛擬井和實際井的參數(shù)進行相關(guān)性分析，選擇其中相關(guān)性較低的參數(shù)（均方根振幅、波阻抗和密度）進行擬合，確定儲層巖石類型（圖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ù)測井曲線與地震屬性值交會圖版，確定對巖性與物性較為敏感的曲線類型，并擬合得到虛擬井的物性分布。

圖5 瓊東南盆地深水峽谷沉積體測井參數(shù)-測井參數(shù)與測井參數(shù)-地震屬性交會圖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ū)分層建模策略

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

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

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

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

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

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

4 相控儲層建模

4.1 沉積相模型隨機模擬

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

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

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

Zone2（II2—III 氣組） 選取與Zone1 相同的方法進行數(shù)據(jù)分析，同樣利用Cell visitation 引入均方根振幅屬性體進行截斷高斯建模，截斷邊界根據(jù)圖2 中砂巖邊界（砂厚約為20 m）進行確定，其中截斷邊界內(nèi)部濁積水道、天然堤與峽谷、遠洋細粒泥巖的比例為5:2:3，截斷邊界外三者的比例為1:4:5。

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

4.2 物性模型隨機模擬

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

表2 瓊東南盆地深水峽谷沉積體分區(qū)分層位建模參數(shù)統(tǒng)計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ì)含量模型

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

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

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

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

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

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

5 結(jié) 論

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

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