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

    東祁連山北緣斷裂帶基于深度學習的密集臺陣地震事件快速檢測與定位研究1

    2022-06-01 08:45:46楊少博王炳文張海江
    震災防御技術 2022年1期
    關鍵詞:檢測

    楊少博 王炳文 高 級 張海江

    1)中國科學技術大學地球和空間科學學院, 合肥 230026

    2)安徽蒙城地球物理國家野外科學觀測研究站, 合肥 230026

    引言

    地震事件檢測與定位是很多地震學研究的基礎。對于活動斷層區(qū)域,可利用完備、可靠的地震事件檢測與定位結果更好地了解區(qū)域地震活動性、刻畫斷層幾何形態(tài)(Ross 等,2019a,2020)。近年來,密集臺陣的大規(guī)模布設有利于更多的微地震活動監(jiān)測(Meng 等,2018;Li 等,2018),使需處理的地震數(shù)量迅速增加。人工智能在地震數(shù)據(jù)處理方面的成功應用提供了一系列快速、可靠的自動化地震數(shù)據(jù)處理方法,包括地震事件檢測(Perol 等,2018;Ross 等,2018;Yang 等,2021)、到時拾取(Zhu 等,2019;Mousavi 等,2020;Xiao 等,2021)、震相關聯(lián)(McBrearty 等,2019;Ross 等,2019b)等,這些自動化方法能夠達到甚至超越人工處理的精度,大幅度提高了數(shù)據(jù)處理效率。

    Wang 等(2020)在Raton Basin(位于美國新墨西哥州和科羅拉多州的邊界)布設了密集臺陣,研究該區(qū)域誘發(fā)地震活動性。該研究將連續(xù)數(shù)據(jù)劃分為多個12 s 的數(shù)據(jù)段,使用PhaseNet(Zhu 等,2019)拾取每個數(shù)據(jù)段中的P、S 波到時,使用REAL 算法(Zhang 等,2019)進行震相關聯(lián)并初定位,使用VELEST 算法(Kissling 等,1994)優(yōu)化定位結果,使用相對定位算法hypoDD(Waldhauser 等,2000)進行重定位。相比該地區(qū)之前的地震目錄,該流程定位了更多的地震,更好地刻畫了斷層的幾何形態(tài)。Liu 等(2020)基于PhaseNet 到時構建了2019 年加州Ridgecrest 地震序列目錄。趙明等(2021b)對PhaseNet 做遷移學習,進而用于構建長寧地震的前震目錄。Ross 等(2018)使用基于深度神經(jīng)網(wǎng)絡的Generalized Phase Detection 算法進行地震事件檢測及拾取,使用PhaseLink(Ross 等,2019b)進行震相關聯(lián),使用傳統(tǒng)定位方法獲得該區(qū)域的地震目錄。上述研究結果均表明將人工智能地震數(shù)據(jù)處理方法與傳統(tǒng)定位方法相結合,能夠快速、可靠地獲得某地區(qū)的地震目錄。

    將上述地震事件檢測與定位流程應用于東祁連山斷裂帶的密集臺陣數(shù)據(jù)處理。由于密集臺陣數(shù)據(jù)量大,在缺少GPU 加速的情況下,使用PhaseNet 拾取連續(xù)數(shù)據(jù)需耗費大量時間。因此,首先使用一個輕量的神經(jīng)網(wǎng)絡CNNDetector(Yang 等,2021)進行地震事件檢測,然后使用PhaseNet 對檢測到的事件進行到時拾取,最后使用傳統(tǒng)方法進行地震定位和震級計算。相比中國地震臺網(wǎng)正式地震目錄,獲得了更多的地震數(shù)據(jù)。

    1 研究區(qū)域

    祁連山位于青藏高原塊體、塔里木塊體和阿拉善塊體的銜接帶上,其構造形變受青藏高原擠壓作用及周緣活動斷裂的控制。始新世末-漸新世早期,伴隨著印度板塊和歐亞板塊的碰撞,祁連山地區(qū)發(fā)生構造變形與隆升;漸新世晚期-中新世早期,祁連山地區(qū)構造穩(wěn)定,地形起伏變小,形成大面積夷平面;中新世中晚期(10~17 Ma)約為斷裂帶起始活動時間,斷裂帶以 NNE-NE 向的擠壓逆沖和地殼縮短增厚為變形特征,造成斷裂帶沿線山脈的快速隆升;中更新世以來,斷裂帶轉變?yōu)樽笮呋瑸橹骷婢吣鏇_運動,調整著青藏高原東北緣地殼物質向東運動,最終被斷裂帶尾端的六盤山逆沖褶皺帶吸收(張培震等,2004;戚幫申,2014;郭鵬,2019)。

    晚第四紀以來,在印度板塊向歐亞板塊NNE 向持續(xù)的構造推擠過程中,青藏高原東北緣受北側阿拉善地塊和東側鄂爾多斯塊體的阻擋,在NE 向或NEE 向擠壓作用下,構造變形十分強烈。祁連山吸收了青藏高原NE 向的擠壓變形,張培震等(2004)計算得到青藏高原北部邊界的阿爾金山、祁連山共吸收了青藏15%~17% 的地殼縮短;祁連-海原斷裂帶作為青藏高原東北緣主要的邊界左旋走滑斷層,在調節(jié)高原東北緣地殼物質相對于阿拉善地塊的向東運動中起重要作用。

    本研究區(qū)位于北祁連山東部,主要包含皇城-雙塔斷裂、冷龍嶺斷裂。冷龍嶺斷裂與海原斷裂、老虎山斷裂、毛毛山斷裂、金強河斷裂和托萊山斷裂構成長約 1 000 km 的祁連-海原斷裂帶。冷龍嶺斷裂早期活動以擠壓逆沖為主,自中更新世以來,斷裂活動性質以左旋走滑為主,兼有正斷層分量,晚第四紀時期主要表現(xiàn)為左旋走滑運動,現(xiàn)今滑動速率約為4~6 mm/yr (郭鵬等,2017;郭鵬,2019)。

    從幾何分段上可將皇城-雙塔斷裂帶劃分為西、中、東段,其中西段為皇城盆地斷裂段,在晚更新世以來活動性弱,中、東段分別為上寺斷裂段和冬青頂斷裂段,活動性均較強。上寺斷裂段以疊瓦狀逆沖斷構造為特征,冬青頂斷裂段分為南、北枝,南枝為正斷層,北枝為逆斷層(侯康明,1998),1927 年古浪地震發(fā)生于東段。

    2 研究數(shù)據(jù)

    為確定研究區(qū)域精細的斷裂結構,并監(jiān)測祁連山斷裂帶附近的微地震活動性,于2019-07-20-2019-08-21 在祁連山斷裂帶布設了面狀密集臺陣。該臺陣包含6 條測線,每條測線由40 個短周期儀器EPS-2-M6Q組成(圖1),儀器頻帶范圍0.2~150 Hz,采樣率200 Hz。測線之間的距離約10 km,測線內臺間距從測線兩端至中間逐漸由數(shù)公里縮小至數(shù)百米,整個臺陣覆蓋范圍約為3 000 km2(50 km×60 km)。對采集數(shù)據(jù)進行預處理,包括去均值、去線性化及重采樣至100 Hz。圖1 中藍色三角表示臺站,四個白色三角形為編號1 001、1 002、1 030、2016 的臺站,黑色圓點表示歷史地震事件(莘海亮,2020),黑色虛線表示斷裂帶,HSF 表示皇城-雙塔斷裂,LLLF 表示冷龍嶺斷裂,兩個紅色五角星為編號317 和編號623 的地震事件。

    圖1 研究區(qū)域臺站及歷史地震事件分布Fig. 1 Distribution of stations and historical earthquakes in the study area

    3 數(shù)據(jù)處理流程和結果

    參考Wang 等(2020)的數(shù)據(jù)處理流程,若直接使用PhaseNet 進行事件檢測與到時拾取,需將24 h 連續(xù)數(shù)據(jù)切割成5 760 個15 s 的數(shù)據(jù)片段,在沒有GPU 的服務器上處理240 個臺站24 h 的數(shù)據(jù)需耗費約36 h。相比PhaseNet,CNNDetector 僅用于檢測地震事件,其網(wǎng)絡結構較簡單,檢測效率更高。在同臺服務器中,處理240 個臺站24 h 的數(shù)據(jù)僅需耗費約4 h,因此首先使用CNNDetector 檢測地震事件,然后使用PhaseNet 拾取到時,再使用REAL 初定位,最后使用hypoDD 進行重定位。

    3.1 地震事件檢測

    Yang 等(2021)提出基于卷積神經(jīng)網(wǎng)絡的地震事件檢測算法(CNNDetector),此算法綜合多個臺站的三分量地震記錄判斷1 個時窗內是否存在地震事件,其檢測精度高、誤檢率小,并在多個區(qū)域測試中表現(xiàn)出較好的泛化能力。進行祁連山密集臺陣數(shù)據(jù)處理時,為檢測到更多的微小地震,將240 個臺站按區(qū)域分為12 組,對數(shù)據(jù)進行3~30 Hz 的帶通濾波后,使用CNNDetector 對每組臺站分別進行檢測。設置閾值時,首先需滿足能夠檢測到大多數(shù)已知地震事件的要求,然后嘗試降低閾值,直至出現(xiàn)虛假檢測結果,由此將閾值設置為0.5。在同一時窗內,若有2 組以上的CNNDetector 檢測值超過閾值,則認為其為地震事件。在此標準下,共檢測到1 036 個地震事件,其中2 個低閾值的地震事件如圖2 所示,其地震位置由圖1 紅色五角星表示。

    圖2 CNNDetector 檢測出的2 個地震事件波形Fig. 2 Waveforms of two seismic events detected by CNNDetector

    3.2 到時拾取

    由于臺站數(shù)量多,手動拾取檢測到的1 036 個地震事件需耗費大量時間,因此使用自動到時拾取算法進行拾取。PhaseNet 是目前應用較多的基于深度學習的到時自動拾取算法(Zhu 等,2019),其訓練數(shù)據(jù)是北加州地震數(shù)據(jù)中心超過30 年記錄到的約700 000 個地震波形,其對低信噪比的數(shù)據(jù)也具有一定拾取精度,且在不同地區(qū)具有較好的泛化能力(趙明等,2021a)。將P 波和S 波拾取閾值均設為0.3,共拾取到74 585個P 波到時和75 191 個S 波到時。地震事件317 的部分到時拾取結果如圖3 所示(圖1 中白色三角形為對應的接收臺站位置),由圖3 可知,部分信噪比較低的信號也能被準確拾取。

    圖3 PhaseNet 到時拾取結果和不同臺站接收到的事件317 的波形圖Fig. 3 Examples of PhaseNet arrival time picking results and waveforms of event 317 received at different stations

    3.3 震相關聯(lián)及初定位

    手動拾取到時時會考慮臺站間信號的關聯(lián)性,從而判斷待拾取的是地震信號還是干擾信號。但自動拾取不會考慮,因此在震相關聯(lián)前的時距曲線中有很多偏離P、S 波時距曲線的點。因此對于自動拾取的震相到時,做進一步的震相關聯(lián)是十分必要的。這里我們使用基于網(wǎng)格搜索的REAL 算法(Zhang 等,2019)進行震相關聯(lián)及初定位。

    選取臺陣中心點處的USTClitho 2.0(Han 等,2022)速度模型作為震相關聯(lián)及初定位使用的速度模型(圖4)。REAL 的搜索范圍為0.6°×0.6°×35 km,搜索網(wǎng)格尺寸為0.05°×0.05°×2 km,最終保留了550 個滿足以下條件的地震事件(圖5(a)):超過10 個P 波到時、5 個S 波到時、20 個P 波加S 波到時且臺站間隙角<240°,其中震相走時殘差需<0.5 s。震相關聯(lián)后剩余43 152 個P 波到時及42 414 個S 波到時。震相關聯(lián)后的時距曲線如圖5(b)所示,異常點已剔除,其中灰色圓點為關聯(lián)前的震相散點,紅點為關聯(lián)后剩余的震相散點。初定位后的殘差分布如圖5(c)所示。

    圖4 地震定位使用的一維速度模型Fig. 4 1D velocity models used for earthquake location

    圖5 REAL 震相關聯(lián)及初定位結果Fig. 5 Seismic phase association and earthquake location results by the REAL algorithm

    3.4 重定位

    為進一步提高地震定位精度,使用相對定位算法hypoDD(Waldhauser 等,2000)進行重定位。hypo-DD 算法使用地震對的相對到時數(shù)據(jù)進行地震定位,由于震中距遠大于地震對的距離,2 個地震近似擁有共同路徑,因此可減小速度模型不準確帶來的相對定位誤差,且可更好地約束地震的相對位置。重定位時使用USTClitho2.0 速度模型,選取的地震對距離<10 km,相互關聯(lián)的地震數(shù)目超過8 個。重定位后剩余517 個地震事件,地震分布如圖6 中黑色圓點所示,其中灰色圓點為2013-2016 年重定位的歷史地震分布(莘海亮,2020),灰色虛線為斷裂帶。重定位的迭代過程及殘差分布如圖7 所示。

    圖6 HypoDD 定位結果Fig. 6 HypoDD location results

    圖7 HypoDD 迭代過程及走時殘差分布Fig. 7 HypoDD iterative process and residual distribution

    3.5 震級計算

    參考Wang 等(2020)研究方法計算震級:

    式中,ML為 震級,A為振幅,d為震中距。

    對事件波形數(shù)據(jù)去儀器響應、濾波至1~50 Hz 后,對于每個地震,在所有存在到時的臺站上計算震級,取所有臺站計算震級的中位數(shù)為最終的震級,如圖8(a)所示,2 個地震目錄中共同地震事件的震級對比如圖8(b)所示,由此可知,2 種震級差異較小,說明式(1)計算的震級是可靠的。

    圖8 震級分布與對比Fig. 8 Magnitude distribution and comparison

    4 討論與結論

    在祁連山密集臺陣觀測期間,中國地震臺網(wǎng)正式地震目錄中共有39 個位于臺陣內的地震事件,而密集臺陣共監(jiān)測到550 個地震事件,約是地震臺網(wǎng)地震數(shù)量的15 倍,尤其監(jiān)測到了更多小于0 級的地震事件。通過對比發(fā)現(xiàn),中國地震臺網(wǎng)正式地震目錄中的39 個地震事件全部被密集臺陣監(jiān)測。因此通過密集臺陣監(jiān)測,可大大提高微地震監(jiān)測能力,可監(jiān)測到-2 級的微震事件。

    由于數(shù)據(jù)量大,為減少人工拾取工作量,提高數(shù)據(jù)處理效率,使用最新發(fā)展的基于人工智能的地震數(shù)據(jù)處理方法。在缺少GPU 加速的情況下,首先使用CNNDetector 進行地震事件檢測,然后使用PhaseNet 拾取P 波和S 波到時,大大節(jié)約了計算時間。在人工智能拾取到時的基礎上,結合傳統(tǒng)定位方法REAL 及hypo-DD,精確地定位出517 個地震,這表明基于深度學習進行地震事件檢測與到時拾取,并在此基礎上進行地震定位是可靠、高效的。密集臺陣重定位的地震位置與研究區(qū)域歷史地震活動具有較好的一致性,均由西至東呈現(xiàn)出4 個條帶狀分布。最西側的地震簇和冷龍嶺斷裂相關,走向與斷裂走向類似,呈近NW-SE 向。與地表斷層跡線相比,最西側的地震簇位于斷層東北側,垂向剖面顯示該斷層向EN 向傾斜。與最西側的地震簇相比,東側的3 個地震條帶走向和皇城-雙塔斷裂、冷龍嶺斷裂走向不一致,而是呈NNW-NNE 向,與最東側的斷層走向一致。地震定位結果表明,研究區(qū)域除地表顯現(xiàn)的NW-SE/NEE-SWW 向的斷層外,還存在活躍的NNW-NNE 向隱伏活躍斷層。這些隱伏活躍斷層與地表顯現(xiàn)的斷層呈一定斜交關系,且更活躍。與中國地震臺網(wǎng)地震目錄給出的地震定位相比,基于密集臺陣給出的地震定位表現(xiàn)出了更明顯的線性結構,體現(xiàn)了地震定位的高精度,可刻畫出斷層面的三維幾何形態(tài)。

    致謝 感謝國家重點研發(fā)計劃(2018YFC1504102)的資助。感謝中國地震臺網(wǎng)中心國家地震科學數(shù)據(jù)中心(http://data.earthquake.cn)提供的數(shù)據(jù)支撐,感謝莘海亮博士提供的歷史地震目錄。感謝審稿人提供的寶貴建議。

    猜你喜歡
    檢測
    QC 檢測
    “不等式”檢測題
    “一元一次不等式”檢測題
    “一元一次不等式組”檢測題
    “幾何圖形”檢測題
    “角”檢測題
    “有理數(shù)的乘除法”檢測題
    “有理數(shù)”檢測題
    “角”檢測題
    “幾何圖形”檢測題
    中文字幕av成人在线电影| 日本av免费视频播放| 午夜福利网站1000一区二区三区| 九九在线视频观看精品| 内射极品少妇av片p| 黄色一级大片看看| 精品人妻一区二区三区麻豆| 啦啦啦啦在线视频资源| 久久久久性生活片| 国产精品久久久久久久久免| 青春草国产在线视频| 在现免费观看毛片| videossex国产| 丰满少妇做爰视频| 一级二级三级毛片免费看| 大码成人一级视频| 男女下面进入的视频免费午夜| 日韩欧美 国产精品| 人妻夜夜爽99麻豆av| 大话2 男鬼变身卡| 交换朋友夫妻互换小说| 亚洲av欧美aⅴ国产| 婷婷色麻豆天堂久久| 亚洲熟女精品中文字幕| 日本色播在线视频| 亚洲欧洲国产日韩| 国产精品人妻久久久影院| 免费久久久久久久精品成人欧美视频 | 亚洲精品亚洲一区二区| 青春草亚洲视频在线观看| 最近2019中文字幕mv第一页| 欧美日本视频| 91aial.com中文字幕在线观看| 欧美少妇被猛烈插入视频| 中文乱码字字幕精品一区二区三区| 丝袜喷水一区| 我的老师免费观看完整版| 国产精品国产三级国产av玫瑰| 91精品国产国语对白视频| 久久久久久久大尺度免费视频| 我的女老师完整版在线观看| 爱豆传媒免费全集在线观看| 在线 av 中文字幕| 午夜免费观看性视频| 国产精品一区二区性色av| 18禁裸乳无遮挡动漫免费视频| 一区二区三区免费毛片| 国产一区有黄有色的免费视频| 国产熟女欧美一区二区| 日韩av免费高清视频| 男人和女人高潮做爰伦理| av在线老鸭窝| 交换朋友夫妻互换小说| 欧美一区二区亚洲| 内地一区二区视频在线| 亚洲欧美一区二区三区黑人 | 国产av国产精品国产| 欧美精品一区二区大全| 亚洲精品久久久久久婷婷小说| 在线观看免费日韩欧美大片 | 日韩视频在线欧美| 一个人看视频在线观看www免费| 亚洲第一区二区三区不卡| 国产精品国产三级国产av玫瑰| 美女cb高潮喷水在线观看| 国产一区有黄有色的免费视频| 成年美女黄网站色视频大全免费 | 亚洲美女视频黄频| 国产精品久久久久成人av| 又黄又爽又刺激的免费视频.| 国产高清不卡午夜福利| 亚洲国产精品专区欧美| 五月天丁香电影| 日韩视频在线欧美| 久久99精品国语久久久| 久久6这里有精品| 纯流量卡能插随身wifi吗| 人妻制服诱惑在线中文字幕| 最新中文字幕久久久久| 久久国内精品自在自线图片| 国产成人精品一,二区| 久久人妻熟女aⅴ| 色婷婷av一区二区三区视频| 久久影院123| 最后的刺客免费高清国语| 婷婷色av中文字幕| 男人添女人高潮全过程视频| 久久国产精品大桥未久av | 一本—道久久a久久精品蜜桃钙片| 午夜福利影视在线免费观看| 国产欧美日韩精品一区二区| 国产精品成人在线| 亚洲高清免费不卡视频| 观看免费一级毛片| 午夜福利在线在线| 久久6这里有精品| 亚洲av在线观看美女高潮| 九九在线视频观看精品| 国产精品欧美亚洲77777| 国产精品99久久99久久久不卡 | 一级毛片aaaaaa免费看小| av网站免费在线观看视频| 国产精品欧美亚洲77777| 国产乱来视频区| 欧美+日韩+精品| 久久精品久久久久久噜噜老黄| 一区二区三区免费毛片| 久久久精品免费免费高清| 久久久午夜欧美精品| 午夜免费鲁丝| 最近最新中文字幕大全电影3| 亚洲av日韩在线播放| 一边亲一边摸免费视频| 啦啦啦啦在线视频资源| 你懂的网址亚洲精品在线观看| 国产一区亚洲一区在线观看| 少妇 在线观看| 99re6热这里在线精品视频| 国产一区二区在线观看日韩| 免费观看的影片在线观看| 男女无遮挡免费网站观看| 国产永久视频网站| 欧美zozozo另类| 亚洲av在线观看美女高潮| 精品熟女少妇av免费看| 男女国产视频网站| 久久久久久人妻| 欧美区成人在线视频| 精品午夜福利在线看| 久久精品夜色国产| 成人亚洲精品一区在线观看 | 国产91av在线免费观看| 中文天堂在线官网| 久久国产亚洲av麻豆专区| 最近最新中文字幕大全电影3| 亚洲精品aⅴ在线观看| 国产精品爽爽va在线观看网站| 久热久热在线精品观看| 亚洲第一区二区三区不卡| 久久久精品94久久精品| 欧美高清性xxxxhd video| 高清午夜精品一区二区三区| 舔av片在线| 亚洲人成网站在线播| 国产av码专区亚洲av| 2018国产大陆天天弄谢| 久久99热6这里只有精品| 网址你懂的国产日韩在线| 日韩大片免费观看网站| 免费看av在线观看网站| 成人高潮视频无遮挡免费网站| 精品久久久久久电影网| 日本欧美国产在线视频| 99re6热这里在线精品视频| 久久99热这里只有精品18| av不卡在线播放| 国产又色又爽无遮挡免| 久久久成人免费电影| 精品一区二区三区视频在线| 看免费成人av毛片| 2018国产大陆天天弄谢| 国产欧美日韩一区二区三区在线 | 黄色一级大片看看| 亚州av有码| 久久人人爽人人片av| 国产一区二区三区综合在线观看 | 国产精品国产av在线观看| 精品久久久久久久久亚洲| 亚洲高清免费不卡视频| 高清在线视频一区二区三区| 一级二级三级毛片免费看| 欧美另类一区| 毛片一级片免费看久久久久| 免费看日本二区| 精品人妻一区二区三区麻豆| 精品熟女少妇av免费看| 一区二区av电影网| 丰满迷人的少妇在线观看| 亚洲国产色片| 国产欧美日韩一区二区三区在线 | 免费黄网站久久成人精品| 成人午夜精彩视频在线观看| 九九在线视频观看精品| 亚洲国产av新网站| 久久青草综合色| 少妇人妻精品综合一区二区| 日韩视频在线欧美| 99精国产麻豆久久婷婷| 国产精品av视频在线免费观看| 丰满迷人的少妇在线观看| 亚洲精品第二区| 免费少妇av软件| 亚洲成色77777| 哪个播放器可以免费观看大片| 国产精品精品国产色婷婷| 精品久久久精品久久久| 美女中出高潮动态图| 欧美日韩视频精品一区| 欧美精品亚洲一区二区| 看十八女毛片水多多多| 一级二级三级毛片免费看| videos熟女内射| 亚洲精品日韩在线中文字幕| 91精品国产国语对白视频| 在线观看免费视频网站a站| 亚洲成人中文字幕在线播放| 国产精品一区www在线观看| 亚洲最大成人中文| 久久ye,这里只有精品| 亚洲不卡免费看| 在线观看av片永久免费下载| 亚洲精品国产av成人精品| 在线播放无遮挡| 色吧在线观看| 国产欧美日韩一区二区三区在线 | 欧美区成人在线视频| 久久久久人妻精品一区果冻| av在线老鸭窝| 少妇被粗大猛烈的视频| 少妇熟女欧美另类| 交换朋友夫妻互换小说| 亚洲精品国产av成人精品| 欧美人与善性xxx| 亚洲精品乱久久久久久| 九九爱精品视频在线观看| 22中文网久久字幕| 插阴视频在线观看视频| 日韩,欧美,国产一区二区三区| 久久人人爽人人爽人人片va| 国产日韩欧美亚洲二区| 免费观看a级毛片全部| 国产精品嫩草影院av在线观看| 国产毛片在线视频| 一级毛片电影观看| 国产色婷婷99| av免费在线看不卡| 99热这里只有精品一区| 亚洲激情五月婷婷啪啪| 女性被躁到高潮视频| 亚洲中文av在线| 激情五月婷婷亚洲| 国产成人精品一,二区| 最黄视频免费看| 狂野欧美激情性bbbbbb| 波野结衣二区三区在线| 日本av手机在线免费观看| 成年女人在线观看亚洲视频| 国产成人免费观看mmmm| 午夜精品国产一区二区电影| 欧美日韩在线观看h| 热99国产精品久久久久久7| 国产精品一区二区三区四区免费观看| 夫妻午夜视频| 国产色婷婷99| 18禁裸乳无遮挡动漫免费视频| 久热久热在线精品观看| 国产精品成人在线| 久久久久国产精品人妻一区二区| 亚洲欧美一区二区三区黑人 | 中文在线观看免费www的网站| 自拍偷自拍亚洲精品老妇| 国产精品久久久久成人av| 岛国毛片在线播放| 亚洲精品久久午夜乱码| 国产精品不卡视频一区二区| 国产亚洲91精品色在线| 老师上课跳d突然被开到最大视频| 欧美精品一区二区免费开放| 国产视频内射| 男女边摸边吃奶| 国语对白做爰xxxⅹ性视频网站| 欧美+日韩+精品| 亚洲欧洲日产国产| 久久国产精品男人的天堂亚洲 | 亚洲丝袜综合中文字幕| 我的女老师完整版在线观看| 精品亚洲乱码少妇综合久久| 99热国产这里只有精品6| 少妇 在线观看| 寂寞人妻少妇视频99o| 大片免费播放器 马上看| 99re6热这里在线精品视频| 夫妻午夜视频| 又黄又爽又刺激的免费视频.| 少妇丰满av| 国产乱来视频区| 欧美精品一区二区大全| av女优亚洲男人天堂| 九九在线视频观看精品| 最后的刺客免费高清国语| 黄色配什么色好看| 成人影院久久| 夜夜爽夜夜爽视频| 亚洲av男天堂| 亚洲精品乱码久久久v下载方式| 亚洲美女黄色视频免费看| 下体分泌物呈黄色| 在线看a的网站| 免费观看无遮挡的男女| 精品人妻一区二区三区麻豆| 少妇裸体淫交视频免费看高清| 久久久精品94久久精品| 中文乱码字字幕精品一区二区三区| 日韩视频在线欧美| 国产精品福利在线免费观看| 好男人视频免费观看在线| 久久国产精品男人的天堂亚洲 | 日本vs欧美在线观看视频 | 91精品伊人久久大香线蕉| 亚洲精品中文字幕在线视频 | 综合色丁香网| 欧美日韩国产mv在线观看视频 | 久久久久久久久大av| 中文乱码字字幕精品一区二区三区| 亚洲精华国产精华液的使用体验| 哪个播放器可以免费观看大片| 18禁裸乳无遮挡免费网站照片| 日韩一区二区视频免费看| av播播在线观看一区| 国产免费一区二区三区四区乱码| 亚洲第一av免费看| 在线观看免费日韩欧美大片 | 春色校园在线视频观看| 80岁老熟妇乱子伦牲交| 亚洲精品自拍成人| 亚洲成人手机| 人妻 亚洲 视频| 亚洲欧美成人综合另类久久久| 狂野欧美白嫩少妇大欣赏| kizo精华| 免费人妻精品一区二区三区视频| 老司机影院成人| 黑丝袜美女国产一区| 91在线精品国自产拍蜜月| 在线播放无遮挡| a级毛片免费高清观看在线播放| 精品视频人人做人人爽| 午夜激情久久久久久久| 国产成人精品一,二区| 日韩大片免费观看网站| 久久久久久久亚洲中文字幕| 日本vs欧美在线观看视频 | 美女脱内裤让男人舔精品视频| 欧美zozozo另类| 国产成人freesex在线| 日韩强制内射视频| 亚洲精华国产精华液的使用体验| 99九九线精品视频在线观看视频| 伦精品一区二区三区| 男女免费视频国产| 嫩草影院入口| 久久国产精品大桥未久av | 九九爱精品视频在线观看| 国产一区二区在线观看日韩| 亚洲美女搞黄在线观看| 黄片wwwwww| 美女视频免费永久观看网站| 最黄视频免费看| 亚洲国产成人一精品久久久| 一区二区av电影网| 乱码一卡2卡4卡精品| 高清日韩中文字幕在线| 尤物成人国产欧美一区二区三区| 亚洲av男天堂| 亚洲第一av免费看| 色哟哟·www| 午夜免费观看性视频| 国产精品三级大全| 少妇的逼水好多| av一本久久久久| 免费看不卡的av| 麻豆国产97在线/欧美| 大又大粗又爽又黄少妇毛片口| 久久久成人免费电影| 精品亚洲乱码少妇综合久久| 各种免费的搞黄视频| 国产日韩欧美在线精品| 亚洲国产欧美人成| 国产伦理片在线播放av一区| 91在线精品国自产拍蜜月| 国产欧美亚洲国产| 18禁在线无遮挡免费观看视频| 性高湖久久久久久久久免费观看| 美女内射精品一级片tv| 精品熟女少妇av免费看| 丝袜脚勾引网站| 日韩一本色道免费dvd| 免费黄色在线免费观看| 街头女战士在线观看网站| 干丝袜人妻中文字幕| 黄片无遮挡物在线观看| 国产欧美日韩一区二区三区在线 | 汤姆久久久久久久影院中文字幕| 婷婷色av中文字幕| 搡女人真爽免费视频火全软件| 亚洲精品一区蜜桃| 免费不卡的大黄色大毛片视频在线观看| 精品亚洲乱码少妇综合久久| 中国三级夫妇交换| 亚洲三级黄色毛片| 成人特级av手机在线观看| 亚洲精华国产精华液的使用体验| 亚洲精品久久午夜乱码| 国产精品一区二区在线观看99| 亚洲国产最新在线播放| 1000部很黄的大片| 三级国产精品片| 国产精品成人在线| 九九爱精品视频在线观看| 亚洲综合色惰| 欧美一级a爱片免费观看看| 有码 亚洲区| kizo精华| 啦啦啦视频在线资源免费观看| 女性生殖器流出的白浆| 中文字幕精品免费在线观看视频 | 妹子高潮喷水视频| 日韩伦理黄色片| 成年av动漫网址| 亚洲国产av新网站| 成人综合一区亚洲| 91久久精品国产一区二区三区| 亚洲国产精品999| 人人妻人人添人人爽欧美一区卜 | 偷拍熟女少妇极品色| 亚洲国产日韩一区二区| 国语对白做爰xxxⅹ性视频网站| 久久人妻熟女aⅴ| 国产人妻一区二区三区在| 日韩av在线免费看完整版不卡| 成年免费大片在线观看| 26uuu在线亚洲综合色| av线在线观看网站| 国产亚洲欧美精品永久| 欧美老熟妇乱子伦牲交| 天美传媒精品一区二区| 国产男女超爽视频在线观看| 国产精品熟女久久久久浪| 久久久亚洲精品成人影院| 色视频www国产| 精品人妻视频免费看| 国产爽快片一区二区三区| 国产黄色免费在线视频| 51国产日韩欧美| 日产精品乱码卡一卡2卡三| 免费看av在线观看网站| 日本wwww免费看| tube8黄色片| 有码 亚洲区| 91狼人影院| 卡戴珊不雅视频在线播放| 高清午夜精品一区二区三区| 亚洲欧美一区二区三区国产| av播播在线观看一区| 蜜桃久久精品国产亚洲av| 亚洲av男天堂| 久久久精品免费免费高清| 在线观看美女被高潮喷水网站| 国产欧美亚洲国产| 春色校园在线视频观看| 亚洲av国产av综合av卡| 干丝袜人妻中文字幕| 日本wwww免费看| 日韩视频在线欧美| 国产精品女同一区二区软件| 免费av不卡在线播放| 毛片女人毛片| 亚洲真实伦在线观看| 成年av动漫网址| 亚洲自偷自拍三级| 欧美高清成人免费视频www| 天美传媒精品一区二区| 99视频精品全部免费 在线| 男人爽女人下面视频在线观看| 这个男人来自地球电影免费观看 | 久久99热6这里只有精品| 欧美精品一区二区免费开放| 日本黄色片子视频| av播播在线观看一区| 国产乱来视频区| 中文字幕av成人在线电影| 在现免费观看毛片| 国产高清不卡午夜福利| 欧美成人一区二区免费高清观看| 亚洲中文av在线| 男男h啪啪无遮挡| 免费大片黄手机在线观看| 久久6这里有精品| 女人久久www免费人成看片| 一级毛片久久久久久久久女| 国产又色又爽无遮挡免| 激情五月婷婷亚洲| 在线观看免费视频网站a站| 国产男人的电影天堂91| 免费黄频网站在线观看国产| 男女下面进入的视频免费午夜| 国产v大片淫在线免费观看| 日本午夜av视频| 亚洲av不卡在线观看| 黄片wwwwww| 亚洲精品乱码久久久v下载方式| 草草在线视频免费看| 亚洲图色成人| 久久久久久伊人网av| 国产亚洲av片在线观看秒播厂| 99久久中文字幕三级久久日本| 97在线人人人人妻| 亚洲av.av天堂| 国产一区有黄有色的免费视频| 黑人高潮一二区| 美女中出高潮动态图| 多毛熟女@视频| 国产精品人妻久久久久久| 亚洲国产精品一区三区| 中国美白少妇内射xxxbb| 亚洲国产精品一区三区| 亚洲,欧美,日韩| 校园人妻丝袜中文字幕| 国产欧美亚洲国产| 91精品伊人久久大香线蕉| 大陆偷拍与自拍| 日韩免费高清中文字幕av| 午夜福利在线观看免费完整高清在| 国产精品嫩草影院av在线观看| 一个人看的www免费观看视频| 国产成人a∨麻豆精品| 国产一区二区三区综合在线观看 | 嫩草影院新地址| 99国产精品免费福利视频| 日韩成人av中文字幕在线观看| 欧美人与善性xxx| 国产av国产精品国产| 三级经典国产精品| 美女福利国产在线 | 一本色道久久久久久精品综合| 国产日韩欧美在线精品| 欧美精品国产亚洲| h视频一区二区三区| 伦理电影免费视频| 国产成人一区二区在线| 天美传媒精品一区二区| 欧美激情极品国产一区二区三区 | 狠狠精品人妻久久久久久综合| 视频区图区小说| 国产日韩欧美亚洲二区| 亚洲av福利一区| 成人午夜精彩视频在线观看| 丰满迷人的少妇在线观看| 国产成人freesex在线| 亚洲高清免费不卡视频| 九九在线视频观看精品| 性色av一级| 免费不卡的大黄色大毛片视频在线观看| 高清视频免费观看一区二区| 美女cb高潮喷水在线观看| 深夜a级毛片| 男人添女人高潮全过程视频| 国产免费视频播放在线视频| 亚洲精品亚洲一区二区| videos熟女内射| 欧美日韩视频高清一区二区三区二| 午夜福利在线观看免费完整高清在| 蜜桃亚洲精品一区二区三区| 国产男女内射视频| 久久久久久久久久久免费av| 人妻系列 视频| 免费不卡的大黄色大毛片视频在线观看| 亚洲在久久综合| 精品一区二区三区视频在线| 精品国产三级普通话版| 亚洲欧美成人精品一区二区| 久久国内精品自在自线图片| 看非洲黑人一级黄片| 国产在线男女| 日韩免费高清中文字幕av| 久久国产乱子免费精品| av专区在线播放| 国模一区二区三区四区视频| 欧美精品一区二区免费开放| 美女中出高潮动态图| 国产成人a∨麻豆精品| 国产精品国产三级专区第一集| 亚洲国产毛片av蜜桃av| 亚洲欧美日韩卡通动漫| 亚洲精品亚洲一区二区| 一个人免费看片子| 久久6这里有精品| 午夜日本视频在线| 免费黄频网站在线观看国产| av在线播放精品| 国产永久视频网站| 中文字幕人妻熟人妻熟丝袜美| 18禁在线播放成人免费| 1000部很黄的大片| 国产在线一区二区三区精| 一级毛片久久久久久久久女| 91久久精品电影网| 新久久久久国产一级毛片| 精品少妇久久久久久888优播| 欧美精品人与动牲交sv欧美| 久久精品国产亚洲av天美| 美女xxoo啪啪120秒动态图| 1000部很黄的大片| 精品亚洲乱码少妇综合久久| 色视频在线一区二区三区| 中文字幕av成人在线电影| 免费黄频网站在线观看国产| 日日摸夜夜添夜夜爱| 亚洲av中文av极速乱| 欧美一级a爱片免费观看看| 大片免费播放器 马上看| 欧美日韩在线观看h| 久久久精品免费免费高清| 日本一二三区视频观看| 男女下面进入的视频免费午夜|