——習(xí)近平總書記在致中國科學(xué)院建院70周年賀信中作出的“兩加快一努力”重要指示要求
——中國科學(xué)院辦院方針
語音播報(bào)
地震震相數(shù)據(jù)是研究地球內(nèi)部結(jié)構(gòu)的重要資料。隨著近些年來地震觀測的不斷增多,地震記錄的數(shù)量飛速增長,通過人工的方式已經(jīng)無法及時(shí)高效地從這些海量數(shù)據(jù)中提取震相走時(shí)信息。盡管研究人員提出了許多經(jīng)典的自動(dòng)拾取算法,然而這些方法往往對噪聲組成復(fù)雜的數(shù)據(jù)適應(yīng)性差且拾取精度無法達(dá)到人工拾取水平。因此,在進(jìn)行層析成像等研究工作中,研究人員仍需通過人工或半自動(dòng)的方式拾取大量的震相,這限制了相關(guān)研究工作的效率和數(shù)據(jù)的利用程度。如何對海量地震數(shù)據(jù)進(jìn)行快速、高精度的自動(dòng)拾取,是影響地震學(xué)家對地球深部精細(xì)結(jié)構(gòu)進(jìn)行研究的重要問題。
針對上述問題,中國科學(xué)院地質(zhì)與地球物理研究所地球與行星物理重點(diǎn)實(shí)驗(yàn)室副研究員王建團(tuán)隊(duì)與中國科學(xué)院大學(xué)電子與通信工程學(xué)院博士劉暢合作,以計(jì)算機(jī)視覺邊界檢測領(lǐng)域前沿的深度學(xué)習(xí)算法為基礎(chǔ),提出了基于深度學(xué)習(xí)的端到端地震震相自動(dòng)拾取算法,命名為PickNet(圖1)。它將歸一化的原始波形記錄作為數(shù)據(jù),輸出一個(gè)類脈沖序列,序列的最大值對應(yīng)的時(shí)刻即為震相的到時(shí)(圖2)。此前,震相拾取問題被看作回歸問題或者是語義分割問題,基于深度學(xué)習(xí)算法的PickNet算法將其看作邊界檢測問題,將地震波形序列映射為尖銳的到時(shí)脈沖序列,具有更高的拾取精度和泛化能力。
他們在日本Hi-net臺(tái)網(wǎng)的實(shí)際地震數(shù)據(jù)上檢驗(yàn)了PickNet算法并與人工拾取結(jié)果進(jìn)行比較。從拾取精度方面看,在PickNet和人類專家共同拾取的震相中,73.91%(85.41%)的初至P波偏差在0.05秒(0.10秒)之內(nèi),60.75%(77.47%)的初至S波偏差在0.10秒(0.20秒)之內(nèi)(圖3)。從拾取數(shù)目方面看,PickNet獲得了約8倍于臺(tái)網(wǎng)中心提供的地震震相,對地震數(shù)據(jù)的利用更加充分。為了進(jìn)一步驗(yàn)證PickNet方法可以有效服務(wù)于地球深部結(jié)構(gòu)研究,他們將PickNet拾取的震相直接用于層析成像研究,僅使用了299個(gè)地震獲得的成像結(jié)果,與此前使用數(shù)千個(gè)地震得到的速度結(jié)構(gòu)較一致(圖4)。研究人員還利用該算法對中國地震局臺(tái)網(wǎng)、美國地震臺(tái)網(wǎng)記錄的數(shù)據(jù)進(jìn)行了測試,均取得了良好的效果。
PickNet方法具有以下優(yōu)點(diǎn):
(1)準(zhǔn):拾取精度接近人工拾取結(jié)果,拾取后的震相數(shù)據(jù)能直接用于研究地球內(nèi)部結(jié)構(gòu)。
(2)快:用一塊Nvidia Geforce GTX 1080Ti GPU約4分鐘能拾取20多萬條地震波形記錄。
(3)多:拾取的走時(shí)數(shù)量是臺(tái)網(wǎng)中心提供的數(shù)量的幾倍。
(4)廣:理論上可適用于拾取其它后續(xù)震相。
該研究表明深度學(xué)習(xí)可以從海量的地震波形數(shù)據(jù)中自動(dòng)高效地挖掘震相數(shù)據(jù),為地震學(xué)家提供了有效的工具,能加速人們對地球內(nèi)部結(jié)構(gòu)的了解。
研究成果發(fā)表于JGR: Solid Earth。
圖1 PickNet算法網(wǎng)絡(luò)結(jié)構(gòu)圖。其中三角形表示上采樣,圓圈與十字的組合表示殘差單元
圖2 PickNet輸入輸出示意圖。(a和b),P波、S波對應(yīng)臺(tái)站記錄(輸入);(c和d),震相到時(shí)脈沖序列(理想輸出);(e和f),PickNet輸出(實(shí)際輸出)
圖3 PickNet在300個(gè)日本區(qū)域地震對應(yīng)的Hi-net臺(tái)網(wǎng)波形記錄上的拾取結(jié)果與JMA提供的人工拾取震相對比。(a)潛在震相數(shù)目(即地震數(shù)目乘以臺(tái)站數(shù)目)、PickNet自動(dòng)拾取數(shù)目、JMA人工拾取數(shù)目的對比。當(dāng)波形中的初至難以拾取時(shí),PickNet輸出會(huì)全部為零,即此時(shí)對該波形記錄做拒識(shí)處理,這對保證所拾取震相的質(zhì)量十分重要;(b)同時(shí)被PickNet和人類專家拾取的震相的到時(shí)偏差對比
圖4 將PickNet拾取的震相數(shù)據(jù)直接用于層析成像反演得到的日本俯沖帶地區(qū)P波(左),S波(右)速度異常。成像結(jié)果揭示了菲律賓海板塊俯沖形態(tài)的變化、火山下方的低速體及其與低頻地震分布之間的聯(lián)系等,這些結(jié)果與前人層析成像結(jié)果有很好的一致性
© 1996 - 中國科學(xué)院 版權(quán)所有 京ICP備05002857號(hào)-1 京公網(wǎng)安備110402500047號(hào) 網(wǎng)站標(biāo)識(shí)碼bm48000002
地址:北京市西城區(qū)三里河路52號(hào) 郵編:100864
電話: 86 10 68597114(總機(jī)) 86 10 68597289(總值班室)
編輯部郵箱:casweb@cashq.ac.cn
© 1996 - 中國科學(xué)院 版權(quán)所有 京ICP備05002857號(hào)-1 京公網(wǎng)安備110402500047號(hào) 網(wǎng)站標(biāo)識(shí)碼bm48000002
地址:北京市西城區(qū)三里河路52號(hào) 郵編:100864
電話: 86 10 68597114(總機(jī)) 86 10 68597289(總值班室)
編輯部郵箱:casweb@cashq.ac.cn
© 1996 - 中國科學(xué)院 版權(quán)所有
京ICP備05002857號(hào)-1
京公網(wǎng)安備110402500047號(hào)
網(wǎng)站標(biāo)識(shí)碼bm48000002
地址:北京市西城區(qū)三里河路52號(hào) 郵編:100864
電話:86 10 68597114(總機(jī))
86 10 68597289(總值班室)
編輯部郵箱:casweb@cashq.ac.cn