壓縮感知與形態(tài)濾波在勘探地震數據處理中的應用研究
發(fā)布時間:2020-05-19 12:09
【摘要】:地震信號處理是勘探地震學的一個重要組成部分。本文詳細地討論了兩種數字信號處理理論,即數學形態(tài)學與壓縮感知,以及它們在勘探地震數據噪音壓制與弱信號提取,重構與插值中的應用?碧降卣饠祿行盘柵c噪音的分離是一個研究已久,卻又沒有完全解決的問題。其原因之一是由于地震數據的復雜性,很難用一種或者幾種屬性將信號與噪音完全分離。傳統(tǒng)的方法利用頻率、波數、振幅等屬性來壓制噪音。但是在某些情況下,信號與噪音在這些屬性上的差異較小,使得傳統(tǒng)方法的應用效果不佳。對于此,本文研究了一種新的地震信號屬性,即形態(tài)尺度,并嘗試從波形差異的角度來分離信號和噪音。針對不同類型的地震噪音,依次發(fā)展了四類形態(tài)濾波方法。具體為:1)建立了形態(tài)尺度與地震信號頻率的關系,利用時間域形態(tài)濾波壓制地震低頻噪音,且保存信號的低頻成分;2)受形態(tài)濾波壓制低頻噪音啟發(fā),將形態(tài)濾波從時間域發(fā)展到空間域,利用規(guī)則噪音在其軌跡方向上的連續(xù)性對其進行壓制(低波數噪音壓制)。結合軌跡追蹤技術,空間域形態(tài)濾波可壓制線性噪音、擬線性噪音、地滾波、外源噪音、多次波等多種相干噪音;3)利用隨機噪音在時間和空間方向上都具有較差相關性的特點,將一維形態(tài)濾波推廣到時間-空間域二維形態(tài)濾波,同時利用隨機噪音與信號在時間和空間方向上形態(tài)尺度差異來壓制隨機噪音;4)針對地震微弱信號的特點,發(fā)展了正則化的非穩(wěn)態(tài)多尺度形態(tài)重構方法。將弱信號識別與提取問題轉化成一個反演問題,并利用整形正則化和共軛梯度法進行求解。地震數據重構技術可以為反演或成像提供一個高密度采樣、規(guī)則的輸入數據。本文研究了壓縮感知框架下的地震數據重構技術,詳細地討論了地震信號的低維/稀疏表示方法。首先,推導了一種新的信號低維表示方法——雙重最小二乘投影法。該方法從主成分分析法出發(fā),解決了傳統(tǒng)低維表示中的“信號估計偏離”問題。由第一重最小二乘投影找到觀測數據在一個信號維度超平面上的最佳估計;由第二重最小二乘投影找到真實信號在另一個信號維度超平面上的最佳估計。雙重最小二乘投影法可以從觀測數據中得到一個真實信號的估計(低維表示)。本文同時也給出了該方法的幾何意義,給予了理論保證。其次,傳統(tǒng)方法利用數學基(例如小波基和曲波基)來對地震數據進行稀疏表示,將地震數據看成一種“圖像”來處理,并沒有利用數據背后的地震波傳播的物理機制與規(guī)律。對于此,本文提出利用Dreamlet來表示地震數據。Dreamlet能自動滿足波動方程,因此可以獲得一個更稀疏的地震數據表示。除此之外,本文推導了一種阻尼Dreamlet稀疏表示方法。在基本Dreamlet稀疏表示中引入一個阻尼算子,能進一步提高Dreamlet的稀疏表示效果。
【圖文】:
圖 2.3 灰值形態(tài)學基本運算示例;圖中黑色的實線為輸入信號,藍色的點線對應于膨脹運算的結果,綠色的點線對應于腐蝕運算的結果,青色的虛線對應于開啟運算的結果,紫色的虛線對應于閉合運算的結果,紅色的小半圓對應于結構元素Fig. 2.3 A demonstration of four basic grayscale morphological operators. The solid black lineis the input signal. The dotted blue and green lines correspond to the dilation and erosiontransformations. The dashed cyan and magenta lines correspond to the opening and closingoperations. The red semicircle denotes structuring element圖 2.3 為灰值形態(tài)學腐蝕、膨脹、開啟、閉合的示意圖。膨脹和腐蝕可以認為是兩個小球(圖中兩個紅色小半圓)在原始數據上從左到右滾動,,小球中心點(圖中紅點)的軌跡即為膨脹和腐蝕的結果。從圖中可以看出,腐蝕和膨脹具有平滑信號的作用。開啟和閉合為腐蝕和膨脹的組合運算,相當于進行二次平滑、去除毛刺的作用。在地震信號處理中,我們一般利用灰值形態(tài)學運算來進行處理。此時, d ( t )為輸入數據, b( )為給定的結構元素。因此,在下文中如不作特殊說明,形態(tài)學即指的是灰值形態(tài)學。
A 0為結構元素的高度。在傳統(tǒng)形態(tài)學中[4],某一個結構元素通常由形狀、寬度、高度三個參數確定。在圖2.4 中,直線型、半圓型和三角型這三種形狀的結構元素的高度均為 1,而寬度均為 9。
【學位授予單位】:中國石油大學(北京)
【學位級別】:博士
【學位授予年份】:2018
【分類號】:P631.44
本文編號:2670869
【圖文】:
圖 2.3 灰值形態(tài)學基本運算示例;圖中黑色的實線為輸入信號,藍色的點線對應于膨脹運算的結果,綠色的點線對應于腐蝕運算的結果,青色的虛線對應于開啟運算的結果,紫色的虛線對應于閉合運算的結果,紅色的小半圓對應于結構元素Fig. 2.3 A demonstration of four basic grayscale morphological operators. The solid black lineis the input signal. The dotted blue and green lines correspond to the dilation and erosiontransformations. The dashed cyan and magenta lines correspond to the opening and closingoperations. The red semicircle denotes structuring element圖 2.3 為灰值形態(tài)學腐蝕、膨脹、開啟、閉合的示意圖。膨脹和腐蝕可以認為是兩個小球(圖中兩個紅色小半圓)在原始數據上從左到右滾動,,小球中心點(圖中紅點)的軌跡即為膨脹和腐蝕的結果。從圖中可以看出,腐蝕和膨脹具有平滑信號的作用。開啟和閉合為腐蝕和膨脹的組合運算,相當于進行二次平滑、去除毛刺的作用。在地震信號處理中,我們一般利用灰值形態(tài)學運算來進行處理。此時, d ( t )為輸入數據, b( )為給定的結構元素。因此,在下文中如不作特殊說明,形態(tài)學即指的是灰值形態(tài)學。
A 0為結構元素的高度。在傳統(tǒng)形態(tài)學中[4],某一個結構元素通常由形狀、寬度、高度三個參數確定。在圖2.4 中,直線型、半圓型和三角型這三種形狀的結構元素的高度均為 1,而寬度均為 9。
【學位授予單位】:中國石油大學(北京)
【學位級別】:博士
【學位授予年份】:2018
【分類號】:P631.44
【參考文獻】
相關期刊論文 前9條
1 李會儉;王潤秋;曹思遠;藥鑫蕊;王芳琳;孫立鵬;;利用多尺度形態(tài)學識別微地震監(jiān)測中的弱信號[J];石油地球物理勘探;2015年06期
2 張華;陳小宏;吳信民;;基于壓縮感知理論與傅立葉變換的地震數據重建(英文)[J];Applied Geophysics;2013年02期
3 吳靜;吳志強;胡天躍;何玉華;王璞;閆桂京;李琳;;基于構建虛同相軸壓制地震層間多次波[J];地球物理學報;2013年03期
4 唐剛;馬堅偉;楊慧珠;;基于學習型超完備字典的地震數據去噪(英文)[J];Applied Geophysics;2012年01期
5 高建軍;陳小宏;李景葉;劉國昌;馬劍;;基于POCS方法指數閾值模型的不規(guī)則地震數據重建(英文)[J];Applied Geophysics;2010年03期
6 張釗;龔育齡;王瀟;李國法;;f-x域濾波方法的應用研究[J];山西能源與節(jié)能;2010年03期
7 劉洋;劉財;王典;李勤學;馮fE;王建民;;時變中值濾波技術在地震隨機噪聲衰減中的應用[J];石油地球物理勘探;2008年03期
8 王潤秋,鄭桂娟,付洪洲,李青;地震資料處理中的形態(tài)濾波去噪方法[J];石油地球物理勘探;2005年03期
9 蔡加銘,周興元,吳律;f-x域算子外推去噪技術研究[J];石油地球物理勘探;1999年03期
本文編號:2670869
本文鏈接:http://www.sikaile.net/kejilunwen/diqiudizhi/2670869.html