地震DNA地層超剝點線識別技術(shù)及應(yīng)用

羅紅梅

(1.中國石油化工股份有限公司勝利油田分公司勘探開發(fā)研究院,山東東營257015;2.中國石油化工股份有限公司勝利油田分公司博士后工作站,山東東營257015)

地震DNA地層超剝點線識別技術(shù)及應(yīng)用

羅紅梅1,2

(1.中國石油化工股份有限公司勝利油田分公司勘探開發(fā)研究院,山東東營257015;2.中國石油化工股份有限公司勝利油田分公司博士后工作站,山東東營257015)

摘要:地層油藏超剝點線的地球物理精細(xì)描述有利于提高含油氣儲層識別的可靠性和鉆井成功率,但在地震描述中根據(jù)地震反射特征追蹤的超剝點線與實際地層油藏的超剝點線往往存在較大的誤差。介紹了地震DNA地層超剝點線識別技術(shù)的方法原理與實現(xiàn)步驟,通過將地震反射數(shù)據(jù)和螞蟻體屬性數(shù)據(jù)轉(zhuǎn)換成相應(yīng)的字符數(shù)據(jù),采用地震DNA搜索方法將對應(yīng)字符所表示的地震數(shù)據(jù)體中的超剝點識別出來。在此基礎(chǔ)上,應(yīng)用聚類算法將超剝點進(jìn)行合理連線得到地層超剝線。模型合成地震記錄的超剝點識別結(jié)果表明,地震DNA地層超剝點識別技術(shù)能夠有效識別出超剝點,與傳統(tǒng)瞬時相位識別的超剝點相比,該方法的識別精度更高。東營凹陷草橋地區(qū)的實際應(yīng)用表明,該技術(shù)可以更加有效地描述超剝點線,實現(xiàn)了該區(qū)低序級不整合超剝邊界的精細(xì)刻畫,比基于瞬時相位屬性描述的超剝線的識別精度明顯提高。

關(guān)鍵詞：地震DNA;字符數(shù)據(jù);超剝點;聚類;超剝線

地層油藏是一種重要的隱蔽油氣藏類型,在世界石油和天然氣的產(chǎn)量和儲量中,43%的石油儲量和30%的天然氣儲量存在于地層圈閉中,但在地震描述中，根據(jù)地震反射特征追蹤的地層尖滅線與實際地層油藏的尖滅線常常存在較大誤差,地層油藏超剝線精細(xì)落實困難,導(dǎo)致鉆探成功率低,僅在新疆的克拉瑪依、準(zhǔn)噶爾、南襄盆地的泌陽和渤海灣盆地等少數(shù)地區(qū)實現(xiàn)了成功鉆探[1-4],如何準(zhǔn)確落實地層油藏的超剝線位置是地層油藏勘探開發(fā)的關(guān)鍵。目前,在地層油藏勘探開發(fā)中,超剝點和超剝線描述技術(shù)發(fā)展較快,其中,地震反射夾角外推技術(shù)是利用不整合與地層的傾角差外推地層尖滅點的預(yù)測技術(shù)[5-6],該技術(shù)工作量大,局限性嚴(yán)重,難以得到廣泛應(yīng)用。利用地震屬性進(jìn)行超剝點和超剝線描述的技術(shù)[7-8],是沿剖面的不整合面或超覆層提取多種地震屬性,通過對比分析確定用于地層圈閉識別的優(yōu)勢屬性,如平均瞬時相位、波峰波谷振幅屬性等,但是由于受地震資料分辨率限制,根據(jù)地震剖面和屬性識別出來的地層超剝點與實際地層的超剝點之間存在著一定的誤差。利用地質(zhì)統(tǒng)計方法進(jìn)行超剝點識別和描述的技術(shù)[9-10]，是對未剝蝕區(qū)的地層厚度數(shù)據(jù)擬合，然后外推確定超剝點位置，但該技術(shù)僅適用于被動大陸邊緣、大型三角洲沉積區(qū)及坳陷型盆地。構(gòu)造地質(zhì)學(xué)中制作復(fù)原橫剖面的地層厚度梯度法可以用于確定超剝點的位置[11],但其精度難以把握。周衛(wèi)紅等[12]、王章青等[13]、謝風(fēng)猛等[14]、胡宇雙等[15]應(yīng)用精細(xì)地層對比、層序地層學(xué)方法確定地層超剝點位置,但其計算結(jié)果往往因人而異。波形分析技術(shù)[16]是通過計算“曼哈頓距離”來識別尖滅點的準(zhǔn)確位置。所謂“曼哈頓距離”是指衡量地震波形與所選取參考子波非相似程度的參數(shù)，在尖滅點位置“曼哈頓距離”為最大值，其他位置則相對較小。地震DNA技術(shù)[17-19]能夠有效提取數(shù)據(jù)體中的構(gòu)造特點,為此,我們提出了地震DNA地層超剝點線識別技術(shù),采用類似于生物遺傳學(xué)基因組合識別的方式識別超剝點,并在此基礎(chǔ)上,利用聚類算法將超剝線提取出來,最后用理論和實際數(shù)據(jù)驗證本文提出的方法有效性。

1方法原理及實現(xiàn)

1.1地震DNA技術(shù)原理

在生物信息學(xué)中,經(jīng)常用到DNA搜索技術(shù)去檢測DNA分子中的堿基對序列,地震DNA技術(shù)就是受這種方法的啟發(fā)而產(chǎn)生的一項新技術(shù)。該技術(shù)已經(jīng)被成功應(yīng)用于地震層位和層序的自動提取[18-19]。這種搜索技術(shù)的核心類似于文本搜索中的正則表達(dá)式。

地震DNA的實現(xiàn)過程主要有兩步。第一步叫做“轉(zhuǎn)化”,也就是將輸入的數(shù)值轉(zhuǎn)化為文本。通過指定每個字符為一個唯一范圍的數(shù)據(jù)值將連續(xù)的數(shù)據(jù)轉(zhuǎn)換為字符。例如:

(1)

式中:字符a,b,c分別對應(yīng)于輸入數(shù)據(jù)范圍[min,-0.1),[-0.1,0.1),[0.1,max],min和max分別對應(yīng)于輸入數(shù)據(jù)的最小值和最大值。這樣就完成了數(shù)據(jù)從數(shù)值到字符的轉(zhuǎn)化過程,也就是將輸入數(shù)據(jù)從模型域轉(zhuǎn)化到DNA域(圖1)。第二步就是設(shè)置搜索模式,這個搜索模式也叫做“基因”。眾所周知,基因是由一系列“核苷酸”排序而成。目的是為了在龐大的DNA庫中搜索出感興趣的“基因”,實現(xiàn)方法是創(chuàng)建一個能夠描述待尋找特征的正則表達(dá)式。圖1中基因為a{nmin,nmax},b{nmin,nmax}和c{nmin,nmax},其中,字符a,b,c是字符轉(zhuǎn)換中的字符,花括號中的數(shù)字表示對應(yīng)字符重復(fù)次數(shù)的取值范圍,花括號的左邊表示最小重復(fù)次數(shù),右邊表示最大重復(fù)次數(shù)。例如一個簡單的字符基因組合的正則表達(dá)式(DNA鏈):a{2,3}b{1,2}c{4}b{0,1}a{3}。如果使用在字符轉(zhuǎn)換中描述的模式描述式(1),則地震數(shù)據(jù)范圍在[min,-0.1)中的轉(zhuǎn)換成字符a,范圍在[-0.1,0.1)中的轉(zhuǎn)換成字符b,范圍在[0.1,max]中的字符轉(zhuǎn)換成字符c。根據(jù)正則表達(dá)式(DNA鏈)的形式,字符a至少重復(fù)2次,至多重復(fù)3次,緊接著字符b至少重復(fù)1次,至多重復(fù)2次,然后字符c重復(fù)4次,接著字符b至少重復(fù)0次,至多重復(fù)1次,最后字符a至少重復(fù)3次,至多重復(fù)5次,即它們的組合有aabccccaaa,aaabccccaaa,aabbccccaaa,aaabbccccaaa,aabccccbaaa,aaabccccbaaa,aabbccccbaaa,aaabbccccbaaa這八種有效匹配方式。正則表達(dá)式搜索還允許字符為零的重復(fù),這意味零重復(fù)也是一個有效的匹配,使用零重復(fù)這個特性,可以創(chuàng)建更復(fù)雜的正則表達(dá)式,讓其中一些子特性在一個有效的匹配中出現(xiàn)或者不出現(xiàn)。

圖1　地震DNA流程結(jié)構(gòu)

地震DNA技術(shù)并不限定輸入數(shù)據(jù),其輸入數(shù)據(jù)可以為地震振幅數(shù)據(jù),也可以是任意一種地震屬性數(shù)據(jù)。由于地震DNA可以在不同屬性之中同時進(jìn)行搜索,這種多屬性的特征給此方法帶來了很大的靈活性。

1.2地層超剝點識別方法

1.2.1方法的提出

傳統(tǒng)的地震DNA搜索主要是應(yīng)用于地震層位和層序的自動提取。因為地震層位是能產(chǎn)生連續(xù)反射的界面,在地震剖面中表現(xiàn)為相同或者相似的波形、波組或者波系結(jié)構(gòu)。即傳統(tǒng)的地震DNA搜索主要是搜索地震數(shù)據(jù)中相似的基因。而地層超剝點的識別與此特性相反。如圖2所示,通過將地震反射振幅數(shù)據(jù)轉(zhuǎn)換成字符數(shù)據(jù),可以發(fā)現(xiàn)超剝點兩側(cè)明顯都具有不同的地震DNA基因序列,或表現(xiàn)為基因類型突變,或表現(xiàn)為某個字符重復(fù)次數(shù)突變?；诖?我們創(chuàng)新性地將地震DNA技術(shù)應(yīng)用于超剝點識別。

圖2　東營凹陷超剝點處的DNA突變示意

1.2.2識別方法及其實現(xiàn)步驟

將地震反射振幅數(shù)據(jù)、相位數(shù)據(jù)或螞蟻體屬性數(shù)據(jù)轉(zhuǎn)換成字符數(shù)據(jù),發(fā)現(xiàn)字符序列即DNA在超剝點處表現(xiàn)為突變,即超剝點兩側(cè)DNA(基因類型或單個字符的重復(fù)次數(shù))是不相同的,基于此特點,建立如圖3所示的地震DNA超剝點識別流程,主要由下述8個步驟組成:

1) 輸入三維地震數(shù)據(jù);

2) 按照Inline或者Crossline方向讀入二維地震剖面數(shù)據(jù);

3) 輸入不整合層位數(shù)據(jù);

4) 地震數(shù)據(jù)預(yù)處理:根據(jù)地震資料的品質(zhì)可選做,旨在提高信噪比和分辨率,如f-x域Cadzow濾波技術(shù)提高地震數(shù)據(jù)的信噪比[20];邊界保持各向異性擴(kuò)散濾波[21],使得超剝現(xiàn)象變得清楚;以及譜反演提高分辨率等[22];

5) 截取不整合面上下一定時窗內(nèi)的地震數(shù)據(jù)或?qū)傩詳?shù)據(jù),并設(shè)定字符轉(zhuǎn)換參數(shù)(字符轉(zhuǎn)換個數(shù)與范圍參數(shù)),將地震數(shù)據(jù)或?qū)傩詳?shù)據(jù)轉(zhuǎn)化為字符數(shù)據(jù);

6) 設(shè)定搜索參數(shù),包括相似度、容忍度和匹配度;

7) 搜索二維剖面的超剝點,記錄超剝點信息;

8) 搜索三維數(shù)據(jù)的超剝點。

其中,字符轉(zhuǎn)換個數(shù)決定圖像的分辨率,而范圍參數(shù)直接決定提取目標(biāo)同相軸的完整性和清晰性。字符轉(zhuǎn)換個數(shù)通常通過不同字符個數(shù)的對比試驗來優(yōu)選。對于范圍參數(shù),如果設(shè)定的值過小,追蹤出來的同相軸中包含的干擾信號就會比較多,無法有效區(qū)分目標(biāo)同相軸和干擾;如果設(shè)定的值過大,會導(dǎo)致目標(biāo)同相軸不完整或不清晰,甚至沒有目標(biāo)同相軸,通常通過對比試驗來優(yōu)選范圍參數(shù)。搜索參數(shù)中的相似度是對原始DNA算法的改進(jìn)之一,該參數(shù)在原始DNA算法中是一種字符串的匹配,也就是給定一個字符串,它能將與該字符串相似的字符串搜索出來，即搜索的是相似字符串。而對于尋找的超剝點就不一樣,超剝點兩側(cè)的字符串類型是不相似的,當(dāng)搜索到突變點時,搜索停止,搜索目的就是尋找突變點。因此,在超剝點識別中該參數(shù)為需要檢測的當(dāng)前層位上下不同字符的個數(shù),也就是表示需要檢測的同相軸個數(shù),可以通過在不整合層位之下字符串類型搭配的情況來尋找突變點。容忍度是對原始DNA算法的改進(jìn)之二,表示橫向上不滿足相似條件的道數(shù)超過該值時被認(rèn)為是超剝點,該參數(shù)的設(shè)定是為了減少誤差,因為在搜索過程中難免會存在誤差,假如字符串類型出現(xiàn)突變,而突變后的字符類型又與之前的一致,這時就得判斷該突變點是超剝點還是因誤差導(dǎo)致的。鑒于此,設(shè)定容忍度來有效控制誤判的發(fā)生。匹配度是對原始DNA算法的改進(jìn)之三,它表示的是相鄰道的字符最大重復(fù)次數(shù)的差值即代表同相軸厚度的變化范圍,一般在同相軸分叉的地方,同相軸的視周期較大,而分叉又是導(dǎo)致字符串突變的原因之一,在字符轉(zhuǎn)換之后,字符的最大重復(fù)次數(shù)與同相軸的視周期息息相關(guān),而在沿著不整合層位向下搜索同相軸的過程中,常常會出現(xiàn)字符串類型相同,但字符重復(fù)次數(shù)不同的情況,這就需要設(shè)定一個控制參數(shù)即匹配度來判斷突變點。

圖3　地震DNA超剝點識別流程

1.3地層超剝線識別方法

地層超剝線是由識別出來的超剝點連接而成,因此首先要對超剝點進(jìn)行合理分類,將屬于同一條超剝線的超剝點聚成一類,從而實現(xiàn)超剝線的提取。超剝點的聚類原則是將合理的超剝點分為一類,所謂合理,就是指超剝點之間的距離和方向滿足一定的最優(yōu)條件。由于超剝點之間的距離和方向?qū)τ诰垲惥哂屑s束意義,因此可以根據(jù)這些約束條件對超剝點進(jìn)行聚類。描述超剝點之間的接近程度用“距離”來度量。兩個超剝點之間的距離越小,表示兩者之間共同點越多,距離越大,共同點越小。在平面上,超剝點都是用線號和道號來表示,所以用歐式距離公式表示:

(2)

同時將兩個類中最近的兩個超剝點距離作為這兩個類的距離,即最短距離原則[23]。

聚類流程如圖4所示。包含5個步驟:

1) 輸入超剝點數(shù)據(jù);

2) 按照距離和方向的要求從第i個超剝點開始判斷,i小于最大線號;

3) 給定一定條件進(jìn)行聚類搜索,然后判斷兩超剝點間的距離是否滿足條件,判斷兩超剝點間的傾角方向是否滿足條件,否就回到步驟2),i=i+1;

4) 記錄滿足條件的超剝點信息;

5) 輸出超剝線。

圖4　超剝點聚類算法流程

2模型試算

圖5為模擬實際地層沉積結(jié)構(gòu)設(shè)計的地質(zhì)模型與相應(yīng)的地震響應(yīng),自下而上發(fā)育5套厚度逐漸減薄的砂體,砂巖速度為2600～2900m/s,泥巖速度為2400～2700m/s,頻率為30Hz。不整合面角度和剝蝕地層的角度見圖5a。圖5b中還給出了5套砂體超剝層試算結(jié)果,很明顯5套砂體實際的超剝點與地震響應(yīng)記錄上的超剝點存在著一定的誤差,而且不難看出,隨著砂體厚度的減薄,這種距離誤差隨之增大(表1)。

圖5　地質(zhì)模型(a)與相應(yīng)的地震響應(yīng)(b)

剝蝕地層厚度/m實際剝蝕點位置/m傳統(tǒng)瞬時相位識別結(jié)果基于地震DNA識別結(jié)果剝蝕點位置/m誤差/m剝蝕點位置/m誤差/m5037045080409.539.530730830100777.047.020113012801501218.088.020137015201501470.0100.010155019003501774.5224.5

針對上述地質(zhì)模型,將正演模擬得到的合成地震記錄作為地震反射振幅數(shù)據(jù)輸入,該數(shù)據(jù)未經(jīng)過任何地震預(yù)處理,采用本文方法進(jìn)行超剝點識別(圖6)。圖6中從上往下依次為地質(zhì)模型的地震響應(yīng)記錄、地震DNA超剝點識別結(jié)果和傳統(tǒng)瞬時相位屬性的識別結(jié)果。對比發(fā)現(xiàn),兩種方法識別的結(jié)果都比較準(zhǔn)確,每一個識別出來的點都有效對應(yīng)于一個砂體的超剝點。為了更加清楚地說明地震DNA識別算法的精度,對識別結(jié)果進(jìn)行了定量統(tǒng)計(表1),進(jìn)一步分析了超剝點識別誤差的大小和造成誤差的原因。

不難發(fā)現(xiàn),無論是傳統(tǒng)瞬時相位還是基于地震DNA識別的結(jié)果都存在著一定的誤差,這是由于在地震剖面上地層超剝點線附近砂層的厚度明顯減薄,同時受地震資料分辨率的限制,地震反射往往提前變?nèi)趸蛘呦А谋?中也能看出隨著剝蝕地層厚度的減小,傳統(tǒng)瞬時相位與基于地震DNA識別結(jié)果的誤差都是增大的。地層厚度越薄,地震資料的分辨率就會越低,識別誤差就會越大。盡管如此,地震DNA地層超剝點識別技術(shù)比傳統(tǒng)瞬時相位識別的精度要高。這是因為地震DNA方法通過字符的轉(zhuǎn)換參數(shù)控制可以有效增加地層超剝點線附近的字符,從而間接增強了地震反射。對比圖6中地震響應(yīng)和搜索結(jié)果可以看出地層超剝點線附近地震反射明顯增強。

圖6　超剝點識別結(jié)果對比

為了進(jìn)一步驗證方法的有效性,建立了更為復(fù)雜的地質(zhì)模型(圖7)進(jìn)行分析。該模型給出了兩個不整合層位和不整合層位下超剝點的實際位置,超剝點的實際位置是以一個寬度范圍給出的。其中,Hor1_1級不整合和Hor2_2級不整合面之下分別發(fā)育7套和5套單層厚度相當(dāng)?shù)纳绑w。對該模型的正演模擬地震記錄進(jìn)行地震DNA超剝點識別,圖8和圖9是分別對Hor1_1級不整合和Hor2_2級不整合層位下超剝點的識別結(jié)果?？梢钥闯?對于1級不整合層位和2級不整合層位之下的每一個砂體,其對應(yīng)的超剝點都被很好地識別出來,識別結(jié)果刻畫模型砂體比較準(zhǔn)確。

圖7　復(fù)雜地質(zhì)模型地震響應(yīng)

圖8　Hor1_1級不整合層位下超剝點的識別結(jié)果

圖9　Hor2_2級不整合層位下超剝點的識別結(jié)果

3實例分析

3.1地震數(shù)據(jù)的預(yù)處理

選取東營凹陷草橋南部三維工區(qū)的實際資料進(jìn)行了地層超剝點、線的識別與描述。在地震振幅數(shù)據(jù)體中,存在斷層、超剝結(jié)構(gòu)等等。斷層處基因會發(fā)生突變,因此,在斷裂系統(tǒng)比較發(fā)育的地區(qū)進(jìn)行地震DNA超剝點線識別之前應(yīng)首先進(jìn)行斷層的提取工作,從而消除斷層的影響,目前斷層的識別方法有很多[24],如相干體技術(shù)[25]、波形分析技術(shù)[26]等,而現(xiàn)今最常用的自動斷層提取方法是螞蟻追蹤算法[27]?？紤]到地震反射同相軸與斷層數(shù)據(jù)體中斷層的相似性,可以利用螞蟻追蹤算法直接對地震振幅數(shù)據(jù)體進(jìn)行地震同相軸追蹤,進(jìn)而提高地震同相軸對尖滅點的分辨能力。

在進(jìn)行超剝點線識別之前,首先對地震數(shù)據(jù)進(jìn)行螞蟻體追蹤預(yù)處理。圖10a為草橋地區(qū)東西向地震振幅數(shù)據(jù)剖面,圖10b為螞蟻追蹤預(yù)處理的螞蟻體屬性剖面,對比可以看出,螞蟻體屬性剖面上同相軸連續(xù)性明顯更高,這充分說明螞蟻追蹤算法對于提高地層連續(xù)性是有效的。以此數(shù)據(jù)體為基礎(chǔ),進(jìn)行地震DNA分析,可以避免斷層的干擾,更好地識別出真實的超剝點線。

圖10　二維地震剖面(a)及螞蟻體屬性剖面(b)

3.2草橋地區(qū)超剝點線識別

首先利用常規(guī)瞬時相位方法對該地區(qū)地層超剝線進(jìn)行識別,再應(yīng)用本文提出的地震DNA地層超剝點線識別技術(shù)進(jìn)行識別。比較兩種方法的識別結(jié)果,以驗證本文方法的正確性和有效性。

從草橋地區(qū)過4口已鉆井的東西向地震剖面來看(圖11),在沙二段底部發(fā)育了一個二級不整合T3(黃色標(biāo)明的位置),為一個沉積不整合面,由于其角度較小,并且上下波阻抗差異小,地震上表現(xiàn)為連續(xù)性較差的中等振幅反射。從利用瞬時相位識別的二級不整合T3之下剝蝕地層邊界結(jié)果來看(圖12a),可以大致識別出幾條屬性相近可以連起來的條帶,圖中畫出的趨勢線就認(rèn)為是地層超剝線的展布。但識別結(jié)果并不太清楚,無法確切、有效地刻畫超剝線的具體位置。進(jìn)一步地,先輸入地震螞蟻體數(shù)據(jù),然后利用地震DNA地層超剝點線識別技術(shù)來識別二級不整合T3之下剝蝕地層的邊界。首先識別出地層超剝點,進(jìn)而利用超剝點聚類算法得到地層超剝線(圖12b),圖中展示的近東西向粗紅線為斷層,其兩側(cè)基因類型為ab-ba或ba-ab。其它顏色的細(xì)線為地層超剝線,刻畫完整,展布特征明顯。分析典型地震剖面的地震DNA識別結(jié)果(圖13)可以看到,每個超剝點左右兩側(cè)的DNA鏈(基因類型和單個字符的重復(fù)次數(shù))與斷層兩側(cè)基因差異較大。因此,地震DNA地層超剝點線識別結(jié)果清晰,每一個識別出來的超剝點均有效對應(yīng)一條超剝線,說明本文所提出方法應(yīng)用于地層超剝點、超剝線的識別非常有效。

圖11　過4口井近東西向地震剖面

圖12　草橋地區(qū)地層超剝線瞬時相位識別(a)與地震DNA識別(b)結(jié)果對比

圖13　草橋地區(qū)H線地震DNA識別結(jié)果

4結(jié)論

本文提出了一種地震DNA地層超剝點線識別與描述方法,該方法將振幅和螞蟻體屬性數(shù)據(jù)轉(zhuǎn)換為字符數(shù)據(jù),通過一個正則表達(dá)式的模式識別來進(jìn)行字符數(shù)據(jù)的搜索以獲得超剝點,并通過聚類提取超剝線。通過模型數(shù)據(jù)和實際資料的處理,以及與傳統(tǒng)的瞬時相位識別方法對比可以看出,該方法可以提高地層超剝點和超剝線描述的精度,尤其對于剝蝕地層厚度小于20m的情況下,識別精度提高明顯。

與傳統(tǒng)的瞬時相位識別超剝點相比,地震DNA方法具有更強的靈活性。地震DNA方法不限定輸入數(shù)據(jù),也就是說可以引入任意一種對識別有幫助的屬性進(jìn)行處理,甚至可以多屬性同時識別。在本文研究中,僅僅使用振幅數(shù)據(jù)和螞蟻體數(shù)據(jù)進(jìn)行處理,下一步可以考慮加入更多屬性同步處理,比如加入相干體和曲率屬性等。同時,地震DNA作為一種靈活性很強的方法技術(shù),在地震資料解釋工作中也具有廣泛的應(yīng)用前景,下一步可以考慮將其應(yīng)用于某些特定地質(zhì)體的自動識別。

參考文獻(xiàn)

[1]周小進(jìn),陳躍,黃繼文,等.庫車坳陷南緣地層超覆圈閉特征及含油氣前景[J].石油實驗地質(zhì),2000,22(1):40-45

ZHOU X J,CHEN Y,HUANG J W,et al.Characteristics of stratigraphic overlap traps and petroleum prospect in the south margin of Kuqa depression[J].Experimental Petroleum Geology,2000,22(1):40-45

[2]陳宗清.川東石炭系地層—構(gòu)造復(fù)合圈閉氣藏的預(yù)測、分類和評價[J].石油學(xué)報,2000,21(3):25-30

CHEN Z Q.Forecast,classification and assesment for gas reservoirs of carboniferous stratigraphic-structural combination traps in eastern Sichuan[J].Acta Petrolei Sinica,2000,21(3):25-30

[3]李丕龍.陸相斷陷盆地油氣地質(zhì)與勘探(卷五)[M].北京:石油工業(yè)出版社,2003:31-39

LI P L.Petroleum geology and exploration of continental fault basin (book five)[M].Beijing:Petroleum Industry Press,2003:31-39

[4]隋風(fēng)貴,趙樂強,林會喜,等.濟(jì)陽坳陷第三系地層油藏形成機(jī)制與組合模式[J].油氣地質(zhì)與采收率,2005,12(6):1-5

SUI F G,ZHAO L Q,LIN H X,et al.The forming mechanism and its combination mode of stratigraphic reservoir in Tertiary in Jiyang Depression[J].Petroleum Geology and Recovery Efficiency,2005,12(6):1-5

[5]蘇朝光,閆昭岷,張營革,等.地層油藏超剝尖滅線夾角定量外推方法模型研究[J].地球物理學(xué)進(jìn)展,2007,22(6):1841-1846

SU C G,YAN Z M,ZHANG Y G,et al.The study of quantitive extrapolation model for overlap-denudation pitching-out line stratigraphic reservoir[J].Progress in Geophysics,2007,22(6):1841-1846

[6]張福利,趙麗平.定量確定地層超覆線位置的地震反射夾角外推法[J].勘探地球物理進(jìn)展,2008,31(4):278-280

ZHANG F L,ZHAO L P.Quantitative determination of overlap line position via extrapolation of angle of seismic reflection[J].Progress in Exploration Geophysics,2008,31(4):278-280

[7]管曉燕,穆星,劉軍.相位分析技術(shù)在地層圈閉描述中的應(yīng)用[J].油氣地質(zhì)與采收率,2007,14(6):56-58

GUAN X Y,MU X,LIU J.Application of the phase analysis technique for stratigraphic trap description[J].Petroleum Geology and Recovery Efficiency,2007,14(6):56-58

[8]魯紅利.哈山南緣淺層超剝帶圈閉特征研究[J].內(nèi)江科技,2012,33(12):150-146

LU H L.Research of pitching-out belt trap in the south shallow of Hasan[J].Neijiang Technology,2012,33(12):150-146

[9]鐘廣法,周祖翼,耿建華.有效厚度擬合定量恢復(fù)原始地層超覆點位置[J].沉積學(xué)報,2003,21(4):614-618

ZHONG G F,ZHOU Z Y,GENG J H.Extrapolating the original positions of removed onlap-offlap points using a piecewise-constructed model of effective stratal thickness[J].Acta Sedimentological Sinica,2003,21(4):614-618

[10]張營革,田建華.沾車地區(qū)館陶組地層超覆油藏地震描述技術(shù)應(yīng)用及效果分析[J].石油物探,2003,42(3):322-327

ZHANG Y G,TIAN J H.Seismic description technology for stratigraphic overlap reservoir of the Guantao Formation and its effect analysis in Zhanche area[J].Geophysical Prospecting for Petroleum,2003,42(3):322-327

[11]江濤,苗洪波,王芳,等.松遼盆地西部盆緣帶地層超覆油藏形成條件[J].大慶石油地質(zhì)與開發(fā),2006,25(4):24-26

JIANG T,MIAO H B,WANG F,et al.Generation condition for overlap reservoir in border belt of western Songliao basin[J].Petroleum Geology and Oilfield Development in Daqing,2006,25(4):24-26

[12]周衛(wèi)紅,謝風(fēng)猛,劉健.東營凹陷西北部地層超覆油藏勘探技術(shù)研究[J].承德石油高等?？茖W(xué)校學(xué)報,2005,7(3):23-26

ZHOU W H,XIE F M,LIU J.Exploration technology research on stratigraphic overlapped oil reservoirs in northwest of Dongying depression[J].Journal of Chengde Petroleum College,2005,7(3):23-26

[13]王章青,楊立戰(zhàn),林勝春.地震綜合解釋技術(shù)在地層超覆圈閉研究中的應(yīng)用[J].江漢石油學(xué)院學(xué)報,2003,25(2):60-61

WANG Z Q,YANG L Z,LIN S C.Application of integrated seismic interpretation technique in studying stratigraphic overlap trap[J].Journal of Jianghan Petroleum Institute,2003,25(2):60-61

[14]謝風(fēng)猛,李彥強.濱南油區(qū)地層超覆油藏的勘探[J].油氣地質(zhì)與采收率,2004,11(6):27-30

XIE F M,LI Y Q.Exploration of the stratigraphic overlap oil reservoirs in Binnan petroliferous province[J].Petroleum Geology and Recovery Efficiency,2004,11(6):27-30

[15]胡玉雙,李黎明.層序地層學(xué)在松遼東部斷陷有利勘探區(qū)帶預(yù)測中的應(yīng)用[J].大慶石油地質(zhì)與開發(fā),1999,18(2):10-12,26

HU Y S,LI L L.Application of sequence stratigraphy to the prediction of favorable exploration areas in the east rift of Songliao Basin[J].Petroleum Geology and Oilfield Development in Daqing,1999,18(2):10-12,26

[16]李海亮,陳啟林,朱允輝,等.波形分析技術(shù)在識別小斷層和剝蝕線中的應(yīng)用[J].新疆石油地質(zhì),2007,28(2):232-234

LI H L,CHEN Q L,ZHU Y H,et al.Application of wave analysis to identifying minor fault and nip-out line[J].Xinjiang Petroleum Geology,2007,28(2):232-234

[17]BAKKE J,GRAMSTAD O,S?NNELAND L.Seismic DNA——a novel visually guided search method using non-local search and multi attribute data sets[J].Expanded Abstract of 82ndAnnual Internat SEG Mtg,2012:589-593

[18]GRAMSTAD O,BAKKE J,S?NNELAND L,et al.Simultaneous extraction of stratigraphic sequences using iterative seismic DNA detection[J].Interpretation,2014,2(4):T167-T176

[19]GRAMSTAD O,BAKKE J,S?NNELAND L.Seismic surface extraction using iterative seismic DNA detection[J].Expanded Abstract of 82ndAnnual Internat SEG Mtg,2012:600-604

[20]崔樹果,朱凌燕,王建花.f-x域Cadzow技術(shù)分塊壓制隨機(jī)噪聲及其應(yīng)用[J].石油物探,2012,51(1):43-50

CUI S G,ZHU L Y,WANG J H.Random noise attenuation with Cadzow technique inf-xdomain and its application[J].Geophysical Prospecting of Petroleum,2012,51(1):43-50

[21]嚴(yán)哲,顧漢明,蔡成國.基于各向異性擴(kuò)散濾波的地震圖像增強處理[J].石油地球物理勘探,2013,48

(3):390-394

YAN Z,GU H M,CAI C G.Seismic image enhancement based on anisotropic diffusion[J].Oil Geophysical Prospecting,2013,48(3):390-394

[22]柴新濤,李振春,韓文功,等.基于LSQR算法的譜反演方法研究[J].石油物探,2012,51(1):11 -18

CHAI X T,LI Z C,HAN W G,et al.Spectral inversion method analysis based on the LSQR algorithm[J].Geophysical Prospecting for Petroleum,2012,51(1):11-18

[23]向東進(jìn),李宏偉,劉小雅.實用多元統(tǒng)計分析[M].武漢:中國地質(zhì)大學(xué)出版社,2005:1-184

XIANG D J,LI H W,LIU X Y.Applied multivariate statistical analysis[M].Wuhan:China University of Geosciences Press,2005:1-184

[24]王鵬,鐘建華,張亞金,等.斷層識別技術(shù)及其在肇州油田的應(yīng)用[J].石油物探,2011,50(5):521-525

WANG P,ZHONG J H,ZHANG Y J,et al.Fault identification technique and its application in Zhaozhou oilfield[J].Geophysical Prospecting for Petroleum,2011,50(5):521-525

[25]孫夕平,楊國權(quán).三維地震相干體技術(shù)在目標(biāo)區(qū)沉積相研究中的應(yīng)用[J].石油物探,2004,43(6):591-594

SUN X P,YANG G Q.Application of 3D seismic coherence cube technology to target area sedimentary facies research[J].Geophysical Prospecting for Petroleum,2004,43(6):591-594

[26]師永民,祁軍,張成學(xué),等.應(yīng)用地震波形分析技術(shù)預(yù)測裂縫的方法探討[J].石油物探,2005,44(2):128-130

SHI Y M,QI J,ZHANG C X,et al.Discussion of the method of applying seismic waveform analysis techniques to predict fracture[J].Geophysical Prospecting for Petroleum,2005,44(2):128-130

[27]嚴(yán)哲.三維地震斷層自動識別與智能解釋[D].武漢:中國地質(zhì)大學(xué),2010

YAN Z.Fault automatic recognition and intelligent interpretation in 3-D seismic data[D].Wuhan:China University of Geosciences,2010

(編輯：朱文杰)

Stratigraphic pitching-out point and line identification technology based on seismic DNA and its application

LUO Hongmei1,2

(1.ResearchInstituteofPetroleumExplorationandDevelopment,SinopecShengliOilfieldCompany,Dongying257015,China;2.PostdoctoralWorkstation,ShengliOilfield,SINOPEC,Dongying257015,China)

Abstract:The geophysical description of pitching-out point and line in stratigraphic reservoir is beneficial to improve the reliability of oil and gas reservoir recognition and success rate of well drilling.The pitching-out point and line is traced by seismic reflection characteristics.There is a large difference between the pitching-out point and line and the actual point and line.In this paper we introduce the method principle and implementation steps of stratigraphic pitching-out point and line identification technology based on seismic DNA.The seismic reflection data and ant-attribute data is converted into the corresponding character data,and then the pitching-out point from the character data indicated corresponding seismic data is identified using seismic DNA.On the basis,stratigraphic pitching-out points are reasonably linked to be the pitching-out line by using clustering.Model simulated records show that stratigraphic pitching-out point and line identification technology based on seismic DNA can effectively identify pitching-out point,and can improve the identification precision compared with the result of conventional instantaneous phase.The application result of actual data in Caoqiao area,Dongying depression shows that this technology can effectively describe pitching-out line and the fine description of the unconformity pitching-out boundary for thin layer is realized.

Keywords:seismic DNA,character data,pitching-out point,cluster,pitching-out line

收稿日期：2015-10-22;改回日期：2016-01-16。

作者簡介：羅紅梅(1973—),女,高級工程師,博士,主要從事地震地質(zhì)數(shù)據(jù)綜合處理與解釋方面的研究。

基金項目：中國石油化工股份有限公司科技開發(fā)部“東營凹陷南部地層超剝帶精細(xì)描述與目標(biāo)評價”項目(P13039)資助。

中圖分類號：P631

文獻(xiàn)標(biāo)識碼：A

文章編號：1000-1441(2016)03-0414-11

DOI:10.3969/j.issn.1000-1441.2016.03.012

This project is financially supported by the Project of SINOPEC Development Planning Department (Grant No.P13039).