桑塔木地區三疊系深部地震薄儲集層預測

汪佳蓓，黃捍東，張文珠

（1.中國地質大學（北京）能源學院，北京 100083；2.中國石油大學（北京）油氣資源與探測國家重點實驗室，北京 102249；3.中國石油大學（北京）提高采收率研究院，北京 102249；4.中海石油深圳分公司研究院，廣東 廣州 510240）

桑塔木地區三疊系深部地震薄儲集層預測

汪佳蓓1，黃捍東2，3，張文珠4

（1.中國地質大學（北京）能源學院，北京 100083；2.中國石油大學（北京）油氣資源與探測國家重點實驗室，北京 102249；3.中國石油大學（北京）提高采收率研究院，北京 102249；4.中海石油深圳分公司研究院，廣東 廣州 510240）

深部隱蔽油氣藏逐步成為現階段油氣勘探的重難點，其具有埋深大、地震信號弱、構造復雜等特征。埋藏深度的增大使得地震反射信息量減少，識別儲集層尤其是薄儲集層的難度加大。文中采用地震相控非線性隨機反演方法，融入測井和層序界面等地質類信息，對桑塔木油田三疊系近5 000 m埋深的薄儲集層進行了鉆前預測，成功刻畫了深層薄儲集層的縱橫向連通性及平面展布特征。通過未參與反演井預測效果分析，少井反演與多井反演效果對比，預測吻合率統計分析和新鉆井與預測結果對比等方面來評價本次地震資料儲集層預測的可信度。預測結果與實鉆情況吻合程度較好，說明該方法在深層薄儲集層預測中有很好的應用前景。

深層；薄儲集層；相控反演；桑塔木油田

0　引言

隨著油氣勘探研究的不斷深入，研究重點已經從尋找淺層、構造油氣藏向深層、巖性油氣藏方向發展［1-3］。深部地層一般埋深超過4 000 m，地質條件復雜，地震反射弱、主頻低、品質差，因此，深部地震儲集層，特別是薄層預測有很大難度。

對于薄砂體儲層預測，基于20世紀80年代提出的Widess準則［4］，能識別大于1/4地震波長的薄儲層。此后許多研究發現，薄層描述可以突破傳統分辨率的極限［5-13］。

塔里木盆地桑塔木油田三疊系地層埋深近5 000m，目的層油組含油砂巖的單層厚度在10 m左右，是典型的深部薄層油藏。常規反演方法在該地區的應用遇到了2個問題：1）砂泥巖聲波測井速度差異小，反演剖面上儲集層特征得不到清晰展現；2）含油砂巖厚度低于地震資料分辨率，無法有效地識別薄層。以上問題導致三疊系油藏的地質模型認識不清，井間油層砂體規模、縱橫向連通性、平面展布及尖滅等不明。

針對以上難題，本文采用地震相控非線性隨機反演方法，對塔里木盆地桑塔木地區三疊系地層進行儲集層預測。

1　地震相控非線性隨機反演方法原理

地震相控非線性隨機反演是在反演的過程中，融入地震相界面和沉積特征等地質類信息，在相模型的約束和指導下，結合鉆探得到的測井資料，采用非線性隨機反演算法，在計算過程中逐道進行外推反演［14］。具體實現過程可以分為地震相控約束外推計算、隨機模擬處理和非線性反演3個部分。

1.1地震相控約束外推計算

地震相是地震反射時間剖面上所表現出的反射波的面貌，是沉積相在地震剖面上的反映［15］。在反演之前劃分出相界面可以用來控制反演的低頻趨勢。反演中通過“相界面”控制的同時，結合“井約束”外推算法，可很好地解決反演所需的低頻趨勢問題，又能克服反演中構建的初始模型精度對最終反演結果的影響等。

實際應用中，首先，利用鉆井、綜合錄井資料編制工區內井的層序劃分與地層界面解釋對比圖；然后，利用這個結果，在對應的連井地震剖面上解釋出層序界面，建立三維層序和相控模型，從而可以將該層序和相控模型運用到反演的約束條件中去，將宏觀模型的特征融入到反演中，為地震相控約束反演奠定約束條件。

由于地下地質情況具有強不確定性，在相控外推計算過程中，可以運用多項式相位時間擬合的方法來建立地震道間的外推關系。具體實現方法如下：在層序和相控模型控制的特定地震時窗范圍約束下，將從測井信息中或井旁地震道中獲取的先驗模型參數向量從井點處，開始往外推計算。

1.2隨機模擬處理

考慮到地質條件的隨機特性，首先利用測井巖石物理統計分析，建立速度、密度等模型參數點之間的統計關系；然后采用高斯隨機函數把模型參數作為區域化變量來進行隨機模擬。

具體需要利用變差函數來描述空間數據場中數據之間的相互關系，從而建立起空間儲層參數點之間的統計相關函數關系式［16-17］。變差函數描述的是不同位置變量的不相似性。

1.3非線性反演

目前廣泛應用于生產的線性反演算法，計算速度快，穩定性強。然而，在實際情況中遇到的反演問題基本都屬于非線性的，如果直接轉成線性問題來求解，會導致最終的反演結果精度低，反映不出薄層特征，計算結果也容易陷入局部極小等問題［18-19］。非線性反演算法直接求解非線性問題，解的精度高，但穩定性差。本文采用非線性最優化算法，求解過程采用截斷奇異值分解法并加入約束條件來提高解的穩定性與精度。

若假定密度為常數，則波阻抗反演就直接轉換為速度反演［20］。先定義目標函數為式（1），需要求解m使目標函數在最小二乘意義下取得極小值。

式中：m為速度；uobs為實際地震數據；umbd為模型響應。

然后，將式（1）在初始模型響應mi處按Taylor公式展開，忽略二次以上的高階項得：

式中：D為f（m）的二次偏導數矩陣。

令m=mi+1，則有：

將式（3）對mi+1求一階導數，并令該導數為0，可得迭代公式為

在求取D-1時往往存在計算量大、矩陣病態等問題，為了增強解的穩定性，本文運用了截斷奇異值分解法，避免直接求逆存在。同時在迭代過程中，通過加入最大最小值、參數包絡等約束，進一步提高解的穩定性和精度。通過式（4）迭代計算，可得到最終反演速度體。

2　實際應用

2.1地質背景

桑塔木地區位于塔里木盆地塔北隆起輪南低凸起的斜坡帶上，三疊系地層是主力油氣產層之一，包含3套含油儲集層，分別為上三疊統黃山街組T1油組、中三疊統克拉瑪依組T2油組和T3油組，埋深為4 600～5 000 m，除T3油組含油薄砂體儲集層的厚度有部分大于10 m外，其他T1油組和T2油組含油砂層的厚度均較薄，一般在3～5 m。該區沉積體系已經進行了大量研究，整體為辮狀河三角洲沉積模式［21-22］，但通過地震地質重新研究和認識，發現該區還存在部分曲流河三角洲沉積模式。該模式位于克拉瑪依組，橫向變化大，在地震剖面上表現為強振幅、中低頻、中—連續性特征。

由于目標薄儲集層中單砂體厚度很薄，只能在測井曲線上看到測井響應，而在地震剖面上無法看到單砂體地震響應。因此，本文嘗試采用地震相控非線性隨機反演方法對塔里木油田三疊系進行儲集層預測，期望在反演剖面中能夠預測出多套薄層砂體，解決該區深層薄儲集層預測的難題。

2.2速度反演效果分析

2.2.1未參與反演井預測結果

LN23井約束反演的反演速度連井剖面見圖1，LN23井為參與反演的約束井，JF121井（距LN23井2.9 km）和JF123井（距LN23井9.5 km）為未參與反演的驗證井?？梢钥闯?，T1，T2油組砂層厚度大致分布在2.5～6.5 m，T3油組砂層厚度主要分布在7～14 m。反演結果中的3個油組特征明顯，與該區完鉆井吻合。JF121井和JF123井上3個油組的厚、薄儲集層均得到了較好的預測，與測井解釋成果吻合程度較高，說明該反演方法在該區預測效果很好。

2.2.2少井反演與多井反演效果對比

檢驗反演方法預測性的有效方法之一，是不論初始約束井多或少，只要基于相同的地震資料其最終反演結果應該趨于一個方向［23-24］。圖2a為26口約束井（不包括LN48井）的反演剖面，圖2b為46口約束井（包括LN48井）的反演剖面。通過對比發現，兩者在地層結構、巖性組合等方面非常相近，只是在一些微觀細節處存在差異。這說明參與反演的約束井的數量對反演結果的影響不大，對初始模型的依賴程度也不高，反演結果穩定可靠。

2.2.3預測吻合率統計分析

工區內，選取未參與反演的50口井的綜合測井解釋及巖性錄井資料，與利用反演剖面預測的砂體厚度進行對比統計。50口井預測結果中有6口井不符合實鉆情況，吻合率為88%。說明地震預測與實際鉆井之間差異性較小，可將該反演成果用于儲集層分布預測研究。

2.2.4新鉆井驗證結果分析

反演工作完成后，在工區北部LN30井東約700 m處鉆了一口新井LN302h，圖3為LN30井至LN302h井在LN302h井鉆探之前的反演剖面。

圖中：1）老井LN30井T3油組1小層井段為4 869～4 875 m的油層向LN302h井方向，速度降低，趨于泥巖速度?？梢悦黠@看到1小層連井處，反演剖面從暖色偏高值區域漸變為冷色偏低值區域，反演剖面上含油砂體的厚度逐漸變薄直至尖滅，盡管砂泥巖聲波測井速度差異小，反演剖面上依然可以看到砂體尖滅處的細微差別，較好地解決了砂泥巖常規反演區分不清的問題。2）LN30井2小層井段4 873.5～4 882.5 m橫向較穩定。LN302h井完鉆結果顯示，在1小層井段4 864.0～4 866 m鉆遇2 m差油層，在2小層4 870.0～4884 m鉆遇14 m水層，2小層之間（即4 866～4 870 m）為4 m的泥巖隔層，預測情況與實際鉆探結果吻合。新井預測結果表明，本文研究方法可用于新區滾動勘探指導新井井位的部署。

2.3儲集層平面特征預測

T1油組1小層砂體厚度見圖4。該小層一直是桑塔木地區主力產層，50多口井的鉆探顯示均有油氣。

從圖4中可以看出，砂體小層分布可分為東部及西北部2個體系：東南部及南部向北西發育的體系可看作東部物源的一支；西部自西向東的2個垛體及自西北部進入工區的小垛體構成了西部體系。相比而言，東部體系占絕對優勢，尤其是進入研究區中部的一支，進一步分為多個分支的大片砂體。砂體展布具有典型的曲流河三角洲平原及三角洲前緣沉積特點，符合該層實際沉積特征，為井位部署和后期鉆探提供依據。

綜上所述，相控反演的預測結果很好地突顯了儲集層的平面分布形態、厚度展布特征及沉積相特征，為本區深部巖性油氣藏的勘探部署提供了重要的依據。

3　結論

1）地震相控非線性隨機反演方法很好地融合了沉積層序等地質類的信息，反演的砂體預測結果更加符合區域沉積背景，解釋性和預測性更強。

2）速度反演的效果對比分析表明，該反演方法不過分依賴于參與反演的約束井數量和初始模型，約束井參與反演前后對比圖非常一致，且新鉆井的實鉆結果和預測結果吻合度很好，說明反演方法的可靠程度高，該方法優于常規反演方法，反演結果能為深部油氣藏的勘探部署提供可靠依據。

3）深部薄層預測一直是儲集層預測的難點，本文的反演方法清晰預測出了近5 000 m地層處的薄儲集層，為今后在深層、薄層以及巖性尖滅體等隱蔽油氣藏的勘探提供了解決思路。

［1］郝芳，鄒華耀，方勇.隱蔽油氣藏研究的難點和前沿［J］.地學前緣，2005，12（4）：581-488.

［2］張厚福，徐兆輝.從油氣藏研究的歷史論地層-巖性油氣藏勘探［J］.巖性油氣藏，2008，20（1）：114-123.

［3］李慶忠.巖性油氣藏地震勘探若干問題討論（Ⅰ）［J］.巖性油氣藏，2008，20（2）：1-5.

［4］WIDESS M.How thin is a thin bed［J］.Geophysics，1973，38（6）：1176-1180.

［5］KALLWEIT R，WOOD L.The limits of resolution of zero-phase wavelets［J］.Geophysics，1982，47（7）：1035-1046.

［6］ZENG H.How thin is a thin bed an alternative perspective［J］.The Leading Edge，2009，28（10）：1192-1197.

［7］ZENG H，BACKUS M.Interpretive advantages of 90°phase wavelets （Part 1）：seismic applications［J］.Geophysis，2005，70（3）：7-15.

［8］ZENG H，BACKUS M.Interpretive advantages of 90°phase wavelets （Part 2）：seismic applications［J］.Geophysis，2005，70（3）：17-24.

［9］COGGINS J.Thinner than expected？Don′t blame the seismic［J］.The Leading Edge，2011，30（8）：858-860.

［10］徐偉慕，郭平，胡天躍.薄互層調諧與分辨率分析［J］.石油地球物理勘探，2013，48（5）：750-757.

［11］楊飛，黃誠，李鵬飛，等.深層砂泥巖薄互層的地震識別方法［J］.重慶大學學報，2014，37（9）：136-143.

［12］黃捍東，張如偉.地震非線性隨機反演方法在陸相薄砂巖儲層預測中的應用［J］.石油學報，2009，30（3）：340-386.

［13］張志偉，王春生，林雅平，等.地震相控非線性隨機反演在阿姆河盆地A區塊碳酸鹽巖儲層預測中的應用［J］.石油地球物理勘探，2011，46（2）：304-310.

［14］黃捍東，羅群，付艷，等.地震相控非線性隨機反演研究與應用［J］.石油地球物理勘探，2007，42（6）：694-698.

［15］曹卿榮，李佩，孫凱，等.應用地震屬性分析技術刻畫河道砂體［J］.巖性油氣藏，2007，19（2）：93-96.

［16］慎國強，孟憲軍，王玉梅，等.隨機地震反演方法及其在埕北35井區的應用［J］.石油地球物理勘探，2004，39（1）：75-81.

［17］高金玉，馬明福，沙雪梅，等.隨機模擬方法在阿爾及利亞Oued Mya盆地的應用［J］.巖性油氣藏，2009，21（3）：86-89.

［18］劉全新，高建虎，董雪華.儲層預測中的非線性反演方法［J］.巖性油氣藏，2007，19（1）：81-85.

［19］王輝，張玉芬.基于模型的疊前數據多參數非線性反演［J］.巖性油氣藏，2008，20（2）：108-113.

［20］姜秀娣，魏修成，黃捍東，等.部分疊加資料縱橫波速度測井約束反演［J］.地球物理學進展，2009，24（1）：254-262.

［21］呂雪雁，朱筱敏，申銀民，等.塔里木盆地臺盆區三疊系層序地層研究和有利勘探區預測［J］.石油勘探與開發，2002，29（1）：32-35.

［22］龐雯，鄭俊茂.輪南地區三疊系辮狀河三角洲沉積儲層特征［J］.西南石油大學學報（自然科學版），2008，30（1）：57-62.

［23］劉文嶺.地震約束儲層地質建模技術［J］.石油學報，2008，29（1）：64-68.

［24］袁曉宇，李映濤，張哨楠，等.塔西南玉北地區鷹山組縫洞型儲層裂縫檢測技術及應用［J］.石油實驗地質，2015，37（3）：394-401.

（編輯楊會朋）

Seismic thin reservoir prediction in deep zones:a case study of Triassic Formation in Tarim Basin

WANG Jiabei1，HUANG Handong2，3，ZHANG Wenzhu4
（1.School of Energy Resources，China University of Geosciences，Beijing 100083，China；2.State Key Laboratory of Petroleum Resources and Prospecting，China University of Petroleum，Beijing 102249，China；3.EOR Research Institute，China University of Petroleum，Beijing 102249，China；4.Research Institute，Shenzhen Branch of CNOOC Ltd.，Guangzhou 510240，China）

As the development of oil and gas exploration，people are gradually focusing on the prospecting of reservoir in deep zones. Seismic reflection information will decrease as the depth increases.The identification of reservoir especially thin-layer reservoir becomes much more difficult.In this paper，the seismic phase-controlled non-linear random inversion method is used to predict the Triassic Formation in Sangtamu Area，where depth is nearly 5000 m，before drilling by taking logging，geological information as constraints.The lateral connectivity and plane distribution characteristics in deep thin reservoir are successfully portrayed.Through the contrast of the inversion results in different constraints such as no drilled-well constraint and less well constraints，speed inversion profile characteristics are analyzed.With the new drilling information to authenticate the inversion results，which shows that the prediction coincides with the real drilling situation and the seismic prediction method has good creditability and applicability in the process ofdeep reservoir exploration and development.

deep layer；thin reservoir；phase-controlled inversion；Sangtamu Oilfield

國家重點基礎研究發展計劃（“973”計劃）項目“深部有效碎屑巖儲層成因機制與發育模式”（2011CB201104）

TE132.1+4；P631

A

10.6056/dkyqt201604009

2015-10-10；改回日期：2016-05-07。

汪佳蓓，女，1986年生，在讀博士研究生，主要從事復雜油氣藏儲層預測研究。E-mail：lemon357@126.com。

引用格式：汪佳蓓，黃捍東，張文珠.桑塔木地區三疊系深部地震薄儲集層預測［J］.斷塊油氣田，2016，23（4）：447-450.

WANG Jiabei，HUANG Handong，ZHANG Wenzhu.Seismic thin reservoir prediction in deep zones：a case study of Triassic Formation in Tarim Basin［J］.Fault-Block Oil&Gas Field，2016，23（4）：447-450.