數(shù)字處理ch32子波提取與子波整形反褶積

《數(shù)字處理ch32子波提取與子波整形反褶積》由會員分享，可在線閱讀，更多相關(guān)《數(shù)字處理ch32子波提取與子波整形反褶積（42頁珍藏版）》請?jiān)谘b配圖網(wǎng)上搜索。

1、單擊此處編輯母版標(biāo)題樣式,*,單擊此處編輯母版文本樣式,第二級,第三級,第四級,第五級,3.1,反褶積的概念,3.2,最佳維納濾波,3.3,最小平方反褶積,3.4,預(yù)測反褶積,3.5,子波提取與子波整形反褶積,3.6,同態(tài)反褶積,3.7,地表一致性反褶積,3.8,反,Q,濾波及譜白化,第三章反褶積,3.5,子波提取與子波整形反褶積,希爾伯特變換法或,Wold,Kolmogorov,公式,海上及陸地勘探有差異,單位圓外無根,要求有高質(zhì)量的聲波、密度測井資料及井旁地震記錄,自,學(xué)并查相關(guān)文獻(xiàn)完成作業(yè),零相位子波,最小相位子波,3.5,子波提取與子波整形反褶積,非因果信號,零相位子波比最小相位子波有

2、更高的分辨率！,據(jù)你所知都有哪些相位的信號？,3.5,子波提取與子波整形反褶積,3.5,子波提取與子波整形反褶積,子波與反子波時域分布關(guān)系,子波 反子波,n,n,n,n,非因果的,3.5,子波提取與子波整形反褶積,？,如何獲得？,找同學(xué)上來求！,3.5,子波提取與子波整形反褶積,？,當(dāng)期望輸出與子波具有相同的相位（最大）延遲時，誤差能量最?。?3.5,子波提取與子波整形反褶積,子波,期望輸出,整形反褶積算子,褶積輸出,最佳結(jié)果,子波整形反褶積前后資料的比較,3.5,子波提取與子波整形反褶積,掃描信號的相關(guān)得到接近零相位的子波,存在漂移的原因？,最小相位,零相位,消除了漂移，匹配良好,消除了漂移

3、，匹配良好,3.1,反褶積的概念,3.2,最佳維納濾波,3.3,最小平方反褶積,3.4,預(yù)測反褶積,3.5,子波提取與子波整形反褶積,3.6,同態(tài)反褶積,3.7,地表一致性反褶積,3.8,反,Q,濾波及譜白化,第三章反褶積,一復(fù)賽譜,Cepstrum,和,D,特征系統(tǒng),3.6,同態(tài)反褶積,一復(fù)賽譜,Cepstrum,和,D,特征系統(tǒng),兩邊取對數(shù),記為,傅氏反變換,該詞來自英文,Cepstrum,，,Spectrum,的反寫,由褶積變?yōu)榧臃ǎ?3.6,同態(tài)反褶積,3.6,同態(tài)反褶積,3.6,同態(tài)反褶積,兩邊取對數(shù),傅氏反變換,泰勒展開,復(fù)賽譜,3.6,同態(tài)反褶積,對復(fù)賽譜作高通濾波，通過逆,D,

4、系統(tǒng)，可得反射系數(shù)，實(shí)現(xiàn)同態(tài)反褶積。,同態(tài)反褶積前,同態(tài)反褶積后,3.1,反褶積的概念,3.2,最佳維納濾波,3.3,最小平方反褶積,3.4,預(yù)測反褶積,3.5,子波提取與子波整形反褶積,3.6,同態(tài)反褶積,3.7,地表一致性反褶積,3.8,反,Q,濾波及譜白化,第三章反褶積,3.7,地表一致性反褶積,Surface Consistent,Deconvolution,3.7,地表一致性反褶積,常假定為最小相位,3.7,地表一致性反褶積,地表一致性反褶積效果比較,(a),普通反褶積,(b),地表一致性反褶積,3.7,地表一致性反褶積,（,分辨率提高、振幅均勻,）,3.1,反褶積的概念,3.2,最

5、佳維納濾波,3.3,最小平方反褶積,3.4,預(yù)測反褶積,3.5,子波提取與子波整形反褶積,3.6,同態(tài)反褶積,3.7,地表一致性反褶積,3.8,反,Q,濾波及譜白化,第三章反褶積,二、譜白化,一、反,Q,濾波,3.8,反,Q,濾波及譜白化,由于地層的吸收作用，地震波經(jīng)地層傳播后，能量被衰減損耗，頻率變低。特別是深層，分辨率大大下降。,一、反,Q,濾波,為恢復(fù)地震波原來的能量，提高分辨率，處理時必須做吸收補(bǔ)償，即,Q,補(bǔ)償,。從濾波的角度講，就是,反,Q,濾波,。,另外，,Q,因子本身也是一種很好的吸收屬性,它還可在油藏描述中得到應(yīng)用。,粘彈性介質(zhì)衰減,地震資料處理的常規(guī)處理步驟之一,Inela

6、stic Attenuation,Q,品質(zhì)因子，,quality factor,反,Q,濾波（,Q,補(bǔ)償）,Q Compensation,振幅衰減主要是幾何（球面）擴(kuò)散與,吸收,。,1,、,Q,因子的物理含義,Q,因子反映的就是能量被吸收損耗的比率，一個波長,內(nèi)，原地震波的能量,E,與傳播中吸收所損耗的能量,E,之比即為,Q,。,彈性能轉(zhuǎn)換成熱能的過程,物理意義？,3.8,反,Q,濾波及譜白化,2,、,Q,值的估算,3.8,反,Q,濾波及譜白化,常有,振幅包絡(luò)法,，頻譜比法，,常,Q,掃描法,和,李氏公式法,等。,一、反,Q,濾波,1,、,Q,因子的物理含義,選一段記錄，用不同的,Q,值從小到

7、大各掃一次做反,Q,濾波，組成許多拼起來的圖，然后處理人員在這些圖上從淺至深選擇一個個,Q,值。,實(shí)際生產(chǎn)中常用,教材中有介紹，自學(xué),常,Q,掃描法：,具有速度資料時可使用,李氏經(jīng)驗(yàn)公式：,3,、反,Q,濾波的實(shí)現(xiàn)方法,大地濾波濾波器的振幅譜為,假定大地濾波過程為最小相位，利用,Hilbert,變換根據(jù)振幅譜求取相位譜：,令,則反,Q,濾波器的頻譜為,時間域的濾波因子如下：,3.8,反,Q,濾波及譜白化,改造后頻譜,Q,補(bǔ)償后，即未衰減前,Q,補(bǔ)償前，即吸收衰減后,具有時變性,，逐點(diǎn)變化,因?yàn)轳薹e因子具有逐點(diǎn)變化的時變性，無法求得精確解，需要做如下近似：,顯然，這種假設(shè)在,Q,時，即沒有衰減的

8、情況下是成立的。,3.8,反,Q,濾波及譜白化,假設(shè)在同一地層中，也就是在一段時間,t,中，,Q,有確定值，濾波因子,h,不變。這樣，就可得到：,注意：對于每一個時間,t,都要在頻率域計(jì)算一個積分，計(jì)算量是巨大的。,因此有必要采取某種措施，以避免積分運(yùn)算！,Dave Hale,（,SEP,，,1982,）將,Q(f,t),用,Taylor,級數(shù)展開，設(shè),G(f),對應(yīng)的時域信號為,g(t),，,則有：,泰勒展開,俞壽朋,高分辨率地震勘探,中對于,Q,提取及影響因素有詳細(xì)的討論。,常,Q,掃描,實(shí)例一,實(shí)例二,實(shí)例二（放大）,二、譜白化,(Spectrum,Whiltening,),譜白化是一種

9、展寬頻譜的方法。,不過它不改變子波的相位譜，,是一種,“,純振幅,”,的濾波過程,。,譜白化處理可以在,頻率域,中完成，也可以在,時間域,里進(jìn)行。,3.8,反,Q,濾波及譜白化,頻率域?qū)崿F(xiàn)方法,時間域?qū)崿F(xiàn)方法,二、譜白化,在頻率域中：,可將地震道振幅譜的各尖峰連成沒有極小值的外包線，,再將其實(shí)行平滑濾波。,再在縱坐標(biāo)上加適量“白噪成分”，然后求一個倒數(shù)，依此倒數(shù)比例放大原來各振幅譜的值，使外包線展平為“白色”的寬頻譜。,在不改原先的相位譜的情況下，做傅氏反變換，即得譜白化結(jié)果。,在時域中進(jìn)行時：,先將頻譜區(qū)間劃分成,3-4,個濾波頻段，每個濾波頻段由梯形組成，它的相鄰兩個邊是互補(bǔ)的,(,加起來

10、為,1,),。,分完濾波頻段之后，用分頻檔濾波的方法將記錄分為四個頻檔的時域形態(tài)。,然后統(tǒng)計(jì)各頻檔的平均振幅，并設(shè)法讓它們乘以不同的放大倍數(shù)，使每個頻檔的平均振幅都互相看齊。,然后再把四個頻檔加起來，即得到譜白化的效果。,二、譜白化,時變譜白化效果的好壞在于合理選擇分頻濾波門的大小，應(yīng)根據(jù)實(shí)際資料的頻率掃描結(jié)果確定！,二、譜白化,第三章作業(yè),1、詳細(xì)推導(dǎo)（329）式和（335）式，并給出每一步推導(dǎo)過程的依據(jù)。,2、設(shè)輸入子波為,(1,，,1/3,),，請根據(jù)最小平方反褶積的方法，求取反子波。,3、請給出Q值與吸收系數(shù)的關(guān)系。若已知頻率為f=40Hz，速度為v=2500m/s，請作出Q與的關(guān)系曲線。,4、子波提取方法綜述。,要求：（1）查閱文獻(xiàn)，并詳細(xì)敘述每一種方法的原理；,（2）注明參考文獻(xiàn)的出處；,（3）手寫，不少于500字。,5、反褶積方法綜述,要求：（1）查閱文獻(xiàn)，并詳細(xì)敘述每一種方法的原理；,（2）詳細(xì)敘述一種課本上未講過的方法；,（3）注明參考文獻(xiàn)的出處；,（4）手寫，不少于1500字。,作業(yè),:,1,、證明,Rxx,和,Rxz,式,2,、做一道維納濾波的題,3,、預(yù)測反褶積的題,4,、子波求取的題,5,、同態(tài)反褶積的題,6,、反褶積方法綜述（期中作業(yè)，找一種課本上未講的方法詳細(xì)敘述）,7,、子波提取方法綜述（期中作業(yè)，找一種課本上未講的方法詳細(xì)敘述）。,