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

    TSP斷層模型數(shù)值模擬

    2015-03-06 11:50:55劉爭平王朝令
    關(guān)鍵詞:反射層云圖斷層

    林 義,劉爭平,王朝令,肖 締

    1.西南交通大學(xué)地球科學(xué)與環(huán)境工程學(xué)院,成都 610031 2.中鐵西北院科學(xué)研究院有限公司,蘭州 730000 3.四川農(nóng)業(yè)大學(xué)土木工程學(xué)院,成都 611830

    ?

    TSP斷層模型數(shù)值模擬

    林 義1,2,劉爭平1,王朝令3,肖 締3

    1.西南交通大學(xué)地球科學(xué)與環(huán)境工程學(xué)院,成都 610031 2.中鐵西北院科學(xué)研究院有限公司,蘭州 730000 3.四川農(nóng)業(yè)大學(xué)土木工程學(xué)院,成都 611830

    隧道在施工開挖中會(huì)遇到各種地質(zhì)問題,其中以斷層和軟弱帶居多,目前隧道地質(zhì)預(yù)報(bào)主要采用TSP(tunnel seismic prediction)系統(tǒng)進(jìn)行。雖然TSP技術(shù)應(yīng)用廣泛,但目前對它的研究工作主要集中于工程應(yīng)用實(shí)例,采用正演模擬方法進(jìn)行的研究較少。筆者采用有限元方法模擬隧道地震波場,采用波場快照與時(shí)間記錄相結(jié)合的方法研究斷層對隧道地震波場傳播的影響,并對含斷層模型的時(shí)間記錄進(jìn)行了反演處理,得到了數(shù)值模型的速度云圖和反射層位圖。數(shù)據(jù)處理結(jié)果表明:采用TSP Win軟件默認(rèn)值處理得到的速度云圖與模型設(shè)定的斷層位置一致;根據(jù)反射層位圖,對異常速度帶的層狀模型來說,P波預(yù)報(bào)的準(zhǔn)確性更高。研究表明,TSP系統(tǒng)具有良好的抗噪性能。通過對工程實(shí)例的處理,驗(yàn)證了數(shù)值模擬所得結(jié)論。

    TSP;TSPwin;波場快照;隧道地震預(yù)報(bào);斷層

    0 引言

    TSP(tunnel seismic prediction)系統(tǒng)是瑞士Amberg公司開發(fā)的、進(jìn)行隧道地質(zhì)超前預(yù)報(bào)的地球物理勘探設(shè)備。研究原理是在隧道內(nèi)[1]采用小藥量的炸藥爆破作為震源,按照地震勘探中震源與檢波器的互換原理,在隧道的邊壁進(jìn)行一系列激發(fā),只有兩個(gè)檢波器接收,這樣可以實(shí)現(xiàn)單炮激發(fā)、多道接收的效果。由震源激發(fā)產(chǎn)生的地震波遇到地層界面、節(jié)理面,特別是斷層破碎帶、溶洞、暗河、巖溶陷落柱、淤泥帶等不良地質(zhì)界面時(shí),將產(chǎn)生反射波?;诘卣鹂碧皆?,檢波器所接收到的反射波傳播速度、延遲時(shí)間、波形、強(qiáng)度和方向等均與界面的性質(zhì)和產(chǎn)狀密切相關(guān)[2],并通過不同振幅、同相軸等波形特征表現(xiàn)出來。TSP主要用于超前預(yù)報(bào)隧道掌子面前方、上方和下方不良地質(zhì)體的性質(zhì)、位置和規(guī)模,最大預(yù)報(bào)距離300~500 m,最高可以分辨直徑為1 m的地質(zhì)體[3]。

    作為工程建設(shè)中的一種實(shí)用技術(shù),TSP具有追求短平快的特點(diǎn),前人已從工程應(yīng)用的角度對TSP進(jìn)行了研究和分析。劉志剛等[1]對TSP基本原理、技術(shù)狀況及其在隧道隧洞超前預(yù)報(bào)中能夠解決的技術(shù)問題、應(yīng)用效果、能達(dá)到的技術(shù)指標(biāo)以及發(fā)展進(jìn)行研究,論述了其存在的問題和解決的途徑。何剛[4]分析了分辨率、頻帶寬度及激發(fā)炸藥藥量之間的相互關(guān)系,提出了對TSP接收系統(tǒng)的改進(jìn)意見,并深入剖析系統(tǒng)測量的誤差因素。李堅(jiān)[5]應(yīng)用TSP對地層分界面、破碎帶、巖溶等進(jìn)行預(yù)報(bào),預(yù)報(bào)結(jié)果得到了后續(xù)施工開挖的驗(yàn)證。此外:Ashida[6]以隧道爆破和TBM(tunnel boring machine)震動(dòng)作為震源,對隧道掌子面前方的地質(zhì)情況做了成像研究,成像結(jié)果與前方實(shí)際地質(zhì)情況一致;Bohlen等[7]研究了隧道掌子面上RSSR現(xiàn)象,即隧道地震波場中的Rayleigh面波在掌子面上轉(zhuǎn)換為S波,遇到反射界面后反射到掌子面上再轉(zhuǎn)換為Rayleigh面波;王朝令等[8]基于有限元法對隧道地震波場進(jìn)行了模擬,在其中討論了頻散、邊界條件等環(huán)節(jié)。

    筆者應(yīng)用數(shù)值模擬方法模擬存在斷層的TSP地震波場,提取其時(shí)間記錄和波場快照,分析存在斷層的情況下的地震波場特征;再應(yīng)用TSP后處理軟件TSPwin對得到的時(shí)間記錄進(jìn)行處理,對比模型圖與提取到反射界面的差異,并研究TSP系統(tǒng)的抗噪性能;最后通過工程實(shí)例驗(yàn)證了研究結(jié)果,以期為TSP隧道地震預(yù)報(bào)提供研究依據(jù)。

    1 TSP系統(tǒng)

    TSPwin是TSP系統(tǒng)的數(shù)據(jù)處理軟件。處理的結(jié)果是形成反映相關(guān)界面或地質(zhì)體反射能量的影像點(diǎn)圖和隧道平面、剖面圖,以供工程技術(shù)人員解譯。數(shù)據(jù)處理流程如圖1所示,主要處理步驟包括帶通濾波、初至拾取、Radon變換、Q值濾波、速度分析、深度偏移和反射界面提取。

    圖1 TSP數(shù)據(jù)處理流程圖Fig.1 Flow chart of TSP data processing

    ①帶通濾波:通過Fourier變換得到地震數(shù)據(jù)的頻譜,然后根據(jù)有效波的頻率范圍設(shè)置濾波窗口,將頻率在窗口之外的干擾波濾除。

    ②初至拾取:提取直達(dá)波到達(dá)檢波器的起跳時(shí)刻點(diǎn)。

    ③Radon變換:Radon變換是將t-x域地震記錄轉(zhuǎn)變?yōu)棣?p/q域數(shù)據(jù)的一種數(shù)學(xué)變換[9]。τ-p/q變換后各種波的運(yùn)動(dòng)學(xué)特征和動(dòng)力學(xué)特征會(huì)發(fā)生變化,使得某些問題在τ-p/q域中解決更為方便[10]。

    ④Q值濾波[11]:Q為衡量地震波在傳播過程中能量衰減的品質(zhì)因子,基于它可以分析傳播介質(zhì)的非彈性情況。

    ⑤速度分析:速度是地震勘探的關(guān)鍵參數(shù),速度分析的主要目的是為水平疊加、偏移等處理步驟提供處理參數(shù)。

    ⑥深度偏移[12]:深度偏移是將檢波器接收到的時(shí)間剖面轉(zhuǎn)換成以空間坐標(biāo)表示的地質(zhì)構(gòu)造剖面。

    ⑦反射界面提?。禾崛〉椒瓷浣缑娴臋E圓形曲線代表界面的反射振幅。其明顯程度表示物理反射界面的波阻抗強(qiáng)度對比,反射振幅越明顯,則反射系數(shù)越大,巖石介質(zhì)的彈性波阻抗分界面屬性差異也就越大。

    2 TSP地震波場數(shù)值模擬

    為了模擬松動(dòng)圈存在時(shí)的TSP地震波場,建立如圖2所示的模型圖。震源采用300 Hz的Ricker子波,這個(gè)頻率與TSP實(shí)際工程數(shù)據(jù)采集的震源主頻相近。通常地震子波是由2至3個(gè)或多個(gè)相位組成的地震脈沖,最常用的Ricker子波由1個(gè)波峰和2個(gè)波谷組成[10]。在本文的數(shù)值模擬研究中,采用0相位的Ricker子波,其表示公式如下:

    (1)

    式中:u表示位移;f0表示主頻;t代表時(shí)刻;t0代表中心時(shí)刻。

    TSP數(shù)據(jù)采集時(shí),采用單個(gè)三分量接收,24炮激發(fā);根據(jù)炮檢互換原理,相當(dāng)于單炮激發(fā),24道三分量檢波器接收。為了方便計(jì)算,在設(shè)置模型時(shí),即采用單炮激發(fā)多道接收,使得計(jì)算模型的觀測系統(tǒng)與TSP實(shí)際采集時(shí)相同。需要說明的是,在模型數(shù)據(jù)處理過程中,速度云圖和反射層圖都是從震源位置起算。

    模型(圖2)周邊加載了黏彈性邊界,它可以吸收邊界反射波[13]。模型的長度為135 m,高度為75 m,隧道高度6 m,有限元網(wǎng)格劃分為0.5 m×0.5 m,偏移距15 m,道間距1 m,采樣間隔62.5 μs,采樣點(diǎn)數(shù)7 218,采樣時(shí)窗451.125 ms。由于隧道開挖過程中存在松動(dòng)圈[14],因此模擬時(shí)在隧道邊壁中設(shè)置了斷層帶,其厚度為1 m。工作面前面斷層帶總共設(shè)置3層,層厚3 m,層間距15 m,從震源到第一、二、三層異常帶的水平距離分別為65、80、95 m,每個(gè)區(qū)域內(nèi)的縱波速度(vP)、橫波速度(vS)以及密度(ρ)如圖2所示。

    圖2 模型圖Fig.2 Model diagram

    作為地震勘探方法在隧道預(yù)報(bào)中的一種應(yīng)用,TSP的分辨率遵循地震勘探分辨率的Knapp準(zhǔn)則,TSP能區(qū)分的厚度是波長的1/2。本文模型中縱波的波長為4 000/300≈13.3,故TSP能區(qū)分6.65 m的厚度。模型中雖然斷層的厚度為3 m,但斷層間距為15 m,故可以分辨。模型中斷層只是作為一個(gè)反射層進(jìn)行勘探,這是因?yàn)門SP探測的目的主要是勘探斷層的位置和數(shù)量,對斷層厚度可以不予區(qū)分。

    a.x分量;b.y分量。圖3 數(shù)值模擬的時(shí)間記錄Fig.3 Time recording of numerical simulation

    由于模擬是在二維中進(jìn)行的,故只能取到x、y兩個(gè)分量。原始時(shí)間剖面如圖3所示。圖3a中的同相軸依次為直達(dá)縱波、直達(dá)橫波、左邊界角點(diǎn)繞射橫波和反射橫波。由于反射縱波能量較小,在圖中沒有顯示。時(shí)間剖面x分量(圖3a)各同相軸起跳點(diǎn)明顯,y分量(圖3b)各同相軸一致性好且較清晰,信噪比和分辨率都比較高。需要指出的是RSSR現(xiàn)象,即隧道內(nèi)沿邊壁傳播的Rayleigh面波在工作面角點(diǎn)位置轉(zhuǎn)換為S波,以體波的形式向前傳播,遇到界面反射回到角點(diǎn)后轉(zhuǎn)換為Rayleigh面波的現(xiàn)象[7]。圖3a中前方界面的反射橫波經(jīng)RSSR轉(zhuǎn)換在隧道內(nèi)產(chǎn)生Rayleigh面波,由于其與橫波速度相差很小,使得60~80 ms時(shí)間段內(nèi)x分量剖面的反射波同相軸混亂,這影響了反射橫波的分辨率,也降低了地震解釋的可靠性。為了降低這種干擾,提高資料解釋的可靠性,觀測系統(tǒng)應(yīng)盡量采用大偏移距,使得RSSR現(xiàn)象在隧道內(nèi)產(chǎn)生Rayleigh面波與前方界面的反射橫波分離開來。

    圖4是模型激發(fā)y分量的波場快照。t=1 ms時(shí),震源開始激發(fā);t=5 ms時(shí)觀察到S波傳播;t=8、10 ms時(shí)波傳播到左邊界;t=12 ms時(shí)可以看到波被左邊界吸收;t=15、18 ms時(shí)繼續(xù)在接收排列上傳播;t=20、22 ms時(shí),波傳播到底邊界,底邊界將其吸收;t=25、28、30 ms時(shí),波傳播到工作面角點(diǎn)位置;t=32、35 ms時(shí),RSSR轉(zhuǎn)換出的S波向前方傳播,以角點(diǎn)為震源產(chǎn)生新的波陣面。此過程是波到達(dá)斷層之前的。

    圖4 模擬y分量的波場快照Fige.4 The y-component snapshot of simulation

    t=38 ms時(shí)波到達(dá)第一個(gè)反射層位置;t=40 ms時(shí)S波開始反射回工作面;t=42 ms時(shí)波通過第一個(gè)反射層,到達(dá)第一個(gè)反射層與第二個(gè)反射層之間的空間位置;t=45 ms時(shí),波傳播到達(dá)第二個(gè)反射層;t=48、50 ms時(shí),波穿過第二個(gè)反射層,此時(shí)從第一個(gè)反射層返回的波又傳播到排列上;t=52、55 ms時(shí),波傳播到第三個(gè)反射層;t=58、60 ms時(shí),穿過第三個(gè)反射層,傳播到右邊界上;t=62、65 ms時(shí),從第一個(gè)反射層返回的波已經(jīng)通過接收排列;t=68、70 ms,不斷有從第二、三反射層上的波向排列上傳播;t=72~80 ms,反射波完全穿過排列,此時(shí)仍有大量的雜波在快照上存在。

    3 模擬數(shù)據(jù)處理

    按照如圖1所示的處理流程應(yīng)用TSPwin對得到的時(shí)間記錄進(jìn)行處理,各處理步驟均采用TSPwin默認(rèn)值進(jìn)行,可以得到如圖5所示速度云圖:橫坐標(biāo)表示預(yù)報(bào)距離,計(jì)算以震源位置起算;縱坐標(biāo)表示以隧道向兩邊展布半徑。圖5a中P波速度云圖的低速帶存在位置為75~110 m,模型(圖2)設(shè)定的斷層帶位置為65~100 m,兩者非常接近。圖5b中S波速度云圖低速帶在80~120 m,與模型(圖2)設(shè)置的斷層帶也比較接近;但所呈現(xiàn)的形態(tài)一致性不高,分析是受到RSSR轉(zhuǎn)換的干擾,使得速度分布云圖出現(xiàn)偏差。綜合兩張圖可知,P波異常帶的速度分布無論是位置還是形態(tài)都非常接近模型設(shè)定的斷層帶區(qū)域,表明采用TSPwin默認(rèn)值處理所得到的P波速度分布云圖與模型設(shè)定非常接近。

    圖5 速度云圖Fig.5 Velocity Contours

    圖6是提取得到的反射層圖。P波反射界面所處位置與模型設(shè)置非常接近,且形態(tài)也一致,預(yù)報(bào)的反射界面正好處于隧道軸線上(圖6a),這正是隧道超前預(yù)報(bào)所需要解決的。S波提取反射界面(圖6b)與模型設(shè)定位置比較接近,但形態(tài)接近程度相對P波較差。對隧道地震超前預(yù)報(bào)來說,用P波進(jìn)行預(yù)報(bào)的準(zhǔn)確性比較高,而S波提取到的反射界面準(zhǔn)確性相對稍差。

    對于如本模型這樣的多層異常帶模型來說,S波沒有P波準(zhǔn)確,P波的速度分布和反射界面都與模型設(shè)定的位置非常吻合;對隧道地震超前預(yù)報(bào)來說,P波預(yù)報(bào)的準(zhǔn)確性更高。

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

    4 TSP抗噪性分析

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

    對加噪聲信號(hào)采用TSPwin進(jìn)行處理,經(jīng)過帶通濾波、初至拾取、反射波提取、P-S波分離、深度偏移和反射層提取,得到圖8所示的速度云圖和圖9的反射層提取圖。

    圖8所示為加噪聲的速度云圖,對比模型(圖2)可以發(fā)現(xiàn),在65 ~100 m區(qū)間內(nèi),P波的速度云圖存在低速區(qū),這與模型設(shè)定相近;相比之下,S波的速度云圖在60~80 m范圍內(nèi)也存在低速區(qū),但與模型的對應(yīng)程度較差。

    圖9為加噪聲數(shù)據(jù)所提取的反射層,對比模型(圖2)可以看到:在70~100 m區(qū)間內(nèi), P波反射界面所處位置與模型設(shè)置相近,且形態(tài)也較一致,但是相比于圖8未加噪聲的P波提取反射層,對應(yīng)性有所降低,且反射層強(qiáng)度也有所下降;S波在80~100 m存在反射層,這接近模型的設(shè)定位置,與圖8未加噪聲的S波提取反射層相比,一致性和對應(yīng)性相近。

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

    a.x分量;b.y分量。圖7 加噪聲的時(shí)間記錄Fig.7 Time recording with additive noise

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

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

    5 實(shí)例

    考慮到隧道超前預(yù)報(bào)方法的施工特點(diǎn),若要驗(yàn)證預(yù)報(bào)結(jié)果的準(zhǔn)確性,只能選取已經(jīng)開挖完的隧道,比較開挖之前的預(yù)報(bào)結(jié)果和開挖后揭示的地質(zhì)情況,對比兩者的差異。依據(jù)常規(guī)TSP觀測排列的布置方法,在隧道邊墻布置地震波信息接收孔,孔深2 m;按約1.5 m的間距布置24個(gè)激發(fā)孔,孔深1.5 m左右,向下傾斜約10°;每個(gè)激發(fā)孔裝填藥量為100 g,激發(fā)孔和接收孔基本保持在同一個(gè)高度。

    圖10 實(shí)例P波速度云圖(a)和反射層圖(b)Fig.10 Real case of P-wave velocity contours(a) and layers(b) of extracting reflection

    圖10是實(shí)例P波的速度云圖和反射層圖。圖10a中,110~140 m存在異常低速區(qū)域。由于是實(shí)例預(yù)報(bào),地震波場比較復(fù)雜,因此提取的反射層比較多,分布位置比較分散。其中,在110~140 m,存在比較多的反射層,形態(tài)以豎向反射層為主(圖10b)。它表示本段圍巖節(jié)理裂隙發(fā)育,圍巖破碎,推測該段局部存在節(jié)理密集帶或破碎帶。開挖后,120~135 m段中間偏右存在黃黏土夾層,在拱腰至拱頂向內(nèi)延伸約1.5~3.0 m,是一處破碎帶構(gòu)造,與預(yù)報(bào)結(jié)果相吻合。

    6 結(jié)論與建議

    筆者建立斷層模型進(jìn)行數(shù)值模擬,通過波場快照研究隧道地震波場的傳播特征,對得到的時(shí)間剖面按照TSP的常規(guī)流程進(jìn)行處理。得到如下結(jié)論:

    1)在速度分布云圖中,斷層的速度呈現(xiàn)比周圍速度低的特征,在提取的反射層圖中斷層呈現(xiàn)豎向條狀的特征。

    2)對數(shù)值模擬所得的時(shí)間記錄添加噪聲,結(jié)果表明TSP具有良好的抗噪性能。

    3)工程實(shí)例表明P波預(yù)報(bào)是比較準(zhǔn)確的,為工程中應(yīng)用TSP進(jìn)行隧道地震預(yù)報(bào)的數(shù)據(jù)處理和資料解釋提供理論依據(jù)。

    由于隧道施工對TSP時(shí)效性的要求較高,采集環(huán)境復(fù)雜,且地質(zhì)體多樣,文中的結(jié)論是基于垂直多斷層的情況,因此需要研究不同情況下的TSP預(yù)報(bào)效果,并開展三維隧道地震波場數(shù)值模擬下的TSP預(yù)報(bào)效果研究,這是需要進(jìn)一步研究的工作。

    [1] 劉志剛,劉秀峰.TSP (隧道地震勘探) 在隧道隧洞超前預(yù)報(bào)中的應(yīng)用與發(fā)展[J].巖石力學(xué)與工程學(xué)報(bào),2003,22(8):1399-1402. Liu Zhigang,Liu Xiufeng. TSP Application and Deve-lopment in Tunnel Lead Forecast[J]. Chinese Journal of Rock Mechanics and Engineering,2003,22(8):1399-1402.

    [2] Alimoradi A, Moradzadeh A, Naderi R, et al.Prediction of Geological Hazardous Zones in Front of a Tunnel Face Using TSP-203 and Artificial Neural Networks[J]. Tunnelling and Underground Space Technology, 2008, 23(6): 711-717.

    [3] 楊峰.TSP隧道超前預(yù)報(bào)中典型地質(zhì)異常體的射線追蹤數(shù)值模擬研究[D].成都:西南交通大學(xué),2007. Yang Feng. The Ray Tracing Numerical Simulation Study of Typical Geological Anomalous Body in Tunnel Advance Prediction[D]. Chengdu: Southwest Jiaotong University, 2007.

    [4] 何剛.TSP-203 系統(tǒng)在隧道超前地質(zhì)預(yù)報(bào)中的應(yīng)用研究[D].長沙:中南大學(xué),2005. He Gang. The Research of Application About TSP-203 in Tunnel Geological Predicatio[D]. Changsha:Central South University, 2005.

    [5] 李堅(jiān).TSP法在鐵路客運(yùn)專線隧道超前地質(zhì)預(yù)報(bào)工作中的應(yīng)用前景[J].鐵道勘察,2006,31(6):45-49. Li Jian. Application Prospect of the TSP Method in the Lead Geological Predication of the Dedicated-Passenger Railway Tunnels[J]. Railway Survey, 2006, 31(6): 45-49.

    [6] Ashida Y. Seismic Imaging Ahead of a Tunnel Face with Three-Component Geophones[J]. International Journal of Rock Mechanics and Mining Sciences, 2001,38(6): 823-831.

    [7] Bohlen T, Ulrich L,Wolfgang R, et al. Rayleigh-to-Shear Wave Conversion at the Tunnel Face-From 3D-FD Modeling to Ahead-of-Drill Exploration[J].Geophysics, 2007, 72(6): 67-79.

    [8] 王朝令, 劉爭平, 黃云艷,等.隧道地震預(yù)報(bào)波場的有限元數(shù)值模擬[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版,2014,44(4): 1369-1381. Wang Zhaoling,Liu Zhengping, Huang Yunyan, et al. Wavefield Modeling Based on the Finite Element Method for the Tunnel Seismic Prediction[J]. Journal of Jilin University: Earth Science Editon, 2014,44(4):1369-1381.

    [9] 王朝令,劉爭平.τ-p變換在隧道反射地震超前預(yù)報(bào)波場分離中應(yīng)用的數(shù)值模擬研究[J].地球物理學(xué)進(jìn)展,2012,27(5): 2216-2225. Wang Zhaoling, Liu Zhengping. The Modelling Study ofτ-pTransform Applied in the Tunnel Seismic Prediction[J]. Progress in Geophysics, 2012, 27(5):2216-2225.

    [10] 王朝令.隧道地震超前預(yù)報(bào)中波場分離與反演方法的數(shù)值模擬研究[J].成都:西南交通大學(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.

    [11] 徐彥, 蘇有錦,秦嘉政.Q值研究動(dòng)態(tài)[J].地震研究,2004,27(7):385-389. Xu Yan, Su Youjin, Qin Jiazheng. Development ofQ-Value Studies[J]. Journal of Seismological Research, 2004, 27(7): 385-389.

    [12] 王有新, 史政軍,李少英.深度偏移的發(fā)展與應(yīng)用概況[J].石油地球物理勘探,1994,29(4): 526-531,534. Wang Youxin, Shi Zhengjun, Li Shaoying. Development of Depth Migration and Its Application[J]. Oil Geophysical Prospecting, 1994.29(4):526-531,534.

    [13] 王朝令,劉爭平.黏彈性邊界條件在ANSYS有限元波場模擬中的實(shí)現(xiàn)[J].大地測量與地球動(dòng)力學(xué),2012,32(2): 28-31. Wang Zhaoling, Liu Zhengping. Realization of Vis-coelastic Boundary Condition in Wave Field Simulation with Ansys Finite Element Software[J]. Journal of Geodesy and Geodynamics, 2012, 32(2): 28-31.

    [14] 宋宏偉, 王闖,賈穎絢.用地質(zhì)雷達(dá)測試圍巖松動(dòng)圈的原理與實(shí)踐[J].中國礦業(yè)大學(xué)學(xué)報(bào),2002,31(4):43-46. Song Hongwei, Wang Chuang, Jia Yingxuan. Priciple of Measuring Broken Rock Zone Around Underground Roadway with GPR and Its Application[J]. Journal of China University of Mining& Technology, 2002, 31(4): 43-46.

    Numerical Simulation of TSP Fault Model

    Lin Yi1,2, Liu Zhengping1, Wang Zhaoling3,Xiao Di3

    1.FacultyofGeosciencesandEnvironmentalEngineering,SouthwestJiaotongUniversity,Chengdu610031,China2.NorthwestResearchInstituteCo,.LtdofC.R.E.C,Lanzhou730000,China3.CollegeofCivilEngineering,SichuanAgriculturalUniversity,Chengdu611830,China

    During tunnel excavation, a variety of geological disasters might be encountered, such as faults, caves, et.al. Tunnel seismic prediction (TSP) is adopted to mitigate the possible damages. Although TSP technology is used widely, the research about TSP is currently focused on its engineering application cases. We use the finite element method to simulate the tunnel seismic wave field, employ wave field snapshots and time recording method on the impact of faults on the characteristics of the propagation of tunnel seismic wave field, and inversely process the time record of model containing the fault. The digital model of the velocity scattered image and the reflection interface position are obtained, and the fault position from velocity scattered image processed with the default values set by using TSPwin is agreed to the one from the model. In respect to the layered model for an abnormal velocity zone, P-wave is more precise. The system of TSP is strong for its feature of anti-noise. The numerical simulation is verified finally through the process and analysis of the engineering cases.

    TSP;TSPwin;snapshot;tunnel seismic prediction;fault

    10.13278/j.cnki.jjuese.201506305.

    2015-02-28

    國家自然科學(xué)基金項(xiàng)目(40874051,51508358)

    林義(1980--),男,博士研究生,高級工程師,主要從事隧道地震超前預(yù)報(bào)研究,E-mail:766708597@qq.com

    王朝令(1980--),男,副教授,博士,主要從事隧道地震波場模擬研究,E-mail:wong8010@gmail.com。

    10.13278/j.cnki.jjuese.201506305

    P631.4

    A

    林義,劉爭平,王朝令,等. TSP斷層模型數(shù)值模擬.吉林大學(xué)學(xué)報(bào):地球科學(xué)版,2015,45(6):1870-1878.

    Lin Yi, Liu Zhengping, Wang Zhaoling, et al. Numerical Simulation of TSP Fault Model.Journal of Jilin University:Earth Science Edition,2015,45(6):1870-1878.doi:10.13278/j.cnki.jjuese.201506305.

    猜你喜歡
    反射層云圖斷層
    反射層對NaI(Tl)閃爍體探測器探測效率的影響
    科技視界(2021年10期)2021-05-20 01:46:42
    基于有限元某壓水堆反射層組件模態(tài)分析的建模方法
    四川建筑(2020年3期)2020-07-18 01:56:24
    成都云圖控股股份有限公司
    黃強(qiáng)先生作品《雨后松云圖》
    名家名作(2017年3期)2017-09-15 11:13:37
    基于TV-L1分解的紅外云圖超分辨率算法
    隨鉆雙感應(yīng)測井儀背景影響研究
    云圖青石板
    斷層破碎帶壓裂注漿加固技術(shù)
    河南科技(2014年18期)2014-02-27 14:14:52
    關(guān)于錨注技術(shù)在煤巷掘進(jìn)過斷層的應(yīng)用思考
    河南科技(2014年7期)2014-02-27 14:11:06
    核電廠SCIENCE程序包反射層功能擴(kuò)展研究
    中國核電(2014年3期)2014-02-18 03:12:35
    免费在线观看影片大全网站| 精品午夜福利视频在线观看一区| 精品久久久久久久人妻蜜臀av| 女同久久另类99精品国产91| 亚洲国产看品久久| 欧美日韩乱码在线| 日本撒尿小便嘘嘘汇集6| 日韩欧美国产在线观看| 日韩中文字幕欧美一区二区| 精品一区二区三区av网在线观看| 欧美日韩精品网址| av欧美777| 91av网站免费观看| 国产单亲对白刺激| 香蕉久久夜色| 美女免费视频网站| 日日干狠狠操夜夜爽| xxxwww97欧美| 亚洲国产日韩欧美精品在线观看 | 欧美性猛交╳xxx乱大交人| 大香蕉久久成人网| 国产在线精品亚洲第一网站| 在线观看免费视频日本深夜| 国产精品免费一区二区三区在线| 午夜精品在线福利| 欧美一级a爱片免费观看看 | 久久久久久久久久黄片| 男人的好看免费观看在线视频 | 精品国产国语对白av| 日本成人三级电影网站| 国产私拍福利视频在线观看| 欧美色视频一区免费| 最近最新中文字幕大全电影3 | videosex国产| 久久这里只有精品19| 中文字幕久久专区| www.自偷自拍.com| 中文字幕另类日韩欧美亚洲嫩草| 特大巨黑吊av在线直播 | 免费人成视频x8x8入口观看| 亚洲国产精品久久男人天堂| 又紧又爽又黄一区二区| 丰满的人妻完整版| 制服诱惑二区| 波多野结衣巨乳人妻| 亚洲成人久久性| 亚洲第一欧美日韩一区二区三区| 日韩精品中文字幕看吧| 久久99热这里只有精品18| 久久久久久久精品吃奶| 中文在线观看免费www的网站 | 两个人免费观看高清视频| 18禁美女被吸乳视频| 欧美丝袜亚洲另类 | 男女那种视频在线观看| 久久婷婷成人综合色麻豆| 99热只有精品国产| 操出白浆在线播放| 亚洲精品av麻豆狂野| 每晚都被弄得嗷嗷叫到高潮| 亚洲人成77777在线视频| 女警被强在线播放| 黄色毛片三级朝国网站| 午夜影院日韩av| 黄色视频,在线免费观看| 国产97色在线日韩免费| 欧美性长视频在线观看| 国内揄拍国产精品人妻在线 | 黑人操中国人逼视频| 国产亚洲精品第一综合不卡| 久久久久久久精品吃奶| 黄色a级毛片大全视频| 亚洲一区二区三区色噜噜| 欧美黑人欧美精品刺激| 后天国语完整版免费观看| 一级毛片女人18水好多| 99国产极品粉嫩在线观看| 国产一区二区三区在线臀色熟女| 无人区码免费观看不卡| 波多野结衣巨乳人妻| 久久欧美精品欧美久久欧美| 国内毛片毛片毛片毛片毛片| 女人被狂操c到高潮| 欧美 亚洲 国产 日韩一| 久久国产乱子伦精品免费另类| 久久人妻福利社区极品人妻图片| 中出人妻视频一区二区| 欧美成人午夜精品| 成人免费观看视频高清| 特大巨黑吊av在线直播 | 午夜精品久久久久久毛片777| 欧美最黄视频在线播放免费| 精品人妻1区二区| 欧美性猛交╳xxx乱大交人| 国产成人欧美在线观看| 国产亚洲av高清不卡| 国产私拍福利视频在线观看| 中文字幕精品亚洲无线码一区 | 国产欧美日韩精品亚洲av| 久9热在线精品视频| 国产激情欧美一区二区| 欧美又色又爽又黄视频| 久久精品国产亚洲av高清一级| 国产亚洲欧美精品永久| 亚洲精品在线美女| 国内毛片毛片毛片毛片毛片| 久久午夜综合久久蜜桃| 男人舔奶头视频| 99国产综合亚洲精品| 国产成人欧美在线观看| 国产精品免费视频内射| 国产一区二区在线av高清观看| 三级毛片av免费| 好看av亚洲va欧美ⅴa在| 欧美在线黄色| 男人操女人黄网站| 黄色毛片三级朝国网站| 一a级毛片在线观看| 国产精品久久视频播放| 日韩大尺度精品在线看网址| 在线视频色国产色| 青草久久国产| 久久精品国产综合久久久| 97人妻精品一区二区三区麻豆 | 国产伦一二天堂av在线观看| 国产精品久久久人人做人人爽| 一区二区三区精品91| 精品国产乱子伦一区二区三区| 91字幕亚洲| 国内毛片毛片毛片毛片毛片| 两个人看的免费小视频| 欧美日本亚洲视频在线播放| 最近最新中文字幕大全免费视频| 精品久久久久久久久久免费视频| 男人舔女人的私密视频| 亚洲国产精品成人综合色| 999久久久精品免费观看国产| www.999成人在线观看| 后天国语完整版免费观看| 精品卡一卡二卡四卡免费| 一级a爱片免费观看的视频| 中亚洲国语对白在线视频| 观看免费一级毛片| 2021天堂中文幕一二区在线观 | 精品国产亚洲在线| 19禁男女啪啪无遮挡网站| 国产成人av教育| 一本大道久久a久久精品| 久久亚洲精品不卡| 免费在线观看成人毛片| 熟女电影av网| 精品久久久久久久毛片微露脸| 老司机福利观看| 麻豆成人av在线观看| www.999成人在线观看| 两人在一起打扑克的视频| 免费无遮挡裸体视频| 88av欧美| 欧美一级毛片孕妇| 少妇 在线观看| 深夜精品福利| 日本撒尿小便嘘嘘汇集6| 成年版毛片免费区| 淫秽高清视频在线观看| 精品久久久久久久久久久久久 | 又大又爽又粗| 国产伦一二天堂av在线观看| 免费av毛片视频| 亚洲成人久久爱视频| 搞女人的毛片| 久久久久久久久久黄片| 特大巨黑吊av在线直播 | 国产亚洲欧美98| 国产黄片美女视频| 国产亚洲欧美在线一区二区| 国产av一区在线观看免费| 黄色 视频免费看| 国产亚洲精品久久久久5区| 久久久久精品国产欧美久久久| 欧美亚洲日本最大视频资源| cao死你这个sao货| 成人一区二区视频在线观看| 国产主播在线观看一区二区| 成年免费大片在线观看| 久久婷婷成人综合色麻豆| 香蕉丝袜av| 久久久久精品国产欧美久久久| 久久精品人妻少妇| 久久精品成人免费网站| 中文字幕精品免费在线观看视频| 免费在线观看完整版高清| 97超级碰碰碰精品色视频在线观看| 国产精品99久久99久久久不卡| 长腿黑丝高跟| 欧美中文综合在线视频| 亚洲精品一区av在线观看| 一级片免费观看大全| 成人一区二区视频在线观看| 亚洲激情在线av| 老司机午夜十八禁免费视频| 国产成人欧美| a级毛片在线看网站| 国产精品电影一区二区三区| 精品国产美女av久久久久小说| 久99久视频精品免费| 手机成人av网站| 国产一区二区激情短视频| 国产亚洲欧美精品永久| 国产爱豆传媒在线观看 | 日韩欧美在线二视频| 大型黄色视频在线免费观看| 高清毛片免费观看视频网站| 日本熟妇午夜| 久久久国产精品麻豆| 国产久久久一区二区三区| 日本免费一区二区三区高清不卡| 国产成人精品无人区| 免费观看精品视频网站| 国产免费av片在线观看野外av| 中文字幕精品亚洲无线码一区 | 国产高清有码在线观看视频 | 欧美日本视频| 国产黄a三级三级三级人| 欧美成人免费av一区二区三区| 999精品在线视频| 搡老妇女老女人老熟妇| 国产精品野战在线观看| 久久精品国产亚洲av香蕉五月| 女警被强在线播放| 久久精品91蜜桃| 午夜亚洲福利在线播放| 亚洲国产欧美网| 国产一卡二卡三卡精品| 国产单亲对白刺激| 99国产极品粉嫩在线观看| a在线观看视频网站| 中文资源天堂在线| 国产一区二区在线av高清观看| svipshipincom国产片| 日韩精品青青久久久久久| 视频在线观看一区二区三区| 日韩精品中文字幕看吧| 免费观看人在逋| 精品国内亚洲2022精品成人| 日韩欧美一区二区三区在线观看| 精品午夜福利视频在线观看一区| 国产亚洲精品第一综合不卡| 热99re8久久精品国产| 国产精品亚洲一级av第二区| 俄罗斯特黄特色一大片| 波多野结衣巨乳人妻| 久久中文看片网| www国产在线视频色| 国产精品综合久久久久久久免费| a在线观看视频网站| 欧美 亚洲 国产 日韩一| 国产黄a三级三级三级人| 在线播放国产精品三级| 国产又色又爽无遮挡免费看| 国产区一区二久久| 97碰自拍视频| 99精品久久久久人妻精品| 欧美精品亚洲一区二区| 男女之事视频高清在线观看| 中文资源天堂在线| 午夜福利免费观看在线| 日日摸夜夜添夜夜添小说| 国产欧美日韩一区二区三| 精品欧美一区二区三区在线| 婷婷精品国产亚洲av在线| tocl精华| 熟女电影av网| 色播亚洲综合网| 欧美乱码精品一区二区三区| 久久久久久久午夜电影| 国产激情欧美一区二区| www.999成人在线观看| 日韩欧美在线二视频| 最好的美女福利视频网| 女警被强在线播放| 久久久久久九九精品二区国产 | 国产又黄又爽又无遮挡在线| 国产高清激情床上av| 国产亚洲欧美98| 久久久久久人人人人人| 午夜精品在线福利| 精品不卡国产一区二区三区| 母亲3免费完整高清在线观看| 国产伦一二天堂av在线观看| 极品教师在线免费播放| 久久精品国产综合久久久| 久久精品影院6| 成人18禁在线播放| 两个人免费观看高清视频| 狂野欧美激情性xxxx| 天天躁夜夜躁狠狠躁躁| 黄频高清免费视频| 国产伦一二天堂av在线观看| 99国产综合亚洲精品| 可以在线观看毛片的网站| 国产欧美日韩一区二区三| 成人三级做爰电影| 亚洲精品一区av在线观看| 可以在线观看的亚洲视频| 中文亚洲av片在线观看爽| 欧美成人性av电影在线观看| 少妇熟女aⅴ在线视频| 精品久久久久久久毛片微露脸| 看免费av毛片| 国产黄a三级三级三级人| 九色国产91popny在线| 久久婷婷人人爽人人干人人爱| 亚洲一区高清亚洲精品| 无遮挡黄片免费观看| 国产极品粉嫩免费观看在线| 91字幕亚洲| 久9热在线精品视频| 一区二区三区激情视频| 给我免费播放毛片高清在线观看| 亚洲人成网站在线播放欧美日韩| 日韩大码丰满熟妇| 亚洲人成网站在线播放欧美日韩| 十八禁人妻一区二区| 国产精品98久久久久久宅男小说| 日本 欧美在线| 又黄又爽又免费观看的视频| 国产成人精品久久二区二区91| 又黄又粗又硬又大视频| 香蕉丝袜av| 精品国产国语对白av| 久久久久久国产a免费观看| 欧美一级a爱片免费观看看 | 欧美日韩一级在线毛片| 日韩欧美 国产精品| 午夜精品在线福利| 悠悠久久av| 国产av又大| 动漫黄色视频在线观看| 久久午夜亚洲精品久久| 日韩欧美 国产精品| 亚洲色图 男人天堂 中文字幕| 熟女电影av网| 香蕉国产在线看| 国产视频一区二区在线看| 中亚洲国语对白在线视频| 欧美不卡视频在线免费观看 | 国产亚洲欧美精品永久| 国产主播在线观看一区二区| 国产片内射在线| 国产精品久久久人人做人人爽| √禁漫天堂资源中文www| 日韩精品免费视频一区二区三区| 美女扒开内裤让男人捅视频| 少妇被粗大的猛进出69影院| 中文在线观看免费www的网站 | 亚洲五月婷婷丁香| 在线观看www视频免费| 18禁裸乳无遮挡免费网站照片 | 69av精品久久久久久| 欧美日韩瑟瑟在线播放| 老鸭窝网址在线观看| 色老头精品视频在线观看| 日日干狠狠操夜夜爽| 天天一区二区日本电影三级| 啪啪无遮挡十八禁网站| 国产精品99久久99久久久不卡| 亚洲国产精品成人综合色| 色老头精品视频在线观看| 国内毛片毛片毛片毛片毛片| 亚洲最大成人中文| av超薄肉色丝袜交足视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产午夜精品久久久久久| 免费在线观看日本一区| 久久久久久久午夜电影| 91九色精品人成在线观看| 亚洲狠狠婷婷综合久久图片| 久久久久免费精品人妻一区二区 | 99精品在免费线老司机午夜| 午夜福利免费观看在线| 成人永久免费在线观看视频| 51午夜福利影视在线观看| 不卡av一区二区三区| 最近最新中文字幕大全免费视频| 亚洲精品一卡2卡三卡4卡5卡| 夜夜爽天天搞| 亚洲 欧美 日韩 在线 免费| 久久久国产成人免费| 丝袜在线中文字幕| 国产极品粉嫩免费观看在线| 最好的美女福利视频网| 人妻久久中文字幕网| 亚洲片人在线观看| xxx96com| 琪琪午夜伦伦电影理论片6080| 成人精品一区二区免费| 制服丝袜大香蕉在线| 欧美日韩瑟瑟在线播放| 久久精品成人免费网站| 丝袜在线中文字幕| 精品一区二区三区四区五区乱码| 国产在线观看jvid| 黄色丝袜av网址大全| 国产欧美日韩一区二区三| 日本一本二区三区精品| 欧美激情极品国产一区二区三区| 亚洲午夜理论影院| 两个人视频免费观看高清| 精品福利观看| 91国产中文字幕| 欧美日韩黄片免| 别揉我奶头~嗯~啊~动态视频| 18禁观看日本| 国产区一区二久久| 国内揄拍国产精品人妻在线 | 精品国产一区二区三区四区第35| 国产高清videossex| 国产一卡二卡三卡精品| 亚洲精品av麻豆狂野| 日韩视频一区二区在线观看| 欧美性长视频在线观看| 亚洲精品久久国产高清桃花| 国产成人欧美| av片东京热男人的天堂| 国内揄拍国产精品人妻在线 | 男人舔女人下体高潮全视频| 91大片在线观看| 老汉色av国产亚洲站长工具| 日日爽夜夜爽网站| 俄罗斯特黄特色一大片| 精品福利观看| 男人操女人黄网站| 2021天堂中文幕一二区在线观 | 亚洲国产欧洲综合997久久, | 国产高清有码在线观看视频 | 伊人久久大香线蕉亚洲五| 97超级碰碰碰精品色视频在线观看| 美女大奶头视频| 欧美日韩亚洲国产一区二区在线观看| 在线免费观看的www视频| 观看免费一级毛片| 亚洲av日韩精品久久久久久密| av电影中文网址| 啪啪无遮挡十八禁网站| 国产一区二区三区视频了| 无人区码免费观看不卡| 波多野结衣高清无吗| av福利片在线| 999精品在线视频| 性色av乱码一区二区三区2| 大香蕉久久成人网| 国产精品国产高清国产av| 黑丝袜美女国产一区| 中文字幕av电影在线播放| 18禁国产床啪视频网站| 成人18禁高潮啪啪吃奶动态图| 日韩 欧美 亚洲 中文字幕| 精品欧美一区二区三区在线| 国产精品久久久av美女十八| 99热6这里只有精品| 黄色 视频免费看| 久久九九热精品免费| 又紧又爽又黄一区二区| 久久精品91无色码中文字幕| 亚洲国产精品合色在线| 99国产精品一区二区蜜桃av| www日本在线高清视频| 欧美zozozo另类| 亚洲色图av天堂| 久久精品91无色码中文字幕| 亚洲专区国产一区二区| 中国美女看黄片| 国产视频一区二区在线看| 久久久久国内视频| 精品不卡国产一区二区三区| 在线观看午夜福利视频| 一边摸一边做爽爽视频免费| 免费看十八禁软件| www日本在线高清视频| 精品国内亚洲2022精品成人| 少妇被粗大的猛进出69影院| 高清在线国产一区| 国产91精品成人一区二区三区| 日韩欧美国产在线观看| 一进一出好大好爽视频| 最近最新中文字幕大全电影3 | 男女床上黄色一级片免费看| 国产久久久一区二区三区| 午夜福利欧美成人| 成熟少妇高潮喷水视频| 亚洲精品国产区一区二| 一本久久中文字幕| 熟女电影av网| 香蕉丝袜av| 99国产精品一区二区三区| 久久久久九九精品影院| 韩国精品一区二区三区| 久久伊人香网站| 久久久水蜜桃国产精品网| 亚洲专区字幕在线| 午夜福利视频1000在线观看| 少妇粗大呻吟视频| 午夜福利成人在线免费观看| 99精品欧美一区二区三区四区| 热re99久久国产66热| 69av精品久久久久久| 色哟哟哟哟哟哟| 欧美中文日本在线观看视频| 久久久久久久久久黄片| 成人永久免费在线观看视频| 在线观看免费日韩欧美大片| 久久人人精品亚洲av| 午夜久久久久精精品| or卡值多少钱| 亚洲精品一卡2卡三卡4卡5卡| av福利片在线| 亚洲国产看品久久| 国产精品电影一区二区三区| 99国产精品一区二区蜜桃av| 午夜视频精品福利| 熟女少妇亚洲综合色aaa.| 国产成年人精品一区二区| 国产成+人综合+亚洲专区| 亚洲精品久久国产高清桃花| av天堂在线播放| 国产高清videossex| a级毛片a级免费在线| 一二三四在线观看免费中文在| bbb黄色大片| 日韩欧美一区视频在线观看| 高潮久久久久久久久久久不卡| 91av网站免费观看| 亚洲中文字幕一区二区三区有码在线看 | 变态另类丝袜制服| 亚洲国产精品合色在线| 欧美日韩福利视频一区二区| 久久精品夜夜夜夜夜久久蜜豆 | 日本熟妇午夜| 亚洲一卡2卡3卡4卡5卡精品中文| 美女午夜性视频免费| 亚洲 欧美一区二区三区| 大型av网站在线播放| 久久久久久人人人人人| 国产成人啪精品午夜网站| 国产精品九九99| 欧美不卡视频在线免费观看 | 久久久久久久久中文| 亚洲国产欧美一区二区综合| 男人舔奶头视频| 天堂影院成人在线观看| 亚洲性夜色夜夜综合| 久久国产精品影院| 成年免费大片在线观看| 国产精品二区激情视频| 不卡av一区二区三区| 国产熟女午夜一区二区三区| av欧美777| 国产成人精品久久二区二区免费| 这个男人来自地球电影免费观看| 国产野战对白在线观看| 久久热在线av| 国产伦人伦偷精品视频| 黄色 视频免费看| 亚洲男人的天堂狠狠| 久久久久久大精品| 亚洲熟女毛片儿| 少妇的丰满在线观看| 亚洲熟妇熟女久久| 国产精品1区2区在线观看.| 很黄的视频免费| 精品乱码久久久久久99久播| 亚洲va日本ⅴa欧美va伊人久久| 别揉我奶头~嗯~啊~动态视频| 欧美黑人精品巨大| 欧美日本视频| tocl精华| 欧美黑人精品巨大| 中文资源天堂在线| 香蕉国产在线看| 男人舔女人下体高潮全视频| 国产成人精品久久二区二区免费| 久久久久九九精品影院| 18禁黄网站禁片免费观看直播| 丁香六月欧美| 日韩大码丰满熟妇| 美女免费视频网站| 精品久久久久久久人妻蜜臀av| 丝袜人妻中文字幕| 欧美av亚洲av综合av国产av| 99精品在免费线老司机午夜| 久久天躁狠狠躁夜夜2o2o| 免费看a级黄色片| 高潮久久久久久久久久久不卡| 啪啪无遮挡十八禁网站| 免费看a级黄色片| 啦啦啦观看免费观看视频高清| 欧美在线一区亚洲| 国产av在哪里看| 久久久久久久久免费视频了| 亚洲国产精品成人综合色| 国产av在哪里看| 一级毛片女人18水好多| 久久天躁狠狠躁夜夜2o2o| 精品久久久久久,| 久久 成人 亚洲| 大香蕉久久成人网| 男女下面进入的视频免费午夜 | 国产一区二区在线av高清观看| 在线观看66精品国产| 成年女人毛片免费观看观看9| 国产精品日韩av在线免费观看| 国产亚洲精品久久久久久毛片| 婷婷丁香在线五月| 欧美日韩中文字幕国产精品一区二区三区| 亚洲中文字幕日韩|