基于地震等時格架的傾角導向儲層靜態(tài)建模方法

樂 靖，王 暉，范廷恩，趙衛(wèi)平，樊鵬軍

（中海油研究總院，北京100028）

基于地震等時格架的傾角導向儲層靜態(tài)建模方法

樂 靖，王 暉，范廷恩，趙衛(wèi)平，樊鵬軍

（中海油研究總院，北京100028）

為了建立等時性更好的地層格架，提高井網(wǎng)稀疏時儲層參數(shù)的預測精度，提出了基于地震等時格架的傾角導向儲層靜態(tài)建模方法。以地震傾角體為導向劃分垂向網(wǎng)格層并建立三維地層網(wǎng)格模型，利用敏感地震屬性計算儲層參數(shù)的空間變化規(guī)律，再結(jié)合井筒信息預測儲層參數(shù)值，建立三維儲層參數(shù)模型。與傳統(tǒng)構(gòu)造建模方法相比，該方法的優(yōu)點是采用了傾角導向來建立等時地層格架模型。該方法基于地震屬性求取的轉(zhuǎn)換函數(shù)來建立儲層參數(shù)模型，有別于傳統(tǒng)地質(zhì)統(tǒng)計學采取變差函數(shù)建模的方法，克服了井網(wǎng)稀疏條件下建模精度低的缺陷。模型測試與海上油田的實際應用表明，與傳統(tǒng)地質(zhì)統(tǒng)計學方法相比，該方法建立了更等時的三維地層網(wǎng)格模型，提高了井網(wǎng)稀疏時儲層參數(shù)建模的精度。

傾角體；傾角導向；地震驅(qū)動；等時格架；儲層靜態(tài)建模

儲層靜態(tài)建模實現(xiàn)了對地下地質(zhì)體的三維定量表征，是油藏評價、動態(tài)模擬和剩余油預測的技術(shù)基礎(chǔ)，其與地下實際情況的吻合程度直接影響油田的開發(fā)效果［1］。因此，如何使地質(zhì)模型更貼近地下的真實情況，即提高儲層靜態(tài)建模的精度，一直是國內(nèi)外學者研究的重點之一。

構(gòu)造建模是儲層靜態(tài)建模的主要手段之一，1989年MALLET［2］采用離散光滑插值方法實現(xiàn)了構(gòu)造建模；針對復雜構(gòu)造建模的需求，1998年EULER［3］提出了表征斷層和層面空間拓撲關(guān)系的解決方案；楊欽等［4］將Delaunary三角剖分算法應用于構(gòu)造建模；魏嘉等［5］提出了適用于復雜地質(zhì)體的曲面擬合技術(shù)，并發(fā)展了適應地層沉積特點的三維地層網(wǎng)格剖分方法；孫波等［6］進一步完善了復雜地質(zhì)界面的三維重構(gòu)方法，使其能更好地表征復雜地質(zhì)體的結(jié)構(gòu)。三維構(gòu)造建模技術(shù)較為成熟，隨著層序地層學的發(fā)展，結(jié)合層序地層建立表征地層尖滅、剝蝕和超覆等接觸關(guān)系的層面框架模型成為研究重點。在明確地層接觸關(guān)系之后，根據(jù)不同的地層接觸關(guān)系分別采用平行于頂面、平行于底面、在頂?shù)组g等比例插值這3種方法來垂向劃分層面框架模型。董建華等［7］以地震解釋的層序界面作為頂?shù)准s束，通過傾角體的引導進行層間的垂向劃分，得到了等時的層序界面。屬性模型表征了地下儲層及流體，其表征精度直接影響井位部署和開發(fā)效果，因此屬性建模方法也是國內(nèi)外學者的研究重點。HALDORSON等［8］于1984年提出了基于井數(shù)據(jù)的儲層隨機建模方法；1992年ARAKTINGI等［9］聯(lián)合地震和測井資料進行了屬性建模研究；此后，YANG等［10］和TJOLSEN等［11］研究了采用地震反演約束屬性建模的方法；DOYEN等［12］基于貝葉斯理論實現(xiàn)了地震約束的儲層參數(shù)預測；李緒宣等［13］基于地震屬性的空間變化代表了儲層的空間變化規(guī)律這一特點，利用與儲層參數(shù)相關(guān)性高的地震屬性來取代變差函數(shù)，提出了地震驅(qū)動建模方法，提高了井網(wǎng)稀疏時儲層參數(shù)預測的精度；郭智等［14］將地震反演自然伽馬和地質(zhì)約束條件加入到建模過程中，提高了井間儲層參數(shù)預測精度；胡勇等［15］采用地震反演數(shù)據(jù)約束儲層參數(shù)建模，以地震正演來驗證地質(zhì)模型的可靠性，提高了建模精度。

當層間存在與頂、底面不一致的起伏變化，尤其頂、底界面為不整合類型時，上述方法難以得到等時的層面框架模型，降低了構(gòu)造模型的精度。地震驅(qū)動建模方法提高了井網(wǎng)稀疏時儲層參數(shù)預測的精度，但是沒有解決三維地層網(wǎng)格模型不等時的問題。本文在地震驅(qū)動建模方法的基礎(chǔ)上，利用地震傾角體來引導垂向網(wǎng)格的劃分，建立等時的三維地層網(wǎng)格模型，以井數(shù)據(jù)為已知點，利用敏感地震屬性的變化規(guī)律來估計井間儲層參數(shù)，形成基于地震等時格架的傾角導向儲層靜態(tài)建模方法。理論模型的數(shù)據(jù)測試和海上油田的實際應用表明，該方法能進一步提高建模精度。

1 基本原理

1．1 傾角導向三維地層網(wǎng)格模型建立方法

利用相干估計傾角掃描計算方法計算傾角體［1617］。以傾角體為導向建立三維地層網(wǎng)格模型，選擇各小層種子點的時間（或深度）和局部傾角，根據(jù)不同點之間傾角值的關(guān)系計算出小層下一點的時間（或深度）值，采用這種逐點遞推的方式完成小層的導向追蹤。主要過程為：①在斷層模型的基礎(chǔ)上，以地震解釋的層面作為關(guān)鍵層面；②以關(guān)鍵層面作為頂、底趨勢面，以頂、底面間的傾角體為導向，對其內(nèi)部的小層進行追蹤內(nèi)插；③根據(jù)儲層表征的需要，在內(nèi)插小層基礎(chǔ)上，以層間的傾角體為導向，采用逐點遞推計算的方式從起始點開始追蹤小層內(nèi)的網(wǎng)格層，最后建立三維地層網(wǎng)格模型。

傾角導向小層追蹤的具體算法為：

式中：Zi，j，k表示（i，j）位置處第k層的深度；ΔX和ΔY分別為X和Y方向上的采樣間隔；Xdipi，j，k和Ydipi，j，k分別為X方向和Y方向（i，j）位置處第k層的視傾角。

1．2 基于地震等時格架的儲層靜態(tài)建模方法

反映相同地質(zhì)體的地震屬性和儲層參數(shù)聯(lián)系緊密［18］，因此對儲層參數(shù)敏感的地震屬性可以揭示儲層參數(shù)的變化［19－21］，再聯(lián)合已知井就可以定量預測儲層參數(shù)值?；诘卣鸬葧r格架的傾角導向儲層靜態(tài)建模方法的原理就基于此，如圖1所示。已知點1，2，3處的地震屬性（S1，Sn，S2）和點1，3處的測井數(shù)據(jù)（L1，L2），根據(jù)地震屬性的空間變化規(guī)律預測點2處的儲層參數(shù)值。在不等時的地層網(wǎng)格里提取數(shù)據(jù)時，地震屬性和測井數(shù)據(jù)之間的線性相關(guān)性差，預測的儲層參數(shù)值精度低。因此需要傾角導向建立的三維地層網(wǎng)格模型來約束，非等時地層內(nèi)的測井數(shù)據(jù)和地震屬性不參與該地層的建模，這樣建模結(jié)果更準確。

根據(jù)地震屬性的空間相關(guān)性：

圖1 基于地震等時格架傾角導向儲層靜態(tài)建模原理示意

式中：Sc代表未知點的地震屬性；Si代表已知井的井旁道地震屬性；ωi代表不同井點處的轉(zhuǎn)換函數(shù)。井點處的地震數(shù)據(jù)和測井數(shù)據(jù)分別從不同尺度反映儲層變化情況，當?shù)卣饘傩耘c儲層參數(shù)之間具有較好線性相關(guān)性時，有：

式中：Lc為未知點的儲層參數(shù)值；Li為已知井的儲層參數(shù)值。

由（2）式、（3）式可知，通過地震屬性求取反映儲層參數(shù)空間變化規(guī)律的轉(zhuǎn)換函數(shù)，結(jié)合已知井點的測井數(shù)據(jù)，就可以預測出未知點的儲層參數(shù)值，不需要求取變差函數(shù)。其中，轉(zhuǎn)換函數(shù)的求取是實現(xiàn)基于地震等時格架傾角導向儲層靜態(tài)建模的關(guān)鍵。

采用多道最小化均方誤差來求取轉(zhuǎn)換函數(shù)ωi。以兩口井為例介紹多道最小化均方誤差計算轉(zhuǎn)換函數(shù)的方法。這兩口井分別對應兩個轉(zhuǎn)換函數(shù)ω1，j和ω2，j（j＝0，1，…，m）。k時刻的誤差為：

平方誤差為：

總的平方誤差為：

當E2最小時有：

E2對ω1，j和ω2，j分別求偏導，得：

利用互相關(guān)和相關(guān)計算來化簡方程并求解，得到轉(zhuǎn)換函數(shù)ω1，j和ω2，j。多口井時上述過程同理擴展。

1．3 技術(shù)流程

在三維地層網(wǎng)格模型不等時的情況下，井震之間的線性相關(guān)性變差，相當于用地震屬性去預測非同一時期的儲層參數(shù)，預測精度下降。因此，需要通過傾角導向建立一個等時的三維地層網(wǎng)格模型，在其約束下，提取的地震屬性和測井數(shù)據(jù)才能反映同一時期沉積的儲層，敏感地震屬性才能真正反映儲層參數(shù)的變化，從而取代變差函數(shù)，有效避免變差函數(shù)在井較少或井網(wǎng)稀疏時難以得到穩(wěn)定結(jié)果的缺點［13］。因此，基于地震等時格架的傾角導向儲層靜態(tài)建模對井間儲層參數(shù)預測的精度要高于地質(zhì)統(tǒng)計學建模方法。主要步驟（圖2）如下：

1）根據(jù)斷層數(shù)據(jù)建立斷層模型，基于地震數(shù)據(jù)計算地震傾角體，在斷層模型基礎(chǔ)上，利用傾角體引導建立三維地層網(wǎng)格模型；

2）對測井數(shù)據(jù)進行質(zhì)量檢查，確保無誤后將建模用到的測井數(shù)據(jù)按照厚度加權(quán)的方法粗化到三維網(wǎng)格中；

3）確定需要建立的儲層參數(shù)模型，尋找對該參數(shù)最敏感的地震屬性，將地震屬性數(shù)據(jù)投影到三維地層網(wǎng)格模型中，為建模做好準備；

4）在一個等時的小層內(nèi)提取層內(nèi)的地震屬性與井數(shù)據(jù)，通過地震屬性計算轉(zhuǎn)換函數(shù)，利用該函數(shù)從井出發(fā)估計層內(nèi)未知點的儲層參數(shù)值；

5）根據(jù)步驟4）對下一個等時小層內(nèi)每個網(wǎng)格的儲層參數(shù)值進行預測，不斷循環(huán)，最終建立三維儲層參數(shù)模型。

圖2 基于地震等時格架的傾角導向儲層靜態(tài)建模技術(shù)流程

2 模型試驗

采用由一個背斜和一個向斜組成的褶皺構(gòu)造模型對本文方法進行測試。沿褶皺構(gòu)造方向劃分為100個網(wǎng)格，步長為10m；走向方向劃分為5個網(wǎng)格，步長為20m?？v向由19個小層組成，劃分為550個網(wǎng)格，步長為1m。設(shè)定理論波阻抗對理論孔隙度有很高的敏感性，將波阻抗作為基于地震等時格架傾角導向儲層靜態(tài)建模的地震數(shù)據(jù)。給每個小層填充孔隙度參數(shù)和波阻抗參數(shù)，從而得到理論孔隙度三維模型和理論波阻抗三維模型（圖3）。參與孔隙度建模的井數(shù)據(jù)從理論孔隙度模型中抽取，一共提取了4口井數(shù)據(jù)。

2．1 傾角導向三維地層網(wǎng)格模型的建立

在采用50Hz雷克子波與理論模型褶積正演得到的地震數(shù)據(jù)中，將解釋為褶皺頂、底的層位作為關(guān)鍵層面。在地震剖面中拾取穩(wěn)定連續(xù)的波峰為頂界面，拾取穩(wěn)定連續(xù)的波谷作為底界面，對每條測線都進行解釋?；诘卣饠?shù)據(jù)分別計算Inline方向和Xline方向的視傾角體，其中Xline方向沿走向方向，所以其視傾角值都為0。以解釋的頂、底層位作為趨勢面，基于Inline方向和Xline方向計算的傾角體進行導向追蹤細分，得到細分后的網(wǎng)格層。地震反射同相軸具有等時意義，將細分后的網(wǎng)格層與地震同相軸進行疊合對比，分析其等時性。與地震同相軸一致的網(wǎng)格層等時性好；與地震同相軸趨勢不一致，發(fā)生斜交的網(wǎng)格層不等時。圖4a和圖4b分別給出了采用傾角導向方法和線性等分方法細分的小層與地震剖面疊合結(jié)果。由圖4可以看到，傾角導向細分的網(wǎng)格層等時性很好，而線性等分的網(wǎng)格層在局部不等時。

圖3 理論三維模型示意

圖4 細分小層與地震剖面疊合結(jié)果

2．2 基于地震等時格架的儲層靜態(tài)模型的建立

波阻抗從與網(wǎng)格無關(guān)的反演中得到，將波阻抗和井數(shù)據(jù)重采樣到傾角導向的三維網(wǎng)格中，基于波阻抗的空間變化來預測孔隙度，建立三維孔隙度模型。通過對比預測孔隙度與理論孔隙度來評價預測精度。圖5給出了理論波阻抗模型、理論孔隙度模型以及不同地層模型約束的孔隙度模型剖面。由圖5可以看出，基于地震等時格架傾角導向儲層靜態(tài)建模的預測孔隙度分布和理論孔隙度基本一致，預測精度很高。線性等分方法建立的三維地層網(wǎng)格模型存在局部穿時，在其約束下地震驅(qū)動建立的孔隙度模型在地層模型穿時位置預測孔隙度值與理論孔隙度值差異較大（圖5中黑色框線內(nèi)），預測精度較低。

圖5 理論模型與不同地層模型約束的預測孔隙度模型剖面

圖6 檢驗井處基于不同地層模型的預測孔隙度與理論孔隙度對比

采用不參與建模的檢驗井來檢驗孔隙度模型的預測精度。將檢驗井設(shè)置在x＝200m處，對比該檢驗井的預測孔隙度與理論孔隙度。圖6給出了基于不同地層模型的預測孔隙度與理論孔隙度。由圖6可以看出，在檢驗井處，地震等時格架傾角導向儲層靜態(tài)建模預測的孔隙度與理論孔隙度的相對誤差在5%以內(nèi)，兩者幾乎相同；線性等分地層格架地震驅(qū)動模型的預測孔隙度與理論孔隙度的相對誤差局部超過了10%。這是因為線性等分三維地層網(wǎng)格模型存在局部不等時，導致了孔隙度模型預測精度下降。基于地震等時格架的傾角導向儲層靜態(tài)建模建立了等時的三維地層網(wǎng)格模型，建立的孔隙度模型準確地預測了理論孔隙度值，建模的精度有了進一步提高。

3 實際應用

A油田位于渤海灣地區(qū)，為海上稠油油田，主要目的層為新近系明化鎮(zhèn)組Nm I—Nm IV，總共有4個油組。儲層沉積類型為河流相沉積，儲層橫向變化快，開發(fā)井井距在350～450m，相對陸上油田井距較大。傳統(tǒng)地質(zhì)統(tǒng)計學建模方法對井間儲層的預測存在較大的不確定性，難以精細表征井間儲層參數(shù)的變化，無法滿足油田中后期開發(fā)調(diào)整的需求。采用傾角導向建立一個等時的三維地層網(wǎng)格模型，在其約束下，基于敏感地震屬性計算轉(zhuǎn)換函數(shù)來表征儲層參數(shù)的變化［22］。在井數(shù)據(jù)的基礎(chǔ)上，結(jié)合轉(zhuǎn)換函數(shù)估計儲層參數(shù)，建立三維地質(zhì)模型。

選擇A油田解釋了孔隙度的46口井參與孔隙度模型的建立。分析哪種地震屬性對孔隙度的變化最敏感，以便確定采用此種地震屬性進行基于地震等時格架的傾角導向儲層靜態(tài)建模。結(jié)果表明，目的層段低速度、低密度對應砂巖，高速度、高密度對應泥巖。利用波阻抗可以區(qū)分砂、泥巖，而泥巖是低孔隙度，砂巖是高孔隙度，所以高孔隙度對應低阻抗，低孔隙度對應高阻抗。46口井的數(shù)據(jù)都體現(xiàn)了這一特點。因此采用油田的地震反演波阻抗數(shù)據(jù)來進行基于地震等時格架的傾角導向儲層靜態(tài)建模。其中，A10井作為檢驗井，不參與建模，而用于檢驗模型與實際情況的吻合程度。

基于地震數(shù)據(jù)計算Inline方向和Xline方向的傾角體。選擇目的層段穩(wěn)定、連續(xù)的地震反射同相軸，解釋8個關(guān)鍵層位。這些關(guān)鍵層位都是等時界面，作為頂、底面約束，分別利用傾角導向方法和線性等分方法追蹤網(wǎng)格層，形成三維網(wǎng)格，如圖7所示。從圖7可以看到，傾角導向方法基于地震反射同相軸進行追蹤細分來建立三維地層網(wǎng)格模型，得到的網(wǎng)格層面與地震反射同相軸一致，具有良好的等時性。線性等分方法通過等比例內(nèi)插得到內(nèi)部的網(wǎng)格層，內(nèi)插出來的網(wǎng)格層與地震反射同相軸存在交叉現(xiàn)象，所建立的三維地層網(wǎng)格模型不完全等時。

圖7 采用不同方法建立的三維地層網(wǎng)格模型

將參與建模的測井數(shù)據(jù)和地震屬性粗化到三維網(wǎng)格中。在網(wǎng)格化的等時小層內(nèi)，提取已知井的井旁地震屬性，采用多道最小化均方誤差方法計算反映孔隙度空間變化的轉(zhuǎn)換函數(shù)。基于計算得到的轉(zhuǎn)換函數(shù)，從小層內(nèi)已知井的孔隙度出發(fā)，估計該小層內(nèi)未知點的孔隙度，最終建立三維孔隙度模型?；诓ㄗ杩沟目臻g變化規(guī)律對孔隙度進行預測，所以孔隙度模型的空間變化與波阻抗的空間變化一致（圖8）。從模型提取的Nm31小層等時切片（圖9）清晰地反映了河流相儲層的平面展布特征，表明基于地震等時格架傾角導向儲層靜態(tài)建模建立的三維地層網(wǎng)格更等時?？紫抖饶Ｐ偷姆直媛拭黠@高于地震反演波阻抗，因為基于地震等時格架傾角導向儲層靜態(tài)建模綜合了測井數(shù)據(jù)垂向分辨率高的優(yōu)勢，所以對細節(jié)的刻畫比波阻抗更精細，對砂體的空間展布特征描述更清晰。

圖8 波阻抗模型（a）與基于地震等時格架傾角導向儲層靜態(tài)建模預測的孔隙度模型（b）

圖9 Nm31小層波阻抗等時切片（a）與基于地震等時格架傾角導向儲層靜態(tài)建模預測的孔隙度等時切片（b）

用檢驗井數(shù)據(jù)來驗證地震等時格架傾角導向儲層靜態(tài)建模方法對儲層參數(shù)的預測精度。A10井作為檢驗井不參與建模，沿檢驗井的井軌跡提取模型孔隙度，并與測井孔隙度進行對比，評價兩種孔隙度之間的相對誤差。相對誤差越小，表明建模方法的預測精度越高。

圖10給出了評價井A10井基于不同方法預測得到的模型孔隙度與測井孔隙度的對比結(jié)果。從圖10可以看到，地震等時格架傾角導向儲層靜態(tài)建模方法的模型孔隙度變化趨勢與測井孔隙度相同，兩者孔隙度值的相對誤差小于10%。線性等分劃分網(wǎng)格方法的模型孔隙度與測井孔隙度吻合程度要差一些，12%的孔隙度值相對誤差大于10%。配置協(xié)克里金序貫高斯方法模擬的孔隙度值在局部與測井孔隙度值差異較大，28%的孔隙度值的相對誤差大于10%。因此，配置協(xié)克里金序貫高斯方法對孔隙度的預測精度最差。這是因為配置協(xié)克里金序貫高斯方法以井數(shù)據(jù)為儲層參數(shù)預測的核心，敏感地震屬性僅起到輔助作用；而基于地震等時格架傾角導向儲層靜態(tài)建模方法能夠充分利用敏感地震屬性來表征儲層參數(shù)的空間變化，地震屬性與井數(shù)據(jù)都是儲層參數(shù)預測的核心。因此基于地震等時格架傾角導向儲層靜態(tài)建模方法對儲層的表征精度更高，比地質(zhì)統(tǒng)計學方法具有優(yōu)越性，采用傾角導向的方法自動追蹤形成的網(wǎng)格層，比常規(guī)方法建立的網(wǎng)格層更接近于等時?；诘卣鸬葧r格架傾角導向儲層靜態(tài)建模方法比常規(guī)方法建立的地層模型約束的地震驅(qū)動建模精度更高，在井網(wǎng)稀疏的條件下能更好地表征儲層的非均質(zhì)性，為油田開發(fā)井網(wǎng)部署和開發(fā)方案的調(diào)整優(yōu)化提供了堅實的技術(shù)支持。

圖10 評價井A10井基于不同方法預測得到的模型孔隙度與測井孔隙度對比

4 結(jié)論

當層間存在與頂、底面不一致的起伏變化時，常規(guī)方法劃分的垂向網(wǎng)格層等時性較差。嘗試采用地震傾角體來引導垂向網(wǎng)格層的劃分，建立等時的三維地層網(wǎng)格模型，更好地反映了儲層的層序結(jié)構(gòu)和沉積特征，對提高屬性建模的精度具有現(xiàn)實意義。在此基礎(chǔ)上，利用敏感地震屬性來求取轉(zhuǎn)換函數(shù)，以井數(shù)據(jù)為已知點，采用轉(zhuǎn)換函數(shù)來估計井間儲層參數(shù)，不再采用變差函數(shù)來預測儲層參數(shù)的變化，降低了井網(wǎng)稀疏時地質(zhì)統(tǒng)計學方法對井間儲層參數(shù)預測的不確定性，形成了一種基于地震等時格架的傾角導向儲層靜態(tài)建模方法。

理論模型測試和實際應用表明，傾角導向三維地層網(wǎng)格模型比常規(guī)方法建立的三維地層網(wǎng)格模型更接近等時，提高了屬性建模的精度。在井網(wǎng)稀疏的條件下，與傳統(tǒng)地質(zhì)統(tǒng)計學方法相比，基于地震等時格架的傾角導向儲層靜態(tài)建模方法能更精細地描述儲層的非均質(zhì)性，降低了井間儲層預測的不確定性。該方法在地震屬性與儲層參數(shù)線性相關(guān)性好的地區(qū)具有較好的推廣應用價值。

［1］ 吳勝和，李宇鵬．儲層地質(zhì)建模的現(xiàn)狀與展望［J］．海相油氣地質(zhì)，2007，12（3）：53－60

WU S H，LI Y P．Reservoir modeling：current situation and development prospect［J］．Marine Origin Petroleum Geology，2007，12（3）：53－60

［2］ MALLET J L．Discrete smooth interpolation［J］．ACMTransactions on Graphics，1989，8（2）：121－144

［3］ EULER N．A new tool to seal 3Dearth model：a cut with constraint［J］．Expanded Abstracts of 68thAnnual Internat SEG Mtg，1998：321－324

［4］ 楊欽，徐永安，陳其明，等．三維約束Delaunay三角化的研究［J］．計算機輔助設(shè)計與圖形學學報，2000，12（8）：590－594

YANG Q，XU Y A，CHEN Q M，et al．Research on 3D constrained delaunay triangulation［J］．Journal of Computer－Aided Design ＆Computer Graphics，2000，12（8）：590－594

［5］ 魏嘉，唐杰，岳承祺，等．三維地質(zhì)構(gòu)造建模技術(shù)研究［J］．石油物探，2008，47（4）：319－326

WEI J，TANG J，YUE C Q，et al．Study of 3－D geological structure model building［J］．Geophysical Prospecting for Petroleum，2008，47（4）：319－326

［6］ 孫波，劉大安．復雜地質(zhì)界面三維重構(gòu)與評價方法［J］．巖石力學與工程學報，2015，34（3）：556－564

SUN B，LIU D A．Three－dimensional reconstruction of complex geological interfaces and its evaluation method［J］．Chinese Journal of Rock Mechanics and Engineering，2015，34（3）：556－564

［7］ 董建華，范廷恩，高云峰，等．真地層切片的拾取及應用［J］．石油地球物理勘探，2010，45（增刊1）：150－153

DONG J H，F(xiàn)AN T E，GAO Y F，et al．Pickup and ap－plication of the true stratal slice［J］．Oil Geophysical Prospecting，2010，45（S1）：150－153

［8］ HALDORSON H，LAKE L A．New approach to shale management in field scale simulation models［J］．Society of Petroleum Engineers Journal，1984，24（4）：447－457

［9］ 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［J］．SPE Technical Conference ＆Exhibition，1992：913－926

［10］ YANG C T，CHOPRA A K．Integrated geostatistical reservoir description using petrophysical，geological，and seismic data for Yacheng 13－1gas field［J］．SPE Technical Conference ＆Exhibition，1995：357－372

［11］ TJOLSEN C B，JOHNSEN G，HALVORSEN G，et al．Seismic data can improve the stochastic facies modeling［J］．SPE Formation Evaluation，1996，11（3）：141－146

［12］ DOYEN P M，PSAILA D E，DEN BOER L D，et al．Reconciling data at seismic and well log scales in 3D earth modeling［J］．SPE Technical Conference ＆Exhibition，1997：465－474

［13］ 李緒宣，胡光義，范廷恩，等．基于地震驅(qū)動的海上油氣田儲層地質(zhì)建模方法［J］．中國海上油氣，2011，23（3）：143－147

LI X X，HU G Y，F(xiàn)AN T E，et al．A method of geological reservoir modeling based on seismic drive for offshore fields［J］．China Offshore Oil and Gas，2011，23（3）：143－147

［14］ 郭智，孫龍德，賈愛林，等．辮狀河相致密砂巖氣藏三維地質(zhì)建模［J］．石油勘探與開發(fā)，2015，42（1）：76－83

GUO Z，SUN L D，JIA A L，et al．3Dgeological modeling for tight sand gas reservoir of braided river facies［J］．Petroleum Exploration and Development，2015，42（1）：76－83

［15］ 胡勇，于興河，李勝利，等．應用地震正反演技術(shù)提高地質(zhì)建模精度［J］．石油勘探與開發(fā)，2014，41（2）：190－197

HU Y，YU X H，LI S L，et al．Improving the accuracy of geological model by using seismic forward and inversion technique［J］．Petroleum Exploration and Development，2014，41（2）：190－197

［16］ MARFURT K J，KIRLIN R L，F(xiàn)ARMER S L．3Dseismic attributes using a semblance－based coherency algorithm［J］．Geophysics，1998，63（4）：1150－1165

［17］ MARFURT K J．Robust estimates of 3Dreflector dip and azimuth［J］．Geophysics，2006，71（4）：29－40

［18］ 劉偉，尹成，王敏，等．河流相砂泥巖薄互層基本地震屬性特征研究［J］．石油物探，2014，53（4）：468－476

LIU W，YIN C，WANG M，et al．Study on the characteristics of basic seismic attributes in fluvial sand－shale interbedded layers［J］．Geophysical Prospecting for Petroleum，2014，53（4）：468－476

［19］ 張營革．能量半衰時屬性在濁積巖儲層預測中的應用研究［J］．石油物探，2013，52（6）：662－668

ZHANG Y G．The application of energy half－time attribute to turbidite reservoir prediction［J］．Geophysical Prospecting for Petroleum，2013，52（6）：662－668

［20］ 畢俊鳳，顧漢明，劉書會，等．河道砂體地震響應特征及影響因素分析——以墾東1地區(qū)館陶組上段河道砂為例［J］．石油物探，2013，52（1）：97－103

BI J F，GU H M，LIU S H，et al．Analysis on seismic response characteristics of channel sands and its influence factors：case study of Upper Guantao Formation in KD1area［J］．Geophysical Prospecting for Petroleum，2013，52（1）：97－103

［21］ 唐湘蓉，蔡涵鵬，賀振華．地震波高頻信息在薄層砂體預測中的應用［J］．石油物探，2012，51（3）：244－250

TANG X R，CAI H P，HE Z H．Thin－bed sand body prediction based on seismic wave high－frequency information［J］．Geophysical Prospecting for Petroleum，2012，51（3）：244－250

［22］ 張建寧，韓文功．東營凹陷濁積砂體疊后地震屬性與儲層物性相關(guān)性研究［J］．石油物探，2012，51（2）：204－212

ZHANG J N，HAN W G．Correlative analysis on poststack seismic attributes and physical properties of reservoir in turbidite sandstone of Dongying Depression［J］．Geophysical Prospecting for Petroleum，2012，51（2）：204－212

（編輯：陳 杰）

A method of dip steering reservoir static modeling based on seismic isochronal stratigraphic framework

LE Jing，WANG Hui，F(xiàn)AN Ting’en，ZHAO Weiping，F(xiàn)AN Pengjun
（CNOOC Research Institute，Beijing100028，China）

When dividing vertical layer is not isochronal and the oil field has sparse well pattern，the accuracy of reservoir parameter prediction in reservoir static modeling will be reduced．Especially in the offshore oil field with a little drilling，it is difficult to establish the reservoir parameter model with high accuracy based on geological statistics method．In order to establish a better isochronous stratigraphic framework and improve the accuracy of reservoir parameter prediction，the method of dip steering reservoir static modeling based on seismic isochronal stratigraphic framework is studied．For the method we establish three－dimensional strata model by dividing vertical layers under the steering of seismic dip cubes，predict reservoir spatial variation from sensitive seismic attributes and combine with well data to establish three dimensional model of reservoir parameters．Comparing with conventional structural modeling methods，stratigraphic framework model under the steering of seismic dip cubes is established in the method．At the same time，the transfer function calculated from seismic attributes to establish reservoir parameter model，which is different from geological statistics method that takes variation function to build reservoir model．And the method can improve the accuracy of reservoir model when the oil field has sparse well pattern．The practical application of an offshore oil field with sparse well pattern shows the method establishes a more isochronous stratigraphic framework than conventional modeling method and further improves the accuracy of the geological modeling．

dip cube，angle guidance，seismic drive，isochronal stratigraphic framework，reservoir static modeling

P631

A

1000－1441（2017）03－0449－10

10．3969／j．issn．1000－1441．2017．03．015

樂靖，王暉，范廷恩，等．基于地震等時格架的傾角導向儲層靜態(tài)建模方法［J］．石油物探，2017，56（3）：449－458

LE Jing，WANG Hui，F(xiàn)AN Ting’en，et al．A method of dip steering reservoir static modeling based on seismic isochronal stratigraphic framework［J］．Geophysical Prospecting for Petroleum，2017，56（3）：449－458

2016－04－11；改回日期：2016－07－26。

樂靖（1984—），男，碩士，工程師，主要從事地震綜合解釋方法研究。

國家科技重大專項“海上開發(fā)地震關(guān)鍵技術(shù)及應用研究”（2011ZX05024－001）資助。

This research is financially supported by the National Science and Technology Major Project of China（Grant No．2011ZX05024－001）．