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

    堤防CT成像的數(shù)值模擬研究

    2021-03-28 11:49:17潘紀順劉宇鋒李長征
    人民黃河 2021年3期
    關(guān)鍵詞:堤防梯形射線

    潘紀順 劉宇鋒 李長征

    摘 要:簡述以射線理論為基礎(chǔ)的地震CT的基本原理、數(shù)學方法、圖像重構(gòu)技術(shù)、ART算法和SIRT算法,在現(xiàn)有井間地震CT的基礎(chǔ)上研究堤防CT,利用MATLAB編程對堤防進行地震斷面CT的數(shù)值模擬,通過ART和SIRT重建方法進行反演,分析得出二者在堤防CT應(yīng)用中的優(yōu)缺點及適用條件。結(jié)果表明:ART算法和SIRT算法都得出了較好的反演結(jié)果,提高了探測效率和精度,較好地驗證了堤防CT成像技術(shù)的正確性和可行性。

    關(guān)鍵詞:堤防;地震CT;ART算法;SIRT算法

    Abstract: The basic principle, mathematical method, image reconstruction technology, ART algorithm and SIRT algorithm of seismic CT based on ray theory were briefly introduced. Based on the existing cross well CT, the embankment CT was studied and the seismic cross-section CT numerical simulation was carried out by MATLAB programming. The advantages and disadvantages of ART and SIRT reconstruction methods in the application of embankment CT were analyzed. Finally, good inversion results were obtained, which improved the accuracy and verified the correctness and feasibility of embankment CT imaging technology.

    Key words: embankment; seismic CT; ART algorithm; SIRT algorithm

    1 引 言

    長久以來,每到汛期,洪災(zāi)威脅對兩岸居民的安全和生活造成了嚴重的影響。抵御洪水泛濫、保護人民生活和工農(nóng)業(yè)的生產(chǎn),堤防已經(jīng)成為一道重要的屏障。我國江河堤防質(zhì)量良莠不齊,為了有效防控洪水,堤防隱患探測工作已經(jīng)變得十分重要。如何高效、準確地確定隱患的分布位置和大小,已經(jīng)成為堤防隱患探測迫切需要解決的問題 [1]。

    20世紀80年代以來,我國開始采用物探方法[2]開展堤防隱患探測的研究工作,如高密度電法[3]、地質(zhì)雷達法[4]、瞬變電磁法[5]等。物探方法是根據(jù)堤防土體與隱患之間的地球物理特性差異來判斷是否存在堤防隱患,在堤防質(zhì)量檢測中起到了重要的作用。

    目前,地震CT技術(shù)[6]已經(jīng)被多次應(yīng)用于地球物理勘探的不同工況中[7]。閻生存等[8]對大壩進行ART算法成像研究,推動了CT技術(shù)在大壩隱患檢測中的發(fā)展[9];胡明順等[10]提出井地地震CT探測的方法;李世民等[11]運用地震CT技術(shù)對隧道前方地質(zhì)情況進行超前預(yù)報;靳洪曉等[12]運用地面地震CT技術(shù)對淺層地質(zhì)情況進行勘探研究。但是目前地震CT技術(shù)應(yīng)用在堤防上的案例較少,堤防CT和一般的井間地震CT不同,堤防是一個近似規(guī)則的梯形結(jié)構(gòu)[13],兩邊均有一定坡度,分為迎水面和背水面。因此,堤防CT模型應(yīng)該是一個梯形模型。

    堤防模型的建立,除了考慮近似于梯形結(jié)構(gòu)的地形影響外,還需要考慮堤防的構(gòu)筑材料對波速的影響。堤防主要有土堤、鋼筋混凝土堤、土石混合堤等,其中土堤最為常見,在某些缺乏土料的地區(qū)常采用土石混合堤。堤防的上層一般是不透水覆蓋層,下層一般為透水的砂礫石層和細砂層。堤防的土一般為壤土類土和黏土類土[14],其中可能夾雜著少量的砂土。這些堤防中的不同構(gòu)筑材料都對波速有著不同程度的影響。

    堤防的迎水面和背水面工作條件較好,在采用地震CT進行堤防隱患探測時,容易實現(xiàn)地震波的激發(fā)和接收。本文利用基于MATLAB編寫的程序進行堤防CT的數(shù)值模擬,反演堤防中隱患的位置和大小,從而驗證在堤防上進行地震CT的準確性和可行性。

    2 堤防CT基本原理

    2.1 離散圖像重建原理

    如圖1所示,在Y方向沿梯形左側(cè)(迎水面)安置系列震源、Y方向沿梯形右側(cè)(背水面)安置系列檢波器。

    震源產(chǎn)生的地震波被檢波器采集后,以電信號的形式傳輸至主機。通過拾取地震波初至時間,可實現(xiàn)反演這個梯形剖面的速度v(x,y)或慢度f(x,y)=1/v(x,y)。設(shè)第i個地震波的走時為τi,傳播路徑的長度為Li,則

    式(1)中df代表弧長微元,速度和傳播路徑都是未知條件,而時間τi是試驗采集到的數(shù)據(jù),為已知的。在堤防中地下介質(zhì)近似相同,速度場可以看作是不變的,射線的路徑Li看作近似直線。但是,在實際的堤防中地下介質(zhì)情況是不同的,射線路徑也往往是彎曲的。假設(shè)圖像重建區(qū)的網(wǎng)格數(shù)目為M,單元格的慢度依次為f1、f2、…、fM。第i條射線的傳播時間為

    通過式(3)計算得出一個速度場的分布矩陣,根據(jù)獲得的數(shù)據(jù)進行反演成像。在計算反演矩陣A時往往是沒有規(guī)則的,且其中的射線只通過了反演區(qū)域的一部分。在實際工程應(yīng)用中常用迭代法求解式(3),重建探測區(qū)域的速度場。

    2.2 ART算法

    目前,反演的方法很多,主要有以下幾種:矩陣法、付氏變換法、卷積反投影法、代數(shù)重建法(簡稱ART法)[15]。ART算法[16]適用性較廣,直射線和彎曲射線均能使用。

    ART算法的基本思想是先給反演區(qū)域中的每個像元都賦予一個初始時間值tij,一般為0;然后將投影結(jié)果沿射線方向反投影回原像元,同時不斷地對反演區(qū)域進行修正;最終在滿足要求時停止修正。假設(shè)fq,ij為第q輪迭代時第i條射線上對于第j像元的慢度(估計值),則第i條射線的理論走時可以表示為

    所有射線修正完后,若計算得出走時差的模ΔTq<ε(ε為給定允許誤差),則停止迭代;否則,再進行第q+1輪的迭代,從i=1開始。

    2.3 SIRT算法

    SIRT算法[17]是由Gilbert首先提出的。該算法與ART算法不同,典型的ART算法用第j條射線的修改值來計算第j+1條射線的修改值,而SIRT算法是把第k輪迭代的所有射線的修改值保存到一個指定空間,在這一輪結(jié)束后求得某種平均修改值Δf-,然后再對每個像元的慢度進行修正,即f(k+1)=f(k)+Δf-,并用于下一輪迭代。

    3 軟件介紹

    軟件基于MATLAB編寫,反演方法使用ART算法和SIRT算法。MATLAB軟件可以用于矩陣運算、函數(shù)和數(shù)據(jù)繪制、三維成圖和創(chuàng)建圖形用戶界面等[18],將數(shù)據(jù)進行矩陣化處理,運用MATLAB進行計算,極大地提高了工作效率。

    軟件主要分為兩個部分,分別為規(guī)則的梯形堤防模型和有淤背區(qū)的堤防模型。有淤背區(qū)的堤防模型是在規(guī)則梯形軟件的基礎(chǔ)上設(shè)計一個淤背區(qū)。

    梯形模型軟件的主要功能是建立一個堤防梯形模型,對此模型進行正演和反演,得出地震波在堤防中的速度分布圖,最后找出堤防隱患的位置,從而提高堤防隱患探測的效率和數(shù)據(jù)解釋的可視化程度。軟件的運行過程如下:第1步,建立一個網(wǎng)格化的梯形模型,在此模型上設(shè)置激發(fā)點和采集點;第2步,分別設(shè)置實際的背景探測速度和假設(shè)的背景探測速度,通過計算得到每條射線經(jīng)過像元的長度,并計算出每條射線的理論時間和實際時間;第3步,將結(jié)果利用ART算法和SIRT算法進行反演,當滿足條件時,停止迭代;第4步,計算結(jié)果成圖,得出探測區(qū)速度分布圖,從圖中的速度分布情況可以看出隱患異常體所在的位置。

    4 數(shù)值模擬

    堤防觀測系統(tǒng)具備對稱性,因此建立一個梯形堤防數(shù)值模型。異常區(qū)與堤防填筑料的平均速度差對成像質(zhì)量的影響很大,速度差越大,成像質(zhì)量越好。地震波通過堤防時,如果遇裂縫、空洞等異常體,速度會降低。假設(shè)規(guī)則梯形探測區(qū)背景速度為1 000 m/s,模型中異常體像元的速度為900 m/s,如圖2(a)所示,設(shè)計一個規(guī)則的二維梯形模型,網(wǎng)格區(qū)域為5×6,網(wǎng)格劃分為1×1。梯形左側(cè)為迎水面,右側(cè)為背水面。激發(fā)點分布在梯形的左側(cè),采集點分布在梯形的右側(cè),設(shè)置激發(fā)點和采集點間距為0.5 m或1 m,重建區(qū)網(wǎng)格數(shù)為26。

    淤背區(qū)對堤防工程具有重要的作用,能夠增強堤防的抗?jié)B和抗震穩(wěn)定性。因此,設(shè)計重建區(qū)網(wǎng)格圖形模型為有淤背區(qū)的堤防,右側(cè)有一個高2 m、寬1 m的淤背區(qū)。基于MATLAB使用ART算法和SIRT算法進行反演,最大迭代次數(shù)均為20次。

    采用圖像距離測量值來評價圖像重建效果,即歸一化均方根距離測量值d、歸一化平均絕對距離測量值r和最壞情況距離測量值e。d=0表示重建圖像忠實地再現(xiàn)原始圖像,d值愈大表示兩者偏差愈大;r=0說明沒有誤差,r增大,說明誤差增大;e為原始圖像和重建圖像之間的最大平均慢度差。建立3類模型,第1類為含低速體堤防,第2類為層狀堤防,第3類為有淤背區(qū)的層狀堤防,進行對比分析。

    4.1 含低速體堤防模型

    圖2(a)為激發(fā)點和采集點的間距為1 m的含低速體模型,一共25條射線。圖3(a)為激發(fā)點和采集點的間距為0.5 m的含低速體模型,一共81條射線。

    圖2和圖3中的迭代次數(shù)均設(shè)置為20次。激發(fā)點與采集點的間距越小,經(jīng)過每一個像元的射線條數(shù)越多。由表1和圖3可知,在迭代10次后,SIRT算法結(jié)果趨于穩(wěn)定,ART算法的圖像距離測量值還在增大;迭代20次后ART算法的結(jié)果才趨于穩(wěn)定。由此可見,對于堤防含低速體模型,SIRT算法比ART算法收斂速度更快,且收斂性更好。

    4.2 層狀堤防模型

    層狀堤防ART算法和SIRT反演結(jié)果見圖4。圖4(a)為激發(fā)點和采集點的間距為0.5 m的層狀堤防模型,一共81條射線。探測區(qū)的背景速度為1 000 m/s,模型中像元(層狀)的速度為900 m/s。

    由圖4和表2可知,SIRT算法的結(jié)果好于ART算法,且收斂速度更快。ART算法在第10次迭代后,d值還在增大,而SIRT算法d值在逐漸減小。圖4中(c)和(d)迭代次數(shù)均為20次,都可以明顯看出異常體的位置和大小。

    4.3 有淤背區(qū)的層狀堤防模型

    有淤背的層狀堤防模型反演結(jié)果見圖5和表3。圖5(a)為激發(fā)點和采集點的間距為0.5 m的有淤背區(qū)的層狀堤防模型,共78條射線。探測區(qū)背景速度為1 000 m/s,異常體的速度為900 m/s,在淤背區(qū)有輕微繞射,速度為800 m/s。

    由圖5和表3可知,SIRT算法的收斂速度大于ART算法的,且受到淤背區(qū)的干擾較小,上下層都能明顯反演出異常體的位置和大小。

    5 結(jié) 語

    (1)由于地震波傳播實際路徑并不是一條直線,因此不能完全重建出真實速度場,但根據(jù)速度場相鄰區(qū)域的差異可以識別出異常體在堤防中的空間分布及形態(tài)。

    (2)當?shù)螖?shù)相同時,SIRT算法重建圖像的距離測量值比ART算法的小。當距離測量值相近時,SIRT算法重建的迭代次數(shù)少于ART算法的??梢奡IRT算法比ART算法的迭代收斂性要好,收斂速度快,而且對投影數(shù)據(jù)誤差的敏感度小。

    (3)震源和檢波器的排布設(shè)計直接影響射線分布的密度、均勻性,從而影響圖像重建的精度和清晰度。因此,在實際應(yīng)用中,要選擇合適的震源和檢波器的排布,以獲得真實可靠的重建圖像。

    利用MATLAB編程軟件進行堤防CT的隱患探測研究,通過數(shù)值模擬,ART算法和SIRT算法反演都得出了較好的結(jié)果,驗證了地震CT進行堤防隱患探測的正確性和可行性。

    將地震CT方法應(yīng)用在堤防隱患探測上,提高了隱患探測的效率和精度,并且減少了鉆孔費用。在實際工作中,如果在觀測系統(tǒng)設(shè)計、數(shù)據(jù)采集和數(shù)據(jù)處理三個階段選擇合適的方法,那么堤防CT成像技術(shù)將取得較好的探測結(jié)果。下一步還需要進行實際堤防隱患探測試驗,與實際鉆孔資料對比,進一步驗證堤防CT成像技術(shù)的實用性。

    參考文獻:

    [1] 周華敏,肖國強,周黎明,等.堤防隱患物探技術(shù)研究現(xiàn)狀與展望[J].長江科學院院報,2019,36(12):164-168.

    [2] 欒明龍.堤防隱患的地球物理模型及其探測與監(jiān)測的思考[J].北京水利,2005(5):36-38.

    [3] 宋朝陽,王銳,李長征,等.高密度電法探測堤防隱患研究[J].人民黃河,2020,42(7):104-106.

    [4] 武桂芝,張寶森,李春江,等.陣列地質(zhì)雷達在黃河堤防隱患探測中的應(yīng)用[J].人民黃河,2020,42(8):113-116.

    [5] 鄧洪亮,謝向文,郭玉松,等.黃河下游堤防工程隱患探測技術(shù)與應(yīng)用[J].地球物理學進展,2008,23(3):936-941.

    [6] 王旭,張晨,由靜,等.地震波CT在隧道地質(zhì)災(zāi)害探測中的應(yīng)用[J].大壩與安全,2010(1):49-52.

    [7] 羅炬,李志海,王海濤.地震層析成像研究方法綜述[J].內(nèi)陸地震,2011,25(4):311-320.

    [8] 閻生存,李珍照,薛桂玉,等.大壩CT中ART算法及其改進探討[J].武漢大學學報(工學版),2001,34(4):29-34.

    [9] 余志雄,薛桂玉,周洪波,等.大壩CT技術(shù)研究概況與進展[J].巖石力學與工程學報,2004,23(8):1394-1397.

    [10] 胡明順,潘冬明,李兵,等.井地地震CT成像數(shù)值模擬研究與應(yīng)用[J].地球物理學進展,2009,24(4):1399-1404.

    [11] 李世民,張曉培,牛建軍,等.地震CT技術(shù)在隧道工程中的應(yīng)用[J].探礦工程(巖土鉆掘工程),2004(9):63-65.

    [12] 靳洪曉,趙永貴,李勤,等.地面地震CT在淺層勘探中的應(yīng)用[J].工程地質(zhì)學報,2000,8(2):239-243.

    [13] 孫東亞,姚秋玲,趙進勇,等.堤防工程建設(shè)技術(shù)進展[J].中國防汛抗旱,2009,19(6):34-37.

    [14] 楊代齊.水利堤防土方填筑方法探討[J].水利技術(shù)監(jiān)督,2015,23(3):92-95.

    [15] 楊文采.地球物理反演的理論與方法[M].北京:地質(zhì)出版社,1997:83-100.

    [16] 駱循,朱介壽.ART算法及其在地球物理反演中的應(yīng)用[J].CT理論與應(yīng)用研究,1987(1):8-16.

    [17] 楊利普,徐志萍,李德慶,等.SIRT法層析成像在某核電站巖溶區(qū)的應(yīng)用[J].CT理論與應(yīng)用研究,2014(2):6-14.

    [18] 段成龍,閻長虹,許寶田,等.跨孔地震CT技術(shù)在地鐵工程施工溶洞探測方面的應(yīng)用[J].地質(zhì)論評,2013,59(6):1242-1248.

    【責任編輯 許立新】

    猜你喜歡
    堤防梯形射線
    玩轉(zhuǎn)梯形
    幼兒100(2023年37期)2023-10-23 11:38:58
    “直線、射線、線段”檢測題
    梯形達人
    一類變延遲中立型微分方程梯形方法的漸近估計
    『直線、射線、線段』檢測題
    梯形
    啟蒙(3-7歲)(2017年6期)2017-11-27 09:34:55
    2016年河南省己建成堤防長度
    赤石脂X-射線衍射指紋圖譜
    中成藥(2017年3期)2017-05-17 06:09:16
    河南省2014年已治理堤防長度
    河南省2014年已治理達標堤防長度
    你懂的网址亚洲精品在线观看| 亚洲激情五月婷婷啪啪| 国产伦精品一区二区三区视频9| 亚洲av电影在线观看一区二区三区| 国产精品熟女久久久久浪| 中文字幕制服av| 亚洲性久久影院| 高清黄色对白视频在线免费看 | 免费看av在线观看网站| 伊人久久国产一区二区| 成人美女网站在线观看视频| 七月丁香在线播放| 国产一区二区在线观看av| 极品人妻少妇av视频| 日本色播在线视频| 欧美日韩在线观看h| 免费人妻精品一区二区三区视频| 男人舔奶头视频| 国产精品偷伦视频观看了| 国产精品福利在线免费观看| 美女主播在线视频| 国产精品国产三级国产av玫瑰| 国产亚洲午夜精品一区二区久久| 9色porny在线观看| av卡一久久| 亚洲国产色片| 女的被弄到高潮叫床怎么办| 99热这里只有是精品50| 日韩一区二区视频免费看| 免费观看在线日韩| 免费久久久久久久精品成人欧美视频 | 日本91视频免费播放| h日本视频在线播放| 欧美日韩在线观看h| 人人澡人人妻人| 国产精品国产三级国产专区5o| 国产av一区二区精品久久| 成人毛片60女人毛片免费| 一本—道久久a久久精品蜜桃钙片| 精品卡一卡二卡四卡免费| 亚洲人成网站在线播| 一个人看视频在线观看www免费| 欧美精品一区二区免费开放| 久久6这里有精品| 妹子高潮喷水视频| 蜜桃在线观看..| 国产极品天堂在线| 偷拍熟女少妇极品色| 秋霞在线观看毛片| 国产成人freesex在线| 国产真实伦视频高清在线观看| 青春草亚洲视频在线观看| 久久精品夜色国产| 熟女av电影| av在线观看视频网站免费| 日韩强制内射视频| 国产精品伦人一区二区| .国产精品久久| 日本wwww免费看| 男男h啪啪无遮挡| 成年av动漫网址| 成人亚洲精品一区在线观看| 色吧在线观看| 精品人妻熟女av久视频| 中文资源天堂在线| 99热国产这里只有精品6| 啦啦啦视频在线资源免费观看| 欧美成人精品欧美一级黄| 欧美激情国产日韩精品一区| 如日韩欧美国产精品一区二区三区 | 22中文网久久字幕| 国产精品欧美亚洲77777| 精品人妻熟女av久视频| 内射极品少妇av片p| 如日韩欧美国产精品一区二区三区 | 亚洲精品自拍成人| 国产精品99久久久久久久久| 秋霞伦理黄片| 亚洲成色77777| 亚洲一区二区三区欧美精品| 免费观看a级毛片全部| 九色成人免费人妻av| 久久97久久精品| 一级黄片播放器| 午夜免费观看性视频| 大话2 男鬼变身卡| 一本—道久久a久久精品蜜桃钙片| av在线播放精品| 日韩一区二区三区影片| 精品亚洲乱码少妇综合久久| 伦精品一区二区三区| 在线观看av片永久免费下载| 2018国产大陆天天弄谢| 久久 成人 亚洲| 性色av一级| 久久久精品免费免费高清| 在线看a的网站| 少妇人妻 视频| 亚洲av电影在线观看一区二区三区| 在线 av 中文字幕| 三上悠亚av全集在线观看 | 中文精品一卡2卡3卡4更新| 一级毛片我不卡| 国产欧美日韩一区二区三区在线 | www.av在线官网国产| 日产精品乱码卡一卡2卡三| 九色成人免费人妻av| 一区二区三区免费毛片| 99热6这里只有精品| av在线app专区| 女性生殖器流出的白浆| 婷婷色综合大香蕉| 青青草视频在线视频观看| 人妻系列 视频| 赤兔流量卡办理| 午夜福利网站1000一区二区三区| av免费在线看不卡| av国产久精品久网站免费入址| 狠狠精品人妻久久久久久综合| 色5月婷婷丁香| 亚洲av国产av综合av卡| 草草在线视频免费看| 日本av手机在线免费观看| 亚洲精品中文字幕在线视频 | 极品人妻少妇av视频| 日日撸夜夜添| 99九九在线精品视频 | 美女国产视频在线观看| 好男人视频免费观看在线| 成年av动漫网址| 五月开心婷婷网| 中文欧美无线码| 亚洲欧洲日产国产| 精品人妻熟女av久视频| 日韩欧美一区视频在线观看 | 精品卡一卡二卡四卡免费| 免费黄网站久久成人精品| 免费看光身美女| 久久热精品热| 少妇的逼好多水| 老司机亚洲免费影院| 国产精品麻豆人妻色哟哟久久| 波野结衣二区三区在线| 最新中文字幕久久久久| 国产av国产精品国产| 精品人妻熟女av久视频| 又爽又黄a免费视频| 伦精品一区二区三区| 国产精品福利在线免费观看| 成人影院久久| 热99国产精品久久久久久7| 亚洲精品乱码久久久v下载方式| 国产精品99久久久久久久久| 男女边摸边吃奶| 亚洲精品国产色婷婷电影| 少妇人妻精品综合一区二区| 国产色爽女视频免费观看| 日本黄色日本黄色录像| 日日啪夜夜撸| 中文精品一卡2卡3卡4更新| 午夜91福利影院| 人人妻人人添人人爽欧美一区卜| 日日摸夜夜添夜夜爱| 久久青草综合色| 国产 精品1| 老熟女久久久| 18+在线观看网站| 久久久久久久久久久久大奶| 大码成人一级视频| 国模一区二区三区四区视频| 久久精品久久久久久久性| 久久久国产欧美日韩av| 免费人成在线观看视频色| 精品少妇黑人巨大在线播放| 亚洲av二区三区四区| 亚洲av欧美aⅴ国产| 极品人妻少妇av视频| 午夜福利视频精品| 91精品国产九色| 国产毛片在线视频| 亚洲精品亚洲一区二区| 一边亲一边摸免费视频| 热99国产精品久久久久久7| 99热这里只有是精品50| 制服丝袜香蕉在线| 亚洲国产欧美在线一区| 国语对白做爰xxxⅹ性视频网站| 国产精品久久久久久久久免| 在线观看www视频免费| 久久久精品免费免费高清| 五月伊人婷婷丁香| 最近2019中文字幕mv第一页| 国产在线免费精品| 自拍欧美九色日韩亚洲蝌蚪91 | 国产精品成人在线| 日本欧美国产在线视频| 午夜福利,免费看| 欧美日韩在线观看h| videossex国产| 久久毛片免费看一区二区三区| 亚洲天堂av无毛| 夜夜骑夜夜射夜夜干| 免费看不卡的av| 亚洲怡红院男人天堂| 狂野欧美激情性bbbbbb| 久久人妻熟女aⅴ| 国产精品欧美亚洲77777| 欧美 日韩 精品 国产| 女人久久www免费人成看片| av免费在线看不卡| 国产欧美日韩一区二区三区在线 | 国国产精品蜜臀av免费| 丁香六月天网| 国产午夜精品久久久久久一区二区三区| 精品国产国语对白av| 欧美变态另类bdsm刘玥| 各种免费的搞黄视频| 99热全是精品| 高清黄色对白视频在线免费看 | 国产成人精品久久久久久| 久久久久久久精品精品| 亚洲人成网站在线观看播放| 亚洲综合色惰| 欧美精品亚洲一区二区| 亚洲第一区二区三区不卡| 综合色丁香网| 久热这里只有精品99| 亚洲怡红院男人天堂| 最后的刺客免费高清国语| 久久精品国产a三级三级三级| 久久久国产精品麻豆| 99热国产这里只有精品6| 成人漫画全彩无遮挡| 色视频www国产| 日韩精品免费视频一区二区三区 | 在线观看三级黄色| 婷婷色综合www| 蜜桃久久精品国产亚洲av| 三级经典国产精品| 亚洲激情五月婷婷啪啪| 高清av免费在线| 国内精品宾馆在线| 少妇的逼水好多| 国产永久视频网站| 大码成人一级视频| 亚洲精品国产av成人精品| videos熟女内射| 欧美区成人在线视频| 十分钟在线观看高清视频www | 国产精品麻豆人妻色哟哟久久| 日本黄大片高清| 99久久精品热视频| 一级毛片我不卡| av在线老鸭窝| 黑丝袜美女国产一区| 日本vs欧美在线观看视频 | 男女无遮挡免费网站观看| av网站免费在线观看视频| 多毛熟女@视频| 97精品久久久久久久久久精品| 99久久精品一区二区三区| 国产女主播在线喷水免费视频网站| 97在线视频观看| a级毛色黄片| 日本黄色日本黄色录像| 精品国产一区二区久久| 性色avwww在线观看| av在线老鸭窝| 天堂中文最新版在线下载| 国产亚洲91精品色在线| 狂野欧美激情性bbbbbb| 亚洲久久久国产精品| 国内揄拍国产精品人妻在线| 建设人人有责人人尽责人人享有的| 制服丝袜香蕉在线| 性色avwww在线观看| 亚洲在久久综合| 欧美精品国产亚洲| 久久久国产欧美日韩av| 精品一区二区免费观看| 中文字幕免费在线视频6| 日韩av不卡免费在线播放| 亚洲精品国产av成人精品| 欧美日本中文国产一区发布| 国产精品女同一区二区软件| 国产精品人妻久久久久久| 亚洲精品乱久久久久久| 国产91av在线免费观看| 欧美日韩精品成人综合77777| 亚洲欧美一区二区三区国产| 熟女电影av网| 精品卡一卡二卡四卡免费| 麻豆成人av视频| 中文精品一卡2卡3卡4更新| 亚洲不卡免费看| av播播在线观看一区| 亚洲国产色片| 欧美日韩一区二区视频在线观看视频在线| 亚洲高清免费不卡视频| 久久久久久久久久久免费av| 国产av国产精品国产| 中文资源天堂在线| 久久影院123| 亚洲三级黄色毛片| 777米奇影视久久| 极品少妇高潮喷水抽搐| 日韩不卡一区二区三区视频在线| 中文字幕制服av| 久久久亚洲精品成人影院| 热99国产精品久久久久久7| 性高湖久久久久久久久免费观看| 天美传媒精品一区二区| 一个人看视频在线观看www免费| 亚洲欧洲精品一区二区精品久久久 | 久久午夜福利片| 免费大片黄手机在线观看| 人人澡人人妻人| 免费看不卡的av| 精品久久久久久久久av| 一级毛片aaaaaa免费看小| 丰满人妻一区二区三区视频av| 黑人巨大精品欧美一区二区蜜桃 | 亚洲一区二区三区欧美精品| 97精品久久久久久久久久精品| 国产免费一级a男人的天堂| 夜夜骑夜夜射夜夜干| av福利片在线观看| 一二三四中文在线观看免费高清| 国产精品.久久久| av国产精品久久久久影院| 午夜福利,免费看| 亚洲欧洲精品一区二区精品久久久 | 97在线视频观看| 亚洲三级黄色毛片| 日本猛色少妇xxxxx猛交久久| 久久久a久久爽久久v久久| 菩萨蛮人人尽说江南好唐韦庄| 国产91av在线免费观看| 国产精品秋霞免费鲁丝片| 晚上一个人看的免费电影| 如日韩欧美国产精品一区二区三区 | 在现免费观看毛片| 国产成人精品无人区| 夜夜骑夜夜射夜夜干| 人体艺术视频欧美日本| 一本大道久久a久久精品| 最近中文字幕2019免费版| 亚洲高清免费不卡视频| 美女脱内裤让男人舔精品视频| 久久午夜福利片| 人人妻人人澡人人爽人人夜夜| 91精品国产国语对白视频| 国产在线免费精品| 美女内射精品一级片tv| 熟女电影av网| 老司机影院成人| 制服丝袜香蕉在线| 人妻少妇偷人精品九色| 婷婷色麻豆天堂久久| 午夜激情福利司机影院| 亚洲精品,欧美精品| 99热全是精品| 欧美日韩在线观看h| 免费看av在线观看网站| 少妇被粗大猛烈的视频| 99久久人妻综合| 99久久精品热视频| 免费黄频网站在线观看国产| 精品午夜福利在线看| 丝袜喷水一区| 国产精品人妻久久久影院| 一级爰片在线观看| 午夜久久久在线观看| 精品一品国产午夜福利视频| 欧美另类一区| 麻豆成人av视频| 色视频在线一区二区三区| 99国产精品免费福利视频| 久久久久久久精品精品| 三级经典国产精品| 最近手机中文字幕大全| 亚洲av日韩在线播放| 麻豆成人午夜福利视频| 婷婷色麻豆天堂久久| 日韩一区二区三区影片| 免费观看a级毛片全部| 亚洲国产成人一精品久久久| 精品久久久久久电影网| 国产日韩欧美在线精品| 黑人猛操日本美女一级片| 亚洲人成网站在线观看播放| 国产综合精华液| 99国产精品免费福利视频| 日本爱情动作片www.在线观看| 国产欧美亚洲国产| 久久午夜福利片| 一级毛片黄色毛片免费观看视频| 亚洲国产精品成人久久小说| 日日啪夜夜撸| 女的被弄到高潮叫床怎么办| 国内揄拍国产精品人妻在线| 久久久国产欧美日韩av| 有码 亚洲区| 久久久欧美国产精品| 欧美性感艳星| 性高湖久久久久久久久免费观看| 午夜福利影视在线免费观看| 两个人的视频大全免费| 国产精品国产三级国产专区5o| 亚洲美女黄色视频免费看| 如何舔出高潮| 日日撸夜夜添| av专区在线播放| 亚洲精品国产色婷婷电影| 国产毛片在线视频| 男男h啪啪无遮挡| 熟女电影av网| 熟女人妻精品中文字幕| 欧美日韩在线观看h| 亚洲国产成人一精品久久久| 女人精品久久久久毛片| 成年女人在线观看亚洲视频| 我要看日韩黄色一级片| 亚洲中文av在线| 精品久久久精品久久久| 纵有疾风起免费观看全集完整版| 99热这里只有是精品在线观看| 亚洲天堂av无毛| 老熟女久久久| 一区二区av电影网| 男人狂女人下面高潮的视频| 两个人的视频大全免费| 欧美日韩一区二区视频在线观看视频在线| 国产男人的电影天堂91| 久久精品国产亚洲av天美| 亚洲精品aⅴ在线观看| 我的女老师完整版在线观看| 亚洲真实伦在线观看| 一本一本综合久久| 免费av不卡在线播放| 日本av手机在线免费观看| 一级毛片aaaaaa免费看小| 亚洲综合色惰| 男人狂女人下面高潮的视频| 精品人妻偷拍中文字幕| 日本欧美国产在线视频| 亚洲成人av在线免费| 一个人看视频在线观看www免费| 日韩一区二区三区影片| 一级毛片黄色毛片免费观看视频| 男男h啪啪无遮挡| 国产在线免费精品| 校园人妻丝袜中文字幕| 国产毛片在线视频| 成人18禁高潮啪啪吃奶动态图 | 久久精品国产亚洲网站| 亚洲精品亚洲一区二区| 国产成人aa在线观看| 欧美精品人与动牲交sv欧美| 国内精品宾馆在线| 肉色欧美久久久久久久蜜桃| 欧美日韩国产mv在线观看视频| 三级经典国产精品| 黄色毛片三级朝国网站 | 国产黄片视频在线免费观看| av视频免费观看在线观看| 九色成人免费人妻av| 免费少妇av软件| 一本—道久久a久久精品蜜桃钙片| 久久久久久久久久久丰满| 日韩精品有码人妻一区| 亚洲精品国产av蜜桃| 观看av在线不卡| 日韩不卡一区二区三区视频在线| 亚洲欧洲日产国产| 久久久精品免费免费高清| 99热网站在线观看| 欧美精品亚洲一区二区| 插逼视频在线观看| 午夜免费鲁丝| a级一级毛片免费在线观看| 最近最新中文字幕免费大全7| 妹子高潮喷水视频| 男女啪啪激烈高潮av片| 一区二区三区免费毛片| 99国产精品免费福利视频| 国产精品免费大片| 久久久a久久爽久久v久久| 哪个播放器可以免费观看大片| 国产片特级美女逼逼视频| 日韩熟女老妇一区二区性免费视频| 97在线视频观看| 国产一区有黄有色的免费视频| 久久精品国产自在天天线| 亚洲欧美日韩另类电影网站| 在线天堂最新版资源| 欧美成人精品欧美一级黄| 美女大奶头黄色视频| 夜夜骑夜夜射夜夜干| 欧美最新免费一区二区三区| av国产精品久久久久影院| a级片在线免费高清观看视频| 大码成人一级视频| 精品亚洲乱码少妇综合久久| 亚洲综合色惰| 精品亚洲成a人片在线观看| 中国三级夫妇交换| 免费黄频网站在线观看国产| 一本—道久久a久久精品蜜桃钙片| 欧美成人精品欧美一级黄| 22中文网久久字幕| 久久精品久久精品一区二区三区| 中文资源天堂在线| 成人18禁高潮啪啪吃奶动态图 | 超碰97精品在线观看| 免费高清在线观看视频在线观看| 麻豆成人午夜福利视频| 国产精品一区二区性色av| 免费观看的影片在线观看| 免费观看在线日韩| 春色校园在线视频观看| 高清午夜精品一区二区三区| 大片免费播放器 马上看| 看十八女毛片水多多多| 在线观看免费高清a一片| 伦理电影大哥的女人| 91精品国产国语对白视频| 不卡视频在线观看欧美| 免费观看性生交大片5| 伊人久久精品亚洲午夜| av在线播放精品| 国产在线免费精品| 免费观看无遮挡的男女| 人妻制服诱惑在线中文字幕| 午夜免费鲁丝| 亚洲国产日韩一区二区| 亚洲,一卡二卡三卡| 久久久国产精品麻豆| 毛片一级片免费看久久久久| 久久久国产精品麻豆| 日韩欧美精品免费久久| 国产av精品麻豆| 一级毛片我不卡| 日本黄大片高清| 乱系列少妇在线播放| 亚洲精品视频女| 性色av一级| 69精品国产乱码久久久| 国产91av在线免费观看| 99久久精品热视频| 91成人精品电影| 欧美精品国产亚洲| 免费人成在线观看视频色| 精品久久久久久久久av| 国产伦精品一区二区三区四那| 亚洲av福利一区| 在线观看免费视频网站a站| √禁漫天堂资源中文www| 午夜免费男女啪啪视频观看| 99九九线精品视频在线观看视频| 欧美人与善性xxx| av国产精品久久久久影院| 精品酒店卫生间| 插阴视频在线观看视频| 寂寞人妻少妇视频99o| 国产老妇伦熟女老妇高清| 22中文网久久字幕| 国产精品一二三区在线看| 九九在线视频观看精品| 91在线精品国自产拍蜜月| 久久女婷五月综合色啪小说| 夜夜骑夜夜射夜夜干| 成人毛片60女人毛片免费| 国产中年淑女户外野战色| 亚洲美女黄色视频免费看| 国产91av在线免费观看| av免费在线看不卡| 一级毛片aaaaaa免费看小| 亚洲欧洲精品一区二区精品久久久 | 女人久久www免费人成看片| 在线观看免费高清a一片| 国产熟女欧美一区二区| 国产淫片久久久久久久久| 91午夜精品亚洲一区二区三区| 日韩一区二区三区影片| 久久国产乱子免费精品| 亚洲美女搞黄在线观看| 最黄视频免费看| 精品99又大又爽又粗少妇毛片| 草草在线视频免费看| 观看免费一级毛片| 97在线视频观看| 人妻一区二区av| 啦啦啦中文免费视频观看日本| 国产黄片美女视频| 久久精品久久精品一区二区三区| 伊人久久精品亚洲午夜| 99久久综合免费| 日韩欧美 国产精品| 性色av一级| 一区二区av电影网| 国产日韩欧美亚洲二区| 高清av免费在线| 色5月婷婷丁香| 热re99久久国产66热| 国产成人精品无人区| 麻豆乱淫一区二区| 欧美日韩精品成人综合77777| 国产亚洲av片在线观看秒播厂| 2018国产大陆天天弄谢| 国产毛片在线视频| 一区二区三区免费毛片| 交换朋友夫妻互换小说| 妹子高潮喷水视频| 日日撸夜夜添|