• <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è)
    在线看a的网站| 免费观看av网站的网址| 18禁观看日本| 久久久精品国产亚洲av高清涩受| 老汉色∧v一级毛片| 777米奇影视久久| 电影成人av| 久久免费观看电影| 婷婷色综合大香蕉| 哪个播放器可以免费观看大片| 午夜福利影视在线免费观看| 91aial.com中文字幕在线观看| 两个人免费观看高清视频| 国产精品香港三级国产av潘金莲 | 最新的欧美精品一区二区| 99久久人妻综合| 男人爽女人下面视频在线观看| 亚洲 欧美一区二区三区| 婷婷色麻豆天堂久久| 青春草视频在线免费观看| 国产免费视频播放在线视频| 卡戴珊不雅视频在线播放| 亚洲三区欧美一区| 99久久人妻综合| 亚洲精品一区蜜桃| 国产日韩一区二区三区精品不卡| 国产精品 国内视频| 精品国产一区二区三区四区第35| 亚洲国产成人一精品久久久| 亚洲成av片中文字幕在线观看 | 久久免费观看电影| 欧美日韩亚洲高清精品| av天堂久久9| 十分钟在线观看高清视频www| 青草久久国产| 精品一区二区三区四区五区乱码 | 高清黄色对白视频在线免费看| 婷婷成人精品国产| 一二三四在线观看免费中文在| 天天躁夜夜躁狠狠躁躁| 免费在线观看完整版高清| 午夜日本视频在线| 国产97色在线日韩免费| 99九九在线精品视频| 国产欧美日韩一区二区三区在线| 精品人妻熟女毛片av久久网站| freevideosex欧美| www.熟女人妻精品国产| 久久青草综合色| 亚洲精品乱久久久久久| 91精品三级在线观看| 国产又爽黄色视频| 美女午夜性视频免费| www.av在线官网国产| 免费女性裸体啪啪无遮挡网站| 婷婷色麻豆天堂久久| 街头女战士在线观看网站| 日韩中文字幕视频在线看片| 亚洲精品国产色婷婷电影| 大片免费播放器 马上看| 亚洲欧美精品综合一区二区三区 | 亚洲人成77777在线视频| 亚洲精华国产精华液的使用体验| 最近中文字幕2019免费版| 在线观看免费日韩欧美大片| 中文字幕色久视频| 国产精品女同一区二区软件| 不卡av一区二区三区| 国产成人aa在线观看| 一二三四在线观看免费中文在| 国产欧美日韩一区二区三区在线| 日韩精品免费视频一区二区三区| 大香蕉久久网| 国产高清不卡午夜福利| 久久久久久久久久久久大奶| 如日韩欧美国产精品一区二区三区| 看非洲黑人一级黄片| 久久久久久伊人网av| av免费在线看不卡| 国产成人免费观看mmmm| 18禁裸乳无遮挡动漫免费视频| 国产男人的电影天堂91| 一边亲一边摸免费视频| 国产综合精华液| 久久久久久免费高清国产稀缺| 亚洲精华国产精华液的使用体验| 制服丝袜香蕉在线| 汤姆久久久久久久影院中文字幕| 日本欧美国产在线视频| 99热全是精品| 婷婷色综合www| 久久精品久久久久久久性| 日韩中文字幕欧美一区二区 | 在线精品无人区一区二区三| 波野结衣二区三区在线| 久久精品亚洲av国产电影网| 亚洲成色77777| 一本大道久久a久久精品| 欧美 亚洲 国产 日韩一| 久久久国产欧美日韩av| 亚洲欧美一区二区三区黑人 | 天天躁夜夜躁狠狠躁躁| 国产av精品麻豆| 中文字幕色久视频| 欧美日韩精品网址| 久热这里只有精品99| 超色免费av| 国产综合精华液| 18禁观看日本| 国产有黄有色有爽视频| 亚洲伊人久久精品综合| 国产成人a∨麻豆精品| 国产亚洲精品第一综合不卡| 国产成人精品无人区| 可以免费在线观看a视频的电影网站 | 中文字幕人妻丝袜一区二区 | 国产精品国产三级专区第一集| 免费人妻精品一区二区三区视频| 亚洲欧洲日产国产| 午夜免费鲁丝| 国产精品久久久久久精品电影小说| 韩国av在线不卡| 国产av一区二区精品久久| 叶爱在线成人免费视频播放| freevideosex欧美| 日韩大片免费观看网站| 黑人欧美特级aaaaaa片| 国语对白做爰xxxⅹ性视频网站| 热re99久久精品国产66热6| 永久免费av网站大全| 婷婷色麻豆天堂久久| 一级片'在线观看视频| 高清不卡的av网站| 热99国产精品久久久久久7| 久久久精品区二区三区| 久久婷婷青草| 日韩在线高清观看一区二区三区| 黄色配什么色好看| 美女脱内裤让男人舔精品视频| 天堂中文最新版在线下载| 国产熟女午夜一区二区三区| 国产精品免费大片| 国产1区2区3区精品| 日韩免费高清中文字幕av| 人妻系列 视频| 一区二区三区激情视频| 国产成人精品一,二区| 国产成人精品一,二区| 国产 一区精品| 日韩中文字幕视频在线看片| 日本欧美国产在线视频| 午夜免费观看性视频| 1024视频免费在线观看| 欧美少妇被猛烈插入视频| 丰满饥渴人妻一区二区三| 伦精品一区二区三区| 免费在线观看黄色视频的| 亚洲国产最新在线播放| 国产精品麻豆人妻色哟哟久久| 久久精品国产亚洲av高清一级| 熟女少妇亚洲综合色aaa.| 777久久人妻少妇嫩草av网站| a级片在线免费高清观看视频| 久久99蜜桃精品久久| 日本爱情动作片www.在线观看| 国产精品久久久久久久久免| 日日啪夜夜爽| 免费看不卡的av| 日韩一本色道免费dvd| 观看美女的网站| 天堂俺去俺来也www色官网| 国产av一区二区精品久久| 在线观看免费日韩欧美大片| 久久久国产一区二区| 亚洲欧美精品综合一区二区三区 | 2022亚洲国产成人精品| 午夜福利视频精品| 男人添女人高潮全过程视频| 只有这里有精品99| 嫩草影院入口| 亚洲精品久久成人aⅴ小说| 国产黄频视频在线观看| 99久久综合免费| 日韩电影二区| 99香蕉大伊视频| 免费在线观看视频国产中文字幕亚洲 | 满18在线观看网站| 国产一区二区 视频在线| 日韩制服丝袜自拍偷拍| 大话2 男鬼变身卡| 久久久久久久久免费视频了| 久久精品国产鲁丝片午夜精品| 亚洲av欧美aⅴ国产| 久久精品久久精品一区二区三区| 久久久久久久亚洲中文字幕| 国产精品 国内视频| 啦啦啦视频在线资源免费观看| 少妇精品久久久久久久| 亚洲国产精品一区二区三区在线| 国产成人av激情在线播放| 精品人妻一区二区三区麻豆| 最新的欧美精品一区二区| 色播在线永久视频| 免费看不卡的av| 美女国产视频在线观看| 久久久久国产精品人妻一区二区| 美女高潮到喷水免费观看| 中文字幕色久视频| 丝袜在线中文字幕| av.在线天堂| 我的亚洲天堂| 欧美av亚洲av综合av国产av | 久久青草综合色| 99久久人妻综合| 国产精品一区二区在线观看99| 久久韩国三级中文字幕| 日韩欧美精品免费久久| www.自偷自拍.com| 天天躁夜夜躁狠狠久久av| 亚洲精品av麻豆狂野| 日本欧美视频一区| 看非洲黑人一级黄片| 黄色视频在线播放观看不卡| 亚洲人成77777在线视频| 亚洲精品国产色婷婷电影| 女人高潮潮喷娇喘18禁视频| 丝袜人妻中文字幕| 99九九在线精品视频| 亚洲国产精品999| 可以免费在线观看a视频的电影网站 | 国产成人aa在线观看| 午夜福利在线免费观看网站| 久久ye,这里只有精品| 国产无遮挡羞羞视频在线观看| 久久午夜综合久久蜜桃| 只有这里有精品99| 国产激情久久老熟女| 国产黄色视频一区二区在线观看| 亚洲精品美女久久久久99蜜臀 | 久久久久久久大尺度免费视频| 黑丝袜美女国产一区| 久久久精品94久久精品| 中文字幕亚洲精品专区| 香蕉丝袜av| 日韩av免费高清视频| 美女国产高潮福利片在线看| 巨乳人妻的诱惑在线观看| 日韩视频在线欧美| 黄色一级大片看看| 精品一区二区免费观看| 麻豆乱淫一区二区| 中文字幕另类日韩欧美亚洲嫩草| 国产高清国产精品国产三级| 久久久精品区二区三区| 1024香蕉在线观看| 少妇人妻久久综合中文| 国产97色在线日韩免费| 亚洲婷婷狠狠爱综合网| 免费少妇av软件| 又大又黄又爽视频免费| 欧美精品亚洲一区二区| 黄色配什么色好看| 满18在线观看网站| 建设人人有责人人尽责人人享有的| 久久狼人影院| 永久免费av网站大全| 久久人人爽av亚洲精品天堂| 久久久国产欧美日韩av| 男人操女人黄网站| 欧美精品亚洲一区二区| 肉色欧美久久久久久久蜜桃| 久久97久久精品| 亚洲视频免费观看视频| 大片电影免费在线观看免费| 伊人久久国产一区二区| 青草久久国产| 最近2019中文字幕mv第一页| 女性生殖器流出的白浆| 制服人妻中文乱码| 精品国产一区二区三区四区第35| 日韩中字成人| 久久久精品区二区三区| 纵有疾风起免费观看全集完整版| 午夜日韩欧美国产| 国产精品二区激情视频| 有码 亚洲区| 日本91视频免费播放| 国产熟女午夜一区二区三区| 亚洲欧美一区二区三区国产| 永久网站在线| 亚洲欧美一区二区三区久久| 日本黄色日本黄色录像| 青春草亚洲视频在线观看| 国产又爽黄色视频| 国产成人免费观看mmmm| 男人舔女人的私密视频| 欧美日韩视频精品一区| 亚洲精品第二区| 国产福利在线免费观看视频| 亚洲精品国产av成人精品| 在线精品无人区一区二区三| 中文字幕最新亚洲高清| 午夜免费观看性视频| 国产精品人妻久久久影院| 国产欧美日韩一区二区三区在线| 国产精品成人在线| 欧美国产精品va在线观看不卡| 午夜免费男女啪啪视频观看| 免费观看a级毛片全部| 日本vs欧美在线观看视频| 多毛熟女@视频| 晚上一个人看的免费电影| 日本免费在线观看一区| 久久久久国产一级毛片高清牌| 少妇熟女欧美另类| videosex国产| 欧美中文综合在线视频| 女的被弄到高潮叫床怎么办| 人体艺术视频欧美日本| 黑人猛操日本美女一级片| 久久久久久人人人人人| 肉色欧美久久久久久久蜜桃| 成人18禁高潮啪啪吃奶动态图| 老熟女久久久| 边亲边吃奶的免费视频| 亚洲色图综合在线观看| 美女午夜性视频免费| 国产精品一国产av| 丰满少妇做爰视频| 午夜福利一区二区在线看| 欧美bdsm另类| 亚洲欧美中文字幕日韩二区| 亚洲成av片中文字幕在线观看 | 在线观看美女被高潮喷水网站| 涩涩av久久男人的天堂| 国产极品天堂在线| 丰满少妇做爰视频| 日韩在线高清观看一区二区三区| 久久这里只有精品19| 97精品久久久久久久久久精品| 看免费成人av毛片| 国产一区亚洲一区在线观看| 成人二区视频| 满18在线观看网站| 国精品久久久久久国模美| 精品视频人人做人人爽| 日韩伦理黄色片| 亚洲欧美中文字幕日韩二区| 国产成人精品福利久久| 国产男女内射视频| 女人被躁到高潮嗷嗷叫费观| 男女免费视频国产| 黄色怎么调成土黄色| 亚洲精品一区蜜桃| 亚洲av中文av极速乱| 午夜免费鲁丝| 免费少妇av软件| 欧美日韩视频精品一区| 春色校园在线视频观看| 久久国产精品男人的天堂亚洲| 亚洲激情五月婷婷啪啪| 黄色毛片三级朝国网站| 美女中出高潮动态图| 午夜av观看不卡| 一本久久精品| 国产一区二区三区综合在线观看| 亚洲情色 制服丝袜| 亚洲精华国产精华液的使用体验| 国产精品一区二区在线不卡| 亚洲精华国产精华液的使用体验| 久久午夜综合久久蜜桃| 国产伦理片在线播放av一区| 制服人妻中文乱码| 深夜精品福利| 丰满迷人的少妇在线观看| 久久亚洲国产成人精品v| 久久热在线av| 精品人妻熟女毛片av久久网站| xxxhd国产人妻xxx| 久久鲁丝午夜福利片| 日日摸夜夜添夜夜爱| 看免费成人av毛片| 国产成人a∨麻豆精品| 90打野战视频偷拍视频| 国产精品久久久久久av不卡| 国产欧美亚洲国产| 男男h啪啪无遮挡| 午夜福利在线免费观看网站| 精品久久久精品久久久| 亚洲第一区二区三区不卡| 在线 av 中文字幕| 午夜福利网站1000一区二区三区| 男人爽女人下面视频在线观看| 亚洲av欧美aⅴ国产| 一级a爱视频在线免费观看| 亚洲欧美精品自产自拍| 人妻少妇偷人精品九色| 18在线观看网站| 高清欧美精品videossex| 久久鲁丝午夜福利片| 有码 亚洲区| 亚洲天堂av无毛| 9191精品国产免费久久| 国产又色又爽无遮挡免| 亚洲av福利一区| 国产精品久久久久久久久免| 纵有疾风起免费观看全集完整版| 国产日韩一区二区三区精品不卡| videos熟女内射| 两性夫妻黄色片| 亚洲av免费高清在线观看| 亚洲av日韩在线播放| 美女国产高潮福利片在线看| 1024视频免费在线观看| 免费看av在线观看网站| 免费播放大片免费观看视频在线观看| 日韩免费高清中文字幕av| 看十八女毛片水多多多| 人妻人人澡人人爽人人| 永久免费av网站大全| 免费黄色在线免费观看| 国产 一区精品| 自线自在国产av| 色吧在线观看| 一级爰片在线观看| 免费久久久久久久精品成人欧美视频| 亚洲欧美成人精品一区二区| 色视频在线一区二区三区| 日韩熟女老妇一区二区性免费视频| 天天躁日日躁夜夜躁夜夜| 中文字幕人妻丝袜制服| 97人妻天天添夜夜摸| 亚洲欧美一区二区三区黑人 | 男女午夜视频在线观看| 亚洲精品一二三| 日韩制服骚丝袜av| 热re99久久精品国产66热6| 国产黄频视频在线观看| 久久人人97超碰香蕉20202| 亚洲伊人久久精品综合| 男人添女人高潮全过程视频| 久久午夜综合久久蜜桃| 五月伊人婷婷丁香| 欧美日韩亚洲高清精品| 美女国产视频在线观看| 青春草视频在线免费观看| 国产片内射在线| 久久av网站| 在线观看免费高清a一片| 日本免费在线观看一区| 日韩 亚洲 欧美在线| 欧美精品一区二区免费开放| 老熟女久久久| 国产男女内射视频| 深夜精品福利| 精品久久久精品久久久| 成人毛片a级毛片在线播放| 天堂俺去俺来也www色官网| 免费在线观看黄色视频的| 久久久a久久爽久久v久久| 国产男女内射视频| 高清视频免费观看一区二区| 亚洲国产最新在线播放| 伊人亚洲综合成人网| 亚洲精品在线美女| 成人国产麻豆网| 国产精品麻豆人妻色哟哟久久| 丰满乱子伦码专区| 看免费成人av毛片| 我要看黄色一级片免费的| 1024视频免费在线观看| 女的被弄到高潮叫床怎么办| 亚洲精品乱久久久久久| 久久久欧美国产精品| 日日爽夜夜爽网站| 免费观看在线日韩| 国产乱来视频区| 欧美日韩一级在线毛片| 丁香六月天网| 久久午夜综合久久蜜桃| 久久精品人人爽人人爽视色| 欧美日韩av久久| 国产 一区精品| videos熟女内射| 男女高潮啪啪啪动态图| 黄色配什么色好看| 国产精品久久久久成人av| 91午夜精品亚洲一区二区三区| 蜜桃国产av成人99| 十分钟在线观看高清视频www| av视频免费观看在线观看| 久久青草综合色| 免费观看av网站的网址| 国产乱来视频区| 日韩制服骚丝袜av| 两个人免费观看高清视频| 国产av国产精品国产| 国产av码专区亚洲av| 亚洲五月色婷婷综合| 欧美bdsm另类| 精品一品国产午夜福利视频| 中文字幕av电影在线播放| 中文字幕另类日韩欧美亚洲嫩草| 麻豆乱淫一区二区| 国产麻豆69| 美女午夜性视频免费| 我的亚洲天堂| 欧美xxⅹ黑人| 美女视频免费永久观看网站| xxx大片免费视频| 亚洲欧美成人精品一区二区| 亚洲欧美一区二区三区久久| 国产极品天堂在线| 啦啦啦啦在线视频资源| 欧美精品一区二区免费开放| 国产精品久久久久久精品古装| 日本免费在线观看一区| 日韩伦理黄色片| 中文字幕亚洲精品专区| 性少妇av在线| 免费少妇av软件| 伦理电影免费视频| 精品少妇黑人巨大在线播放| 国产免费视频播放在线视频| 一级毛片黄色毛片免费观看视频| 亚洲av免费高清在线观看| 精品人妻偷拍中文字幕| 午夜91福利影院| 满18在线观看网站| 欧美日韩av久久| 大话2 男鬼变身卡| 一级毛片我不卡| 亚洲经典国产精华液单| 国产在视频线精品| 久久99一区二区三区| 亚洲四区av| 久久精品国产亚洲av高清一级| 欧美成人午夜精品| 久久99一区二区三区| 春色校园在线视频观看| 成人毛片60女人毛片免费| 黄网站色视频无遮挡免费观看| 香蕉国产在线看| 老熟女久久久| 亚洲综合色网址| 一区福利在线观看| 中文字幕最新亚洲高清| 国产成人一区二区在线| av国产精品久久久久影院| 日韩熟女老妇一区二区性免费视频| 丰满少妇做爰视频| 交换朋友夫妻互换小说| 国产日韩欧美视频二区| 免费观看a级毛片全部| av在线app专区| av网站在线播放免费| 久久久欧美国产精品| 91在线精品国自产拍蜜月| 久久久久久久久久久免费av| 欧美人与性动交α欧美软件| 国产免费福利视频在线观看| 最近的中文字幕免费完整| 日韩av不卡免费在线播放| 热re99久久精品国产66热6| 人人妻人人澡人人看| tube8黄色片| 午夜久久久在线观看| 精品一品国产午夜福利视频| 日本猛色少妇xxxxx猛交久久| 国产精品国产av在线观看| 啦啦啦中文免费视频观看日本| 欧美日韩综合久久久久久| 国产片特级美女逼逼视频| 亚洲精品,欧美精品| 久久人人97超碰香蕉20202| 爱豆传媒免费全集在线观看| 在线观看免费视频网站a站| 国产亚洲欧美精品永久| 亚洲欧美成人精品一区二区| 9色porny在线观看| 亚洲精品美女久久久久99蜜臀 | 久久久a久久爽久久v久久| 中文字幕色久视频| 欧美黄色片欧美黄色片| 精品少妇内射三级| 午夜福利网站1000一区二区三区| 免费看av在线观看网站| 少妇人妻久久综合中文| 成人午夜精彩视频在线观看| 亚洲 欧美一区二区三区| av.在线天堂| 久久人人爽人人片av| 国产熟女午夜一区二区三区| 日本91视频免费播放| 99热网站在线观看| 久热这里只有精品99| 黑人猛操日本美女一级片| 亚洲中文av在线| 不卡av一区二区三区| 男女啪啪激烈高潮av片| 久久久久久久大尺度免费视频| 爱豆传媒免费全集在线观看| 国产又色又爽无遮挡免| 男男h啪啪无遮挡| 成人毛片60女人毛片免费| 亚洲av综合色区一区| 99久久中文字幕三级久久日本| 日韩电影二区| 美女xxoo啪啪120秒动态图| 午夜免费男女啪啪视频观看| 国产综合精华液| 黑丝袜美女国产一区| 欧美精品亚洲一区二区| 男人操女人黄网站| 国产免费又黄又爽又色|