近震定位及其應(yīng)用本科畢業(yè)論文-全文預(yù)覽
【正文】 tn)**2t(1)**2)/ coef(istn,4)=t(istn)+t(1) delt(istn) =(x(istn)**2x(1)**2+y(istn)**2y(1)**2)/ end do !least squarroot method to get the solution ! dt(1) de(1,1) de(1,2) de(1,3) de(1,4) de(1,5) ! dt(2) de(2,1) de(2,2) de(2,3) de(2,4) de(2,5) ! . ! . ! . ! dt(n) de(n,1) de(n,2) de(n,3) de(n,4) de(n,5) ! de(1,1) de(2,1) ... de(n,1) ! de(1,2) de(2,2) ... de(n,2) ! de(1,3) de(2,3) ... de(n,3) ! de(1,4) de(2,4) ... de(n,4) ! de(1,5) de(2,5) ... de(n,5) do i=1,4 do j=1,4 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,4 solut(i,5)=0 do k=1,nstn solut(i,5)=solut(i,5)+coef(k,i)*delt(k) end do end do call gaussian(solut,4,5) xx=solut(1,5) yy=solut(2,5) v0=sqrt(solut(3,5)) t0=solut(4,5)/solut(3,5)h=0！print *,xx,yy do istn=1,nstnh=h+sqrt(v0*v0*(t(7)t0)*(t(7)t0)(x(7)xx)**2(y(7)yy)**2) end do h=h/nstn call xy2lp(xx,yy,la0,ph0,lai,phi) print *,39。) nstn v= h=do istn=1,nstn lai=x(istn) phi=y(istn) call lp2xy(lai,phi,la0,ph0,xx,yy) t(istn)=sqrt(xx**+yy**+h**)/v write(2,39。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?！ ∮捎诒救怂接邢蓿酒撐募俺绦蛴写M一步改進，請各位老師批評指正。謹(jǐn)致以衷心的感謝和崇高的敬意。解決此問題的基本思想就是將非線性問題化為線性問題，然后利用最小二乘法原理迭代，全主元消元法等求解，通過反復(fù)修定，以得到震源參數(shù)的最佳值。同樣采用上述反演方法定位后，可確定該地震事件發(fā)生的位置為（,）。桂林理工大學(xué)本科畢業(yè)設(shè)計由圖可看出反演后的結(jié)果。下圖分別是利用geotool軟件在選定地震事件前后，通過屏幕抓圖得到的圖像。地震發(fā)生后，立即在震源區(qū)布設(shè)流動觀測地震臺站（網(wǎng)）是修正主震震源深度精度的有效方法。，震源深度誤差隨著震中距的增大和地震波速 度的增大而增大。即：震中位置: （x，y）發(fā)震時刻: T=W/ V2.同理，再由求出震源的深度。當(dāng)?shù)卣鹋_站數(shù)目n小于5時，上述方程組為不定方程，有多解；大于5時，方程組為超定方程，通過變換 （A’為A的轉(zhuǎn)置）使得超定方程組變?yōu)檎ǚ匠探M，從而計算出震源參數(shù)，這就是最小二乘法原理。式（）為非線性方程組，為便于直接求解，首先應(yīng)將其線性化，變?yōu)榫€性方程組。 基本原理 波速V未知的初定為方法此方法通常用在臺數(shù)n4，震中距小于170km，已知直達縱波P到時，而且臺站方位和震中距分布合理的情況下。4 反演計算對于反演模型，主要目的就是利用地震波在各接收臺站的到時，結(jié)合臺站經(jīng)緯度，通過數(shù)學(xué)推導(dǎo)而確定震源參數(shù)。將直角坐標(biāo)值，化為經(jīng)緯度，的換算公式為： () 正演計算如果已知震源及臺站坐標(biāo)；設(shè)震源坐標(biāo),發(fā)生時刻，均勻介質(zhì)中波速；有n個臺站，其中任一臺站坐標(biāo)。在計算程序中，坐標(biāo)原點確定在接收到該次地震信號各臺站經(jīng)緯度的平均值0及0上，x軸向東，y軸向北。為此需要將經(jīng)緯度換算成局部的直角坐標(biāo)系。為了對此問題進行研究，有必要建立近震定位的正反演模型。由臺網(wǎng)提供的到時數(shù)據(jù)主要是指P 波的到時，因為各種續(xù)至波的相位通常是難以識別的。通過反復(fù)修定，以得到震源參數(shù)的最佳值。3. Pn/Sn震相Pn/Sn分別為地震波在地殼底部的莫霍界面產(chǎn)生的折射縱、橫波，~~。 近震的主要震相由于近震的震中距很小，其地震波的傳播主要局限于地殼及莫霍面附近，主要震相有Pg/Sg（）、P*/S*、Pn/Sn、P*P/S*S、P*S/S*P、PmP/SmS及PmS/SmP等。有時也將中源地震和深源地震統(tǒng)稱為深震。對于天然地震，有下述分類：（一）按震源深度分：1．淺源地震：震源深度小于60公里的天然地震稱為淺震；也稱正常深度地震。5．發(fā)震時刻
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1