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

    青海北部中小地震視應(yīng)力特征研究

    2022-01-26 05:09:04劉文邦李啟雷張朋濤屠泓為青海省地震局青海西寧810001
    地震工程學(xué)報(bào) 2022年1期
    關(guān)鍵詞:拐角震級(jí)應(yīng)力場(chǎng)

    劉文邦, 李啟雷, 張朋濤, 屠泓為(青海省地震局, 青海 西寧 810001)

    0 引言

    青海北部位于青藏高原東北緣,在NEE向的阿爾金斷裂帶、NWW向祁連山斷裂帶和NWW向東昆侖斷裂帶等3條巨型左旋走滑斷裂所圍限的活動(dòng)塊體內(nèi)部。該地區(qū)構(gòu)造變形強(qiáng)烈,發(fā)育了大量褶皺、逆沖和走滑斷裂,是高原隆升和變形最敏感的地區(qū)之一[1]。有研究表明,自1997年西藏瑪尼7.5級(jí)地震以來(lái),中國(guó)大陸圍繞青藏高原巴顏喀拉地塊邊界發(fā)生了昆侖山口西8.1級(jí)地震、汶川8.0級(jí)地震及玉樹(shù)7.1級(jí)地震等8次7級(jí)以上地震,然而,近年來(lái)一些6級(jí)以上地震未發(fā)生在巴顏喀拉地塊及鄰區(qū),而是發(fā)生在青藏地塊周緣的邊界構(gòu)造上[2]。作為巴顏喀拉塊體的相鄰構(gòu)造區(qū),青海北部發(fā)生了2014年10月2日青海烏蘭5.1級(jí)、2015年11月23日青海祁連5.2級(jí)、2016年1月21日青海門(mén)源6.4級(jí)、2019年9月16日甘肅甘州5.0級(jí)、2019年10月28日甘肅夏河5.7級(jí)等地震,地震活動(dòng)顯著增強(qiáng)。青海北部地區(qū)的地震危險(xiǎn)性有待進(jìn)一步研究。

    地震孕育過(guò)程中可能伴隨著地殼應(yīng)力變化,研究其變化與地震的關(guān)系是揭示地震前兆、探索具有明確物理意義的地震預(yù)報(bào)方法的途徑之一。近年來(lái),視應(yīng)力是一個(gè)與震源動(dòng)力學(xué)過(guò)程相關(guān)的動(dòng)力學(xué)參數(shù),可以反映區(qū)域應(yīng)力場(chǎng)的強(qiáng)弱,在地殼應(yīng)力場(chǎng)監(jiān)視和地震預(yù)報(bào)中得到廣泛應(yīng)用[3-8]。吳忠良等[9]分析了視應(yīng)力在地震危險(xiǎn)性及地震預(yù)測(cè)中的應(yīng)用。易桂喜等[10]發(fā)現(xiàn)蘆山7.0級(jí)地震前龍門(mén)山斷裂帶南段發(fā)生高視應(yīng)力地震。李艷娥等[11]研究了日本MW9.1地震前后破裂區(qū)地震視應(yīng)力時(shí)空變化,發(fā)現(xiàn)震前數(shù)年裂區(qū)存在大面積視應(yīng)力高值。劉紅桂等[12]研究發(fā)現(xiàn)在云南地區(qū),高視應(yīng)力地震的發(fā)生可作為預(yù)測(cè)該地區(qū)未來(lái)發(fā)生中強(qiáng)地震的一個(gè)參考指標(biāo)。楊志高等[13]認(rèn)為水庫(kù)地區(qū)蓄水前后視應(yīng)力標(biāo)度率存在顯著變化。鄭建常和張博等[14-15]通過(guò)視應(yīng)力等震源動(dòng)力學(xué)參數(shù)對(duì)震群類(lèi)型和發(fā)震過(guò)程進(jìn)行了分析。

    本文利用青海北部地區(qū)的數(shù)字化地震波資料,反演計(jì)算了2010年以來(lái)青海北部地區(qū)中小地震視應(yīng)力。通過(guò)分析研究及預(yù)測(cè)效能檢驗(yàn),得到一些有意義的結(jié)果,對(duì)青海北部地震視應(yīng)力的預(yù)測(cè)指示作用有一定的認(rèn)識(shí)。

    1 計(jì)算方法

    根據(jù)視應(yīng)力(σapp)定義[16]:

    (1)

    式中:ES為地震輻射能量;M0為地震矩;μ為剪切能量(對(duì)于地殼介質(zhì) ,μ取3×104MPa)。由上式可知,計(jì)算地震視應(yīng)力的關(guān)鍵是先計(jì)算地震波輻射能量和地震矩。 通過(guò)對(duì)數(shù)字波形資料進(jìn)行震源譜反演分析,可計(jì)算得到地震輻射能量和地震矩。本文具體計(jì)算步驟參考文獻(xiàn)[17]。

    對(duì)中小地震,震源譜符合Brune圓盤(pán)模型[18],震源譜可表示為:

    (2)

    式中:Ω0為震源譜零頻極限值;fC為拐角頻率。當(dāng)Ω0和fC給定時(shí),式(2) 即確定。因此,確定Ω0和fC就成為震源譜計(jì)算的主要內(nèi)容。地震矩M0可以根據(jù)式(3)求得:

    (3)

    式中:ρ為地殼介質(zhì)密度(可取2.71 kg/m3);v為波速(P波取6.1 km/s,S波取3.5 km/s);d為震源距;Ω0為震源譜零頻極限值;R為輻射因子。

    本文所用的波形資料都為速度記錄,三分向速度記錄變換到頻率域,利用下式可得到速度譜V(f),即:

    (4)

    地震能量ES可由對(duì)速度譜的平方積分求得,即:

    (5)

    由于地震計(jì)頻帶寬度的有限性,在計(jì)算過(guò)程中考慮低頻和高頻補(bǔ)償,即:

    (6)

    式中:Ω0為零頻極限;β為S波速度;f1、f3分別為擬合震源譜過(guò)程中選擇的最低頻率和最高頻率,中小地震振幅譜符合Brune 模型,低頻段為水平段,高頻段為衰減段,計(jì)算過(guò)程中f1和f2根據(jù)實(shí)際情況在振幅譜的平坦部分選擇,而f3選在振幅譜的衰減段。Ω(f3)對(duì)應(yīng)頻率f3的振幅值,可由將Ω0、fc和f3代入式(2)中計(jì)算理論值代替。由于地震能量由S波攜帶,本文主要計(jì)算S波的地震能量。

    對(duì)某一地震,由于地震輻射花樣及破裂方向的影響,每個(gè)臺(tái)站計(jì)算結(jié)果有所不同,圖1為2020年4月7日青海祁連ML3.7地震QIL、TIJ兩個(gè)臺(tái)站的位移譜。

    圖1 2020年4月7日祁連ML3.7地震QIL、TIJ兩個(gè)臺(tái)位移譜Fig.1 Displacement spectra of QIL and TIJ stations during the Qilian ML3.7 earthquake on April 7,2020

    為了消除個(gè)別臺(tái)站的異常高值對(duì)平均值的影響,在由各個(gè)臺(tái)站的值求平均時(shí)采用Archuleta等的方法[19],即式中,xi為各臺(tái)站的地震矩或地震能量;N為臺(tái)站數(shù);Δx為誤差因子,其意義為當(dāng)x以對(duì)數(shù)坐標(biāo)作圖時(shí)的標(biāo)準(zhǔn)差。

    (7)

    (8)

    2 使用資料

    本文使用了青海測(cè)震臺(tái)網(wǎng)記錄的青海北部2010年01月至2020年5月的ML≥3.0地震數(shù)字化波形資料,考慮到視應(yīng)力計(jì)算要求,選用了震中距在180 km以?xún)?nèi)的臺(tái)站,性噪比較高記錄清晰的臺(tái)站波形資料。對(duì)這些資料進(jìn)行去傾斜及儀器響應(yīng)的校正處理,選用S波段波形進(jìn)行震源譜計(jì)算。圖2為震中、臺(tái)站及次級(jí)塊體分布。

    3 計(jì)算結(jié)果

    利用青海區(qū)域臺(tái)網(wǎng)地震波形資料,計(jì)算了青海北部2010年1月至2020年5月共637個(gè)ML3.0以上地震的震源譜參數(shù),其中ML3.0~3.9地震537個(gè),ML4.0~4.9地震85個(gè),ML≥5.0地震15個(gè)。計(jì)算得到了它們的拐角頻率、地震矩、地震輻射能量及視應(yīng)力等震源動(dòng)力學(xué)參數(shù)。

    圖2 研究區(qū)域2010年1月至2020年5月 ML≥3.0地震震中及臺(tái)站分布Fig.2 Epicenter and station distribution of ML ≥ 3.0 earthquakes in the study area from January 2010 to May 2020

    3.1 震源參數(shù)的標(biāo)度關(guān)系

    通過(guò)分析637個(gè)ML3.0以上地震的震源參數(shù),得到震源參數(shù)與震級(jí)的標(biāo)度關(guān)系,如圖3所示。

    如圖3,震源參數(shù)與震級(jí)的擬合線(xiàn)性趨勢(shì)明顯。其中地震矩、地震輻射能量與震級(jí)均呈正相關(guān)性,拐角頻率與震級(jí)呈負(fù)相關(guān)性,地震矩、地震輻射能量與震級(jí)的相關(guān)性較拐角頻率及視應(yīng)力與震級(jí)相關(guān)性更強(qiáng)。ML3.0~6.4地震的拐角頻率分布在0.23~6.75 Hz,拐角頻率與震級(jí)呈負(fù)相關(guān)。拐角頻率反映震源尺度大小的物理量,研究發(fā)現(xiàn)拐角頻率和震級(jí)存在相關(guān)關(guān)系,地震越大其頻譜的低頻成分越豐富,拐角頻率越低[14]。

    圖3 震源參數(shù)與震級(jí)之間的標(biāo)度關(guān)系Fig.3 Scaling relationship between focal parameters and magnitude

    視應(yīng)力分布在0.003~2.39 MPa,平均值為0.15 MPa,視應(yīng)力與震級(jí)存在較明顯正相關(guān)性,一般來(lái)說(shuō),震級(jí)越大視應(yīng)力越高,但也存在一些低震級(jí)高視應(yīng)力地震。視應(yīng)力與震級(jí)間的變化關(guān)系較復(fù)雜,研究認(rèn)為與斷層的滑動(dòng)類(lèi)型、破裂過(guò)程及介質(zhì)強(qiáng)度等有關(guān),在介質(zhì)較破碎的斷層弱化帶無(wú)法積累較大的應(yīng)力,在堅(jiān)硬的巖石層發(fā)生的地震導(dǎo)致應(yīng)力的集中釋放可以導(dǎo)致區(qū)域構(gòu)造應(yīng)力場(chǎng)的改變而引發(fā)地震,高應(yīng)力的中小地震可以作為預(yù)測(cè)中強(qiáng)地震的重要參考指標(biāo)之一[10-12],下面在研究視應(yīng)力時(shí)空分布時(shí),主要討論高視應(yīng)力地震與中強(qiáng)地震的關(guān)系。

    3.2 青海北部視應(yīng)力隨時(shí)間變化特征

    考慮到視應(yīng)力與震級(jí)之間存在正相關(guān)性,將相差較大的兩個(gè)震級(jí)檔地震放在一起比較,很難區(qū)分震級(jí)大導(dǎo)致的高視應(yīng)力還是高視應(yīng)力地震。鑒于此,在分析視應(yīng)力時(shí)空變化特征時(shí),選取震級(jí)相近的震級(jí)檔的進(jìn)行視應(yīng)力對(duì)比,本文通過(guò)ML3.0~3.5、ML3.6~4.0兩個(gè)震級(jí)檔的視應(yīng)力時(shí)空變化特征,研究青海北部中小地震視應(yīng)力與中強(qiáng)地震的關(guān)系。

    由圖4可知,分震級(jí)檔視應(yīng)力隨時(shí)間變化,ML3.0~3.5震級(jí)檔視應(yīng)力在兩個(gè)時(shí)間段視應(yīng)力出現(xiàn)增強(qiáng)現(xiàn)象:2013—2015年間高視應(yīng)力地震顯著增多,尤其是多次地震視應(yīng)力超過(guò)0.5 MPa,5級(jí)以上地震同步有所增強(qiáng),在柴達(dá)木盆地內(nèi)部發(fā)生了2013年2月12日冷湖5.1級(jí)、2013年6月5日大柴旦5.0級(jí)、2014年10月2日烏蘭5.1級(jí)等地震,在祁連山地震帶發(fā)生了2013年9月20日門(mén)源5.1級(jí)、2015年11月23日祁連5.2級(jí)、2016年1月21日門(mén)源6.4級(jí)等地震;2018—2019年間視應(yīng)力高值異常顯著,發(fā)生了2019年3月28日茫崖5.0級(jí)、2019年9月6日甘州5.0級(jí)、2019年10月28日夏河5.7級(jí)等地震。ML3.6~4.0震級(jí)檔視應(yīng)力變化曲線(xiàn)在2013年底存在突升變化,2018年前后存在高值異常,兩震級(jí)檔視應(yīng)力均存在同步變化。

    3.3 青海北部視應(yīng)力的空間分布特征

    如圖5所示,ML3.0~3.5、ML3.6~4.0震級(jí)檔視應(yīng)力空間分布,兩個(gè)震級(jí)檔高視應(yīng)力異常區(qū)集中在柴達(dá)木中東部的柴達(dá)木北中央斷裂、柴達(dá)木北緣斷裂、宗務(wù)隆山斷裂、鄂拉山斷裂附近,在附近發(fā)生了大柴旦5.0級(jí)、烏蘭5.1級(jí)地震;青海東部的高視應(yīng)力異常區(qū)主要分布在貴德斷裂、拉脊山南緣斷裂、光蓋山—迭山北麓斷裂、托萊山斷裂等斷裂附近,在其附近發(fā)生了祁連5.2級(jí)、門(mén)源6.4級(jí)、甘州5.0級(jí)及夏河5.7級(jí)地震。在空間上,絕大部分5級(jí)地震發(fā)生在高視應(yīng)力異常區(qū)邊緣,也有些地震與異常在空間上存在一定距離。研究認(rèn)為高應(yīng)力集中區(qū)的應(yīng)力釋放后,會(huì)引起應(yīng)力快速調(diào)整并使其它地區(qū)的應(yīng)力水平快速增長(zhǎng)而引發(fā)地震,區(qū)域構(gòu)造應(yīng)力場(chǎng)快速調(diào)整,在應(yīng)力場(chǎng)調(diào)整范圍內(nèi)可能發(fā)生中強(qiáng)地震[12],這與本文得到結(jié)論基本一致,大部分中強(qiáng)地震前屬同構(gòu)造區(qū)域?qū)?yīng)有高視應(yīng)力異常。

    圖4 地震序列視應(yīng)力隨時(shí)間變化及M-T圖Fig.4 Variation of apparent stress with time and M-T diagram of the earthquake sequence

    圖5 研究區(qū)ML3.0~3.5、ML3.6~4.0震級(jí)檔地震視應(yīng)力空間分布Fig.5 Spatial distribution of seismic apparent stress of ML3.0~3.5 and ML3.6~4.0 earthquakes in the study area

    3.4 高視應(yīng)力異常的預(yù)測(cè)效能檢驗(yàn)

    本文針對(duì)青海北部ML3.0~4.9地震高視應(yīng)力異常與5級(jí)以上地震的對(duì)應(yīng)情況,進(jìn)行預(yù)測(cè)效能R值檢驗(yàn)。R=報(bào)對(duì)地震數(shù)/應(yīng)報(bào)地震數(shù)-預(yù)測(cè)占用時(shí)間/統(tǒng)計(jì)時(shí)間段總長(zhǎng)度,其中應(yīng)報(bào)地震數(shù)為發(fā)生在給定統(tǒng)計(jì)區(qū)域、給定統(tǒng)計(jì)時(shí)段內(nèi)滿(mǎn)足震級(jí)條件的強(qiáng)震總數(shù);報(bào)對(duì)地震數(shù)為統(tǒng)計(jì)區(qū)預(yù)測(cè)時(shí)間內(nèi)發(fā)生的強(qiáng)震總數(shù);預(yù)測(cè)占用時(shí)間為所有異常預(yù)測(cè)時(shí)間總和。只有當(dāng)R>R0時(shí),預(yù)測(cè)意義顯著,其中R0表示考慮報(bào)準(zhǔn)率和漏報(bào)率、置信度為97.5%時(shí)的最低R值[20]。

    計(jì)算得到研究區(qū)ML3.0~4.9地震視應(yīng)力平均值為0.14 MPa。本文分別以1 MPa、0.5 MPa和0.3 MPa為視應(yīng)力異常閾值進(jìn)行內(nèi)符檢驗(yàn)。

    (1) 預(yù)測(cè)規(guī)則1:以單次視應(yīng)力大于1 MPa的地震認(rèn)定為異常,共有12次地震視應(yīng)力大于1 MPa,以地震發(fā)生當(dāng)月為異常時(shí)間,共有9組異常,如表1所列。計(jì)算R值,得到的最佳預(yù)測(cè)時(shí)間為210天,R=0.55,R0=0.37。R>R0,通過(guò)顯著性檢驗(yàn),具有預(yù)測(cè)意義。

    表1 青海北部高視應(yīng)力中小地震與后續(xù)5級(jí)以上地震統(tǒng)計(jì)

    (2) 預(yù)測(cè)規(guī)則2:以單次視應(yīng)力大于0.5 MPa的地震認(rèn)定為異常,共33次地震視應(yīng)力大于0.5 MPa,無(wú)法通過(guò)檢驗(yàn)。為了提高異常的信度,對(duì)同一地點(diǎn)或鄰近區(qū)域連續(xù)發(fā)生3次以上視應(yīng)力大于0.5 MPa的地震定為一次異常,以最后一次地震為異常統(tǒng)計(jì)時(shí)間點(diǎn),共確認(rèn)5組異常,計(jì)算R值,得到的最佳預(yù)測(cè)時(shí)間為480天,R=0.45,R0=0.37。R>R0,通過(guò)統(tǒng)計(jì)檢驗(yàn),具有預(yù)測(cè)意義。

    (3) 預(yù)測(cè)規(guī)則3:以單次視應(yīng)力大于0.3 MPa地震為異常,共80次地震視應(yīng)力大于0.3 MPa,同樣無(wú)法通過(guò)檢驗(yàn)。統(tǒng)計(jì)發(fā)現(xiàn)高視應(yīng)力地震年頻次異常顯著,2013年、2014年、2018年及2019年發(fā)生的視應(yīng)力大于0.3 MPa地震頻次均超過(guò)10次,地震頻次明顯偏高,以青海北部發(fā)生視應(yīng)力大于0.3 MPa年頻次大于10次為異常,共確認(rèn)4組異常,計(jì)算R值,得到的最佳預(yù)測(cè)時(shí)間為300天,R=0.48,R0=0.38。R>R0,通過(guò)統(tǒng)計(jì)檢驗(yàn),具有預(yù)測(cè)意義。

    在青海北部地區(qū),發(fā)生視應(yīng)力大于1 MPa的中小地震或較大范圍高視應(yīng)力地震年頻次異常,其實(shí)就是中小地震反映了地殼應(yīng)力場(chǎng)前兆異常的表現(xiàn)。高視應(yīng)力異??赡苁钦鹪磮?chǎng)前兆,如2012年9月27日發(fā)生大柴旦ML4.3地震后,2013年在震中附近發(fā)生了5.0級(jí)地震;也有可能是是區(qū)域構(gòu)造場(chǎng)前兆,有些區(qū)域性質(zhì)的中小地震發(fā)生后,區(qū)域應(yīng)力場(chǎng)快速調(diào)整,在離震中較遠(yuǎn)的構(gòu)造相關(guān)區(qū)域發(fā)生中強(qiáng)地震,如2015年7月26日共和ML3.2地震后,在其東北方向約130 km的托萊山—冷龍嶺斷裂先后發(fā)生了2015年11月23日祁連5.2級(jí)地震和2016年1月21日門(mén)源6.4級(jí)地震。較大范圍高視應(yīng)力地震年頻次異常可能反映了地震增強(qiáng)活動(dòng)及區(qū)域構(gòu)造應(yīng)力場(chǎng)調(diào)整的過(guò)程,2013年、2014年、2018年及2019年高視應(yīng)力地震頻次偏高,5級(jí)地震也相對(duì)活躍。本文通過(guò)不同預(yù)測(cè)規(guī)則下的R值評(píng)分顯著性檢驗(yàn)結(jié)果顯示,高視應(yīng)力中小地震對(duì)研究區(qū)域未來(lái)一段時(shí)間發(fā)生中強(qiáng)地震具有一定預(yù)測(cè)意義,可以發(fā)現(xiàn)視應(yīng)力越高單次地震異常的可信度越高,同一地點(diǎn)發(fā)生多次高視應(yīng)力地震或高視應(yīng)力中小地震的年度高頻次異常,均有利于青海北部地區(qū)發(fā)生中強(qiáng)以上地震,預(yù)測(cè)的優(yōu)勢(shì)時(shí)段為半年至一年。

    4 結(jié)論與討論

    本文計(jì)算了2010年以來(lái)青海北部發(fā)生的ML≥3.0地震的震源參數(shù),通過(guò)分析震源參數(shù)間的標(biāo)度關(guān)系、ML3.0~3.5、ML3.6~4.0兩個(gè)震級(jí)檔視應(yīng)力的時(shí)空特征、高視應(yīng)力中小地震與中強(qiáng)震的對(duì)應(yīng)關(guān)系及預(yù)測(cè)效能檢驗(yàn),主要得出以下結(jié)論:

    (1) 通過(guò)分析震源參數(shù)的標(biāo)度關(guān)系,地震矩、地震輻射能量與震級(jí)存在較好的正相關(guān)性;拐角頻率、視應(yīng)力與震級(jí)也存在明顯的相關(guān)性,拐角頻率隨震級(jí)增大而減小,視應(yīng)力隨震級(jí)增大而增大。

    (2) 視應(yīng)力隨時(shí)間變化:ML3.0~3.5震級(jí)檔在2013—2015年及2018—2019年兩個(gè)時(shí)間段高視應(yīng)力異常顯著,ML3.6~4.0震級(jí)檔在2013年底和2018年前后存在高視應(yīng)力異常,該時(shí)間段5級(jí)地震相對(duì)活躍。

    (3) 分析視應(yīng)力的空間分布特征,高視應(yīng)力集中在柴達(dá)木盆地中東部的柴達(dá)木北中央斷裂、柴達(dá)木北緣斷裂、宗務(wù)隆山斷裂及鄂拉山斷裂附近,青海東部視應(yīng)力異常分布在貴德斷裂、拉脊山南緣斷裂、光蓋山—迭山北麓斷裂及托萊山斷裂等斷裂附近。大部分5級(jí)地震發(fā)生在異常區(qū)邊緣,這可能與應(yīng)力場(chǎng)調(diào)整后發(fā)生中強(qiáng)地震有關(guān)。

    (4) 本文通過(guò)不同預(yù)測(cè)規(guī)則的R值顯著性檢驗(yàn),結(jié)果顯示視應(yīng)力越高地震異常的可信度越高,相鄰區(qū)域發(fā)生多次高視應(yīng)力地震或整個(gè)區(qū)域高視應(yīng)力中小地震的年高頻次異常,均有利于青海北部地區(qū)發(fā)生中強(qiáng)以上地震,預(yù)測(cè)的優(yōu)勢(shì)時(shí)段為半年至一年。

    致謝:中國(guó)地震局地球物理研究所李艷娥副研究員提供了計(jì)算程序,在此表示感謝。

    猜你喜歡
    拐角震級(jí)應(yīng)力場(chǎng)
    拐 角
    基于累積絕對(duì)位移值的震級(jí)估算方法
    地震后各國(guó)發(fā)布的震級(jí)可能不一樣?
    Where Is My Home?
    新震級(jí)國(guó)家標(biāo)準(zhǔn)在大同臺(tái)的應(yīng)用與評(píng)估
    山西地震(2020年1期)2020-04-08 07:34:26
    走過(guò)那一個(gè)拐角
    美文(2017年4期)2017-02-23 14:26:12
    拐角遇到奇跡
    鋁合金多層多道窄間隙TIG焊接頭應(yīng)力場(chǎng)研究
    焊接(2016年9期)2016-02-27 13:05:22
    中國(guó)地震臺(tái)網(wǎng)面波震級(jí)與矩震級(jí)的統(tǒng)計(jì)關(guān)系
    考慮斷裂破碎帶的丹江口庫(kù)區(qū)地應(yīng)力場(chǎng)與水壓應(yīng)力場(chǎng)耦合反演及地震預(yù)測(cè)
    国产精品爽爽va在线观看网站| 成年版毛片免费区| 又黄又爽又免费观看的视频| 国产三级在线视频| 九九热线精品视视频播放| av视频在线观看入口| 可以在线观看的亚洲视频| 久久久久精品国产欧美久久久| 中文字幕高清在线视频| 国产一区二区三区视频了| 91在线观看av| 亚洲av电影不卡..在线观看| 国产欧美日韩精品一区二区| 又爽又黄a免费视频| 久久久久久久久久成人| 自拍偷自拍亚洲精品老妇| 十八禁网站免费在线| 搡老妇女老女人老熟妇| 国产精品伦人一区二区| 一个人看视频在线观看www免费| 欧美绝顶高潮抽搐喷水| 一区二区三区激情视频| 色吧在线观看| 一区二区三区免费毛片| 久久伊人香网站| 国产综合懂色| 色综合婷婷激情| 韩国av一区二区三区四区| 搞女人的毛片| 深夜a级毛片| 精品久久久久久久久av| 88av欧美| 在线观看午夜福利视频| 亚洲av成人精品一区久久| 日本成人三级电影网站| 欧美黄色片欧美黄色片| 99精品久久久久人妻精品| 精品一区二区三区av网在线观看| 亚洲成人免费电影在线观看| 久久人人爽人人爽人人片va | 国产视频内射| 99久久99久久久精品蜜桃| 亚洲 国产 在线| 精品久久久久久久久av| xxxwww97欧美| 欧美激情久久久久久爽电影| 男女做爰动态图高潮gif福利片| 久久6这里有精品| 欧美色欧美亚洲另类二区| 中文字幕高清在线视频| 久久国产乱子免费精品| 岛国在线免费视频观看| 尤物成人国产欧美一区二区三区| 亚洲人成网站在线播放欧美日韩| 亚洲av免费在线观看| 嫁个100分男人电影在线观看| 日本在线视频免费播放| 亚洲无线观看免费| 亚洲国产精品sss在线观看| 深爱激情五月婷婷| 亚洲av电影不卡..在线观看| av福利片在线观看| 深夜精品福利| 亚洲国产精品成人综合色| 国产日本99.免费观看| 亚洲中文日韩欧美视频| 久久国产乱子免费精品| 亚洲五月天丁香| 亚洲第一电影网av| 1000部很黄的大片| 俄罗斯特黄特色一大片| 亚洲欧美清纯卡通| 欧美xxxx性猛交bbbb| 中文字幕久久专区| 中文字幕av成人在线电影| 午夜福利免费观看在线| 91麻豆av在线| 亚洲中文字幕日韩| 国产一区二区亚洲精品在线观看| 久久99热6这里只有精品| 少妇的逼水好多| 老熟妇仑乱视频hdxx| 国产成人福利小说| 日本撒尿小便嘘嘘汇集6| 韩国av一区二区三区四区| 亚洲av成人不卡在线观看播放网| 欧美黄色片欧美黄色片| 国产中年淑女户外野战色| 天堂av国产一区二区熟女人妻| 国产爱豆传媒在线观看| 久久久久久大精品| 亚洲人成网站在线播| 白带黄色成豆腐渣| 日韩欧美三级三区| 国产欧美日韩一区二区精品| 欧美高清性xxxxhd video| 两性午夜刺激爽爽歪歪视频在线观看| 国产主播在线观看一区二区| 午夜福利视频1000在线观看| 人人妻人人看人人澡| 精品人妻1区二区| 亚洲最大成人av| 国产美女午夜福利| 蜜桃久久精品国产亚洲av| 在线播放国产精品三级| 国产伦一二天堂av在线观看| 麻豆成人av在线观看| 91午夜精品亚洲一区二区三区 | 国产伦精品一区二区三区视频9| 亚洲熟妇中文字幕五十中出| 一区二区三区四区激情视频 | 久久九九热精品免费| 久久性视频一级片| 国产真实乱freesex| 成年女人毛片免费观看观看9| 天堂影院成人在线观看| 久久久久久久久中文| av福利片在线观看| 国产av在哪里看| 亚洲av免费在线观看| 一级毛片久久久久久久久女| 欧美乱妇无乱码| 欧美成人免费av一区二区三区| 精品国内亚洲2022精品成人| 欧美不卡视频在线免费观看| 91久久精品电影网| 有码 亚洲区| 久久99热6这里只有精品| 亚洲欧美日韩东京热| 久久精品夜夜夜夜夜久久蜜豆| 午夜福利在线观看吧| 国产成人a区在线观看| 99久久久亚洲精品蜜臀av| 久久久精品欧美日韩精品| 亚洲第一欧美日韩一区二区三区| 国产精品久久久久久久电影| 深爱激情五月婷婷| 中文在线观看免费www的网站| 麻豆久久精品国产亚洲av| 欧美一区二区国产精品久久精品| 亚洲avbb在线观看| 久久精品91蜜桃| 有码 亚洲区| 国产亚洲欧美98| 国产精品国产高清国产av| 高清日韩中文字幕在线| 亚洲成a人片在线一区二区| 国内精品一区二区在线观看| 精品乱码久久久久久99久播| 亚洲一区高清亚洲精品| 久久精品91蜜桃| 午夜福利免费观看在线| av黄色大香蕉| 婷婷精品国产亚洲av在线| 一区二区三区四区激情视频 | 欧美日韩中文字幕国产精品一区二区三区| 99久国产av精品| 国产老妇女一区| 有码 亚洲区| 免费看日本二区| 国产视频一区二区在线看| 国产黄a三级三级三级人| 两人在一起打扑克的视频| 伦理电影大哥的女人| 色吧在线观看| 国产中年淑女户外野战色| 国产亚洲精品综合一区在线观看| 精品久久久久久久久亚洲 | 久久久精品大字幕| 久久国产乱子伦精品免费另类| 亚洲精品影视一区二区三区av| 国产成人a区在线观看| 99国产极品粉嫩在线观看| 丝袜美腿在线中文| 午夜日韩欧美国产| 宅男免费午夜| 国产v大片淫在线免费观看| 99视频精品全部免费 在线| 麻豆国产97在线/欧美| 亚洲精品456在线播放app | 一进一出抽搐gif免费好疼| 极品教师在线免费播放| 在线播放国产精品三级| av天堂在线播放| 亚洲精品影视一区二区三区av| 国产精品电影一区二区三区| 久久久成人免费电影| 99精品久久久久人妻精品| 天天一区二区日本电影三级| 国产一区二区在线av高清观看| 蜜桃亚洲精品一区二区三区| 久久热精品热| 国产精品女同一区二区软件 | 中文亚洲av片在线观看爽| 成年女人永久免费观看视频| 国内久久婷婷六月综合欲色啪| 一进一出抽搐gif免费好疼| 亚洲av第一区精品v没综合| 亚洲久久久久久中文字幕| 国内精品美女久久久久久| 国产精品女同一区二区软件 | 欧美日韩黄片免| 18禁在线播放成人免费| 国产色婷婷99| 久久这里只有精品中国| 欧美色欧美亚洲另类二区| 小说图片视频综合网站| 国产精品爽爽va在线观看网站| 亚洲精品乱码久久久v下载方式| 免费在线观看影片大全网站| 搡老妇女老女人老熟妇| 中文字幕高清在线视频| 日本精品一区二区三区蜜桃| 精品99又大又爽又粗少妇毛片 | 变态另类丝袜制服| 中文在线观看免费www的网站| av在线观看视频网站免费| 99riav亚洲国产免费| 日韩欧美在线乱码| 天堂√8在线中文| 婷婷丁香在线五月| 99国产极品粉嫩在线观看| 俺也久久电影网| 一本综合久久免费| 色哟哟·www| 蜜桃久久精品国产亚洲av| 在线观看av片永久免费下载| 毛片女人毛片| 在线十欧美十亚洲十日本专区| 欧美黑人欧美精品刺激| 日韩欧美在线二视频| 久久热精品热| 亚洲av免费高清在线观看| 熟妇人妻久久中文字幕3abv| 欧美性猛交黑人性爽| 色尼玛亚洲综合影院| 日韩欧美国产一区二区入口| 久久欧美精品欧美久久欧美| 90打野战视频偷拍视频| 日本一本二区三区精品| 男女下面进入的视频免费午夜| 能在线免费观看的黄片| 国产精品久久久久久久电影| 国产黄色小视频在线观看| 亚洲在线自拍视频| 久久99热这里只有精品18| 亚洲欧美激情综合另类| 久久热精品热| 真人一进一出gif抽搐免费| 国产精品一区二区三区四区免费观看 | 村上凉子中文字幕在线| 日本免费a在线| 久久草成人影院| 欧美激情久久久久久爽电影| 亚洲真实伦在线观看| 国产av不卡久久| 精品欧美国产一区二区三| 麻豆一二三区av精品| 日韩人妻高清精品专区| 午夜福利免费观看在线| 精品熟女少妇八av免费久了| 亚洲av成人精品一区久久| 亚洲中文字幕日韩| 免费高清视频大片| 热99re8久久精品国产| 少妇的逼水好多| 又紧又爽又黄一区二区| 日韩欧美三级三区| 18美女黄网站色大片免费观看| 日韩 亚洲 欧美在线| 少妇熟女aⅴ在线视频| 国产高清视频在线观看网站| 国产成人欧美在线观看| 国产视频一区二区在线看| 男女下面进入的视频免费午夜| 欧美色视频一区免费| 国产激情偷乱视频一区二区| 久久精品人妻少妇| 免费大片18禁| 国产精品一区二区免费欧美| 91麻豆精品激情在线观看国产| 极品教师在线视频| av中文乱码字幕在线| 麻豆国产av国片精品| 麻豆成人av在线观看| 国产高清有码在线观看视频| 久久精品国产亚洲av涩爱 | 国产精品不卡视频一区二区 | 欧美日韩国产亚洲二区| 亚洲av电影在线进入| 亚洲人成网站高清观看| 免费av毛片视频| 欧美日韩亚洲国产一区二区在线观看| 亚洲国产精品sss在线观看| 热99re8久久精品国产| 国产极品精品免费视频能看的| 亚洲aⅴ乱码一区二区在线播放| 9191精品国产免费久久| or卡值多少钱| 国产精品电影一区二区三区| 哪里可以看免费的av片| 亚洲av成人精品一区久久| 99久久精品热视频| 亚洲美女搞黄在线观看 | 91午夜精品亚洲一区二区三区 | 可以在线观看毛片的网站| 国产精品久久久久久精品电影| 国产精品久久久久久久久免 | 成人av在线播放网站| 国产老妇女一区| 免费搜索国产男女视频| 老司机深夜福利视频在线观看| 丰满乱子伦码专区| 国产精品野战在线观看| 国产一区二区亚洲精品在线观看| 国产伦精品一区二区三区四那| 深夜精品福利| 欧美一级a爱片免费观看看| 亚洲激情在线av| 亚洲男人的天堂狠狠| 亚洲精品日韩av片在线观看| 首页视频小说图片口味搜索| 亚洲精品一区av在线观看| 精品熟女少妇八av免费久了| 欧美一区二区精品小视频在线| 免费观看的影片在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av熟女| 熟妇人妻久久中文字幕3abv| 国产黄色小视频在线观看| 中文资源天堂在线| 美女xxoo啪啪120秒动态图 | 成人无遮挡网站| 成人高潮视频无遮挡免费网站| 美女xxoo啪啪120秒动态图 | 亚洲五月婷婷丁香| 亚洲中文字幕日韩| 欧美乱色亚洲激情| 69人妻影院| 看片在线看免费视频| 国产精品,欧美在线| 亚洲不卡免费看| 少妇丰满av| 国产精品久久久久久久久免 | 精品一区二区三区视频在线| 我的老师免费观看完整版| 青草久久国产| 日本黄大片高清| 青草久久国产| 国产高潮美女av| 精品久久久久久久末码| 国产亚洲精品综合一区在线观看| 久久99热这里只有精品18| 少妇人妻精品综合一区二区 | .国产精品久久| 丰满人妻一区二区三区视频av| 夜夜夜夜夜久久久久| 免费在线观看影片大全网站| 91久久精品国产一区二区成人| 国产欧美日韩精品一区二区| 久久久成人免费电影| 老司机午夜十八禁免费视频| 亚洲成人免费电影在线观看| 欧美成人性av电影在线观看| 国产精品不卡视频一区二区 | 国产激情偷乱视频一区二区| 91麻豆精品激情在线观看国产| 91麻豆av在线| 午夜福利在线在线| 日韩高清综合在线| 国内精品一区二区在线观看| 黄色视频,在线免费观看| 免费在线观看日本一区| 精品一区二区免费观看| 最新在线观看一区二区三区| 久久久久久九九精品二区国产| 中亚洲国语对白在线视频| 一个人看视频在线观看www免费| 国产精品免费一区二区三区在线| 蜜桃久久精品国产亚洲av| 国产真实伦视频高清在线观看 | 99热这里只有精品一区| 欧美日本亚洲视频在线播放| 1024手机看黄色片| 国产乱人伦免费视频| 俄罗斯特黄特色一大片| 日本一本二区三区精品| 国产蜜桃级精品一区二区三区| 女人被狂操c到高潮| 日韩免费av在线播放| 18美女黄网站色大片免费观看| 99热精品在线国产| 久久精品国产自在天天线| 欧美色欧美亚洲另类二区| 午夜影院日韩av| 性欧美人与动物交配| 97人妻精品一区二区三区麻豆| 精品午夜福利视频在线观看一区| 黄色女人牲交| 亚洲av成人精品一区久久| 久久人人爽人人爽人人片va | 在线观看一区二区三区| 尤物成人国产欧美一区二区三区| 亚洲18禁久久av| 免费人成在线观看视频色| 国产精品日韩av在线免费观看| 18+在线观看网站| 色哟哟·www| 亚洲午夜理论影院| 亚洲av熟女| 波野结衣二区三区在线| 日本一二三区视频观看| 成人特级黄色片久久久久久久| 51午夜福利影视在线观看| 国产精品嫩草影院av在线观看 | 国内揄拍国产精品人妻在线| 亚洲精品影视一区二区三区av| 中文字幕人妻熟人妻熟丝袜美| 亚洲真实伦在线观看| 国产亚洲精品久久久com| 在线播放无遮挡| 久久这里只有精品中国| 国产中年淑女户外野战色| 久久国产精品影院| 欧美性猛交黑人性爽| 一本精品99久久精品77| 亚洲国产精品sss在线观看| 亚洲黑人精品在线| 色综合婷婷激情| 在线免费观看的www视频| aaaaa片日本免费| 亚洲人成网站在线播放欧美日韩| 精品免费久久久久久久清纯| 人妻久久中文字幕网| 悠悠久久av| 免费在线观看成人毛片| 在线观看一区二区三区| 男女那种视频在线观看| 国产一区二区亚洲精品在线观看| 成人永久免费在线观看视频| 最近最新免费中文字幕在线| 国产极品精品免费视频能看的| 能在线免费观看的黄片| 国产午夜精品久久久久久一区二区三区 | 色噜噜av男人的天堂激情| 色尼玛亚洲综合影院| 97超视频在线观看视频| 看片在线看免费视频| 欧美日韩乱码在线| 国产老妇女一区| 久9热在线精品视频| 99国产极品粉嫩在线观看| 国产v大片淫在线免费观看| 色综合亚洲欧美另类图片| 天堂网av新在线| 亚洲av日韩精品久久久久久密| 特级一级黄色大片| 中国美女看黄片| 12—13女人毛片做爰片一| 久99久视频精品免费| 成人毛片a级毛片在线播放| 亚洲欧美日韩高清专用| 亚洲精品在线观看二区| 内地一区二区视频在线| 免费高清视频大片| 欧美激情国产日韩精品一区| 日韩欧美精品免费久久 | 在线免费观看不下载黄p国产 | 91在线精品国自产拍蜜月| 国产老妇女一区| 精品人妻熟女av久视频| 91av网一区二区| 午夜久久久久精精品| 国产一区二区激情短视频| 久久精品国产亚洲av涩爱 | 波多野结衣高清无吗| 蜜桃亚洲精品一区二区三区| 免费电影在线观看免费观看| 麻豆一二三区av精品| 亚洲七黄色美女视频| 国产精品亚洲一级av第二区| 在线免费观看不下载黄p国产 | 精品久久久久久久末码| 午夜福利高清视频| 老司机午夜福利在线观看视频| 99在线人妻在线中文字幕| 午夜免费男女啪啪视频观看 | 噜噜噜噜噜久久久久久91| 又黄又爽又免费观看的视频| 精品欧美国产一区二区三| 日本与韩国留学比较| 日本三级黄在线观看| 在线免费观看的www视频| 欧美高清成人免费视频www| 91久久精品国产一区二区成人| 欧美日韩黄片免| 日本黄色片子视频| 免费电影在线观看免费观看| 国产精品不卡视频一区二区 | 久久热精品热| 桃色一区二区三区在线观看| 日本与韩国留学比较| 校园春色视频在线观看| 欧美精品国产亚洲| 欧美+日韩+精品| 18美女黄网站色大片免费观看| 日韩人妻高清精品专区| 国产精品久久电影中文字幕| 国产黄a三级三级三级人| 亚洲最大成人中文| 好男人在线观看高清免费视频| 精品久久久久久,| 一本精品99久久精品77| 乱人视频在线观看| 99视频精品全部免费 在线| 精品无人区乱码1区二区| 国产主播在线观看一区二区| 国产激情偷乱视频一区二区| 夜夜爽天天搞| 欧美+亚洲+日韩+国产| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久久午夜电影| 网址你懂的国产日韩在线| 成人性生交大片免费视频hd| 欧美bdsm另类| 国产精品野战在线观看| 男女视频在线观看网站免费| 又粗又爽又猛毛片免费看| 亚洲精品乱码久久久v下载方式| 麻豆国产97在线/欧美| 国产精品一区二区三区四区免费观看 | 久久久久久久久大av| 亚洲av中文字字幕乱码综合| 亚洲五月天丁香| 男人舔奶头视频| 亚洲第一电影网av| 麻豆国产97在线/欧美| 色综合婷婷激情| 日本免费a在线| 人妻制服诱惑在线中文字幕| 狠狠狠狠99中文字幕| 国产黄片美女视频| а√天堂www在线а√下载| 白带黄色成豆腐渣| 日韩高清综合在线| 每晚都被弄得嗷嗷叫到高潮| 欧美色视频一区免费| 欧美一区二区精品小视频在线| 简卡轻食公司| 久久精品人妻少妇| 色视频www国产| 中文字幕人成人乱码亚洲影| 国产精品久久电影中文字幕| 亚洲综合色惰| 久久久精品大字幕| 亚洲国产欧美人成| 亚洲无线观看免费| 岛国在线免费视频观看| 变态另类成人亚洲欧美熟女| 国产精品嫩草影院av在线观看 | 国产精品98久久久久久宅男小说| 桃红色精品国产亚洲av| 国产免费一级a男人的天堂| 午夜激情欧美在线| 村上凉子中文字幕在线| 中文字幕免费在线视频6| 日韩欧美一区二区三区在线观看| 99热这里只有是精品50| 国产亚洲欧美在线一区二区| 91狼人影院| av视频在线观看入口| 又紧又爽又黄一区二区| 麻豆国产av国片精品| 久久久久亚洲av毛片大全| 久久草成人影院| 精品一区二区免费观看| 国内精品久久久久精免费| 免费高清视频大片| 中文资源天堂在线| 亚洲黑人精品在线| 一区二区三区激情视频| 日本黄色片子视频| 天堂√8在线中文| 国产精品98久久久久久宅男小说| 亚洲成a人片在线一区二区| 国产高清视频在线观看网站| 婷婷亚洲欧美| 在线看三级毛片| 97超视频在线观看视频| 欧美日韩瑟瑟在线播放| 国产成人av教育| 国产精品免费一区二区三区在线| 亚洲最大成人中文| 国产麻豆成人av免费视频| 亚洲,欧美精品.| 欧美精品啪啪一区二区三区| 午夜激情欧美在线| 美女高潮喷水抽搐中文字幕| 国产精品一区二区三区四区免费观看 | 成年女人永久免费观看视频| 欧美日韩综合久久久久久 | 天堂网av新在线| 一个人观看的视频www高清免费观看| 午夜福利在线在线| 亚洲在线观看片| 欧美黑人巨大hd| 国产av在哪里看| 欧美最黄视频在线播放免费| 国模一区二区三区四区视频| 色在线成人网| 又爽又黄a免费视频| 成人特级av手机在线观看| 国产黄片美女视频| 欧美色欧美亚洲另类二区| АⅤ资源中文在线天堂|