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

    利用波速擾動定量分析Hyposat方法的定位誤差

    2017-11-29 05:14:33于俊誼馬起楊
    華北地震科學 2017年4期
    關(guān)鍵詞:深度模型

    于俊誼,馬起楊

    (浙江省地震局,杭州 310013)

    利用波速擾動定量分析Hyposat方法的定位誤差

    于俊誼,馬起楊

    (浙江省地震局,杭州 310013)

    以浙江省數(shù)字地震臺網(wǎng)的臺站布局為基礎(chǔ),在浙江省內(nèi)以0.1°×0.1°網(wǎng)格點作為虛擬地震震中,利用浙江2015地殼速度模型計算了各臺站Pn、Sn、Pg、Sg震相到時,對速度模型及各震相走時增加適量擾動,進而使用Hyposat進行地震定位。定位結(jié)果與虛擬震源位置對比顯示:①對于震中位置,當速度模型中速度偏差小于0.1 km/s時,所得震中偏差基本在3 km以內(nèi),但在某些區(qū)域會出現(xiàn)無法有效定位的情況;②對于震源深度,與僅使用Pg、Sg(或Pb、Sb)震相定位的結(jié)果相比,當Pn、Sn震相參與定位時,震源深度的計算結(jié)果更接近實際震源深度,且計算結(jié)果的一致性更好。

    Hyposat;定位偏差;震源深度偏差

    0 引言

    地震定位研究是地震學研究的基礎(chǔ)內(nèi)容。隨著地震學家對震相走時理論研究的不斷深入、臺站到時提取精度的提高以及計算機技術(shù)的廣泛應(yīng)用,地震定位的精度得到了不斷的提高?;诓煌牡卣鹫鹣嗉白邥r方程,研究人員開發(fā)了不同的地震定位方法,目前使用較為普遍的定位方法主要有Hypo81[1]、單純型法[2]、Hypo2000[3]、Hyposat、網(wǎng)格搜索法[4]、遺傳算法[5]等地震定位算法。

    測震臺網(wǎng)軟件JOPENS是我國地震臺網(wǎng)地震資料處理的基礎(chǔ)軟件。Hyposat定位方法外掛于該軟件的人機交互模塊MSDP中,已成為各區(qū)域臺網(wǎng)進行地震定位時使用最多的定位方法之一。Hyposat是一個用于地震定位的軟件包,開發(fā)者開發(fā)該軟件包的初衷是為了最大限度地使用觀測到的震相信息進行地震定位。在地震定位過程中使用了震相到時、反方位角及震相的射線參數(shù)等信息。Hyposat可以使用IASP91標準震相列表中的所有震相進行地震定位,也可以使用單臺記錄到不同震相的到時差作為獨立參數(shù)參與定位計算[6]。

    在使用過程中,一些學者發(fā)現(xiàn)Hyposat方法可以較好地確定震中位置,卻不能穩(wěn)定地確定震源深度。特別是在使用直達波震相進行定位時,介質(zhì)波速的小幅改變都會帶來震源深度較大的偏差,甚至出現(xiàn)無法確定震源深度的情況。一些學者認為這種現(xiàn)象是由于定位時使用了與實際情況偏差較大的速度模型所造成的。本研究基于浙江省區(qū)域地震波速度模型,通過對速度模型添加擾動,評估本區(qū)域內(nèi)地殼一維速度模型對Hyposat方法定位誤差的影響。

    1 評估方法

    影響地震定位精度的主要因素包括速度模型、臺站布局、震相的數(shù)量及震相識別精確度。為了評估Hyposat定位方法在浙江省定位誤差的分布情況,本研究在浙江省數(shù)字地震臺網(wǎng)臺站布局(圖1)的基礎(chǔ)上,以給定的速度模型計算幾種常見震相的到時(即理論走時)。由此得到虛擬地震發(fā)震時刻、震中位置、震源深度及各臺站的震相到時。為了模擬人工讀數(shù)誤差,向震相到時中加入了一定的時間擾動,然后用Hyposat方法進行定位。定位過程中適當改變定位速度模型中Pg、Pn震相的波速值(Sg、Sn震相的波速按波速比系數(shù)調(diào)整),最后用定位結(jié)果與理論震源位置的偏差來評估速度模型變化對定位誤差的影響。

    2 數(shù)據(jù)處理

    真實的地下結(jié)構(gòu)是不均勻的,各種定位方法用來定位的速度模型均為近似的真實速度結(jié)構(gòu)。本研究使用浙江2015地殼模型[7]作為理論速度模型計算各種震相的走時。浙江2015地殼速度模型如表1。

    表1 浙江2015地殼速度模型

    在浙江省區(qū)域范圍內(nèi),以經(jīng)緯度0.1°間隔的網(wǎng)格點作為虛擬地震震中,以5、10、15、20 km作為虛擬地震的震源深度,根據(jù)表1給出的速度模型按照式(1)、式(2)分別計算Pg、Sg、Pn、Sn震相的走時,由式(3)可計算得到首波盲區(qū)范圍。

    (1)

    (2)

    (3)

    (4)

    (5)

    (6)

    當震源深度為25 km時,震源處于下地殼,不再出現(xiàn)Pg、Sg震相。這時使用Pb、Sb震相替代Pg、Sg震相進行計算,下地殼地震Pb、Sb震相走時通過解方程組(4)獲得,Pn、Sn震相走時按照式(5)計算,由式(6)可計算得到首波盲區(qū)范圍。其中:Δ為震中距;h為震源深度;h1、h2分別為速度模型中第一層、第二層的厚度;v1、v2、v3為速度模型中第一層、第二層及第三層頂端對應(yīng)的波速;e1、e2為地震波由下地殼經(jīng)過康拉德界面上行至地面時的出射角和入射角。由于人工識別震相時會有一定的讀數(shù)誤差,在理論到時上添加均勻分布的-0.05~0.05 s隨機誤差以模擬讀數(shù)誤差,從而得到用于測試Hyposat方法的到時數(shù)據(jù)。本文后續(xù)所有定位測試均在該到時數(shù)據(jù)的基礎(chǔ)上完成。

    為了分析縱波速度變化對Hyposat方法定位誤差的影響,本文對表1列出速度模型中的波速進行定步長調(diào)整,得到的速度模型如表2,使用調(diào)整后的速度模型進行定位。由于第二層(即康拉德界面與莫霍面之間的部分)的波速變化對Pg的走時不產(chǎn)生影響,對Pn的走時影響較小,為了突出重點及降低工作量,在對表1列出的速度模型進行變化時未對第二層的速度進行調(diào)整。

    表2 用于定位的速度模型

    使用表2列出的5種速度模型分別以Hyposat方法進行定位。通過對比定位結(jié)果與理論震源位置的偏差可對速度模型變化帶來的定位誤差進行定量評估。本研究選擇了震中距在200 km范圍內(nèi)的臺站參與定位,并分別使用了Pg、Sg震相組合及Pg、Sg、Pn、Sn震相組合(震源深度為25 km的情況下,使用Pb、Sb震相代替Pg、Sg震相)進行定位計算。

    3 定位誤差分析

    當震源位于上地殼時,以震源深度為10 km為例。如圖2所示,以理論速度模型進行定位時,定位震中與理論震中偏差均在1 km范圍內(nèi)。隨著定位所使用的速度模型與理論速度模型之間偏差增大,震中偏差也逐漸增大,當定位使用的Pg、Pn波速偏差達到±0.1 km/s時(速度模型1、5),定位震中偏差最大值可達3 km。

    注:圖a~e為使用Pg、Sg震相組合定位震中偏差分布;圖f~j為使用Pn、Pg、Sn、Sg震相組合定位震中偏差分布;定位使用的均為表2所示的速度模型,其中:圖a、f使用了速度模型1,圖b、g使用了速度模型2,圖c、h使用了速度模型3,圖d、i使用了速度模型4,圖e、j使用了速度模型5。 圖2 945個理論深度為10 km地震的定位震中偏差

    當定位使用的速度模型與理論速度模型相差不大時,使用兩種震相組合在網(wǎng)內(nèi)部分定位的震中偏差無明顯差別。在網(wǎng)緣部分,與使用Pg、Sg、Pn、Sn震相組合相比,使用Pg、Sg震相組合得到的震中位置與理論震中位置偏差略小。網(wǎng)緣部分出現(xiàn)的這種現(xiàn)象,作者分析后認為是由浙江臺網(wǎng)接入了周邊省市的共享臺站造成的。由于外省市共享臺站的存在,浙江省邊緣區(qū)域仍屬于網(wǎng)內(nèi)區(qū)域,但Pn、Sn震相需在首波盲區(qū)范圍外才會出現(xiàn),導致Pn、Sn震相記錄的臺站在省界邊緣區(qū)域缺失,從而因震相分布不均勻而出現(xiàn)了上述現(xiàn)象。當定位使用的速度模型偏差達到0.1 km/s時,使用Pg、Sg震相組合得到的震中位置偏差更小,但在某些臺站布局的情況下會出現(xiàn)無法定位的情況(圖3)。

    震源深度為10 km時的定位深度如圖4。僅使用Pg、Sg震相計算得到的震源深度主要分布于3~15 km的范圍內(nèi),而使用Pn、Pg、Sn、Sg震相共同計算時得到的震源深度主要分布于8~11 km的范圍內(nèi)。當定位使用的Pg、Pn波速偏差達到+0.1 km/s時(模型5),在某些區(qū)域的地震使用Pg、Sg震相組合進行定位計算將無法得到有效的震源深度。使用表2所示的5種速度模型進行定位時,與使用Pg、Sg震相組合的定位深度相比,以Pn、Pg、Sn、Sg震相組合定位得到的震源深度偏差更小,且定位結(jié)果的一致性更好。

    注:圖a為使用Pg、Sg震相組合的定位結(jié)果,圖中灰色部分表示未得到有效的定位結(jié)果;圖b為使用Pg、Pn、Sg、Sn震相組合的定位結(jié)果。圖3 使用表2所示速度模型5定位的震中偏差

    注:圖a~e為使用Pg、Sg震相組合定位深度分布;圖f~j為使用Pn、Pg、Sn、Sg震相組合定位深度分布;使用的模型同圖2。圖4 945個理論深度為10 km地震的定位深度

    注:理論震源深度為10 km;圖a為使用Pg、Sg震相組合的定位結(jié)果,圖中灰色部分表示未得到有效的定位結(jié)果;圖b為使用Pg、Pn、Sg、Sn震相組合的定位結(jié)果圖5 使用表2所示速度模型5定位得到的震源深度

    當震源位于下地殼時,震中位置方面與震源深度10 km的情況基本相同(圖6)。震源深度分布見圖6~7。當僅使用Pb、Sb震相進行定位計算時得到的震源深度結(jié)果偏差較大,且結(jié)果分布較為離散(圖7a~7e);當Pn、Sn震相加入定位計算后,得到的震源深度計算結(jié)果偏差較小,且結(jié)果一致性較好(圖7f~7j)。

    當震源位于上地殼的其它深度時(震源深度為5 km、15 km、20 km)定位結(jié)果對比得到的結(jié)論與震源深度為10 km的情況基本相同,本文不再詳細羅列。

    4 結(jié)論

    本文通過設(shè)置虛擬震源,使用浙江省區(qū)域理論速度模型計算各震相走時。對速度模型添加速度擾動后,使用Hyposat方法以擾動后的速度模型進行定位計算。將定位結(jié)果與虛擬震源位置進行對比后發(fā)現(xiàn),在浙江省內(nèi)速度模型及震相組合的選取對定位結(jié)果具有較大影響。具體表現(xiàn)為:

    震中位置方面,當定位使用的速度模型偏差不大時,使用不同震相組合計算得到的震中基本一致;當定位使用的速度模型與實際情況偏差較大時,使用Pg、Sg震相組合得到的震中位置較使用Pg、Sg、Pn、Sn震相組合得到的震中位置偏差略小,但在部分區(qū)域會出現(xiàn)無法得到有效定位的結(jié)果。

    注:理論深度為25 km;圖a、b為震中位置偏差圖;圖c、d為計算得到的震源深度;圖a、c為使用Pg、Sg震相組合的定位結(jié)果;圖b、d為使用Pg、Pn、Sg、Sn震相組合的定位結(jié)果。圖6 使用表2所示速度模型5時的定位結(jié)果

    注:圖a~e為使用Pb、Sb震相組合定位深度分布;圖f~j為使用Pn、Pb、Sn、Sb震相組合定位深度分布;使用的模型同圖2。圖7 945個理論深度為25 km地震的定位深度

    震源深度方面,與使用Pg、Sg震相組合的定位深度相比,當Pn、Sn震相參與定位計算時得到的震源深度偏差會更小,且一致性更好。在浙江省區(qū)域范圍內(nèi),使用Hyposat方法進行定位時應(yīng)盡可能多地添加Pn、Sn震相參與定位,這樣可以使得到的震源位置更為準確、可靠。

    [1] 趙仲和. 多重模型地震定位程序及其在北京臺網(wǎng)的應(yīng)用[J]. 地震學報, 1983, 5(2): 242-254.

    [2] 徐果明. 反演理論及其應(yīng)用[M]. 北京: 地震出版社, 2003: 120-122.

    [3] Klein F W. 85.8-The HYPOINVERSE2000 earthquake location program[J]. International Geophysics, 2003, 81: 1619-1620.

    [4] 趙仲和, 牟磊育. 單臺地震自動定位網(wǎng)格搜索法及其MATLAB試驗[J]. 地震地磁觀測與研究, 2005, 26(4): 1-12.

    [5] 朱新運, 張帆, 于俊誼. 浙江珊溪水庫地震精細定位及構(gòu)造研究[J]. 中國地震, 2010, 26(4): 380-390.

    [6] 彼得·鮑曼. 新地震觀測實踐手冊[M]. 中國地震局監(jiān)測預(yù)報司, 譯. 北京: 地震出版社, 2006: 909-927.

    [7] 朱元清, 宋秀青, 劉雙慶. 中國地震測定參考速度結(jié)構(gòu)研究[M]. 北京: 地震出版社, 2016: 202-215.

    HypocenterLocationErrorAnalysisoftheHyposatMethodUsingSeismicVelocityPerturbation

    YU Jun-yi,MA Qi-yang

    (Zhejiang Eatrhquake Agency, Hangzhou 310013, China)

    Based on the stations of Zhejiang digital seismic network and using 0.1°×0.1° grid points in Zhejiang as the virtual epicenters, we calculated the travel time of Pn, Sn, Pg and Sg seismic phases by 2015yr crust velocity model of Zhejiang province. After adding the random error on the travel time, we used Hyposat software to calculate the hypocenter location by velocity models with different range of disturbances. Comparing the virtual hypocenter with the calculation, we found that when the P-wave velocity deviation reaches within 0.1km/s, the deviation of epicenter is less than 3km. However, even in this case, there are also some improper location results in part of numerical simulations. Comparing with only making use of the Pg and Sg (Pb, Sb) phases, on the condition that Pn, Sn, Pg and Sg (Pb,Sb) phases are simultaneously used, the calculation of focal depth is closer to the virtual given value, and the calculation results are better in consistency.

    Hyposat; location error; focal depth error

    于俊誼,馬起楊.利用波速擾動定量分析Hyposat方法的定位誤差[J]. 華北地震科學,2017,35(4):29-33.

    2017-05-25

    中國地震局星火計劃攻關(guān)項目(XH18019);浙江省地震局局科技項目(2017ZJJ03)

    于俊誼(1982—),男,山東煙臺人,工程師,主要從事地震監(jiān)測方面工作.E-mail:yqlyjy@163.com

    P315.6

    A

    1003-1375(2017)04-0029-05

    10.3969/j.issn.1003-1375.2017.04.005

    猜你喜歡
    深度模型
    一半模型
    深度理解一元一次方程
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    深度觀察
    深度觀察
    深度觀察
    深度觀察
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    av黄色大香蕉| 午夜福利在线观看免费完整高清在| 国产精品无大码| 免费不卡的大黄色大毛片视频在线观看| 免费大片18禁| 韩国av在线不卡| 国产日韩欧美亚洲二区| 精品国产一区二区三区久久久樱花| 日本色播在线视频| 久久综合国产亚洲精品| 另类精品久久| 久久av网站| 日韩人妻高清精品专区| 一级毛片电影观看| 国产精品99久久99久久久不卡 | 日本91视频免费播放| 99国产精品免费福利视频| 精品人妻偷拍中文字幕| 国产高清国产精品国产三级| 国产精品一二三区在线看| 日韩视频在线欧美| 国产精品一区二区性色av| 精品国产乱码久久久久久小说| 国产精品三级大全| 日韩 亚洲 欧美在线| 午夜影院在线不卡| 久久久亚洲精品成人影院| 精品午夜福利在线看| 一区二区三区免费毛片| 国产高清三级在线| 三级国产精品欧美在线观看| 五月玫瑰六月丁香| 欧美日韩av久久| 精品久久久噜噜| 九草在线视频观看| 亚洲怡红院男人天堂| 校园人妻丝袜中文字幕| 中文字幕免费在线视频6| 久久鲁丝午夜福利片| 哪个播放器可以免费观看大片| 一级二级三级毛片免费看| 久久久久精品性色| 精品少妇久久久久久888优播| 麻豆乱淫一区二区| 亚洲国产精品999| 亚洲精品日韩在线中文字幕| 亚洲国产av新网站| 亚洲av成人精品一二三区| a级片在线免费高清观看视频| 午夜视频国产福利| 久久毛片免费看一区二区三区| 精品一品国产午夜福利视频| 熟女人妻精品中文字幕| 下体分泌物呈黄色| 国产视频内射| 欧美 日韩 精品 国产| 亚洲图色成人| av黄色大香蕉| 高清黄色对白视频在线免费看 | 国产美女午夜福利| 色5月婷婷丁香| 免费人妻精品一区二区三区视频| 偷拍熟女少妇极品色| 免费在线观看成人毛片| 亚洲久久久国产精品| 久久99精品国语久久久| 日本av手机在线免费观看| 妹子高潮喷水视频| 亚洲在久久综合| 如日韩欧美国产精品一区二区三区 | 国产一区二区三区av在线| 亚洲av电影在线观看一区二区三区| 最近2019中文字幕mv第一页| 欧美+日韩+精品| 国产午夜精品久久久久久一区二区三区| 狂野欧美激情性xxxx在线观看| 国产 精品1| av专区在线播放| 22中文网久久字幕| 久久人人爽av亚洲精品天堂| 亚洲精品乱码久久久v下载方式| 亚洲av二区三区四区| 草草在线视频免费看| 91久久精品国产一区二区成人| av不卡在线播放| 久久久久久久久久久丰满| 日韩在线高清观看一区二区三区| 亚洲国产精品一区二区三区在线| 91久久精品国产一区二区三区| 国产黄色视频一区二区在线观看| 免费看光身美女| 五月伊人婷婷丁香| 春色校园在线视频观看| 欧美高清成人免费视频www| 2022亚洲国产成人精品| 久久青草综合色| 亚洲国产成人一精品久久久| 成人亚洲欧美一区二区av| 精品久久久久久电影网| 亚洲综合精品二区| 少妇被粗大的猛进出69影院 | 国产午夜精品一二区理论片| 欧美亚洲 丝袜 人妻 在线| 天堂8中文在线网| 美女视频免费永久观看网站| 国产精品99久久久久久久久| 永久网站在线| 国产伦理片在线播放av一区| 亚洲人成网站在线播| 精品人妻一区二区三区麻豆| 精品人妻一区二区三区麻豆| 下体分泌物呈黄色| 亚洲av免费高清在线观看| 老司机影院成人| 免费看日本二区| 日日啪夜夜爽| 一级爰片在线观看| 精品视频人人做人人爽| 午夜福利,免费看| 九草在线视频观看| 久久久午夜欧美精品| 99热国产这里只有精品6| 亚州av有码| 伦理电影免费视频| 久久6这里有精品| 51国产日韩欧美| 在线观看免费视频网站a站| 亚洲真实伦在线观看| 黄色怎么调成土黄色| 免费看日本二区| 黑丝袜美女国产一区| 亚洲精品色激情综合| 乱系列少妇在线播放| 熟女av电影| 日韩精品有码人妻一区| 中文乱码字字幕精品一区二区三区| 18+在线观看网站| 国产 一区精品| 一区二区三区免费毛片| 亚洲丝袜综合中文字幕| 亚洲欧美精品专区久久| 大片免费播放器 马上看| 精品亚洲成a人片在线观看| 久久人妻熟女aⅴ| 日韩强制内射视频| 如何舔出高潮| 国产熟女欧美一区二区| 成人漫画全彩无遮挡| 特大巨黑吊av在线直播| 大话2 男鬼变身卡| 一二三四中文在线观看免费高清| 免费人成在线观看视频色| 蜜臀久久99精品久久宅男| 午夜免费鲁丝| 亚洲不卡免费看| 中文字幕人妻丝袜制服| 一级毛片久久久久久久久女| av在线老鸭窝| 只有这里有精品99| 丰满少妇做爰视频| 80岁老熟妇乱子伦牲交| 99热网站在线观看| 国产无遮挡羞羞视频在线观看| 国产黄片美女视频| 久久久精品94久久精品| 精品卡一卡二卡四卡免费| 一区二区av电影网| 国模一区二区三区四区视频| 男人舔奶头视频| 纯流量卡能插随身wifi吗| 啦啦啦视频在线资源免费观看| 99热国产这里只有精品6| 国产高清不卡午夜福利| a级毛色黄片| 亚洲成人手机| 亚洲情色 制服丝袜| 国产国拍精品亚洲av在线观看| 免费播放大片免费观看视频在线观看| 亚洲成人av在线免费| 日本91视频免费播放| 在线观看免费视频网站a站| 精品人妻偷拍中文字幕| 久久精品久久精品一区二区三区| 亚洲精品一区蜜桃| 男的添女的下面高潮视频| 成人影院久久| 国产乱人偷精品视频| 中文字幕人妻丝袜制服| 亚洲情色 制服丝袜| av线在线观看网站| 91久久精品电影网| 中文字幕av电影在线播放| 新久久久久国产一级毛片| 我的老师免费观看完整版| 亚洲欧洲日产国产| 成年av动漫网址| 亚洲第一区二区三区不卡| 欧美另类一区| 亚洲美女搞黄在线观看| 亚洲精品亚洲一区二区| 精品卡一卡二卡四卡免费| 欧美97在线视频| av在线观看视频网站免费| 久久97久久精品| av线在线观看网站| 午夜福利在线观看免费完整高清在| 2018国产大陆天天弄谢| 晚上一个人看的免费电影| 欧美成人午夜免费资源| 国产淫片久久久久久久久| 久久久久人妻精品一区果冻| 亚洲av成人精品一二三区| 欧美日韩一区二区视频在线观看视频在线| 久久精品国产亚洲av天美| 免费大片18禁| 最近中文字幕2019免费版| av线在线观看网站| 久久青草综合色| 国产精品福利在线免费观看| 在线天堂最新版资源| 国产伦精品一区二区三区四那| 精品亚洲成a人片在线观看| 99九九在线精品视频 | 精品人妻熟女av久视频| 偷拍熟女少妇极品色| 久久国产乱子免费精品| 少妇高潮的动态图| 18禁动态无遮挡网站| 春色校园在线视频观看| 有码 亚洲区| 九色成人免费人妻av| 亚洲精品国产色婷婷电影| 天美传媒精品一区二区| 日日摸夜夜添夜夜爱| 2018国产大陆天天弄谢| 久久久午夜欧美精品| 国产黄片美女视频| 午夜免费观看性视频| 国产免费福利视频在线观看| 国精品久久久久久国模美| 国产免费又黄又爽又色| 欧美精品一区二区免费开放| 下体分泌物呈黄色| 成人亚洲欧美一区二区av| 欧美另类一区| 亚洲第一av免费看| 成人国产麻豆网| 亚洲欧美中文字幕日韩二区| 夫妻性生交免费视频一级片| 人人澡人人妻人| 乱码一卡2卡4卡精品| 在线观看免费日韩欧美大片 | 在线亚洲精品国产二区图片欧美 | 午夜久久久在线观看| 亚洲av综合色区一区| 夜夜看夜夜爽夜夜摸| 丰满迷人的少妇在线观看| 精品亚洲成国产av| 久久久久精品久久久久真实原创| 一级黄片播放器| 欧美区成人在线视频| av天堂久久9| 免费少妇av软件| 一区在线观看完整版| 欧美精品亚洲一区二区| 全区人妻精品视频| 99九九在线精品视频 | 嫩草影院入口| 色5月婷婷丁香| 97超视频在线观看视频| 亚洲精品国产av蜜桃| 天天躁夜夜躁狠狠久久av| 能在线免费看毛片的网站| 尾随美女入室| 97在线视频观看| 久久久a久久爽久久v久久| 亚洲熟女精品中文字幕| 日本黄色片子视频| 最新的欧美精品一区二区| 热99国产精品久久久久久7| 国产一区二区三区综合在线观看 | 国产精品一区二区性色av| 国产免费视频播放在线视频| 亚洲久久久国产精品| 国产欧美另类精品又又久久亚洲欧美| 黑人巨大精品欧美一区二区蜜桃 | 交换朋友夫妻互换小说| 婷婷色综合www| 国产一区二区三区综合在线观看 | 中文精品一卡2卡3卡4更新| 51国产日韩欧美| 精品国产一区二区三区久久久樱花| 黄色怎么调成土黄色| 中文字幕久久专区| 亚洲国产欧美在线一区| 国产精品麻豆人妻色哟哟久久| 又大又黄又爽视频免费| 国产精品久久久久久久久免| 免费高清在线观看视频在线观看| 国产美女午夜福利| 国产成人精品无人区| 午夜老司机福利剧场| 内射极品少妇av片p| 精品久久国产蜜桃| 国产免费又黄又爽又色| 大码成人一级视频| 观看免费一级毛片| 最近最新中文字幕免费大全7| 晚上一个人看的免费电影| 国产精品国产三级专区第一集| 欧美 亚洲 国产 日韩一| 亚洲中文av在线| 在线观看国产h片| 欧美97在线视频| 色吧在线观看| 亚洲精品一区蜜桃| 极品人妻少妇av视频| 免费看日本二区| 亚洲av国产av综合av卡| 久久人人爽人人片av| 亚洲av日韩在线播放| 国语对白做爰xxxⅹ性视频网站| 国产一区二区三区综合在线观看 | 国产亚洲午夜精品一区二区久久| 夜夜爽夜夜爽视频| 国产精品人妻久久久影院| 丝瓜视频免费看黄片| 精品久久久久久久久亚洲| 人妻系列 视频| 久久久亚洲精品成人影院| 国产精品久久久久成人av| 国产白丝娇喘喷水9色精品| 在线观看美女被高潮喷水网站| kizo精华| av天堂久久9| 少妇被粗大的猛进出69影院 | 亚洲精品视频女| 一本一本综合久久| 国产免费一级a男人的天堂| 午夜老司机福利剧场| h日本视频在线播放| 亚洲精品久久午夜乱码| 国产毛片在线视频| 丰满迷人的少妇在线观看| 国产又色又爽无遮挡免| 91久久精品国产一区二区三区| 精品一品国产午夜福利视频| 自拍欧美九色日韩亚洲蝌蚪91 | 伦理电影免费视频| 午夜免费男女啪啪视频观看| av一本久久久久| 久久狼人影院| 国产精品免费大片| 赤兔流量卡办理| 一边亲一边摸免费视频| 美女中出高潮动态图| 免费不卡的大黄色大毛片视频在线观看| 亚洲图色成人| 国产一区二区在线观看av| 中文字幕制服av| 午夜av观看不卡| 三级国产精品片| 日产精品乱码卡一卡2卡三| 亚洲欧美一区二区三区国产| 汤姆久久久久久久影院中文字幕| 国产精品国产三级国产专区5o| 国产精品一区二区性色av| 一本大道久久a久久精品| 人人妻人人澡人人爽人人夜夜| 免费观看无遮挡的男女| 黄色配什么色好看| 国产成人免费无遮挡视频| 国产一区二区三区综合在线观看 | 丰满乱子伦码专区| 日韩成人伦理影院| 天堂8中文在线网| 色5月婷婷丁香| 我要看日韩黄色一级片| 欧美人与善性xxx| 人人澡人人妻人| 亚洲,欧美,日韩| 国产精品女同一区二区软件| 国产在线视频一区二区| 中文字幕av电影在线播放| 在线精品无人区一区二区三| 高清不卡的av网站| 少妇高潮的动态图| 性色avwww在线观看| 成人亚洲欧美一区二区av| 性色avwww在线观看| 欧美日韩一区二区视频在线观看视频在线| 五月天丁香电影| 永久免费av网站大全| 女性被躁到高潮视频| 肉色欧美久久久久久久蜜桃| 日韩亚洲欧美综合| 国内少妇人妻偷人精品xxx网站| 久久久久精品性色| 蜜臀久久99精品久久宅男| 久久99热6这里只有精品| 国产精品无大码| 超碰97精品在线观看| 国产一区二区三区av在线| 欧美 日韩 精品 国产| 观看av在线不卡| 欧美日韩亚洲高清精品| 久久久久久伊人网av| av黄色大香蕉| 五月玫瑰六月丁香| 久久毛片免费看一区二区三区| av一本久久久久| 午夜福利,免费看| 18禁在线播放成人免费| 国产一区二区三区综合在线观看 | 国产精品国产三级国产av玫瑰| 韩国高清视频一区二区三区| 亚洲精品日韩av片在线观看| 亚洲成人一二三区av| 欧美精品人与动牲交sv欧美| av不卡在线播放| 国产精品一区二区性色av| 最近2019中文字幕mv第一页| 汤姆久久久久久久影院中文字幕| 欧美成人精品欧美一级黄| 欧美bdsm另类| 久久久a久久爽久久v久久| 大片电影免费在线观看免费| 熟女av电影| 视频中文字幕在线观看| av福利片在线观看| 亚洲av免费高清在线观看| 在线观看av片永久免费下载| 亚洲av综合色区一区| 免费播放大片免费观看视频在线观看| 丝袜喷水一区| 久久久精品94久久精品| 国产精品99久久久久久久久| 久久 成人 亚洲| 插逼视频在线观看| 日韩 亚洲 欧美在线| 亚洲欧美精品专区久久| 久久久久久久久久久丰满| 五月玫瑰六月丁香| 日韩成人伦理影院| 亚洲成人手机| 欧美精品一区二区大全| 婷婷色av中文字幕| 九九久久精品国产亚洲av麻豆| 亚洲国产毛片av蜜桃av| 一级毛片电影观看| 在线精品无人区一区二区三| 国产免费一区二区三区四区乱码| 国产成人精品久久久久久| 亚洲精品第二区| 日本av免费视频播放| 国产91av在线免费观看| 黑人高潮一二区| a级一级毛片免费在线观看| 青春草亚洲视频在线观看| 久久精品国产亚洲av涩爱| 一区二区三区免费毛片| 国产成人精品无人区| 麻豆精品久久久久久蜜桃| av免费观看日本| 亚洲精品乱码久久久v下载方式| 亚洲精品色激情综合| 午夜福利影视在线免费观看| 51国产日韩欧美| a级片在线免费高清观看视频| 亚洲精品日韩av片在线观看| 一二三四中文在线观看免费高清| 内地一区二区视频在线| 人人妻人人澡人人爽人人夜夜| 久久久久视频综合| 亚洲美女搞黄在线观看| 亚洲欧美精品专区久久| 国产无遮挡羞羞视频在线观看| 亚洲精品亚洲一区二区| 熟女人妻精品中文字幕| 日本免费在线观看一区| 日韩欧美精品免费久久| 黄色日韩在线| 在线观看国产h片| 久久国产精品男人的天堂亚洲 | 婷婷色综合www| 少妇的逼水好多| 国内少妇人妻偷人精品xxx网站| 最后的刺客免费高清国语| 久久精品熟女亚洲av麻豆精品| 亚洲精品自拍成人| 日韩免费高清中文字幕av| 久久久久久久久久人人人人人人| 国产成人免费无遮挡视频| 观看av在线不卡| 亚洲,一卡二卡三卡| 久久国产精品男人的天堂亚洲 | 亚洲不卡免费看| 亚洲欧美精品专区久久| 午夜激情福利司机影院| 亚洲av日韩在线播放| 老司机亚洲免费影院| 国产精品一区二区性色av| 麻豆成人午夜福利视频| 22中文网久久字幕| 最近2019中文字幕mv第一页| 一区二区三区精品91| 免费看日本二区| 成人国产麻豆网| 精品亚洲成国产av| 一级片'在线观看视频| 午夜免费男女啪啪视频观看| 综合色丁香网| 国产精品一区二区在线观看99| 最后的刺客免费高清国语| 精品久久久精品久久久| 亚洲电影在线观看av| 91精品伊人久久大香线蕉| 午夜免费男女啪啪视频观看| 亚洲国产精品一区三区| 久久精品国产亚洲av涩爱| 亚洲,一卡二卡三卡| 久久久久久久久久久免费av| 我要看黄色一级片免费的| 日本黄大片高清| 久久久久精品性色| 久久99精品国语久久久| 亚洲精品中文字幕在线视频 | 最近手机中文字幕大全| 99视频精品全部免费 在线| 亚洲国产精品成人久久小说| av在线观看视频网站免费| 纯流量卡能插随身wifi吗| 丰满人妻一区二区三区视频av| 国产成人精品一,二区| 噜噜噜噜噜久久久久久91| 岛国毛片在线播放| 高清欧美精品videossex| 久久免费观看电影| 免费在线观看成人毛片| 国产极品粉嫩免费观看在线 | 亚洲av在线观看美女高潮| 又粗又硬又长又爽又黄的视频| 在线观看免费视频网站a站| 青青草视频在线视频观看| 女人久久www免费人成看片| 人妻人人澡人人爽人人| 交换朋友夫妻互换小说| 精品少妇久久久久久888优播| 久久久久久久久久久免费av| 亚洲va在线va天堂va国产| 成人毛片60女人毛片免费| 日本黄色片子视频| 亚州av有码| 亚洲国产精品999| 99久久精品国产国产毛片| 久久国产乱子免费精品| 久久精品熟女亚洲av麻豆精品| 91久久精品国产一区二区成人| 久久热精品热| 亚洲精品中文字幕在线视频 | 女的被弄到高潮叫床怎么办| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲精品日本国产第一区| 日韩伦理黄色片| 中文字幕久久专区| 91精品国产国语对白视频| 成年美女黄网站色视频大全免费 | 久久久久久久久久成人| 2021少妇久久久久久久久久久| 国产午夜精品久久久久久一区二区三区| 在线观看人妻少妇| 一级二级三级毛片免费看| 国产有黄有色有爽视频| 建设人人有责人人尽责人人享有的| 亚洲av男天堂| 亚洲性久久影院| 一本一本综合久久| 亚洲精品,欧美精品| 少妇精品久久久久久久| 亚洲欧美日韩卡通动漫| 日韩人妻高清精品专区| videos熟女内射| 青青草视频在线视频观看| 成人毛片60女人毛片免费| 亚洲精品乱码久久久v下载方式| 又粗又硬又长又爽又黄的视频| 成年人免费黄色播放视频 | 美女脱内裤让男人舔精品视频| 国产视频内射| 国产精品人妻久久久影院| 五月开心婷婷网| 最近的中文字幕免费完整| 国产成人91sexporn| 国产毛片在线视频| 五月天丁香电影| 99久久精品热视频| 国产精品99久久久久久久久| 人妻少妇偷人精品九色| 桃花免费在线播放| 高清视频免费观看一区二区| 下体分泌物呈黄色| 日日撸夜夜添| 国产欧美亚洲国产| 久久精品国产亚洲网站| 九九久久精品国产亚洲av麻豆| 高清视频免费观看一区二区| 人妻少妇偷人精品九色| av免费在线看不卡| 黄色一级大片看看| 国产精品嫩草影院av在线观看| 少妇被粗大猛烈的视频| 国产av一区二区精品久久| 久久国内精品自在自线图片| 哪个播放器可以免费观看大片| 亚洲人成网站在线观看播放|