近震定位及其應(yīng)用本科畢業(yè)論文(文件)
【正文】 任取一臺(tái)站(i)數(shù)據(jù)代入此式便得震源深度： 當(dāng)波速V已知的初定為方法因?yàn)樗芯繀^(qū)域位于青藏高原，其地殼的厚度較大，所以近震的深度基本在上地殼的范圍之內(nèi)，其速度大約為v=，在上述方法中的線性方程組就會(huì)減少一個(gè)未知數(shù)，所以（）變化之后 ， （）其中 W=V2T。具體來說，震源深度的精度誤差受震中距、到時(shí)誤差和速度模型（地殼模型）三個(gè)因素制約，而且這些因素對(duì)震源深度的影響是非線性的。這一結(jié)論與許多深源地震的定位結(jié)果相一致。 區(qū)域臺(tái)站及示例震源分布圖1．臺(tái)網(wǎng)記錄的地震事件里該近震記錄根據(jù)美國國家高級(jí)地震系統(tǒng)（ANSS）提供的全球地震事件目錄，可知該事件的位置（,），發(fā)震時(shí)刻為2003年第230天的09:03,該地震事件總有38道記錄。黑色p標(biāo)注處的時(shí)刻是經(jīng)人工判斷后標(biāo)記的該地震道記錄的起跳時(shí)刻，是實(shí)際到時(shí)。而且由于地表的覆蓋層厚度，松散度等等的不同等等都會(huì)引起波至到時(shí)的誤差，所以利用地震資料到時(shí)時(shí)，還要經(jīng)過判斷和人工調(diào)整。將得出其經(jīng)度、緯度、發(fā)震時(shí)刻、深度等信息，按照標(biāo)準(zhǔn)格式寫入origin事件，利用geotool軟件對(duì)該事件就行處理后，其中最小的誤差為0. 1s。 a:近震定位前 b：近震定位后 實(shí)例二近震定位前（a）和定位后（b）桂林理工大學(xué)本科畢業(yè)設(shè)計(jì)本次畢業(yè)設(shè)計(jì)是在王教授悉心指導(dǎo)下，經(jīng)過不斷的學(xué)習(xí)和修改完成的。在此一并表示衷心的感謝。(a)39。) write(2,39。) x(istn)*(), y(istn)*(), t(istn)+6 end do write(2,39。 print *, lai*(),phi*(),v0,t0, hend program testSUBROUTINE Gaussian(A,N1,M1) implicit none integer ::n1,m1,k1,k2,ik,i,j ,l, i1 REAL*8 :: A(n1,m1) REAL*8 :: BMAX,T,EPS EPS=0. DO K1=1,N1 BMAX=0. DO I=K1,N1 IF(BMAXABS(A(I,K1)).) THEN BMAX=ABS(A(I,K1)) L=I END IF 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 RETURNEND SUBROUTINE Gaussiansubroutine lp2xy(lai,phi,la0,ph0,xx,yy) implicit none real*8 :: ph0, la0, phi, lai, ee, a, r1,r2, r3, dl, dp, sinxcos, xx, yy ee= a= sinxcos=sin(ph0)*cos(ph0) R1=sqrt(1ee*sin(ph0)**) r2=a/R1 r3=a*(1ee)/R1**3 dl=laila0 dp=phiph0 xx=r2*dl*cos(ph0) yy=r3*dp+*r2*dl***sinxcos returnend subroutine lp2xysubroutine xy2lp(xx,yy,la0,ph0,lai,phi) implicit none real*8 :: ph0, la0, phi, lai, ee, a, r1,r2, r3, sinxcos, xx, yy ee= a= sinxcos=sin(ph0)*cos(ph0) R1=sqrt(1ee*sin(ph0)**) r2=a/R1 r3=a*(1ee)/R1**3 lai=xx/(r2*cos(ph0))+la0 phi=yy/*xx**2*tan(ph0)/(r2*r3)+ph0 returnend subroutine xy2lp當(dāng)速度給定時(shí)，在求解方程參數(shù)和未知數(shù)的量稍有變化。) la0*(),ph0*(),v,h close(2) ！stop 10 continue !(lai,phi)(x,y) do istn=1,nstn lai=x(istn) phi=y(istn) call lp2xy(lai,phi,la0,ph0,xx,yy) x(istn)=xx y(istn)=yy end do !calculating coefficents for the linear equations do istn=2,nstn coef(istn,1)=x(istn)x(1) coef(istn,2)=y(istn)y(1) coef(istn,3)=(t(is
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1