地震震源机制即时波形逆推系统之发展－频率波数法 - 中央气象局.PDF

地震震源機制即時波形逆推系統之發展－頻率波數 法 趙里 中央研究院地球科學研究所 張建興 甘志文 中央氣象局地震測報中心 摘 要 本研究的最終目標是開發一個計算軟體系統來實現地震發生後震源機制的即時 確定，從而把震源機制解包含在當前的地震預警系統中。為此，在第一年的研究計 畫中所要實現的任務是首先用過去已經發生的 地震為例，建立起利用初動P波波形確 定地震震源機制的電腦程式。目前這一程式已經建立起來。在建立該程式的過程中 完成了以下四項工作： 壹、 選取地震範例 從地震目錄中選取了四個發生在台南地區的三，四級淺源地震作為範例，收集 了這些地震在各種地震監測網的測站上記錄到的波形資料。其中包括中研院地球所 的BATS的寬 頻記錄，氣象局台網(CWBSN)的寬頻記錄，以及TAIGER計畫在中部和南部 所建立的流動台站的寬頻記錄。圖一顯示所選取的地震與測站分佈圖。圖中可見幾 個寬頻台網在台南地區已經能夠提供許多可供選擇的記錄。在一些寬頻台網台站較 少或一些規模較小的地震，也可以通過利用氣象局短週期台網和強地面運動台網來 彌補台站的分佈問題。本研究中僅使用了圖一中顯示的寬頻台站。 貳、 確定一維速度模型 255 建立從地震震央到測站之間的一維分層平均速度模型。因為臺灣地區地殼結構 在側向上變化很大，在不同的地區使用同一個一維速度模型將會為波形的計算帶來 較大的誤差，進而影響到震源機制解的可靠性 。因此，我們採取對不同的震央-測站 路徑採用不同的一維分層速度模型的方法。在本研究中我們選用吳逸民教授在2007 年發表的臺灣地區三維速度結構模型WYM(Wu et al. 2007)。 首先，我們用 2008年TAIGER計畫的炸測記錄分析了WYM三維速度模型的初動波到 時殘差。對於每一個炸點，用有限差分數值模擬的方法計 算WYM三維速度模型中各個 測站上的理論地震圖。然後在實際記錄和理論地震圖上量取初動到時，從而得到WYM 三維模型的到時殘差。因為我們已經掌握了爆炸源精確的時間和地點，所以利用得 到的初動到時殘差可以有效地評估速度模型的優劣。圖二中顯示了用炸點S4得到的 WYM模型的初動到時殘差結果。從殘差分佈可見當震央距在200公里以內時，初動波 到時殘差均小於2秒。因為對初動到時影響最大的是從震源到台站之間的速度結構而 不是震源機制，此外，我們在震源機制的逆推中所依賴的主要是初動的極性和初達 波的波形。因此，殘差小於2秒的初動到時對於震源機制解不會產生影響。 在確定沿特定路徑上的一維平均模型中，對於特定的震央與測站位置，我們在 沿震央到測站的大圓路徑上選取一系列採樣點。然後確定在WYM三維模型中在每個採 樣點處的速度隨深度變化的一維模型。最後將這些採樣點處的一維模型平均，求得 相應的一維分層平均速度模型。該一維模型即代表了從震源到測站的路徑上平均的 一維速度結構。為此我們已經編寫了所需的FORTRAN程式。只要給定一對經緯度座 標，該程式就可以給出一個該兩點間的一維分層平均速度模型。圖三給出了從地震 2007111413到測站TPUB的一維平均模型的過程。 參、 理論地震圖計算 在從前一步中得到的一維分層平均速度模型中計算理論地震圖的程式。為此， 我們使用朱露培教授所編寫的在地震學界已經廣泛應用的基於頻率-波數積分法(F-K) 的FORTRAN程式(Zhu Rivera 2002 )。目前我們已經將該程式進行了必要的修改以 適用於我們的特定目的，並在自己的電腦系統上通過調試運行。使用該程式可以很 方便地計算頻率高達10Hz的理論地震圖。 肆、 建立逆推震源機制的電腦程式 25