• <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
    热re99久久精品国产66热6| 国产aⅴ精品一区二区三区波| 老鸭窝网址在线观看| 国产精品1区2区在线观看. | 午夜激情久久久久久久| 他把我摸到了高潮在线观看 | 国产亚洲精品久久久久5区| 丝袜喷水一区| 欧美 亚洲 国产 日韩一| 国产国语露脸激情在线看| 99精国产麻豆久久婷婷| 在线天堂中文资源库| 人妻一区二区av| 久久国产精品男人的天堂亚洲| 国产在线观看jvid| www.自偷自拍.com| 美女福利国产在线| 国产熟女午夜一区二区三区| 91av网站免费观看| 热99久久久久精品小说推荐| 丝袜喷水一区| 国产淫语在线视频| 国产成人欧美| 精品高清国产在线一区| 国产成人av教育| 电影成人av| 免费观看av网站的网址| 天堂中文最新版在线下载| 日韩欧美免费精品| 男人舔女人的私密视频| 免费在线观看日本一区| 国产单亲对白刺激| 欧美成狂野欧美在线观看| 欧美日韩视频精品一区| 久久久国产成人免费| 国产成人免费观看mmmm| 精品午夜福利视频在线观看一区 | 精品一区二区三区四区五区乱码| 亚洲av日韩精品久久久久久密| 久久亚洲精品不卡| 色在线成人网| 露出奶头的视频| 久久久精品区二区三区| 新久久久久国产一级毛片| 他把我摸到了高潮在线观看 | 午夜精品国产一区二区电影| av在线播放免费不卡| 亚洲精品美女久久久久99蜜臀| 亚洲精品久久成人aⅴ小说| 精品国产亚洲在线| 窝窝影院91人妻| 91国产中文字幕| 男女免费视频国产| 欧美精品人与动牲交sv欧美| 国产av一区二区精品久久| 嫩草影视91久久| 午夜视频精品福利| 国产亚洲精品第一综合不卡| 1024香蕉在线观看| 国产精品av久久久久免费| 这个男人来自地球电影免费观看| 亚洲免费av在线视频| 高潮久久久久久久久久久不卡| 亚洲自偷自拍图片 自拍| 欧美激情久久久久久爽电影 | 别揉我奶头~嗯~啊~动态视频| 亚洲精华国产精华精| 亚洲五月婷婷丁香| 一区二区三区精品91| 99精品久久久久人妻精品| 国产麻豆69| 免费观看a级毛片全部| e午夜精品久久久久久久| 国产精品免费大片| 精品第一国产精品| 成人国语在线视频| 国产黄色免费在线视频| 狠狠婷婷综合久久久久久88av| 精品人妻熟女毛片av久久网站| 精品国产一区二区三区四区第35| 丁香六月天网| 搡老岳熟女国产| 一区二区三区精品91| 人人妻人人爽人人添夜夜欢视频| 久久亚洲精品不卡| 999久久久国产精品视频| 亚洲成人免费电影在线观看| 97人妻天天添夜夜摸| 午夜激情av网站| 久久热在线av| 制服人妻中文乱码| bbb黄色大片| 在线av久久热| 日韩中文字幕视频在线看片| 精品第一国产精品| 一区二区三区乱码不卡18| 女同久久另类99精品国产91| 久久精品国产a三级三级三级| avwww免费| 女性生殖器流出的白浆| 高清黄色对白视频在线免费看| 性高湖久久久久久久久免费观看| 国产又色又爽无遮挡免费看| 好男人电影高清在线观看| 久久久久久人人人人人| av电影中文网址| 又黄又粗又硬又大视频| 一二三四在线观看免费中文在| 欧美日韩亚洲高清精品| 中文字幕制服av| kizo精华| 啦啦啦 在线观看视频| 老司机影院毛片| 国产日韩欧美在线精品| 丝瓜视频免费看黄片| 99国产极品粉嫩在线观看| 国产在线一区二区三区精| 老熟女久久久| 考比视频在线观看| 91成年电影在线观看| 久久久久久久久久久久大奶| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲综合色网址| 18禁裸乳无遮挡动漫免费视频| 亚洲熟女精品中文字幕| 久久免费观看电影| 99精国产麻豆久久婷婷| 成人免费观看视频高清| 美女高潮到喷水免费观看| 性色av乱码一区二区三区2| 99riav亚洲国产免费| 中文字幕色久视频| 色尼玛亚洲综合影院| 91麻豆av在线| 啪啪无遮挡十八禁网站| 欧美激情高清一区二区三区| 叶爱在线成人免费视频播放| 日本a在线网址| 汤姆久久久久久久影院中文字幕| 狠狠狠狠99中文字幕| 性色av乱码一区二区三区2| 一夜夜www| 国产xxxxx性猛交| 欧美变态另类bdsm刘玥| 亚洲,欧美精品.| 少妇裸体淫交视频免费看高清 | 久久精品国产综合久久久| 精品久久久久久久毛片微露脸| 午夜福利,免费看| 女性生殖器流出的白浆| 中文字幕人妻丝袜制服| 国产成人精品久久二区二区91| 欧美日韩成人在线一区二区| 三级毛片av免费| www.精华液| 一区二区三区国产精品乱码| 国产欧美日韩一区二区三区在线| 国产精品国产av在线观看| 五月开心婷婷网| 汤姆久久久久久久影院中文字幕| 久久久国产一区二区| 国产1区2区3区精品| 欧美亚洲日本最大视频资源| 午夜福利在线免费观看网站| 中文字幕最新亚洲高清| 少妇 在线观看| 国产野战对白在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 女人高潮潮喷娇喘18禁视频| 在线观看人妻少妇| 涩涩av久久男人的天堂| 丝袜喷水一区| 欧美黑人欧美精品刺激| 精品久久久久久久毛片微露脸| 亚洲熟女精品中文字幕| 法律面前人人平等表现在哪些方面| 精品卡一卡二卡四卡免费| 亚洲精品在线美女| 欧美午夜高清在线| 国产亚洲精品久久久久5区| 国产精品久久久久久人妻精品电影 | 中文字幕另类日韩欧美亚洲嫩草| 高清视频免费观看一区二区| 母亲3免费完整高清在线观看| 久久久国产一区二区| 国产1区2区3区精品| 激情视频va一区二区三区| 新久久久久国产一级毛片| 久久精品人人爽人人爽视色| 如日韩欧美国产精品一区二区三区| 国产成人av激情在线播放| 热re99久久精品国产66热6| 国产精品免费一区二区三区在线 | 999久久久国产精品视频| 色在线成人网| 免费少妇av软件| 男女无遮挡免费网站观看| 午夜福利一区二区在线看| 女性生殖器流出的白浆| 大型av网站在线播放| 精品一区二区三区四区五区乱码| 天堂8中文在线网| 下体分泌物呈黄色| 精品亚洲成国产av| 午夜福利在线免费观看网站| 99香蕉大伊视频| 正在播放国产对白刺激| 水蜜桃什么品种好| 男女无遮挡免费网站观看| 99国产精品一区二区蜜桃av | 午夜福利免费观看在线| 少妇裸体淫交视频免费看高清 | 91麻豆精品激情在线观看国产 | 在线观看www视频免费| 精品视频人人做人人爽| 免费观看av网站的网址| 色婷婷av一区二区三区视频| 少妇猛男粗大的猛烈进出视频| 在线观看66精品国产| 一边摸一边做爽爽视频免费| 12—13女人毛片做爰片一| 欧美+亚洲+日韩+国产| 80岁老熟妇乱子伦牲交| 亚洲七黄色美女视频| 在线观看免费午夜福利视频| 十八禁人妻一区二区| 天天添夜夜摸| 麻豆国产av国片精品| 怎么达到女性高潮| 999久久久精品免费观看国产| 亚洲国产看品久久| 国产精品.久久久| 又大又爽又粗| 精品久久蜜臀av无| 美女福利国产在线| 欧美日韩亚洲高清精品| 午夜精品国产一区二区电影| 色老头精品视频在线观看| 精品少妇久久久久久888优播| 色老头精品视频在线观看| 免费少妇av软件| 91成年电影在线观看| 久久精品亚洲av国产电影网| 极品人妻少妇av视频| 淫妇啪啪啪对白视频| 国产视频一区二区在线看| 男女无遮挡免费网站观看| 国产伦人伦偷精品视频| 国产一区二区在线观看av| 久久精品国产综合久久久| 国产男靠女视频免费网站| 一个人免费在线观看的高清视频| 人人妻人人澡人人看| 国产亚洲午夜精品一区二区久久| 黄色毛片三级朝国网站| 999久久久精品免费观看国产| 18禁观看日本| av欧美777| 色94色欧美一区二区| 久久九九热精品免费| svipshipincom国产片| 亚洲欧美精品综合一区二区三区| av欧美777| 一级,二级,三级黄色视频| 成人18禁高潮啪啪吃奶动态图| 麻豆乱淫一区二区| 91麻豆精品激情在线观看国产 | 久久九九热精品免费| 精品一区二区三区四区五区乱码| 久久久久久亚洲精品国产蜜桃av| 久久国产亚洲av麻豆专区| 亚洲三区欧美一区| 热99久久久久精品小说推荐| 高清毛片免费观看视频网站 | videos熟女内射| 国产精品自产拍在线观看55亚洲 | 精品国产一区二区三区久久久樱花| a级毛片在线看网站| 老司机午夜十八禁免费视频| 免费在线观看日本一区| 伊人久久大香线蕉亚洲五| 亚洲成av片中文字幕在线观看| 亚洲色图 男人天堂 中文字幕| 青青草视频在线视频观看| 国产高清视频在线播放一区| 19禁男女啪啪无遮挡网站| 视频区图区小说| svipshipincom国产片| 美国免费a级毛片| 美女国产高潮福利片在线看| 中文字幕制服av| 国产精品二区激情视频| 国产成人影院久久av| 高清av免费在线| 亚洲精品在线美女| 视频区欧美日本亚洲| 国产一区二区在线观看av| 男男h啪啪无遮挡| 激情视频va一区二区三区| 少妇精品久久久久久久| 91精品三级在线观看| 51午夜福利影视在线观看| 国产国语露脸激情在线看| 我的亚洲天堂| 国产深夜福利视频在线观看| 91麻豆精品激情在线观看国产 | 久久人妻av系列| 午夜日韩欧美国产| 另类亚洲欧美激情| 久久九九热精品免费| 美女国产高潮福利片在线看| 久久毛片免费看一区二区三区| 久久久精品国产亚洲av高清涩受| 免费日韩欧美在线观看| 1024香蕉在线观看| 精品熟女少妇八av免费久了| 不卡av一区二区三区| 一本一本久久a久久精品综合妖精| 精品久久蜜臀av无| 日韩免费av在线播放| 国产黄色免费在线视频| 精品国内亚洲2022精品成人 | 免费在线观看影片大全网站| 精品免费久久久久久久清纯 | 欧美激情久久久久久爽电影 | 高清毛片免费观看视频网站 | 欧美在线一区亚洲| 日日爽夜夜爽网站| 80岁老熟妇乱子伦牲交| 在线观看免费午夜福利视频| 黄色成人免费大全| 欧美日韩中文字幕国产精品一区二区三区 | 香蕉丝袜av| 国产人伦9x9x在线观看| 少妇粗大呻吟视频| 91av网站免费观看| 国产高清国产精品国产三级| 欧美日韩亚洲高清精品| 夜夜爽天天搞| 亚洲av片天天在线观看| 精品视频人人做人人爽| 国产成人精品无人区| 窝窝影院91人妻| 午夜福利视频在线观看免费| 中文字幕人妻熟女乱码| www日本在线高清视频| 亚洲精品粉嫩美女一区| 国产成人啪精品午夜网站| 日韩中文字幕欧美一区二区| 18禁美女被吸乳视频| 天堂8中文在线网| 99re在线观看精品视频| 日本av免费视频播放| 在线 av 中文字幕| 纯流量卡能插随身wifi吗| 51午夜福利影视在线观看| 久久人妻av系列| videos熟女内射| 一边摸一边抽搐一进一小说 | 国产成人av教育| 亚洲第一欧美日韩一区二区三区 | 操美女的视频在线观看| 我的亚洲天堂| 中文亚洲av片在线观看爽 | 国产日韩欧美亚洲二区| 国产精品自产拍在线观看55亚洲 | 如日韩欧美国产精品一区二区三区| 久久青草综合色| 纯流量卡能插随身wifi吗| 久久久久久久精品吃奶| 国产日韩一区二区三区精品不卡| 三上悠亚av全集在线观看| 性色av乱码一区二区三区2| 黑人欧美特级aaaaaa片| 久久精品aⅴ一区二区三区四区| 欧美成人免费av一区二区三区 | 精品人妻在线不人妻| 美女扒开内裤让男人捅视频| 日韩精品免费视频一区二区三区| 国产成人啪精品午夜网站| 一区福利在线观看| 正在播放国产对白刺激| 男女床上黄色一级片免费看| 日日夜夜操网爽| 久久久久精品人妻al黑| 国产免费视频播放在线视频| 蜜桃在线观看..| 久久久久久亚洲精品国产蜜桃av| 黄色视频在线播放观看不卡| 亚洲成人免费av在线播放| 肉色欧美久久久久久久蜜桃| 亚洲欧洲精品一区二区精品久久久| 男人舔女人的私密视频| 久久av网站| 一区二区av电影网| 老司机深夜福利视频在线观看| av欧美777| 91国产中文字幕| 午夜福利影视在线免费观看| 男女无遮挡免费网站观看| 精品熟女少妇八av免费久了| www日本在线高清视频| 91麻豆精品激情在线观看国产 | 亚洲色图 男人天堂 中文字幕| 免费在线观看日本一区| 亚洲国产毛片av蜜桃av| 国产欧美日韩一区二区精品| 99国产精品一区二区三区| 午夜福利免费观看在线| 正在播放国产对白刺激| 男女无遮挡免费网站观看| 国产av又大| 香蕉久久夜色| 国产91精品成人一区二区三区 | 夜夜夜夜夜久久久久| 国产在线精品亚洲第一网站| 成人国产一区最新在线观看| 精品福利永久在线观看| 午夜久久久在线观看| 9热在线视频观看99| 免费在线观看完整版高清| 国产精品秋霞免费鲁丝片| 精品少妇内射三级| 国产在线观看jvid| 国产欧美日韩一区二区精品| www日本在线高清视频| 99re6热这里在线精品视频| 性少妇av在线| 国产成人欧美| av超薄肉色丝袜交足视频| 黄色丝袜av网址大全| 亚洲色图综合在线观看| 精品国产一区二区三区四区第35| 欧美人与性动交α欧美精品济南到| 少妇 在线观看| 夜夜爽天天搞| 91精品国产国语对白视频| 日韩欧美三级三区| 欧美中文综合在线视频| 精品久久久精品久久久| 国产野战对白在线观看| av福利片在线| 热re99久久精品国产66热6| 母亲3免费完整高清在线观看| 日韩一卡2卡3卡4卡2021年| 久久国产亚洲av麻豆专区| 国产欧美日韩一区二区三| 欧美一级毛片孕妇| 好男人电影高清在线观看| 如日韩欧美国产精品一区二区三区| 黄色成人免费大全| 最近最新中文字幕大全免费视频| 丰满少妇做爰视频| 亚洲精品国产精品久久久不卡| 久久午夜综合久久蜜桃| 成年动漫av网址| 欧美黄色淫秽网站| 99久久人妻综合| 精品国产乱码久久久久久小说| 女同久久另类99精品国产91| 国产成人精品无人区| 国产成人精品久久二区二区免费| 国产欧美日韩精品亚洲av| 久久影院123| 久久精品aⅴ一区二区三区四区| 狂野欧美激情性xxxx| 免费看a级黄色片| av福利片在线| 亚洲精品国产区一区二| 免费看十八禁软件| 欧美日韩精品网址| 男女下面插进去视频免费观看| 夜夜骑夜夜射夜夜干| 老司机影院毛片| 国产99久久九九免费精品| 国产91精品成人一区二区三区 | 黄网站色视频无遮挡免费观看| 电影成人av| 欧美激情久久久久久爽电影 | 黄色片一级片一级黄色片| 三上悠亚av全集在线观看| 久9热在线精品视频| 动漫黄色视频在线观看| 国产欧美日韩一区二区精品| 水蜜桃什么品种好| 国产成人精品久久二区二区91| 久久久久网色| 大型av网站在线播放| 欧美人与性动交α欧美软件| 俄罗斯特黄特色一大片| 脱女人内裤的视频| 中文字幕人妻丝袜一区二区| 精品第一国产精品| 国产91精品成人一区二区三区 | 日韩 欧美 亚洲 中文字幕| 老司机午夜福利在线观看视频 | 天天操日日干夜夜撸| 亚洲午夜理论影院| 国产三级黄色录像| 日日摸夜夜添夜夜添小说| 成年人免费黄色播放视频| 我的亚洲天堂| 欧美黑人欧美精品刺激| 青青草视频在线视频观看| 午夜福利,免费看| 久久久久久久久免费视频了| 1024香蕉在线观看| 熟女少妇亚洲综合色aaa.| av在线播放免费不卡| 一级片'在线观看视频| 好男人电影高清在线观看| 99久久国产精品久久久| 亚洲国产欧美日韩在线播放| 免费观看a级毛片全部| 亚洲国产欧美日韩在线播放| 香蕉国产在线看| 热re99久久精品国产66热6| 青草久久国产| 国产无遮挡羞羞视频在线观看| 国产精品成人在线| 老司机在亚洲福利影院| 久久免费观看电影| 国产视频一区二区在线看| 高清毛片免费观看视频网站 | 欧美精品一区二区大全| 欧美精品亚洲一区二区| 69av精品久久久久久 | cao死你这个sao货| 国产精品久久久久久精品电影小说| 电影成人av| 中文亚洲av片在线观看爽 | 国产精品秋霞免费鲁丝片| 伊人久久大香线蕉亚洲五| 国产成人一区二区三区免费视频网站| 丝瓜视频免费看黄片| 日本撒尿小便嘘嘘汇集6| 无遮挡黄片免费观看| 可以免费在线观看a视频的电影网站| 欧美av亚洲av综合av国产av| av又黄又爽大尺度在线免费看| 久久狼人影院| 国产精品久久电影中文字幕 | 纯流量卡能插随身wifi吗| 免费在线观看黄色视频的| 亚洲 国产 在线| 亚洲伊人久久精品综合| 国产野战对白在线观看| 一本—道久久a久久精品蜜桃钙片| 丁香欧美五月| √禁漫天堂资源中文www| 91麻豆精品激情在线观看国产 | 天堂中文最新版在线下载| 美女国产高潮福利片在线看| 久久免费观看电影| 午夜两性在线视频| 亚洲五月婷婷丁香| 热re99久久精品国产66热6| 中文字幕人妻熟女乱码| 丝瓜视频免费看黄片| 99热网站在线观看| 国产日韩欧美在线精品| 欧美中文综合在线视频| 久久天堂一区二区三区四区| 国产在线精品亚洲第一网站| 99在线人妻在线中文字幕 | 久久精品亚洲熟妇少妇任你| 色精品久久人妻99蜜桃| 999久久久精品免费观看国产| 国产熟女午夜一区二区三区| 午夜视频精品福利| 黄色片一级片一级黄色片| 精品亚洲乱码少妇综合久久| 国产97色在线日韩免费| 丝袜喷水一区| 国产亚洲欧美精品永久| 日韩精品免费视频一区二区三区| 法律面前人人平等表现在哪些方面| 宅男免费午夜| 99精品在免费线老司机午夜| 国产欧美亚洲国产| 热re99久久国产66热| 香蕉丝袜av| www日本在线高清视频| a级毛片在线看网站| 亚洲专区国产一区二区| 91国产中文字幕| 老熟女久久久| 亚洲avbb在线观看| 人人妻人人爽人人添夜夜欢视频| 亚洲中文日韩欧美视频| 免费观看人在逋| 午夜激情av网站| 99香蕉大伊视频| 国产日韩欧美在线精品| 9色porny在线观看| 91成年电影在线观看| 久久久久网色| 亚洲色图 男人天堂 中文字幕| 热99国产精品久久久久久7| 欧美老熟妇乱子伦牲交| 脱女人内裤的视频| 99国产综合亚洲精品| 国产精品98久久久久久宅男小说| 麻豆国产av国片精品| 日本欧美视频一区| 高潮久久久久久久久久久不卡| 国产有黄有色有爽视频| 国产成+人综合+亚洲专区| 亚洲一区二区三区欧美精品| 大型黄色视频在线免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品一二三| 亚洲av成人一区二区三| 久久久久久久精品吃奶| 亚洲欧洲日产国产|