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

    九寨溝MS7.0地震強(qiáng)地震動(dòng)模擬及參數(shù)敏感性分析

    2022-10-10 02:01:24李明會(huì)景月嶺張宇璇鄭司元宛明晨何李浩
    關(guān)鍵詞:九寨溝震動(dòng)臺(tái)站

    李明會(huì) 景月嶺,2 張宇璇 鄭司元 宛明晨 何李浩

    1 合肥工業(yè)大學(xué)土木與水利工程學(xué)院,合肥市屯溪路193號(hào),230001 2 水沙科學(xué)與水利水電工程國家重點(diǎn)實(shí)驗(yàn)室,北京市清華園1號(hào),100089

    2017-08-08四川省九寨溝縣發(fā)生MS7.0地震。關(guān)于九寨溝地震強(qiáng)地面運(yùn)動(dòng)數(shù)值模擬的研究取得了一些成果[1-4],但仍有不足之處,如應(yīng)用復(fù)合震源模型進(jìn)行模擬及其震源參數(shù)敏感性分析的研究較少、隨機(jī)有限斷層法計(jì)算地震動(dòng)分量方向單一、地震峰值加速度擬合精度不高等。為此,本文應(yīng)用復(fù)合震源模型對(duì)九寨溝MS7.0地震強(qiáng)地面運(yùn)動(dòng)數(shù)值模擬展開研究,并探討不同模型參數(shù)對(duì)合成地震動(dòng)的影響規(guī)律,為九寨溝地區(qū)工程場(chǎng)地地震安全評(píng)價(jià)提供參考依據(jù)。

    1 復(fù)合震源模型基本理論

    復(fù)合震源模型是用于計(jì)算寬頻帶強(qiáng)地面運(yùn)動(dòng)的有限源模型[5]。該模型通過描述震源的運(yùn)動(dòng)學(xué)特征,結(jié)合地震波在層狀介質(zhì)中的傳播過程,將廣義反射和透射系數(shù)法[6]得到的格林函數(shù)與震源時(shí)間函數(shù)進(jìn)行卷積,計(jì)算得到寬頻帶地震動(dòng)。模型中,強(qiáng)震震源被認(rèn)為是大量具有恒定應(yīng)力降圓形子震的疊加。首先將斷層面劃分成一系列尺寸相同的正方形子斷層網(wǎng)格,之后用一定數(shù)量和相應(yīng)半徑大小分布服從Frankel自相似模型[7]的圓形子震填充該斷層面,并允許其部分重疊。圓形子震相當(dāng)于一系列小震,在震源破裂過程中相繼被激發(fā),最終釋放的能量在斷層面上不斷被疊加,達(dá)到模擬強(qiáng)震的目的。設(shè)圓形子震半徑為R、數(shù)目為N,則有:

    (1)

    式中,分形維數(shù)D=2,控制斷層面上圓形子震數(shù)量的分布及圓形子震的上升時(shí)間,P為比例常數(shù):

    (2)

    式中,M0為地震矩,Δσ為應(yīng)力降,Rmax為最大圓形子震半徑,Rmin為最小圓形子震半徑。對(duì)式(1)進(jìn)行積分,可得半徑大于R且小于Rmax圓形子震的分布數(shù)量N(R):

    (3)

    (4)

    (5)

    (6)

    式中,β為剪切波的速度。斷層面上滑動(dòng)位移量的空間分布由圓形子震滑動(dòng)位移量疊加得到,第i個(gè)圓形子震的滑動(dòng)位移量ΔUi為[1]:

    (7)

    式中,C為依賴斷層類型的常數(shù),μ為介質(zhì)剪切模量。

    2 地震數(shù)據(jù)與震源參數(shù)

    本文研究區(qū)域(103°~105°E,31°~35°N)位于巴顏喀拉塊體,該塊體近年來地震頻發(fā)。九寨溝地震后,中國地震臺(tái)網(wǎng)中心等機(jī)構(gòu)公布了其震源位置、地震矩、震源深度、斷層節(jié)面的走向、傾角和滑動(dòng)角等震源參數(shù)(表1),其中,根據(jù)美國地質(zhì)調(diào)查局的研究,斷層節(jié)面Ⅱ?yàn)閷?shí)際斷層面[1],因此本文取節(jié)面Ⅱ進(jìn)行研究。

    表1 九寨溝地震震源參數(shù)

    基于孟令媛等[9]對(duì)余震分布情況及斷層面上滑動(dòng)位移的反演結(jié)果,建立長30 km、寬20 km的斷層模型,并將其劃分為150個(gè)2 km×2 km的正方形子斷層。斷層面上圓形子震的分布如圖1所示,其中圓形子震地震矩的總和等于總地震矩。通過式(7)計(jì)算得到斷層面上滑動(dòng)位移量的不均勻分布,如圖2所示。九寨溝地震震中距150 km區(qū)域范圍內(nèi)臺(tái)站信息如表2所示,九寨溝地區(qū)介質(zhì)速度信息(圖3)來自全球地殼模型Crust1.0(https:∥igppweb.ucsd.edu/~gabi/crust1.html),臺(tái)站位置分布如圖4所示。對(duì)所有臺(tái)站觀測(cè)數(shù)據(jù)進(jìn)行基線校準(zhǔn)和濾波處理,帶通濾波范圍為0.1~20 Hz。

    圖1 斷層面上圓形子震分布Fig.1 Distribution of circular subevents on the fault plane

    圖2 斷層面上滑移量的不均勻分布Fig.2 Non-uniform distribution of slip on the fault plane

    表2 臺(tái)站信息

    圖3 九寨溝地區(qū)速度模型Fig.3 Velocity model of the Jiuzhaigou region

    圖4 臺(tái)站分布Fig.4 Distribution of the stations

    3 參數(shù)敏感性分析

    在運(yùn)動(dòng)學(xué)震源模型中,斷層的破裂面積、長度、寬度、埋置深度、走向、傾角和滑移角等是描述斷層特征的參數(shù);地震矩、應(yīng)力降、平均滑動(dòng)時(shí)間、平均破裂速度、破裂方式和初始破裂點(diǎn)位置等是描述平均破裂過程的參數(shù)。這些震源參數(shù)控制著地震動(dòng)強(qiáng)度和分布[10],分析這些參數(shù)的影響規(guī)律可為地震動(dòng)合成精度提供保證。本文首先以無強(qiáng)震觀測(cè)記錄為前提,研究地震矩、應(yīng)力降、破裂速度、子斷層尺寸、震源深度、斷層走向、傾角和滑移角等8個(gè)震源參數(shù)在所有合理取值范圍內(nèi)對(duì)復(fù)合震源模型計(jì)算結(jié)果的影響程度;然后根據(jù)九寨溝地震觀測(cè)數(shù)據(jù)確定適合的震源參數(shù)研究范圍,并進(jìn)行敏感性規(guī)律分析,在其他參數(shù)保持不變的條件下,僅改變單個(gè)參數(shù)變量,利用評(píng)價(jià)指標(biāo)比較分析計(jì)算結(jié)果,總結(jié)震源參數(shù)對(duì)地震動(dòng)計(jì)算的影響規(guī)律。

    3.1 評(píng)價(jià)參數(shù)和評(píng)價(jià)指標(biāo)

    為表征震源參數(shù)對(duì)臺(tái)站NS和EW向峰值加速度(peak ground acceleration,PGA)的影響,用誤差值E表示計(jì)算的PGA相對(duì)于觀測(cè)的PGA的變化:

    E=|AOBS-ASYN|/AOBS×100%

    (8)

    式中,AOBS為觀測(cè)的PGA,ASYN為計(jì)算的PGA。

    地震動(dòng)強(qiáng)度主要用PGA、峰值速度或譜加速度等地震動(dòng)參數(shù)表達(dá)。本文以PGA為例,根據(jù)王德才等[11]的參數(shù)敏感性分析模型,通過計(jì)算模型偏差和模型標(biāo)準(zhǔn)差討論不同震源參數(shù)對(duì)PGA的影響。

    用模型偏差B表征計(jì)算值和觀測(cè)值之間的差異,偏差越大,計(jì)算值越偏離觀測(cè)數(shù)據(jù):

    (9)

    式中,M為總臺(tái)站數(shù)。

    用模型標(biāo)準(zhǔn)差S描述數(shù)據(jù)離散程度,標(biāo)準(zhǔn)差越大,數(shù)據(jù)波動(dòng)性越強(qiáng):

    (10)

    3.2 參數(shù)敏感性分析結(jié)果

    首先選取地震矩、應(yīng)力降、破裂速度、子斷層尺寸、震源深度、斷層走向、傾角和滑移角等8個(gè)震源參數(shù)的最大合理取值范圍,計(jì)算得到不同震源參數(shù)數(shù)值下臺(tái)站NS和EW向的PGA,然后以模型偏差B和標(biāo)準(zhǔn)差S為衡量標(biāo)準(zhǔn),分析8個(gè)震源參數(shù)對(duì)復(fù)合震源模型計(jì)算結(jié)果的影響程度。

    地震矩根據(jù)Hanks等[12]提出的地震矩與矩震級(jí)的關(guān)系,在強(qiáng)震6~8級(jí)范圍內(nèi)取值計(jì)算,間隔為1×1018Nm;應(yīng)力降參考王振宇等[13]的取值范圍10~200 bar,間隔為5 bar;斷層的平均破裂速度約為剪切波速的0.6~0.9倍[14],本文選取此區(qū)間,并設(shè)定區(qū)間間隔為0.05;子斷層尺寸取值分別為1 km×1 km、2 km×2 km、6 km×6 km、8 km×8 km和10 km×10 km;震源深度范圍選定為0~30 km,取值間隔為1 km;斷層走向、傾角和滑移角依據(jù)理論范圍分別取0°~360°、0°~90°和-180°~180°,取值間隔分別為10°、5°和10°。

    8個(gè)震源參數(shù)取值范圍和梯度各不相同,本文通過式(9)和式(10)計(jì)算出相應(yīng)PGA的偏差和標(biāo)準(zhǔn)差,然后將每個(gè)震源參數(shù)下PGA偏差、標(biāo)準(zhǔn)差的最大值及最小值作差,消除AOBS實(shí)測(cè)數(shù)據(jù)的影響,得到每個(gè)震源參數(shù)下PGA偏差和標(biāo)準(zhǔn)差的變化幅度,結(jié)果如圖5所示。由圖5可知,地震矩、應(yīng)力降、破裂速度、子斷層尺寸和震源深度對(duì)PGA偏差和標(biāo)準(zhǔn)差的影響明顯,其中地震矩和應(yīng)力降影響更突出;斷層走向、傾角和滑移角對(duì)PGA偏差和標(biāo)準(zhǔn)差的影響極不明顯。根據(jù)地震學(xué)理論,地震矩直接決定地震釋放的能量大小,地震矩越大,地震波攜帶的能量越多,使得整體加速度計(jì)算值增加;應(yīng)力降表示斷層上某點(diǎn)破裂前后的應(yīng)力強(qiáng)度變化,應(yīng)力降越大,斷層上釋放的應(yīng)力越多;破裂速度決定斷層面圓形子震的發(fā)震時(shí)間,間接影響整個(gè)斷層面的破裂過程;子斷層尺寸的劃分會(huì)對(duì)斷層面上滑移量分布產(chǎn)生很大的影響;震級(jí)相同的地震,震源越淺,造成的破壞越嚴(yán)重。

    圖5 不同震源參數(shù)對(duì)PGA偏差和標(biāo)準(zhǔn)差的影響幅度Fig.5 The influence amplitude of different sourceparameters on PGA deviation and standard deviation

    基于8個(gè)震源參數(shù)的影響規(guī)律,結(jié)合表1不同研究機(jī)構(gòu)和眾多學(xué)者[2-4,8]提供的震源參數(shù),確定九寨溝地震震源參數(shù)研究范圍,進(jìn)一步研究8個(gè)震源參數(shù)對(duì)九寨溝11個(gè)臺(tái)站NS和EW向PGA的影響大小和規(guī)律,取值間隔如圖6所示。

    從圖6看出,地震矩、應(yīng)力降、破裂速度、子斷層尺寸和震源深度對(duì)臺(tái)站PGA的影響大;斷層走向、傾角和滑移角對(duì)臺(tái)站PGA的影響小。說明在應(yīng)用復(fù)合震源模型計(jì)算九寨溝地震時(shí),即使8個(gè)震源參數(shù)取值范圍有所限制,沒有取到理論范圍內(nèi)的所有值,但其對(duì)計(jì)算結(jié)果的影響規(guī)律基本與圖5一致。

    圖6 不同震源參數(shù)對(duì)臺(tái)站PGA的影響Fig.6 Influence of different source parameters on station PGA

    圖7為地震矩、應(yīng)力降、破裂速度、子斷層尺寸和震源深度等5個(gè)震源參數(shù)對(duì)九寨溝地震地震動(dòng)模擬的影響規(guī)律,其中實(shí)心點(diǎn)表示本文取值。由表3和圖7可知,受地震矩影響,PGA的偏差變化范圍超過1,PGA的偏差和標(biāo)準(zhǔn)差隨地震矩的增加而減小,其中PGA的偏差下降趨勢(shì)近似為線性變化,地震矩為7×1018Nm時(shí),PGA的偏差基本為0,同時(shí)PGA的標(biāo)準(zhǔn)差也較小。應(yīng)力降越大,臺(tái)站整體PGA越大,PGA的偏差越小,應(yīng)力降為40 bar時(shí),PGA的偏差接近于0,這與俞瑞芳等[15]應(yīng)用隨機(jī)有限斷層法得到PGA隨著應(yīng)力降增大而逐漸增加的結(jié)論類似。破裂速度越大,PGA的偏差越小,當(dāng)破裂速度等于剪切速度0.8倍(2.96 km/s)時(shí),偏差接近于0,此時(shí)標(biāo)準(zhǔn)差也最小,與傅磊等[16]應(yīng)用隨機(jī)有限斷層法和李宗超等[4]應(yīng)用經(jīng)驗(yàn)格林函數(shù)法在地震模擬過程中選取的破裂速度與剪切速度的比例關(guān)系一致,PGA的偏差變化范圍為0.9。子斷層尺寸越小,PGA的偏差越小,最佳尺寸為2 km×2 km,此時(shí)標(biāo)準(zhǔn)差也為最小值,與Dang等[1]利用隨機(jī)有限斷層法的研究結(jié)論相同。震源深度在20 km處PGA的偏差幾乎為0,同時(shí)標(biāo)準(zhǔn)差最小,隨著震源深度的增加,PGA的偏差增長趨勢(shì)變緩;而且如圖6所示,震源深度為25 km時(shí),大部分臺(tái)站的峰值加速度很小,與劉啟方等[10]運(yùn)用震源位錯(cuò)模型得到的研究結(jié)果相似。

    表3 不同震源參數(shù)PGA偏差和標(biāo)準(zhǔn)差統(tǒng)計(jì)結(jié)果

    圖7 不同震源參數(shù)對(duì)PGA偏差和標(biāo)準(zhǔn)差的影響Fig.7 Influence of different source parameters on PGA deviation and standard deviation

    斷層走向、傾角和滑移角中,PGA的偏差最大變化幅度為0.08,標(biāo)準(zhǔn)差最大變化幅度僅為0.1,對(duì)計(jì)算結(jié)果影響不明顯,與王德才等[11]采用混合方法[17]和隨機(jī)方法[18]研究的結(jié)果一致。

    3.3 地震輸入?yún)?shù)組合

    由上文分析結(jié)果可知,震源參數(shù)整體對(duì)PGA偏差影響幅度大,對(duì)標(biāo)準(zhǔn)差影響幅度相對(duì)較小,因此以PGA的偏差為主要衡量標(biāo)準(zhǔn),標(biāo)準(zhǔn)差為次要標(biāo)準(zhǔn),得到一組與觀測(cè)記錄擬合效果較好的震源參數(shù)。從表1各機(jī)構(gòu)提供的數(shù)據(jù)中選取地震矩、破裂速度、子斷層尺寸、震源深度、斷層走向、傾角和滑移角;應(yīng)力降參考Sun等[3]的結(jié)果,品質(zhì)因子Q(f)和高頻衰減kappa值參考孟令媛等[9]的結(jié)果,具體數(shù)值如表4所示。

    表4 九寨溝地震復(fù)合震源模型輸入的震源物理參數(shù)

    4 地震動(dòng)模擬結(jié)果分析

    根據(jù)式(8)計(jì)算11個(gè)臺(tái)站NS和EW向PGA的相對(duì)誤差,并分別與Dang等[1]和李宗超等[4]的計(jì)算結(jié)果進(jìn)行對(duì)比(圖8)??梢钥闯?,11個(gè)臺(tái)站NS和EW向PGA相對(duì)誤差均小于50%。62SHW站NS和EW向PGA相對(duì)誤差最小,均小于3%;51HSL站NS向、51HSD站NS和EW向PGA誤差較大,且計(jì)算值均大于觀測(cè)值,主要原因可能是臺(tái)站震中距較大,導(dǎo)致傳播過程中散射的復(fù)雜性增加,同時(shí),臺(tái)站局部山體地形引起的場(chǎng)地效應(yīng)也會(huì)增加計(jì)算誤差。11個(gè)臺(tái)站NS和EW向PGA相對(duì)誤差均在20%以內(nèi),說明所有臺(tái)站PGA與實(shí)際數(shù)據(jù)相比誤差小。與Dang等[1]和李宗超等[4]的計(jì)算結(jié)果相比,除李宗超等[4]51JZY站EW向、62ZHQ站NS向和62MXT站NS向外,本文各臺(tái)站PGA相對(duì)誤差均較優(yōu)。為方便直觀地進(jìn)行地震動(dòng)持時(shí)和波形對(duì)比,所有時(shí)程均截取45 s時(shí)間窗(圖9),臺(tái)站加速度記錄數(shù)據(jù)中最小PGA為2.75 cm/s2,為避免地震動(dòng)持時(shí)出現(xiàn)零持時(shí)這種不合理狀態(tài),選取Bolt持時(shí)的閾值為0.5倍的PGA 作為計(jì)算相對(duì)Bolt持時(shí)的閾值[19]。從圖9看出,個(gè)別臺(tái)站誤差較大,但大多數(shù)臺(tái)站地震動(dòng)持時(shí)基本保持一致。如圖10所示,合成反應(yīng)譜整體變化趨勢(shì)與觀測(cè)反應(yīng)譜一致,說明11個(gè)臺(tái)站的峰值加速度、地震動(dòng)持時(shí)和反應(yīng)譜整體擬合效果較好,與觀測(cè)數(shù)據(jù)較為符合。

    圖8 不同臺(tái)站NS和EW向PGA的相對(duì)誤差對(duì)比Fig.8 Comparison of relative errors of PGA in NS and EW directions at different stations

    圖9 不同臺(tái)站NS和EW向加速度時(shí)程Fig.9 Time history of acceleration in NS and EW directions at different stations

    圖10 不同臺(tái)站合成和觀測(cè)加速度反應(yīng)譜Fig.10 Synthetic and observed acceleration response spectra at different stations

    部分臺(tái)站地震動(dòng)三要素計(jì)算結(jié)果精度偏低的可能原因是:1)沒有充分考慮地形因素對(duì)地震動(dòng)的影響,部分臺(tái)站如51HSL、51HSS和51HSD等位于山體上,山體地形的高度及底部延展距離對(duì)斷層動(dòng)力學(xué)破裂過程影響較大,進(jìn)而影響到相應(yīng)的地面地震動(dòng)分布[20],山體的坡角、地震動(dòng)入射角度會(huì)對(duì)山體地震動(dòng)產(chǎn)生抑制作用,導(dǎo)致臺(tái)站PGA變小[21];2)本文在計(jì)算地震動(dòng)時(shí)采用各向同性散射模型,但實(shí)際地震波在地殼介質(zhì)中的傳播是各向異性的,而且隨著傳播路徑的增加,散射的影響也會(huì)逐漸加強(qiáng);3)有些震源參數(shù)的選取進(jìn)行了簡化,在地震動(dòng)模擬過程中假定某些參數(shù)(如應(yīng)力降和破裂速度)為常數(shù),忽略了參數(shù)與時(shí)間和空間的相關(guān)性;4)多數(shù)臺(tái)站低頻段反應(yīng)譜模擬結(jié)果高于觀測(cè)值,可能是由于本文模型在破裂過程中沒有考慮到圓形子震傾角的可變性,或是在合成格林函數(shù)中未納入小尺度地殼非均勻性的散射效應(yīng)[22]。

    PGA是表征地震能量的重要指標(biāo)之一,能夠?yàn)榻ㄖ锟拐鹪O(shè)防提供合理的參考依據(jù)。綜合現(xiàn)有11個(gè)臺(tái)站PGA的觀測(cè)值和計(jì)算值可以得到九寨溝地震PGA 空間分布特征,如圖11所示。由圖11可知,隨著震中距增大,NS和EW向PGA從震源中心呈橢圓形向四周衰減擴(kuò)散,長軸方向大致呈SW走向,這可能與斷層破裂的方向性有關(guān),且2個(gè)方向PGA最大值均出現(xiàn)在51JZB站;在51MXD和51PWM站之間有一個(gè)PGA值未連通的區(qū)域,可能是由于該區(qū)域缺少合適的臺(tái)站記錄所致。總體而言,PGA計(jì)算值和觀測(cè)值的空間分布特征較為一致。

    圖11 九寨溝地震PGA空間分布特征Fig.11 PGA spatial distribution characteristics of Jiuzhaigou earthquake

    5 結(jié) 語

    1)在最大合理變化范圍內(nèi),相對(duì)于斷層走向、傾角和滑移角,地震矩、應(yīng)力降、破裂速度、子斷層尺寸和震源深度對(duì)地震動(dòng)計(jì)算結(jié)果影響較大。

    2)震源參數(shù)敏感性分析規(guī)律為九寨溝地震地震動(dòng)模擬的參數(shù)選取提供了依據(jù),進(jìn)而得到一組合理的震源參數(shù)組合,同時(shí)能夠獲得九寨溝地區(qū)內(nèi)任意一點(diǎn)PGA的變化區(qū)間。這一結(jié)論可以進(jìn)一步推廣到其他缺乏強(qiáng)震觀測(cè)記錄的場(chǎng)地,計(jì)算出該區(qū)域強(qiáng)度特征的上限值,為建筑物抗震設(shè)防提供依據(jù)。

    3)復(fù)合震源模型強(qiáng)地震動(dòng)模擬了特定圓形子震在斷層上的破裂過程、斷層面上滑動(dòng)位移的不均勻空間分布以及震源時(shí)間函數(shù),計(jì)算結(jié)果可以再現(xiàn)九寨溝地震的地震動(dòng)特征,模擬區(qū)域臺(tái)站的地震動(dòng)PGA、地震動(dòng)持時(shí)、反應(yīng)譜和PGA分布特征與觀測(cè)值擬合較好,反應(yīng)譜在高頻段擬合精度更高。

    猜你喜歡
    九寨溝震動(dòng)臺(tái)站
    中國科學(xué)院野外臺(tái)站檔案工作回顧
    氣象基層臺(tái)站建設(shè)
    西藏科技(2021年12期)2022-01-17 08:46:38
    震動(dòng)減脂儀可以減肥?
    己亥秋日九寨溝采風(fēng)得句
    岷峨詩稿(2019年4期)2019-11-14 21:25:32
    赴九寨溝道上(外四首)
    岷峨詩稿(2019年4期)2019-04-20 09:02:06
    題九寨溝(外五首)
    岷峨詩稿(2019年4期)2019-04-20 09:01:58
    震后九寨溝縱覽(外四首)
    岷峨詩稿(2019年4期)2019-04-20 09:01:58
    水電工程場(chǎng)地地震動(dòng)確定方法
    振動(dòng)攪拌 震動(dòng)創(chuàng)新
    中國公路(2017年18期)2018-01-23 03:00:38
    基層臺(tái)站綜合觀測(cè)業(yè)務(wù)管理之我見
    西藏科技(2015年6期)2015-09-26 12:12:13
    成人永久免费在线观看视频| 国产亚洲精品第一综合不卡| 色综合站精品国产| 免费看十八禁软件| 天堂影院成人在线观看| www日本在线高清视频| 国产一区二区三区视频了| av中文乱码字幕在线| 久久久久久久久久黄片| 久久久久国内视频| 欧美国产日韩亚洲一区| 黑人巨大精品欧美一区二区mp4| 欧美日本亚洲视频在线播放| 制服丝袜大香蕉在线| bbb黄色大片| 50天的宝宝边吃奶边哭怎么回事| 亚洲 欧美 日韩 在线 免费| 国产蜜桃级精品一区二区三区| 日本a在线网址| 97超级碰碰碰精品色视频在线观看| 国产99久久九九免费精品| 夜夜看夜夜爽夜夜摸| 国产av一区二区精品久久| 狠狠狠狠99中文字幕| 99热6这里只有精品| 国产精品1区2区在线观看.| 精品人妻1区二区| 19禁男女啪啪无遮挡网站| 国产精华一区二区三区| 九色国产91popny在线| av免费在线观看网站| 精品国产一区二区三区四区第35| 色播亚洲综合网| 久久精品国产99精品国产亚洲性色| 人人妻人人看人人澡| 国产高清视频在线播放一区| 久久久国产精品麻豆| 麻豆成人午夜福利视频| 色播在线永久视频| 最好的美女福利视频网| 男女之事视频高清在线观看| 亚洲av电影不卡..在线观看| 又黄又爽又免费观看的视频| 黄色女人牲交| 给我免费播放毛片高清在线观看| 国产v大片淫在线免费观看| 一本综合久久免费| 免费看a级黄色片| 露出奶头的视频| 精品无人区乱码1区二区| 免费一级毛片在线播放高清视频| av天堂在线播放| 久久久国产精品麻豆| bbb黄色大片| 波多野结衣巨乳人妻| 国内少妇人妻偷人精品xxx网站 | 妹子高潮喷水视频| 老汉色av国产亚洲站长工具| 一级作爱视频免费观看| 看黄色毛片网站| 伊人久久大香线蕉亚洲五| 国产一区二区在线av高清观看| 久9热在线精品视频| 精品国产乱子伦一区二区三区| 国产精品亚洲一级av第二区| 长腿黑丝高跟| 久9热在线精品视频| 一级毛片精品| 女人高潮潮喷娇喘18禁视频| 免费女性裸体啪啪无遮挡网站| 十分钟在线观看高清视频www| 国产免费男女视频| 精华霜和精华液先用哪个| 性色av乱码一区二区三区2| 成年女人毛片免费观看观看9| 国产伦在线观看视频一区| 日韩精品青青久久久久久| 国产一区二区激情短视频| 久久精品国产亚洲av高清一级| 亚洲人成电影免费在线| 丰满人妻熟妇乱又伦精品不卡| 怎么达到女性高潮| 日本免费一区二区三区高清不卡| 久久中文字幕一级| 午夜成年电影在线免费观看| 亚洲中文字幕日韩| 国产精品九九99| 18禁黄网站禁片午夜丰满| 亚洲国产欧美网| 日本五十路高清| 丝袜美腿诱惑在线| 精品久久久久久成人av| 国产亚洲精品久久久久久毛片| 宅男免费午夜| 美女 人体艺术 gogo| 看黄色毛片网站| 午夜影院日韩av| 露出奶头的视频| 亚洲成国产人片在线观看| 色婷婷久久久亚洲欧美| 久热这里只有精品99| 视频在线观看一区二区三区| 久久 成人 亚洲| 国产精品久久久久久亚洲av鲁大| 欧美黑人巨大hd| 中文在线观看免费www的网站 | 亚洲成av片中文字幕在线观看| 中文字幕人妻熟女乱码| 在线国产一区二区在线| 午夜影院日韩av| 欧美zozozo另类| √禁漫天堂资源中文www| 欧美成人午夜精品| 精品久久久久久久毛片微露脸| 亚洲成a人片在线一区二区| 女性被躁到高潮视频| 草草在线视频免费看| av超薄肉色丝袜交足视频| 亚洲天堂国产精品一区在线| 亚洲 国产 在线| 丝袜美腿诱惑在线| 最新美女视频免费是黄的| 亚洲精品中文字幕在线视频| 叶爱在线成人免费视频播放| 最新美女视频免费是黄的| 一进一出抽搐gif免费好疼| 色综合站精品国产| 精品日产1卡2卡| 国产精品一区二区精品视频观看| 亚洲中文av在线| 桃色一区二区三区在线观看| 亚洲熟女毛片儿| 99在线人妻在线中文字幕| 国内精品久久久久久久电影| 国产乱人伦免费视频| 91麻豆精品激情在线观看国产| 亚洲色图 男人天堂 中文字幕| 国产激情偷乱视频一区二区| 啦啦啦免费观看视频1| 亚洲精品久久国产高清桃花| 久久青草综合色| 精品人妻1区二区| 欧美在线一区亚洲| 亚洲av中文字字幕乱码综合 | 午夜久久久在线观看| 久久香蕉精品热| 成人免费观看视频高清| 一级作爱视频免费观看| 亚洲成国产人片在线观看| 欧美性猛交╳xxx乱大交人| 欧美乱妇无乱码| 欧美另类亚洲清纯唯美| 人人妻人人看人人澡| 久久欧美精品欧美久久欧美| 国产私拍福利视频在线观看| 国产色视频综合| 高潮久久久久久久久久久不卡| 国产爱豆传媒在线观看 | 久久精品人妻少妇| 国产黄a三级三级三级人| 最近最新中文字幕大全免费视频| 首页视频小说图片口味搜索| 99久久99久久久精品蜜桃| 国产av在哪里看| 香蕉丝袜av| 午夜精品久久久久久毛片777| 久久热在线av| 熟女电影av网| 午夜免费观看网址| 男女之事视频高清在线观看| 亚洲五月婷婷丁香| 丰满人妻熟妇乱又伦精品不卡| 国产精品乱码一区二三区的特点| 亚洲真实伦在线观看| 视频区欧美日本亚洲| 免费看日本二区| 国产乱人伦免费视频| 欧美激情 高清一区二区三区| 俺也久久电影网| 怎么达到女性高潮| 波多野结衣高清无吗| 天天添夜夜摸| 亚洲成国产人片在线观看| 婷婷亚洲欧美| 色综合亚洲欧美另类图片| 中亚洲国语对白在线视频| 久久久久精品国产欧美久久久| 精品不卡国产一区二区三区| 亚洲激情在线av| 日本 欧美在线| 老司机午夜福利在线观看视频| 国产又黄又爽又无遮挡在线| 精品久久久久久,| www.www免费av| 久久久久久久久中文| 一a级毛片在线观看| 真人一进一出gif抽搐免费| 国产伦一二天堂av在线观看| 国产欧美日韩一区二区三| 非洲黑人性xxxx精品又粗又长| 色老头精品视频在线观看| 亚洲熟妇熟女久久| 一卡2卡三卡四卡精品乱码亚洲| 每晚都被弄得嗷嗷叫到高潮| 亚洲精华国产精华精| 亚洲av中文字字幕乱码综合 | 久久久久久久久免费视频了| 一a级毛片在线观看| 欧美日韩一级在线毛片| 亚洲三区欧美一区| 亚洲一码二码三码区别大吗| 热re99久久国产66热| 99re在线观看精品视频| 亚洲成a人片在线一区二区| 国产视频一区二区在线看| 91成年电影在线观看| 精品久久久久久久久久久久久 | 麻豆成人午夜福利视频| www日本在线高清视频| 久久久水蜜桃国产精品网| 欧美另类亚洲清纯唯美| 无遮挡黄片免费观看| 男人的好看免费观看在线视频 | 精品久久久久久,| 久久午夜亚洲精品久久| 侵犯人妻中文字幕一二三四区| 欧美性猛交黑人性爽| 国内久久婷婷六月综合欲色啪| 韩国精品一区二区三区| 后天国语完整版免费观看| 岛国在线观看网站| 免费看美女性在线毛片视频| 国产激情欧美一区二区| 免费电影在线观看免费观看| www.精华液| 人人妻,人人澡人人爽秒播| 99riav亚洲国产免费| 9191精品国产免费久久| 亚洲精品美女久久av网站| 午夜久久久在线观看| 日韩精品青青久久久久久| www.www免费av| 一本大道久久a久久精品| 日本精品一区二区三区蜜桃| 欧美国产日韩亚洲一区| 大型黄色视频在线免费观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲片人在线观看| 亚洲国产看品久久| 欧美成人性av电影在线观看| 老司机午夜福利在线观看视频| 午夜免费鲁丝| 精品一区二区三区视频在线观看免费| 精品乱码久久久久久99久播| 亚洲五月婷婷丁香| 怎么达到女性高潮| 精品电影一区二区在线| 亚洲人成77777在线视频| 在线看三级毛片| 国产免费男女视频| 亚洲片人在线观看| 这个男人来自地球电影免费观看| 日本一本二区三区精品| 日韩欧美 国产精品| 精品久久久久久久毛片微露脸| www国产在线视频色| 曰老女人黄片| 在线永久观看黄色视频| av片东京热男人的天堂| 精品电影一区二区在线| 亚洲成av人片免费观看| 黄网站色视频无遮挡免费观看| 精品免费久久久久久久清纯| 好看av亚洲va欧美ⅴa在| 久久久水蜜桃国产精品网| 国产片内射在线| 丁香欧美五月| 午夜老司机福利片| 999久久久精品免费观看国产| 欧美成人性av电影在线观看| 免费在线观看视频国产中文字幕亚洲| 99在线人妻在线中文字幕| 亚洲欧洲精品一区二区精品久久久| 午夜影院日韩av| 免费在线观看黄色视频的| 精品高清国产在线一区| 亚洲国产欧美日韩在线播放| 在线国产一区二区在线| 在线观看www视频免费| 88av欧美| av片东京热男人的天堂| 日本五十路高清| 一级作爱视频免费观看| 老司机靠b影院| 视频在线观看一区二区三区| 一级a爱视频在线免费观看| 18禁美女被吸乳视频| 亚洲av成人一区二区三| 一个人免费在线观看的高清视频| 日日摸夜夜添夜夜添小说| 1024香蕉在线观看| 黄色丝袜av网址大全| 精华霜和精华液先用哪个| 国产日本99.免费观看| 亚洲免费av在线视频| 欧美三级亚洲精品| 一个人观看的视频www高清免费观看 | 免费无遮挡裸体视频| 国内毛片毛片毛片毛片毛片| 亚洲午夜理论影院| 午夜福利免费观看在线| 色播在线永久视频| 久久久精品欧美日韩精品| 国产激情欧美一区二区| 欧美日本视频| 国产精品久久久av美女十八| 成年免费大片在线观看| 人人妻,人人澡人人爽秒播| 十八禁人妻一区二区| 婷婷亚洲欧美| 男女视频在线观看网站免费 | 99久久久亚洲精品蜜臀av| 又大又爽又粗| 久久香蕉精品热| 免费人成视频x8x8入口观看| 岛国在线观看网站| 免费人成视频x8x8入口观看| 老司机在亚洲福利影院| 色精品久久人妻99蜜桃| 国产精品亚洲一级av第二区| 久久亚洲真实| 国产精品99久久99久久久不卡| 亚洲一码二码三码区别大吗| 白带黄色成豆腐渣| 亚洲精华国产精华精| 婷婷亚洲欧美| 亚洲精华国产精华精| 白带黄色成豆腐渣| 久久天躁狠狠躁夜夜2o2o| 老熟妇乱子伦视频在线观看| 男女视频在线观看网站免费 | 久久精品aⅴ一区二区三区四区| 欧美成人免费av一区二区三区| 国产亚洲精品第一综合不卡| 狂野欧美激情性xxxx| 日本五十路高清| 老司机午夜十八禁免费视频| 亚洲美女黄片视频| 很黄的视频免费| 日韩免费av在线播放| 精品熟女少妇八av免费久了| 18禁观看日本| 美女 人体艺术 gogo| 91字幕亚洲| 欧美日韩福利视频一区二区| 成人国产一区最新在线观看| 国产成人av教育| 一进一出抽搐gif免费好疼| 麻豆国产av国片精品| 在线观看www视频免费| 久久性视频一级片| 亚洲中文字幕一区二区三区有码在线看 | 亚洲精品一区av在线观看| 18禁黄网站禁片午夜丰满| 亚洲全国av大片| 99久久国产精品久久久| 男女午夜视频在线观看| 欧美一级毛片孕妇| 男女下面进入的视频免费午夜 | 国产三级在线视频| 国产日本99.免费观看| 日本五十路高清| e午夜精品久久久久久久| 婷婷精品国产亚洲av在线| 久久这里只有精品19| 亚洲 国产 在线| 精品无人区乱码1区二区| 久久久久亚洲av毛片大全| 精品一区二区三区四区五区乱码| 午夜福利18| www国产在线视频色| 久久国产精品男人的天堂亚洲| 国产v大片淫在线免费观看| 啦啦啦免费观看视频1| 一卡2卡三卡四卡精品乱码亚洲| 男人舔女人下体高潮全视频| 国产精品 欧美亚洲| av在线播放免费不卡| 法律面前人人平等表现在哪些方面| 久热这里只有精品99| 国产片内射在线| www.www免费av| 国产亚洲精品一区二区www| 久久中文字幕人妻熟女| 亚洲自偷自拍图片 自拍| 午夜亚洲福利在线播放| 手机成人av网站| 亚洲精华国产精华精| 国产精品亚洲一级av第二区| 久久久国产成人免费| 天天躁夜夜躁狠狠躁躁| 亚洲一区高清亚洲精品| 精品国产一区二区三区四区第35| 天堂影院成人在线观看| 国产av不卡久久| 真人做人爱边吃奶动态| 制服丝袜大香蕉在线| 757午夜福利合集在线观看| 国产伦人伦偷精品视频| 亚洲在线自拍视频| 曰老女人黄片| 91老司机精品| 999久久久精品免费观看国产| 天堂√8在线中文| 国产成人欧美在线观看| 精品福利观看| 成人免费观看视频高清| www.熟女人妻精品国产| 欧美+亚洲+日韩+国产| 这个男人来自地球电影免费观看| av电影中文网址| 久久中文看片网| 欧洲精品卡2卡3卡4卡5卡区| 亚洲三区欧美一区| 一个人观看的视频www高清免费观看 | 很黄的视频免费| 嫩草影视91久久| 午夜福利在线观看吧| 欧美国产精品va在线观看不卡| 欧美性猛交╳xxx乱大交人| 精品国内亚洲2022精品成人| 88av欧美| 观看免费一级毛片| 欧美黑人精品巨大| 人妻久久中文字幕网| 久久久国产欧美日韩av| 久久九九热精品免费| 人人妻人人看人人澡| 免费高清视频大片| 一区二区三区精品91| 最新在线观看一区二区三区| 国产亚洲精品第一综合不卡| av福利片在线| 国产精品国产高清国产av| 黄片小视频在线播放| 国产高清videossex| 久久狼人影院| 免费在线观看成人毛片| 侵犯人妻中文字幕一二三四区| 男女之事视频高清在线观看| 麻豆一二三区av精品| 欧美性猛交黑人性爽| 日韩国内少妇激情av| 1024视频免费在线观看| 色老头精品视频在线观看| 99riav亚洲国产免费| 久9热在线精品视频| 国产精品久久久久久亚洲av鲁大| 老司机午夜福利在线观看视频| 久久婷婷人人爽人人干人人爱| 日本一本二区三区精品| 天天添夜夜摸| 欧美性长视频在线观看| 18禁裸乳无遮挡免费网站照片 | 亚洲av第一区精品v没综合| 变态另类成人亚洲欧美熟女| 在线观看免费午夜福利视频| 美女免费视频网站| 日本一本二区三区精品| 亚洲性夜色夜夜综合| 欧美不卡视频在线免费观看 | 怎么达到女性高潮| 免费看美女性在线毛片视频| 国语自产精品视频在线第100页| 嫁个100分男人电影在线观看| 国产爱豆传媒在线观看 | 久久久久国产精品人妻aⅴ院| 十八禁人妻一区二区| 天天躁夜夜躁狠狠躁躁| 亚洲av电影不卡..在线观看| www.999成人在线观看| 人人妻人人看人人澡| 久久久久久久精品吃奶| 999精品在线视频| 桃色一区二区三区在线观看| 国产精品久久久久久精品电影 | 欧美成人性av电影在线观看| 亚洲中文字幕日韩| 国产日本99.免费观看| 日本 av在线| 脱女人内裤的视频| 亚洲七黄色美女视频| 亚洲国产高清在线一区二区三 | 悠悠久久av| 桃色一区二区三区在线观看| 免费av毛片视频| 黄色a级毛片大全视频| 亚洲狠狠婷婷综合久久图片| 一二三四在线观看免费中文在| 一级作爱视频免费观看| 美女午夜性视频免费| 亚洲在线自拍视频| 中文字幕另类日韩欧美亚洲嫩草| 一级黄色大片毛片| x7x7x7水蜜桃| 一区二区日韩欧美中文字幕| 亚洲 欧美一区二区三区| 叶爱在线成人免费视频播放| 亚洲人成网站在线播放欧美日韩| 热99re8久久精品国产| 午夜a级毛片| 男女那种视频在线观看| 听说在线观看完整版免费高清| 亚洲欧美精品综合一区二区三区| 一进一出抽搐动态| 一区二区日韩欧美中文字幕| 99国产精品99久久久久| 99国产精品一区二区蜜桃av| 身体一侧抽搐| 狂野欧美激情性xxxx| 久久久久久久久久黄片| 国产伦一二天堂av在线观看| 午夜福利视频1000在线观看| av中文乱码字幕在线| 精品久久久久久久久久久久久 | 日韩大尺度精品在线看网址| 男女之事视频高清在线观看| 国产一区二区三区视频了| 久久久国产成人精品二区| 国产精品,欧美在线| netflix在线观看网站| 18美女黄网站色大片免费观看| 伊人久久大香线蕉亚洲五| 日韩欧美在线二视频| 国产精品日韩av在线免费观看| 婷婷亚洲欧美| bbb黄色大片| 男女做爰动态图高潮gif福利片| 国产激情偷乱视频一区二区| 国产精品 国内视频| 嫩草影视91久久| 丰满人妻熟妇乱又伦精品不卡| 91成年电影在线观看| 又紧又爽又黄一区二区| 非洲黑人性xxxx精品又粗又长| av在线播放免费不卡| 国产精品野战在线观看| 亚洲五月天丁香| 国产高清视频在线播放一区| 美女 人体艺术 gogo| 亚洲美女黄片视频| 一本精品99久久精品77| 久久中文字幕一级| 一级a爱视频在线免费观看| 亚洲五月婷婷丁香| 18禁观看日本| 成人18禁在线播放| 国产精品精品国产色婷婷| 一区二区三区激情视频| 99热6这里只有精品| 中文字幕人成人乱码亚洲影| 法律面前人人平等表现在哪些方面| 日韩 欧美 亚洲 中文字幕| 久久午夜综合久久蜜桃| 亚洲精品色激情综合| 国产精品野战在线观看| 午夜福利在线观看吧| 人人妻人人澡欧美一区二区| 精品久久久久久久久久免费视频| 男女之事视频高清在线观看| 午夜视频精品福利| 成人18禁在线播放| 丝袜美腿诱惑在线| 999精品在线视频| 亚洲国产中文字幕在线视频| 亚洲狠狠婷婷综合久久图片| 在线观看www视频免费| 亚洲欧美精品综合久久99| 国产一区二区激情短视频| 一级黄色大片毛片| 欧美中文日本在线观看视频| 啦啦啦 在线观看视频| 一级a爱视频在线免费观看| e午夜精品久久久久久久| 中文字幕人妻丝袜一区二区| 女性被躁到高潮视频| 可以在线观看的亚洲视频| 国产欧美日韩一区二区三| 国产免费男女视频| 可以在线观看的亚洲视频| 午夜免费观看网址| 欧美性长视频在线观看| 18禁黄网站禁片午夜丰满| 亚洲成人国产一区在线观看| 日韩国内少妇激情av| 操出白浆在线播放| 国产精品,欧美在线| 欧美精品啪啪一区二区三区| 三级毛片av免费| 美女高潮喷水抽搐中文字幕| 欧美不卡视频在线免费观看 | 国产熟女xx| 好男人在线观看高清免费视频 | 麻豆一二三区av精品| 午夜免费激情av| 国产在线观看jvid| 中文字幕高清在线视频| a在线观看视频网站| 国产成人系列免费观看| 成人精品一区二区免费| 久久性视频一级片| 国产主播在线观看一区二区| 91九色精品人成在线观看| 麻豆成人午夜福利视频| 欧美午夜高清在线|