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

    TSP隧道地震預(yù)報空洞模型的數(shù)值模擬研究

    2016-11-17 02:32:50王朝令
    大地測量與地球動力學(xué) 2016年11期
    關(guān)鍵詞:模型

    王朝令 楊 茜 陳 章 楊 文

    1 四川農(nóng)業(yè)大學(xué)土木工程學(xué)院,都江堰市建設(shè)路288號,611830

    ?

    TSP隧道地震預(yù)報空洞模型的數(shù)值模擬研究

    王朝令1楊 茜1陳 章1楊 文1

    1 四川農(nóng)業(yè)大學(xué)土木工程學(xué)院,都江堰市建設(shè)路288號,611830

    采用有限元方法模擬隧道地震波場,采用波場快照與時間記錄相結(jié)合的方法,研究空洞對隧道地震波場傳播的影響,并對含空洞模型的時間記錄進(jìn)行處理,得到數(shù)值模型的速度云圖和反射層位圖。經(jīng)數(shù)據(jù)處理表明,采用TSPwin設(shè)定默認(rèn)值處理得到的速度云圖與模型設(shè)定的空洞位置比較一致,在提取反射層位圖上,空洞反射層呈現(xiàn)條帶狀特征,需要結(jié)合速度云圖來確定空洞位置,且P波預(yù)報的準(zhǔn)確性相對較高。對TSP系統(tǒng)的抗噪性進(jìn)行研究,表明其具有良好的抗噪性能。最后通過對工程實例的處理,驗證了數(shù)值模擬所得的結(jié)論。

    TSP;TSPwin;空洞;波場快照;隧道地震預(yù)報

    溶洞內(nèi)填充物與周圍巖體具有明顯的波阻抗差異,當(dāng)?shù)卣鸩▓鰝鞑サ讲ㄗ杩菇缑鏁r,會產(chǎn)生明顯的反射,因此其具有良好的地球物理特征。鐘世航[1]指出,TSP預(yù)報失誤的主要原因在于對探測對象的地質(zhì)特征及地球物理特征認(rèn)識不清。李術(shù)才等[2]對巖溶裂隙水和不良地質(zhì)情況在超前預(yù)報中存在的主要問題進(jìn)行了總結(jié)。孫克國等[3]采用TSP預(yù)報方法對多個巖溶隧道進(jìn)行預(yù)報,總結(jié)出洞穴和暗河的判別規(guī)律。許振浩等[4]研究巖溶隧道內(nèi)斷層破碎帶、地下水充填和無充填型溶洞的地震波反射特性。王朝令[5]基于有限元方法,應(yīng)用數(shù)值模擬研究隧道地震波場的數(shù)值模擬、分離和反演。林義等[6]應(yīng)用數(shù)值模擬方法研究了圍巖松動圈對隧道地震預(yù)報系統(tǒng)的影響。

    本文基于有限元方法模擬空洞存在情形下隧道地震波場的傳播過程,并用TSPwin軟件對模擬結(jié)果進(jìn)行處理分析,得到速度分布云圖和反射層提取圖,將兩者與模型圖進(jìn)行比較,分析TSP預(yù)報的準(zhǔn)確性。然后在數(shù)據(jù)中加入噪聲進(jìn)行分析,研究TSPwin的抗噪性能。最后以實例驗證結(jié)果的準(zhǔn)確性。

    1 空洞模型的數(shù)值模擬

    為模擬空洞模型的隧道地震波場傳播情形,建立模型時需要結(jié)合計算能力,并能滿足分析需要。若模型太大,計算時間長,不利于進(jìn)行修改模擬。為此,建立如圖1所示的模型,模型長135 m,高75 m,在隧道模型的前方設(shè)置一個圓形空洞,其圓心坐標(biāo)為(87,5),半徑為20 m,內(nèi)部填充P波為0,S波為0。檢波器設(shè)置與上述兩個模型相同,都位于隧道的邊壁上,震源激發(fā)采用主頻為300 Hz的零相位Ricker子波。當(dāng)主頻為300 Hz時,網(wǎng)格劃分長度為0.5 m,可以最大程度地減小頻散效應(yīng)[1]。采樣間隔為62.5 μs,點數(shù)為7 128,與TSP的設(shè)置相同。

    圖1 模型圖Fig.1 Model diagram

    TSP數(shù)據(jù)采集時,采用單個三分量接收,24炮激發(fā),根據(jù)炮檢互換原理,相當(dāng)于單炮激發(fā),24道三分量檢波器接收。為方便計算,在數(shù)值模擬時,采用單炮激發(fā)多道接收,使得計算模型的觀測系統(tǒng)與TSP實際采集時的相同。需要說明的是,在模型數(shù)據(jù)處理過程中,速度云圖和提取反射層圖都是從震源位置起算。此模型中,空洞左邊緣距離震源67 m,右邊緣距離震源107 m。

    由于模擬是在二維中進(jìn)行的,只能取到x、y兩個分量。圖2是提取的時間記錄,圖中的直達(dá)縱波、直達(dá)橫波都比較清晰,第一個同相軸的反射橫波是RSSR現(xiàn)象引起的,第二個反射同相軸是前方界面的反射橫波。隧道內(nèi)沿邊壁傳播的Rayleigh面波在工作面角點位置轉(zhuǎn)換為S波,以體波的形式向前傳播,遇到界面反射回到角點后,轉(zhuǎn)換為Rayleigh面波,此即為RSSR現(xiàn)象[2]。當(dāng)?shù)卣鸩ㄓ龅娇斩磿r,x分量相對于y分量來說要復(fù)雜一些。由于半徑為20m的空洞,其直徑大于第一Fresnel帶,因此整個時間剖面存在反射縱橫波,其中直達(dá)橫波和反射橫波能量強,所以比較明顯。

    圖2 數(shù)值模擬的時間剖面Fig.2 Time recording of numerical simulation

    圖3為空洞模型提取的波場快照,T=7 ms顯示震源所激發(fā)的P波、S波開始分離;T=10 ms、15 ms時,波場傳播到接收排列上;T=20 ms、25 ms、30 ms過程中,波場傳播過接收排列,并到達(dá)工作面位置,在角點上產(chǎn)生RSSR現(xiàn)象;T=33 ms、35 ms,經(jīng)RSSR轉(zhuǎn)換出來的橫波與原來的波場產(chǎn)生疊加一起向外傳播,產(chǎn)生新的波場;T=37 ms、40 ms過程中,橫波傳播到空洞位置,此時P波已經(jīng)到達(dá)下邊界;T=42 ms、45 ms,由于空洞的存在,產(chǎn)生繞射現(xiàn)象;T=49 ms、52 ms、55 ms、59 ms時反射波由空洞反射到工作面,可以觀察到明顯的繞射現(xiàn)象;T=63 ms、68 ms,經(jīng)過RSSR轉(zhuǎn)化的Rayleigh面波繼續(xù)沿排列方向傳播;T=72 ms、78 ms時,可以觀察到左邊界有波反射回來,這是由繞射現(xiàn)象產(chǎn)生的波場,同時反射波繼續(xù)沿排列傳播,此現(xiàn)象可以與圖2的時間剖面進(jìn)行對比,形成對應(yīng)關(guān)系,此時波場已經(jīng)變得非常復(fù)雜。

    圖3 模擬y分量的波場快照Fig.3 Y-component snapshot of simulation

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

    按照TSP的處理流程,應(yīng)用TSPwin對得到的時間記錄進(jìn)行處理,各處理步驟均采用TSPwin默認(rèn)值進(jìn)行,得到如圖4所示的速度云圖。圖4中,橫坐標(biāo)表示預(yù)報距離(以震源位置起算),縱坐標(biāo)表示以隧道向兩邊展布半徑,兩者的單位均為m。在計算地震波傳播過程中,如果所建立的速率模型與巖層中的傳播速度一樣,那么任何反射界面的位置都與其傳播距離無關(guān)。與模型進(jìn)行對比可知,低速異常區(qū)的P波速在2 000 m/s左右,S波在1 000 m/s左右,且所在位置與模型設(shè)定的位置非常接近,表明對空洞模型而言,速度分析云圖可以提供比較準(zhǔn)確的定位。

    圖4 空洞模型的速度圖Fig.4 Velocity contours of cave model

    圖5是20 m半徑空洞模型經(jīng)過偏移后處理得到的反射層位分布??梢园l(fā)現(xiàn),圖中反射層位分布比較明顯,模型中的分布位置相符,說明TSP在探測空洞時,反射層所呈現(xiàn)的形態(tài)是以前后界面來表示的,這增加了空洞在提取反射層圖上的設(shè)定難度。此時,結(jié)合速度分布云圖,可以提高準(zhǔn)確性。相比于P波,S波的反射層能量較弱,在提取反射層圖上不太明顯,因此P波的預(yù)報效果更好。

    圖5 提取反射層Fig.5 Layers of extracting reflection

    3 TSP抗噪性分析

    為研究TSP的抗干擾能力,對數(shù)值模擬所提取的時間記錄加入隨機噪聲,得到圖6所示的低信噪比信號。從圖6可以看到,信號湮沒在隨機噪聲中,整個剖面不容易辨識,x分量中除了直達(dá)波比較明顯外,反射波幾乎被噪聲所掩蓋,y分量信噪比稍高。

    圖6 加噪聲的時間記錄Fig.6 Time recording with additive noise

    對加噪聲信號采用TSPwin進(jìn)行處理,經(jīng)過帶通濾波、初至拾取、反射波提取、P-S波分離、深度偏移和反射層提取,得到圖7所示的速度云圖和圖8的反射層提取圖。比較圖7的加噪聲速度云圖與未加噪聲的速度云圖,P波的速度為2 000 m/s,S波的速度為1 000 m/s左右;P波的速度云圖與設(shè)定的位置非常相近,而S波速度云圖相差較大,故P波的抗干擾能力更強,精度更高。

    圖7 加噪聲的速度云圖Fig.7 Velocity contours with additive noise

    圖8 加噪聲的反射層提取Fig.8 Layers of extracting reflection with additive noise

    圖8為加噪聲數(shù)據(jù)所提取的反射層。對比模型圖1可以看到,在70~100 m區(qū)間內(nèi), P波反射界面所處位置與模型設(shè)置相近,且形態(tài)也較一致,但是相比于圖8未加噪聲的P波反射層提取,對應(yīng)性有所降低,且反射層強度也有所下降。圖8(b)為提取的S波反射層,可以看到,在70~100 m之間存在反射層,接近模型的設(shè)定位置,與圖7未加噪聲的S波提取反射層相比,相差較大。因此,就預(yù)報精度而言,P波更好。

    通過對加噪聲的數(shù)據(jù)處理表明,TSP具有良好的抗噪性能。TSP在數(shù)據(jù)采集時會采用在隧道邊壁打孔,用環(huán)氧樹脂將套管固定在接收孔中,然后將傳感器與套管緊密耦合,以此接收地震信號,使得原始時間記錄具有很高的信噪比,因此TSP系統(tǒng)所具有的抗噪性能可以滿足實際工程需要。

    4 工程實例

    考慮到隧道超前預(yù)報方法的施工特點,若要驗證預(yù)報結(jié)果的準(zhǔn)確性,則只能選取已經(jīng)開挖完的隧道,比較開挖之前的預(yù)報結(jié)果和開挖后揭示的地質(zhì)情況,對比兩者的差異。依據(jù)常規(guī)TSP觀測排列的布置方法,在隧道邊墻布置地震波信息接收孔,孔深2 m,按約1.5 m的間距布置24個激發(fā)孔,孔深1.5 m左右,向下傾斜約10°,每個激發(fā)孔裝填藥量為100 g,激發(fā)孔和接收孔基本保持在同一個高度。采用上文所述處理流程中設(shè)定的參數(shù)對此數(shù)據(jù)進(jìn)行處理,得到圖9所示的速度云圖和反射層提取圖。

    圖9 P波速度云圖和反射層圖Fig.9 Velocity contours and extracted reflection of P wave

    圖9(b)反射層圖顯示,在110~140 m段反射層密集,且正負(fù)反射交錯出現(xiàn),推測發(fā)育有溶洞或是破碎帶。根據(jù)圖9(a)速度云圖,由于低速區(qū)呈區(qū)塊狀分布,推測為空洞;隧道開挖結(jié)果顯示,在110 m處溶洞覆蓋掌子面左側(cè)拱腰處有溶洞,溶洞規(guī)模長33.7 m、寬32.4 m、高24 m,溶洞深度超過隧底高程,整個溶洞靠線路右側(cè),并超過右邊墻8 m以上。這說明預(yù)報結(jié)果與開挖結(jié)果非常一致,采用文中的參數(shù)對數(shù)據(jù)進(jìn)行處理,可以比較準(zhǔn)確地預(yù)報前方的地質(zhì)災(zāi)害。隧道地震超前預(yù)報作為一種地球物理勘探方法,其結(jié)果存在多解性,在參照速度云圖和反射層圖時,需要結(jié)合兩者作出判斷,并在施工時進(jìn)行多種準(zhǔn)備,應(yīng)對可能發(fā)生的施工事故。

    5 結(jié) 語

    TSP在工程中應(yīng)用非常廣泛,但目前國內(nèi)外都將其作為一種工程方法進(jìn)行實例預(yù)報研究,采用正演模擬并進(jìn)行處理分析的研究相對較少。本文建立溶洞模型進(jìn)行數(shù)值模擬,通過波場快照研究隧道地震波場的傳播特征,對得到的時間剖面按照TSP的常規(guī)流程進(jìn)行處理,分析空洞模型速度云圖和提取反射層圖的特征。在速度分布云圖中,空洞的速度呈現(xiàn)比周圍速度低的特征,在提取的反射層圖中空洞呈現(xiàn)豎向條狀包圍的特征,因此需要結(jié)合速度分析空洞的位置和形態(tài)。分析表明,P波的準(zhǔn)確度較高。對數(shù)值模擬所得的時間記錄添加噪聲,研究表明,TSP具有良好的抗噪性能。最后通過工程實例表明,P波預(yù)報是比較準(zhǔn)確的。

    本文以空洞模型為例,采用有限元數(shù)值計算模擬了隧道的地震波場,經(jīng)過對時間剖面的反演處理,分析空洞在速度分布云圖和反射層圖中呈現(xiàn)的特征,并研究了TSP系統(tǒng)的抗噪性能,為工程中應(yīng)用TSP進(jìn)行隧道地震預(yù)報的數(shù)據(jù)處理和資料解釋提供理論依據(jù)。

    [1] 鐘世航.TSP作隧道掌子面前方地質(zhì)預(yù)報幾例失誤原因分析[J]. 巖石力學(xué)與工程學(xué)報, 2003,22(增1): 2 443-2 446(Zhong Shihang.Analysis of the Failure Reasons in Survey Forward from Working Face of Tunnels by Means of TSP[J].Chinese Journal of Rock Mechanics and Engineering,2003,22(S1):2 443-2 446)

    [2] 李術(shù)才, 李樹忱, 張慶松,等. 巖溶裂隙水與不良地質(zhì)情況超前預(yù)報研究[J]. 巖石力學(xué)與工程學(xué)報, 2007, 26(2):217-225(Li Shucai, Li Shuchen, Zhang Qingsong, et al.Forecast of Karst-Fractured Groundwater and Defective Geological Conditions[J]. Chinese Journal of Rock Mechenics and Engineering, 2007,26(2):217-225)

    [3] 孫克國, 李術(shù)才, 張慶松,等.TSP在巖溶區(qū)山嶺隧道預(yù)報中的應(yīng)用研究[J]. 山東大學(xué)學(xué)報:工學(xué)版,2008,38(1): 74-79(Sun Keguo, Li Shucai, Zhang Qingsong, et al. Application of the TSP Geological Forecast Method of a Mountain Tunnel in a Karst Zone[J]. Journal of Shandong University:Engineering Science,2008,38(1):74-79)

    [4] 許振浩, 李術(shù)才, 張慶松,等.TSP超前地質(zhì)預(yù)報地震波反射特性研究[J]. 地下空間與工程學(xué)報,2008,4(4): 640-644(Xu Zhenhao, Li Shucai, Zhang Qingsong, et al. Reflection Characteristic of Seismic Wave in TSP Advance Geological Predicion[J].Chinese Journal of Underground Space and Engineering,2008,4(4): 640-644)

    [5] 王朝令.隧道地震超前預(yù)報中波場分離與反演方法的數(shù)值模擬研究[D]. 成都: 西南交通大學(xué),2012(Wang Zhaoling. The Numerical Simulation Research on Separation of Wave Field and Inversion Method in Tunnel Seismic Prediction[D]. Chengdu: Southwest Jiaotong University,2012)

    [6] 林義,劉爭平,王朝令.圍巖松動圈對TSP預(yù)報的影響研究[J]. 地球物理學(xué)進(jìn)展,2015,30(2): 912-919(Lin Yi, Liu Zhengping, Wang Zhaoling, et.al Effects of Loosing Circle in Surrounding Rock of TSP[J].Progress in Geophysics, 2015, 30(2):912-919)

    About the first author:WANG Zhaoling, PhD, associate professor, majors in tunnel seismic prediction, E-mail:wong8010@gmail.com.

    Numerical Simulation of TSP Tunnel Cave Model of Seismic Prediction

    WANGZhaoling1YANGQian1CHENZhang1YANGWen1

    1 College of Civil Engineering, Sichuan Agriculture University, 288 Jianshe Road, Dujiangyan 611830, China

    Tunnel excavation encounters a variety of geological challenges, such as faults, caves, etc. Currently tunnel seismic prediction adopts TSP, which is used widely, but research about TSP is currently focused on its engineering application examples. In this paper, we use the finite element method to simulate the tunnel seismic wave field, employing wave field snapshots and combined time recording methods to study the impact of caves on the characteristics of the propagation of seismic wave field tunnel. We also employ inversion processing of the time record of the model containing the cave. We compare our method with the velocity scattered image and reflection interface position, velocity processing with the default values set by using TSPwin. The results are good for tunnel seismic velocity cave layered model with single prediction, and have a higher prediction accuracy of the P-wave. Finally, we process and analyze engineering examples to verify the numerical simulation the conclusions.

    TSP;TSPwin; cave;snapshot;tunnel seismic prediction

    National Natural Science Foundation of China, No.5150835.

    2015-12-14

    項目來源:國家自然科學(xué)基金(51508358)。

    王朝令,博士,副教授,主要研究方向為隧道地震超前預(yù)報,E-mail:wong8010@gmail.com。

    10.14075/j.jgg.2016.11.012

    1671-5942(2016)011-0998-05

    P315

    A

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機模型
    提煉模型 突破難點
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    3D打印中的模型分割與打包
    精品久久久久久久久亚洲| 久久久国产成人免费| 亚洲欧美日韩无卡精品| 色5月婷婷丁香| 在线免费观看不下载黄p国产| 成年女人永久免费观看视频| 国产av在哪里看| 色5月婷婷丁香| 婷婷色av中文字幕| 天美传媒精品一区二区| 精品国产三级普通话版| 国产成人freesex在线| 亚洲欧美精品综合久久99| 一区二区三区四区激情视频| www日本黄色视频网| 在线观看66精品国产| 夜夜爽夜夜爽视频| 国国产精品蜜臀av免费| 亚洲欧美成人精品一区二区| 国产69精品久久久久777片| 国产黄色小视频在线观看| 精品久久久久久久久亚洲| 午夜福利高清视频| 国产极品精品免费视频能看的| 高清视频免费观看一区二区 | 国产精品无大码| 欧美97在线视频| 亚洲欧美中文字幕日韩二区| 亚洲精品国产成人久久av| 亚洲精品自拍成人| 亚洲成av人片在线播放无| 午夜激情福利司机影院| 淫秽高清视频在线观看| 午夜福利在线观看免费完整高清在| 国产亚洲av片在线观看秒播厂 | 久久99蜜桃精品久久| 欧美日本视频| 成年免费大片在线观看| 18禁在线播放成人免费| 亚洲成人精品中文字幕电影| 高清视频免费观看一区二区 | 免费黄色在线免费观看| 国产高清国产精品国产三级 | 亚洲欧美成人精品一区二区| 国产白丝娇喘喷水9色精品| 一本一本综合久久| 亚洲av成人av| 久久精品人妻少妇| 天天一区二区日本电影三级| 中文字幕免费在线视频6| 亚洲高清免费不卡视频| 日本wwww免费看| 国产av码专区亚洲av| 欧美97在线视频| 99热6这里只有精品| 超碰av人人做人人爽久久| 中文欧美无线码| 精品少妇黑人巨大在线播放 | 亚洲欧美精品自产自拍| 免费人成在线观看视频色| 中文字幕精品亚洲无线码一区| 国产一区二区在线av高清观看| 亚洲精品,欧美精品| www.av在线官网国产| 国产精品久久电影中文字幕| 日韩亚洲欧美综合| 丰满少妇做爰视频| 国产精品,欧美在线| 亚洲精品亚洲一区二区| 免费看光身美女| 乱码一卡2卡4卡精品| 六月丁香七月| 亚洲av成人精品一二三区| 欧美精品一区二区大全| 亚洲最大成人中文| 中文字幕制服av| 久99久视频精品免费| 99久久无色码亚洲精品果冻| 亚洲精品国产成人久久av| 七月丁香在线播放| 国产亚洲5aaaaa淫片| 全区人妻精品视频| 免费av观看视频| 午夜福利成人在线免费观看| 国产麻豆成人av免费视频| 乱人视频在线观看| 亚洲综合色惰| 亚洲真实伦在线观看| 色网站视频免费| 国产午夜精品论理片| 国产高清三级在线| 欧美成人精品欧美一级黄| 成人性生交大片免费视频hd| 亚洲高清免费不卡视频| 美女xxoo啪啪120秒动态图| 亚洲一级一片aⅴ在线观看| 2021天堂中文幕一二区在线观| 一级黄色大片毛片| 精品久久久久久久人妻蜜臀av| 一卡2卡三卡四卡精品乱码亚洲| 成人性生交大片免费视频hd| 成人无遮挡网站| 国产视频首页在线观看| 日韩欧美三级三区| 黄片无遮挡物在线观看| 午夜福利高清视频| 日本免费在线观看一区| 亚洲性久久影院| 国产一区有黄有色的免费视频 | 欧美性猛交╳xxx乱大交人| 国产亚洲91精品色在线| 亚洲国产精品专区欧美| 久久99蜜桃精品久久| 久久久亚洲精品成人影院| 免费观看在线日韩| 黄色日韩在线| 九草在线视频观看| 日韩在线高清观看一区二区三区| 联通29元200g的流量卡| 又爽又黄a免费视频| 亚洲av电影不卡..在线观看| 精品无人区乱码1区二区| 日本欧美国产在线视频| 国模一区二区三区四区视频| 久久99热这里只频精品6学生 | 欧美色视频一区免费| 好男人视频免费观看在线| 免费大片18禁| 久久久精品94久久精品| 国产高清视频在线观看网站| 天堂√8在线中文| 欧美日本亚洲视频在线播放| 1000部很黄的大片| 插逼视频在线观看| 成年版毛片免费区| 日韩,欧美,国产一区二区三区 | 丰满人妻一区二区三区视频av| 久久国产乱子免费精品| 两个人的视频大全免费| 国产老妇女一区| 亚洲av二区三区四区| 国产伦一二天堂av在线观看| 精品久久久久久久人妻蜜臀av| 大又大粗又爽又黄少妇毛片口| 国产高清三级在线| videos熟女内射| 级片在线观看| 国产v大片淫在线免费观看| 美女高潮的动态| 韩国av在线不卡| 日韩 亚洲 欧美在线| 亚洲欧美日韩高清专用| 免费观看性生交大片5| 天美传媒精品一区二区| 狂野欧美激情性xxxx在线观看| 美女大奶头视频| 日韩av在线免费看完整版不卡| 国产一级毛片七仙女欲春2| 夜夜看夜夜爽夜夜摸| 欧美最新免费一区二区三区| 亚洲经典国产精华液单| 国内精品一区二区在线观看| 欧美激情国产日韩精品一区| 欧美+日韩+精品| 麻豆成人av视频| 国产伦一二天堂av在线观看| 国产伦理片在线播放av一区| 小蜜桃在线观看免费完整版高清| 亚洲丝袜综合中文字幕| 欧美97在线视频| 亚洲丝袜综合中文字幕| 日韩一区二区视频免费看| 久久精品熟女亚洲av麻豆精品 | 男人的好看免费观看在线视频| 99久久无色码亚洲精品果冻| 一个人观看的视频www高清免费观看| 久久6这里有精品| 日韩av在线大香蕉| 性插视频无遮挡在线免费观看| 美女内射精品一级片tv| 毛片一级片免费看久久久久| 91久久精品国产一区二区成人| 毛片一级片免费看久久久久| 18禁在线无遮挡免费观看视频| 毛片一级片免费看久久久久| 一区二区三区免费毛片| 日韩av在线大香蕉| 国模一区二区三区四区视频| 久久久久九九精品影院| 男人和女人高潮做爰伦理| 国产精品嫩草影院av在线观看| 国产片特级美女逼逼视频| 精品一区二区三区视频在线| 精品一区二区三区视频在线| 亚洲成人久久爱视频| 最近2019中文字幕mv第一页| 热99re8久久精品国产| 久久精品久久久久久久性| 亚洲精品456在线播放app| 午夜精品国产一区二区电影 | 国产白丝娇喘喷水9色精品| 三级男女做爰猛烈吃奶摸视频| 美女脱内裤让男人舔精品视频| 国产极品精品免费视频能看的| 一区二区三区高清视频在线| 成人漫画全彩无遮挡| 亚洲国产欧洲综合997久久,| 一级毛片电影观看 | 精品欧美国产一区二区三| 熟女人妻精品中文字幕| 秋霞在线观看毛片| 99久国产av精品| 久久精品国产亚洲av天美| 99久久成人亚洲精品观看| 国产成人精品一,二区| 中文字幕人妻熟人妻熟丝袜美| 国产亚洲av片在线观看秒播厂 | 国产精华一区二区三区| 久久国内精品自在自线图片| 水蜜桃什么品种好| 国产视频内射| 如何舔出高潮| 级片在线观看| 看片在线看免费视频| 午夜激情欧美在线| 国产在视频线在精品| 国产午夜精品一二区理论片| 日韩在线高清观看一区二区三区| 搞女人的毛片| 人人妻人人看人人澡| 美女被艹到高潮喷水动态| 亚洲国产日韩欧美精品在线观看| 不卡视频在线观看欧美| 99视频精品全部免费 在线| 久久久久久九九精品二区国产| 国产午夜精品久久久久久一区二区三区| av免费观看日本| 日韩大片免费观看网站 | 禁无遮挡网站| 能在线免费观看的黄片| 久热久热在线精品观看| 亚洲国产精品久久男人天堂| 啦啦啦观看免费观看视频高清| 日日摸夜夜添夜夜添av毛片| 一个人看的www免费观看视频| 麻豆久久精品国产亚洲av| 村上凉子中文字幕在线| 美女国产视频在线观看| 国产黄a三级三级三级人| 成人漫画全彩无遮挡| av国产久精品久网站免费入址| 午夜爱爱视频在线播放| 久久草成人影院| 久久久久久久亚洲中文字幕| 国产亚洲精品av在线| 三级经典国产精品| 22中文网久久字幕| 亚洲成人中文字幕在线播放| 亚洲美女视频黄频| 99在线人妻在线中文字幕| 狠狠狠狠99中文字幕| 最近手机中文字幕大全| 免费在线观看成人毛片| 国产精品永久免费网站| 久久久欧美国产精品| .国产精品久久| 久久欧美精品欧美久久欧美| 国产伦精品一区二区三区视频9| 国产一区二区在线av高清观看| 亚洲av熟女| 国产淫片久久久久久久久| 18禁动态无遮挡网站| 精品国内亚洲2022精品成人| 国产探花极品一区二区| 三级国产精品欧美在线观看| 91av网一区二区| 久久人人爽人人爽人人片va| 亚洲欧美日韩卡通动漫| 国产淫片久久久久久久久| 久久久久久久国产电影| av专区在线播放| av天堂中文字幕网| 欧美又色又爽又黄视频| 色吧在线观看| 校园人妻丝袜中文字幕| 久久久久网色| 蜜桃久久精品国产亚洲av| 国产在线男女| 插逼视频在线观看| 国产 一区精品| 精品久久久久久久久av| 日日干狠狠操夜夜爽| 午夜免费激情av| 如何舔出高潮| 波多野结衣巨乳人妻| 国产老妇女一区| 深夜a级毛片| 国产亚洲午夜精品一区二区久久 | 天堂av国产一区二区熟女人妻| 三级经典国产精品| 又粗又爽又猛毛片免费看| 亚洲伊人久久精品综合 | 亚洲国产精品成人久久小说| 亚洲精品乱码久久久v下载方式| 亚洲人成网站在线观看播放| 精品国产露脸久久av麻豆 | 日日摸夜夜添夜夜添av毛片| 国产免费福利视频在线观看| eeuss影院久久| 舔av片在线| 男人的好看免费观看在线视频| 久久久久网色| 免费黄色在线免费观看| 高清毛片免费看| 国产伦在线观看视频一区| 我的老师免费观看完整版| 男女下面进入的视频免费午夜| or卡值多少钱| 1024手机看黄色片| 精品99又大又爽又粗少妇毛片| 午夜视频国产福利| 全区人妻精品视频| 久久亚洲精品不卡| 在线观看美女被高潮喷水网站| 精品人妻熟女av久视频| 国产av一区在线观看免费| 国产伦理片在线播放av一区| 1000部很黄的大片| 高清日韩中文字幕在线| 欧美zozozo另类| 欧美3d第一页| 一级爰片在线观看| 老司机影院毛片| 少妇被粗大猛烈的视频| 久久综合国产亚洲精品| 欧美成人免费av一区二区三区| 九九热线精品视视频播放| 欧美xxxx性猛交bbbb| 国产大屁股一区二区在线视频| 国产精品乱码一区二三区的特点| 男人舔奶头视频| 99在线人妻在线中文字幕| 成人无遮挡网站| 内地一区二区视频在线| 国产成人freesex在线| 亚洲真实伦在线观看| 国产成人福利小说| 99久久精品国产国产毛片| 直男gayav资源| 欧美激情久久久久久爽电影| 2022亚洲国产成人精品| 久久久久久国产a免费观看| 大香蕉97超碰在线| 国产又黄又爽又无遮挡在线| 网址你懂的国产日韩在线| 国产一区二区在线观看日韩| 一级毛片久久久久久久久女| 三级国产精品片| 97超碰精品成人国产| 久久精品人妻少妇| 少妇的逼好多水| 1024手机看黄色片| 国产免费又黄又爽又色| 欧美日韩一区二区视频在线观看视频在线 | 黄色日韩在线| 日韩精品有码人妻一区| 久久6这里有精品| 特大巨黑吊av在线直播| 不卡视频在线观看欧美| 精品人妻视频免费看| 亚洲最大成人av| 麻豆乱淫一区二区| 国产精品人妻久久久久久| 亚洲精品乱码久久久久久按摩| 国产精品三级大全| 国产精品美女特级片免费视频播放器| 国产淫语在线视频| 男的添女的下面高潮视频| 日本猛色少妇xxxxx猛交久久| 日韩亚洲欧美综合| 欧美97在线视频| av在线蜜桃| 精品久久国产蜜桃| 1024手机看黄色片| 国产精品电影一区二区三区| 91久久精品电影网| 久久久久久久午夜电影| 非洲黑人性xxxx精品又粗又长| 国产又色又爽无遮挡免| 亚洲真实伦在线观看| 美女xxoo啪啪120秒动态图| 欧美丝袜亚洲另类| 欧美变态另类bdsm刘玥| 国产三级中文精品| 国产精品爽爽va在线观看网站| 免费黄色在线免费观看| 国产精品不卡视频一区二区| 日本与韩国留学比较| 在线观看美女被高潮喷水网站| 亚洲精品456在线播放app| 久久久久国产网址| 婷婷六月久久综合丁香| 午夜老司机福利剧场| 国产一区有黄有色的免费视频 | 亚洲天堂国产精品一区在线| 最近中文字幕2019免费版| 在线免费观看的www视频| 在线免费观看不下载黄p国产| 国内精品一区二区在线观看| av免费观看日本| 久久99热这里只频精品6学生 | 日韩在线高清观看一区二区三区| 自拍偷自拍亚洲精品老妇| 日本爱情动作片www.在线观看| 亚洲av不卡在线观看| 国产精品.久久久| 免费电影在线观看免费观看| 1024手机看黄色片| 亚洲精华国产精华液的使用体验| 亚洲自拍偷在线| 青青草视频在线视频观看| 欧美极品一区二区三区四区| 欧美日韩综合久久久久久| 97超视频在线观看视频| 精品人妻熟女av久视频| 中文资源天堂在线| 男人舔女人下体高潮全视频| 最后的刺客免费高清国语| 啦啦啦韩国在线观看视频| 亚洲激情五月婷婷啪啪| 18禁在线播放成人免费| 成年av动漫网址| 永久网站在线| 嫩草影院精品99| 久久久久久久国产电影| 中文亚洲av片在线观看爽| 精品久久久久久久久av| 成人鲁丝片一二三区免费| 久久久久久久久久久免费av| 黑人高潮一二区| 国产精品国产三级国产专区5o | 欧美性猛交╳xxx乱大交人| 午夜福利成人在线免费观看| 国产熟女欧美一区二区| 97在线视频观看| 蜜桃亚洲精品一区二区三区| 亚洲成人av在线免费| 亚洲欧洲日产国产| 免费看日本二区| 一级黄色大片毛片| 国产精品三级大全| 男人狂女人下面高潮的视频| 色网站视频免费| 麻豆国产97在线/欧美| 一级av片app| 国产精品久久久久久久电影| 尾随美女入室| 亚洲美女视频黄频| 国产精品久久久久久精品电影小说 | 能在线免费看毛片的网站| 亚洲三级黄色毛片| 久久国产乱子免费精品| 男插女下体视频免费在线播放| 欧美变态另类bdsm刘玥| 欧美bdsm另类| 久久久久网色| 国产成人freesex在线| 亚州av有码| 干丝袜人妻中文字幕| 国产熟女欧美一区二区| 噜噜噜噜噜久久久久久91| 99久久人妻综合| 91午夜精品亚洲一区二区三区| 岛国在线免费视频观看| 啦啦啦观看免费观看视频高清| 国产精品无大码| 亚洲四区av| 国产免费男女视频| 亚洲精品色激情综合| 啦啦啦啦在线视频资源| 精品国产露脸久久av麻豆 | 欧美精品一区二区大全| 狠狠狠狠99中文字幕| 欧美97在线视频| 少妇被粗大猛烈的视频| 精品一区二区三区视频在线| 天天一区二区日本电影三级| 国产亚洲91精品色在线| 国产精品久久久久久久久免| 国产免费福利视频在线观看| 一个人免费在线观看电影| 亚洲av成人精品一区久久| 国产 一区精品| 精品一区二区三区视频在线| 日本熟妇午夜| 国产一区二区在线av高清观看| 一卡2卡三卡四卡精品乱码亚洲| 国产爱豆传媒在线观看| 男女视频在线观看网站免费| 国产午夜精品一二区理论片| 国产成人91sexporn| 美女高潮的动态| 三级国产精品欧美在线观看| 成人午夜高清在线视频| 天堂av国产一区二区熟女人妻| 美女国产视频在线观看| 亚洲av电影不卡..在线观看| 久久久午夜欧美精品| 在现免费观看毛片| 国产精品国产高清国产av| 色综合站精品国产| 乱码一卡2卡4卡精品| av免费观看日本| 亚洲欧美一区二区三区国产| 在线观看66精品国产| 我要搜黄色片| 国产欧美日韩精品一区二区| 色综合亚洲欧美另类图片| 啦啦啦啦在线视频资源| 中文天堂在线官网| 成人无遮挡网站| 亚洲欧美日韩东京热| 久久久精品94久久精品| 国产色爽女视频免费观看| 国产不卡一卡二| 国产视频内射| 边亲边吃奶的免费视频| 日韩高清综合在线| 亚洲国产精品合色在线| 亚洲精品国产成人久久av| 国产高清三级在线| 国内揄拍国产精品人妻在线| 成年免费大片在线观看| 边亲边吃奶的免费视频| 一级毛片我不卡| 白带黄色成豆腐渣| 超碰av人人做人人爽久久| 日韩一区二区视频免费看| 美女黄网站色视频| 亚洲不卡免费看| 久久国产乱子免费精品| 最近2019中文字幕mv第一页| 午夜福利在线观看免费完整高清在| av卡一久久| 伦精品一区二区三区| a级毛片免费高清观看在线播放| 午夜久久久久精精品| 午夜a级毛片| 嘟嘟电影网在线观看| 婷婷色麻豆天堂久久 | 在线a可以看的网站| 日日干狠狠操夜夜爽| 看黄色毛片网站| 床上黄色一级片| 国产日韩欧美在线精品| 草草在线视频免费看| 欧美成人免费av一区二区三区| 2021天堂中文幕一二区在线观| 成人av在线播放网站| 国产免费福利视频在线观看| 国模一区二区三区四区视频| 女人久久www免费人成看片 | 国产又色又爽无遮挡免| 国产一区二区在线av高清观看| 久久久精品大字幕| 在线免费观看的www视频| 久久久久久久久久久丰满| 哪个播放器可以免费观看大片| 国产精品国产三级国产专区5o | 国产亚洲精品久久久com| 国产老妇女一区| 午夜福利在线观看吧| 亚洲av二区三区四区| 日韩av在线大香蕉| www.av在线官网国产| 日本爱情动作片www.在线观看| 色5月婷婷丁香| 美女xxoo啪啪120秒动态图| 国产精品人妻久久久久久| 成人特级av手机在线观看| 麻豆一二三区av精品| 男女边吃奶边做爰视频| 国产探花在线观看一区二区| 精品午夜福利在线看| 男女那种视频在线观看| 日韩大片免费观看网站 | 伦理电影大哥的女人| 成人午夜高清在线视频| 国产毛片a区久久久久| 欧美zozozo另类| 白带黄色成豆腐渣| 国产成人a∨麻豆精品| 亚洲成色77777| 亚洲国产最新在线播放| 精品国产露脸久久av麻豆 | www.色视频.com| 国产久久久一区二区三区| 久久这里只有精品中国| 国产在线一区二区三区精 | 一区二区三区四区激情视频| 成人高潮视频无遮挡免费网站| 亚洲欧美日韩卡通动漫| 国产黄片美女视频| 18禁在线无遮挡免费观看视频| 亚洲av二区三区四区| 成人无遮挡网站| 久久精品91蜜桃| 1024手机看黄色片| 91久久精品电影网| 亚洲人与动物交配视频| 日产精品乱码卡一卡2卡三| 韩国高清视频一区二区三区| 国产一区二区三区av在线| 国产精品蜜桃在线观看| 精品久久久久久久久亚洲|