• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    考慮到時誤差的地震定位算法及其在四川地區(qū)2001—2008年地震定位的應(yīng)用

    2012-12-06 09:13:56萬永革盛書中程萬正張致偉曄趙曉燕卜玉菲薛志芳劉金麗
    地震地質(zhì) 2012年1期
    關(guān)鍵詞:臺網(wǎng)臺站震源

    萬永革 盛書中 程萬正 張致偉 武 曄趙曉燕 卜玉菲 薛志芳 劉金麗

    1)防災(zāi)科技學(xué)院,河北三河燕郊開發(fā)區(qū) 065201

    2)四川省地震局,成都 610041

    3)甘肅省地震局,蘭州 730000

    4)河北省地震局秦皇島中心臺,秦皇島 066100

    5)河北省曲周縣實驗中學(xué),邯鄲 057250

    考慮到時誤差的地震定位算法及其在四川地區(qū)2001—2008年地震定位的應(yīng)用

    萬永革1)盛書中1)程萬正2)張致偉2)武 曄1)趙曉燕1)卜玉菲3)薛志芳4)劉金麗5)

    1)防災(zāi)科技學(xué)院,河北三河燕郊開發(fā)區(qū) 065201

    2)四川省地震局,成都 610041

    3)甘肅省地震局,蘭州 730000

    4)河北省地震局秦皇島中心臺,秦皇島 066100

    5)河北省曲周縣實驗中學(xué),邯鄲 057250

    精確的地震位置對于地震活動性、地震層析成像和地殼應(yīng)力場反演具有相當(dāng)重要的意義,對于地震速報也具有重要的應(yīng)用價值。將觀測到時的不確定性、臺站高程、地震震源深度進(jìn)行約束的同時,根據(jù)反演理論給出了地震震源位置精確估計和誤差估計的方法。該算法聯(lián)合考慮Pg波、Sg波、Pn波和Sn波的到時進(jìn)行反演,數(shù)據(jù)量的增加可以增強地震位置的準(zhǔn)確性,并可同時應(yīng)用于地方震和區(qū)域地震。采用模擬數(shù)據(jù)對該地震定位算法進(jìn)行檢驗發(fā)現(xiàn),該算法在觀測數(shù)據(jù)的不確定性不等時明顯優(yōu)于其他方法。將該算法應(yīng)用于四川地區(qū)2001—2008年間的地震定位,得到的地震位置更加符合地震的叢集性并集中于斷裂帶附近。這些結(jié)果為四川地區(qū)的地震活動性、斷層構(gòu)造以及地震層析成像研究打下了基礎(chǔ),并且為汶川地震之前的地震活動前兆研究也提供了有益幫助。

    地震定位 不確定性 臺站高程 四川地區(qū)

    0 引言

    地震定位是地震學(xué)中最基本的問題之一,對于研究地震活動構(gòu)造、地球內(nèi)部結(jié)構(gòu)、震源的幾何構(gòu)造等地震學(xué)基本問題都有重要意義。另外,基于快速準(zhǔn)確的地震定位的地震速報,對于震后的減災(zāi)、救災(zāi)工作也具有重要的實際意義。因此,地震學(xué)家一直在不斷改進(jìn)或提出新的定位方法。

    地震定位算法有很多種(Lee et al.,1975;趙仲和,1983)。一些研究者采用S波與P波的到時差和虛波速度確定地震的位置,如Inglada法。該方法的優(yōu)點在于不需要給定初始震中位置,但采用了與實際相差較大的均勻地殼模型(徐果明等,1982)。最常用的地震定位算法為Geiger法(Geiger,1912)。這種方法可采用較為復(fù)雜的地殼模型,甚至三維不均勻速度結(jié)構(gòu)(Wu et al.,2008),其好處是能夠根據(jù)與實際相近的地殼模型反演震源位置,但必須給出與真實震源位置相差不大的初始震源位置和發(fā)震時刻。還有的研究者是采用遺傳算法(萬永革等,1995,1997)和單純形法(趙珠等,1994),這類方法的好處是不需要求理論走時對震源參數(shù)的偏導(dǎo)數(shù),缺點是搜索較慢,且難以估計參數(shù)的不確定性。還有一些研究者采用相對定位法,如主事件定位法(Spence et al.,1980;周仕勇等,1999)、雙差定位法(Waldhauser et al.,2000)。主事件定位法必須給出較為準(zhǔn)確的主事件位置,并且其他事件與主事件的空間距離要比震源到臺站的距離小很多。雙差定位法可以得到事件之間的相對位置,其絕對位置的準(zhǔn)確性也能得到保證,但該方法不能對零散地震進(jìn)行定位。

    牟磊育等(2006)研究了原始的Geiger法,給出了聯(lián)合采用Inglada法和Geiger法對區(qū)域臺網(wǎng)地震到時數(shù)據(jù)的精確定位,編寫了MATLAB程序。雖然該方法是一種較為理想的區(qū)域臺網(wǎng)定位法,但該程序沒有利用Pn波到時,沒有考慮各種震相到時拾取準(zhǔn)確性的差別,并且不能給出定位參數(shù)的不確定范圍。本文在他們工作的基礎(chǔ)上增加了Pn波到時的利用,并且考慮了不同震相到時拾取的不同誤差,利用反演理論給出了估計震源參數(shù)的標(biāo)準(zhǔn)差。將我們研制的方法應(yīng)用于四川地區(qū)2001—2008年的地震定位,得到了較為精確的地震位置。

    1 定位算法及其改進(jìn)

    1.1 考慮觀測誤差的Geiger法反演地震震源參數(shù)

    4個參數(shù)描述了發(fā)震時間和震源位置,我們把這些參數(shù)稱為模型,定義模型矢量

    采用Geiger法是通過考慮目標(biāo)位置的擾動把問題線性化

    這里m0是猜想的最佳位置,可采用Inglada法求出粗略的震源發(fā)震時間和位置(徐果明等,1982)。m是離m0很小距離的新的位置??梢杂锰├占墧?shù)展開式的第1項來近似表示在m處所預(yù)測的時間

    為了使這些殘差最小,我們試圖用下列方式求解模型的修正量Δm

    這里G 是偏導(dǎo)數(shù)Gij=?ti/?mj的矩陣,i=1,2,…,n,j=1,…,4。 為了得到位置的調(diào)整值 Δm,用標(biāo)準(zhǔn)的最小二乘法得到方程(4)的最佳擬合。下一步,我們使m0處于m0+Δm的位置,重復(fù)這一過程,直到定位收斂。假如初始猜想的位置離實際位置不是太遠(yuǎn),迭代過程一般能相當(dāng)快地收斂。

    考慮觀測誤差時,式(6)可以寫為

    其中,δd為觀測誤差,對于Pg波,Sg波和Pn波可以有不同的誤差,Cd為觀測誤差的協(xié)方差矩陣。利用最大似然解估計公式可得(Jackson et al.,1985;Shen et al.,1996,萬永革等,2004,2008a,b)

    解的協(xié)方差矩陣可表示為

    該矩陣對角線元素的開方即為解的標(biāo)準(zhǔn)差。

    平均殘差定義為

    其中m為數(shù)據(jù)個數(shù),即方程的個數(shù)。

    1.2 Pn,Sn波走時及其偏導(dǎo)數(shù)

    設(shè)地球平層模型的地殼層的序列編號為l,最下一層為地幔層m,地殼層厚度為dl,速度為Vl,入射角為θl,地幔層的速度為Vm。震源位于j層中,dz為震源層中震源距該層頂部的距離。根據(jù)該地球模型,只有震源層的入射角滿足時才會出現(xiàn)首波,在其他層中的入射角為

    臨界觀測到首波的距離為

    震中距為Δ的首波的走時為

    Pn波走時的偏導(dǎo)數(shù)

    1.3 臺站高程的校正

    在平原或高原地區(qū),臺站高程相差不大,對定位沒有太大影響。但如果在地形起伏比較大的地區(qū),臺站高程往往相差較大,如果不考慮臺站高程會導(dǎo)致定位結(jié)果誤差較大。牟磊育等(2006)的程序因用于平原地區(qū)而沒有采用臺站高程校正。本研究在計算某一臺站走時時,采用地殼模型第1層厚度增加臺站高程值的方法來校正臺站高程的影響,使得定位結(jié)果更為準(zhǔn)確。

    1.4 地震震源深度的約束

    其他定位方法一般沒有對深度進(jìn)行約束。只是限定如果深度不滿足一定的范圍(如空中地震或深度特別大地震),則強制給出某一深度值。這樣處理的結(jié)果會導(dǎo)致震源參數(shù)不滿足最小殘差的要求。我們在每一次迭代過程中判斷修正后的震源深度是否在合理的范圍之內(nèi),如果不合理,只將參數(shù)增加到合理深度的邊界,進(jìn)行下一步迭代。下一步迭代還有可能使得震源深度向較為合理的方向進(jìn)展。最終得到殘差最小且深度在合理范圍內(nèi)的震源參數(shù)。

    1.5 地震深度的精確求解

    由于地震震相數(shù)據(jù)觀測數(shù)據(jù)有限,并且地震臺站大多位于同一水平面上,對深度的分辨率往往不高。為了使深度有更好的約束,我們分別在初始水平位置的基礎(chǔ)上設(shè)置每一層的中點作為深度的初始深度進(jìn)行迭代求解,選擇殘差最小的解為最優(yōu)解。

    2 數(shù)值實驗

    為了測試以上的定位方法,假設(shè)一個由9個臺站組成的虛擬臺網(wǎng),虛擬臺網(wǎng)之形狀為邊長300km的正方形(圖1)。虛擬臺網(wǎng)下的速度模型按照地球物理較為全面的觀測結(jié)果得到的Crust5.1(Mooney et al.,1998)給定。具體參數(shù)見表 1。

    圖1 虛擬臺站及虛擬地震位置Fig.1 Locations of the virtual stations and earthquake.臺站的編號同表1

    假設(shè)在虛擬臺網(wǎng)內(nèi)發(fā)生一個地震,發(fā)震時間 0時刻,位于橫軸 75.071km,縱軸50km處,深度為19km(圖1)。為簡單起見,虛擬臺網(wǎng)的高程全部設(shè)為0,根據(jù)虛擬臺網(wǎng)內(nèi)臺站的位置坐標(biāo)、速度模型、震中位置和發(fā)震時間,計算出虛擬臺網(wǎng)內(nèi)各臺站的Pg、Sg和Pn的到時,分別隨機產(chǎn)生以0.5,2.5,1.5幅度的隨機干擾,得到模擬 Pg、Sg和 Pn的到時(表2),然后采用我們修訂的程序進(jìn)行定位計算。初始位置X,Y,Z,t分別設(shè)定為50,30,18,0。當(dāng)采用相同的到時拾取誤差時,得到的地震震源 X,Y,Z,t分別為 76.0866,50.1371,18.1681,0.1544;采用與虛擬隨機擾動相同的誤差進(jìn)行反演得到的 X,Y,Z,t分別為 75.6551km、50.1251km、17.6933km、0.1366s,得到的該參數(shù)的標(biāo)準(zhǔn)差為 1.4504km、1.5413km、4.2750km 和0.1958s。采用不同的到時拾取誤差得到的結(jié)果更加接近虛擬設(shè)置的地震位置,可見引用相對到時標(biāo)準(zhǔn)差的重要性。

    3 2001—2008年的四川地區(qū)地震定位實例

    采用我們的程序?qū)λ拇ǖ貐^(qū)2001—2008年的地震進(jìn)行了研究,對震相進(jìn)行核對,得到1,831個地震的位置(圖2)。原來所有地震震相的平均殘差減少為原來的0.35%。定位后的位置更加集中在斷裂帶上,并具有更強的叢集性,反映了小震分布的基本特征。重新定位后的小震叢集區(qū)主要為龍門山斷裂帶、鮮水河斷裂帶及其西南區(qū)域、馬邊-鹽津斷裂帶附近及華鎣山斷裂帶北部和東南區(qū)域。很多研究者采用雙差定位法來確定地震的叢集特性(Waldhauser et al.,2000),其重定位是將地震組成地震簇而精確確定地震之前的相對位置,由于去除了不能組成地震簇的離散地震,因此能反映地震的叢集特性從而確定斷層的形狀(萬永革等,2008a;王福昌等,2008)。而本研究沒有去除離散地震,也明顯表現(xiàn)出地震的叢集特性,說明了研究結(jié)果的客觀性。趙珠等(1997)采用該地區(qū)早期較長時間的地震資料確定的地震位置也表現(xiàn)為地震叢集在龍門山斷裂帶附近的特性,與本文的研究結(jié)果有類似之處,表明了本文處理方法的正確性。

    在深度上,地震目錄列出的某些震源深度達(dá)60km左右(已達(dá)不易發(fā)生地震的上地幔),重新定位后的震源深度主要集中在8km以上的上地殼(圖3),下地殼的地震較少。這與Takagi等(1997)的結(jié)果類似。另外,由于這期間大部分為小震,均發(fā)生在淺部的4~8km的范圍內(nèi),深部缺少地震,似乎可以稱之為深部“地震空區(qū)”,這也許是地震預(yù)測研究人員要注意的現(xiàn)象。

    4 結(jié)論與討論

    與其他地震定位方法對比,我們的方法有如下優(yōu)勢:1)可以利用Pn波到時數(shù)據(jù);2)可以對數(shù)據(jù)給予不同的觀測誤差,從而確定地震定位的不確定性;3)考慮了臺站高程;4)可以約束震源深度的搜索范圍;5)通過對各個層中分別設(shè)置震源來去除由于地殼結(jié)構(gòu)各層速度差別較大而造成非最優(yōu)解的情況。實驗結(jié)果表明,在數(shù)據(jù)具有較大誤差的情況下,我們的程序仍然可以得到較為精確的震源位置。將該算法用于2001年至汶川地震之前的四川地區(qū)地震記錄,得到了殘差更小的震源位置。

    圖2 定位的地震震源位置Fig.2 The hypocentral locations of earthquakes located.

    將該算法用于2001年至汶川地震之前的四川地區(qū)1831個小震震源精確定位,結(jié)果更為精確地顯示了地震的叢集性和在斷裂帶附近發(fā)生的特性,并且震源深度集中在淺部4~8km的范圍內(nèi)。這些結(jié)果為進(jìn)一步研究汶川地震之前的地震活動特征及分析地震前兆現(xiàn)象打下了基礎(chǔ)。

    本文提出了考慮P波、S波和Pn波到時拾取誤差來進(jìn)行地震定位的方法。實際上地震波到時的讀取和地震定位由不同的人員進(jìn)行,僅根據(jù)地震觀測報告無法得到確切的到時讀取的誤差。因此,本文采用直達(dá)P波到時的精度為0.5、S波到時拾取誤差為1.5s、Pn波到時的誤差為1.0s的原則來進(jìn)行定位。更為嚴(yán)格的地震定位需要給定不同臺站不同拾取震相不同到時的誤差,這樣才能更好地利用到時數(shù)據(jù)的拾取誤差。因此,為了更好確定地震位置,建議在地震觀測報告中增加震相拾取的誤差參數(shù)。

    圖3 1831個地震的深度分布Fig.3 The depth distribution of the 1831 earthquakes.

    我們編寫的MATLAB程序可向感興趣的讀者提供。

    致謝 審稿人對本文提出了建設(shè)性修改意見,特此致謝。

    牟磊育,趙仲和,張偉,等.2006.用INGLADA與GEIGER方法實現(xiàn)近震精定位[J].中國地震,22(3):294—302.

    MOU Lei-yu,ZHAO Zhong-he,ZHANGWei,et al.2006.Using the INGLADA and GEIGER method to realize regional earthquake location accurately[J].Earthquake Research in China,22(3):294—302(in Chinese).

    萬永革,李鴻吉.1995.遺傳算法在確定震源位置中的應(yīng)用[J].地震地磁觀測與研究,16(6):1—7.

    WAN Yong-ge,LI Hong-ji.1995.The preliminary study on the seismic hypocenter location using genetic algorithms[J].Seismological and Geomagnetic Observation and Research,16(6):1—7(in Chinese).

    萬永革,李清河,李鴻吉,等.1997.用遺傳算法確定三維橫向不均勻介質(zhì)中的近震震源位置[J].西北地震學(xué)報,19(2):7—14.

    WAN Yong-ge,LIQing-he,LIHong-ji,et al.1997.Hypocentral locations of the near earthquakes in 3-D lateral heterogeneousmedium determined by the genetic algorithms[J].Northwestern Seismological Journal,19(2):7—14(in Chinese).

    萬永革,王敏,沈正康,等.2004.利用GPS和水準(zhǔn)測量資料反演2001年昆侖山口西8.1級地震的同震滑動分布[J].地震地質(zhì),26(3):393—404.

    WAN Yong-ge,WANG Min,SHEN Zheng-kang,et al.2004.Co-seismic slip distribution of the 2001 west of Kunlun Mountain Pass earthquake inverted by GPSand leveling data[J].Seismology and Geology,26(3):393—404(in Chinese).

    萬永革,沈正康,刁桂苓,等.2008a.利用小震分布和區(qū)域應(yīng)力場確定大震斷層面參數(shù)方法及其在唐山地震序列中的應(yīng)用[J].地球物理學(xué)報,51(3):793—804.

    WAN Yong-ge,SHEN Zheng-kang,DIAO Gui-ling,et al.2008a.An algorithm of fault parameter determination using distribution of small earthquakes and parameters of regional stress field and its application to Tangshan earthquake sequence[J].Chinese JGeophys,51(3):569—583(in Chinese).

    萬永革,沈正康,王敏,等.2008b.根據(jù)GPS和InSAR數(shù)據(jù)反演2001年昆侖山口西地震同震破裂分布[J].地球物理學(xué)報,51(4):1074—1084.

    WAN Yong-ge,SHEN Zheng-kang,WANG Min,et al.,2008b.Coseismic slip distribution of the 2001 West Kunlun Mountain Pass earthquake constrained using GPSand InSAR data[J].Chinese JGeophysics,51(4):753—764.

    萬永革,沈正康,盛書中,等.2009.2008年汶川大地震對周圍斷層的影響[J].地震學(xué)報,31(2):128—139.

    WAN Yong-ge,SHEN Zheng-kang,SHENG Shu-zhong,etal.2009.The influence of2008Wenchuan earthquake on surrounding faults[J].Acta Seismologica Sinica,31(2):128—139(in Chinese).

    王福昌,萬永革,胡順田.2008.粒子群算法在主震斷層面參數(shù)估計中的應(yīng)用[J].地震研究,31(2):149—154.

    WANG Fu-chang,WAN Yong-ge,HU Shun-tian.2008.Application of particle swarm optimization to the estimation of mainshock fault plane parameters[J].Journal of Seismological Research,31(2):149—154(in Chinese).

    徐果明,周惠蘭.1982.地震學(xué)原理[M].北京:科學(xué)出版社.

    XU Guo-ming,ZHOU Hui-lan.1982.Principle of Seismology[M].Science Press,Beijing(in Chinese).

    趙珠,丁志峰,易桂喜,等.1994.西藏地震定位:一種使用單純形優(yōu)化的非線性方法[J].地震學(xué)報,16(2):212—219.

    ZHAO Zhu,DING Zhi-feng,YIGui-xi,et al.1994.Location of Tibetan earthquakes:A nonlinear approach by a simplex optimized technique[J].Acta Seismologica Sinica,16(2):212—219(in Chinese).

    趙珠,范軍,鄭斯華,等.1997.龍門山斷裂帶地殼速度結(jié)構(gòu)和震源位置的精確修定[J].地震學(xué)報,19(6):615—622.

    ZHAO Zhu,F(xiàn)AN Jun,ZHENG Si-hua,et al.1997.Crustal structure and accurate hypocenter determination along the Longmenshan Fault zone[J].Acta Seismologica Sinica,19(6):615—622(in Chinese).

    趙仲和.1983.多重模型地震定位程序及其在北京臺網(wǎng)的應(yīng)用[J].地震學(xué)報,5(2):242—254.

    ZHONG Zhong-he.1983.An earthquake location program withmultiple velocitymodel and its application in the Beijing seismic network[J].Acta Seismologica Sinica,5(2):242—254(in Chinese).

    周仕勇,許忠淮,韓京,等.1999.主地震定位方法分析以及1997年新疆伽師震群高精度定位[J].地震學(xué)報,21(3):258—265.

    ZHOU Shi-yong,XU Zhong-huai,HAN Jing,et al.1999.Analysis on the master eventmethod and precise location of 1997 Jiashi strong earthquake swarm in western China[J].Acta Seismologica Sinica,21(3):258—265(in Chinese).

    Geiger L.1912.Probability method for the determination of earthquake epicenters from arrival time only [J].Bull St Louis Univ,8:60—71.

    Jackson D D,Matsu'ura M.1985.A Bayesian approach to nonlinear inversion[J].JGeophys Res,90(B1):581—591.

    Lee W H K,Lahr JC.1975.A computer program for determining hypocenter,magnitude,and firstmotion pattern of local earthquakes[J].U SGeol Surv Open-file Rep,75—311.

    Mooney W D,Laske G,Masters G.1998.A global crustalmodel at5°x 5°[J].JGeophys Re,103:727—747.

    Shen Z K,Jackson D,Ge B.1996.Crustal deformation across and beyond the Los Angeles Basin from geodeticmeasurements[J].JGeophys Res,101(B12):27957—27980.

    SpenceW.1980.Relative epicenter determination using P-wave arrival-time differences[J].Bull Seism Soc Amer,70(1):171—183.

    Takagi A,Hasegawa A,Umino N.1977.Seismic activity in the Northeastern Janpan Arc[J].JPhys Earth,25(Suppl):95—104.

    Waldhauser F,Ellsworth W L.2000.A double-difference earthquake location algorithm:Method and application to the Northern Hayward Fault,California[J].Bull Seism Soc Amer,90(6):1353—1368.

    Wu Y M,Chang C H,Zhao L,et al.2008.A comprehensive relocation of earthquakes in Taiwan from 1991 to 2005[J].Bull Seism Soc Amer,98(3):1471—1481.

    EARTHQUAKE LOCATION METHODW ITH ARRIVAL TIME UNCERTAINTY CONSIDERED AND ITS APPLICATION TO LOCATION OF EARTHQUAKES FROM 2001 TO 2008 IN SICHUAN AREA

    WAN Yong-ge1)SHENG Shu-zhong1)CHENGWan-zheng2)ZHANG Zhi-wei2)WU Ye1)ZHAO Xiao-yan1)BU Yu-fei3)XUE Zhi-fang4)LIU Jin-li5)

    1)Institute of Disaster Prevention,Yanjiao,Hebei 065201,China

    2)Sichuan Earthquake Administration,Chengdu 610041,China

    3)Gansu Earthquake Administration,Lanzhou 730000,China

    4)Qinhuangdao Center Station,Hebei Earthquake Administration,Qinhuangdao 066100,China

    5)Quzhou Experimental High School,Hebei Province,Quzhou 057250,China

    The precisely located earthquake catalogue is important to seismicity,seismic tomography and crustal stress inversion studies.It also has greatapplication value in rapid report of an earthquake that just occurred.By considering the arrival time uncertainty,and the constraints on station elevation and seismic depth,we propose a relatively accuratemethod to estimate hypocentral location and its uncertainty based on inversion theory.Ourmethod can combine the arrival times of Pg wave,Sg wave,Pn wave and Sn wave in hypocenter location,so it increases the location accuracy by involvingmore data;and it can be also used in local and regional earthquake location simultaneously.In order to test our location method,we located earthquakes by using the simulated data with different uncertainty of Pg,Sg,Pn,Sn arrivals.The result shows that the location determined by using ourmethod ismore accurate than that by using othermethod.We apply it to earthquakes occurring in the period from 2001 to 2008 in Sichuan area,and obtained amore clustered hypocentral distribution convergent to the fault zones.The result provides a solid foundation for studies of seismicity,geometry of the active faults and seismic tomography in Sichuan region.It is also helpful to study the seismicity precursors before theWenchuan earthquake.

    earthquake location,uncertainty,station elevation,Sichuan area

    P315.5

    A

    0253-4967(2012)01-0001-10

    10.3969/j.issn.0253-4967.2012.01.001

    2011-05-03收稿,2011-07-30改回。

    國家自然科學(xué)基金(41074072、40874022)、中國地震局教師科研基金(20100101)、中央高?;究蒲袠I(yè)務(wù)費專項(ZY20110101)和地震行業(yè)專項(200808053,201208009)共同資助。

    萬永革,男,1967年生,2001年在中國地震局地球物理研究所獲得博士學(xué)位,2002—2006年在中國地震局地質(zhì)研究所進(jìn)行博士后合作研究,研究員,主要從事地震學(xué)和地球動力學(xué)方面的研究,電話:010-61597607;E-mail:wanyg217217@vip.sina.com;wanyongge8@gmail.com。

    猜你喜歡
    臺網(wǎng)臺站震源
    中國科學(xué)院野外臺站檔案工作回顧
    氣象基層臺站建設(shè)
    西藏科技(2021年12期)2022-01-17 08:46:38
    地球物理臺網(wǎng)儀器維修信息管理的研究與實現(xiàn)
    地震研究(2021年1期)2021-04-13 01:05:30
    推進(jìn)報臺網(wǎng)深度融合 做強區(qū)級融媒體中心
    傳媒評論(2018年10期)2019-01-17 01:33:16
    震源的高返利起步
    西藏地震應(yīng)急流動臺網(wǎng)淺析
    西藏科技(2016年9期)2016-09-26 12:21:39
    MRG9000媒資衛(wèi)士在市級電視臺全臺網(wǎng)的應(yīng)用
    可控震源地震在張掖盆地南緣逆沖斷裂構(gòu)造勘探中的應(yīng)用
    基層臺站綜合觀測業(yè)務(wù)管理之我見
    西藏科技(2015年6期)2015-09-26 12:12:13
    同步可控震源地震采集技術(shù)新進(jìn)展
    国产成人a∨麻豆精品| 少妇人妻久久综合中文| 国产男女内射视频| 亚洲婷婷狠狠爱综合网| 简卡轻食公司| 亚洲av综合色区一区| 精品亚洲成国产av| 黄色一级大片看看| 国内少妇人妻偷人精品xxx网站| 国内揄拍国产精品人妻在线| 国产精品精品国产色婷婷| 男人和女人高潮做爰伦理| 蜜桃久久精品国产亚洲av| 麻豆成人av视频| 男女免费视频国产| 九九久久精品国产亚洲av麻豆| 免费大片18禁| a级毛色黄片| 老女人水多毛片| 久久毛片免费看一区二区三区| 国产免费福利视频在线观看| 亚洲精品乱久久久久久| 美女福利国产在线 | 人妻 亚洲 视频| 免费观看av网站的网址| 黄片无遮挡物在线观看| 亚洲熟女精品中文字幕| 久久精品久久久久久噜噜老黄| 国产伦精品一区二区三区视频9| 欧美日韩视频精品一区| 视频区图区小说| 一级毛片 在线播放| 亚洲怡红院男人天堂| 成人国产av品久久久| 亚洲天堂av无毛| 久久久久精品久久久久真实原创| 九九久久精品国产亚洲av麻豆| 国产成人91sexporn| 日本wwww免费看| 人人妻人人添人人爽欧美一区卜 | 国产有黄有色有爽视频| 99久久精品热视频| 天天躁日日操中文字幕| 亚洲美女搞黄在线观看| 人人妻人人看人人澡| 日韩欧美 国产精品| 一级毛片 在线播放| 在线观看免费高清a一片| 美女高潮的动态| 日韩欧美一区视频在线观看 | 日韩制服骚丝袜av| 天美传媒精品一区二区| 狂野欧美激情性bbbbbb| 欧美极品一区二区三区四区| 国产色婷婷99| 亚洲人成网站在线观看播放| av专区在线播放| 国产精品.久久久| 久久久久久久久久久免费av| 熟女电影av网| 国产乱人视频| 精品国产三级普通话版| 少妇的逼好多水| 成人18禁高潮啪啪吃奶动态图 | 国产高清三级在线| 高清视频免费观看一区二区| 你懂的网址亚洲精品在线观看| freevideosex欧美| 久久久成人免费电影| 舔av片在线| 国模一区二区三区四区视频| 亚洲成人手机| 99久久精品国产国产毛片| 免费观看性生交大片5| 一区二区三区免费毛片| 亚洲精品第二区| 国产中年淑女户外野战色| 视频中文字幕在线观看| 中国国产av一级| 人妻夜夜爽99麻豆av| 国产免费福利视频在线观看| 大话2 男鬼变身卡| 永久网站在线| 99视频精品全部免费 在线| 美女脱内裤让男人舔精品视频| videos熟女内射| 麻豆精品久久久久久蜜桃| 高清欧美精品videossex| 日本av手机在线免费观看| 在线观看一区二区三区| av播播在线观看一区| 日韩制服骚丝袜av| 日本午夜av视频| 日本黄大片高清| 亚洲第一区二区三区不卡| 久久av网站| 91狼人影院| 国产欧美日韩精品一区二区| 亚洲精品一区蜜桃| 热re99久久精品国产66热6| 99久久人妻综合| 成年女人在线观看亚洲视频| 欧美xxⅹ黑人| 成人毛片a级毛片在线播放| 亚洲欧美精品专区久久| 爱豆传媒免费全集在线观看| av在线播放精品| 下体分泌物呈黄色| 新久久久久国产一级毛片| 亚洲最大成人中文| 免费观看a级毛片全部| 久久久久视频综合| 嫩草影院入口| 欧美成人精品欧美一级黄| 夜夜爽夜夜爽视频| av.在线天堂| 亚洲精品日韩av片在线观看| 国产国拍精品亚洲av在线观看| 91精品伊人久久大香线蕉| 青春草国产在线视频| 国产在线男女| 精品国产露脸久久av麻豆| 哪个播放器可以免费观看大片| 丰满人妻一区二区三区视频av| 久久久久精品久久久久真实原创| 国产高清有码在线观看视频| 亚洲av综合色区一区| 欧美变态另类bdsm刘玥| 波野结衣二区三区在线| 免费高清在线观看视频在线观看| 一级片'在线观看视频| 最近中文字幕高清免费大全6| 国产中年淑女户外野战色| 久久精品久久久久久久性| 尤物成人国产欧美一区二区三区| 一个人免费看片子| 国产高清不卡午夜福利| 精品人妻视频免费看| 国产男人的电影天堂91| 高清不卡的av网站| 妹子高潮喷水视频| 色5月婷婷丁香| 免费观看av网站的网址| 一区二区av电影网| 欧美成人a在线观看| 青春草国产在线视频| 亚洲精品456在线播放app| 久久影院123| 久久97久久精品| videossex国产| 激情 狠狠 欧美| 色哟哟·www| 26uuu在线亚洲综合色| 视频区图区小说| 日本与韩国留学比较| 亚洲美女黄色视频免费看| 搡女人真爽免费视频火全软件| 精品人妻视频免费看| h日本视频在线播放| 网址你懂的国产日韩在线| 一级爰片在线观看| 国产伦精品一区二区三区四那| 国产成人午夜福利电影在线观看| 小蜜桃在线观看免费完整版高清| 大陆偷拍与自拍| 成年美女黄网站色视频大全免费 | 亚洲精品一二三| 欧美日本视频| 日本黄色日本黄色录像| 大又大粗又爽又黄少妇毛片口| 性色av一级| 嘟嘟电影网在线观看| av黄色大香蕉| 亚洲丝袜综合中文字幕| 国产伦理片在线播放av一区| 三级国产精品欧美在线观看| kizo精华| 久久久久国产精品人妻一区二区| 久久久色成人| 热99国产精品久久久久久7| 成年免费大片在线观看| 欧美 日韩 精品 国产| 99热这里只有是精品50| 伦精品一区二区三区| 亚洲三级黄色毛片| 乱系列少妇在线播放| 男人和女人高潮做爰伦理| 99热这里只有是精品在线观看| 久久久久久久大尺度免费视频| 极品少妇高潮喷水抽搐| 国产亚洲一区二区精品| 精品亚洲成国产av| 欧美老熟妇乱子伦牲交| av福利片在线观看| 欧美老熟妇乱子伦牲交| tube8黄色片| 亚洲精品456在线播放app| 97精品久久久久久久久久精品| 国产视频内射| 精品亚洲成a人片在线观看 | 久久99热这里只频精品6学生| 激情 狠狠 欧美| 国产黄频视频在线观看| 日韩中文字幕视频在线看片 | 日本一二三区视频观看| 精品久久久久久久久av| 夜夜看夜夜爽夜夜摸| 波野结衣二区三区在线| 国产乱人视频| 国产精品秋霞免费鲁丝片| 国产真实伦视频高清在线观看| 日韩av不卡免费在线播放| 成人亚洲欧美一区二区av| 精品熟女少妇av免费看| 亚洲天堂av无毛| 黑人猛操日本美女一级片| 简卡轻食公司| a级毛色黄片| 成人二区视频| 一级av片app| 丝瓜视频免费看黄片| 欧美人与善性xxx| 九九爱精品视频在线观看| 国产爽快片一区二区三区| 男女无遮挡免费网站观看| 国产精品国产av在线观看| 午夜福利在线观看免费完整高清在| 少妇丰满av| av一本久久久久| 男女边吃奶边做爰视频| 国产欧美另类精品又又久久亚洲欧美| 国精品久久久久久国模美| 秋霞在线观看毛片| 欧美丝袜亚洲另类| 日本与韩国留学比较| 永久免费av网站大全| 国内少妇人妻偷人精品xxx网站| 少妇被粗大猛烈的视频| 99热网站在线观看| 国产欧美另类精品又又久久亚洲欧美| av一本久久久久| 亚州av有码| 国产亚洲午夜精品一区二区久久| av国产精品久久久久影院| 男人爽女人下面视频在线观看| 日本爱情动作片www.在线观看| 亚洲欧洲日产国产| 亚洲国产欧美在线一区| 黄片wwwwww| 我的老师免费观看完整版| 久久久成人免费电影| 日本色播在线视频| 日本av免费视频播放| 亚洲电影在线观看av| 一区在线观看完整版| 国产一区二区在线观看日韩| 男女边摸边吃奶| 亚洲成人中文字幕在线播放| 亚洲一级一片aⅴ在线观看| 国产一区二区三区av在线| 亚洲精品色激情综合| 舔av片在线| 国产免费一区二区三区四区乱码| 51国产日韩欧美| 国产精品久久久久久精品古装| 欧美变态另类bdsm刘玥| 赤兔流量卡办理| 干丝袜人妻中文字幕| 中文在线观看免费www的网站| 99re6热这里在线精品视频| 欧美日韩一区二区视频在线观看视频在线| 免费人妻精品一区二区三区视频| 久久99热这里只有精品18| 亚洲国产欧美人成| 欧美成人a在线观看| 一级黄片播放器| 99久久精品一区二区三区| 身体一侧抽搐| 日韩伦理黄色片| 在线观看三级黄色| 啦啦啦视频在线资源免费观看| 亚洲精品成人av观看孕妇| 免费黄频网站在线观看国产| 亚洲av成人精品一区久久| 国产精品久久久久成人av| 国产黄色视频一区二区在线观看| 免费看av在线观看网站| 亚洲国产最新在线播放| 2021少妇久久久久久久久久久| 特大巨黑吊av在线直播| 又黄又爽又刺激的免费视频.| 亚洲美女搞黄在线观看| 乱系列少妇在线播放| 少妇熟女欧美另类| 国产淫语在线视频| 精品国产一区二区三区久久久樱花 | 男女啪啪激烈高潮av片| av免费在线看不卡| 97在线视频观看| 亚洲,一卡二卡三卡| 亚洲国产精品成人久久小说| 老熟女久久久| 亚洲久久久国产精品| 亚洲精品亚洲一区二区| 一本—道久久a久久精品蜜桃钙片| 欧美性感艳星| 简卡轻食公司| 五月天丁香电影| 在线观看人妻少妇| 成人高潮视频无遮挡免费网站| 男人舔奶头视频| 国产精品人妻久久久影院| 亚洲国产欧美在线一区| 久久午夜福利片| 国内精品宾馆在线| 亚洲国产av新网站| 爱豆传媒免费全集在线观看| 精品久久国产蜜桃| 一二三四中文在线观看免费高清| 亚洲精品456在线播放app| 国产精品蜜桃在线观看| 日韩电影二区| 精品少妇黑人巨大在线播放| 人妻制服诱惑在线中文字幕| 91久久精品电影网| 日日撸夜夜添| 色吧在线观看| 99九九线精品视频在线观看视频| 99久国产av精品国产电影| 91精品国产九色| 高清在线视频一区二区三区| 小蜜桃在线观看免费完整版高清| 人体艺术视频欧美日本| 国产伦在线观看视频一区| 日本-黄色视频高清免费观看| 亚洲熟女精品中文字幕| 欧美精品国产亚洲| 亚洲av电影在线观看一区二区三区| 五月天丁香电影| 一级爰片在线观看| 成人亚洲欧美一区二区av| 久久女婷五月综合色啪小说| 欧美精品一区二区免费开放| 精品午夜福利在线看| 最近中文字幕高清免费大全6| 国产乱来视频区| 久久久精品免费免费高清| 日韩中字成人| 赤兔流量卡办理| 国产淫语在线视频| 欧美日韩综合久久久久久| 一区在线观看完整版| 我的女老师完整版在线观看| 天堂8中文在线网| 干丝袜人妻中文字幕| 免费少妇av软件| 91精品一卡2卡3卡4卡| 亚洲精品aⅴ在线观看| 91久久精品国产一区二区三区| 日韩强制内射视频| 欧美成人一区二区免费高清观看| 久久精品夜色国产| 久久99热6这里只有精品| 特大巨黑吊av在线直播| 五月开心婷婷网| 亚洲精品日韩av片在线观看| 亚洲av电影在线观看一区二区三区| 亚洲欧美一区二区三区国产| 国产精品麻豆人妻色哟哟久久| 高清午夜精品一区二区三区| 97超视频在线观看视频| 午夜激情久久久久久久| 欧美日韩国产mv在线观看视频 | 久久久久精品久久久久真实原创| 亚洲欧美日韩卡通动漫| 国产久久久一区二区三区| 人妻制服诱惑在线中文字幕| 国产男女超爽视频在线观看| 日本午夜av视频| 国产免费一区二区三区四区乱码| 欧美成人精品欧美一级黄| 日本欧美视频一区| 国产又色又爽无遮挡免| 高清黄色对白视频在线免费看 | 大片电影免费在线观看免费| 在线免费十八禁| 亚洲精品中文字幕在线视频 | 我要看黄色一级片免费的| 亚洲高清免费不卡视频| 黄色日韩在线| 人妻制服诱惑在线中文字幕| 亚洲最大成人中文| 熟女av电影| 精品人妻视频免费看| av一本久久久久| 久久精品国产亚洲网站| 联通29元200g的流量卡| h视频一区二区三区| 日韩欧美 国产精品| 亚洲欧美精品自产自拍| 亚洲国产欧美在线一区| av国产久精品久网站免费入址| 99热这里只有精品一区| 亚洲,欧美,日韩| 国产久久久一区二区三区| 亚洲性久久影院| 精品人妻偷拍中文字幕| 青春草亚洲视频在线观看| 97热精品久久久久久| 免费看日本二区| 久久精品国产亚洲av涩爱| 午夜福利网站1000一区二区三区| 亚洲av成人精品一二三区| 丰满乱子伦码专区| 一本久久精品| 国模一区二区三区四区视频| 国产精品麻豆人妻色哟哟久久| 身体一侧抽搐| 伦理电影免费视频| 在线免费观看不下载黄p国产| 一级av片app| 乱系列少妇在线播放| 亚洲av中文字字幕乱码综合| 久久女婷五月综合色啪小说| 午夜激情久久久久久久| 国产一区亚洲一区在线观看| 男男h啪啪无遮挡| 欧美精品亚洲一区二区| 免费看av在线观看网站| 夜夜骑夜夜射夜夜干| av黄色大香蕉| 寂寞人妻少妇视频99o| 欧美bdsm另类| 观看美女的网站| 亚洲av不卡在线观看| 国产在线视频一区二区| 成人国产麻豆网| 精品一区在线观看国产| 麻豆精品久久久久久蜜桃| 精品人妻一区二区三区麻豆| 高清视频免费观看一区二区| 成人黄色视频免费在线看| 天天躁日日操中文字幕| 久久久久性生活片| 亚洲国产日韩一区二区| 视频区图区小说| 永久网站在线| 亚洲第一区二区三区不卡| 最黄视频免费看| 国产男人的电影天堂91| 大片电影免费在线观看免费| 热99国产精品久久久久久7| 国产欧美日韩精品一区二区| 国产成人精品福利久久| 三级经典国产精品| 久久久久久人妻| 久久久久久久久久成人| 亚洲激情五月婷婷啪啪| 国产亚洲欧美精品永久| 热99国产精品久久久久久7| 岛国毛片在线播放| 国产伦精品一区二区三区视频9| 少妇猛男粗大的猛烈进出视频| 久久久久久久亚洲中文字幕| 97超碰精品成人国产| 国产伦在线观看视频一区| 在线观看免费日韩欧美大片 | 国产女主播在线喷水免费视频网站| 高清av免费在线| 欧美zozozo另类| 国产美女午夜福利| 欧美+日韩+精品| 亚洲国产欧美在线一区| 国产黄片视频在线免费观看| 日本爱情动作片www.在线观看| 欧美成人午夜免费资源| 国产永久视频网站| av国产久精品久网站免费入址| av视频免费观看在线观看| 国产真实伦视频高清在线观看| 卡戴珊不雅视频在线播放| 国产精品爽爽va在线观看网站| 麻豆成人午夜福利视频| h日本视频在线播放| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲国产精品专区欧美| 99久久精品热视频| 18禁裸乳无遮挡免费网站照片| 男女无遮挡免费网站观看| 亚洲电影在线观看av| 多毛熟女@视频| 91精品国产九色| 美女中出高潮动态图| 亚洲精品日韩av片在线观看| 啦啦啦视频在线资源免费观看| 乱码一卡2卡4卡精品| 日韩国内少妇激情av| 黑人猛操日本美女一级片| 国产 一区 欧美 日韩| freevideosex欧美| 日日摸夜夜添夜夜爱| 欧美 日韩 精品 国产| 久久久亚洲精品成人影院| 国产精品精品国产色婷婷| 欧美日韩亚洲高清精品| 亚洲国产高清在线一区二区三| 高清av免费在线| 亚洲欧美日韩另类电影网站 | 精品国产露脸久久av麻豆| 国产精品三级大全| 成人一区二区视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 美女视频免费永久观看网站| 只有这里有精品99| 亚洲精品久久午夜乱码| 不卡视频在线观看欧美| 亚洲av.av天堂| 黄色视频在线播放观看不卡| 网址你懂的国产日韩在线| 国国产精品蜜臀av免费| 一级黄片播放器| 精品国产露脸久久av麻豆| 久久久精品94久久精品| 免费在线观看成人毛片| 国产淫语在线视频| 又大又黄又爽视频免费| 狂野欧美白嫩少妇大欣赏| 岛国毛片在线播放| 国产精品国产三级国产av玫瑰| 欧美bdsm另类| 亚洲美女搞黄在线观看| freevideosex欧美| 亚洲欧美日韩另类电影网站 | 国产精品不卡视频一区二区| 最近最新中文字幕免费大全7| 欧美极品一区二区三区四区| av国产免费在线观看| 亚洲精品乱码久久久v下载方式| 波野结衣二区三区在线| 高清不卡的av网站| 国产精品一区二区在线观看99| 日本vs欧美在线观看视频 | 国产黄色免费在线视频| 久久久久国产精品人妻一区二区| 亚洲真实伦在线观看| 久久精品久久久久久久性| 大话2 男鬼变身卡| 狂野欧美白嫩少妇大欣赏| 日韩av在线免费看完整版不卡| a级毛色黄片| 国产成人精品久久久久久| 六月丁香七月| 久久av网站| 人人妻人人添人人爽欧美一区卜 | 日韩欧美 国产精品| 乱系列少妇在线播放| 国产有黄有色有爽视频| 一本一本综合久久| 国产爱豆传媒在线观看| 国内精品宾馆在线| 免费av不卡在线播放| 一区二区三区免费毛片| 国产免费一级a男人的天堂| 狂野欧美激情性bbbbbb| av在线观看视频网站免费| 最近2019中文字幕mv第一页| 22中文网久久字幕| 99久久精品国产国产毛片| 久久久久久久久久人人人人人人| 丰满少妇做爰视频| 亚洲欧美精品自产自拍| 舔av片在线| 欧美人与善性xxx| 美女xxoo啪啪120秒动态图| 在线观看免费视频网站a站| 色哟哟·www| 日本欧美视频一区| 肉色欧美久久久久久久蜜桃| 国产精品熟女久久久久浪| 91久久精品电影网| 我要看黄色一级片免费的| 99国产精品免费福利视频| 少妇猛男粗大的猛烈进出视频| 全区人妻精品视频| 狂野欧美白嫩少妇大欣赏| 国产精品一区二区性色av| 王馨瑶露胸无遮挡在线观看| 精品国产露脸久久av麻豆| 在线免费十八禁| 精品人妻偷拍中文字幕| 超碰97精品在线观看| 亚洲天堂av无毛| 亚洲av在线观看美女高潮| 中文字幕av成人在线电影| 熟女人妻精品中文字幕| 乱码一卡2卡4卡精品| 91久久精品电影网| 色视频www国产| 亚洲欧美成人综合另类久久久| tube8黄色片| 一级a做视频免费观看| 精品少妇黑人巨大在线播放| 春色校园在线视频观看| 免费看光身美女| 欧美xxxx黑人xx丫x性爽| 亚洲欧美精品自产自拍| 久久久久精品性色| 国产精品福利在线免费观看| 欧美极品一区二区三区四区| 美女国产视频在线观看| 最近最新中文字幕大全电影3| 国国产精品蜜臀av免费| 亚洲丝袜综合中文字幕| 在线观看免费高清a一片| 色综合色国产|