——以廣東和河北地區(qū)為例"/>
  • <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Pn-Pg震相組合測(cè)定震源深度方法的應(yīng)用研究
    ——以廣東和河北地區(qū)為例

    2016-08-05 08:57:00劉雙慶
    物探化探計(jì)算技術(shù) 2016年3期

    謝 靜, 劉雙慶

    (天津市地震局,天津 300201)

    ?

    Pn-Pg震相組合測(cè)定震源深度方法的應(yīng)用研究
    ——以廣東和河北地區(qū)為例

    謝靜, 劉雙慶

    (天津市地震局,天津300201)

    摘要:在直達(dá)波射線走時(shí)算法改良的基礎(chǔ)上,利用天然地震初至震相Pn、Pg進(jìn)行組合的方法,對(duì)區(qū)域地震震源深度進(jìn)行重新定位,結(jié)果顯示,廣東臺(tái)網(wǎng)觀測(cè)報(bào)告里震源深度和這里新方法定位深度總體上吻合較好。河北臺(tái)網(wǎng)觀測(cè)報(bào)告里震源深度普遍偏淺3 km~6 km,新定位結(jié)果顯示唐山老震區(qū)還存在下地殼地震。Pn-Pg震相組合法對(duì)唐山老震區(qū)的震源深度定位效果一般更深些,同時(shí)也進(jìn)一步反映了斷層構(gòu)造信息。

    關(guān)鍵詞:震源深度; Pn-Pg震相組合法; 華北平原地震帶; 震相報(bào)告

    0引言

    震源深度是描述震源的最基本參數(shù)之一,在研究地震活動(dòng)與斷層的關(guān)系[1]、地球動(dòng)力學(xué)意義上的板塊運(yùn)動(dòng)和殼幔構(gòu)造[2-4]、地震成核和震源破裂[9]及地震視應(yīng)力隨深度變化特征[10]等方面,獲得準(zhǔn)確的震源深度是不可忽視的工作。同時(shí),震源深度這個(gè)重要的參數(shù)也是目前最難確定的參數(shù)之一[5]。目前,省級(jí)日常地震定位工作中,主要使用MSDP軟件中掛載的單純形法、LocalSat、HypoSat方法進(jìn)行絕對(duì)定位;遺傳算法等非線性搜索定位程序。但每種方法給出的效果差異較大,特別是對(duì)單純形法的深度定位結(jié)果大多數(shù)人普遍反映偏淺(圖1),而且給出的深度空間分布與當(dāng)?shù)氐刭|(zhì)構(gòu)造沒(méi)有較好的關(guān)聯(lián)性。因此,研究人員對(duì)各省局給出的地震震相報(bào)告進(jìn)行了校對(duì),并利用主事件、雙差定位等相對(duì)定位方法對(duì)小區(qū)域地震群進(jìn)行深度定位,結(jié)合層析成像資料勾畫相應(yīng)的發(fā)震構(gòu)造[8,11],但相對(duì)定位法要求震源之間距離遠(yuǎn)小于震中距且定位結(jié)果還存在系統(tǒng)性偏差現(xiàn)象。研究人員用深度震相,如sPg,sPmP,sPn進(jìn)行深度測(cè)定,該方法雖然原理上較其他方法在結(jié)果敏感度方面較高,但因這些震相是后續(xù)震相,信噪比很低,識(shí)別困難,所以也沒(méi)有得到有效應(yīng)用[7]。作者在前人提出的算法基礎(chǔ)上[6],進(jìn)一步優(yōu)化了直達(dá)波走時(shí)算法,對(duì)省局的地震震相報(bào)告用Pn-Pg震源深度重定位,以進(jìn)一步分析震相報(bào)告中震源深度信息的可靠性,并對(duì)一些顯著地震構(gòu)造的震源深度分布進(jìn)行討論。

    1.Pn-Pg震相組合法

    1.1Pn-Pg震相組合法的靈敏度優(yōu)勢(shì)

    Pn-Pg震相組合深度定位方法由朱元清等[6]提出。當(dāng)震中位置確定時(shí),不同震源深度到同一臺(tái)站的Pg、Pn震相走時(shí)是不同的,因此可以通過(guò)不同臺(tái)站的這兩個(gè)震相及已有震中位置進(jìn)行深度測(cè)定。該方法既使用了初至震相到時(shí)提取精度高的優(yōu)點(diǎn),同時(shí)又利用了Pn、Pg震相離源射線反向?qū)ι疃茸兓舾械奶攸c(diǎn),而且不受發(fā)震時(shí)刻測(cè)定誤差的影響,從而可以有效提高深度定位精度。其中,Pn、Pg震相組合相對(duì)于只用直達(dá)Pg、Sg震相組合,其提高深度變化對(duì)走時(shí)差異影響的靈敏度分析如下:

    圖1 華北地區(qū)的地震深度分布隨臺(tái)網(wǎng)產(chǎn)出時(shí)間的變化Fig.1 The distribution of seismic depth in north China changes with network output time

    為簡(jiǎn)化理論推導(dǎo),考慮一層地殼模型,假定震中距為Δ,震源深度為h,發(fā)震時(shí)刻為T0,Tpg、Tsg、Tpn分別是直達(dá)波P波、S波和首波的到時(shí),Vpg、Vsg是直達(dá)波P波、S波的速度,則走時(shí)方程如式(1)~式(3)所示。

    (1)

    (2)

    (3)

    對(duì)式(1)~式(3)分別求關(guān)于h的偏導(dǎo)數(shù),將式(2)的求導(dǎo)結(jié)果與式(1)的求導(dǎo)結(jié)果相減,以及將式(3)的求導(dǎo)結(jié)果與式(1)的求導(dǎo)結(jié)果相減,將二者的差求商得式(4)。

    (4)

    假定Vpg=1.73×VSg=6km/s,Vpn=8 km/s,則κ值見(jiàn)表1。從表1中可以看出,用Pn、Pg震相組合相對(duì)Pg、Sg震相組合,深度變化對(duì)兩者走時(shí)差的影響在2倍以上——對(duì)于Pn處于盲區(qū)的情形,可以通過(guò)沿Moho面滑移來(lái)處理[6]。震源深度越淺震中距越大時(shí),這里的震相組合方法對(duì)深度變化的敏感性更大。

    表1 不同震中距不同深度時(shí)的靈敏度κ分布

    1.2Pg射線追蹤的算法優(yōu)化

    兩點(diǎn)間的射線追蹤,即確定震源與接收點(diǎn)之間的射線路徑和走時(shí),是諸如地震定位、走時(shí)反演、合成理論地震圖、地震層析成像等許多地震學(xué)問(wèn)題中的基本問(wèn)題,也是最耗時(shí)的問(wèn)題,因此,田玥[12],提出了一種快速準(zhǔn)確的水平層狀介質(zhì)中兩點(diǎn)間的射線追蹤方法該方法對(duì)任意復(fù)雜的水平層狀速度結(jié)構(gòu)、任意大小的入射角都能穩(wěn)定、快速、準(zhǔn)確地收斂。

    給定震源和接收點(diǎn)的位置及速度模型,確定直達(dá)波和反射波的射線路徑,射線路徑由射線參數(shù)p唯一確定,p可以由震中距方程直接求得。

    (5)

    1)將射線參數(shù)p轉(zhuǎn)換為新的參數(shù)q。不直接求解射線參數(shù)p,而引入新參數(shù)q來(lái)代替p。q是最高速度層(其等效厚度為hM)中射線走過(guò)的水平距離(圖2),則p與q之間有簡(jiǎn)單的轉(zhuǎn)換關(guān)系為式(6)。

    (6)

    將式(6)代入式(5)得式(7)。

    (7)

    式中,εk=vk/vM。

    2)用牛頓法可求解關(guān)于q的震中距方程,對(duì)于給定的震中距Δ0,方程(7)可寫成關(guān)于q的非線性方程:

    圖2 入射角接近90°時(shí)的射線路徑Fig.2 The ray path when the incident angle is closed to 90 °

    (8)

    用牛頓法求解方程(8)的迭代公式為:

    (9)

    3)將得到的新參數(shù)q轉(zhuǎn)換為射線參數(shù)p。作者按上述思路,重新編寫了水平層狀介質(zhì)直達(dá)波走時(shí)計(jì)算程序,并進(jìn)行了程序可靠性檢驗(yàn)。給出式(5)、式(7)兩種方法的Pg射線追蹤算法的數(shù)值算例,速度結(jié)構(gòu)見(jiàn)表2,震中距50 km,3種震源深度所得的計(jì)算結(jié)果,用牛頓法求解p和求解q的迭代收斂情況比較如圖3所示??梢钥闯?,在一般情況下,直接求解p收斂較慢,需要迭代20次左右,然而求解q收斂非常迅速,只需要迭代2次~3次,并且入射角越接近90°時(shí),p收斂得越慢(迭代32次),q收斂得越快(迭代2次)。

    表2 分層速度模型

    2區(qū)域臺(tái)網(wǎng)震相資料震源深度重定位

    2.1廣東臺(tái)網(wǎng)震源深度重定位

    中國(guó)地震局“十五”期間的測(cè)震臺(tái)網(wǎng)項(xiàng)目中,向各省局推廣的Jopens系統(tǒng)(MSDP軟件為其子模塊)及相關(guān)配置由廣東省的相關(guān)技術(shù)負(fù)責(zé)人進(jìn)行研發(fā),他們對(duì)相關(guān)的技術(shù)細(xì)致掌握相對(duì)到位,資料產(chǎn)出的可靠性更高。特別是地殼模型中的華南走時(shí)表模型是針對(duì)該區(qū)域的,因此多方面的不確定性因素相對(duì)較少,重定位前、后震源深度應(yīng)該相差較小。因此這里首先利用廣東省的震相報(bào)告及華南地殼模型進(jìn)行深度重定位,分析Pn-Pg震相組合方法的適用性。對(duì)2009年—2013年廣東臺(tái)網(wǎng)的所有ML2.5以上地震進(jìn)行震源深度重定位,結(jié)果見(jiàn)圖4。直方圖顯示新計(jì)算的結(jié)果主值在10 km~14 km,原結(jié)果主值在9 km~12 km,新方法中個(gè)別未能測(cè)定深度的結(jié)果未畫出,直方圖主值分布相差2 km左右,但總體上,兩者震源深度集中區(qū)吻合較好。

    圖3 用牛頓法求解p和q的收斂性比較Fig.3 Comparison the convergence of using Newton method to solve p and q(a)震源深度為16 km時(shí)兩種方法的收斂情況;(b)震源深度為20 km時(shí)兩種方法的收斂情況;(c)震源深度為29 km時(shí)兩種方法的收斂情況

    圖4 廣東省2009年-2013年ML2.5以上地震深度重定位對(duì)比Fig.4 Comparison seismic depth after relocation which magnitude are greater than 2.5 in Guangdong province from 2009 to 2013

    2.2河北臺(tái)網(wǎng)震源深度重定位

    河北臺(tái)網(wǎng)是首都圈地區(qū)的老臺(tái)網(wǎng),臺(tái)網(wǎng)中部區(qū)域位于華北平原沉積區(qū),部分沉積厚度在5 km以上。該地區(qū)存在1966年邢臺(tái)Ms7.2、1976年唐山Ms7.8兩個(gè)老震區(qū),在天津西南鄰區(qū)2006年曾發(fā)生文安Ms5.1地震,天津東部寧河1976年曾發(fā)生Ms6.9地震,天津西北鄰區(qū)1679年曾發(fā)生三河—平谷Ms8.0地震。河北臺(tái)網(wǎng)的臺(tái)間距相對(duì)全國(guó)臺(tái)網(wǎng)相對(duì)較小,平均約45 km。作者抽取了2009年1月至2011年12月ML2.0以上河北臺(tái)網(wǎng)的華北平原地震帶內(nèi)的震相報(bào)告進(jìn)行震源深度重定位。定位結(jié)果如圖5、圖6所示:①圖5(b)和圖5(d)對(duì)比分析可見(jiàn)原始深度過(guò)于集中于5 km~10 km,重定位后震源深度分布較為均勻且過(guò)渡平緩(圖6);②直方圖顯示河北臺(tái)網(wǎng)觀測(cè)報(bào)告震源深度主值在6 km,并且特別的集中于6 km,新定位結(jié)果主值在9 km~12 km范圍內(nèi),在21 km以下的下地殼,作者仍獲取到一些地震深度結(jié)果,而這些震源較深的地震多與老震區(qū)分布有關(guān);③由于研究地區(qū)的沉積蓋層比較厚,同時(shí)結(jié)合近幾年在測(cè)震臺(tái)網(wǎng)編目的經(jīng)驗(yàn),利用MSDP掛載的單純形法定位,研究地區(qū)幾乎沒(méi)有出現(xiàn)過(guò)震源深度大于15 km的情況(圖1),這可以從圖5(b)震源深度和經(jīng)度關(guān)系圖和圖6直方圖的原始深度分布中可以看出,因此作者認(rèn)為6 km的優(yōu)勢(shì)分布深度可能不能反映研究地區(qū)的真實(shí)震源深度特征。

    3認(rèn)識(shí)及討論

    天然地震震源初始破裂深度的確定一直是一個(gè)棘手的難題,其客觀上受地震震源位于臺(tái)站下方若干公里的地殼內(nèi)部,震源深度以下沒(méi)有臺(tái)站,無(wú)法實(shí)現(xiàn)立體觀測(cè)所導(dǎo)致。雖然現(xiàn)實(shí)物理?xiàng)l件限制了其觀測(cè)效果,但震源發(fā)射出來(lái)的地震波因其傳到各個(gè)臺(tái)站的特殊運(yùn)動(dòng)屬性,對(duì)其進(jìn)行適當(dāng)?shù)慕M合取舍,可以壓制敏感度低的測(cè)定組合,提高震源深度測(cè)定的效果。這里使用的Pn-Pg組合方法即是其中的一種,該方法原理直觀,實(shí)現(xiàn)便捷。

    圖6 華北平原地震帶震源深度重定位前后對(duì)比直方圖Fig.6 Histogram contrast before and after   seismic depth relocation in north China plain

    研究表明,該方法定位出的震源深度總體上比原報(bào)告中的結(jié)果偏深,廣東地區(qū)偏深2 km,河北地區(qū)偏深3 km~6 km,新計(jì)算的結(jié)果在深度分布上過(guò)渡比較平緩且未出現(xiàn)明顯的深度分布異常高的 “易震”層位。

    當(dāng)然這里的方法對(duì)各向異性比較突出的地區(qū)的深度定位效果不太好,對(duì)臺(tái)站高程相差較大,地殼橫向速度變化較大的區(qū)域,也沒(méi)有找到一個(gè)比較好的通用處理技術(shù)。由于該方法是在固定震中位置的前提下對(duì)深度進(jìn)行定位,如果臺(tái)網(wǎng)密度不高,震中定位精度不夠的地震事件,該方法對(duì)其深度定位結(jié)果的改善也比較有限。

    參考文獻(xiàn):

    [1]GALDEANO C.S.,CASADO C.L.,DELGADO J,et al.Shallow seismicity and active faults in the Betic Cordillera:a preliminary approach to seismic sources associated with specific faults[J].Tectonophysics,1995,248(3-4):293-302.

    [2]GLENNON M.A.,CHEN,W.P..Rupture of deep-focus earthquakes in the north-western Pacific and their implications on seismogenesis[J].Geophy,J.Inter.,1995,120(3):706-720

    [3]鄭斯華.青藏高原地震的震源深度及其構(gòu)造意義[J].中國(guó)地震,1995,11(2):99-106.

    ZHENG S H.The earthquake focal depth of Qinghai Tibet Plateau and its tectonic significance[J].Earthquake Research in China,1995,11(2):99-106.(In Chinese)

    [4]張國(guó)民,汪素云,李麗,等.中國(guó)大陸地震震源深度及其構(gòu)造含義[J].科學(xué)通報(bào),2002,47(9):663-668.

    ZHANG G M,WANG S Y,LI L,et al.The earthquake focal depth and its tectonic implications in the Chinese mainland[J].Chinese Science Bulletin,2002,47(9):663-668.(In Chinese)

    [5]高原,周蕙蘭,鄭斯華,等.測(cè)定震源深度的意義的初步討論[J].中國(guó)地震,1997,13(4):321-329.

    GAO Y,ZHOU H L,ZHENG S H,et al.A preliminary discussion on the significance of determining focal depth[J].Earthquake Research in China,1997,13(4):321-329.(In Chinese)

    [6]朱元清,石耀霖.一種確定地震深度的新方法[J].地震地磁觀測(cè)與研究,1990,11(2):4-12.

    ZHU Y Q,SHI Y L.A new method for determining focal depth [J].Seismological and Geomagnetic Observation and Research,1990,11(2):4-12.(In Chinese)

    [7]崇加軍,倪四道,曾祥方.sPL,一個(gè)近距離確定震源深度的震相[J].地球物理學(xué)報(bào),2010,53(11):2620-2630.

    CONG J J,NI S D,ZENG X F.sPL,an effective seismic phase for determining focal depth at near distance[J].Chinese Journal of Geophysics,2010,53(11):2620-2630.(In Chinese)

    [8]黃媛,楊建思,張?zhí)熘?2003年新疆巴楚—伽師地震序列的雙差法重新定位研究[J].地球物理學(xué)報(bào),2006,49(1):162-169.

    HUANG Y,YANG J S,ZHANG T Z.Relocation of the Ba chu-Jia shi,Xinjiang earthquake sequence in 2003 using the double-difference location algorithm[J].Chinese Journal of Geophysics,2006,49(1):162-169.(In Chinese)

    [9]許力生,陳運(yùn)泰.用經(jīng)驗(yàn)格林函數(shù)方法從長(zhǎng)周期數(shù)字波形資料中提取共和地震的震源時(shí)間函數(shù)[J].地震學(xué)報(bào),1996,18(2):156-169.

    XU L S,CHEN Y T.Using the empirical Green function method to extract the source time function of Gonghe earthquake from the long period digital waveform data[J].Acta seimologica Sinica,1996,18(2):156-169.(In Chinese)

    [10]吳忠良.震源深度與矩心深度的分布特征及其與震源機(jī)制類型之間的關(guān)系[J].中國(guó)地震,2002,18(4):337-345.

    WU Z L.The distribution features of focal depth and centroid depth and the relationship between the type of focal mechanism[J].Earthquake Research in China,2002,18(4):337-345.(In Chinese)

    [11]朱艾斕,徐錫偉,于貴華,等.玉樹地震序列重新定位及其地震構(gòu)造研究[J].地學(xué)前緣,2012,19(4):8-14.

    ZHU A L,XU X W,YU G H,et al.Reloaction of the Yushu Ms7.1 earthquake sequence and investigation on its seismotectonics[J].Earth Sceince Frontiers,2012,19(4):8-14.(In Chinese)

    [12]田玥,陳曉非.水平層狀介質(zhì)中的快速兩點(diǎn)間射線追蹤方法[J].地震學(xué)報(bào),2005,27(2):147-154.

    TIAN Y,CHRN X F.A Rapid and Accurate two-point Ray Tracing Method in Horizontally Layered Velocity Model[J].Acta seimologica Sinica,2005,27(2):147-154.(In Chinese)

    收稿日期:2015-06-04改回日期:2015-08-03

    基金項(xiàng)目:星火計(jì)劃(XH14002Y);中國(guó)地震局測(cè)震青年骨干項(xiàng)目(20130202)

    作者簡(jiǎn)介:謝靜( 1981- ),女,碩士,主要研究方向?yàn)樘烊坏卣穑珽-mail:juliaxj@163.com。

    文章編號(hào):1001-1749(2016)03-0361-06

    中圖分類號(hào):P 631.4

    文獻(xiàn)標(biāo)志碼:A

    DOI:10.3969/j.issn.1001-1749.2016.03.11

    Studying on the application of Pn-Pg seismic phases combination method to determine the regional hypocenter depth——take Guangdong and Hebei as an example

    XIE Jing, LIU Shuang-qing

    (Earthquake Administration of Tianjin Municipality ,Tianjin300201,China)

    Abstract:Based on the improved direct wave ray travel time algorithm, in this paper the regional earthquake focal depths are relocated by using the combination of first arrival phases of Pn and Pg. The relocation indicates that the focal depths in Guangdong network's observation reports are consistent with focal depths given by the new method. But the focal depths in Hebei network's observation report are generally shallower than 3-6 km. And the new results indicate that part of earthquakes also exist in lower crustal in Tangshan seismic belt. The new hypocenter depths yielded by the phase Pn-Pg combination method are generally deeper in Tangshan seismic belt, and also the results can reflect the seismic fault information in the further extension.

    Key words:focal depth; the method of Pn-Pg phases combination; seismic belt in North China plain; phases report

    久久青草综合色| 国产成人系列免费观看| 国产一区二区三区在线臀色熟女 | 男女午夜视频在线观看| 午夜福利在线免费观看网站| 人人妻人人澡人人看| 丰满迷人的少妇在线观看| 欧美日韩av久久| 国产野战对白在线观看| 精品国产乱码久久久久久男人| 国产欧美日韩综合在线一区二区| 日日夜夜操网爽| 妹子高潮喷水视频| 亚洲av美国av| 欧美精品一区二区免费开放| 搡老乐熟女国产| 一边摸一边抽搐一进一出视频| 亚洲少妇的诱惑av| 欧美激情高清一区二区三区| 亚洲国产欧美在线一区| 国产精品国产高清国产av | 久久国产精品大桥未久av| 精品国产乱码久久久久久小说| 淫妇啪啪啪对白视频| 日韩一区二区三区影片| kizo精华| 亚洲熟女精品中文字幕| 免费一级毛片在线播放高清视频 | 国产97色在线日韩免费| 欧美国产精品va在线观看不卡| 免费看a级黄色片| 女人久久www免费人成看片| 亚洲天堂av无毛| 人人妻人人添人人爽欧美一区卜| 纯流量卡能插随身wifi吗| 一二三四社区在线视频社区8| 久久精品国产亚洲av高清一级| 午夜福利免费观看在线| 丁香六月欧美| 看免费av毛片| 老司机午夜十八禁免费视频| 成人国产av品久久久| 多毛熟女@视频| 9色porny在线观看| 国产成人精品久久二区二区免费| av不卡在线播放| 九色亚洲精品在线播放| 亚洲精品中文字幕在线视频| 自线自在国产av| 国产日韩欧美在线精品| 国产亚洲欧美在线一区二区| 极品教师在线免费播放| 99riav亚洲国产免费| 国产色视频综合| 欧美黄色淫秽网站| 一区二区三区乱码不卡18| 午夜免费成人在线视频| 90打野战视频偷拍视频| 交换朋友夫妻互换小说| 极品教师在线免费播放| 欧美激情久久久久久爽电影 | 男男h啪啪无遮挡| 丝袜美足系列| 一边摸一边抽搐一进一出视频| 国产精品久久久久成人av| 大片免费播放器 马上看| 99久久99久久久精品蜜桃| 久久久久网色| 久久久精品区二区三区| 国产深夜福利视频在线观看| 国产淫语在线视频| 国产欧美日韩精品亚洲av| 久久人妻熟女aⅴ| 超色免费av| 操美女的视频在线观看| 成人影院久久| 两性午夜刺激爽爽歪歪视频在线观看 | 丝袜美腿诱惑在线| 日韩三级视频一区二区三区| 色94色欧美一区二区| 一级片免费观看大全| 国产欧美日韩一区二区精品| 亚洲成国产人片在线观看| 国产成人av激情在线播放| 一本综合久久免费| 热re99久久精品国产66热6| tube8黄色片| 一本一本久久a久久精品综合妖精| 精品国产乱码久久久久久男人| 免费人妻精品一区二区三区视频| 狂野欧美激情性xxxx| 午夜福利视频精品| 极品少妇高潮喷水抽搐| 母亲3免费完整高清在线观看| 亚洲精品国产色婷婷电影| 黄网站色视频无遮挡免费观看| 国产伦理片在线播放av一区| 亚洲国产av新网站| 777久久人妻少妇嫩草av网站| 欧美变态另类bdsm刘玥| 日本av免费视频播放| 精品国内亚洲2022精品成人 | 亚洲第一av免费看| 一级,二级,三级黄色视频| av天堂久久9| 建设人人有责人人尽责人人享有的| 亚洲久久久国产精品| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品亚洲精品国产色婷小说| 男女无遮挡免费网站观看| 免费看十八禁软件| 色播在线永久视频| 天天影视国产精品| 两性夫妻黄色片| 免费一级毛片在线播放高清视频 | 伦理电影免费视频| 好男人电影高清在线观看| 99国产综合亚洲精品| 亚洲国产欧美在线一区| 亚洲三区欧美一区| 免费观看人在逋| 大片电影免费在线观看免费| 飞空精品影院首页| 欧美日韩中文字幕国产精品一区二区三区 | 美女福利国产在线| 精品一品国产午夜福利视频| 十八禁人妻一区二区| 国产99久久九九免费精品| 99久久人妻综合| 男女之事视频高清在线观看| 最近最新中文字幕大全电影3 | 欧美在线黄色| 麻豆国产av国片精品| 亚洲欧美一区二区三区久久| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日本中文国产一区发布| 一本一本久久a久久精品综合妖精| 亚洲三区欧美一区| 日韩一区二区三区影片| 精品视频人人做人人爽| 性高湖久久久久久久久免费观看| av超薄肉色丝袜交足视频| 大片电影免费在线观看免费| 91精品国产国语对白视频| 少妇 在线观看| 免费少妇av软件| 国产精品自产拍在线观看55亚洲 | 两性夫妻黄色片| 免费av中文字幕在线| 色婷婷久久久亚洲欧美| 99国产极品粉嫩在线观看| 啦啦啦 在线观看视频| 母亲3免费完整高清在线观看| 老司机午夜福利在线观看视频 | 啦啦啦免费观看视频1| 天天添夜夜摸| 黄色片一级片一级黄色片| 国产成人av激情在线播放| 9热在线视频观看99| 女人高潮潮喷娇喘18禁视频| 亚洲精品粉嫩美女一区| 欧美精品一区二区免费开放| 搡老乐熟女国产| 国产精品免费视频内射| 在线观看免费视频网站a站| 少妇裸体淫交视频免费看高清 | 丰满人妻熟妇乱又伦精品不卡| av不卡在线播放| 国产野战对白在线观看| 国产成人精品在线电影| 免费在线观看黄色视频的| 国产99久久九九免费精品| 国产精品久久久久久精品古装| 亚洲精品一二三| 亚洲精华国产精华精| 高清av免费在线| 欧美国产精品va在线观看不卡| 一个人免费看片子| 叶爱在线成人免费视频播放| 精品久久久久久电影网| 成人精品一区二区免费| 亚洲全国av大片| 国产在线免费精品| 成人黄色视频免费在线看| 成人影院久久| 夜夜骑夜夜射夜夜干| 999精品在线视频| 欧美老熟妇乱子伦牲交| 亚洲av欧美aⅴ国产| 可以免费在线观看a视频的电影网站| av电影中文网址| 精品亚洲成a人片在线观看| 国产精品一区二区在线观看99| 亚洲色图 男人天堂 中文字幕| 成人18禁高潮啪啪吃奶动态图| 国产精品免费视频内射| 男女无遮挡免费网站观看| 一夜夜www| 日本撒尿小便嘘嘘汇集6| 美女福利国产在线| 国产欧美日韩一区二区三| 国产精品秋霞免费鲁丝片| 激情在线观看视频在线高清 | 久久中文字幕一级| 久久影院123| 丝瓜视频免费看黄片| 母亲3免费完整高清在线观看| 精品一区二区三区av网在线观看 | 亚洲少妇的诱惑av| 久久国产亚洲av麻豆专区| 日韩一区二区三区影片| 日韩 欧美 亚洲 中文字幕| 午夜福利在线免费观看网站| 操出白浆在线播放| 看免费av毛片| 999精品在线视频| 亚洲熟妇熟女久久| 亚洲人成电影免费在线| 亚洲av欧美aⅴ国产| 国产又色又爽无遮挡免费看| 久久天堂一区二区三区四区| 欧美精品人与动牲交sv欧美| 国产单亲对白刺激| 久久精品亚洲精品国产色婷小说| 丁香六月天网| 天天添夜夜摸| 精品久久蜜臀av无| 一级黄色大片毛片| 精品少妇黑人巨大在线播放| 99热国产这里只有精品6| 午夜福利视频精品| 久久人人97超碰香蕉20202| 波多野结衣av一区二区av| 欧美日韩亚洲综合一区二区三区_| 夜夜爽天天搞| 国产精品香港三级国产av潘金莲| 久久久国产精品麻豆| 视频区图区小说| 亚洲国产欧美网| 涩涩av久久男人的天堂| 久久精品国产亚洲av高清一级| 欧美另类亚洲清纯唯美| 变态另类成人亚洲欧美熟女 | 亚洲av日韩精品久久久久久密| 热99re8久久精品国产| 亚洲中文日韩欧美视频| 女人高潮潮喷娇喘18禁视频| 日本a在线网址| 多毛熟女@视频| 欧美日韩黄片免| 国产一区二区三区在线臀色熟女 | 精品人妻熟女毛片av久久网站| 国产成人免费观看mmmm| 妹子高潮喷水视频| a在线观看视频网站| 免费在线观看视频国产中文字幕亚洲| 欧美日韩中文字幕国产精品一区二区三区 | a级毛片黄视频| 国产三级黄色录像| 50天的宝宝边吃奶边哭怎么回事| 热99久久久久精品小说推荐| 欧美另类亚洲清纯唯美| 亚洲少妇的诱惑av| 日本wwww免费看| 一级黄色大片毛片| 电影成人av| 99久久99久久久精品蜜桃| 深夜精品福利| 黄色丝袜av网址大全| tube8黄色片| 国产极品粉嫩免费观看在线| 国产一区二区三区视频了| 欧美在线一区亚洲| 老司机亚洲免费影院| 国产亚洲欧美在线一区二区| 欧美久久黑人一区二区| 91麻豆精品激情在线观看国产 | 日韩一卡2卡3卡4卡2021年| 91老司机精品| 精品一品国产午夜福利视频| 欧美成人免费av一区二区三区 | 欧美变态另类bdsm刘玥| 婷婷丁香在线五月| 丁香六月天网| 国产区一区二久久| 极品少妇高潮喷水抽搐| 波多野结衣一区麻豆| 亚洲欧美精品综合一区二区三区| 亚洲av欧美aⅴ国产| 精品久久久精品久久久| 国产一区二区三区在线臀色熟女 | 国产高清国产精品国产三级| 少妇被粗大的猛进出69影院| 免费av中文字幕在线| 欧美精品高潮呻吟av久久| 国产欧美日韩综合在线一区二区| 欧美日韩亚洲高清精品| 午夜精品久久久久久毛片777| 国产日韩欧美视频二区| 一级片免费观看大全| av不卡在线播放| 天天添夜夜摸| 在线观看舔阴道视频| 国产欧美亚洲国产| 久久青草综合色| 国产91精品成人一区二区三区 | 蜜桃在线观看..| av网站在线播放免费| 国产日韩欧美亚洲二区| 日韩一卡2卡3卡4卡2021年| 久久久久久久大尺度免费视频| 又大又爽又粗| 国产精品亚洲一级av第二区| 亚洲精品自拍成人| 激情在线观看视频在线高清 | 久久精品91无色码中文字幕| 久久99热这里只频精品6学生| 国产极品粉嫩免费观看在线| 大片免费播放器 马上看| 激情视频va一区二区三区| 女人久久www免费人成看片| 一级片'在线观看视频| 国产精品98久久久久久宅男小说| 国产一区有黄有色的免费视频| 男女之事视频高清在线观看| 中文字幕最新亚洲高清| 亚洲国产成人一精品久久久| 国产有黄有色有爽视频| 久久久精品94久久精品| 国产精品久久久久久精品古装| 久久精品人人爽人人爽视色| tocl精华| 国产免费视频播放在线视频| 国产高清视频在线播放一区| 亚洲第一青青草原| 人人妻人人澡人人看| av免费在线观看网站| 国产免费现黄频在线看| www.精华液| 老熟妇仑乱视频hdxx| 欧美在线一区亚洲| 午夜91福利影院| 90打野战视频偷拍视频| 王馨瑶露胸无遮挡在线观看| 黑丝袜美女国产一区| 国产成人免费无遮挡视频| 国产精品二区激情视频| 高潮久久久久久久久久久不卡| 欧美激情 高清一区二区三区| 久久久久精品国产欧美久久久| 亚洲国产看品久久| 午夜久久久在线观看| 国产淫语在线视频| 免费av中文字幕在线| 日韩欧美一区视频在线观看| 在线播放国产精品三级| 国产91精品成人一区二区三区 | 精品国产一区二区三区四区第35| 2018国产大陆天天弄谢| av视频免费观看在线观看| 午夜成年电影在线免费观看| 亚洲视频免费观看视频| 黄片播放在线免费| 久久久水蜜桃国产精品网| 色视频在线一区二区三区| 一本一本久久a久久精品综合妖精| 最新美女视频免费是黄的| 少妇被粗大的猛进出69影院| 精品第一国产精品| 国产一区二区三区综合在线观看| 日日摸夜夜添夜夜添小说| 99国产精品免费福利视频| 日本av手机在线免费观看| www.999成人在线观看| 黑人欧美特级aaaaaa片| 精品少妇黑人巨大在线播放| 91精品国产国语对白视频| 日韩熟女老妇一区二区性免费视频| 免费一级毛片在线播放高清视频 | 成年人午夜在线观看视频| 国产成人精品久久二区二区91| 久久久久网色| 国产福利在线免费观看视频| 亚洲成人手机| 一进一出好大好爽视频| av视频免费观看在线观看| 人人妻人人澡人人爽人人夜夜| 人妻久久中文字幕网| 亚洲专区中文字幕在线| 久久免费观看电影| 国产成人精品久久二区二区免费| kizo精华| 精品久久蜜臀av无| 12—13女人毛片做爰片一| 色精品久久人妻99蜜桃| www日本在线高清视频| av有码第一页| 日本wwww免费看| 在线观看66精品国产| 久久久国产精品麻豆| 午夜成年电影在线免费观看| 国产日韩欧美亚洲二区| 后天国语完整版免费观看| 午夜激情av网站| 国产有黄有色有爽视频| 啦啦啦在线免费观看视频4| 久久久精品免费免费高清| 久久久久久亚洲精品国产蜜桃av| 亚洲国产欧美日韩在线播放| 母亲3免费完整高清在线观看| 咕卡用的链子| 欧美日本中文国产一区发布| 久久久精品94久久精品| 亚洲欧美色中文字幕在线| 国产成人一区二区三区免费视频网站| 国产日韩欧美视频二区| 满18在线观看网站| 黑人操中国人逼视频| 一二三四社区在线视频社区8| 国产精品久久电影中文字幕 | 黄色视频,在线免费观看| 亚洲精品国产一区二区精华液| 国产成人欧美| 国产亚洲欧美在线一区二区| 国产精品欧美亚洲77777| 人人妻,人人澡人人爽秒播| 1024视频免费在线观看| 亚洲av成人一区二区三| 一二三四社区在线视频社区8| 高清av免费在线| 日本五十路高清| 久久久国产成人免费| 在线观看免费视频日本深夜| 老司机午夜十八禁免费视频| 757午夜福利合集在线观看| 久久精品人人爽人人爽视色| 色视频在线一区二区三区| 亚洲精品粉嫩美女一区| 夜夜骑夜夜射夜夜干| 交换朋友夫妻互换小说| 欧美日韩黄片免| 免费看十八禁软件| 国产老妇伦熟女老妇高清| 一个人免费在线观看的高清视频| 大片电影免费在线观看免费| 久久99热这里只频精品6学生| 国产一卡二卡三卡精品| 免费不卡黄色视频| 欧美 日韩 精品 国产| 色精品久久人妻99蜜桃| av不卡在线播放| 丝袜美腿诱惑在线| 精品卡一卡二卡四卡免费| 久久久久精品人妻al黑| 久久影院123| 国产免费av片在线观看野外av| 午夜激情久久久久久久| 国产一区有黄有色的免费视频| 一区二区av电影网| 国产精品免费视频内射| 交换朋友夫妻互换小说| 咕卡用的链子| 人妻一区二区av| 正在播放国产对白刺激| 免费看a级黄色片| 成年动漫av网址| 久9热在线精品视频| 又紧又爽又黄一区二区| 啪啪无遮挡十八禁网站| 美国免费a级毛片| 女人久久www免费人成看片| 国产99久久九九免费精品| 深夜精品福利| 9色porny在线观看| 他把我摸到了高潮在线观看 | 丰满饥渴人妻一区二区三| a级毛片在线看网站| 激情在线观看视频在线高清 | 国产黄频视频在线观看| 一级a爱视频在线免费观看| 亚洲av日韩在线播放| 我要看黄色一级片免费的| 国产男女超爽视频在线观看| 国产免费视频播放在线视频| 国产在线免费精品| av国产精品久久久久影院| 99re6热这里在线精品视频| 久久久久久久久免费视频了| 黄片小视频在线播放| 国产日韩一区二区三区精品不卡| 久久久欧美国产精品| 在线播放国产精品三级| 亚洲国产毛片av蜜桃av| 国产日韩一区二区三区精品不卡| 久久国产精品大桥未久av| 黄色视频在线播放观看不卡| 成年人黄色毛片网站| 国产在线一区二区三区精| 国产亚洲欧美精品永久| 亚洲久久久国产精品| 午夜福利影视在线免费观看| 欧美日韩成人在线一区二区| 女人高潮潮喷娇喘18禁视频| 成人影院久久| 成在线人永久免费视频| 精品少妇一区二区三区视频日本电影| 日本wwww免费看| 日韩欧美一区视频在线观看| 国产1区2区3区精品| 亚洲欧美一区二区三区黑人| 精品一区二区三区视频在线观看免费 | 色在线成人网| 亚洲欧洲日产国产| 在线永久观看黄色视频| 精品第一国产精品| 亚洲熟妇熟女久久| 丝袜美足系列| 亚洲精品自拍成人| 亚洲伊人色综图| 在线播放国产精品三级| 国产一区二区三区在线臀色熟女 | 热99国产精品久久久久久7| 大香蕉久久成人网| 无遮挡黄片免费观看| 欧美久久黑人一区二区| 国产亚洲精品久久久久5区| 成人特级黄色片久久久久久久 | 天天躁夜夜躁狠狠躁躁| 国产精品免费大片| 美女高潮到喷水免费观看| 美女国产高潮福利片在线看| 别揉我奶头~嗯~啊~动态视频| 日韩 欧美 亚洲 中文字幕| 国产99久久九九免费精品| 久久99一区二区三区| 亚洲精品国产区一区二| 精品午夜福利视频在线观看一区 | 久久午夜综合久久蜜桃| 久久中文字幕一级| 久久久久久免费高清国产稀缺| 一本—道久久a久久精品蜜桃钙片| 午夜福利,免费看| 一本大道久久a久久精品| 欧美精品一区二区免费开放| 日本欧美视频一区| 亚洲色图综合在线观看| 欧美精品一区二区大全| www.自偷自拍.com| 正在播放国产对白刺激| 国内毛片毛片毛片毛片毛片| 超碰97精品在线观看| 精品午夜福利视频在线观看一区 | 9色porny在线观看| 女性生殖器流出的白浆| 无限看片的www在线观看| 欧美在线一区亚洲| xxxhd国产人妻xxx| av天堂在线播放| 少妇的丰满在线观看| 日韩欧美一区二区三区在线观看 | 制服诱惑二区| 国产深夜福利视频在线观看| 久久精品国产99精品国产亚洲性色 | 国产精品免费视频内射| 午夜精品久久久久久毛片777| 在线观看人妻少妇| 正在播放国产对白刺激| 亚洲av美国av| 午夜日韩欧美国产| 久久久久久久精品吃奶| 婷婷成人精品国产| 一区二区av电影网| 国产不卡一卡二| 人人澡人人妻人| 老司机影院毛片| 曰老女人黄片| 在线观看免费日韩欧美大片| 嫩草影视91久久| 欧美另类亚洲清纯唯美| 91精品三级在线观看| 窝窝影院91人妻| 亚洲av日韩在线播放| 欧美日韩黄片免| 国产精品亚洲av一区麻豆| 另类亚洲欧美激情| 久久 成人 亚洲| 99国产精品一区二区蜜桃av | 免费看a级黄色片| 一级毛片电影观看| 亚洲精品国产精品久久久不卡| 在线十欧美十亚洲十日本专区| 国产精品偷伦视频观看了| 亚洲 国产 在线| 中文字幕色久视频| 久久久欧美国产精品| 美女福利国产在线| 在线看a的网站| 欧美日韩黄片免| 国产成人欧美在线观看 | 天堂中文最新版在线下载| 又紧又爽又黄一区二区| 日本一区二区免费在线视频| 99re在线观看精品视频| 91老司机精品| 免费观看人在逋| 一边摸一边做爽爽视频免费| 热99re8久久精品国产| avwww免费| 啦啦啦中文免费视频观看日本| 国产av精品麻豆| 黄色成人免费大全| 男人操女人黄网站| 天天躁狠狠躁夜夜躁狠狠躁| 天堂8中文在线网| 欧美日韩亚洲综合一区二区三区_|