【文章內容簡介】 對于第一個臺，（i=1）可以寫出 ()將()式中i分別代以2，3，…，n與()式相減得 （）共有n1個線性方程式。()式以消去了深度參數(shù)z，只有x、y、T、V四個未知數(shù)。由于n1一般都大于4，故()式所表示的線性方程個數(shù)多于未知數(shù)個數(shù)，即構成所謂的超定方程組…… …… …… （）為了表達簡潔起見，將此超定方程組用矩陣形式表達 ()其中，對于n個臺站的數(shù)據,其中U=V2， W=V2T。當?shù)卣鹋_站數(shù)目n小于5時，上述方程組為不定方程，有多解；大于5時，方程組為超定方程，通過變換 （A’為A的轉置）使得超定方程組變?yōu)檎ǚ匠探M，從而計算出震源參數(shù)，這就是最小二乘法原理。只要解出此方程組()，相應震源參數(shù)也可求得。即：震中位置: （x，y）介質的波速: V＝發(fā)震時刻: T=W/U.又由： 任取一臺站(i)數(shù)據代入此式便得震源深度： 當波速V已知的初定為方法因為所研究區(qū)域位于青藏高原，其地殼的厚度較大，所以近震的深度基本在上地殼的范圍之內，其速度大約為v=，在上述方法中的線性方程組就會減少一個未知數(shù)，所以（）變化之后 ， （）其中 W=V2T。只要解出此方程組()，相應震源參數(shù)也可求得。即：震中位置: （x，y）發(fā)震時刻: T=W/ V2.同理，再由求出震源的深度。這兩種方法的主要流程為： 震源計算機定位流程圖 模型反演結果 為了更好的檢驗程序，在正演t的結果上均加了6s，下面為各個反演方法所得到的結果：，兩個模型反演出的地震參數(shù)分別為： mod1:震中(,) t0= v= h= mod2:震中(,) t0= v= h=（V=）時，兩個模型反演出的地震參數(shù)分別為： mod1:震中(,) t0= h= mod2:震中(,) t0= h= 誤差分析從數(shù)據的反演來看，深度的誤差很大，這主要是因為震源深度對地震波到時的讀取誤差反應極其敏感，是地震學中最難準確測定的參數(shù)之一，各種方法對于震源深度的估計都具相當程度的不確定性，影響著人們對震源過程的認識。具體來說，震源深度的精度誤差受震中距、到時誤差和速度模型（地殼模型）三個因素制約，而且這些因素對震源深度的影響是非線性的。，震源深度的誤差隨著震中距或臺站位置的增大和走時殘差的增大而增大。，震源深度誤差隨著震中距的增大和地震波速 度的增大而增大。，走時殘差一定時，越淺的地震，定位誤差可能越大。這一結論與許多深源地震的定位結果相一致。盡管震源深度是地震學中很難精確測定的參數(shù)之一，各種方法測定震源深度的結果不同，但可以將不同的測定震源深度的結果進行對比分析，可能會改善對震源深度的測定精度。地震發(fā)生后，立即在震源區(qū)布設流動觀測地震臺站（網）是修正主震震源深度精度的有效方法。 5 實例剖析 所研究區(qū)域位于青藏高原南迦巴瓦地段。 區(qū)域臺站及示例震源分布圖1．臺網記錄的地震事件里該近震記錄根據美國國家高級地震系統(tǒng)（ANSS）提供的全球地震事件目錄，可知該事件的位置（,），發(fā)震時刻為2003年第230天的09:03,該地震事件總有38道記錄。經過上述方法反演后得到的位置為：（,），反演的結果與實際經緯度分別相差（,）。下圖分別是利用geotool軟件在選定地震事件前后，通過屏幕抓圖得到的圖像。其中，紅色的p標注的時刻是geotool軟件算出來的該地震記錄道的起跳時刻，是理論到時。黑色p標注處的時刻是經人工判斷后標記的該地震道記錄的起跳時刻，是實際到時。，a為定位前的地震記錄圖像，b為定位后的圖像。由圖可看出反演后的結果。反演后理論到時和實際的起跳時刻仍然存在幾秒的誤差，因為模型畢竟是理想狀態(tài)下的情況，而實際的資料并不是那么完美，在實際操作中很難得到極其精確的臺站到時，一方面存在人為因素，如時間拾取的精度；另一方面也存在機器設備的物理因素。而且由于地表的覆蓋層厚度，松散度等等的不同等等都會引起波至到時的誤差，所以利用地震資料到時時，還要經過判斷和人工調整。在該地震事件中，各個記錄道起跳刻的最小誤差為2s，最大的誤差為5s。桂林理工大學本科畢業(yè)設計論文 a： 近震定位前 b： 近震定位后圖 實例一近震定位前（a）和定位后（b）2．臺站網記錄的地震事件里沒有近震記錄時 當臺站網記錄的地震事件里沒有近震記錄時，利用上述反演方法進行地震的定位，可確定該地震事件發(fā)生的經度和緯度為：（,）。將得出其經度、緯度、發(fā)震時刻、深度等信息，按照標準格式寫入origin事件，利用geotool軟件對該事件就行處理后，其中最小的誤差為0. 1s。 a:近震定位前 b：近震定位后 實例二近震定位前（a）和定位后（b）為了進一步驗證此種定位的正確與否，再加一例輔證。同樣采用上述反演方法定位后，可確定該地震事件發(fā)生的位置為（,）。，其最小的誤差為0. 09s，最大的誤差為3s。 a:近震定位前 b：近震定位后