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

    基于CAP方法的震源機制研究

    2010-09-05 12:57:33趙凌云陳俊華申學(xué)林魏貴春
    長江科學(xué)院院報 2010年5期
    關(guān)鍵詞:秭歸臺站震源

    趙凌云,鄧 津,陳俊華,申學(xué)林,戴 苗,魏貴春

    (1.湖北省地震局,武漢 430071;2.甘肅省地震局,蘭州 730000)

    基于CAP方法的震源機制研究

    趙凌云1,鄧 津2,陳俊華1,申學(xué)林1,戴 苗1,魏貴春1

    (1.湖北省地震局,武漢 430071;2.甘肅省地震局,蘭州 730000)

    利用來源于湖北區(qū)域臺網(wǎng)共6個臺的寬頻帶數(shù)字地震記錄,采用CAP法(Cut and Paste Method)反演了2008年11月22日湖北 秭歸Ms4.1級地震和09月27日Ms3.2中小震的震源機制解。結(jié)果顯示,2次地震的震源機制解有一些差別,但是均為逆沖兼走滑分量性質(zhì)斷層,主壓應(yīng)力P軸主要集中在NE向,與區(qū)域主壓應(yīng)力方向有一定一致性。分析認(rèn)為秭歸地震可能是仙女山斷裂受庫水滲透和載荷變化影響引發(fā)的水庫構(gòu)造地震。通過與P波初動地震震源機制結(jié)果對比發(fā)現(xiàn),該地震的震源機制的反演結(jié)果有較好的一致性,結(jié)果較可信。

    秭歸地震;震源機制解;CAP法(Cut and Paste Method)

    地震的震源機制從一定角度描述了震源的性質(zhì)及其破裂過程,因此,及時確定地震震源機制,對于地震本身的研究、孕震機理的解釋及震后應(yīng)力的分布,具有十分重要的意義。對于震源機制理論和方法的研究,一直是地震學(xué)研究的熱點,并取得了大量重要成果。1980年,Helmberger等利用地震P波(P&Pn1)模擬[1],通過理論和實際波形的對比,采用格點搜索的方法,得到了研究地震的震源機制。姚振興(Yao and Harkrider,1983)等通過廣義射線理論[2],利用30°~90°范圍內(nèi)的長周期P波資料測定中強地震震源過程的地震矩張量反演方法。陳運泰、吳忠良等(Chen Y T,et al,1991;陳運泰等,1992;吳忠良等,1994)利用近震源寬頻帶記錄資料進(jìn)行了矩張量反演研究[3]。近年來,隨著區(qū)域數(shù)字化地震計的改造,獲取了大量高質(zhì)量寬頻地震數(shù)據(jù),因此利用近震地震波形反演震源機制成為可能。本文擬采用Cut and Paste(CAP)[4-7]方法,探討了該地區(qū)中小地震的發(fā)震特征及區(qū)域構(gòu)造應(yīng)力場特征。

    1 方法原理部分

    利用Zhao等[4-8]提出的CAP(Cut and Paste Method;Zhao and Helmberger,1994,1996)方法,通過分割波形記錄為Pnl和Snl部分,分別賦予不同的權(quán)重,通過格點搜索的方法進(jìn)行地震震源機制的反演?;驹恚涸O(shè)u(t)是臺站記錄到的去除儀器響應(yīng)后的地震波形,s(t)是相應(yīng)的理論波形。一個雙力偶震源的理論合成位移s(t)可以表示為

    式中,i=1,2,3時分別對應(yīng)垂直走向滑動斷層、垂直傾向滑動斷層和45°傾向滑動斷層3種最基本的斷層類型。M0為標(biāo)量地震矩,Gi為格林函數(shù),Ai為震源的輻射花樣信息,θ為臺站方位角,φ,δ,λ分別代表斷層的走向、傾角、滑動角。

    采用頻率波數(shù)方法(F K)[7,8],計算各震中距的格林函數(shù),由格林函數(shù)得到合成地震圖后,把合成與觀測的地震數(shù)據(jù)做互相關(guān)。

    由互相關(guān)函數(shù)確定時間偏移:

    在反演過程中,以合成地震位移與觀測地震位移一致作為判斷標(biāo)準(zhǔn):

    通過計算合成地震圖,對比P波、pP波和sP波的理論與觀測數(shù)據(jù)的到時,準(zhǔn)確定出震源深度;根據(jù)不同震相的幅度差異,調(diào)整震源參數(shù),使得波形能夠較好地吻合觀測數(shù)據(jù)確定震源機制解。

    由于要求的震源機制解的未知參數(shù)不多,且0°≤φi≤360°,0°≤δ≤90°,-180°≤λ≤180°,因此求解方程(2)可定義一個誤差目標(biāo)函數(shù)來衡量s和u的差異,直接采用網(wǎng)格搜索的方法得到最佳的震源機制解。本研究中,使用經(jīng)震中距矯正后的絕對誤差值作為誤差目標(biāo)函數(shù)[7](Zhao and Helmberger,1996),定義為

    式中:r為震中距;r0為選定的參考震中距;p為指數(shù)因子,一般而言,體波p=1,面波p=0.5。

    2 觀測資料和地殼模型選取

    2.1 觀測臺站分布

    選用湖北省地震臺網(wǎng)的6個臺站,在0.05~20 Hz的頻率范圍內(nèi),這6個臺站儀器的幅頻特性曲線都是平直的,較好地獲得了臺站信息資料(表1)。

    表1 臺站表Table 1 Coordinates of stations used in this study

    圖1 地震震中與臺站分布圖Fig.1 Seism ic epicenter and stations arrangement

    2.2 地殼速度模型的選取

    通過比較陳學(xué)波、廖武林、李強[9]等對三峽地區(qū)模型的研究成果,選取走時殘差最小的地殼模型,經(jīng)實際驗證綜合分析得到最接近于真實的速度模型,從地表到40 km深度處分為8層,建立初始模型計算理論格林函數(shù),見表2。

    3 計算結(jié)果與分析

    在反演過程中,將理論波形與旋轉(zhuǎn)后的觀測波形,濾波后得到Pnl部分(0.05~0.2 Hz)和Snl部分(0.02~0.05 Hz)進(jìn)行擬合[10-12],全空間中進(jìn)行格點搜索震源參數(shù),搜索得到不同深度上的震源機制和誤差,調(diào)整震源參數(shù),使得波形能夠較好的吻合觀測數(shù)據(jù),得到最佳的震源機制解。

    本文利用CAP方法計算了湖北省秭歸Ms4.1與Ms3.2級地震的震源機制,結(jié)果可靠,增加了該區(qū)中小地震震源機制解的基礎(chǔ)資料。

    表2 湖北秭歸地區(qū)地殼速度模型Table 2 Crustal velocity m odel for Zigui area,Hubei Province

    3.1 震源深度

    震源深度始終是一個很難確定的一個參數(shù),往往地震目錄中給定的深度可能與實際深度存在較大的差異。通過F K方法,計算基于上述速度模型和不同震中距的格林函數(shù)時,震源深度對理論地震圖的形態(tài)將產(chǎn)生直接影響,并在一定程度上影響震源機制的反演結(jié)果。為了減小震源深度誤差對震源機制反演結(jié)果的影響,反演過程中采用格點法對震源深度進(jìn)行了搜索,以波形擬合誤差最小的震源深度和相應(yīng)的震源機制解作為最佳擬合結(jié)果。圖2是地震事件2008年11月22日的例子,圖2(a)是震源深度在7 km附近時的位移譜擬合結(jié)果;圖2(b)是誤差隨著深度的分布,可以看出震源深度的變化對擬合誤差具有明顯的影響,但在最佳深度附近的震源機制變化不大。根據(jù)誤差最小判定,本次地震事件的最佳震源深度為7 km,波形擬合為最好,最符合實際震源機制結(jié)果。

    比較地震觀測報告中給出的震源深度(表3中H1)和震源機制反演過程中根據(jù)波形擬合得到的最佳震源深度(H2)的比較結(jié)果,發(fā)現(xiàn)兩者定位深度比較接近,說明此2次地震集中在5~10 km范圍內(nèi),而后者由于震級太小,誤差隨深度的變化結(jié)果并不是很好,地震深度的可信度可以為判斷本次地震可能是仙女山斷裂受庫水滲透和載荷變化影響引發(fā)的水庫構(gòu)造地震提供依據(jù)。

    3.2 結(jié)果分析

    本文共計算了2次地震的震源參數(shù)與相應(yīng)的震源機制解(如表3所示)。根據(jù)表中地震的2個斷層面以及P軸、T軸和B軸方向,分析得到:這2次地震的震源機制以逆沖為主兼有走滑分量;主壓應(yīng)力P軸主要集中在NE向,這與該區(qū)主要的北東向地震構(gòu)造基本一致,初步確定該震與北東向構(gòu)造活動密切有關(guān),需要指出的是,P波初動資料所得震源機制解中有2組節(jié)面,而極震區(qū)等震線長軸約為北西西走向,那么震中處應(yīng)當(dāng)同時存在北西西走向的斷裂構(gòu)造。分析認(rèn)為本次地震可能是仙女山斷裂受庫水滲透和載荷變化影響引發(fā)的水庫構(gòu)造地震。

    圖2 CAP方法計算的地震事件2008112216震源機制解Fig.2 Focalmechanism solution results of the earthquakes on Nov.22 2008 w ith the CAPmethod

    表3 湖北 秭歸Ms4.1與Ms3.2地震的震源機制解Table 3 Focalmechanism solution results of the Ms 4.1 and Ms 3.2 earthquakes in Hubei Zigui

    4 結(jié)論與討論

    4.1 結(jié) 論

    (1)本文根據(jù)波形擬合得到的最佳震源深度(H2)的比較結(jié)果主要集中分布在10 km附近,與李強[9]等通過雙差定位對三峽地區(qū)對深度的統(tǒng)計的結(jié)論一致(88.46%的地震深度分布在0~10 km范圍內(nèi),其中42.55%的震源深度在6 km以內(nèi),平均深度為6.59 km),說明定位震源深度都比較好且本文選擇的深度合理,對震源機制解結(jié)果沒有起到影響作用。

    (2)4.1級主震的最佳雙力偶解為節(jié)面I走向213°,傾角81°,滑動角110°;節(jié)面Ⅱ走向326°,傾角22°,滑動角25°;3.2級地震的最佳雙力偶解節(jié)面Ⅰ走向211°,傾角27°,滑動角109°;節(jié)面Ⅱ走向10°,傾角65°,滑動角81°。這2次地震的震源機制解有一些差別,但均為逆斷兼有走滑性質(zhì),主壓應(yīng)力P軸主要集中在NE向,與區(qū)域主壓應(yīng)力方向有一定一致性,分析認(rèn)為本次地震可能是仙女山斷裂受庫水滲透和載荷變化影響引發(fā)的水庫構(gòu)造地震。

    (3)通過與P波初動地震震源機制結(jié)果(節(jié)面Ⅰ走向213°,傾角81°,滑動角110°)對比發(fā)現(xiàn),該地震的震源機制的反演結(jié)果有較好的一致性,結(jié)果較可信。

    (4)本文給出2008年09月27日秭歸Ms3.2地震,2008年11月22日秭歸Ms4.1的發(fā)震構(gòu)造是北東向的逆斷兼有走滑性質(zhì)斷層,這與該區(qū)主要的北東向地震構(gòu)造基本一致,初步確定該地震與北東向構(gòu)造活動密切有關(guān),需要指出的是,P波初動資料所得震源機制解中有2組節(jié)面,而極震區(qū)等震線長軸約為北西西走向,那么震中處應(yīng)當(dāng)同時存在北西西走向的斷裂構(gòu)造[13]。

    4.2 討 論

    本文是首次在三峽地區(qū)利用數(shù)字波形資料反演震源參數(shù),盡管是初步分析,但也得到了比較有益的結(jié)果。隨著數(shù)字化地震臺網(wǎng)的改造,臺網(wǎng)密度增大、信噪比提高,記錄了越來越多的的中小地震,因此對中小地震震源機制的研究將會起到越來越重要的作用和良好的應(yīng)用前景。湖北三峽地區(qū)應(yīng)力狀態(tài)可根據(jù)多個中小地震震源機制解詳細(xì)確定應(yīng)力場變化的細(xì)節(jié)特征,有待今后的進(jìn)一步研究。

    參考文獻(xiàn):

    [1] HELMBERGER D V,ENGEN GR.Modeling the long pe riod body waves from shallow earthquakes at regional ran ges[J].Bul1.Seism.Soc.Am.,1980,70:1699-1714.

    [2] 姚振興,鄭天愉,溫聯(lián)星,等.用P波波形資料反演中強地震矩張量的方法[J].地球物理學(xué)報,1994,37(1):36-44.(YAO Zhen xin,ZHEN Tian yun,WEN Lian xing,et a1.Moment tensor inversion of the great or moderate earthquakes with P wave records[J].Chinese Journal of Geophysics.1994,37(1):36-44.(in Chi nese))

    [3] 吳忠良,陳運泰,倪江川,等.近震源寬頻帶記錄的地震矩張量反演[J].地震學(xué)報,1994,(2):141-152.(WU Zhong liang,CHEN Yun tai,NI Jiang chuan,et a1.Moment Tensor Inversion of broadband records local earthquake[J].Acta Seismologica Sinica.1994,(2):141-152.(in Chinese))

    [4] ZHAO L S,HELMBERGER D V.Source estimation from broadband regional seismograms[J].Bull.Seis.Soc.A mer,1994,84(1):91-104.

    [5] ZHU L,HELMBERGER D V.Advancement in source es timation techniques using broadband regional seismograms[J].Bull Seism Soc Am,1996,86:1634-1641.

    [6] ZHU L,TAN Y,HELMBERGER D V.Calibration of the Tibetan Plateau using regional seismic waveforms[J].Pa geoph,2006,163:1193-1213.

    [7] TAN Y,ZHU L,HELMBERGER DV,et a1.Locating and modeling regional earthquakes with two stations[J].J Geophys.Res.,2006,111(B01):306-314.

    [8] ZHU and Rivera.A note on the dynamic and static dis placements from a point source in multilayered media[J].Geophys.J.Int.,2002,148:619-627.

    [9] 張愛民,李強.雙差定位法在三峽庫區(qū)地震定位中的應(yīng)用[J].大地測量與地球動力學(xué),2006,26(3):73-77.(ZHANG Ai min,LIQiang.Application of double difference locationmethod to earthquake in Three Gorges Reservoir area[J].Journal of Geodesy and Geodynam ics,2006,26(3):73-77.(in Chinese))

    [10]張愛萍,倪四道,楊曉勇,等.2004年東烏珠穆沁旗地震震源參數(shù)研究[J].地震,2008,28(3):61-68.(ZHANG Ai ping,NISi dao,YANG Xiao yong,et al.Source parameters of 2004 Xilin Gol League earthquake[J].Earthquake,2008,28(3):61-68.(in Chinese))

    [11]黃建平,倪四道,傅容珊,等.綜合近震及遠(yuǎn)震波形反演2006文安地震(Mw5.1)的震源機制解[J].地球物理學(xué)報,2009,52(1):120-130.(HUANG Jian ping,NISi dao,F(xiàn)U Rong shan,et a1.Source mechanism of the 2006 Mw5.1 Wenan Earthquake determined from a joint inversion of local and teleseismic broadband wave form data[J].Chinese Journal of Geophysics,2009,52(1):120-130.(in Chinese))

    [12]呂 堅,鄭 勇,倪四道,等.2005年11月26日九江-瑞昌Ms5.7、Ms4.8地震的震源機制解與發(fā)震構(gòu)造研究[J].地球物理學(xué)報,2008,51(1):158-164.(Lv Jian,ZHENG Yong,NISi dao,et al.Focalmechanisms and seismogenic structures of the Ms 5.7 and Ms4.8 Jiu jiang Ruichang earthquakes of Nov.26,2005[J].Chi nese Journal of Geophysics,2008,51(1):158-164.(in Chinese))

    [13]楊福平,嚴(yán)尊國.湖北省秭歸ML3.6地震震源特征[J].大地測量與地球動力學(xué),2002,22(4):92-96.(YANG Fu ping,YAN Zun guo.Source characteristics of ML3.6 earthquake in Zigui,Hubei[J].Journal of Ge odesy and Geodynamics,2002,22(4):92-96(in Chi nese))

    (編輯:劉運飛)

    Study on Source M echanism of M s 4.1 and M s 3.2 Earthquakes in Hubei Ziguiw ith CAP M ethod

    ZHAO Ling yun1,DENG Jin2,CHEN Jun hua1,SHEN Xue lin1,DAIMiao1,WEIGui chun1
    (1.Earthquake Administration of Hubei Province,Wuhan 430071,China;2.Earthquake Administration of Gansu Province,Lanzhou 730000,China)

    On the basis of the records of6 stations in Hubei regional seismic network,we obtained the focalmecha nisms of the Ms 4.1 and Ms3.2 earthquakes in Ziguion Nov.22,2008 and Sep.27,2008,respectivelywith the Cut and Paste(CAP)method.The result shows that earthquakes have strike slip and strike slip faultmechanisms and the orientations ofmaximum compressional stresses aremainly in NE.It is certainly consistentwith the orientation of the regionalmaximum compressional stresses.We deduce the earthquake is the reservoir induced earthquake by water permeation and dynamic changes load.The comparison shows that the focusmechanism of two events presen ted in this paper is consistentwith that of small moderate earthquakes determined by P wave first motion.

    Zigui earthquake;focusmechanism solution;CAPmethod

    P315.6

    A

    1001-5485(2010)05-0081-04

    2010 03 01

    趙凌云(1981 ),女,碩士,山西五寨人,主要從事數(shù)字地震波理論及應(yīng)用研究,(電話)027 87868209(電子信箱)zhaolingyun002@126.com。

    猜你喜歡
    秭歸臺站震源
    中國科學(xué)院野外臺站檔案工作回顧
    氣象基層臺站建設(shè)
    西藏科技(2021年12期)2022-01-17 08:46:38
    震源的高返利起步
    匠人:在秭歸打造一艘傳統(tǒng)龍舟
    中國三峽(2016年6期)2017-01-15 13:59:18
    秭歸端午節(jié)日文化的十大類型
    中國三峽(2016年6期)2017-01-15 13:59:13
    原味·秭歸端午
    中國三峽(2016年6期)2017-01-15 13:59:02
    秭歸端午:楚鄉(xiāng)遺俗至今留
    中國三峽(2016年6期)2017-01-15 13:59:02
    可控震源地震在張掖盆地南緣逆沖斷裂構(gòu)造勘探中的應(yīng)用
    基層臺站綜合觀測業(yè)務(wù)管理之我見
    西藏科技(2015年6期)2015-09-26 12:12:13
    同步可控震源地震采集技術(shù)新進(jìn)展
    亚洲国产精品sss在线观看 | 深夜精品福利| 激情在线观看视频在线高清| 国产区一区二久久| 久久伊人香网站| 久久午夜综合久久蜜桃| 欧美日韩亚洲国产一区二区在线观看| 50天的宝宝边吃奶边哭怎么回事| 黄色a级毛片大全视频| 在线国产一区二区在线| www.熟女人妻精品国产| 午夜亚洲福利在线播放| a级片在线免费高清观看视频| 欧美午夜高清在线| 男女床上黄色一级片免费看| 亚洲精品一卡2卡三卡4卡5卡| 精品少妇一区二区三区视频日本电影| 999久久久精品免费观看国产| 亚洲色图av天堂| 一级毛片女人18水好多| 80岁老熟妇乱子伦牲交| 日本五十路高清| 搡老熟女国产l中国老女人| 中文字幕精品免费在线观看视频| 国产熟女午夜一区二区三区| 丰满的人妻完整版| 亚洲欧美日韩无卡精品| 黑人巨大精品欧美一区二区mp4| 精品卡一卡二卡四卡免费| 日韩精品中文字幕看吧| 看黄色毛片网站| 国产激情欧美一区二区| 亚洲欧美精品综合一区二区三区| 免费在线观看影片大全网站| 精品人妻在线不人妻| 久久中文看片网| 亚洲性夜色夜夜综合| 亚洲aⅴ乱码一区二区在线播放 | 在线观看免费日韩欧美大片| 亚洲国产精品一区二区三区在线| 18禁黄网站禁片午夜丰满| 亚洲人成电影免费在线| 欧美大码av| 韩国av一区二区三区四区| 国产精品免费视频内射| 亚洲精品一区av在线观看| 国产蜜桃级精品一区二区三区| 日本黄色日本黄色录像| 久久久久久久精品吃奶| 免费搜索国产男女视频| 夜夜夜夜夜久久久久| 99热国产这里只有精品6| 成人免费观看视频高清| 午夜福利在线观看吧| 18禁美女被吸乳视频| 一二三四社区在线视频社区8| 久久性视频一级片| 制服人妻中文乱码| a级毛片黄视频| 天天躁狠狠躁夜夜躁狠狠躁| 长腿黑丝高跟| 国产精品偷伦视频观看了| 99久久精品国产亚洲精品| 国产成人影院久久av| 久99久视频精品免费| 性少妇av在线| 欧美大码av| 国产亚洲精品综合一区在线观看 | 精品久久久久久久久久免费视频 | 视频区图区小说| 精品日产1卡2卡| 看免费av毛片| 成年版毛片免费区| 日本五十路高清| 免费av毛片视频| 亚洲av日韩精品久久久久久密| 在线看a的网站| 久久久国产成人免费| 亚洲成国产人片在线观看| 久久久久久久久免费视频了| 少妇 在线观看| 婷婷精品国产亚洲av在线| a级片在线免费高清观看视频| 午夜亚洲福利在线播放| 97超级碰碰碰精品色视频在线观看| 在线观看舔阴道视频| 中文亚洲av片在线观看爽| 午夜影院日韩av| 亚洲精品一二三| 国产亚洲精品久久久久5区| 桃红色精品国产亚洲av| 中出人妻视频一区二区| 男男h啪啪无遮挡| 在线观看66精品国产| 乱人伦中国视频| 男人舔女人下体高潮全视频| 久久久久久久午夜电影 | 男女之事视频高清在线观看| 久热这里只有精品99| 免费人成视频x8x8入口观看| 亚洲成a人片在线一区二区| 最新美女视频免费是黄的| 在线播放国产精品三级| 欧美 亚洲 国产 日韩一| 亚洲 欧美 日韩 在线 免费| 国产深夜福利视频在线观看| 亚洲激情在线av| 亚洲九九香蕉| 99久久国产精品久久久| 中文字幕最新亚洲高清| 天天影视国产精品| 久久精品国产亚洲av香蕉五月| 免费人成视频x8x8入口观看| a在线观看视频网站| 久久久国产成人精品二区 | 十八禁网站免费在线| 极品教师在线免费播放| 人人妻人人澡人人看| √禁漫天堂资源中文www| 99久久99久久久精品蜜桃| 久热爱精品视频在线9| 国产男靠女视频免费网站| 97碰自拍视频| 国产一区二区激情短视频| 69av精品久久久久久| 99久久国产精品久久久| 在线av久久热| 亚洲av片天天在线观看| 日韩免费av在线播放| 岛国视频午夜一区免费看| 最近最新中文字幕大全免费视频| 最新在线观看一区二区三区| 免费观看精品视频网站| 亚洲五月色婷婷综合| 亚洲成人精品中文字幕电影 | 亚洲精品久久成人aⅴ小说| 午夜免费观看网址| 国产精品久久视频播放| 午夜成年电影在线免费观看| 亚洲午夜理论影院| 精品国产一区二区三区四区第35| 在线观看午夜福利视频| 91在线观看av| 免费在线观看影片大全网站| 另类亚洲欧美激情| 久久中文字幕一级| 男人操女人黄网站| 怎么达到女性高潮| 一二三四在线观看免费中文在| 免费看a级黄色片| 亚洲精品美女久久久久99蜜臀| 夜夜爽天天搞| 老司机福利观看| 亚洲av成人一区二区三| 久热爱精品视频在线9| 国产精品久久久av美女十八| 国产蜜桃级精品一区二区三区| 国产成人精品在线电影| 香蕉国产在线看| 亚洲avbb在线观看| 欧美国产精品va在线观看不卡| 色婷婷久久久亚洲欧美| 亚洲国产精品sss在线观看 | 欧美一区二区国产精品久久精品| 国产精品爽爽va在线观看网站| 欧美成人一区二区免费高清观看| 亚洲国产日韩欧美精品在线观看| 神马国产精品三级电影在线观看| 亚洲av成人av| 欧美丝袜亚洲另类 | 亚洲第一区二区三区不卡| 精品久久久久久久久av| 别揉我奶头 嗯啊视频| 国产午夜福利久久久久久| 成年女人永久免费观看视频| 亚洲人成网站在线播| 在线天堂最新版资源| 亚洲精品粉嫩美女一区| 蜜桃久久精品国产亚洲av| 亚洲电影在线观看av| 久久久久久国产a免费观看| 久久久国产成人精品二区| 免费av不卡在线播放| 国产国拍精品亚洲av在线观看| 麻豆av噜噜一区二区三区| 在线观看66精品国产| 亚洲内射少妇av| 男人狂女人下面高潮的视频| 三级国产精品欧美在线观看| 91在线精品国自产拍蜜月| 一卡2卡三卡四卡精品乱码亚洲| 亚洲中文字幕日韩| 日韩欧美免费精品| 免费人成在线观看视频色| 亚洲美女黄片视频| 国产亚洲精品久久久com| 91九色精品人成在线观看| 搡女人真爽免费视频火全软件 | 天堂√8在线中文| 欧美中文日本在线观看视频| 国产免费一级a男人的天堂| 国产精品女同一区二区软件 | 精品一区二区三区视频在线观看免费| 搡老妇女老女人老熟妇| 亚洲av中文字字幕乱码综合| 久久久成人免费电影| 国产美女午夜福利| 欧美日韩综合久久久久久 | 啪啪无遮挡十八禁网站| 色综合婷婷激情| 热99在线观看视频| 免费在线观看成人毛片| 不卡一级毛片| 天美传媒精品一区二区| 男女那种视频在线观看| 亚洲精品日韩av片在线观看| 无人区码免费观看不卡| 国产亚洲精品久久久久久毛片| 舔av片在线| 亚洲男人的天堂狠狠| 身体一侧抽搐| 啪啪无遮挡十八禁网站| 国产伦精品一区二区三区四那| 国产真实伦视频高清在线观看 | 精品久久久久久久久久久久久| 麻豆久久精品国产亚洲av| 亚洲性夜色夜夜综合| 亚洲男人的天堂狠狠| 免费在线观看影片大全网站| 黄色视频,在线免费观看| 午夜久久久久精精品| av在线观看视频网站免费| 欧美成人一区二区免费高清观看| 国产精品乱码一区二三区的特点| 九色成人免费人妻av| 亚洲欧美精品综合久久99| 一进一出抽搐gif免费好疼| 啪啪无遮挡十八禁网站| 欧美乱色亚洲激情| 国产一区二区亚洲精品在线观看| 日本三级黄在线观看| 99国产极品粉嫩在线观看| 男女视频在线观看网站免费| 老司机深夜福利视频在线观看| 国产成年人精品一区二区| 欧美中文日本在线观看视频| 午夜福利在线在线| 精品欧美国产一区二区三| 国产毛片a区久久久久| 国产一区二区在线av高清观看| 精品99又大又爽又粗少妇毛片 | 变态另类丝袜制服| 日本一本二区三区精品| 国产蜜桃级精品一区二区三区| 色视频www国产| 免费在线观看影片大全网站| 美女cb高潮喷水在线观看| 国产探花极品一区二区| 极品教师在线免费播放| 精品国内亚洲2022精品成人| 欧美在线黄色| 少妇人妻精品综合一区二区 | 18禁在线播放成人免费| 中文字幕av在线有码专区| 黄片小视频在线播放| 国产单亲对白刺激| 亚洲精品在线观看二区| 午夜精品久久久久久毛片777| h日本视频在线播放| 国产精品精品国产色婷婷| 久久久久久久亚洲中文字幕 | 观看免费一级毛片| 搡女人真爽免费视频火全软件 | 国产精品影院久久| 久久午夜亚洲精品久久| 国产视频一区二区在线看| 超碰av人人做人人爽久久| 成人一区二区视频在线观看| 亚洲国产欧洲综合997久久,| 国产色爽女视频免费观看| 又粗又爽又猛毛片免费看| 日韩高清综合在线| h日本视频在线播放| 深爱激情五月婷婷| 亚洲国产日韩欧美精品在线观看| 亚洲无线在线观看| a级一级毛片免费在线观看| 丰满人妻一区二区三区视频av| 老司机福利观看| 亚洲精品一区av在线观看| 日韩亚洲欧美综合| 两性午夜刺激爽爽歪歪视频在线观看| 听说在线观看完整版免费高清| 亚洲av成人精品一区久久| 男女做爰动态图高潮gif福利片| 欧美色欧美亚洲另类二区| 悠悠久久av| 亚洲精华国产精华精| 免费在线观看日本一区| 一二三四社区在线视频社区8| 高清在线国产一区| 国产免费av片在线观看野外av| 亚洲一区高清亚洲精品| 国产av在哪里看| 青草久久国产| 变态另类成人亚洲欧美熟女| 亚洲黑人精品在线| 国内精品一区二区在线观看| 欧美性猛交╳xxx乱大交人| 国产精品野战在线观看| 亚洲国产欧美人成| 麻豆一二三区av精品| 亚洲国产欧洲综合997久久,| 亚洲欧美日韩高清在线视频| 免费高清视频大片| 精品一区二区三区av网在线观看| 草草在线视频免费看| 色综合欧美亚洲国产小说| 国产三级在线视频| 久久久久性生活片| 99国产极品粉嫩在线观看| 亚洲成人久久爱视频| 成人特级黄色片久久久久久久| 91久久精品国产一区二区成人| 国产成+人综合+亚洲专区| 五月玫瑰六月丁香| 国产精品一区二区免费欧美| 亚洲精品色激情综合| 欧美日韩亚洲国产一区二区在线观看| 久久国产乱子伦精品免费另类| 可以在线观看的亚洲视频| 1000部很黄的大片| 久久久色成人| 国产免费一级a男人的天堂| 男插女下体视频免费在线播放| 日本成人三级电影网站| 国模一区二区三区四区视频| 黄色日韩在线| 日韩欧美精品免费久久 | 久99久视频精品免费| 三级国产精品欧美在线观看| 亚洲国产色片| 日本免费一区二区三区高清不卡| 欧美日韩福利视频一区二区| 国产精品综合久久久久久久免费| 成人性生交大片免费视频hd| 又粗又爽又猛毛片免费看| 成人性生交大片免费视频hd| 91午夜精品亚洲一区二区三区 | 午夜免费男女啪啪视频观看 | 免费电影在线观看免费观看| 中文在线观看免费www的网站| 一区福利在线观看| 亚洲不卡免费看| 3wmmmm亚洲av在线观看| 精品日产1卡2卡| 97人妻精品一区二区三区麻豆| 午夜激情欧美在线| 欧美日本亚洲视频在线播放| 日本黄色片子视频| 亚洲自拍偷在线| 欧美黄色淫秽网站| www日本黄色视频网| 欧美中文日本在线观看视频| 国产一区二区在线观看日韩| 色av中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 99热这里只有精品一区| 九九久久精品国产亚洲av麻豆| 久久99热这里只有精品18| 国产视频内射| 国产白丝娇喘喷水9色精品| 国产av不卡久久| 久久久久久久亚洲中文字幕 | 欧美成人a在线观看| 日日摸夜夜添夜夜添av毛片 | 1024手机看黄色片| 亚洲中文字幕一区二区三区有码在线看| 嫩草影院入口| 国产黄片美女视频| 十八禁人妻一区二区| 国产白丝娇喘喷水9色精品| 一夜夜www| 欧美绝顶高潮抽搐喷水| 国产伦人伦偷精品视频| 欧美精品啪啪一区二区三区| 免费人成视频x8x8入口观看| 国产精品98久久久久久宅男小说| 国产精品久久视频播放| 伦理电影大哥的女人| 一个人观看的视频www高清免费观看| 又黄又爽又免费观看的视频| 搡老妇女老女人老熟妇| 少妇高潮的动态图| 免费搜索国产男女视频| h日本视频在线播放| 91久久精品国产一区二区成人| 内射极品少妇av片p| 偷拍熟女少妇极品色| 韩国av一区二区三区四区| 变态另类丝袜制服| 99久久99久久久精品蜜桃| 露出奶头的视频| 看黄色毛片网站| 变态另类丝袜制服| 亚洲成av人片在线播放无| 露出奶头的视频| 日日摸夜夜添夜夜添av毛片 | 性欧美人与动物交配| 91av网一区二区| 欧美3d第一页| 男女下面进入的视频免费午夜| 精品99又大又爽又粗少妇毛片 | 性色avwww在线观看| 国产老妇女一区| 国产亚洲av嫩草精品影院| 亚洲电影在线观看av| 免费观看的影片在线观看| 三级国产精品欧美在线观看| 午夜激情福利司机影院| 久久人人精品亚洲av| 91麻豆av在线| 亚洲av.av天堂| 国产大屁股一区二区在线视频| 高清日韩中文字幕在线| 成人一区二区视频在线观看| 舔av片在线| 久久精品国产亚洲av涩爱 | 国产精品一区二区性色av| 别揉我奶头~嗯~啊~动态视频| 成人国产一区最新在线观看| 国产精品伦人一区二区| 国产午夜精品论理片| 国内毛片毛片毛片毛片毛片| 亚洲av五月六月丁香网| 免费观看的影片在线观看| 日韩欧美 国产精品| 精品无人区乱码1区二区| 精品久久久久久,| 欧美另类亚洲清纯唯美| 欧美黄色片欧美黄色片| 国产在线男女| 男女之事视频高清在线观看| 在线a可以看的网站| 美女cb高潮喷水在线观看| 日韩欧美 国产精品| 精品无人区乱码1区二区| 亚洲国产精品sss在线观看| 99热这里只有精品一区| 欧美黑人欧美精品刺激| 国产在线男女| 欧美成人性av电影在线观看| 亚洲成人免费电影在线观看| 欧美色视频一区免费| 亚洲专区中文字幕在线| 中出人妻视频一区二区| 久久中文看片网| 国产欧美日韩精品一区二区| 午夜亚洲福利在线播放| 婷婷精品国产亚洲av| 一本一本综合久久| 国产av在哪里看| 中出人妻视频一区二区| 国产精品一及| 中文字幕高清在线视频| 性色av乱码一区二区三区2| 白带黄色成豆腐渣| 欧美乱妇无乱码| 亚洲av成人不卡在线观看播放网| 国产精品伦人一区二区| 国产精品一及| 中文字幕高清在线视频| 日韩av在线大香蕉| 国产精品一区二区性色av| 亚洲欧美激情综合另类| 午夜福利欧美成人| 久久国产精品人妻蜜桃| 久久久久久大精品| 99久久九九国产精品国产免费| 日本 欧美在线| 国产精品一区二区性色av| 在线看三级毛片| av欧美777| 午夜精品一区二区三区免费看| 亚洲av日韩精品久久久久久密| 日本熟妇午夜| 欧美乱妇无乱码| 国产av在哪里看| 日韩欧美在线二视频| 免费高清视频大片| 熟妇人妻久久中文字幕3abv| 俺也久久电影网| 国产三级中文精品| 欧美高清性xxxxhd video| 亚洲国产欧洲综合997久久,| 国产伦在线观看视频一区| 欧美xxxx黑人xx丫x性爽| 老熟妇仑乱视频hdxx| 久久草成人影院| 黄片小视频在线播放| 在线国产一区二区在线| 欧美zozozo另类| 亚洲一区二区三区色噜噜| 日本黄色片子视频| 久久久成人免费电影| 淫妇啪啪啪对白视频| 综合色av麻豆| 免费av不卡在线播放| 久99久视频精品免费| 五月伊人婷婷丁香| 国产色爽女视频免费观看| 久久久精品欧美日韩精品| 国产黄片美女视频| 国产美女午夜福利| 国内揄拍国产精品人妻在线| 美女被艹到高潮喷水动态| а√天堂www在线а√下载| 欧美日韩综合久久久久久 | 婷婷丁香在线五月| 亚洲,欧美精品.| 欧美日本视频| 久久天躁狠狠躁夜夜2o2o| 亚洲成人久久性| 午夜激情欧美在线| 九九久久精品国产亚洲av麻豆| 国产av麻豆久久久久久久| 久久久久久久精品吃奶| 国产高清视频在线播放一区| 欧美极品一区二区三区四区| 国产精品一区二区三区四区久久| 欧美日韩国产亚洲二区| 男人的好看免费观看在线视频| 97碰自拍视频| 亚洲美女搞黄在线观看 | 亚洲av.av天堂| 9191精品国产免费久久| 国产成人影院久久av| 亚洲精品一区av在线观看| 午夜久久久久精精品| 90打野战视频偷拍视频| 麻豆成人av在线观看| 老司机午夜十八禁免费视频| 高清日韩中文字幕在线| 91狼人影院| 青草久久国产| 一级a爱片免费观看的视频| 欧美一级a爱片免费观看看| 超碰av人人做人人爽久久| 成年女人看的毛片在线观看| 日本黄色片子视频| 国产成年人精品一区二区| 特级一级黄色大片| 国产欧美日韩一区二区精品| 少妇丰满av| 真人做人爱边吃奶动态| 69人妻影院| 欧美不卡视频在线免费观看| 18禁裸乳无遮挡免费网站照片| 国产极品精品免费视频能看的| 美女 人体艺术 gogo| 色在线成人网| 欧美高清性xxxxhd video| 99在线人妻在线中文字幕| 欧美不卡视频在线免费观看| 51国产日韩欧美| 亚洲综合色惰| 国产亚洲精品av在线| 国产精华一区二区三区| 国产男靠女视频免费网站| 久久国产精品人妻蜜桃| 嫩草影院新地址| 国产免费一级a男人的天堂| 亚洲人成网站在线播| 综合色av麻豆| 亚洲国产高清在线一区二区三| 成年女人毛片免费观看观看9| 日韩欧美三级三区| 成人高潮视频无遮挡免费网站| 国产精品美女特级片免费视频播放器| 久久婷婷人人爽人人干人人爱| 亚洲国产日韩欧美精品在线观看| 毛片女人毛片| 男女下面进入的视频免费午夜| 色综合婷婷激情| 免费观看的影片在线观看| eeuss影院久久| 人人妻人人看人人澡| 亚洲五月婷婷丁香| 日韩中字成人| 精品久久久久久久人妻蜜臀av| a在线观看视频网站| 18禁黄网站禁片免费观看直播| a级一级毛片免费在线观看| 午夜老司机福利剧场| 成人三级黄色视频| 亚洲三级黄色毛片| 色精品久久人妻99蜜桃| 日韩欧美精品v在线| 丁香六月欧美| 俺也久久电影网| 日韩中字成人| 亚洲午夜理论影院| 在现免费观看毛片| 男女下面进入的视频免费午夜| 9191精品国产免费久久| 男人和女人高潮做爰伦理| 国产一区二区在线观看日韩| 国产精品野战在线观看| a在线观看视频网站| 3wmmmm亚洲av在线观看| 99久久99久久久精品蜜桃| 精品欧美国产一区二区三| 久久99热这里只有精品18| 国产探花极品一区二区| 国产成人av教育|