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