汪佳蓓，黃捍東，張文珠

（1.中國(guó)地質(zhì)大學(xué)（北京）能源學(xué)院，北京 100083；2.中國(guó)石油大學(xué)（北京）油氣資源與探測(cè)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 102249；3.中國(guó)石油大學(xué)（北京）提高采收率研究院，北京 102249；4.中海石油深圳分公司研究院，廣東 廣州 510240）

桑塔木地區(qū)三疊系深部地震薄儲(chǔ)集層預(yù)測(cè)

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

（1.中國(guó)地質(zhì)大學(xué)（北京）能源學(xué)院，北京 100083；2.中國(guó)石油大學(xué)（北京）油氣資源與探測(cè)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 102249；3.中國(guó)石油大學(xué)（北京）提高采收率研究院，北京 102249；4.中海石油深圳分公司研究院，廣東 廣州 510240）

深部隱蔽油氣藏逐步成為現(xiàn)階段油氣勘探的重難點(diǎn)，其具有埋深大、地震信號(hào)弱、構(gòu)造復(fù)雜等特征。埋藏深度的增大使得地震反射信息量減少，識(shí)別儲(chǔ)集層尤其是薄儲(chǔ)集層的難度加大。文中采用地震相控非線性隨機(jī)反演方法，融入測(cè)井和層序界面等地質(zhì)類(lèi)信息，對(duì)桑塔木油田三疊系近5 000 m埋深的薄儲(chǔ)集層進(jìn)行了鉆前預(yù)測(cè)，成功刻畫(huà)了深層薄儲(chǔ)集層的縱橫向連通性及平面展布特征。通過(guò)未參與反演井預(yù)測(cè)效果分析，少井反演與多井反演效果對(duì)比，預(yù)測(cè)吻合率統(tǒng)計(jì)分析和新鉆井與預(yù)測(cè)結(jié)果對(duì)比等方面來(lái)評(píng)價(jià)本次地震資料儲(chǔ)集層預(yù)測(cè)的可信度。預(yù)測(cè)結(jié)果與實(shí)鉆情況吻合程度較好，說(shuō)明該方法在深層薄儲(chǔ)集層預(yù)測(cè)中有很好的應(yīng)用前景。

深層；薄儲(chǔ)集層；相控反演；桑塔木油田

0　引言

隨著油氣勘探研究的不斷深入，研究重點(diǎn)已經(jīng)從尋找淺層、構(gòu)造油氣藏向深層、巖性油氣藏方向發(fā)展［1-3］。深部地層一般埋深超過(guò)4 000 m，地質(zhì)條件復(fù)雜，地震反射弱、主頻低、品質(zhì)差，因此，深部地震儲(chǔ)集層，特別是薄層預(yù)測(cè)有很大難度。

對(duì)于薄砂體儲(chǔ)層預(yù)測(cè)，基于20世紀(jì)80年代提出的Widess準(zhǔn)則［4］，能識(shí)別大于1/4地震波長(zhǎng)的薄儲(chǔ)層。此后許多研究發(fā)現(xiàn)，薄層描述可以突破傳統(tǒng)分辨率的極限［5-13］。

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

針對(duì)以上難題，本文采用地震相控非線性隨機(jī)反演方法，對(duì)塔里木盆地桑塔木地區(qū)三疊系地層進(jìn)行儲(chǔ)集層預(yù)測(cè)。

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

地震相控非線性隨機(jī)反演是在反演的過(guò)程中，融入地震相界面和沉積特征等地質(zhì)類(lèi)信息，在相模型的約束和指導(dǎo)下，結(jié)合鉆探得到的測(cè)井資料，采用非線性隨機(jī)反演算法，在計(jì)算過(guò)程中逐道進(jìn)行外推反演［14］。具體實(shí)現(xiàn)過(guò)程可以分為地震相控約束外推計(jì)算、隨機(jī)模擬處理和非線性反演3個(gè)部分。

1.1地震相控約束外推計(jì)算

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

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

由于地下地質(zhì)情況具有強(qiáng)不確定性，在相控外推計(jì)算過(guò)程中，可以運(yùn)用多項(xiàng)式相位時(shí)間擬合的方法來(lái)建立地震道間的外推關(guān)系。具體實(shí)現(xiàn)方法如下：在層序和相控模型控制的特定地震時(shí)窗范圍約束下，將從測(cè)井信息中或井旁地震道中獲取的先驗(yàn)?zāi)Ｐ蛥?shù)向量從井點(diǎn)處，開(kāi)始往外推計(jì)算。

1.2隨機(jī)模擬處理

考慮到地質(zhì)條件的隨機(jī)特性，首先利用測(cè)井巖石物理統(tǒng)計(jì)分析，建立速度、密度等模型參數(shù)點(diǎn)之間的統(tǒng)計(jì)關(guān)系；然后采用高斯隨機(jī)函數(shù)把模型參數(shù)作為區(qū)域化變量來(lái)進(jìn)行隨機(jī)模擬。

具體需要利用變差函數(shù)來(lái)描述空間數(shù)據(jù)場(chǎng)中數(shù)據(jù)之間的相互關(guān)系，從而建立起空間儲(chǔ)層參數(shù)點(diǎn)之間的統(tǒng)計(jì)相關(guān)函數(shù)關(guān)系式［16-17］。變差函數(shù)描述的是不同位置變量的不相似性。

1.3非線性反演

目前廣泛應(yīng)用于生產(chǎn)的線性反演算法，計(jì)算速度快，穩(wěn)定性強(qiáng)。然而，在實(shí)際情況中遇到的反演問(wèn)題基本都屬于非線性的，如果直接轉(zhuǎn)成線性問(wèn)題來(lái)求解，會(huì)導(dǎo)致最終的反演結(jié)果精度低，反映不出薄層特征，計(jì)算結(jié)果也容易陷入局部極小等問(wèn)題［18-19］。非線性反演算法直接求解非線性問(wèn)題，解的精度高，但穩(wěn)定性差。本文采用非線性最優(yōu)化算法，求解過(guò)程采用截?cái)嗥娈愔捣纸夥ú⒓尤爰s束條件來(lái)提高解的穩(wěn)定性與精度。

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

式中：m為速度；uobs為實(shí)際地震數(shù)據(jù)；umbd為模型響應(yīng)。

然后，將式（1）在初始模型響應(yīng)mi處按Taylor公式展開(kāi)，忽略二次以上的高階項(xiàng)得：

式中：D為f（m）的二次偏導(dǎo)數(shù)矩陣。

令m=mi+1，則有：

將式（3）對(duì)mi+1求一階導(dǎo)數(shù)，并令該導(dǎo)數(shù)為0，可得迭代公式為

在求取D-1時(shí)往往存在計(jì)算量大、矩陣病態(tài)等問(wèn)題，為了增強(qiáng)解的穩(wěn)定性，本文運(yùn)用了截?cái)嗥娈愔捣纸夥ǎ苊庵苯忧竽娲嬖?。同時(shí)在迭代過(guò)程中，通過(guò)加入最大最小值、參數(shù)包絡(luò)等約束，進(jìn)一步提高解的穩(wěn)定性和精度。通過(guò)式（4）迭代計(jì)算，可得到最終反演速度體。

2　實(shí)際應(yīng)用

2.1地質(zhì)背景

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

由于目標(biāo)薄儲(chǔ)集層中單砂體厚度很薄，只能在測(cè)井曲線上看到測(cè)井響應(yīng)，而在地震剖面上無(wú)法看到單砂體地震響應(yīng)。因此，本文嘗試采用地震相控非線性隨機(jī)反演方法對(duì)塔里木油田三疊系進(jìn)行儲(chǔ)集層預(yù)測(cè)，期望在反演剖面中能夠預(yù)測(cè)出多套薄層砂體，解決該區(qū)深層薄儲(chǔ)集層預(yù)測(cè)的難題。

2.2速度反演效果分析

2.2.1未參與反演井預(yù)測(cè)結(jié)果

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

2.2.2少井反演與多井反演效果對(duì)比

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

2.2.3預(yù)測(cè)吻合率統(tǒng)計(jì)分析

工區(qū)內(nèi)，選取未參與反演的50口井的綜合測(cè)井解釋及巖性錄井資料，與利用反演剖面預(yù)測(cè)的砂體厚度進(jìn)行對(duì)比統(tǒng)計(jì)。50口井預(yù)測(cè)結(jié)果中有6口井不符合實(shí)鉆情況，吻合率為88%。說(shuō)明地震預(yù)測(cè)與實(shí)際鉆井之間差異性較小，可將該反演成果用于儲(chǔ)集層分布預(yù)測(cè)研究。

2.2.4新鉆井驗(yàn)證結(jié)果分析

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

圖中：1）老井LN30井T3油組1小層井段為4 869～4 875 m的油層向LN302h井方向，速度降低，趨于泥巖速度?？梢悦黠@看到1小層連井處，反演剖面從暖色偏高值區(qū)域漸變?yōu)槔渖椭祬^(qū)域，反演剖面上含油砂體的厚度逐漸變薄直至尖滅，盡管砂泥巖聲波測(cè)井速度差異小，反演剖面上依然可以看到砂體尖滅處的細(xì)微差別，較好地解決了砂泥巖常規(guī)反演區(qū)分不清的問(wèn)題。2）LN30井2小層井段4 873.5～4 882.5 m橫向較穩(wěn)定。LN302h井完鉆結(jié)果顯示，在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的泥巖隔層，預(yù)測(cè)情況與實(shí)際鉆探結(jié)果吻合。新井預(yù)測(cè)結(jié)果表明，本文研究方法可用于新區(qū)滾動(dòng)勘探指導(dǎo)新井井位的部署。

2.3儲(chǔ)集層平面特征預(yù)測(cè)

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

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

綜上所述，相控反演的預(yù)測(cè)結(jié)果很好地突顯了儲(chǔ)集層的平面分布形態(tài)、厚度展布特征及沉積相特征，為本區(qū)深部巖性油氣藏的勘探部署提供了重要的依據(jù)。

3　結(jié)論

1）地震相控非線性隨機(jī)反演方法很好地融合了沉積層序等地質(zhì)類(lèi)的信息，反演的砂體預(yù)測(cè)結(jié)果更加符合區(qū)域沉積背景，解釋性和預(yù)測(cè)性更強(qiáng)。

2）速度反演的效果對(duì)比分析表明，該反演方法不過(guò)分依賴(lài)于參與反演的約束井?dāng)?shù)量和初始模型，約束井參與反演前后對(duì)比圖非常一致，且新鉆井的實(shí)鉆結(jié)果和預(yù)測(cè)結(jié)果吻合度很好，說(shuō)明反演方法的可靠程度高，該方法優(yōu)于常規(guī)反演方法，反演結(jié)果能為深部油氣藏的勘探部署提供可靠依據(jù)。

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

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

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

［3］李慶忠.巖性油氣藏地震勘探若干問(wèn)題討論（Ⅰ）［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］徐偉慕，郭平，胡天躍.薄互層調(diào)諧與分辨率分析［J］.石油地球物理勘探，2013，48（5）：750-757.

［11］楊飛，黃誠(chéng)，李鵬飛，等.深層砂泥巖薄互層的地震識(shí)別方法［J］.重慶大學(xué)學(xué)報(bào)，2014，37（9）：136-143.

［12］黃捍東，張如偉.地震非線性隨機(jī)反演方法在陸相薄砂巖儲(chǔ)層預(yù)測(cè)中的應(yīng)用［J］.石油學(xué)報(bào)，2009，30（3）：340-386.

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

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

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

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

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

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

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

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

［21］呂雪雁，朱筱敏，申銀民，等.塔里木盆地臺(tái)盆區(qū)三疊系層序地層研究和有利勘探區(qū)預(yù)測(cè)［J］.石油勘探與開(kāi)發(fā)，2002，29（1）：32-35.

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

［23］劉文嶺.地震約束儲(chǔ)層地質(zhì)建模技術(shù)［J］.石油學(xué)報(bào)，2008，29（1）：64-68.

［24］袁曉宇，李映濤，張哨楠，等.塔西南玉北地區(qū)鷹山組縫洞型儲(chǔ)層裂縫檢測(cè)技術(shù)及應(yīng)用［J］.石油實(shí)驗(yàn)地質(zhì)，2015，37（3）：394-401.

（編輯楊會(huì)朋）

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

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃（“973”計(jì)劃）項(xiàng)目“深部有效碎屑巖儲(chǔ)層成因機(jī)制與發(fā)育模式”（2011CB201104）

TE132.1+4；P631

A

10.6056/dkyqt201604009

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

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

引用格式：汪佳蓓，黃捍東，張文珠.桑塔木地區(qū)三疊系深部地震薄儲(chǔ)集層預(yù)測(cè)［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.