SP測井小波變換在單井層序劃分中的應用

賀 聰，吳遠東，吉利明，張明震，梁曉飛

(1.甘肅省油氣資源研究重點實驗室/中國科學院油氣資源研究重點實驗室，甘肅 蘭州 730000；2.中國科學院大學，北京 100049；3.中國科學院地質與地球物理研究所，北京 100029)

?

SP測井小波變換在單井層序劃分中的應用

賀 聰1，2，3，吳遠東1，2，吉利明1，3，張明震1，梁曉飛1，2

(1.甘肅省油氣資源研究重點實驗室/中國科學院油氣資源研究重點實驗室，甘肅 蘭州 730000；2.中國科學院大學，北京 100049；3.中國科學院地質與地球物理研究所，北京 100029)

利用單井自然電位(SP)測井曲線經coif5小波變換，結合露頭剖面傳統層序地層劃分標準，探討了小波變換曲線與鄂爾多斯盆地南部A井三疊系陸相延長組三級層序界面的對應關系。結果顯示，SP測井曲線經coif5小波變換后的d12小波(SP_coif5_d12)極大值可能對應湖退—湖侵旋回轉換面(即層序界面)；而小波曲線的極小值對應最大湖泛面，SP_coif5_d12小波可以作為湖盆的三級層序劃分指標，說明coif5小波變換極值點位置檢測到的SP測井信號與層序界面具有對應關系。該方法非常簡捷，在陸相湖盆層序地層研究中具有實際應用價值，尤其適于缺乏巖心和地震等資料的單井研究。

自然電位測井；小波變換；層序地層；鄂爾多斯盆地；延長組

0 引 言

層序地層學理論及其技術方法在陸相盆地層序分析和油氣勘探開發中的適用性和有效性已經被學界廣泛接受，其關鍵是識別和劃分不同成因的界面以及不同級次的基準面旋回[1-2]。近年來，小波分析被引入到層序地層的研究中，利用測井信號小波變換被認為是實現定量劃分層序的有效途徑[3-7]，該技術的焦點是小波變換得到的信號與各級次分辨率層序的對應關系。利用露頭、巖心、測井和地震等資料，根據旋回沉積學和層序地層學理論劃分層序的方法已經很成熟，其結果通常也更易被人接受。然而獲取全井巖心和地震資料的成本很高，對于某些缺乏露頭、巖心和地震資料的井而言，傳統方法顯然有一定的缺陷。相比而言，利用測井信號進行小波變換劃分層序地層則具有簡單、快速和低成本的優勢。利用鄂爾多斯盆地南部A井的SP曲線進行coif5小波變換，探討其與三級層序地層界面的對應關系，并將研究區附近露頭剖面具有高可信度的傳統層序地層劃分方案[8]作為對比依據，來判斷和校驗該方法和指標的可靠性。

1 地質背景

1.1 研究區地質概況

鄂爾多斯盆地是中國中部非常重要的含油氣盆地，是華北克拉通最穩定的一部分，在中生代之前屬于華北古陸伸向秦祁海域的臺地邊緣區。晚三疊世開始，由于華北地臺西緣受到自西向東擠壓而形成坳陷型盆地，白堊紀晉西地區隆起后形成獨立的盆地[9]，現今被分為6個一級地質構造單元。印支期華北克拉通逐漸萎縮，進入鄂爾多斯內陸盆地演化階段，形成了一套東北部較淺、西南部較深的，總厚度超過1 000 m的河湖相—三角洲相碎屑沉積[10]，是盆地內重要的油氣產層。延長組在鄂爾多斯盆地四周都有物源區，在盆地中形成明顯的由河流沉積、三角洲沉積、半深湖沉積所組成的環狀相帶，使延長組經歷了內陸湖泊產生、發展乃至消亡的完整過程。利用研究區測井資料，根據凝灰巖和泥頁巖標志層自下而上將延長組劃分為長10—長1共10個油層組，其中盆地南部有些地區延長組頂部為長4+5，局部地區僅存長3部分地層[9]。

1.2 延長組層序地層特征

在鄂爾多斯盆地延長組的層序地層研究工作中，應用旋回沉積學和層序地層學理論，根據露頭、巖心、測井和地震等資料已獲得了大量成果[8，11-13]。延長組頂部和底部都是區域不整合面，屬于一個二級地層層序，代表了28.4 Ma，記錄了晚三疊世鄂爾多斯湖盆整個演化階段[8]。Zou等[8]根據金鎖關露頭資料并結合盆地其他地方的多個露頭剖面識別出延長組發育4個大尺度湖侵—湖退旋回(圖1)。延長組這4個湖侵—湖退地層旋回分別對應4個三級層序(SQ1、SQ2、SQ3和SQ4)，每個層序的時間跨度平均為7 Ma，最大湖泛面主要位于長9、長7、長4+5和長1層段中(圖1)。

圖1 金鎖關剖面延長組綜合柱狀圖(據Zou et al.， 2010，修改)

延長組頂部和底部層序界面是區域性不整合面，內部的層序界面主要表現為湖退—湖侵旋回轉換面，通常是大型河道底部下切面[8]。

2 小波變換及其劃分層序的方法

小波變換是一種工具，把數據、函數或運算符分割成不同頻率的成分，然后再用分解的方法(通過濾波器實現)去研究對應尺度下的各個成分。關于小波變換的原理可以參考文獻[14]。信號中的突變部分和奇異點等不規則部分通常包含重要的信息，應用小波變換可以檢測出信號中突變點的位置、類型以及變化的幅度，如將信號的突變或瞬變過渡得比較陡峭或者平穩，在小波多尺度變化上就表現為最大(小)值的情況[15]。常見的小波主要有haar小波、db小波、bior小波、coif小波、sym小波、morlet小波以及mexican-hat小波等，其中只有前5種小波可以用于離散小波變換，而haar小波在實際的信號處理與應用中具有較大限制，bior小波主要用于信號與圖像的重構，coif小波具有比db小波及其改進小波(sym小波)更好的對稱性，因此，選取coif小波開展相關研究。coif小波表示為coifN(N=1、2、3、4、5)，N越大處理后的信號去噪效果和平滑度越好，故選用coif5小波，其小波函數、尺度函數和分解濾波器信息見圖2。

圖2 coif5小波函數、尺度函數及分解濾波器圖

小波變換具有多分辨率的特點，能體現時域、頻域的局部特征，可以模擬由粗到細的層序單元人工劃分方法[5]。測井信號提供的關于地層的各種巖性物理參數具有縱向分辨率高的特點，是一定時間序列內各種沉積事件的物質記錄，能夠真實、連續地反映所測地層的沉積特征，對研究地層多級別旋回性及識別地層信息多分辨率突變具有優勢[7]。測井信號經小波變換處理后可以獲得清晰的頻率結構，探測到各頻率段之間能反映地質環境發生突變的區域，將多個不同周期(尺度)沉積旋回疊加的測井曲線分解成各自周期獨立的沉積旋回，以尺度的形式展示出來，形成特定的界面響應，根據多尺度因子下周期性震蕩特征與各級層序界面的對應關系劃分層序[6]。但是目前尚未見報道測井曲線小波變換劃分層序地層的評價標準，根據傳統層序地層劃分方案來鑒定其結果的可靠程度仍然是最為有效的指標。

3 結果與討論

3.1 測井曲線小波變換

在以往的小波變換劃分層序地層中，一般采用GR測井曲線。但研究區位于鄂爾多斯盆地南部，離火山口很近，GR值受濁流沉積、凝灰質火山巖等干擾嚴重，因此，GR測井不適宜進行該地區層序地層的劃分。而SP測井是另外一個最常用來劃分層序地層的測井曲線，因此選取SP曲線作為劃分層序地層的原始信號。利用Matlab R2013a軟件自帶的小波分析工具箱對A井延長組的SP測井信號基于coif5小波多尺度離散變換結果如圖3。其中S代表自然電位測井原始信號(SP，單位為mV)，a代表低頻信號，d代表高頻信號。原始信號S第1次分解獲得低頻信號a1和高頻信號d1，第2次將低頻信號a1分解為低頻信號a2和高頻信號d2，依次將信號進行12次分解(圖3a)，最后獲得高頻信號d1～d12和低頻信號a12，這13個信號之和即為原始信號S(圖3b)。完全分解模式圖中的橫坐標從左到右代表由淺到深方向，坐標值為自然電位測井采樣點序數，采樣點總數為9 200，對應該井段總深度460 m(圖3b)。

3.2 層序界面劃分與討論

根據前人研究的經驗[3-4，6-7]，通常利用測井曲線小波變換所得到的高頻信號的低頻部分信號與各級次的高分辨率層序地層相對應。將SP曲線coif5小波分解后的各級信號(圖3)分別與研究區金鎖關剖面的三級層序劃分方案(圖1)[8]對比，主要比較小波變換曲線的極大(小)值位置和三級層序界面，結果發現第12次分解獲得的低頻部分信號d12與該方案較吻合(圖4)。金鎖關露頭剖面距離A井僅10 km左右(圖1)，兩者沉積環境和層序發育一致，因此，具有可比性。從圖4中可以看出，SP小波變換曲線(SP_coif5_d12)的2個最大值處對應2個湖退—湖侵旋回轉換面(內部層序界面)，分別是長82頂部和長62頂部；而3個最小值處對應3個最大洪泛面，分別是長9下部、長7下部和長4+5頂部。在長10—長4+5地層中，長10—長82和長81—長62是2個完整的湖侵—湖退旋回，對應2個完整的三級層序，而長61—長4+5是1個湖侵半旋回，對應半個三級層序，這個劃分結果與金鎖關露頭剖面層序劃分結果完全一致。除此之外，根據圖4中長7中下部GR曲線異常高值段可以判斷其為熱頁巖層段，即延長組區域性標志層張家灘油頁巖，是已經公認的最大湖泛面(凝縮層段)發育位置，與小波變換后識別出的長7段最大湖泛面位置也很吻合。由此可見，可以將SP曲線經coif5小波變換后的d12小波(SP_coif5_d12)作為研究區三級層序劃分標準，其中SP_coif5_d12小波的極大值處對應三級層序的層序界面(湖退—湖侵旋回轉換面)，而極小值處對應最大湖泛面。

圖3 SP曲線coif5小波變換

圖4 A井延長組(長10—長4+5)層序地層劃分[8]

SP測井小波變換后的曲線極大(小)值處與層序界面(最大湖泛面)的對應關系，可能與小波變換識別信號的奇異點(突變點)作用密切相關。在信號檢測中，突變點的位置可以是小波變換的過零點，也可以是極值點，但一般來說，根據極值點做檢測比根據過零點做檢測更準確[16]。根據此次研究的結果，層序界面突變點的位置明顯是與極值點對應。SP_coif5_d12小波與研究區三級層序界面表現出一致性，說明對應小波變換識別出的信號奇異點可能就是三級層序界面，因此，可以將其作為湖盆三級層序的劃分指標。

4 結 論

(1)SP測井曲線經coif5小波變換劃分鄂爾多斯陸相湖盆南部延長組層序地層具有可行性，SP_coif5_d12小波的極大值處對應三級層序的層序界面(湖退—湖侵旋回轉換面)，極小值處對應最大湖泛面，說明coif5小波變換極值點位置檢測到的自然電位測井信號突變點與層序界面有對應關系。

(2)SP測井曲線經coif5小波變換后，頻繁突變的信號曲線變得平滑，突出顯示突變強度最大的信息，使得層序界面的識別變得簡便、快速，在陸相湖盆層序地層劃分中具有實用價值，尤其是在巖心和地震等資料缺乏的單井應用中更是可以降低研究成本。

[1] 鄭榮才，彭軍，吳朝容. 陸相盆地基準面旋回的級次劃分和研究意義[J]. 沉積學報， 2001， 19(2)： 249-255.

[2] 李景哲，張金亮，李存磊，等.層序地層學在應用中的兩點誤區[J].特種油氣藏， 2013， 20(6)： 52-55，143.

[3] 鄭希民，郭彥如，劉化清，等. 應用自然伽馬測井曲線小波分析劃分陸相坳陷盆地三級層序的方法——以鄂爾多斯盆地延長組為例[J]. 天然氣地球科學， 2006， 17(5)： 672-676.

[4] 李相博，郭彥如，劉化清，等. 淺談小波分析在鄂爾多斯盆地延長組層序地層劃分中的應用[J]. 天然氣地球科學， 2006， 17(6)： 779-782.

[5] 朱劍兵，紀友亮，趙培坤，等. 小波變換在層序地層單元自動劃分中的應用[J]. 石油勘探與開發， 2005， 32(1)： 84-86.

[6] 王志坤，王多云，宋廣壽，等. 測井信號小波分析在高分辨率層序地層劃分中的應用[J]. 大慶石油學院學報， 2005， 29(6)： 17-20，120.

[7] 余繼峰，李增學. 測井數據小波變換及其地質意義[J]. 中國礦業大學學報， 2003， 32(3)： 127-130.

[8] Zou C N，Zhang X Y，Luo P， et al. Shallow-lacustrine sand-rich deltaic depositional cycles and sequence stratigraphy of the upper triassic yanchang formation， Ordos Basin， China[J]. Basin Research， 2010， 22(1)： 108-125.

[9] 楊俊杰. 鄂爾多斯盆地構造演化與油氣分布規律[M]. 北京： 石油工業出版社， 2002： 200-228.

[10] 楊華，鄧秀芹. 構造事件對鄂爾多斯盆地延長組深水砂巖沉積的影響[J]. 石油勘探與開發， 2013， 40(5)： 513-520.

[11] Yang M H，Liu C Y，Zheng M L， et al. Sequence framework of two different kinds of margins and their response to tectonic activity during the middle-late triassic， Ordos Basin[J]. Science in China Series D-earth Sciences， 2007， 50(S2)： 203-216.

[12] Zhao J X，Shen Z J，Li L， et al. Sequence filling characteristics and distribution laws of large continental depressed lake Basin： a case study of yanchang formation， Ordos Basin[J]. Acta Petrologica Sinica， 2011， 27(8)： 2318-2326.

[13] 楊華，陳洪德，付金華. 鄂爾多斯盆地晚三疊世沉積地質與油藏分布規律[M]. 北京： 科學出版社， 2012： 19-56.

[14] 馬朋善，高秀花，汪桂春，等. 小波分析在地震資料分頻處理中的應用[J]. 特種油氣藏， 2007， 14(3)： 35-38，106.

[15] 周偉. MATLAB小波分析高級技術[M]. 西安： 西安電子科技大學出版社， 2005： 125.

[16] 劉偉，曹思遠. 基于小波變換的信號奇異性檢測在層位識別中的應用[J]. 石油地球物理勘探， 2010， 45(4)： 467，530-533，624.

編輯 黃華彪

20150807；改回日期：20150928

中國科學院戰略性先導科技專項(B類)“頁巖氣勘探開發基礎理論與關鍵技術”(XDB10010103)；國家自然科學基金“鄂爾多斯盆地三疊系延長組湖相烴源巖藻類母質及其屬性研究”(41172131)；甘肅省重點實驗室專項“鄂爾多斯盆地延長組烴源巖顯微組分特征及其對頁巖含氣性的影響”(1309RTSA041)

賀聰(1988-)，男，2012年畢業于中國礦業大學(北京)地質工程專業，現為中國科學院礦物學、巖石學、礦床學專業在讀博士研究生，主要從事石油與天然氣地質研究

蘇奧(1989-)，男，工程師，2011年畢業于中國地質大學(武漢)石油工程專業，2014年畢業于該校礦產普查與勘探專業，獲碩士學位，現主要從事盆地流體地質油氣成藏及地球化學研究工作。

10.3969/j.issn.1006-6535.2015.06.006

TE122.3

A

1006-6535(2015)06-0025-05