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

    基于臺陣資料測量高鐵信號視速度的方法

    2022-10-11 08:06:12殷常陽溫景充石永祥寧杰遠
    北京大學學報(自然科學版) 2022年5期
    關(guān)鍵詞:面波波場臺站

    殷常陽 溫景充 石永祥 寧杰遠,2,?

    基于臺陣資料測量高鐵信號視速度的方法

    殷常陽1溫景充1石永祥1寧杰遠1,2,?

    1.北京大學地球與空間科學學院, 北京 100871; 2.河北紅山地球物理國家野外科學觀測研究站, 邢臺 054000; ?通信作者, E-mail: njy@pku.edu.cn

    針對運行于高架橋上高鐵引發(fā)的震源, 基于對移動組合源形成的波場傳播的理論分析, 提出一種利用多臺互相關(guān)測量高鐵波場視速度的方法。然后, 生成分層介質(zhì)條件下的高鐵波場合成數(shù)據(jù), 計算信號視速度, 并通過與理論結(jié)果的對比, 驗證方法的正確性。

    高鐵; 移動組合源; 干涉場; 視速度

    隨著中國東部經(jīng)濟的發(fā)展, 對地下空間的利用需求越來越迫切。因此, 自動監(jiān)測路基病害, 實現(xiàn)地震預測和預警尤為重要。然而, 中國東部地震震源相對缺乏, 地震觀測資料信噪比低, 嚴重地制約了地震學手段對該區(qū)域地下結(jié)構(gòu)的探測與監(jiān)測。

    高鐵產(chǎn)生的振動是一種全新的優(yōu)質(zhì)地震震源, 能夠激發(fā)穩(wěn)定的、重復性強的地震波場。探索利用高鐵這類移動組合震源進行地下結(jié)構(gòu)探測與監(jiān)測相關(guān)的新理論、新方法和新技術(shù), 既能提供滿足覆蓋要求的進行地下結(jié)構(gòu)探測與監(jiān)測的優(yōu)質(zhì)綠色震源, 也能解決觀測資料信噪比低的問題。

    觀測結(jié)果表明, 高鐵震源不同于地震震源和爆炸源, 具有明顯的分立譜特征[1–4]。與此相呼應的是, 高鐵波場有干涉場特點, 甚至遠臺信號先于近臺信號到達臺站[5], 這與傳統(tǒng)震源的波場明顯不同。要正確地進行地下結(jié)構(gòu)探測, 必須針對高鐵地震波場特征, 改進已有方法或?qū)ふ倚路椒ā?/p>

    雖然高鐵震源反演方面的研究已取得重要進展[6–8], 但是利用高鐵地震波場觀測資料正確地反演地下結(jié)構(gòu)仍然面臨考驗[9–10], 其中的關(guān)鍵問題是如何正確地進行信號提取和應用。王曉凱等[13">[11–12]利用信號分析方法, 從高鐵地震記錄中實現(xiàn)信號提取。但是, 單臺法提取的地震信號容易受到噪聲的干擾。溫景充等[13]基于模擬記錄, 討論高鐵地震面波相速度頻散曲線的提取和修正方法, 為利用高鐵振動的面波信號研究地下結(jié)構(gòu)奠定了基礎, 但是其所用方法為雙臺法, 仍然容易受到噪聲的干擾。

    基于臺陣資料的視速度測量是傳統(tǒng)地震學的經(jīng)典方法[14–16]。該方法測量地震波在地表的傳播速度, 在信噪比高的情況下, 可以得到確定性結(jié)果。對于面波, 得到的是面波相速度和群速度, 避免了地震波傳播路徑偏離大圓弧對計算結(jié)果的影響。對于體波, 可以得到相應的視速度, 為進一步的分析和計算提供高質(zhì)量的基礎資料。

    對于在高架橋上勻速行駛的高鐵, 如果把固定的等間隔橋墩視為點源, 則每個橋墩點源的震源時間函數(shù)大致相同, 依次等間隔延遲激發(fā), 那么將在特定的頻率形成穩(wěn)定干涉場。石永祥等[17]基于理論分析發(fā)現(xiàn), 在特定的頻率, 高鐵波場在特定角度以球面波方式傳播。在臺站間距較小時, 可以將高鐵波場的傳播視為平面波傳播來處理。

    針對上述特點, 本文設計基于多臺法測量不同頻率下視速度的方法。當震相為面波時, 測量結(jié)果即為面波的相速度。最后, 基于合成數(shù)據(jù), 利用本文提出的方法計算視速度。通過與已知速度結(jié)果的對比, 驗證方法的正確性。

    1 多臺法測量視速度

    1.1 高鐵波場特征

    將原點橋墩處的震源時間函數(shù)的傅里葉變換記為(), 則坐標為g的接收點的豎直方向頻域地震記錄(g,)表示為各個橋墩產(chǎn)生波場的疊加:

    圖1 高鐵波場示意圖

    Fig. 1 Schematic diagram of high-speed rail wavefield

    1.2 多臺法測量視速度原理

    從 1.1 節(jié)得到結(jié)論: 在特定的頻率, 可以認為等間距橋墩高鐵地震波場是主要能量沿某一特定方向傳播的平面波, 描述其傳播行為的表達式為

    式(5)表示兩個臺站的互相關(guān)結(jié)果, 相位差為

    利用式(6)可以求取波數(shù), 其中

    這里有 3 個未知數(shù)k,kk, 所以需要 3 個以上臺站對進行求解, 即有

    求解式(7)就可以得到k,kk。然后, 求出的模:

    利用式(8)的結(jié)果, 即可得到波速:

    求解式(10), 得到kk, 求出的模:

    1.3 相位多值問題及解決方法

    式(14)有唯一解的條件是 3 個臺站不在同一直線上。臺站 1 與臺站 2 之間也可以計算相位差, 得到一個與式(14)線性相關(guān)的方程。實際上, 也可以利用更多臺站對的結(jié)果, 形成一個超定方程, 但這樣做不一定有利于方程的正確求解, 因為相位具有非唯一性。

    相位的非唯一性是利用互相關(guān)方法求取兩個臺站間面波相位差的時候遇到的共性問題。例如, 用傳統(tǒng)的雙臺法提取頻散曲線時, 就會得到多條相位相差 2π 的頻散曲線, 且頻率越高, 相鄰頻散曲線的相速度差就越小, 僅通過單一頻率的結(jié)果無法分辨哪個是正確的相速度值。在提取天然遠震面波或背景噪聲相速度時, 頻段是比較完整的, 可以通過低頻長周期的相速度值確定選取哪一條頻散曲線。但是, 對于高鐵地震記錄, 只能得到橋頻及其倍頻附近一段頻率范圍的結(jié)果, 因此難以通過低頻結(jié)果進行約束。針對這一問題, 可以考慮的解決方法有以下兩種。

    第一種方法是盡量減小臺間距, 基于密集臺陣進行觀測, 使得臺間距小于最小波長。這就相當于將相位差限制在 2π 以內(nèi)。但是, 這樣做會使低頻部分的誤差變大, 因為面波方法一般要求臺間距不小于波長的一半[16]。

    第二種方法是對相位進行整數(shù)個 2π 的校正, 將每組臺站對計算出的相位差加上 2π。比如, 整數(shù)取值為?4~5, 則共有 10 組不同的相位, 如果考慮式(10), 方程組中包含個臺站對, 則需要計算的方程數(shù)量為 10, 計算量極大, 且篩選困難。因此, 本文考慮使用式(14), 每次解方程只包含兩個臺站對, 利用不同的臺站對解出兩個結(jié)果, 對兩個結(jié)果進行比較后篩選正確值。作為測試, 設置 4 個臺站, 位置坐標分別為(0m, 0m)、(80m, 150m)、(120m, ?120m)和(200m, 30m)。模擬沿方向傳播的平面瑞利波, 給定的相速度在 300~600m/s 范圍內(nèi), 則能夠直接計算出臺站間的相位差。圖 2(a)和(b)分別展示用前 3 個臺站計算出來的相速度結(jié)果及用后 3個臺站算出來的相速度結(jié)果, 二者僅在模型值(虛線)處重合。利用這一特點, 就可以篩選出正確的相速度值。想要得到誤差更小的測量結(jié)果, 可以考慮利用更多的臺站組合, 每個臺站組合包括 3 個臺站。對于不同臺站組合的篩選結(jié)果, 利用不同臺站正確重合的特點取交集。

    2 合成數(shù)據(jù)結(jié)果

    2.1 合成數(shù)據(jù)的生成

    本文基于接近實際情況的橋墩間距和高鐵速度設計合成數(shù)據(jù), 高鐵運行和臺站位置模型如圖 3 所示。高鐵速度=80m/s, 橋墩間距=30m??紤]合成數(shù)據(jù)在特定時窗內(nèi)的完整性, 將鐵路長度設為9000m, 共 301 個橋墩。臺站 1 的位置坐標為(4500m,2000m), 臺站 2 的位置坐標為(4550m, 2000m), 臺站 3 的位置坐標為(4500m, 2050m), 臺站 4 的位置坐標為(4550m, 2050m)。

    高鐵線路的橋梁段主要分布在我國東部平原地區(qū)。因此, 本文針對近地表為松散沉積的情況給出介質(zhì)結(jié)構(gòu)模型(圖 4(a)): 第一層密度1=2.6kg/m3, P波速度1=800m/s, S 波速度1=350m/s, 厚度1=20 m; 第二層密度2=3.0kg/m3, P 波速度2=1000m/s, S 波速度2=450m/s, 厚度2=50m; 第三層密度3=3.3kg/m3, P 波速度3=1200m/s, S 波速度3=750m/s, 向下無限延伸。由于臺站布設在地表, 求取視速度時, 可以利用面波的相速度結(jié)果, 根據(jù)上述模型計算得到瑞利波理論頻散曲線(圖 4(b))。

    在生成模擬數(shù)據(jù)時, 本文沿用 Yao 等[16]給出的震源時間函數(shù)以及格林函數(shù)計算方法。本文只關(guān)注格林函數(shù)中的33分量, 因此只計算瑞利波的格林函數(shù), 即

    (a)前3個臺站計算結(jié)果的頻散曲線; (b)后3個臺站計算結(jié)果的頻散曲線?;疑珜嵕€為計算的相速度值, 黑色虛線為參考值

    虛線表示高速鐵路, 由橋墩架起, 橋墩的間距為30m, 共有 301 個橋墩; 三角形為布設在地表的 4 個臺站

    其中,()為振幅項,為場源距離,()為瑞利波相速度。計算得到 4 個臺站的垂向地震記錄如圖 5 所示。

    2.2 數(shù)值計算結(jié)果

    基于上述 4 個臺站的模擬地震記錄, 計算得到瑞利波相速度, 經(jīng)過篩選后的結(jié)果見圖 6。由于本文選取的橋墩間距為 30m, 列車行駛速度為 80m/s,對應的橋頻約為 2.67Hz, 兩倍和三倍橋頻分別為5.33Hz 和 8Hz。如圖 6(a)所示, 我們在橋頻附近的2.36~3.16Hz 范圍內(nèi)篩選出的相速度值與模型值一致。如圖 6(b)和(c)所示, 在兩倍橋頻附近, 相速度值在 5.15Hz 左右出現(xiàn)跳變, 分別在 4.28~5.08Hz 和5.19~6.11Hz 范圍內(nèi)得到正確的相速度值。圖 6(d)中, 在三倍橋頻附近的 7.39~8.51Hz 范圍內(nèi)能夠得到正確的相速度值。

    圖4 速度模型(a)和參考頻散曲線(b)

    圖5 4個臺站的 z 方向模擬地震記錄

    (a)一倍橋頻; (b)和(c)兩倍橋頻; (d)三倍橋頻?;疑€是參考頻散曲線, 黑色曲線是基于臺陣法得出的頻散曲線, 灰色虛線是橋頻及其倍頻, 黑色虛線是頻率范圍的端點值

    我們知道, 相速度本身在 5.15Hz 附近并沒有跳變, 所以這個跳變是由計算方法引起的。主要原因是計算兩倍橋頻附近的頻散曲線時相位跳躍了2π, 計算過程中忽略周期跳躍導致數(shù)值不穩(wěn)定。

    2.3 誤差分析及討論

    上述結(jié)果顯示, 多臺法能夠較有效地獲取高鐵振動產(chǎn)生波場的面波相速度值。但是, 出現(xiàn)一個明顯的問題: 相速度結(jié)果在兩倍橋頻附近出現(xiàn)跳變, 導致一小段頻率范圍內(nèi)的結(jié)果不可用。結(jié)合測試過程中臺站分布和高鐵地震波場傳播的規(guī)律, 我們認為這一誤差是由以下原因造成的。由于每一組計算涉及 3 個臺站和兩個臺站對, 若其中一對臺站的連線與某個頻率的面波主要傳播方向垂直, 從理論上講, 這兩個臺站間的相位差趨于零。在模擬數(shù)據(jù)中, 會因為數(shù)據(jù)的截斷產(chǎn)生誤差; 應用到真實數(shù)據(jù)中時, 數(shù)據(jù)本身的噪聲帶來的相位測量誤差也不容忽視。對于理論上相位相等的這兩個臺站, 產(chǎn)生的相對誤差就比較大, 從而使得相速度的測量結(jié)果不準確。解決這一的方法是利用更多的臺站數(shù)據(jù), 基于最小二乘法, 降低不確定性因素的影響。

    3 結(jié)論

    本文分析了以等間距橋墩為震源的高鐵地震信號的疊加干涉特征, 提出利用多臺站相位差提取波場視速度的方法。針對相位計算結(jié)果存在多值性的問題, 本文提出通過比較不同臺站組合計算結(jié)果的方法來獲得正確的相速度值。與雙臺法相比, 本文提出的利用多臺提取相速度的方法無需對計算值進行校正, 能夠適應更復雜的情況?;谀M的瑞利波記錄進行測試的結(jié)果表明, 所得相速度在一定頻率范圍內(nèi)與參考模型值接近, 從而證明了本文方法的有效性。

    [1]徐善輝, 郭建, 李培培, 等. 京津高鐵列車運行引起的地表振動觀測與分析. 地球物理學進展, 2017, 32(1): 421?425

    [2]蔣一然, 鮑鐵釗, 寧杰遠, 等. 高架橋下方高鐵地震信號頻譜特征研究. 北京大學學報(自然科學版), 2019, 55(5): 829?838

    [3]蔣一然, 梁萱, 寧杰遠, 等. 高鐵地震 4D地頻圖及其可用性研究. 北京大學學報(自然科學版), 2019, 55(5): 850?858

    [4]包乾宗, 許杰, 許明瑞. 高鐵地震信號時頻特征對比分析. 北京大學學報(自然科學版), 2019, 55(5): 805?812

    [5]溫景充, 寧杰遠, 張獻兵. 高鐵地震波場特點的理論分析. 北京大學學報(自然科學版), 2019, 55(5): 813–822

    [6]溫景充, 寧杰遠. 高鐵波場等效震源時間函數(shù)反演方法研究. 北京大學學報(自然科學版), 2019, 55(5): 823?828

    [7]張固瀾, 何承杰, 李勇, 等. 高鐵地震震源子波時 間函數(shù)及驗證. 地球物理學報, 2019, 62(6): 2344? 2354

    [8]胡光輝, 張固瀾, 孫思宇, 等. 高鐵地震信號直達波波形反演建模. 地球物理學報, 2019, 62(12): 4817?4824

    [9]Quiros D A, Brown L D, Kim D. Seismic interfero-metry of railroad induced ground motions: body and surface wave imaging. Geophysical Journal Interna-tional, 2016, 205(1): 301–313

    [10]Brenguier F, Boué P, Ben-Zion Y, et al. Train traffic as a powerful noise source for monitoring active faults with seismic interferometry. Geophysical Re-search Letters, 2019, 46(16): 9529?9536

    [11]王曉凱, 陳文超, 溫景充, 等. 高鐵震源地震信號的擠壓時頻分析應用. 地球物理學報, 2019, 62(6): 2328?2335

    [12]王曉凱, 陳建友, 陳文超, 等. 高鐵震源地震信號的稀疏化建模. 地球物理學報, 2019, 62(6): 2336? 2343

    [13]溫景充, 石永祥, 寧杰遠. 高鐵地震面波相速度頻散曲線提取. 地球物理學報, 2021, 64(9): 3246?3256

    [14]Forsyth D W, Webb S, Dorman L M, et al. Phase velocities of Rayleigh waves in the MELT experiment on the East Pacfic Rise. Science, 1998, 280: 1235–1238

    [15]Ritzwoller M H, Lin F C, Shen W. Ambient noise tomography with a large seismic array. Comptes Ren-dus Geosciences, 2011, 343(8): 558?570

    [16]Yao H, Van der Hilst R D, De Hoop M V. Surface-wave array tomography in SE Tibet from ambient seismic noise and two-station analysis —I. phase velocity maps. Geophysical Journal International, 2006, 166: 732?744

    [17]石永祥, 溫景充, 寧杰遠. 高鐵震源地下介質(zhì)成像的理論分析. 中國科學: 地球科學, 2022, 52(5): 893?902

    An Array-Based Method for Measuring Apparent Velocity of the High-Speed Rail Seismic Wavefield

    YIN Changyang1, WEN Jingchong1, SHI Yongxiang1, NING Jieyuan1,2,?

    1. School of Earth and Space Sciences, Peking University, Beijing 100871; 2. Hebei Hongshan National Geophysical Observation and Research Station, Xingtai 054000; ? Corresponding author, E-mail: njy@pku.edu.cn

    For the seismic source generated by the high-speed rail with viaduct, a theoretical method is proposed to measure the apparent velocity of the high-speed rail wavefield generated by a moving source set by using cross-correlation of the waves recorded by multiple stations. Furthermore, synthetic data of the high-speed rail wave-fields in layered media is generated, and the apparent velocity of a phase in it is calculated. The correctness of the method is verified compared with the theoretical results.

    high-speed rail; moving source combination; interference field; apparent velocity

    10.13209/j.0479-8023.2022.046

    2021–05–30;

    2022–05–22

    河北省地震科技星火計劃(DZ20200827054)和國家自然科學基金(41874071)資助

    猜你喜歡
    面波波場臺站
    中國科學院野外臺站檔案工作回顧
    氣象基層臺站建設
    西藏科技(2021年12期)2022-01-17 08:46:38
    gPhone重力儀的面波頻段響應實測研究
    地震研究(2021年1期)2021-04-13 01:04:56
    自適應相減和Curvelet變換組合壓制面波
    彈性波波場分離方法對比及其在逆時偏移成像中的應用
    交錯網(wǎng)格與旋轉(zhuǎn)交錯網(wǎng)格對VTI介質(zhì)波場分離的影響分析
    地震學報(2016年1期)2016-11-28 05:38:36
    基于Hilbert變換的全波場分離逆時偏移成像
    基層臺站綜合觀測業(yè)務管理之我見
    西藏科技(2015年6期)2015-09-26 12:12:13
    旋轉(zhuǎn)交錯網(wǎng)格VTI介質(zhì)波場模擬與波場分解
    MDOS平臺臺站級使用方法及技巧
    精品久久久噜噜| 一进一出抽搐动态| 亚洲婷婷狠狠爱综合网| 免费观看的影片在线观看| 国产单亲对白刺激| 内地一区二区视频在线| 色5月婷婷丁香| 亚洲精品乱码久久久久久按摩| 日韩欧美三级三区| 国产成人a∨麻豆精品| 麻豆av噜噜一区二区三区| 天堂av国产一区二区熟女人妻| 在线观看免费视频日本深夜| 亚州av有码| videossex国产| 成人亚洲精品av一区二区| 久久久久免费精品人妻一区二区| 日产精品乱码卡一卡2卡三| 国产成人福利小说| 日本欧美国产在线视频| 久久人妻av系列| 变态另类成人亚洲欧美熟女| 国产午夜精品久久久久久一区二区三区| 有码 亚洲区| 一本一本综合久久| 成年av动漫网址| 国产免费男女视频| 午夜福利在线在线| 丰满乱子伦码专区| 人妻少妇偷人精品九色| 日韩视频在线欧美| 女同久久另类99精品国产91| 丰满的人妻完整版| 亚洲av第一区精品v没综合| 欧美3d第一页| 国产精品国产高清国产av| 美女大奶头视频| 欧美一区二区亚洲| 韩国av在线不卡| 看黄色毛片网站| 成人毛片a级毛片在线播放| 国产一区亚洲一区在线观看| 久久午夜亚洲精品久久| 亚洲av成人精品一区久久| 熟女电影av网| 久久精品影院6| 国内精品宾馆在线| 亚洲人成网站高清观看| av.在线天堂| 丰满人妻一区二区三区视频av| 97热精品久久久久久| 日本在线视频免费播放| 网址你懂的国产日韩在线| 亚洲欧美日韩无卡精品| 18禁在线播放成人免费| 国产视频内射| 欧美日韩国产亚洲二区| 亚洲欧美精品专区久久| 青青草视频在线视频观看| 国产亚洲欧美98| 熟妇人妻久久中文字幕3abv| 男女做爰动态图高潮gif福利片| 69av精品久久久久久| 女的被弄到高潮叫床怎么办| av福利片在线观看| 久久精品国产清高在天天线| 听说在线观看完整版免费高清| 亚洲国产欧美在线一区| 中文欧美无线码| 国产乱人视频| 国产三级中文精品| 精品久久久久久久久亚洲| 又粗又爽又猛毛片免费看| 少妇人妻一区二区三区视频| 亚洲美女视频黄频| 亚洲一级一片aⅴ在线观看| 男女边吃奶边做爰视频| 亚洲aⅴ乱码一区二区在线播放| 国产淫片久久久久久久久| 免费看日本二区| 国产精品乱码一区二三区的特点| 美女国产视频在线观看| 午夜福利视频1000在线观看| 亚洲人成网站在线观看播放| 久久国产乱子免费精品| 久久6这里有精品| 亚洲欧美精品综合久久99| 黄色一级大片看看| 国产伦精品一区二区三区视频9| 嫩草影院入口| 成人一区二区视频在线观看| 国产精品福利在线免费观看| 一区二区三区免费毛片| 十八禁国产超污无遮挡网站| 只有这里有精品99| 少妇裸体淫交视频免费看高清| 中出人妻视频一区二区| 白带黄色成豆腐渣| www.色视频.com| 性插视频无遮挡在线免费观看| a级毛色黄片| 成人二区视频| 高清毛片免费看| 天堂影院成人在线观看| 亚洲欧美成人精品一区二区| 岛国毛片在线播放| 国产精品综合久久久久久久免费| eeuss影院久久| 国产男人的电影天堂91| 欧美zozozo另类| 国产久久久一区二区三区| 欧美区成人在线视频| 国产精品人妻久久久影院| 可以在线观看的亚洲视频| 麻豆一二三区av精品| 2021天堂中文幕一二区在线观| 亚洲av第一区精品v没综合| 久久精品国产亚洲av天美| 亚洲中文字幕日韩| 国产精品三级大全| 麻豆av噜噜一区二区三区| 国内精品一区二区在线观看| 国产男人的电影天堂91| 非洲黑人性xxxx精品又粗又长| 特级一级黄色大片| 国产精品久久久久久亚洲av鲁大| 中文欧美无线码| 草草在线视频免费看| 少妇裸体淫交视频免费看高清| 99热6这里只有精品| 亚洲人成网站高清观看| 欧美变态另类bdsm刘玥| 成人亚洲欧美一区二区av| 成年版毛片免费区| 成人三级黄色视频| 99九九线精品视频在线观看视频| 少妇人妻精品综合一区二区 | 97超碰精品成人国产| 免费观看在线日韩| 搡女人真爽免费视频火全软件| 免费看美女性在线毛片视频| 国产黄色视频一区二区在线观看 | 老师上课跳d突然被开到最大视频| 亚洲av成人精品一区久久| 99精品在免费线老司机午夜| 日韩在线高清观看一区二区三区| 人妻夜夜爽99麻豆av| 久久久成人免费电影| 日本免费a在线| 欧美丝袜亚洲另类| 久久久久久久久中文| 亚洲中文字幕一区二区三区有码在线看| 日韩一本色道免费dvd| 精品久久久久久久久亚洲| 超碰av人人做人人爽久久| 一边摸一边抽搐一进一小说| 日韩 亚洲 欧美在线| 亚洲美女搞黄在线观看| 欧美日韩国产亚洲二区| 91久久精品电影网| 91麻豆精品激情在线观看国产| 中国国产av一级| 亚洲成人久久爱视频| 级片在线观看| www.av在线官网国产| 少妇被粗大猛烈的视频| 日韩欧美一区二区三区在线观看| 一边亲一边摸免费视频| 97热精品久久久久久| 国产精品久久电影中文字幕| 欧美激情久久久久久爽电影| 日韩欧美 国产精品| 亚洲真实伦在线观看| 深爱激情五月婷婷| 久久国内精品自在自线图片| 性欧美人与动物交配| 午夜视频国产福利| 国产免费一级a男人的天堂| 国产精品久久久久久精品电影| 日韩强制内射视频| 国产真实乱freesex| 国产探花在线观看一区二区| 黄色一级大片看看| 久久99精品国语久久久| 久久久久久久久大av| 日韩欧美精品免费久久| 国产精品蜜桃在线观看 | 日本-黄色视频高清免费观看| 久久久成人免费电影| 免费搜索国产男女视频| 蜜桃亚洲精品一区二区三区| 国产精品不卡视频一区二区| 日韩人妻高清精品专区| av免费在线看不卡| 激情 狠狠 欧美| 亚洲成人久久爱视频| 精品久久久久久久久久免费视频| 国产亚洲精品久久久久久毛片| 亚洲成人中文字幕在线播放| 久久人人爽人人爽人人片va| 中国美女看黄片| 久久久久久久久久久丰满| av天堂中文字幕网| 只有这里有精品99| 偷拍熟女少妇极品色| 欧美色欧美亚洲另类二区| 日韩大尺度精品在线看网址| av在线蜜桃| 亚洲精品日韩av片在线观看| 99在线人妻在线中文字幕| 国产人妻一区二区三区在| 观看免费一级毛片| 免费观看在线日韩| 精品久久久久久久久av| 秋霞在线观看毛片| 免费一级毛片在线播放高清视频| 日本熟妇午夜| 亚洲人成网站在线播放欧美日韩| 悠悠久久av| 日本黄色视频三级网站网址| 国产av在哪里看| 91麻豆精品激情在线观看国产| av视频在线观看入口| 美女 人体艺术 gogo| 亚洲人成网站高清观看| 国产精品久久久久久久电影| 久久这里只有精品中国| 69人妻影院| 在线免费观看不下载黄p国产| 特级一级黄色大片| 男女下面进入的视频免费午夜| 69av精品久久久久久| a级毛色黄片| 尾随美女入室| 亚洲av中文av极速乱| 亚洲欧洲日产国产| 久久韩国三级中文字幕| 亚洲国产精品久久男人天堂| 老熟妇乱子伦视频在线观看| 久久99热这里只有精品18| 一本久久中文字幕| 亚洲av成人精品一区久久| 网址你懂的国产日韩在线| 久久婷婷人人爽人人干人人爱| 日韩av在线大香蕉| 亚洲第一电影网av| 身体一侧抽搐| 亚洲av免费在线观看| 国产片特级美女逼逼视频| 99久久精品一区二区三区| 亚洲av中文av极速乱| 天美传媒精品一区二区| 一卡2卡三卡四卡精品乱码亚洲| 熟妇人妻久久中文字幕3abv| 91av网一区二区| 久久韩国三级中文字幕| 久久精品国产亚洲av天美| 在线天堂最新版资源| 亚洲欧洲日产国产| 久久99热6这里只有精品| 老女人水多毛片| 看片在线看免费视频| 亚洲成人av在线免费| 国产成人91sexporn| 亚洲精品久久国产高清桃花| 国内久久婷婷六月综合欲色啪| 99久久九九国产精品国产免费| 只有这里有精品99| 免费看美女性在线毛片视频| 别揉我奶头 嗯啊视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产不卡一卡二| 亚洲av免费高清在线观看| 看黄色毛片网站| 1000部很黄的大片| 成人漫画全彩无遮挡| 亚洲av不卡在线观看| 女同久久另类99精品国产91| 免费观看a级毛片全部| 亚洲欧美日韩卡通动漫| 一级av片app| 国产精品国产高清国产av| 看非洲黑人一级黄片| 国产黄色小视频在线观看| 性欧美人与动物交配| 一个人看的www免费观看视频| 亚洲美女视频黄频| 欧美xxxx黑人xx丫x性爽| 啦啦啦韩国在线观看视频| 久久99蜜桃精品久久| 一区二区三区四区激情视频 | 国产三级在线视频| 69av精品久久久久久| 国产综合懂色| 不卡一级毛片| 我要看日韩黄色一级片| av在线老鸭窝| 97人妻精品一区二区三区麻豆| 麻豆国产av国片精品| 黄色一级大片看看| 久久久久久大精品| 久久久a久久爽久久v久久| 日本色播在线视频| 精品欧美国产一区二区三| 欧美最黄视频在线播放免费| 国产精品综合久久久久久久免费| 成年女人永久免费观看视频| 国产视频内射| 精品不卡国产一区二区三区| 久久久精品94久久精品| 热99re8久久精品国产| 搞女人的毛片| 色综合亚洲欧美另类图片| 成人av在线播放网站| 久久亚洲国产成人精品v| 亚洲欧洲国产日韩| 国产精品嫩草影院av在线观看| 日韩成人伦理影院| 久久中文看片网| 欧美成人a在线观看| 精品99又大又爽又粗少妇毛片| 黄色配什么色好看| 成人国产麻豆网| 欧美zozozo另类| 老司机福利观看| 成年女人看的毛片在线观看| 久久久久久国产a免费观看| 日本撒尿小便嘘嘘汇集6| 婷婷亚洲欧美| 男女做爰动态图高潮gif福利片| 久久国内精品自在自线图片| 中文字幕av在线有码专区| 免费电影在线观看免费观看| 91精品国产九色| 一级黄色大片毛片| 亚洲av成人av| 国产视频内射| 成人永久免费在线观看视频| 只有这里有精品99| 卡戴珊不雅视频在线播放| 一卡2卡三卡四卡精品乱码亚洲| 看十八女毛片水多多多| 国产黄a三级三级三级人| 一级av片app| 高清午夜精品一区二区三区 | 99久久人妻综合| 天堂中文最新版在线下载 | 26uuu在线亚洲综合色| 免费观看的影片在线观看| 国产精品野战在线观看| 精华霜和精华液先用哪个| 老司机影院成人| 欧美日韩国产亚洲二区| 在线观看一区二区三区| 亚洲精品久久国产高清桃花| 国产精品伦人一区二区| 人妻夜夜爽99麻豆av| 97热精品久久久久久| 日韩成人伦理影院| 97超视频在线观看视频| 可以在线观看的亚洲视频| 不卡视频在线观看欧美| 婷婷色综合大香蕉| 午夜免费激情av| 中文字幕久久专区| 欧美+亚洲+日韩+国产| 欧美激情在线99| 天堂网av新在线| 日产精品乱码卡一卡2卡三| 夜夜爽天天搞| 亚洲国产精品国产精品| 身体一侧抽搐| 亚洲最大成人av| 国产精品久久久久久久久免| 国产精品1区2区在线观看.| 天堂√8在线中文| 午夜精品国产一区二区电影 | 亚洲精品日韩在线中文字幕 | av福利片在线观看| a级毛片免费高清观看在线播放| 国产久久久一区二区三区| 日产精品乱码卡一卡2卡三| 午夜老司机福利剧场| 亚洲成av人片在线播放无| 国产精品久久久久久精品电影| 成人性生交大片免费视频hd| 亚洲丝袜综合中文字幕| 欧美一区二区精品小视频在线| 91久久精品国产一区二区三区| 简卡轻食公司| 久久久欧美国产精品| 国产精品免费一区二区三区在线| 久久人人爽人人片av| 亚洲性久久影院| 亚洲国产欧洲综合997久久,| 能在线免费观看的黄片| h日本视频在线播放| 精品久久久久久久久av| 免费黄网站久久成人精品| 国产精品久久久久久久电影| 国产一区二区三区在线臀色熟女| 日本三级黄在线观看| 日本-黄色视频高清免费观看| 青青草视频在线视频观看| 日本三级黄在线观看| 久久精品人妻少妇| 色播亚洲综合网| 日韩一区二区视频免费看| 九九久久精品国产亚洲av麻豆| 狂野欧美白嫩少妇大欣赏| 国产精品乱码一区二三区的特点| 亚洲国产欧美在线一区| 国产精品蜜桃在线观看 | 长腿黑丝高跟| 高清毛片免费看| 大型黄色视频在线免费观看| 99国产精品一区二区蜜桃av| 精品久久久久久久人妻蜜臀av| 菩萨蛮人人尽说江南好唐韦庄 | 性色avwww在线观看| 亚洲性久久影院| 激情 狠狠 欧美| 国产乱人偷精品视频| 亚洲自拍偷在线| 国产片特级美女逼逼视频| 在线天堂最新版资源| 美女xxoo啪啪120秒动态图| 99热这里只有精品一区| 午夜福利高清视频| 97热精品久久久久久| 亚洲欧美日韩东京热| 91久久精品国产一区二区三区| 久久久久久久亚洲中文字幕| 搡女人真爽免费视频火全软件| 熟女人妻精品中文字幕| 国产极品精品免费视频能看的| 日韩成人av中文字幕在线观看| 中文字幕av在线有码专区| 国产一区二区在线av高清观看| 免费观看精品视频网站| 久久久精品欧美日韩精品| 日韩国内少妇激情av| 18禁在线无遮挡免费观看视频| av.在线天堂| 日韩 亚洲 欧美在线| 亚洲乱码一区二区免费版| 亚洲人成网站在线播| 91麻豆精品激情在线观看国产| 青春草视频在线免费观看| 久久久久久大精品| 国产毛片a区久久久久| 性色avwww在线观看| 91久久精品电影网| 一级av片app| 精品一区二区免费观看| 日韩视频在线欧美| 激情 狠狠 欧美| 欧美激情久久久久久爽电影| 亚洲四区av| 神马国产精品三级电影在线观看| 青春草国产在线视频 | 校园人妻丝袜中文字幕| 插逼视频在线观看| 国产精品久久久久久精品电影小说 | 免费一级毛片在线播放高清视频| 嫩草影院精品99| 亚洲国产精品合色在线| 午夜福利视频1000在线观看| 久久欧美精品欧美久久欧美| 国产91av在线免费观看| 久久久久久久久久黄片| 久久久久久久久中文| 日日干狠狠操夜夜爽| 在线天堂最新版资源| 欧美xxxx性猛交bbbb| 免费观看a级毛片全部| 日韩欧美国产在线观看| a级毛片免费高清观看在线播放| 日本熟妇午夜| 久久国内精品自在自线图片| 亚洲真实伦在线观看| 亚洲久久久久久中文字幕| av又黄又爽大尺度在线免费看 | 国产精品一区二区三区四区免费观看| 18禁裸乳无遮挡免费网站照片| 精品人妻偷拍中文字幕| 成人亚洲精品av一区二区| 亚州av有码| 国产亚洲精品av在线| av卡一久久| 边亲边吃奶的免费视频| 极品教师在线视频| 亚洲自偷自拍三级| 麻豆成人av视频| 在线观看午夜福利视频| ponron亚洲| 欧美成人精品欧美一级黄| 久久中文看片网| 国产精品人妻久久久久久| 夜夜爽天天搞| 天天躁夜夜躁狠狠久久av| 晚上一个人看的免费电影| 亚州av有码| 久久久久免费精品人妻一区二区| 一级毛片我不卡| 成人毛片60女人毛片免费| 一级av片app| www.色视频.com| 国产 一区 欧美 日韩| 99久久中文字幕三级久久日本| 欧美另类亚洲清纯唯美| 自拍偷自拍亚洲精品老妇| 久久精品夜色国产| 99国产极品粉嫩在线观看| 午夜激情欧美在线| 国产免费一级a男人的天堂| 给我免费播放毛片高清在线观看| 夫妻性生交免费视频一级片| 国产成人a∨麻豆精品| 三级经典国产精品| 嫩草影院入口| 精品不卡国产一区二区三区| 日韩欧美精品免费久久| 高清午夜精品一区二区三区 | 天堂中文最新版在线下载 | 久久热精品热| 蜜桃久久精品国产亚洲av| 床上黄色一级片| 岛国毛片在线播放| 成人av在线播放网站| 3wmmmm亚洲av在线观看| 免费人成在线观看视频色| 嫩草影院新地址| 三级国产精品欧美在线观看| 如何舔出高潮| 亚洲内射少妇av| 国产精品一区二区三区四区免费观看| 一区二区三区四区激情视频 | 好男人视频免费观看在线| 中文亚洲av片在线观看爽| 青春草国产在线视频 | 精品久久久久久成人av| 97在线视频观看| 人人妻人人澡人人爽人人夜夜 | 91久久精品国产一区二区三区| 三级经典国产精品| 美女 人体艺术 gogo| 国产精品久久久久久精品电影| 亚洲精华国产精华液的使用体验 | 日韩av不卡免费在线播放| 亚洲国产精品久久男人天堂| 狂野欧美激情性xxxx在线观看| 国产老妇女一区| 久久99热6这里只有精品| 在线免费观看的www视频| 国产免费男女视频| 麻豆国产av国片精品| 在线观看一区二区三区| 小蜜桃在线观看免费完整版高清| 成人无遮挡网站| av在线观看视频网站免费| 午夜福利在线在线| 亚洲第一电影网av| 欧美性猛交╳xxx乱大交人| 夜夜看夜夜爽夜夜摸| 成人鲁丝片一二三区免费| 亚洲熟妇中文字幕五十中出| 国产精品国产高清国产av| 日韩 亚洲 欧美在线| 日韩欧美一区二区三区在线观看| 尾随美女入室| 国产精品嫩草影院av在线观看| 看非洲黑人一级黄片| 午夜视频国产福利| 国产成人影院久久av| 深夜a级毛片| 天堂影院成人在线观看| 日本免费a在线| 国产高潮美女av| 日本三级黄在线观看| 黄色视频,在线免费观看| 午夜福利在线观看免费完整高清在 | 亚洲欧美精品自产自拍| 成熟少妇高潮喷水视频| 综合色av麻豆| 天堂中文最新版在线下载 | 男的添女的下面高潮视频| 国产精品永久免费网站| 日韩强制内射视频| 国产精品1区2区在线观看.| 99热这里只有是精品在线观看| 亚洲国产欧美在线一区| 伦理电影大哥的女人| 少妇熟女欧美另类| 亚洲精品色激情综合| 色5月婷婷丁香| 国产片特级美女逼逼视频| 九九久久精品国产亚洲av麻豆| 91狼人影院| 国模一区二区三区四区视频| av在线蜜桃| 国产69精品久久久久777片| 婷婷色综合大香蕉| 欧美极品一区二区三区四区| 婷婷六月久久综合丁香| 国产久久久一区二区三区| 久久99热这里只有精品18| 色尼玛亚洲综合影院| 国产蜜桃级精品一区二区三区| 精品久久久久久久久久免费视频| 观看美女的网站| 国产乱人视频| 又黄又爽又刺激的免费视频.| 日日摸夜夜添夜夜爱| 亚洲真实伦在线观看|