• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    利用多重濾波法提取面波群速度及其高斯濾波參數(shù)取值

    2019-08-28 01:03:36蔣嬋君王有學(xué)王海燕
    桂林理工大學(xué)學(xué)報 2019年2期
    關(guān)鍵詞:群速度面波臺站

    蔣嬋君,王有學(xué),熊 彬,王海燕

    (1.桂林理工大學(xué) a.地球科學(xué)學(xué)院; b.廣西隱伏金屬礦產(chǎn)勘查重點實驗室,廣西 桂林 541006;2.桂林理工大學(xué)博文管理學(xué)院,廣西 桂林 541006)

    天然地震面波層析成像是探測地殼和上地幔結(jié)構(gòu)的最常用方法之一,其前提是從面波數(shù)據(jù)中提取高質(zhì)量的相速度或者群速度。目前,提取面波群速度的方法有移動窗分析法[1]、多重濾波法[2]和連續(xù)小波變換[3],其中多重濾波法是提取面波群速度最常用的方法[4-9]。

    Landisman等[1]提出了在時間域采用矩形窗作為時間窗口的移動窗分析法, Dziewonski等[2]提出了在頻率域采用高斯濾波器的多重濾波法。如果將移動窗分析法中的矩形窗改為高斯窗,那么時間域移動窗分析法與頻率域多重濾波法是等價的。目前,多重濾波法作為最常用的提取群速度的方法,其效果在很大程度上取決于高斯濾波參數(shù)α,因此需要選擇合適的高斯濾波參數(shù)值,前人對此作了大量的研究: Inston等[10]提出了最佳帶寬濾波; Cara[11]提出了最佳濾波方案; Nyman等[12]提出了均等顯示方法; Herrmann等[13]給出了不同震中距的高斯濾波參數(shù)的參考取值; Kolínsky[14]研究發(fā)現(xiàn),高斯濾波參數(shù)與信號周期之間存在準線性變化的規(guī)律,但僅適用于該研究中所采用的天然地震數(shù)據(jù);朱良保等[15]提出了線性時間分辨方法選擇高斯濾波參數(shù);陳浩朋等[16]提出了不同震中距的高斯濾波參數(shù)按照周期分段進行取值。本文將采用Dziewonski等[2]、Herrmann等[13]、陳浩朋等[16]分別提出的高斯濾波參數(shù)α取值,并且在陳浩朋等[16]的周期分段取值基礎(chǔ)上作了調(diào)整改進,利用多重濾波法對合成瑞利面波數(shù)據(jù)和高質(zhì)量實際數(shù)據(jù)進行群速度的提取,分析高斯濾波參數(shù)α的不同取值對多重濾波法的影響。

    1 多重濾波法原理

    1.1 單臺法

    中心頻率為ωn的高斯帶通濾波器Hn(ω)為

    (1)

    其中:α為高斯濾波參數(shù);ωn為中心頻率;ωln為中心頻率的下限頻率;ωun為中心頻率的上限頻率[2-3,16-17]。

    令BAND=(ωun-ωn)/ωn=(ωn-ωln)/ωn,β=αBAND2。 其中BAND為相對帶寬,β為截斷閾值。 實際計算中必須對高斯函數(shù)進行截斷, 原則上在截斷后不能引起明顯的吉普斯效應(yīng),β一般取值為3,截斷處的能量相對于能量峰值下降了30 DB。當高斯濾波參數(shù)為常數(shù)時,BAND也為常數(shù)時,此時稱為常相對帶寬[15]。

    設(shè)某一臺站記錄的天然地震面波數(shù)據(jù)為f(t), 對f(t)進行傅里葉變換得到f(ω); 在頻率域, 用中心頻率為ωn的高斯帶通濾波器Hn(ω)對f(ω)進行濾波得到Fn(ω); 然后將Fn(ω)進行傅里葉反變換轉(zhuǎn)換到時間域, 并對其進行希爾伯特變換得到幅值包絡(luò), 最大振幅的到時即為該頻率ωn的群速度波包的到時tn; 假設(shè)面波沿大圓路徑傳播, 根據(jù)震中距L和群速度波包到時tn, 即可計算出該中心頻率ωn的群速度Un=L/tn。 重復(fù)上述步驟, 直到求出所有待求頻率(周期)的群速度值[17]。

    1.2 雙臺法

    假設(shè)面波沿大圓路徑,當兩個臺站與震中近似在同一條大圓弧上時(圖1),可以用雙臺法提取兩個臺站間的群速度。

    設(shè)這兩個臺站記錄的天然地震面波數(shù)據(jù)分別為f1(t)、f2(t),其互相關(guān)函數(shù)f(t)為

    (2)

    圖1 雙臺法示意圖Fig.1 Diagram for two-station method

    f(t)相當于第1個臺站激發(fā)第2個臺站接收的地震記錄。 此時, 用單臺法的步驟提取f(t)的群速度[16]。

    單臺法和雙臺法只是面波的數(shù)據(jù)來源不同,單臺法適用于震中在研究區(qū)域內(nèi),雙臺法適用于震中在研究區(qū)域外。在提取群速度的原理和步驟上,單臺法和雙臺法沒有區(qū)別,故而本文均以單臺法為例。

    2 高斯濾波參數(shù)

    2.1 理論瑞利面波數(shù)據(jù)

    為了分析高斯濾波參數(shù)α的取值對多重濾波法的影響,采用模態(tài)疊加方法[13]計算理論瑞利面波數(shù)據(jù):選擇地球速度結(jié)構(gòu)AK135模型[18],地震震源參數(shù)見表1,5個接收臺站,其震中距分別為1 000、2 000、3 000、4 000和8 000 km,方位角45°,數(shù)據(jù)采樣間隔為1 s。震中距為3 000 km的理論瑞利面波數(shù)據(jù)及其振幅譜見圖2,提取的周期范圍為4~220 s。

    2.2 高斯濾波參數(shù)取值

    在多重濾波法中,選擇合適的高斯濾波參數(shù)α可以減少面波群速度測量誤差。

    表1 地震震源參數(shù)

    圖2 震中距3 000 km的理論瑞利面波數(shù)據(jù)(a)及其振幅譜(b)Fig.2 Synthetic seismogram of Rayleigh wave at distance of 3 000 km(a) and its amplitude spectra(b)

    Dziewonski等[2]給出了α=50.3;Herrmann等[13]給出了4~100 s內(nèi)不同震中距下高斯濾波參數(shù)α的參考值(表2,文獻[13]并沒有明確給出震中距3 000 km時的高斯濾波參數(shù)α參考值, 但可以根據(jù)震中距2 000 km和4 000 km的α參考值插值求取)。 采用上述兩種不同的常相對帶寬取值方案, 對不同震中距的理論瑞利面波數(shù)據(jù)進行多重濾波, 其結(jié)果分別見圖3a1—a5、 圖3b1—b5。 采用Dziewonski等[2]、Herrmann等[13]的高斯濾波參數(shù)α取值提取的群速度曲線在周期較大時均出現(xiàn)起伏,與理論群速度曲線擬合程度差;并且震中距越小,群速度曲線在周期較大時的起伏越劇烈,與理論群速度曲線擬合程度越差。這說明,高斯濾波參數(shù)α的選取若只考慮震中距的大小,在周期較大時很難取得較好的效果,因此,當周期范圍跨度較大時,高斯濾波參數(shù)α取值不僅應(yīng)考慮震中距,還應(yīng)考慮周期。

    表2 不同方法的高斯濾波參數(shù)α取值

    陳浩朋等[16]給出了不同震中距不同周期分段的高斯濾波參數(shù)α取值(表2),對不同震中距的理論瑞利面波數(shù)據(jù)進行多重濾波,其結(jié)果見圖3(c1—c5)。 除了震中距1 000 km以外,采用陳浩朋等[16]的高斯濾波參數(shù)α取值提取的群速度曲線在周期較大時沒有出現(xiàn)明顯的起伏,且與理論群速度曲線擬合較好,這是因為陳浩朋等[16]的高斯濾波參數(shù)α取值同時考慮了震中距和周期范圍。

    但當震中距為2 000、3 000、4 000和8 000 km時,尤其是當震中距為2 000、3 000 km時,采用陳浩朋等[16]的高斯濾波參數(shù)α取值提取的群速度曲線在45~60 s范圍內(nèi)與理論群速度曲線擬合相對較差,這說明陳浩朋等[16]的周期分段點60 s,取值過大,經(jīng)過試驗后將周期分段點由60 s調(diào)整為45 s(表2), 改進后提取的群速度曲線見圖3(d2—d5),改進后提取的群速度曲線在45~60 s范圍內(nèi)與理論群速度曲線擬合得更好。

    對于不同震中距,在周期小于45 s范圍內(nèi),盡管Dziewonski等[2]、Herrmann等[13]、陳浩朋等[16]的高斯濾波參數(shù)α取值不同(表2),但提取的群速度曲線均無明顯區(qū)別(圖3a1—a5, b1—b2, c1—c5),這說明在周期小于45 s范圍內(nèi),高斯濾波參數(shù)α取值范圍較大,不應(yīng)局限于某一數(shù)值,改進后的取值范圍見表2。

    另外,當震中距為1 000 km,即震中距較小時,在周期大于45 s范圍內(nèi),采取不同的高斯濾波參數(shù)α取值提取的群速度曲線均與理論群速度曲線擬合差(圖3a1,b1,c1), 這說明提取的周期不應(yīng)超過45 s(表2)。

    3 實際數(shù)據(jù)

    本文使用南嶺西部地區(qū)流動臺站觀測數(shù)據(jù)作為試驗數(shù)據(jù), 地震事件信息見表3, 共12個臺站(表4), 臺站的地震記錄見圖4, 地震記錄質(zhì)量高。 以JH01臺站記錄為例, 采用4種不同高斯濾波參數(shù)α取值方案(表2)對其進行多重濾波, 結(jié)果如圖5所示。

    采用Dziewonski等、 Herrmann等給出的常相對帶寬取值提取的群速度曲線在周期45 s后均出現(xiàn)一定程度的起伏; 采用陳浩朋等給出的周期分段α取值提取的群速度曲線在60 s后起伏消失, 但在45~60 s范圍內(nèi)仍有一定程度的起伏; 采用改進后的周期分段α取值提取的群速曲線在45~60 s范圍內(nèi)起伏也消失。 由此可見, 高斯濾波參數(shù)α取值失當會造成多重濾波法提取的群速度曲線在周期大于45 s范圍內(nèi)出現(xiàn)假象起伏, 這在一定程度上會影響后續(xù)的速度反演精度。 對12個臺站提取的群速度曲線取平均值(圖6), 在4種不同的α取值方案中, 改進后的周期分段α取值提取的平均群速度曲線最光滑。

    表3 地震事件信息

    表4 寬頻帶流動臺站信息

    圖4 地震事件不同臺站垂直分量記錄Fig.4 Vertical components record of different stations for the same earthquake event

    圖5 對JH01臺站記錄采用不同α取值提取的群速度頻散Fig.5 Group-velocity measurements record of JH01 station by MFT using different α values

    4 結(jié)論與建議

    采用Dziewonski等、Herrmann等的常相對帶寬取值、陳浩朋等的周期分段取值、改進的周期分段取值,分別對不同震中距的理論瑞利面波數(shù)據(jù)和高質(zhì)量的實際數(shù)據(jù)進行多重濾波,并對這4種不同α取值方案的提取結(jié)果進行對比,結(jié)果表明:

    (1)當震中距為1 000 km時,即震中距較小時,提取的周期不應(yīng)超過45 s。

    圖6 采用不同α取值提取的所有臺站的平均群速度頻散Fig.6 Average group-velocity measurements of all stations by MFT using different α values

    (2)當周期范圍跨度較大時,高斯濾波參數(shù)α取值不僅應(yīng)考慮震中距,而且還應(yīng)考慮周期。

    (3)對于不同震中距:在周期小于45 s時,不同的高斯濾波參數(shù)α取值對群速度的提取無明顯影響,此時α取值范圍較大,不應(yīng)局限于某一數(shù)值;但在周期大于45 s時,不同的高斯濾波參數(shù)α取值對群速度的提取影響較大,此時α的取值應(yīng)謹慎。

    (4)在周期大于45 s時,采用Dziewonski等、Herrmann等的常相對帶寬取值提取的群速度頻散曲線有明顯的起伏;采用陳浩朋等提出的周期分段α取值提取的群速度曲線雖然避免了周期大于60 s后的假象起伏,但在45~60 s范圍內(nèi)仍有一定程度的假象起伏;將陳浩朋等的周期分段點60 s調(diào)整為45 s后,避免了45~60 s范圍內(nèi)的假象起伏。因此,在4種不同的α取值方案中,本文改進后的周期分段α取值效果最好。

    猜你喜歡
    群速度面波臺站
    中國科學(xué)院野外臺站檔案工作回顧
    VTI介質(zhì)地震波群速度三維擴展各向異性線性近似表征
    氣象基層臺站建設(shè)
    西藏科技(2021年12期)2022-01-17 08:46:38
    gPhone重力儀的面波頻段響應(yīng)實測研究
    地震研究(2021年1期)2021-04-13 01:04:56
    激光調(diào)制摻鉺光纖中光速及其在高靈敏光纖傳感領(lǐng)域的潛在應(yīng)用研究
    自適應(yīng)相減和Curvelet變換組合壓制面波
    物質(zhì)波的波速與頻率公式中的能量
    小議超光速
    基層臺站綜合觀測業(yè)務(wù)管理之我見
    西藏科技(2015年6期)2015-09-26 12:12:13
    MDOS平臺臺站級使用方法及技巧
    梅河口市| 磐安县| 辛集市| 吴桥县| 乐清市| 洪湖市| 偃师市| 星座| 焉耆| 赫章县| 福清市| 阆中市| 庄浪县| 栖霞市| 嘉定区| 华宁县| 石棉县| 肥西县| 井陉县| 扎兰屯市| 普兰店市| 大连市| 尖扎县| 瓮安县| 公安县| 宜兰市| 阿拉善右旗| 阿拉善左旗| 平凉市| 康定县| 分宜县| 孟州市| 扎鲁特旗| 斗六市| 郁南县| 清徐县| 澎湖县| 柘荣县| 丰都县| 元江| 洮南市|