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

    利用最小二乘配置識別川滇地區(qū)跨斷層形變異常的研究

    2016-11-17 02:36:16李瑞莎
    大地測量與地球動力學(xué) 2016年11期
    關(guān)鍵詞:測段魯?shù)?/a>蘆山

    賈 鵬 張 希 路 珍 李瑞莎

    1 中國地震局第二監(jiān)測中心,西安市西影路316號,710054

    ?

    利用最小二乘配置識別川滇地區(qū)跨斷層形變異常的研究

    賈 鵬1張 希1路 珍1李瑞莎1

    1 中國地震局第二監(jiān)測中心,西安市西影路316號,710054

    將最小二乘配置方法的高斯型經(jīng)驗協(xié)方差函數(shù)模型應(yīng)用于區(qū)域跨斷層形變曲線的異常輔助識別。通過川滇地區(qū)跨斷層形變對蘆山7.0級、康定6.3級、魯?shù)?.5級等地震的前兆異常判定曲線,說明最小二乘配置內(nèi)插擬合可以進行區(qū)域跨斷層水準形變曲線的趨勢模擬和輔助查找區(qū)間異常特征,對于跨斷層形變曲線異常識別有一定意義。

    最小二乘配置;川滇地區(qū);跨斷層形變;擬合;前兆異常

    川滇地區(qū)位于南北地震帶中南段,發(fā)育有鮮水河斷裂帶、龍門山斷裂帶、安寧河斷裂帶、紅河斷裂帶、小江斷裂帶等,跨斷層形變場地主要布設(shè)在這些大小斷裂上,對整個川滇地區(qū)的斷裂活動有很好的監(jiān)測作用??鐢鄬有巫儽O(jiān)測對強震中短期(即數(shù)月至1 a左右)前兆有明顯顯示,對地震預(yù)測有一定的意義[1]。2008年和2013年在鮮水河斷裂、龍門山斷裂與安寧河斷裂的“Y”型交匯區(qū),接連發(fā)生汶川M8.0和蘆山M7.0強震,2014年發(fā)生康定M6.3地震,并且在川滇交界處還發(fā)生魯?shù)镸6.5地震,分布在川滇地區(qū)斷裂帶上的跨斷層場地測量記錄到了相關(guān)形變數(shù)據(jù)。以往的研究以識別觀測曲線異常,尋找加速、轉(zhuǎn)折、突跳等直觀異常[2]為主,由于實際監(jiān)測時間并不長(只有30~40 a),對所在斷裂帶構(gòu)造活動總體趨勢及其動態(tài)變化,包括突跳、加速、轉(zhuǎn)折等異常指標,僅從觀測數(shù)據(jù)上看也會受到環(huán)境和其他因素的影響,不好直接判定。

    最小二乘配置是根據(jù)已知點信號、協(xié)方差及其與待估算點的協(xié)方差關(guān)系而獲得的待估算點的無偏最優(yōu)估計,綜合了平差、推估和濾波。江在森、張希等[3-5]對該方法進行過深入的研究和探討,將協(xié)方差經(jīng)驗函數(shù)進行了簡化,并且張希等將一維時間域內(nèi)的推估內(nèi)插進行了驗證。武艷強等[6-7]將其用于GPS連續(xù)站資料分析,驗證了反映時序變化特征的可行性。賈鵬[8]等將時域最小二乘配置的擬合和外推也應(yīng)用到定點形變觀測特征曲線尋找及其異常識別。本文將最小二乘配置的擬合應(yīng)用至跨斷層形變測線中,輔助尋找識別震前尤其是蘆山和魯?shù)榈卣鸬那罢桩惓!?/p>

    1 模型論述

    假設(shè)待內(nèi)插區(qū)域有m0個已知點觀測(或計算)值,設(shè)為L=(g1,g2,…,gm0)T,其中每個點值gi的中誤差值為mgi,i=1,2,…,m0。t為要濾波的已知點信號,n為觀測誤差向量,s是待估算點信號,t和n都是中心分布的。那么最小二乘配置的基本方程為:

    (1)

    則:

    (2)

    (3)

    將某一坐標為(x,y)的待估算點的值表示為g,而用c(a,b)表示變量a、b間協(xié)方差,得:

    (4)

    上述各式中,Cst、Ctt均根據(jù)同一個高斯型經(jīng)驗協(xié)方差函數(shù)確定:

    (5)

    為保證f(0)大于0,定義fr(0)=αfL(0)。0<α≤0.2,令α>0意為必須濾波。

    確定參數(shù)是最小二乘配置實現(xiàn)的關(guān)鍵,但一般情況下很難得到可靠性較好的協(xié)方差圖形,故根據(jù)具體地區(qū)測點分布情況來確定參數(shù)k,即首先確定擬合量在整個區(qū)域的相關(guān)距離S(即超出這一距離,則點間協(xié)方差值接近于零),而參數(shù)[1]

    (6)

    設(shè)dij(i,j=1,2,…,m0)為任兩點間的距離,而

    分別定義為最小相鄰點距、平均相鄰點距、最大相鄰點距、最大點距??扇。?/p>

    (7)

    在文獻[8]中筆者已經(jīng)對模型參數(shù)k進行了探討和確定,本文根據(jù)公式(6-7.2), 使用程序自動計算的缺省值為參數(shù)k值。由于篇幅有限,在此僅選取4個川滇跨斷層場地測段(如圖1,2):虛墟、湯家坪、龍燈壩3個基線場地和爾烏水準場地,其中虛墟和龍燈壩場地在鮮水河斷裂上,湯家坪場地在則木河斷裂上,爾烏水準場地在安寧河斷裂上。圖2中斷層箭頭所指為上盤,其他圖或文字中箭頭所指亦為上盤。對于所選的場地以初次測得兩基線L(水準H)點距(高)差為初始值記為0,利用第n期測的高差值與初始高差值之差作為第n期累積形變ΔL(ΔH)為縱軸,橫軸為時間,獲得跨斷層基線L(水準H)點的動態(tài)變化曲線。對于斷層來講,不同性質(zhì)的活動斷層一般在理想狀況下,ΔL 和ΔH隨時間推移基本都是持續(xù)增加或減少,但是考慮到實際因素的影響,比如氣溫、降雨、人為干擾、測量誤差以及地震等,ΔL和ΔH非持續(xù)增加或減少,而是在一定范圍內(nèi)出現(xiàn)波動和折返。排除這些干擾后所獲取的波動和折返的單期或多期變化異常形態(tài)則是我們所需要的。

    圖1 川滇地區(qū)跨斷層形變場地Fig.1 Cross-fault sites in Sichuan-Yunan area

    圖2 選取測跨斷層形變場地簡圖Fig.2 Simple geological map of selected cross-fault site

    2 區(qū)域跨斷層測段實例應(yīng)用

    南北地震帶南段的川滇地區(qū),主要以鮮水河斷裂為主,在2013-04-20蘆山M7.0地震前數(shù)月出現(xiàn)了多個場地異常,包括侏倭、虛墟、龍燈壩等基線場地,以大幅壓縮、壓性變化為主,震后轉(zhuǎn)折呈恢復(fù)減弱跡象。如圖3,跨斷層短基線虛墟(B←A)測段數(shù)據(jù)獲得始于1981年,整體趨勢在2012年前主要以拉張上行為主,這期間有小幅上下波動變化。從整體上看,在2013年蘆山地震前一年有比較明顯的大幅快速壓縮變化,異常較為明顯。在2013年9月前后,該基線再次連續(xù)壓縮,在恢復(fù)過程中再次發(fā)生了康定M6.3地震。

    圖3 跨斷層基線場地虛墟(B←A)測段變化曲線Fig.3 Variation curve of cross-fault baseline site Xuxu(B←A)

    圖4 跨斷層基線場地虛墟(B←A)測段原曲線(上圖黑色)與擬合值對比(上圖紅色)、絕對值差異曲線(下圖黑色)及二倍均方差(下圖藍色虛線)Fig.4 Original curve of baseline in Xuxu (top black) compared with fitting curve (top red), absolute value curve of the difference (bottom solid) with the double mean square error (bottom blue dashed)

    由圖4,通過最小二乘配置的擬合計算,該場地B←A測段30多年的趨勢特征被很好地擬合出來。擬合趨勢與原測段觀測曲線相比,實際觀測曲線2013年快速壓縮明顯;而通過絕對值與二倍均方差殘差可以看出,蘆山地震前的曲線趨勢異常與擬合值差值很大,可以確定為異常變化;在2013年年底前,擬合值殘差部分同樣小幅超過了二倍均方差,這部分異常變化同樣能反映出2014年康定地震的前期應(yīng)變積累異常。

    圖5 跨斷層短基線場地龍燈壩(D←A)測段變化曲線Fig.5 Variation curve of cross-fault baseline site Longdengba(D←A)

    龍燈壩場地(D←A)基線與虛墟場地相隔不遠,屬鮮水河斷裂中北段,從1984 年對龍燈壩場地觀測開始,至今已連續(xù)觀測了31a,累計觀測數(shù)據(jù)240多期。該基線也在蘆山地震前一年出現(xiàn)了大幅壓縮變化,不同的是該測線整體年變幅度波動較大。2000年、2003年以及蘆山地震恢復(fù)后也出現(xiàn)過較大的壓縮變化,該場地距離康定地震較近。但僅從原始曲線看,除了蘆山地震前壓縮極為明顯,2000年、2003年以及2013年末的異常變化都是存在的,我們利用最小二乘配置進行擬合看是否對應(yīng)為異常情況。

    圖6 跨斷層基線場地龍燈壩(D←A)測段原曲線(上圖黑色)與擬合值對比(上圖紅色)、絕對值差異曲線(下圖黑色)及二倍均方差(下圖藍色虛線)Fig.6 Original curve of baseline in Longdengba(top black) compared with fitting curve (top red), absolute value curve of the difference (bottom solid) with the double mean square error (bottom blue dashed)

    圖6紅色曲線即為最小二乘配置擬合的結(jié)果,整體趨勢平穩(wěn),2004年趨勢由平轉(zhuǎn)折拉張,2013年蘆山地震后的趨勢上行加速,在蘆山地震前出現(xiàn)明顯的趨勢差異,且異常明顯,擬合值和原始值之差也超過二倍均方差幅度。2002年、2006年有小段變化差異超過二倍均方差,但是幅度非常?。?009年小幅度的超出二倍均方差的異常存在,不能排除是汶川地震后調(diào)整的影響。此外,2013年末的一段超出二倍均方差的異常,即在蘆山地震后,基線恢復(fù)后的轉(zhuǎn)折壓縮,時間段上比2009年長,應(yīng)為康定地震前的應(yīng)變積累所致。

    圖7 跨斷層短基線場地湯家坪(B←A)測段變化曲線Fig.7 Variation curve of cross-fault short baseline site Tangjiaping(B←A)

    位于則木河斷裂上的湯家坪與鮮水河斷裂的上面兩處場地不同,在2010年或2011年以來持續(xù)正斷、拉張加速。2013年后再次拉張,2014年轉(zhuǎn)平下行,魯?shù)榈卣鸷筠D(zhuǎn)折。單看測段原基線整體趨勢,變化波動較為繁雜。

    利用最小二乘配置擬合后(如圖8紅色曲線),顯示出湯家坪場地基線變化趨勢。2010年該場地基線以壓縮為主,而2010年到2014年魯?shù)榈卣鹎俺霈F(xiàn)了大幅的拉張變化,魯?shù)榈卣鸷筅厔蒉D(zhuǎn)折。從二倍均方差差異曲線來看,2010~2012年比較集中。而2005年年初也出現(xiàn)了比較大的差異,不能排除與2004-12-26印尼M8.9強震的影響有一定關(guān)系[9]。

    圖8 跨斷層基線場地湯家坪(B←A)測段原曲線(上圖黑色)與擬合值對比(上圖紅色)、絕對值差異曲線(下圖黑色)及二倍均方差(下圖藍色虛線)Fig.8 Original curve of baseline in Tangjiaping (top black) compared with fitting curve (top red), absolute value curve of the difference (bottom solid) with the double mean square error(bottom blue dashed)

    圖9 跨斷層短水準場地爾烏(3←1)測段變化曲線Fig.9 Variation curve of cross-fault short leveling site Erwu(3←1)

    然而,對于跨斷層形變來說,測段曲線異常變化復(fù)雜,例如位于安寧河與則木河斷裂交匯中段的爾烏水準場地(3←1)測段,在汶川地震前后有拉張后的壓縮變化,在蘆山地震前也有一段時期的快速拉張變化。我們利用最小二乘配置擬合和殘差看是否在震前與震后有比較明顯的趨勢差異。

    圖10 跨斷層短水準場地爾烏(3←1)測段原曲線(上圖黑色)與擬合值對比(上圖紅色)、絕對值差異曲線(下圖黑色)及二倍均方差(下圖藍色虛線)Fig.10 Original curve of leveling in Erwu (top black) compared with fitting curve (top red), absolute value curve of the difference (bottom solid) with the double mean square error (bottom blue dashed)

    通過最小二乘配置擬合該測段的整體趨勢可以看出,在汶川地震后,2010年開始整個爾烏水準(3←1)測段是一個拉張的過程,有一個較為明顯的趨勢性異常,原觀測曲線與擬合的二倍均方差差異比較大的也是2010~2012年這一段。另外,爾烏水準場地距離蘆山地震和魯?shù)榈卣鸬恼鹬芯喽疾凰憬J山和魯?shù)榈卣鸷蟛⑽椿謴?fù)到2010年前的原有趨勢,鑒于截至本稿前該場地的異常趨勢依然存在,筆者認為不能排除其對應(yīng)其他前兆異常信號,是否能夠進一步確定是其他地震的前兆異常,要根據(jù)未來的監(jiān)測數(shù)據(jù)和實際場地情況甚至結(jié)合其他測震手段綜合而定。筆者也對川滇地區(qū)的其他場地進行了計算擬合,總體效果良好,限于篇幅,不再贅述。因此,利用最小二乘配置擬合,能夠檢驗跨斷層形變偏離已有數(shù)學(xué)模型所能描述的部分信息,對輔助判定區(qū)域跨斷層形變前兆異常有一定的意義。

    3 結(jié) 語

    筆者嘗試利用最小二乘配置進行區(qū)域跨斷層形變曲線特征模擬和尋找異常特征,通過內(nèi)插區(qū)域內(nèi)所有已知點的相關(guān)性,反映其隨時間變化的趨勢性。川滇地區(qū)斷層形變復(fù)雜,同一條斷層不同地震的孕震、發(fā)震機理甚至遠距離超強大震的影響,也給前兆異常查找和識別帶來了很大困難,很多時候還是一種經(jīng)驗和探索。例如蘆山地震前鮮水河北段短期異常十分突出[10],而地震卻發(fā)生在高應(yīng)力積累的形變閉鎖區(qū)的龍門山斷裂帶;同樣是龍門山斷裂帶上相近的汶川和蘆山兩個大震,前者沒有明顯異常前兆而后者卻有。通過濾波和擬合,可以較好地反映曲線周期,從而根據(jù)該場地的周期和短時變化偏離數(shù)學(xué)模型的部分輔助識別異常。

    [1] 江在森,丁平,王雙緒,等.中國西部大地形變監(jiān)測與地震預(yù)測[M].北京:地震出版社,2001(JiangZaisen,DingPing,WangShuangxu,etal.GeodeticDeformationMonitoringandEarthquakePredictioninWestChina[M].Beijing:SeismologicalPress, 2001)

    [2] 楊曉東,李寧,劉立煒,等. 2013年甘肅岷縣漳縣Ms6.6地震跨斷層形變異常研究[J].地震研究,2014,37(4):578-587(YangXiaodong,LiNing,LiuLiwei,etal.ResearchonCross-faultDeformationAbnormityofMinxian-ZhangxianMs6.6EarthquakeinGansuin2013. [J].JournalofSeismologicalResearch,2014,37(4):578-587)

    [3] 張希,江在森,張四新.借助最小二乘配置整體解算地殼視應(yīng)變場[J].地殼形變與地震,1998,18(2):57-62(ZhangXi,JiangZaisen,ZhangSixin,etal.WholeCalculationofTheCrustalVisualStrainFieldwithTheLeastSquareCollocation[J].CrustalDeformationandEarthquake,1998,18(2):57-62)

    [4] 張希,江在森.用最小二乘配置獲得地形變應(yīng)變場動態(tài)圖像的幾個問題研究[J].地殼形變與地震,1999,19(3):32-39(ZhangXi,JiangZaisen.StudyonSomeQuestionsofDynamicPicturesofCrustalDeformationandStrainFieldsObtainedbytheLeastSquareCollocation[J].CrustalDeformationandEarthquake,1999,19(3):32-39)

    [5] 江在森,張希.華北地區(qū)近期地殼水平運動與應(yīng)力應(yīng)變場特征[J].地球物理學(xué)報,2000,43(5):657-665(JiangZaisen,ZhangXi.CharacteristicsofRecentHorizontalMovementandStrain-StressFiledintheCrustofNorthChina[J].ChineseJournalofGeophysics,2000,43(5):657-665)

    [6] 武艷強,江在森,楊國華.最小二乘配置方法在提取GPS時間序列信息中的應(yīng)用[J].國際地震動態(tài),2007(7):99-103(WuYanqiang,JiangZaisen,YangGuohua.TheApplicationofLeastSquareCollocationinObtainingInformationfromGPSTimeSeries[J].RecentDevelopmentsinWorldSeismology,2007(7):99-103)

    [7] 武艷強,黃立人.時間序列處理的新插值方法[J].大地測量與地球動力學(xué),2004,24(4):43-47(WuYanqiang,HuangLiren.ANewInterpolationMethodinTimeSeriesAnalyzing[J].JournalofGeodesyandGeodynamics,2004,24(4):43-47)

    [8] 賈鵬,張希,唐紅濤,等.時域最小二乘配置求定點形變特征曲線與異常識別的應(yīng)用研究[J].地震研究,2013,36(1):57-62(JiaPeng,Zhangxi,TangHongtao,etal.TheapplicationofLeastSquareCollocationinFittingPointingDeformationCurvesandIdentifyingAnomalyCharacteristics[J].JournalofSeismologicalResearch,2013,36(1):57-62)

    [9] 薛富平,王雙緒.南北地震帶斷層形變與汶川8.0級地震研究[J].國際地震動態(tài),2009(4):40-41(XueFuping,WangShuangxu.ResearchonFaultsDeformationinNorth-SouthSeismicBeltandWenchuanM8.0Earthquake[J].RecentDevelopmentsinWorldSeismology,2009(4):40-41)

    [10]方穎,張晶,江在森,等.用跨斷層形變資料分析鮮水河斷裂西北段的運動特征[J].地球物理學(xué)報,2015,58(5):1 645-1 653(FangYing,ZhangJing,JiangZaisen,etal.MovementCharacteristicsofThnorthwestSegmentoftheXianshuiheFaultZoneDerivedfromCross-FaultDeformationData[J].ChineseJournalofGeophysics,2015,58(5):1 645-1 653)

    About the first author:JIA Peng, engineer,majors in crustal deformation and earthquake prediction, E-mail:45097027@qq.com.

    The Application of Least Square Collocation in Identifying Anomaly Deformation Curves of Cross-Fault in Sichuan-Yunnan Region

    JIAPeng1ZHANGXi1LUZhen1LIRuisha1

    1 Crustal Monitoring and Application Center,CEA,316 Xiying Road,Xi’an 710054,China

    The Gauss-function model, based on least square collocation, is used in anomaly assisted identification of regional cross-fault deformation characteristic curves. The precursory anomaly judgment examples of regional deformation of cross-fault in Lushan Ms7.0, Kangding Ms6.3, Ludian Ms6.5 and other earthquakes, show that least square collocation can be fitted in regional cross-fault deformation and helps in identifying anomaly characteristics. The method has certain significance for the deformation of cross-fault identification.

    least square collocation;Sichuan-Yunnan region;deformation of cross-fault;fitting;anomaly

    The Spark Program of Earthquake Sciences of CEA,No.XH15065Y;Basic Research Program for Natural Science of Shaanxi Province,No.2014JM2-4039;Seismic Regime Tracking Project of CEA,No.2016010205.

    2015-12-15

    項目來源:中國地震局地震科技星火計劃(XH15065Y);陜西省自然科學(xué)基礎(chǔ)研究計劃(2014JM2-4039);中國地震局震情跟蹤定向工作任務(wù)(2016010205)。

    賈鵬,工程師,主要從事地殼形變與地震預(yù)報研究,E-mail:45097027@qq.com。

    10.14075/j.jgg.2016.11.007

    1671-5942(2016)011-0972-05

    P315

    A

    猜你喜歡
    測段魯?shù)?/a>蘆山
    靈武跨斷層水準測量及相關(guān)研究
    西部高原地區(qū)流域水準網(wǎng)施加重力異常改正的必要性
    DINI03電子水準儀數(shù)據(jù)預(yù)處理程序?qū)崿F(xiàn)
    地礦測繪(2016年4期)2017-01-05 01:31:29
    高鐵線路水準基點網(wǎng)復(fù)測穩(wěn)定性分析方法研究
    讓我輕輕地把魯?shù)閾肀Вü?jié)選)
    春回蘆山
    目擊
    蘆山地震公路地質(zhì)災(zāi)害調(diào)查及評估
    從魯?shù)榈卣鹂赐话l(fā)事件應(yīng)急管理機制的進步
    大壩與安全(2014年4期)2014-02-28 02:37:04
    2014年云南魯?shù)?.5級地震烈度分布與房屋震害特征
    地震研究(2014年4期)2014-02-27 09:31:09
    91麻豆av在线| 国产精品野战在线观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲电影在线观看av| 午夜免费成人在线视频| 亚洲精品色激情综合| 99久久无色码亚洲精品果冻| aaaaa片日本免费| 中文字幕最新亚洲高清| 免费一级毛片在线播放高清视频| 19禁男女啪啪无遮挡网站| 亚洲avbb在线观看| 色尼玛亚洲综合影院| 黄色a级毛片大全视频| 又黄又粗又硬又大视频| 精品一区二区三区视频在线观看免费| 热99re8久久精品国产| 免费看十八禁软件| 亚洲精品中文字幕一二三四区| 国产高清有码在线观看视频 | 91在线观看av| 在线国产一区二区在线| 可以在线观看毛片的网站| 丰满人妻熟妇乱又伦精品不卡| 久久久国产成人精品二区| 成人亚洲精品av一区二区| 可以在线观看的亚洲视频| 99riav亚洲国产免费| 国产精华一区二区三区| bbb黄色大片| 精品日产1卡2卡| 精品日产1卡2卡| 日本一本二区三区精品| 国产成人精品久久二区二区免费| 色婷婷久久久亚洲欧美| 动漫黄色视频在线观看| 久久精品人妻少妇| 亚洲精品久久成人aⅴ小说| 在线永久观看黄色视频| 午夜两性在线视频| 中出人妻视频一区二区| 非洲黑人性xxxx精品又粗又长| 免费在线观看黄色视频的| 日本撒尿小便嘘嘘汇集6| 熟女少妇亚洲综合色aaa.| 精品久久久久久久久久免费视频| 午夜日韩欧美国产| 成人国产综合亚洲| 国产乱人伦免费视频| 亚洲国产欧美日韩在线播放| 在线视频色国产色| 国产成年人精品一区二区| 极品教师在线免费播放| 久久亚洲真实| 一级作爱视频免费观看| 精品日产1卡2卡| 啦啦啦观看免费观看视频高清| 国产亚洲精品综合一区在线观看 | 国产亚洲精品av在线| 欧美亚洲日本最大视频资源| 女生性感内裤真人,穿戴方法视频| 国产熟女午夜一区二区三区| 一级毛片精品| 人人妻人人澡欧美一区二区| 国产成人精品久久二区二区91| 特大巨黑吊av在线直播 | 大香蕉久久成人网| 中文字幕最新亚洲高清| 18禁观看日本| 亚洲九九香蕉| 亚洲国产精品久久男人天堂| 成年版毛片免费区| 少妇被粗大的猛进出69影院| 这个男人来自地球电影免费观看| 免费在线观看日本一区| 999久久久国产精品视频| 琪琪午夜伦伦电影理论片6080| 人人妻人人澡人人看| 亚洲中文日韩欧美视频| 日本五十路高清| 成在线人永久免费视频| 亚洲精品国产一区二区精华液| 亚洲av五月六月丁香网| 婷婷六月久久综合丁香| 男人操女人黄网站| 国产精品野战在线观看| 啦啦啦免费观看视频1| 国产成人精品久久二区二区91| 国产成人啪精品午夜网站| 国语自产精品视频在线第100页| 久久精品国产亚洲av香蕉五月| 中亚洲国语对白在线视频| 正在播放国产对白刺激| 国产99白浆流出| 老汉色av国产亚洲站长工具| bbb黄色大片| 91九色精品人成在线观看| 波多野结衣av一区二区av| 搡老妇女老女人老熟妇| 制服人妻中文乱码| 一进一出抽搐gif免费好疼| 人妻久久中文字幕网| 美国免费a级毛片| 在线视频色国产色| 在线观看免费日韩欧美大片| 黄色毛片三级朝国网站| 久久国产精品男人的天堂亚洲| 精品免费久久久久久久清纯| 精品第一国产精品| 亚洲av成人不卡在线观看播放网| 大型av网站在线播放| 亚洲第一电影网av| 国内少妇人妻偷人精品xxx网站 | av在线播放免费不卡| 18美女黄网站色大片免费观看| 国产人伦9x9x在线观看| 久久精品人妻少妇| 男女视频在线观看网站免费 | 国产一级毛片七仙女欲春2 | 中文亚洲av片在线观看爽| 国产精品二区激情视频| 国产97色在线日韩免费| 久久国产精品影院| 久久香蕉国产精品| 亚洲成国产人片在线观看| 激情在线观看视频在线高清| 国产av不卡久久| 中文字幕精品亚洲无线码一区 | 亚洲中文av在线| 日韩一卡2卡3卡4卡2021年| 黑人操中国人逼视频| 成人三级黄色视频| 给我免费播放毛片高清在线观看| 亚洲专区字幕在线| 18禁国产床啪视频网站| 成年免费大片在线观看| 熟女电影av网| 91成年电影在线观看| 少妇的丰满在线观看| 国产精品98久久久久久宅男小说| 亚洲男人天堂网一区| 欧美乱妇无乱码| 夜夜躁狠狠躁天天躁| 色老头精品视频在线观看| 日本黄色视频三级网站网址| 欧美一区二区精品小视频在线| 草草在线视频免费看| 久久中文字幕人妻熟女| 午夜福利视频1000在线观看| 午夜福利成人在线免费观看| 亚洲国产欧美网| 91字幕亚洲| 国产av不卡久久| 欧美成狂野欧美在线观看| 50天的宝宝边吃奶边哭怎么回事| 精品国产超薄肉色丝袜足j| 久久精品国产亚洲av香蕉五月| 日本免费a在线| 很黄的视频免费| 欧美一级a爱片免费观看看 | 免费在线观看影片大全网站| 一进一出好大好爽视频| 亚洲精华国产精华精| 国产激情久久老熟女| 欧美国产精品va在线观看不卡| 91成人精品电影| 国产午夜精品久久久久久| 叶爱在线成人免费视频播放| 美国免费a级毛片| 久久香蕉精品热| www.自偷自拍.com| 午夜免费观看网址| 久久亚洲真实| 中出人妻视频一区二区| 欧美乱码精品一区二区三区| 欧美一区二区精品小视频在线| 天堂动漫精品| 国产精华一区二区三区| 国产av一区在线观看免费| 欧美激情久久久久久爽电影| 丁香欧美五月| 久久中文字幕人妻熟女| 欧美丝袜亚洲另类 | 男女午夜视频在线观看| 国产精品自产拍在线观看55亚洲| 国产视频内射| 亚洲一区二区三区不卡视频| 12—13女人毛片做爰片一| 亚洲成av片中文字幕在线观看| 最近最新免费中文字幕在线| 欧美色欧美亚洲另类二区| 男女床上黄色一级片免费看| 久久精品91无色码中文字幕| 一个人观看的视频www高清免费观看 | 狠狠狠狠99中文字幕| 很黄的视频免费| 欧洲精品卡2卡3卡4卡5卡区| 美国免费a级毛片| 亚洲熟女毛片儿| 99re在线观看精品视频| 久久精品国产亚洲av高清一级| 国产单亲对白刺激| 超碰成人久久| 国产又爽黄色视频| 老鸭窝网址在线观看| 精品久久久久久成人av| 99国产精品一区二区三区| 国产男靠女视频免费网站| 成人亚洲精品一区在线观看| x7x7x7水蜜桃| 99久久99久久久精品蜜桃| 午夜福利在线在线| 久久久精品欧美日韩精品| 国产成人一区二区三区免费视频网站| 热re99久久国产66热| 麻豆久久精品国产亚洲av| 欧美黑人欧美精品刺激| 成人三级做爰电影| 国产午夜福利久久久久久| 日韩欧美免费精品| 一级毛片高清免费大全| 中文字幕人妻丝袜一区二区| 精品欧美国产一区二区三| 精品第一国产精品| 999精品在线视频| 亚洲精品美女久久av网站| 美女 人体艺术 gogo| 国产精品99久久99久久久不卡| 啦啦啦 在线观看视频| 欧美日韩亚洲国产一区二区在线观看| 午夜免费成人在线视频| 国产真实乱freesex| 又黄又爽又免费观看的视频| 国产真人三级小视频在线观看| 18禁美女被吸乳视频| 国产成人av激情在线播放| 波多野结衣av一区二区av| 久久午夜亚洲精品久久| 亚洲一区中文字幕在线| 精品国产乱码久久久久久男人| 人人澡人人妻人| 中文字幕人成人乱码亚洲影| 久久久久久亚洲精品国产蜜桃av| 国产99白浆流出| 大型av网站在线播放| 99国产精品99久久久久| 欧美黑人欧美精品刺激| 欧美成人一区二区免费高清观看 | 日本 欧美在线| 久久久久久久午夜电影| 欧美激情 高清一区二区三区| 国产成年人精品一区二区| 国产精品九九99| 99精品欧美一区二区三区四区| 欧美激情极品国产一区二区三区| 一个人免费在线观看的高清视频| 久久久国产精品麻豆| www.www免费av| 香蕉国产在线看| 国产亚洲欧美98| 亚洲专区字幕在线| 美女午夜性视频免费| 亚洲aⅴ乱码一区二区在线播放 | 久久国产精品人妻蜜桃| 在线观看www视频免费| 波多野结衣高清作品| 欧美色视频一区免费| 国产1区2区3区精品| 伊人久久大香线蕉亚洲五| 成在线人永久免费视频| 久久久久国产一级毛片高清牌| 国产精品99久久99久久久不卡| 亚洲精品美女久久久久99蜜臀| 午夜老司机福利片| 国产免费男女视频| avwww免费| 国产精品久久视频播放| 国产成人欧美在线观看| 亚洲黑人精品在线| 亚洲欧美精品综合久久99| 窝窝影院91人妻| 丝袜在线中文字幕| 欧美日韩中文字幕国产精品一区二区三区| 在线观看免费午夜福利视频| 好看av亚洲va欧美ⅴa在| 亚洲第一av免费看| 在线观看午夜福利视频| 欧美一级a爱片免费观看看 | 亚洲一区二区三区色噜噜| 一个人免费在线观看的高清视频| 夜夜夜夜夜久久久久| 国产精品国产高清国产av| 黄片小视频在线播放| 性色av乱码一区二区三区2| 国产成人欧美在线观看| 黄色成人免费大全| 国产成人精品无人区| 欧美日本视频| 亚洲人成77777在线视频| 人妻丰满熟妇av一区二区三区| 久久精品国产亚洲av香蕉五月| 亚洲第一欧美日韩一区二区三区| 精品久久久久久,| 无人区码免费观看不卡| 欧美+亚洲+日韩+国产| 真人一进一出gif抽搐免费| 女同久久另类99精品国产91| 丝袜人妻中文字幕| 欧美精品亚洲一区二区| 中亚洲国语对白在线视频| 欧美国产精品va在线观看不卡| 黑丝袜美女国产一区| 久久国产亚洲av麻豆专区| 亚洲国产欧洲综合997久久, | 国产亚洲欧美在线一区二区| 啦啦啦韩国在线观看视频| 亚洲第一青青草原| 久久国产精品人妻蜜桃| 88av欧美| 高潮久久久久久久久久久不卡| 在线永久观看黄色视频| 中文字幕最新亚洲高清| 国产精品免费一区二区三区在线| 中文在线观看免费www的网站 | 1024视频免费在线观看| 婷婷精品国产亚洲av| 成人午夜高清在线视频 | 一区二区三区精品91| 亚洲精品一卡2卡三卡4卡5卡| 最近最新中文字幕大全免费视频| 国产免费男女视频| 久久精品国产亚洲av高清一级| 国产精品亚洲美女久久久| videosex国产| 人人妻人人看人人澡| 久久久久久国产a免费观看| 欧美一级a爱片免费观看看 | 午夜激情福利司机影院| 神马国产精品三级电影在线观看 | 色综合婷婷激情| 国产一区二区在线av高清观看| 黄片播放在线免费| www.999成人在线观看| 欧美大码av| 精品国产美女av久久久久小说| 美女扒开内裤让男人捅视频| 久久婷婷人人爽人人干人人爱| xxxwww97欧美| 久久国产乱子伦精品免费另类| av片东京热男人的天堂| 久久欧美精品欧美久久欧美| 99久久精品国产亚洲精品| 日韩精品青青久久久久久| bbb黄色大片| 91字幕亚洲| 国产极品粉嫩免费观看在线| 此物有八面人人有两片| 免费在线观看亚洲国产| 日韩成人在线观看一区二区三区| 色综合欧美亚洲国产小说| 18禁黄网站禁片免费观看直播| 成年女人毛片免费观看观看9| 国产亚洲精品第一综合不卡| 亚洲欧美精品综合久久99| 美女免费视频网站| 国产精品久久电影中文字幕| 女同久久另类99精品国产91| 午夜免费激情av| 国产精品亚洲一级av第二区| 日韩大码丰满熟妇| 国产一级毛片七仙女欲春2 | 99re在线观看精品视频| √禁漫天堂资源中文www| 999久久久国产精品视频| 日本a在线网址| 亚洲第一电影网av| 99国产极品粉嫩在线观看| 一本久久中文字幕| 男女视频在线观看网站免费 | 亚洲av日韩精品久久久久久密| 成人三级做爰电影| 999久久久精品免费观看国产| 日日干狠狠操夜夜爽| 99国产精品99久久久久| 亚洲男人天堂网一区| 国产片内射在线| 色在线成人网| 欧美日韩一级在线毛片| av免费在线观看网站| 亚洲欧美一区二区三区黑人| 人妻久久中文字幕网| 2021天堂中文幕一二区在线观 | 伊人久久大香线蕉亚洲五| 18美女黄网站色大片免费观看| 少妇熟女aⅴ在线视频| 国产一区二区激情短视频| 久久国产精品人妻蜜桃| 国产欧美日韩一区二区三| 99久久综合精品五月天人人| 婷婷精品国产亚洲av在线| av福利片在线| 欧美成人午夜精品| cao死你这个sao货| 中文字幕精品免费在线观看视频| 亚洲人成77777在线视频| 国产精品野战在线观看| 美女高潮喷水抽搐中文字幕| 精品欧美一区二区三区在线| 亚洲精品一区av在线观看| 成年人黄色毛片网站| 超碰成人久久| 嫩草影院精品99| 午夜日韩欧美国产| 香蕉av资源在线| 国产单亲对白刺激| 欧美亚洲日本最大视频资源| 亚洲精华国产精华精| 国产精品亚洲美女久久久| 成人亚洲精品一区在线观看| 欧美最黄视频在线播放免费| 亚洲成av片中文字幕在线观看| 99热这里只有精品一区 | 很黄的视频免费| 国产高清视频在线播放一区| 午夜福利视频1000在线观看| 变态另类成人亚洲欧美熟女| 亚洲av熟女| 又黄又粗又硬又大视频| 欧美激情高清一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 看片在线看免费视频| 波多野结衣高清无吗| 欧美日韩福利视频一区二区| 在线观看66精品国产| 免费人成视频x8x8入口观看| 99国产精品一区二区三区| 听说在线观看完整版免费高清| 我的亚洲天堂| 欧美激情 高清一区二区三区| 国产视频内射| 亚洲电影在线观看av| 亚洲黑人精品在线| 麻豆一二三区av精品| 精品国产美女av久久久久小说| 欧美另类亚洲清纯唯美| 国产av在哪里看| 欧美日韩福利视频一区二区| 老司机靠b影院| 天天添夜夜摸| 久久精品国产综合久久久| 男女床上黄色一级片免费看| 久久国产精品人妻蜜桃| 国产精品,欧美在线| 人妻久久中文字幕网| 91成人精品电影| 精品无人区乱码1区二区| 男人舔奶头视频| 无人区码免费观看不卡| 欧美+亚洲+日韩+国产| 香蕉国产在线看| 成人特级黄色片久久久久久久| 51午夜福利影视在线观看| 9191精品国产免费久久| 国产精品爽爽va在线观看网站 | 波多野结衣高清无吗| 非洲黑人性xxxx精品又粗又长| 国产高清视频在线播放一区| 日韩有码中文字幕| 动漫黄色视频在线观看| 给我免费播放毛片高清在线观看| 在线观看www视频免费| 久久九九热精品免费| 国产午夜福利久久久久久| 嫁个100分男人电影在线观看| 亚洲国产精品合色在线| 日本a在线网址| 成人国产综合亚洲| 亚洲黑人精品在线| 变态另类成人亚洲欧美熟女| 亚洲精品久久国产高清桃花| 黑人巨大精品欧美一区二区mp4| 欧美成人一区二区免费高清观看 | 亚洲一卡2卡3卡4卡5卡精品中文| 老汉色∧v一级毛片| 国产av又大| 日本黄色视频三级网站网址| 久久精品91无色码中文字幕| 国产黄片美女视频| 满18在线观看网站| 一本综合久久免费| 俄罗斯特黄特色一大片| 亚洲av片天天在线观看| 久久香蕉国产精品| 怎么达到女性高潮| 日本精品一区二区三区蜜桃| 久久国产精品人妻蜜桃| 99热只有精品国产| 啦啦啦观看免费观看视频高清| 亚洲精品色激情综合| 精品一区二区三区av网在线观看| 美女免费视频网站| 人妻久久中文字幕网| 美女午夜性视频免费| 亚洲avbb在线观看| 欧美黑人精品巨大| 久久久久精品国产欧美久久久| 久久亚洲精品不卡| 中文字幕最新亚洲高清| 亚洲自拍偷在线| 欧美一级毛片孕妇| 久久婷婷人人爽人人干人人爱| 亚洲欧洲精品一区二区精品久久久| 亚洲成av片中文字幕在线观看| 日日爽夜夜爽网站| 精品久久久久久久毛片微露脸| 男男h啪啪无遮挡| 免费高清在线观看日韩| 成人精品一区二区免费| 免费高清视频大片| 午夜久久久在线观看| 99re在线观看精品视频| 欧美一级毛片孕妇| 男人舔奶头视频| 99久久国产精品久久久| 在线观看一区二区三区| 色av中文字幕| 一区二区日韩欧美中文字幕| 亚洲第一电影网av| 久久天堂一区二区三区四区| 制服人妻中文乱码| 亚洲男人天堂网一区| 高潮久久久久久久久久久不卡| АⅤ资源中文在线天堂| 国产激情欧美一区二区| 国产亚洲欧美在线一区二区| 午夜日韩欧美国产| 极品教师在线免费播放| 欧美乱色亚洲激情| 最近最新中文字幕大全免费视频| 精品福利观看| 在线观看午夜福利视频| 久久久久久免费高清国产稀缺| 后天国语完整版免费观看| 久久婷婷成人综合色麻豆| 最近在线观看免费完整版| 12—13女人毛片做爰片一| 中文字幕精品亚洲无线码一区 | 久久这里只有精品19| 他把我摸到了高潮在线观看| 91九色精品人成在线观看| 男人的好看免费观看在线视频 | а√天堂www在线а√下载| 999久久久精品免费观看国产| 久久草成人影院| 亚洲精品美女久久av网站| 激情在线观看视频在线高清| 99久久国产精品久久久| 国产欧美日韩一区二区三| 搡老熟女国产l中国老女人| 一区二区三区精品91| 欧美av亚洲av综合av国产av| 免费在线观看日本一区| 久久亚洲真实| 国产av在哪里看| 国产精品一区二区免费欧美| 亚洲七黄色美女视频| 99久久久亚洲精品蜜臀av| 欧洲精品卡2卡3卡4卡5卡区| a级毛片a级免费在线| 亚洲三区欧美一区| 亚洲第一欧美日韩一区二区三区| 国产午夜精品久久久久久| 黄色成人免费大全| 最新在线观看一区二区三区| 亚洲精品美女久久久久99蜜臀| 激情在线观看视频在线高清| 精品一区二区三区四区五区乱码| 香蕉av资源在线| 超碰成人久久| 国产精品美女特级片免费视频播放器 | 国产色视频综合| 久久久久久人人人人人| e午夜精品久久久久久久| avwww免费| 夜夜躁狠狠躁天天躁| 国产精品爽爽va在线观看网站 | 久久国产精品男人的天堂亚洲| 国产成人啪精品午夜网站| 亚洲男人天堂网一区| 精品高清国产在线一区| 日韩欧美免费精品| 黄色成人免费大全| 久久久国产成人免费| 91成人精品电影| 日本在线视频免费播放| 天堂√8在线中文| 亚洲国产日韩欧美精品在线观看 | 欧美丝袜亚洲另类 | 国产成人av教育| 91成人精品电影| 久久久久精品国产欧美久久久| 免费观看人在逋| 国产爱豆传媒在线观看 | 妹子高潮喷水视频| 亚洲精品国产一区二区精华液| 亚洲人成77777在线视频| 亚洲中文字幕一区二区三区有码在线看 | 两个人视频免费观看高清| 嫩草影视91久久| 精品一区二区三区av网在线观看| 99久久无色码亚洲精品果冻| 国产亚洲精品第一综合不卡| 窝窝影院91人妻| 嫁个100分男人电影在线观看| 久久狼人影院|