【正文】 ：發(fā)生地震的時刻，一般用字母Ο或T0表示。與震中相對的地球直徑的另一端稱為對震中；或稱震中對點。理論上將震源看成一個點,而實際上是一個區(qū)。相對于體波而言，在彈性界面附近還存在著另一類波動，從能量上而說他們只分布在彈性分界面附近，稱為面波。震源愈深，相對誤差愈??；h300km時，誤差小于其深度的5%。計算機的近震定位方法主要分為速度未知的初定為方法和速度已知的初定為方法。在我國，李善邦先生最初使用方位角和最小二乘進行地震的觀測和定位。利用地震記錄進行定位始源于歐洲和日本。如今地震發(fā)生越來越頻繁，2010年也被“尊稱”為國際地震年。桂林理工大學本科畢業(yè)設計主要表現(xiàn)為：大型建筑物破壞，普通民房破壞，山崩地裂，人畜傷亡。由于天然地震具有很大的能量，它所產(chǎn)生的地震波可以穿透很大的深度，傳播很遠的距離，而且當遇到地球深部的波阻抗差異界面時產(chǎn)生反射波，其中，天然地震的中的近震資料為研究地球深部結構提供了一個重要的信息資源和途徑。1961年，博爾特和威爾莫合作，改進了計算方案，并首先在ISS使用，隨后，國際地震中央局與美國海岸和大地測量局先后使用。其中近震的定位方法主要有石川法、和達法、高橋法、外心方位角法、假定發(fā)震時刻定位、等時量板法等十幾種方法。一般淺源地震（h100km）的震源深度誤差為深度值的10%左右。 地震波的分類在彈性介質內(nèi)傳播的縱波和橫波，由于存在整個彈性空間，因而這些波統(tǒng)稱為體波。 地震學的基本名詞和概念1．震源：地球內(nèi)部發(fā)生地震的地方稱為震源（或稱震源區(qū)）。同時，地面上受破壞最嚴重的地區(qū)叫極震區(qū)，理論上震中區(qū)和極震區(qū)是相同的，實際上由于地表局部對質條件的影響，極震區(qū)不一定是震中區(qū)。它可用線距離表示，也可用地心所張的角來表示。 地震的分類為研究方便，按震動的性質，可分為天然地震、人工地震及脈動三類。已記錄到的最深地震的震源深度約700公里。 3. 遠震：震中距大于1000公里的地震。應當指出，由于康拉德界面并不是全球普遍存在的界面，在一些地區(qū)觀測不到。3 正演計算（模型試算）測定震源位置及發(fā)震時刻最常用的方法是將非線性問題化為線性問題，然后利用最小二乘法原理迭代，全主元消元法等求解。為此，我們需要知道臺網(wǎng)各個臺站的坐標，以及至少五個臺的地震波到時。為了進行參數(shù)的求解，一般要對地殼平均速度做出假設。 基本原理 坐標變換對于近震及地方震的計算，采用平面坐標比球面坐標方便。因而，直角坐標系的原點是浮動的，每次都選來盡可能的靠近震中區(qū)，以減少變換過程中帶來的誤差。因此，必須化為經(jīng)緯度，輸出。Modd2也是如此。在此主要利用前面正演模型所得數(shù)據(jù)，作為反演之用，進而推導出震源參數(shù)。待求的震中位置為x,y,震源深度z,發(fā)震時刻為T（）；按照均勻地殼模型，可列出下列方程： 震源計算機定位原理圖 （） i=1,2,…，n.式中為地震波在地殼中的平均傳播速度，因區(qū)域不同而異。由于n1一般都大于4，故()式所表示的線性方程個數(shù)多于未知數(shù)個數(shù)，即構成所謂的超定方程組…… …… …… （）為了表達簡潔起見，將此超定方程組用矩陣形式表達 ()其中，對于n個臺站的數(shù)據(jù),其中U=V2， W=V2T。只要解出此方程組()，相應震源參數(shù)也可求得。，震源深度的誤差隨著震中距或臺站位置的增大和走時殘差的增大而增大。盡管震源深度是地震學中很難精確測定的參數(shù)之一，各種方法測定震源深度的結果不同，但可以將不同的測定震源深度的結果進行對比分析，可能會改善對震源深度的測定精度。經(jīng)過上述方法反演后得到的位置為：（,），反演的結果與實際經(jīng)緯度分別相差（,）。，a為定位前的地震記錄圖像，b為定位后的圖像。在該地震事件中，各個記錄道起跳刻的最小誤差為2s，最大的誤差為5s。 a:近震定位前 b：近震定位后 實例二近震定位前（a）和定位后（b）為了進一步驗證此種定位的正確與否，再加一例輔證。論文 6 結束語通過本次畢業(yè)設計，使我了解了近震震源參數(shù)定位的方法和原理，從正演和反演的理論知識落腳到實際資料的處理。王教授淵博的學識，豐富的實踐經(jīng)驗，高瞻遠矚、敏銳的科學眼光，將是我永遠學習的楷模；他樂觀的心態(tài)、謙遜的為人、樸實的生活態(tài)度，令我深深敬佩；他嚴謹細致、一絲不茍的作風將是我以后學習、工作中的榜樣，他循循善誘的教導和不拘一格的思路給予我無盡的啟迪。祝愿他們身體健康，工作順利，事業(yè)上取得更大成功。) filename !read the picks for all stations open(1,file=filename) read(1,*) nstn x0=0。(I2)39。()39。即do istn=2,nstn coef(istn,1)=x(istn)x(1) coef(istn,2)=y(istn)y(1) coef(istn,3)=(t(istn)+t(1))*v**2delt(istn)=(x(istn)**2x(1)**2+y(istn)**2y(1)**2)/(t(istn)**2t(1)**2)*v**2/ end dodo i=1,3 do j=1,3 solut(i,j)=0 do k=1,nstn solut(i,j)=solut(i,j)+coef(k,i)*coef(k,j) end do end do end do do i=1,3 solut(i,4)=0 do k=1,nstn solut(i,4)=solut(i,4)+coef(k,i)*delt(k) end do end do