【正文】 法是將非線性問題化為線性問題，然后利用最小二乘法原理迭代，全主元消元法等求解。 圖 近震震相及其射線路徑示意圖 （ α c、 β c分別為康拉德界面及莫霍面的臨界角） 1. Pg/Sg 震相 由于地球淺部存在強 烈的橫向及縱向上的速度變化，尤其是有沉積該層地區(qū)的垂向梯度的影響，地震波射線在地殼淺部 10 公里的范圍內(nèi)發(fā)生廻折而到達地表，并在震中附近數(shù)公里的范圍，這種波束與直達波，被稱之為廻折或潛波（ Diving Wave），對于縱、橫波來說，分別命名為 Pg和 Sg。 2．中源地震：震源深度在 60 公里至 300 公里之間的地震稱為中源地震。一般用字母 ⊿ 表示。為便于對地震的了解研究，現(xiàn)將其基本知識介紹如下。一般淺源地震（ h100km）的震源深度誤差為深度值的 10%左右。 1961年，博爾特和威爾莫合作，改進了計算方案，并首先在 ISS使用，隨后，國際地震中央局與美國海岸和大地測量局先后使用。主要表現(xiàn)為：大型建筑物破壞，普通民房破壞，山崩地裂，人畜傷亡。一次強烈地震的發(fā)生，常伴隨著地面變形和地層錯動 ，其破壞力是相當(dāng)大的?！秶H地震匯編》（ ISS）最先采用最小二乘法進行計算修訂震中。如果有合理的臺網(wǎng)分布和適用的走時表，對于Δ 100km的區(qū)域性事件，定位精度為 ? 1km左右；對于 1001000km的近震，可達 ? 2至5km；遠震為 ? 5— 10km。面波又分為瑞雷波和勒夫波 (如右圖 所示 )。 4．震中距離：在地面上，從震中到任一點沿大圓弧測量的距離稱為震中距離。論文 3 地震，有下述分類： （一）按震源深度分： 1．淺源地震：震源深度小于 60 公里的天然地震稱為淺震；也稱正常深度地震。 近震的主要震相 由于近震的震中距很小，其地震波的傳播主要局限于地殼及莫霍面附近，主要震相有 Pg/Sg（ SP/ ）、 P*/S*、 Pn/Sn、 P*P/S*S、 P*S/S*P、 PmP/SmS 及 PmS/SmP 等，如圖 所示。 桂林理工大學(xué)本科畢業(yè)設(shè)計如果地震發(fā)生在臺網(wǎng)以外，這種答案是不確定的。但是，在作具體計算時，只用到其中一部分臺站，也就是只需要把這部分臺站的經(jīng)緯度（球面坐標(biāo)）化為平面坐標(biāo)。 圖 正演計算示意圖 桂林理工大學(xué)本科畢業(yè)設(shè)計在此主要利用前面正演 模型所得數(shù)據(jù)，作為反演之用，進而推導(dǎo)出震源參數(shù)。 ()式以消去了深度參數(shù) z，只有 x、 y、 T、 V 四個未知數(shù)。 只要解出此方程組 ()，相應(yīng)震源參數(shù)也可求得。 ，走時殘差一定時，越淺的地 震，定位誤差可能越大。下圖分別是利用 geotool 軟件在選定地震事件前后，通過屏幕抓圖得到的圖像。桂林理工大學(xué)本科畢業(yè)設(shè)計 a:近震定位前 b：近震定位后 圖 實例二近震定位前（ a）和定位后（ b） 桂林理工大學(xué)本科畢業(yè)設(shè)計在此一并表示衷心的感謝。y0=0 do istn=1,nstn read(1,*) x(istn),y(istn), t(istn) x(istn)=x(istn)*atan()/45 y(istn)=y(istn)*atan()/45 x0=x0+x(istn) y0=y0+y(istn) end do la0=x0/nstn ph0=y0/nstn close(1) goto 10 !calculating sythetic data la0=99*atan()/45 ph0=27*atan()/45 open(2,file=39。()39。論文 25 end subroutine xy2lp 當(dāng)速度給定時，在求解方程參數(shù)和未知數(shù)的量稍有變化。lai phi v0 t0 h39。) nstn 桂林理工大學(xué)本科畢業(yè)設(shè)計論文 18 參考文獻 [1] Seth Stein,Michael Wysession. An Introduction to Seismology Earthquakes,and Earth Publishing, 2020. [2] 李善邦，中國地震 . 地震出版社 ,1981. [3] 單娜琳 ,程志平 ,劉云禎 .工程地震勘探 . 冶金工業(yè)出版社， 2020. [4] 姜枚 ,王有學(xué) ,錢輝 .造山的高原 — 青藏高原及其鄰區(qū)的寬頻地震探測與地殼上地幔結(jié)構(gòu) .地質(zhì)出版社， 2020. [5] 熊章強 ,方根顯 .淺層地震勘探 .地震出版社， 2020. [6] 大港油田科技叢書編委會 .地震勘探資料處理和解釋技術(shù) .石油工業(yè)出版社， 1999. [7] 時振梁 ,張少泉 ,趙榮國等 .地震工作手冊 .地震出版社， 1990. [8] 傅淑芳 ,劉寶城 ,李文藝 .地震學(xué)教程（上、下冊） . 地震出版社， 1980. [9] 彭國倫 .Fortran 95 程序設(shè)計 .中國電力出版社， 2020. 桂林理工大學(xué)本科畢業(yè)設(shè)計 本次畢業(yè)設(shè)計 是在 王 教授悉心指導(dǎo)下，經(jīng)過不斷的學(xué)習(xí)和修改完成的。 a:近震定位前 b：近震定位后 圖 實例二近震定位前（ a）和定位后（ b） 桂林理工大學(xué)本科畢業(yè)設(shè)計由圖 可看出反演后的結(jié)果 。 桂林理工大學(xué)本科畢業(yè)設(shè)計論文 12 的估計都具相當(dāng)程度的不確定性，影響著人們對震源過程的認識。39。待求的震中位置為 x,y,震源深度 z,發(fā)震時刻為 T（圖 ）；按照均勻地殼模型，可列出下列方程： 圖 震源計算機定位原理圖 22222 )()()( iii TTVzyyxx ?????? （ ） i=1,2,?， n. 式中 V 為地震波在地殼中的平均傳播速度，因區(qū)域不同而異。 桂林理工大學(xué)本科畢業(yè)設(shè)計論文 6 于是，地面上任一點（其經(jīng)緯度為 ? 、 ? ）相對于坐標(biāo)原點 O（ ? 0,? 0）的直角坐標(biāo) x 、 y 的換算公式為： ?????????RnxRyRxmn2t an)(c o s)(02000?????? （ ） 其中參數(shù) : 2/1022 )s in1( ?e aR n ?? 2/30222)s in1( )1( ?e eaR m ? ?? 而參數(shù) nR 和 mR 中的常數(shù)又為： 2e =222aba ? = (地球偏心率 ) a = km (地球長半徑 ) b = (地球短半徑 ) 經(jīng)換算為平面坐