三參數小波多尺度螞蟻追蹤技術研究與應用

樂友喜， 趙 迎*， 張會娟， 劉兵卿， 黃健良

(1.中國石油大學(華東) 地球科學與技術學院，青島 266580；2.中海油研究總院，北京 100027)

?

三參數小波多尺度螞蟻追蹤技術研究與應用

樂友喜1， 趙 迎1*， 張會娟1， 劉兵卿1， 黃健良2

(1.中國石油大學(華東) 地球科學與技術學院，青島 266580；2.中海油研究總院，北京 100027)

螞蟻追蹤技術是近年發展起來的一項新型三維地震解釋技術，可實現三維地震資料中斷層的精細解釋。由于斷層具有不同的尺度和規模，該技術只是對原始地震數據所能刻畫的所有尺度的斷層綜合在一起進行追蹤，不具有針對不同級別的斷層進行變尺度分析的能力。針對這一問題，將三參數小波多尺度分析應用到螞蟻屬性體的提取和解釋中，提出了一種三參數小波多尺度螞蟻追蹤技術，并應用該技術對實際地震資料進行了斷層精細解釋。結果表明，大尺度螞蟻追蹤屬性體較好地刻畫了大尺度斷層的展布及其特征；小尺度螞蟻屬性體對非連續性細節信息顯示更豐富，能夠反映出原始數據體不易識別的一些微小斷層。因此該技術可有效提高斷層解釋的精度和可靠性，能夠滿足不同勘探開發階段對斷裂系統解釋的要求。

三參數小波； 螞蟻追蹤； 多尺度； 斷層解釋

0 引言

在三維地震資料解釋中，斷裂系統的解釋是其中的一項重要內容，其解釋的準確程度，對于沉積盆地內構造斷裂系統的研究，有利油氣區帶的識別，局部小斷裂的精細識別與描述，圈閉邊界封堵性的研究，鉆井井位的正確設計等具有非常重要的影響[1]。常規的地震資料斷層解釋周期長，且存在較強的主觀性和經驗性，低級序斷層識別也比較困難。隨著新技術、新方法的不斷發展，以相干體技術、方差體技術等為代表的一系列地震屬性處理技術的出現，在一定程度上提高了斷層解釋的精度和速度，但在具體應用過程中，會受到地震分辨率等多種因素的制約，存在一定的局限性[2]。因此地震、地質工作者一直都期望能研究出一套更客觀、精準的斷裂系統解釋技術，基于蟻群算法的螞蟻追蹤技術給人們帶來了滿意的答案，它能夠克服以往對斷層解釋的主觀性和經驗性，有效提高斷層解釋的精度、速度和對地質、構造的認知程度，然而地下斷層具有多級尺度的響應，該技術卻不具有多尺度分析和解釋斷層的性能。2005年高靜懷等[3]提出一類新的分析小波—三參數小波，該小波比Morlet小波、最佳匹配地震子波等小波具有更好的時域局部化特性，這里將三參數小波多尺度分析應用到螞蟻屬性體的提取和解釋中，提出了三參數小波多尺度螞蟻追蹤技術，并針對南堡凹陷部分區域的三維地震資料進行了斷裂系統的精細解釋。

1 螞蟻追蹤技術的基本原理和工作流程

1.1 基本原理

20世紀90年代初，意大利學者Dorigo等[4-5]受自然界中真實蟻群的集體覓食行為的啟發首先提出蟻群算法，它是一種基于群體的隨機搜索模擬進化算法。螞蟻在覓食的運動過程中，能夠在經過的路徑上留下一種稱之為外激素的物質，在運動過程中并能感知這種物質，以此來指導自己的運動方向，因此由大量螞蟻組成的蟻群集體行為便表現出一種信息正反饋現象：某一路徑上走過的螞蟻越多，留下的信息就越多，后來者選擇該路徑的概率就越大[6]。Petrel軟件自動構造解釋模塊下的螞蟻追蹤技術以蟻群算法為原理，在地震數據體中散播大量的螞蟻，當有螞蟻發現滿足預設斷裂條件的斷裂痕跡時將釋放某種信號，召集其他區域的螞蟻集中在該斷裂處對其進行追蹤，直到完成該斷裂的追蹤和識別，而其他不滿足斷裂條件的斷裂痕跡將不再進行標注，最終通過運算獲得一個低噪音、具有清晰斷裂痕跡的螞蟻屬性體。

1.2 工作流程

螞蟻追蹤技術進行斷裂系統解釋的工作流程主要包括4個步驟：

1)地震資料預處理。主要采用構造平滑處理、混沌處理、作方差體等邊緣檢測手段增強地震數據在空間上的不連續性，并可通過降低噪音來任意限定地震數據體[7]。

2)利用螞蟻算法生成螞蟻屬性體。螞蟻追蹤技術創立了一種全新的斷裂系統屬性，在預先設定的地震體內突出具有方位特征的斷裂，然后進行運算并產生螞蟻屬性體。該步驟是斷裂系統解釋的核心工作，其成功的關鍵在于參數選擇[8]。

3)提取斷片，同時進行驗證和編輯。

4)建立最終的斷裂解釋模型。

2 三參數小波變換基本原理

高靜懷等[3,9]由匹配地震子波的物理小波出發，提出了一類新的分析小波——三參數小波，該小波是對地震子波模擬公式的修改而得到的，有三個可調參數，對信號做小波分析時有很大的自由度，能夠很好地匹配給定的地震子波或有效信號，與最佳匹配地震子波或其他小波相比，具有更好的時域局部化性質，其時域表達式為式(1)：

φ(t;Λ)=e-τ(t-β)2{p(Λ)[cos(σt)-k(Λ)]+iq(Λ)sin(σt)}

(1)

式中：矢量Λ=(σ,t,β)表示參數σ、τ、β的一組集合；σ為分析小波的調制頻率；τ為能量衰減因子；β為能量延遲因子，(σ,τ,β∈R且σ,τ≥0)。通過對式(1)作傅里葉變換，得其頻域形式為式(2)。

(2)

利用容許條件和分析小波的歸一化條件，可以解得p(Λ)、q(Λ)和k(Λ)。

(3)

(4)

(5)

由式(1)構造的小波函數即為三參數小波，將三參數小波作為基本小波，對于任意給定的信號s(t)∈L2(R)，L2表示平方可積函數空間，s(t)關于分析小波φ(t)的三參數小波變換為：

(6)

式中：a為尺度因子；b為平移因子，a,b∈R且a≠0；φ*表示取復共軛。

尺度因子a對應于三參數小波的主頻。根據Meyersetal[10]提出的方法，可得到尺度因子與小波主頻的關系式：

(7)

式中：f是Fourier頻率；ω*為中心頻率(即平均頻率)；a為尺度因子。若使式(7)對三參數小波成立，則須取能量延遲因子β=0。根據定義，令β=0，可計算出三參數小波的平均頻率為式(8)。

(8)

在應用三參數小波變換進行時頻分析的過程中，可通過參數σ、τ、β調節時頻譜的時間和頻率分辨率。由圖1可知，隨著σ減小，τ增大，時頻譜的時間分辨率增高，頻率分辨率降低；反之，則時頻譜的時間分辨率降低，頻率分辨率增高。

3 三參數小波多尺度螞蟻追蹤技術

小波變換具有很好的時頻局部化性質，它的多分辨率分析和優良的“數學顯微鏡”特性，能豐富解釋信息，提高解釋精度，有利于對油藏進行精細解釋和描述。小波變換具有的良好性質，使其在地震資料分析中得到了廣泛地應用，這些應用一般分為兩大類：①利用小波多分辨率對地震信號進行小波變換，再進行分頻去噪，以提高地震資料的分辨率及信噪比；②在小波域中進行地震屬性分析與解釋方法研究，將研究結果直接用于斷層或構造解釋、儲層參數預測等。例如，小波多尺度相干分析技術、三維多尺度體曲率分析技術等。

圖2比較了不同時頻分析方法的時頻譜效果，從圖2可以看出，三參數小波變換與其他方法相比，具有更好的時域局部化性質，這里選取三參數小波做小波多尺度螞蟻追蹤技術研究。與小波多尺度相干體的提取思想相似，首先采用Al-Dossary等[11]提出的保構造濾波方法，對原始3D地震數據進行預處理，有效地去除噪聲干擾，同時能保留縱橫向邊緣信息；然后對預處理后的數據進行頻譜分析，確定地震資料的頻帶范圍，在此基礎上選取若干個合適的分析尺度，對地震數據做三參數小波變換分頻處理，獲取不同尺度性質的分頻數據體；再逐個對分頻體提取螞蟻屬性體，最后對多尺度螞蟻屬性體進行對比分析。

圖1 參數σ，τ，β對時頻譜時間和頻率分辨率的影響

圖2 不同時頻分析方法的時頻譜效果對比

4 實際應用效果分析

南堡凹陷位于渤海灣盆地黃驊坳陷北部，是在華北板塊基底上以新生代為主要時期發育的“北斷南超”型斷陷盆地。內部斷層十分發育且走向多變,以北東東、東西向為主,存在部分北東、北西向斷層[12]。因此該區構造復雜，斷層走向多樣，剖面解釋以及空間組合都存在較大的困難。針對南堡凹陷三維地震資料，這里應用三參數小波多尺度螞蟻追蹤技術，對部分區域的斷裂系統進行了精細解釋。

圖3(a)為南堡凹陷3號構造帶三維地震資料Inline1780測線的原始地震剖面，從圖3(a)可以看出，部分斷層的斷點有顯示但較模糊，斷面解釋可靠性差，斷層組合難度較大。經三參數小波變換分頻處理后分別得到主頻為15 Hz、30 Hz、50 Hz的地震剖面(如圖3(b)，3(c)，3(d)所示)：①當主頻為15 Hz時(圖3(b))，地震分辨率較低，信噪比較高，圖3(b)中橢圓框內所示的大尺度斷層其斷點較為清晰，斷面上下貫通，解釋的可靠性較高；②當主頻為30 Hz時(圖3(c))，雖然大尺度斷層的深層斷點變得模糊，但是中淺層小尺度斷層其斷點變得更為清晰；③當主頻為50 Hz時(圖3(d))，地震分辨率較高，小尺度斷層非常發育，大尺度斷層已不能可靠解釋。

在三參數小波多尺度分頻數據體的基礎上，對分頻數據體分別提取螞蟻屬性體，然后對多尺度螞蟻屬性體進行對比分析。從多尺度螞蟻屬性體的水平切片上可以了解斷裂的平面展布狀態、主要發育規律，從而對研究區的構造特征形成一個基本認識。圖4(a)為原始地震數據的螞蟻屬性體在2 750 ms處的水平切片，圖4(b)、4(c)、4(d)分別為三參數小波分頻數據體的螞蟻屬性體在2 750 ms處的水平切片，圖中黑色線條表示連續性很差、邊緣變化劇烈的地方。從圖4可以看出，大尺度螞蟻屬性體的水平切片(圖4(b))中線條較稀疏、延續長度較長，表示大尺度螞蟻追蹤結果主要追蹤了大斷層，較大尺度斷層的發育和展布規律較為清晰，不僅弱化了一些原始數據體中可能由噪聲引起的偽邊緣信息，而且一些原始數據體的螞蟻追蹤結果中(圖4(a))不夠清晰或者受其他細小邊緣的干擾，導致展布刻畫不清晰的主要大斷層在大尺度螞蟻屬性體的水平切片中得以很好地描述，如圖中西南角部分。隨著尺度的降低，越來越多的小尺度斷層或低級序斷層在螞蟻體水平切片中被反映出來，其形態和展布也比較清晰直觀，而且原始數據體中沒有識別出的一些小斷裂，在小尺度螞蟻追蹤結果中(圖4(d))中得以清晰地反應出來。綜合分析多個不同尺度下的螞蟻追蹤結果，不僅可以清晰刻畫主要大斷層的空間發育與展布，也能清晰揭示小尺度斷裂信息，從而提高了斷層解釋的精度和可靠性，能夠滿足不同勘探開發階段對斷裂系統解釋的要求。

圖3 Inline1780原始地震剖面及三參數小波分頻剖面

圖4 螞蟻屬性體的水平切片對比圖(2 750 ms)

5 結論

小波變換具有很好的時頻局部化性質，它的多分辨率分析和優良的"數學顯微鏡"特性，能豐富解釋信息，提高解釋精度，三參數小波比Morlet小波、最佳匹配地震子波等小波具有更好的時頻聚焦性；螞蟻追蹤技術能夠克服對斷層解釋的主觀性和經驗性，有效提高斷層解釋的精度、速度和對地質、構造的認知程度。將三參數小波變換與螞蟻追蹤技術有效地結合起來，充分發揮兩者的優勢，實現對不同級別的斷層的變尺度分析，極大地提高斷層解釋的精度和可靠性。

[1] 張磊.螞蟻體追蹤技術在潛西地區構造解釋中的運用[J].江漢石油科技，2013，23(3)：24-27. ZHANG L. Application of ant tracking technology in structural interpretation in Qianxi area. Jianghan Petroleum Science and Technology, 2013, 23(3): 24-27. (In Chinese)

[2] 張繼標，戴俊生，馮建偉，等.螞蟻追蹤技術在大程莊地區斷裂自動解釋中的應用[J].石油天然氣學報，2012，34(5)：53-57. ZHANG J B, DAI J SH, FENG J W, et al. Application of ant tracking technology in automatic fault interpretation in Dachengzhuang area. Journal of Oil and Gas Technology, 2012, 34(5): 53-57. (In Chinese)

[3] 萬濤，高靜懷，陳文超，等.三參數小波及其在薄互層地震資料分析中的應用[J].煤田地質與勘探，2005，33(4)：69-72. WAN T, GAO J H, CHEN W C, et al. TP wavelet and its application to analysis of seismic data of thin interbeds. Coal Geology & Exploration, 2005, 33(4): 69-72. (In Chinese)

[4] DORIGO M, MANIEZZO V, COLORNI A. Positive Feedback as a Search Strategy[R]. Milan: Milan Politecnico di Milano, 1991.

[5] DORIGO M, GAMBARDELLA M. Ant colony system: a cooperative learning approach to the traveling sales-man problem[J]. IEEE Transactions on Evolutionary Computation, 1997, 1(1): 53-66.

[6] 雷筱珍，賴萬欽.一種改進的并行蟻群優化求解算法[J].甘肅聯合大學學報:自然科學版，2009，23(1)：67-71. LEI X ZH, LAI W Q. An improved parallel ant colony optimization algorithm. Journal of Gansu Lianhe University, 2009, 23(1): 67-71. (In Chinese )

[7] 唐琪凌，蘇波，王迪，等.螞蟻算法在斷裂系統解釋中的應用[J].特種油氣藏，2009，16(6)：30-33. TANG Q L, SU B, WANG D, et al. Application of ant colony algorithm in fault system interpretation. Special Oil and Gas Reservoirs, 2009, 16(6): 30-33. (In Chinese)

[8] 張繼標，戴俊生，趙力彬，等.基于螞蟻算法的斷裂自動解釋技術在黃玨南地區的應用[J].中國石油大學學報:自然科學版，2011，35(6)：14-20. ZHANG J B, DAI J S, ZHAO L B, et al. Application of automatic fault interpretation technology based on ant colony algorithm in south Huang-jue area. Journal of China University of Petroleum, 2011, 35(6): 14-20. (In Chinese)

[9] 高靜懷，萬濤，陳文超，等.三參數小波及其在地震資料分析中的應用[J].地球物理學報，2006，49(6)：1802-1812. GAO J H, WAN T, CHEN W C, et al. Three parameter wavelet and its application to seismic data processing. Chinese Journal of Geophysics, 2006, 49(6): 1802-1812. (In Chinese)

[10]MEYERS S D, KELLY B G, O’BRIEN J J. An introduction to wavelet analysis in oceanography and meteorology: With application to the dispersion of Yanai waves[J]. Monthly Weather Review, 1993, 121(10): 2858-2866.

[11]Al-DOSSARY S, WANG Y E, ARAMCO S. Structure preserving smoothing for 3D seismic attributes[C]. SEG 81thannual meeting, Society of Exploration Geophysicists, Las Vegas, 2011, 1004-1008.

[12]周天偉，周建勛，董月霞，等.渤海灣盆地南堡凹陷新生代斷裂系統形成機制[J].中國石油大學學報:自然科學版，2009，33(1)：12-17. ZHOU T W, ZHOU J X, DONG Y X, et al. Formation mechanism of cenozoic fault system of Nanpu sag in Bohai bay basin. Journal of China University of Petroleum, 2009, 33(1): 12-17. (In Chinese)

Study and application of three-parameter wavelet multi-scaleant tracking technology

YUE You-xi1, ZHAO Ying1*, ZHANG Hui-juan1, LIU Bing-qing1, HUANG Jian-liang2

(1.School of Geosciences,China University of Petroleum, Qingdao 266580,China;2.CNOOC Research Institute,Beijing 100027,China)

Ant tracking technology is a newly developed technology for 3D seismic data interpretation, and fine fault interpretation can be well accomplished by using this technology. Since faults have different scales and sizes, and this technology can only synthetically track all scale of faults together which are depicted by original seismic data. So this method does not have the variable scale analysis ability for different levels of faults. To solve this problem, we apply three-parameter wavelet multi-scale analysis to the extraction and interpretation of ant tracking attributes, and propose a three-parameter wavelet multi-scale ant tracking technology. Applications of this method to real fault system interpretation show that large-scale ant tracking attributes mainly depicte the distribution and characteristic of large-scale faults, while smaller-scale ant tracking attributes display more abundant non-continuous detail information and reflect some minor faults which can not be identified easily from original seismic data. Therefore, this technology can effectively improve accuracy and reliability of fault interpretation and meet the requirements of fracture system interpretation in different stages of exploration and development.

three parameters wavelet; ant tracking; multi-scale; fault interpretation

2014-10-08改回日期：2015-04-15

樂友喜(1966-)，男，教授，博士，主要研究方向為地震儲層預測、信號處理與分析、模式識別、開發地震等，E-mail：yueyouxi@upc.edu.cn。

*通信作者：趙迎(1990-)，男，碩士，主要研究方向為地震儲層預測及信號處理與分析，E-mail：zhaoying_211@126.com。

1001-1749(2015)05-0610-06

P 631.4

A

10.3969/j.issn.1001-1749.2015.05.11