【正文】 論文 23 END DO IF() STOP 4444 IF() THEN DO J=K1,M1 T=A(L,J) A(L,J)=A(K1,J) A(K1,J)=T END DO END IF T=1./A(K1,K1) K2=K1+1 DO J=K2,M1 A(K1,J)=A(K1,J)*T DO I=K2,N1 A(I,J)=A(I,J)A(I,K1)*A(K1,J) END DO END DO END DO DO IK=2,N1 I=M1IK I1=I+1 DO J=I1,N1 A(I,M1)=A(I,M1)A(I,J)*A(J,M1) END DO END DO RETURN END SUBROUTINE Gaussian subroutine lp2xy(lai,phi,la0,ph0,xx,yy) implicit none 桂林理工大學(xué)本科畢業(yè)設(shè)計(jì)()39。(a)39。謹(jǐn)致以衷心的感謝和崇高的敬意。同樣采用上述反演方法定位后，可確定該地震事件發(fā)生的位置為（ ,）。而且由于地表的覆蓋層厚度，松散度等等的不同等等都會(huì)引起波至到時(shí)的誤差，所以利用地震資料到時(shí)時(shí)，還要經(jīng)過判斷和人工調(diào)整。 區(qū)域臺(tái)站及示例震源分布圖 1．臺(tái)網(wǎng)記錄的地震事件里該近震記錄 根據(jù)美國國 家高級(jí)地震系統(tǒng)（ ANSS）提供的全球地震事件目錄，可知該事件的位置（ ,），發(fā)震時(shí)刻為 2020 年第 230 天的 09:03,該地震事件總有38 道記錄，震中距范圍在 。 ，震源深度的誤差隨著震中距或臺(tái)站位置的增大和走時(shí)殘差的增大而增大。 只要解出此方程組 ()，相應(yīng)震源參數(shù)也可求得。為此將（ ） 式展開得 )2( 22 22222222TTTTV zyyyyxxxx ii iii ??? ?????? () i=1,2,? ,n. 對(duì)于第一個(gè)臺(tái)，（ i=1） 可以寫出 )2( 22 2121222112212TTTTV zyyyyxxxx ??? ?????? () 將 ()式中 i分別代以 2， 3，?， n 與 ()式相減得 桂林理工大學(xué)本科畢業(yè)設(shè)計(jì)實(shí)際計(jì)算表明，只要輸入數(shù)據(jù)正確，初定的震源位置及發(fā)震時(shí)刻足夠精確，就可供地震速報(bào)之用。因此，必須化為經(jīng)緯度 ? ， ? 輸出。為此需要 將經(jīng)緯度換算成局部的直角坐標(biāo)系。由臺(tái)網(wǎng)提供的到時(shí)數(shù)據(jù)主要是指 P 波的到時(shí)，因?yàn)楦鞣N續(xù)至波的相位通常是難以識(shí)別的。一般來說， Pn 震相在莫霍面的臨界反射以遠(yuǎn)地區(qū)總是已初至波的形式出現(xiàn)，與其它折射波一樣，波至呈現(xiàn)為線性分布。 2. 近震：震中距小于 1000 公里的地震。 地震的分類 為研究方便，按震動(dòng)的性質(zhì)，可分為天然地震 、人工地震及脈動(dòng)三類。同時(shí)，地面上受破壞最嚴(yán)重的地區(qū)叫極震區(qū)，理論上震中區(qū)和極震區(qū)是相同的，實(shí)際上由于地表局部對(duì)質(zhì)條件的影響，極震區(qū)不一定是震中區(qū)。 在彈性介質(zhì)內(nèi)傳播的縱波和橫波，由于存在整個(gè)彈性空間，因而這些波統(tǒng)稱為體波。計(jì)算機(jī)的近震定位方法主要分為速度未知的初定為方法和速度已知的初定為方法。 利用地震記錄進(jìn)行定位始源于歐洲和日本。論文 目 次 摘要 ?????????????????????????????????Ⅰ Abstract??????????????????????????????? Ⅱ 1 前言 ??????????????????????????????? 1 2 近震的理論基礎(chǔ) ??????????????????????????? 2 地震學(xué)的基本名詞和概念???????????????????? 2 地震的分類?????????????????????????? 2 近震的主要震相 ???????????????????????? 3 3 正演計(jì)算（模型試算） ??????????????????????? 5 基本原理??????????????????????????? 5 坐標(biāo)變換???????????????????????? 5 正演計(jì)算???????????????????????? 6 正演模型的建立???????????????????????? 7 4 反演計(jì)算 ????????????????????????? ????? 8 基本原理??????????????????????????? 8 當(dāng)速度 V 未知的初定方法????????????????? 8 當(dāng)速度 V 已知的初定方法????????????????? 10 模型反演結(jié)果????????????????????????? 11 誤差分析??????????????????????????? 11 5 實(shí)例剖析 ?????????????????????????????? 13 6 結(jié)束語 ??????????????????????????????? 17 參考文獻(xiàn) ??????????????????????????????? 18 附錄 近震定位程序 ?????????????????????????? ? 19 桂林理工大學(xué)本科畢業(yè)設(shè)計(jì)任何事情都具有雙面性，地震對(duì)我們的生命安全造成了很大的威脅的同時(shí)也為專家學(xué)者們的研究地球提供了豐富的資料。 1953年，我國開始采用大量觀測(cè)數(shù)據(jù)修訂震中。發(fā)震時(shí)刻的誤差為 1/10— 1/2s. 桂林理工大學(xué)本科畢業(yè)設(shè)計(jì)理論上將震源看成一個(gè)點(diǎn) ,而實(shí)際上是一個(gè)區(qū)。 5．發(fā)震時(shí)刻：發(fā)生地震的時(shí)刻，一般用字母 Ο 或 T0表示。已記錄到的最深地震 的震源深度約 700 公里。論文 4 2. P*/S*震相 P*/S*分別為地震波在地殼內(nèi)部的康拉德界面產(chǎn)生的折射縱、橫波，它們的速度分別為 ~~。此方法對(duì)近震經(jīng)緯度的確定比較準(zhǔn)確的。對(duì)于正演模型，其 基本思路是通過設(shè)定震源初始參數(shù)（經(jīng)度、緯度、高程、發(fā)震時(shí)刻及均勻介質(zhì)中的波速）計(jì)算出地震波從震源到各臺(tái)站的走時(shí)，進(jìn)而求出各臺(tái)站接收到地震波的時(shí)刻；這些數(shù)據(jù)就作為反演的數(shù)據(jù)來源。 圖 坐標(biāo)轉(zhuǎn)換示意圖 桂林理工大學(xué)本科畢業(yè)設(shè)計(jì) Modd2 也是如此。其中任一臺(tái)站的坐標(biāo)為（ xi,yi） 相應(yīng)的直達(dá)波到時(shí)為 Ti。論文 10 BAX? ? BAAXA 39。論文 11 這兩種方法的主要流程為： 臺(tái) 站 經(jīng) 緯 度 轉(zhuǎn) 成 直 角 坐 標(biāo) 值計(jì)算方程組系數(shù)矩陣 和常數(shù)矩陣 AB求解方程組 A X = B依方程組的解進(jìn)一步求得震源參數(shù) ( x , y , z , T )作出震源與臺(tái)站相對(duì)位置圖形直角坐標(biāo)值轉(zhuǎn)回經(jīng)緯度表示輸出 結(jié)束 開始 圖 震源計(jì)算機(jī)定位流程圖 模型反演結(jié)果 為了更好的檢驗(yàn)程序，在正演 t 的結(jié)果上均加了 6s，下面為各個(gè)反演方法所得到的結(jié)果： V未知 時(shí)，兩個(gè)模型反演出的地震參數(shù)分別為： mod1:震中 (,) t0= v= h= mod2:震中 (,) t0= v= h= V作為已知條件（ V=）時(shí)，兩個(gè)模型反演出