近震定位及其應(yīng)用本科畢業(yè)論文-文庫吧資料
2025-06-25 13:05本頁面
【正文】 ,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。) x(istn)*(), y(istn)*(), t(istn)+6 end do write(2,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。) 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。(a)39?！ ∮捎诒救怂接邢?，本篇論文及程序有待進(jìn)一步改進(jìn)，請各位老師批評指正。在此一并表示衷心的感謝。謹(jǐn)致以衷心的感謝和崇高的敬意。本次畢業(yè)設(shè)計是在王教授悉心指導(dǎo)下，經(jīng)過不斷的學(xué)習(xí)和修改完成的。解決此問題的基本思想就是將非線性問題化為線性問題，然后利用最小二乘法原理迭代，全主元消元法等求解，通過反復(fù)修定，以得到震源參數(shù)的最佳值。 a:近震定位前 b：近震定位后 實例二近震定位前（a）和定位后（b）桂林理工大學(xué)本科畢業(yè)設(shè)計同樣采用上述反演方法定位后，可確定該地震事件發(fā)生的位置為（,）。將得出其經(jīng)度、緯度、發(fā)震時刻、深度等信息，按照標(biāo)準(zhǔn)格式寫入origin事件，利用geotool軟件對該事件就行處理后，其中最小的誤差為0. 1s。桂林理工大學(xué)本科畢業(yè)設(shè)計而且由于地表的覆蓋層厚度，松散度等等的不同等等都會引起波至到時的誤差，所以利用地震資料到時時，還要經(jīng)過判斷和人工調(diào)整。由圖可看出反演后的結(jié)果。黑色p標(biāo)注處的時刻是經(jīng)人工判斷后標(biāo)記的該地震道記錄的起跳時刻，是實際到時。下圖分別是利用geotool軟件在選定地震事件前后，通過屏幕抓圖得到的圖像。 區(qū)域臺站及示例震源分布圖1．臺網(wǎng)記錄的地震事件里該近震記錄根據(jù)美國國家高級地震系統(tǒng)（ANSS）提供的全球地震事件目錄，可知該事件的位置（,），發(fā)震時刻為2003年第230天的09:03,該地震事件總有38道記錄。地震發(fā)生后，立即在震源區(qū)布設(shè)流動觀測地震臺站（網(wǎng)）是修正主震震源深度精度的有效方法。這一結(jié)論與許多深源地震的定位結(jié)果相一致。，震源深度誤差隨著震中距的增大和地震波速 度的增大而增大。具體來說，震源深度的精度誤差受震中距、到時誤差和速度模型（地殼模型）三個因素制約，而且這些因素對震源深度的影響是非線性的。即：震中位置: （x，y）發(fā)震時刻: T=W/ V2.同理，再由求出震源的深度。即：震中位置: （x，y）介質(zhì)的波速: V＝發(fā)震時刻: T=W/U.又由： 任取一臺站(i)數(shù)據(jù)代入此式便得震源深度： 當(dāng)波速V已知的初定為方法因為所研究區(qū)域位于青藏高原，其地殼的厚度較大，所以近震的深度基本在上地殼的范圍之內(nèi)，其速度大約為v=，在上述方法中的線性方程組就會減少一個未知數(shù)，所以（）變化之后 ， （）其中 W=V2T。當(dāng)?shù)卣鹋_站數(shù)目n小于5時，上述方程組為不定方程，有多解；大于5時，方程組為超定方程，通過變換 （A’為A的轉(zhuǎn)置）使得超定方程組變?yōu)檎ǚ匠探M，從而計算出震源參數(shù)，這就是最小二乘法原理。()式以消去了深度參數(shù)z，只有x、y、T、V四個未知數(shù)。式（）為非線性方程組，為便于直接求解，首先應(yīng)將其線性化，變?yōu)榫€性方程組。其中任一臺站的坐標(biāo)為（xi,yi）相應(yīng)的直達(dá)波到時為Ti。 基本原理 波速V未知的初定為方法此方法通常用在臺數(shù)n4，震中距小于170km，已知直達(dá)縱波P到時，而且臺站方位和震中距分布合理的情況下。其具體方法如下。4 反演計算對于反演模型，主要目的就是利用地震波在各接收臺站的到時，結(jié)合臺站經(jīng)緯度，通過數(shù)學(xué)推導(dǎo)而確定震源參數(shù)。 正演模型建立假定某地震事件發(fā)生后被8個臺站檢測到() ，這8個臺站的經(jīng)緯度分別是：（， ）、（， ）、（， ）（， ）、（， ）、（， ）（， ）、（， ）通過給定兩個模型：Mod1的震中位置為：（， ）Mod2的震中位置為：（， ）， 臺站分布圖正演的結(jié)果：Mod1：t= 、 、Mod2：t= 、對于mod1而言，震中在其中心，它與（，）臺站相距20，有換算關(guān)系10 ，且速度v=。將直角坐標(biāo)值，化為經(jīng)緯度，的換算公式為： () 正演計算如果已知震源及臺站坐標(biāo)；設(shè)震源坐標(biāo),發(fā)生時刻，均勻介質(zhì)中波速；有n個臺站，其中任一臺站坐標(biāo)。但這樣得到的震中位置仍是局部
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1