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

    基于重力觀測的巴顏喀拉塊體強(qiáng)震震級模擬分析*

    2021-12-10 05:17:54楊彥明張小艷王旭東楊紅櫻
    地震科學(xué)進(jìn)展 2021年11期
    關(guān)鍵詞:同震重力儀面波

    熊 峰 楊彥明※ 龐 博 張小艷 王旭東 楊紅櫻

    1) 內(nèi)蒙古自治區(qū)地震局,呼和浩特 010010

    2) 北京交通大學(xué),北京 100044

    引言

    重力觀測是研究地球內(nèi)部構(gòu)造的重要手段,其固體潮連續(xù)觀測數(shù)據(jù)表征了與潮汐變化有關(guān)的地殼運(yùn)動,以及由于地球各圈層物質(zhì)遷移而導(dǎo)致的重力變化,也是研究地殼構(gòu)造運(yùn)動的重要物理量[1]。地震孕育過程中,地球系統(tǒng)將發(fā)生一系列的物理化學(xué)變化,通過對重力固體潮曲線進(jìn)行分析,提取數(shù)據(jù)中的非潮汐變化,進(jìn)一步識別關(guān)于地球物理的地震異常信息[2]。重力對于同震變化具有敏感性,通過對同震形變波的定量研究,可以估計出地震震源的基本參數(shù)[3-6]。由于地震導(dǎo)致的地面運(yùn)動的唯一性,地震儀器記錄地面運(yùn)動的位移、速度以及加速度值。重力儀與地震儀類似,表征地面運(yùn)動的加速度,能夠記錄到地震破裂和傳播過程中所產(chǎn)生的地震波,對1 s以上的地震波信號均有響應(yīng),對遠(yuǎn)震波形記錄有較好的一致性。利用重力儀可記錄到P波、S波及面波的特性,對比分析PET相對重力儀與JCZ-1甚寬帶地震儀記錄到的地震波信息,可以揭示同震形變波形中蘊(yùn)含地球動力特征的物理機(jī)制。

    巴顏喀拉塊體為青藏高原現(xiàn)今地殼運(yùn)動最為強(qiáng)烈的地區(qū)之一,周緣斷裂活動強(qiáng)烈,強(qiáng)震活動頻繁,我國大陸7級以上強(qiáng)震幾乎都發(fā)生在巴顏喀拉塊體周緣斷裂上[7-11]。內(nèi)蒙古烏加河地震臺PET相對重力儀于2008年正式運(yùn)行,清晰完整地記錄到2008年3月21日新疆于田MS7.3、2010年4月14日青海玉樹MS7.1、2013年4月20日四川蘆山MS7.0、2014年2月12日新疆于田MS7.3、2017年8月8日四川九寨溝MS7.0、2021年5月22日青?,敹郙S7.4地震,全部發(fā)生在巴顏喀拉塊體周緣或內(nèi)部斷裂上,震中距為1000—2500 km。本文通過計算2020年以前發(fā)生在巴顏喀拉塊體周圍的5次7.0級以上地震,模擬分析得出一般震級模型、地應(yīng)力波模型、S波P波幅度比模型,理論震級與實(shí)測震級相關(guān)系數(shù)R值為0.90—0.99。2021年5月22日青?,敹喟l(fā)生MS7.4地震后,再次對比3種震級模型的差異,計算理論震級與實(shí)際震級的誤差及R值,進(jìn)一步驗證模型的可靠性。通過對同震形變波的定量研究,可認(rèn)識到在同一力源作用下,同震波形特征的差異性,這對理解形變前兆觀測的物理本質(zhì)是有意義的[3,12]。利用測震儀與重力儀震例回溯、檢驗,進(jìn)而探討地震的破裂與傳播過程,促進(jìn)地震學(xué)與重力學(xué)學(xué)科之間的相互借鑒與融合,揭示同震形變波形蘊(yùn)含地球動力特征的物理機(jī)制,挖掘能夠反映強(qiáng)地震發(fā)生的典型要素。

    1 資料基本情況

    1.1 數(shù)據(jù)與臺站

    烏加河地震臺地處鄂爾多斯北緣,位于陰山和狼山構(gòu)造帶交匯地區(qū)的河套斷陷帶。觀測山洞為花崗巖,巖體完整,巖性致密均勻,數(shù)據(jù)質(zhì)量好、信噪比高。山洞布設(shè)寬頻帶數(shù)字地震儀和PET相對重力儀。該重力儀分辨率為0.01×10?8m/s2,精度為1×10?8m/s2,數(shù)據(jù)采樣率達(dá)到1 Hz,可記錄到固體潮的微小變化和地震波[2,11]。

    本次研究對象為2008—2021年中國大陸發(fā)生的6次7.0級以上地震,烏加河地震臺PET相對重力儀和寬頻帶數(shù)字地震儀都記錄到完整的地震波數(shù)據(jù)。圖1為臺站與地震分布圖。

    圖1 臺站與地震分布圖Fig. 1 Distribution map of seismic station and earthquakes

    1.2 數(shù)據(jù)處理

    圖2給出了2017年8月8日21:19:46烏加河地震臺JCZ-1甚寬帶地震儀和PET相對重力儀記錄到的四川阿壩九寨溝MS7.0地震波形。此次地震震中位于巴顏喀拉地塊的東北部、東昆侖斷裂帶東段塔藏斷裂附近[7],距離烏加河地震臺1 013.18 km。由于重力儀記錄的是加速度,將地震儀記錄的速度值求導(dǎo)后,得到加速度數(shù)據(jù),將速度、加速度與重力儀記錄的秒采樣重力殘差進(jìn)行對比分析,研究發(fā)現(xiàn)二者相關(guān)系數(shù)均值最高可達(dá)0.9996以上[1]。重力儀和地震儀記錄的速度值求導(dǎo)后,數(shù)據(jù)波形記錄的P波初動方向、P波與S波到時具有較高的一致性,同時在持續(xù)時間和振幅等震相特征上也高度相似。在形變觀測研究中,該研究成果為識別地震異常提供了重要參考,同時,在重力日常觀測和數(shù)據(jù)處理工作中,對波形出現(xiàn)階變、突跳等異常變化時的性質(zhì)判定提供理論標(biāo)準(zhǔn)[2]。

    圖2 地震儀與重力儀地震波形對比Fig. 2 Comparison of seismograph and gravimeter seismic waveforms

    2 同震波形分析

    根據(jù)地震永久形變引起的變化理論,對同震響應(yīng)以及震后變化的物理機(jī)制進(jìn)行分析。地震類型和震級的不同,計算和分析所對應(yīng)的同震響應(yīng),可以進(jìn)一步確定由于破壞性地震引起的形變同震響應(yīng)以及震后影響的空間特征[13]。由于同震波形中攜帶了地震破裂以及地震波傳播的信息,因此,利用同震信息能夠評估地震震源參數(shù),進(jìn)一步分析同震波形的物理特征,將有助于提高對短臨前兆動力特性的認(rèn)識[14-15]。

    2.1 同震形變波的基本物理特征

    受地震波激發(fā),傾斜和重力等形變儀器均可記錄到地震的同震響應(yīng),其特征主要為面波延遲時間、振幅及震時形變波持續(xù)時間等[16]。地震面波主要有Rayleigh(瑞雷波)和Love(勒夫波)兩種類型。Love在水平方向運(yùn)動,是SH型的面波;Rayleigh的質(zhì)點(diǎn)運(yùn)動是一個橢圓。面波能量強(qiáng),只局限在自由表面附近擾動,其傳播距離較遠(yuǎn)。通過對面波進(jìn)行分析,可測定震級,利用面波最大振幅比、震時形變波的同震形變階躍均可確定震源機(jī)制解[6,17]。因此,震時形變波和地震震源之間關(guān)系的確定,為進(jìn)一步識別地震的短臨前兆特征提供參考,是地震短臨預(yù)測的重要途徑之一[3,14]。

    面波延遲時間為測點(diǎn)觀測到的面波初動與發(fā)震時刻之間的時間間隔,與面波速度和震中距相關(guān)聯(lián)[16],圖3a為面波延遲時間與震中距的關(guān)系。選取烏加河PET相對重力儀記錄到的2008—2021年在中國大陸中西部發(fā)生的20次5—8級地震的同震響應(yīng)事件作為研究對象,手動提取面波到時及面波最大幅度,計算得出面波到時與發(fā)震時刻的走時差,也就是面波延遲時間。通過擬合震中距與面波延遲時間,得出二者相關(guān)關(guān)系公式:

    圖3 面波延遲時間與震中距關(guān)系曲線及面波延遲時間與地應(yīng)力波模型關(guān)系曲線Fig. 3 Relationship between surface wave delay time and epicentral distance and relationship between surface wave delay time and in-situ stress wave model

    式中,x為震中距,y為面波延遲時間,R=0.8377,為相關(guān)系數(shù),面波延遲時間與震中距回歸計算后呈線性分布特征,面波延遲時間隨震中距增加而增大,正相關(guān)特征顯著。面波延遲時間范圍為0—20 min[14],主要由震中距決定,震中距越大面波到達(dá)時間越遲。表1為烏加河地震臺PET相對重力儀記錄到的2008—2021年在中國大陸中西部發(fā)生的20次5—8級地震的同震波形參數(shù)。

    表1 重力儀記錄到的地震波形參數(shù)Table 1 Seismic waveform parameters recorded by gravimeter

    2.2 地應(yīng)力波模型

    地應(yīng)力波主要是地殼深部介質(zhì)黏彈或流變特性的反映。地傾斜、應(yīng)變和重力等觀測記錄到的一些短周期前兆異常變化,在空間上均具有可傳播特征,表現(xiàn)為地應(yīng)力波傳播過程的動態(tài)響應(yīng)。由于地應(yīng)力波較窄的傳播途徑和快速的衰減等原因,國內(nèi)外對地應(yīng)力波傳播問題尚未深入研究[18]。

    圖3b為面波延遲時間與應(yīng)力波模型關(guān)系曲線。根據(jù)地震面波延遲時間和震中距的關(guān)系模型,模型增加lg(D/D0)項進(jìn)行分析,設(shè)置D為震中距,D0為地應(yīng)力半波長,約290 km[12,18]。計算震中距與地應(yīng)力半波長的比值后取對數(shù),得出地應(yīng)波模型后,通過擬合lg(D/D0)與面波延時,得出相關(guān)關(guān)系公式:

    式中,x為lg(D/D0),y為面波延遲時間,R=0.8556,為相關(guān)系數(shù)。面波延遲時間與lg(D/D0)回歸分析后呈線性分布特征,面波延遲時間隨震中距增加而增大,正相關(guān)特征顯著。面波延遲時間范圍為0—20 min,主要由震中距決定,震中距越大面波到達(dá)時間越遲,與式(1)計算結(jié)果差別不明顯。同震形變波在不同測點(diǎn)記錄到的最大響應(yīng)幅值以及響應(yīng)波持續(xù)時間與震級均存在一定程度的關(guān)系,面波延遲時間與地震面波的傳播速度以及震中距相關(guān)[17]。

    3 同震形變震級模型研究

    3.1 一般性關(guān)系模型

    應(yīng)用震級和變形幅度的關(guān)系模型[3,5]進(jìn)行擬合分析,獲得同震形變波最大響應(yīng)幅值和地震震級以及震中距之間的關(guān)系模型。震級增大,震時形變波幅度增大,形變幅度與震級存在正相關(guān)的關(guān)系。同時,同震形變幅度隨震中距的增加呈衰減特征,即相同震級的地震,震中距越遠(yuǎn)同震形變幅度越小。

    一般性關(guān)系震級模型公式為:

    式中,M為震級,A為振幅,D為震中距。

    由于振幅受震級和震中距共同約束,研究分析選取樣本震級和震中距偏差不能太大,為排除相關(guān)因素干擾,選取2008—2021年巴顏喀拉塊體周圍發(fā)生的6次7級以上地震進(jìn)行研究,手動提取烏加河PET重力儀記錄到的同震形變面波振幅和周期,對震級、面波振幅及震中距作為基本物理量與震級模型進(jìn)行擬合,求出振幅與震中距系數(shù)及常數(shù),得出同震形變震級公式:

    式中,M為震級,變量lgA為面波振幅的對數(shù)值,變量lgD為震中距的對數(shù)值。通過震級模擬分析,震級與面波振幅的對數(shù)值加震中距的對數(shù)值呈正相關(guān)特征,相關(guān)系數(shù)R值為0.92,震級增大,同震形變面波幅度增大。震級模擬計算得出的標(biāo)準(zhǔn)誤差為0.10,這一誤差水平與利用烏加河地震臺寬頻帶地震計測得的震級誤差基本相當(dāng),即PET連續(xù)重力儀記錄到的同震形變波形,通過提取面波振幅,可以估算得出震級大小或震中距。表2為模型震例參數(shù),圖4為一般性關(guān)系震級模型與震級擬合圖。

    圖4 一般性關(guān)系模型與震級擬合圖Fig. 4 Fitting diagram of general relationship model and magnitude

    表2 模型震例參數(shù)表Table 2 Earthquake parameters for the model

    3.2 地應(yīng)力波模型

    為檢驗同震形變觀測中場地響應(yīng)能力的大小和差異,應(yīng)用應(yīng)力波概念建立同震形變震級模型再次研究同震形變波形中各物理量間的關(guān)系,回歸分析震級模型系數(shù)及參數(shù),得出同震形變震級公式。

    震級模型公式為:

    式中,M為震級,變量A為振幅,變量D為震中距,D0約等于290 km,為地應(yīng)力波半波長。通過對2008—2021年巴顏喀拉塊體周圍發(fā)生的6次7級以上地震進(jìn)行研究,手動提取烏加河PET重力儀記錄到的同震形變震級、面波振幅、震中距作為基本物理量,計算面波振幅的對數(shù)、震中距與地應(yīng)力半波長的比值的對數(shù),回歸分析得出地應(yīng)波概念下同震形變震級模型公式為:

    通過震級模擬分析,震級與面波振幅的對數(shù)值加震中距的對數(shù)值,呈正相關(guān)特征,相關(guān)系數(shù)R值為0.92,標(biāo)準(zhǔn)誤差為0.10,與一般關(guān)系模型擬合結(jié)果一致,公式差別僅常數(shù)C不同,即PET連續(xù)重力儀記錄到的同震形變波形,通過提取面波振幅,加入地應(yīng)力波概念模型,估算得出震級大小與一般關(guān)系模型估算震級區(qū)別不大,應(yīng)用一般關(guān)系模型估算震級較方便快捷。圖5為地應(yīng)力波模型與與震級擬合圖。

    圖5 地應(yīng)力波模型與震級擬合圖Fig. 5 Fitting diagram of in-situ stress wave model and magnitude

    重力同震面波幅度衰減因地震不同,衰減的形態(tài)與速度也有較大差別,目前因地震樣本較少,無法深入研究每次地震衰減的差異性及衰減系數(shù)的差異性。

    3.3 P波S波振幅比模型

    應(yīng)用P波和S波的振幅比計算地震的震源機(jī)制解是地震學(xué)中一個常用的方法,本文通過重力儀可記錄到同震P波和S波的特性,研究了同震形變波形的震級反演,進(jìn)一步揭示同震形變波形蘊(yùn)含地球動力特征的物理機(jī)制。

    震級模型公式為:

    式中,M為震級,AS為S波振幅,AP為P波振幅,D為震中距。通過對2008—2021年巴顏喀拉塊體周圍發(fā)生的6次7級以上地震進(jìn)行研究,把震級、同震P波振幅、同震S波振幅、震中距作為基本物理量,計算同震S波振幅與同震P波振幅比值的對數(shù)、震中距與地應(yīng)力半波長的比值的對數(shù)。經(jīng)過與震級擬合回歸分析出振幅與震中距系數(shù)及常數(shù),得出同震形變震級公式為:

    通過震級模擬分析,震級與同震S波振幅及同震P波振幅比值取對數(shù)加震中距與地應(yīng)力半波長的比值的對數(shù)值,呈正相關(guān)特征,相關(guān)系數(shù)R值為0.99,標(biāo)準(zhǔn)誤差為0.04,震級相關(guān)系數(shù)與標(biāo)準(zhǔn)誤差得到最優(yōu)結(jié)果,即應(yīng)用同震P波振幅與同震S波振幅進(jìn)行震級反演與實(shí)際震級誤差最小。PET相對重力儀同震P波振幅與同震S波振幅中蘊(yùn)含了豐富的地震破裂及地震波傳播過程的信息,有益于重力學(xué)科與測震學(xué)科相互融合研究,相互佐證排除異常干擾。圖6為同震P波同震S波模型與震級擬合圖。

    圖6 P波S波幅度比模型與震級擬合圖Fig. 6 Fitting diagram of P-wave and S-wave amplitude radio model and magnitude

    4 結(jié)論與認(rèn)識

    本文應(yīng)用烏加河地震臺連續(xù)重力觀測資料,分析了重力儀與地震儀的波形相關(guān)關(guān)系,重力儀和地震儀記錄的速度值求導(dǎo)后,數(shù)據(jù)波形記錄的P波初動方向、P波與S波到時具有較高的一致性,同時在持續(xù)時間和振幅等震相特征上也高度相似。通過研究PET重力儀的同震形變物理特征,模擬分析中國大陸發(fā)生的6次7級以上地震的同震P波、S波及面波與震中距、震級的相關(guān)關(guān)系,得出一般震級模型、地應(yīng)力波模型和P波S波幅度比模型公式。將3種模型公式得出的理論震級與中國地震臺網(wǎng)中心測定的面波震級進(jìn)行對比,通過計算誤差與相關(guān)系數(shù)R值,討論了同震變化的物理機(jī)制。結(jié)論如下:

    (1)應(yīng)用震級與變形幅度間的一般性關(guān)系模型進(jìn)行擬合,得到同震形變最大響應(yīng)幅度與地震震級和震中距之間的關(guān)系模型。震級增大,震時形變波幅度增大,變形幅度與震級存在正相關(guān)的關(guān)系。同時,同震形變幅度隨震中距的增加呈衰減特征,即相同震級的地震,震中距越遠(yuǎn)同震形變幅度越小。計算PET相對重力儀數(shù)據(jù)得出震級與烏加河地震臺寬頻帶地震計記錄到的地震波形數(shù)據(jù)測定的面波震級相比,擬合所得震級與其偏差較小,兩者誤差基本相當(dāng)。

    (2)應(yīng)用地應(yīng)力波概念建立同震形變震級模型研究同震形變波形中各物理量間的關(guān)系,模擬分析震級模型系數(shù)及參數(shù),擬合計算得出的地應(yīng)力波概念同震形變震級公式與反演計算得出的一般關(guān)系模型公式結(jié)果一致,公式差別僅常數(shù)不同,應(yīng)用一般關(guān)系模型估算震級較方便快捷。

    (3)應(yīng)用同震形變P波振幅、同震S波振幅、震中距及地應(yīng)力半波長概念作為基本物理量進(jìn)行震級擬合與實(shí)際震級相比誤差最小,PET相對重力儀同震P波振幅與同震S波振幅中蘊(yùn)含了豐富的地震破裂及地震波傳播過程的信息,有益于重力學(xué)科與測震學(xué)科相互融合研究,相互佐證排除異常干擾。

    通過研究發(fā)生在巴顏喀拉塊體周緣斷裂上6次7級以上強(qiáng)震,6次地震最大面波周期其中4次為12 s,最大面波周期地震為2021年5月22日青?,敹喟l(fā)生MS7.4強(qiáng)震,面波周期為14 s。通過3種震級模型,計算理論震級與實(shí)際震級的誤差及R值,理論震級與實(shí)測震級標(biāo)準(zhǔn)誤差達(dá)到0.04,相關(guān)系數(shù)R值為0.92—0.99,進(jìn)一步驗證了模型的可靠性。重力觀測同震形變波形蘊(yùn)含地球動力特征,同震形變波形中記錄到的面波、P波和S波的振幅比可以計算地震的震級,對震級和震中距進(jìn)行一定的估計預(yù)測。通過對重力同震形變波的定量研究,可認(rèn)識到在同一力源作用下,同震波特征的差異性,這對理解形變觀測的物理本質(zhì)是有意義的。利用測震儀與重力儀震例回溯、檢驗,進(jìn)而探討地震的破裂與傳播過程,促進(jìn)地震學(xué)與重力學(xué)學(xué)科之間的相互借鑒與融合,揭示同震形變波形蘊(yùn)含地球動力特征的物理機(jī)制,挖掘能夠反映強(qiáng)地震發(fā)生的典型要素,這對預(yù)測研究有著及其重要的意義。

    猜你喜歡
    同震重力儀面波
    工程化原子重力儀綜述
    gPhone重力儀的面波頻段響應(yīng)實(shí)測研究
    地震研究(2021年1期)2021-04-13 01:04:56
    自適應(yīng)相減和Curvelet變換組合壓制面波
    基于組合滑??刂频慕^對重力儀兩級主動減振設(shè)計
    云南思茅大寨井水位地震同震響應(yīng)特征分析*
    地震研究(2018年4期)2018-11-23 02:29:36
    蘆山地震前后介質(zhì)波速變化與GPS應(yīng)變場相關(guān)性研究?
    蘆山Ms7.0地震引起的水位同震響應(yīng)特征分析
    四川地震(2014年2期)2014-12-02 04:16:30
    淺析工程勘探的面波勘探方法
    河南科技(2014年8期)2014-02-27 14:07:40
    CG-5重力儀彈簧形變對測量的影響
    川滇地區(qū)鉆孔四分量應(yīng)變儀記錄的同震應(yīng)變階分析1
    日韩欧美国产在线观看| 能在线免费观看的黄片| 久久精品久久久久久噜噜老黄 | 国产精品美女特级片免费视频播放器| 搡女人真爽免费视频火全软件| 国产精品一区二区三区四区免费观看| 亚洲欧美日韩高清在线视频| 精品久久久久久成人av| 亚洲国产色片| 狠狠狠狠99中文字幕| 国产不卡一卡二| 国内少妇人妻偷人精品xxx网站| 欧美日韩国产亚洲二区| 国产欧美日韩精品一区二区| 色尼玛亚洲综合影院| 可以在线观看的亚洲视频| 国产精品久久久久久av不卡| 久久久精品94久久精品| 亚洲一级一片aⅴ在线观看| 在线观看一区二区三区| 亚洲自拍偷在线| 国产老妇女一区| 熟妇人妻久久中文字幕3abv| 啦啦啦韩国在线观看视频| 天堂影院成人在线观看| 欧美丝袜亚洲另类| 日韩高清综合在线| 免费看美女性在线毛片视频| 亚洲自偷自拍三级| 国内揄拍国产精品人妻在线| 亚洲精品日韩在线中文字幕 | 久久久成人免费电影| 国产av在哪里看| 91午夜精品亚洲一区二区三区| 美女脱内裤让男人舔精品视频 | 我的老师免费观看完整版| 久久久欧美国产精品| 亚洲人成网站高清观看| 国产乱人偷精品视频| 女人被狂操c到高潮| 禁无遮挡网站| 搡老妇女老女人老熟妇| 爱豆传媒免费全集在线观看| 久久精品国产亚洲av香蕉五月| 国内少妇人妻偷人精品xxx网站| 在线天堂最新版资源| 色综合色国产| av在线老鸭窝| 在线免费观看不下载黄p国产| 国产成人a∨麻豆精品| 一级毛片久久久久久久久女| 亚洲精品色激情综合| 男女下面进入的视频免费午夜| 成人毛片a级毛片在线播放| 九草在线视频观看| 国产精品野战在线观看| 精品久久国产蜜桃| 全区人妻精品视频| 内射极品少妇av片p| 成人欧美大片| 亚洲一区高清亚洲精品| 老司机影院成人| 日本欧美国产在线视频| 性色avwww在线观看| 舔av片在线| 亚洲国产欧美在线一区| 久久精品夜夜夜夜夜久久蜜豆| 亚洲真实伦在线观看| 亚洲第一区二区三区不卡| 九草在线视频观看| 成人高潮视频无遮挡免费网站| 日日摸夜夜添夜夜添av毛片| 老师上课跳d突然被开到最大视频| 免费看日本二区| 国产精品久久久久久精品电影小说 | 亚洲欧洲国产日韩| 国内精品美女久久久久久| 伦理电影大哥的女人| 99热精品在线国产| 别揉我奶头 嗯啊视频| 国产视频内射| 欧美成人一区二区免费高清观看| 我要看日韩黄色一级片| 床上黄色一级片| 国产亚洲91精品色在线| 日本色播在线视频| 只有这里有精品99| 一进一出抽搐gif免费好疼| 国产精品久久久久久久电影| 成人亚洲欧美一区二区av| 伊人久久精品亚洲午夜| 久久久欧美国产精品| 99久久精品一区二区三区| 国产日韩欧美在线精品| 亚洲第一电影网av| 久久草成人影院| 乱系列少妇在线播放| 精品一区二区三区视频在线| 高清毛片免费观看视频网站| 亚洲七黄色美女视频| 一个人观看的视频www高清免费观看| 天美传媒精品一区二区| 成人特级av手机在线观看| 亚洲最大成人av| 一夜夜www| 69av精品久久久久久| 嫩草影院新地址| 一夜夜www| 亚洲aⅴ乱码一区二区在线播放| 嫩草影院精品99| 中国国产av一级| 午夜福利在线在线| 在线天堂最新版资源| 免费av观看视频| 日本成人三级电影网站| 亚洲最大成人中文| 日本熟妇午夜| 国语自产精品视频在线第100页| 夫妻性生交免费视频一级片| 国产熟女欧美一区二区| 国产极品精品免费视频能看的| 人妻夜夜爽99麻豆av| 久久热精品热| 欧美精品一区二区大全| 永久网站在线| 国产精品.久久久| 夜夜看夜夜爽夜夜摸| 亚洲电影在线观看av| 人妻夜夜爽99麻豆av| 12—13女人毛片做爰片一| 22中文网久久字幕| 成人毛片a级毛片在线播放| av在线蜜桃| 国产三级中文精品| 最好的美女福利视频网| 色噜噜av男人的天堂激情| 国产精品久久久久久精品电影| 欧美日韩一区二区视频在线观看视频在线 | 人人妻人人澡欧美一区二区| 国产av在哪里看| 中文欧美无线码| 亚洲国产精品sss在线观看| 91午夜精品亚洲一区二区三区| 人人妻人人澡欧美一区二区| 少妇人妻一区二区三区视频| 日本与韩国留学比较| 国产午夜精品久久久久久一区二区三区| 人人妻人人澡欧美一区二区| eeuss影院久久| avwww免费| 久久久久久久久久黄片| 久久久久久大精品| 日韩成人av中文字幕在线观看| 国产高清有码在线观看视频| 一边摸一边抽搐一进一小说| 欧美xxxx黑人xx丫x性爽| 男人和女人高潮做爰伦理| 国产淫片久久久久久久久| 亚洲国产色片| 晚上一个人看的免费电影| 国产日韩欧美在线精品| 国内精品久久久久精免费| 大又大粗又爽又黄少妇毛片口| 大又大粗又爽又黄少妇毛片口| 国产成人午夜福利电影在线观看| 欧美高清成人免费视频www| av在线亚洲专区| 一区二区三区免费毛片| 久久久久久久久大av| 久久精品影院6| 国产视频首页在线观看| 天美传媒精品一区二区| 丝袜美腿在线中文| 老司机影院成人| 特大巨黑吊av在线直播| 国产一级毛片在线| 亚洲在线自拍视频| 精品久久久久久久人妻蜜臀av| 内射极品少妇av片p| av福利片在线观看| 国产熟女欧美一区二区| 成人毛片60女人毛片免费| 日本五十路高清| 欧美日韩精品成人综合77777| 99热只有精品国产| 亚洲国产欧美人成| 一个人免费在线观看电影| 午夜老司机福利剧场| 亚洲成人久久爱视频| 91精品国产九色| 天堂影院成人在线观看| 波多野结衣巨乳人妻| 亚洲欧美日韩东京热| 欧美最新免费一区二区三区| 91午夜精品亚洲一区二区三区| 自拍偷自拍亚洲精品老妇| eeuss影院久久| 爱豆传媒免费全集在线观看| 一级毛片久久久久久久久女| videossex国产| 可以在线观看的亚洲视频| 国产精品久久久久久久久免| 蜜桃亚洲精品一区二区三区| 国产 一区 欧美 日韩| 少妇的逼水好多| 成人高潮视频无遮挡免费网站| 成人毛片60女人毛片免费| 自拍偷自拍亚洲精品老妇| 亚洲av男天堂| 中文字幕熟女人妻在线| 久久韩国三级中文字幕| 国产精品一区www在线观看| 午夜福利成人在线免费观看| 99热精品在线国产| 国产精品一区二区三区四区免费观看| 小说图片视频综合网站| 久久精品人妻少妇| 少妇被粗大猛烈的视频| 长腿黑丝高跟| 搡老妇女老女人老熟妇| 午夜精品在线福利| 亚洲欧美精品综合久久99| 久久人妻av系列| 久久鲁丝午夜福利片| 国产午夜福利久久久久久| 国产午夜精品久久久久久一区二区三区| 国产成人91sexporn| 综合色av麻豆| 欧美在线一区亚洲| 97热精品久久久久久| 夫妻性生交免费视频一级片| 免费观看在线日韩| 最好的美女福利视频网| 久久精品国产亚洲av香蕉五月| 国产精品三级大全| 欧美激情在线99| 久久99热6这里只有精品| 成人综合一区亚洲| 成年版毛片免费区| 国内精品美女久久久久久| 久久99精品国语久久久| 国产伦精品一区二区三区视频9| 久久久久久国产a免费观看| 少妇裸体淫交视频免费看高清| 色播亚洲综合网| 国产精品一二三区在线看| 老师上课跳d突然被开到最大视频| 精品国内亚洲2022精品成人| 亚洲欧美成人精品一区二区| 熟女电影av网| av在线观看视频网站免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久精品大字幕| 久久亚洲国产成人精品v| 精品久久久噜噜| 欧美变态另类bdsm刘玥| 久久精品91蜜桃| 五月伊人婷婷丁香| 久久久久久大精品| 性色avwww在线观看| 国产精品伦人一区二区| 少妇的逼水好多| av福利片在线观看| 男人舔女人下体高潮全视频| 久久久精品欧美日韩精品| 日韩 亚洲 欧美在线| 精品人妻一区二区三区麻豆| 黄色一级大片看看| 大型黄色视频在线免费观看| 久久精品国产清高在天天线| 偷拍熟女少妇极品色| 国产精品久久久久久精品电影小说 | 日韩成人伦理影院| 国语自产精品视频在线第100页| 亚洲精品成人久久久久久| 我要看日韩黄色一级片| 精品无人区乱码1区二区| 九九爱精品视频在线观看| 免费av观看视频| 男人舔奶头视频| 精品久久国产蜜桃| 国产片特级美女逼逼视频| 啦啦啦观看免费观看视频高清| 国产精品一区二区在线观看99 | 天堂√8在线中文| 99久久精品热视频| 在线观看免费视频日本深夜| 一本一本综合久久| 欧美日韩一区二区视频在线观看视频在线 | 亚洲精品亚洲一区二区| 免费av毛片视频| 亚洲高清免费不卡视频| 中出人妻视频一区二区| av在线亚洲专区| 国产精品久久久久久亚洲av鲁大| 亚洲av免费在线观看| 少妇猛男粗大的猛烈进出视频 | 国产人妻一区二区三区在| 国产av一区在线观看免费| 欧美变态另类bdsm刘玥| 国产免费一级a男人的天堂| 99热网站在线观看| 听说在线观看完整版免费高清| 插阴视频在线观看视频| 亚洲av电影不卡..在线观看| 久久久欧美国产精品| 国产高清不卡午夜福利| 日产精品乱码卡一卡2卡三| 黑人高潮一二区| 最后的刺客免费高清国语| 欧美高清性xxxxhd video| 国产av麻豆久久久久久久| 亚洲激情五月婷婷啪啪| 九色成人免费人妻av| 国产黄片视频在线免费观看| 亚洲一区高清亚洲精品| 亚洲无线观看免费| 国语自产精品视频在线第100页| 亚洲真实伦在线观看| 男女边吃奶边做爰视频| 成人午夜高清在线视频| 最近的中文字幕免费完整| 中文字幕久久专区| 一进一出抽搐gif免费好疼| 爱豆传媒免费全集在线观看| 美女内射精品一级片tv| 久久久久久久久久久丰满| 变态另类成人亚洲欧美熟女| 3wmmmm亚洲av在线观看| 免费观看人在逋| 国产日本99.免费观看| 一级av片app| 亚洲精华国产精华液的使用体验 | 狂野欧美激情性xxxx在线观看| 人妻少妇偷人精品九色| 成人美女网站在线观看视频| 1000部很黄的大片| 国产成年人精品一区二区| 97热精品久久久久久| 天堂影院成人在线观看| 亚洲人成网站在线播放欧美日韩| 国产精品久久久久久精品电影| 嫩草影院入口| 欧美成人精品欧美一级黄| 国产成人freesex在线| 美女大奶头视频| 69人妻影院| 亚洲国产精品国产精品| 成人特级黄色片久久久久久久| 哪里可以看免费的av片| 久久久国产成人免费| 午夜免费男女啪啪视频观看| 久久99蜜桃精品久久| 免费人成视频x8x8入口观看| 高清在线视频一区二区三区 | 久久久国产成人精品二区| 成人高潮视频无遮挡免费网站| 精品国产三级普通话版| 能在线免费看毛片的网站| 久久人妻av系列| 亚洲电影在线观看av| 久久久久久国产a免费观看| 99在线视频只有这里精品首页| 联通29元200g的流量卡| 欧美日韩一区二区视频在线观看视频在线 | 99久国产av精品| 日韩国内少妇激情av| 国产不卡一卡二| 精品久久久久久久末码| 久久人妻av系列| 欧美精品一区二区大全| 一级二级三级毛片免费看| 欧美色欧美亚洲另类二区| 日日摸夜夜添夜夜添av毛片| 国产一区二区三区av在线 | 日韩一区二区视频免费看| 神马国产精品三级电影在线观看| 精品免费久久久久久久清纯| 蜜桃亚洲精品一区二区三区| 女同久久另类99精品国产91| 亚洲熟妇中文字幕五十中出| 韩国av在线不卡| 色综合色国产| 婷婷六月久久综合丁香| 国产精品久久视频播放| 天堂网av新在线| 边亲边吃奶的免费视频| 成人漫画全彩无遮挡| 校园春色视频在线观看| 国产精品一区二区在线观看99 | 国产淫片久久久久久久久| 日韩大尺度精品在线看网址| 国产三级中文精品| 深爱激情五月婷婷| 99热这里只有精品一区| 欧美三级亚洲精品| 97超视频在线观看视频| 亚洲精品国产av成人精品| 久久这里有精品视频免费| 最好的美女福利视频网| 亚洲aⅴ乱码一区二区在线播放| 最近最新中文字幕大全电影3| 亚洲精品日韩在线中文字幕 | 日韩在线高清观看一区二区三区| 国产高潮美女av| 国国产精品蜜臀av免费| 美女高潮的动态| 国产毛片a区久久久久| 国产一区亚洲一区在线观看| 99热这里只有是精品在线观看| 婷婷精品国产亚洲av| 久久精品国产清高在天天线| 麻豆乱淫一区二区| 日韩,欧美,国产一区二区三区 | 日韩中字成人| 欧美成人一区二区免费高清观看| 国产av麻豆久久久久久久| 欧美日韩精品成人综合77777| 一本久久中文字幕| 91精品一卡2卡3卡4卡| 欧美bdsm另类| 国模一区二区三区四区视频| 男女那种视频在线观看| 国产三级中文精品| 亚洲国产精品成人综合色| 级片在线观看| 国产三级在线视频| 国产成人a∨麻豆精品| 国产综合懂色| 在线观看av片永久免费下载| 狠狠狠狠99中文字幕| 久久久久久久亚洲中文字幕| 国产亚洲5aaaaa淫片| 精品久久久久久久末码| 天堂av国产一区二区熟女人妻| 欧美日韩精品成人综合77777| 老司机影院成人| 久久精品国产亚洲av涩爱 | 99久国产av精品| 国产精品综合久久久久久久免费| 精品免费久久久久久久清纯| 久久久国产成人精品二区| 天堂av国产一区二区熟女人妻| 26uuu在线亚洲综合色| 午夜激情欧美在线| 插逼视频在线观看| 深爱激情五月婷婷| 久久精品影院6| 午夜精品国产一区二区电影 | 亚洲无线观看免费| 免费不卡的大黄色大毛片视频在线观看 | 国产精品人妻久久久久久| 日韩制服骚丝袜av| 久久99热6这里只有精品| 亚洲一区二区三区色噜噜| av视频在线观看入口| 99久久精品一区二区三区| 久久99热这里只有精品18| 韩国av在线不卡| 在线天堂最新版资源| 欧美日韩乱码在线| 亚洲在线自拍视频| 日韩欧美精品v在线| 老司机影院成人| 精品少妇黑人巨大在线播放 | 国产亚洲精品av在线| 高清日韩中文字幕在线| 乱系列少妇在线播放| 久久精品国产亚洲av天美| 亚洲精品乱码久久久久久按摩| 色播亚洲综合网| 久久亚洲国产成人精品v| 亚州av有码| 99久国产av精品| 久久欧美精品欧美久久欧美| 美女内射精品一级片tv| 欧美高清成人免费视频www| 亚洲欧美精品自产自拍| 国产综合懂色| 亚洲成人精品中文字幕电影| 美女cb高潮喷水在线观看| 国产成人freesex在线| 成年女人永久免费观看视频| 欧洲精品卡2卡3卡4卡5卡区| 国产av在哪里看| 好男人在线观看高清免费视频| 波野结衣二区三区在线| 成人特级av手机在线观看| 久久人妻av系列| 亚洲成人中文字幕在线播放| 久久精品国产亚洲网站| 亚洲欧美精品自产自拍| 免费人成视频x8x8入口观看| 国产亚洲精品av在线| 国产精品人妻久久久影院| 国产久久久一区二区三区| 2022亚洲国产成人精品| 国产亚洲av嫩草精品影院| 免费看光身美女| 日韩强制内射视频| 免费观看精品视频网站| 丰满的人妻完整版| 欧美丝袜亚洲另类| 边亲边吃奶的免费视频| 成人无遮挡网站| www.av在线官网国产| 国产精品女同一区二区软件| 欧美潮喷喷水| 国产精品一及| 亚洲不卡免费看| 精华霜和精华液先用哪个| 99热这里只有是精品在线观看| 国产精品综合久久久久久久免费| 给我免费播放毛片高清在线观看| 少妇高潮的动态图| 日本av手机在线免费观看| 久久久久久久久中文| 久久人人爽人人爽人人片va| 亚洲天堂国产精品一区在线| eeuss影院久久| 日本黄色片子视频| 男女啪啪激烈高潮av片| 亚洲国产色片| 欧美日韩一区二区视频在线观看视频在线 | 国产乱人偷精品视频| 精品人妻视频免费看| 少妇人妻精品综合一区二区 | 亚洲av男天堂| 亚洲美女搞黄在线观看| www日本黄色视频网| 韩国av在线不卡| 国产精品综合久久久久久久免费| 99久久精品国产国产毛片| 免费无遮挡裸体视频| 青青草视频在线视频观看| 波多野结衣巨乳人妻| 变态另类丝袜制服| 亚洲欧美日韩高清在线视频| 99热网站在线观看| 日本在线视频免费播放| 亚洲成人中文字幕在线播放| 久久人人爽人人片av| 精品人妻偷拍中文字幕| 99久国产av精品国产电影| 女同久久另类99精品国产91| 国产成人a∨麻豆精品| 国产成人一区二区在线| 国产精品爽爽va在线观看网站| 国产色婷婷99| 久久中文看片网| 国产精品永久免费网站| avwww免费| 精品日产1卡2卡| 成人三级黄色视频| 99热全是精品| 亚洲av中文av极速乱| 亚洲成人久久爱视频| 色哟哟哟哟哟哟| 精品欧美国产一区二区三| 国产精品99久久久久久久久| 亚洲精品国产成人久久av| 少妇的逼水好多| 乱人视频在线观看| 国产精品不卡视频一区二区| 哪里可以看免费的av片| 大又大粗又爽又黄少妇毛片口| 一边亲一边摸免费视频| 一级毛片aaaaaa免费看小| 精品午夜福利在线看| 真实男女啪啪啪动态图| 亚洲av免费在线观看| 真实男女啪啪啪动态图| 亚洲18禁久久av| 99热精品在线国产| 色综合站精品国产| 午夜免费男女啪啪视频观看| 久久久精品94久久精品| av专区在线播放| 色哟哟·www| 午夜免费男女啪啪视频观看| 丰满乱子伦码专区| 在线国产一区二区在线| 九色成人免费人妻av| 成人午夜精彩视频在线观看| 国产成人福利小说| 成人欧美大片| АⅤ资源中文在线天堂| 99久久精品一区二区三区| 久久久久久久午夜电影| 欧美高清成人免费视频www| 99九九线精品视频在线观看视频| 欧美日韩国产亚洲二区| а√天堂www在线а√下载| 亚洲欧美日韩东京热| 久久6这里有精品| 婷婷亚洲欧美| 国模一区二区三区四区视频| 亚洲电影在线观看av| 亚洲欧美日韩高清专用| 中文字幕av成人在线电影| 韩国av在线不卡| 亚洲成人中文字幕在线播放| 欧美另类亚洲清纯唯美| 99热网站在线观看| 天堂网av新在线| 狂野欧美激情性xxxx在线观看| 亚洲性久久影院| 欧美日韩乱码在线| 天天一区二区日本电影三级| 色播亚洲综合网| 看免费成人av毛片| 国内精品宾馆在线| 日产精品乱码卡一卡2卡三| 啦啦啦啦在线视频资源| 18禁黄网站禁片免费观看直播|