• <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年已治理達標堤防長度
    2018国产大陆天天弄谢| 建设人人有责人人尽责人人享有的| 日日爽夜夜爽网站| 色哟哟·www| 久久人人爽人人爽人人片va| 黄片无遮挡物在线观看| 国产精品福利在线免费观看| 欧美精品人与动牲交sv欧美| 中文字幕人妻丝袜制服| 国产免费福利视频在线观看| 全区人妻精品视频| 亚洲成人av在线免费| 国产欧美亚洲国产| 黄色配什么色好看| 欧美激情极品国产一区二区三区 | 男女啪啪激烈高潮av片| 精品人妻偷拍中文字幕| 三级国产精品欧美在线观看| 国产黄片美女视频| 国产午夜精品一二区理论片| tube8黄色片| 国产亚洲精品久久久com| 99九九线精品视频在线观看视频| 亚洲精品国产成人久久av| 久久久精品94久久精品| 多毛熟女@视频| 久久久久久久久大av| 少妇丰满av| 一区在线观看完整版| 久久国产乱子免费精品| 亚洲美女黄色视频免费看| 久久99一区二区三区| 午夜福利在线观看免费完整高清在| 少妇高潮的动态图| 不卡视频在线观看欧美| 欧美高清成人免费视频www| 日韩一区二区视频免费看| 99久久精品一区二区三区| 人妻 亚洲 视频| 久久这里有精品视频免费| 色5月婷婷丁香| 日韩成人伦理影院| 国产毛片在线视频| 内射极品少妇av片p| 精品人妻熟女毛片av久久网站| 高清视频免费观看一区二区| 另类亚洲欧美激情| 国产又色又爽无遮挡免| 免费看不卡的av| 国产又色又爽无遮挡免| 麻豆成人av视频| 国产在线视频一区二区| 伦精品一区二区三区| 欧美日本中文国产一区发布| 看免费成人av毛片| 国产男女超爽视频在线观看| 女性生殖器流出的白浆| 成人亚洲精品一区在线观看| .国产精品久久| 亚洲成人av在线免费| 国产精品久久久久久久电影| 三级经典国产精品| 久久久国产精品麻豆| av在线老鸭窝| 免费在线观看成人毛片| 少妇猛男粗大的猛烈进出视频| 观看av在线不卡| 久久鲁丝午夜福利片| 国产一级毛片在线| 啦啦啦中文免费视频观看日本| 国产精品成人在线| 久久久精品免费免费高清| 中文乱码字字幕精品一区二区三区| av免费在线看不卡| 国产精品久久久久久久电影| 久久99热6这里只有精品| 亚洲精品国产色婷婷电影| 汤姆久久久久久久影院中文字幕| 男女无遮挡免费网站观看| 国内精品宾馆在线| 五月玫瑰六月丁香| 这个男人来自地球电影免费观看 | a级一级毛片免费在线观看| 青春草视频在线免费观看| 男人添女人高潮全过程视频| 观看免费一级毛片| 国产亚洲av片在线观看秒播厂| 亚洲国产毛片av蜜桃av| 国产高清三级在线| 大又大粗又爽又黄少妇毛片口| 欧美精品高潮呻吟av久久| 国产亚洲av片在线观看秒播厂| 老女人水多毛片| 内地一区二区视频在线| 久久精品国产亚洲av天美| 如何舔出高潮| 嘟嘟电影网在线观看| 国产亚洲一区二区精品| 亚洲不卡免费看| 日本欧美视频一区| 97在线人人人人妻| 成年av动漫网址| 久久精品国产亚洲av天美| 伊人久久国产一区二区| 国产高清三级在线| 欧美精品高潮呻吟av久久| 久久久国产欧美日韩av| av在线app专区| 成年av动漫网址| 久久久久久人妻| 十八禁高潮呻吟视频 | 国产黄频视频在线观看| 亚洲精品日韩av片在线观看| 亚洲欧美成人综合另类久久久| 国产精品秋霞免费鲁丝片| 一级毛片 在线播放| 亚洲内射少妇av| 久久久久人妻精品一区果冻| 天美传媒精品一区二区| 热re99久久国产66热| 最后的刺客免费高清国语| 狂野欧美激情性bbbbbb| 亚洲av不卡在线观看| av在线app专区| 国产av一区二区精品久久| 人妻一区二区av| 中文字幕人妻丝袜制服| 男人舔奶头视频| 国产熟女午夜一区二区三区 | 国产欧美日韩综合在线一区二区 | 中文字幕免费在线视频6| 亚洲欧美清纯卡通| 大香蕉97超碰在线| 一级黄片播放器| 国产乱人偷精品视频| 中文字幕久久专区| 亚洲国产精品999| 黑丝袜美女国产一区| 国内少妇人妻偷人精品xxx网站| 五月开心婷婷网| 一级毛片电影观看| 一区二区三区乱码不卡18| 日韩欧美精品免费久久| 久久人人爽人人爽人人片va| 亚洲av欧美aⅴ国产| 精品酒店卫生间| 99久久精品一区二区三区| 丝瓜视频免费看黄片| 国产精品99久久99久久久不卡 | 伊人久久精品亚洲午夜| 你懂的网址亚洲精品在线观看| 亚洲精品乱码久久久久久按摩| 最新中文字幕久久久久| 国产精品免费大片| 观看美女的网站| 热99国产精品久久久久久7| 国产精品国产三级国产专区5o| 18禁动态无遮挡网站| 亚洲精品国产av蜜桃| 国产色婷婷99| 五月开心婷婷网| 纵有疾风起免费观看全集完整版| 日韩人妻高清精品专区| 亚洲综合色惰| 精品久久国产蜜桃| 免费观看性生交大片5| 性色avwww在线观看| 国产无遮挡羞羞视频在线观看| 一级爰片在线观看| 国产欧美亚洲国产| 高清黄色对白视频在线免费看 | 欧美xxxx性猛交bbbb| 乱人伦中国视频| 日韩在线高清观看一区二区三区| 黄色日韩在线| 午夜日本视频在线| 午夜av观看不卡| 亚洲av福利一区| 国产av国产精品国产| 一区二区三区乱码不卡18| 麻豆乱淫一区二区| av线在线观看网站| 国产成人一区二区在线| 国产黄片美女视频| 日本与韩国留学比较| 亚洲一级一片aⅴ在线观看| 亚洲情色 制服丝袜| 欧美 亚洲 国产 日韩一| 亚洲精品中文字幕在线视频 | 欧美xxⅹ黑人| 成人黄色视频免费在线看| 又粗又硬又长又爽又黄的视频| 免费观看无遮挡的男女| 在线 av 中文字幕| 国产爽快片一区二区三区| 18+在线观看网站| 久久99一区二区三区| 99热国产这里只有精品6| 国产日韩一区二区三区精品不卡 | av免费在线看不卡| 色视频在线一区二区三区| 国产片特级美女逼逼视频| 亚洲美女搞黄在线观看| 尾随美女入室| 亚洲精品aⅴ在线观看| 亚洲美女视频黄频| 青春草亚洲视频在线观看| 国产精品久久久久久精品电影小说| 99热这里只有精品一区| 另类亚洲欧美激情| 国产免费福利视频在线观看| 成人黄色视频免费在线看| av天堂久久9| 国产亚洲精品久久久com| 久久久久久人妻| 内射极品少妇av片p| 日韩一本色道免费dvd| 各种免费的搞黄视频| 成人无遮挡网站| 特大巨黑吊av在线直播| a级片在线免费高清观看视频| 亚洲欧美精品自产自拍| 搡女人真爽免费视频火全软件| 人人澡人人妻人| 亚洲欧美中文字幕日韩二区| 丁香六月天网| h日本视频在线播放| 只有这里有精品99| 国产男女内射视频| 最近手机中文字幕大全| 纯流量卡能插随身wifi吗| 亚洲欧美成人精品一区二区| 黑人巨大精品欧美一区二区蜜桃 | 国产精品秋霞免费鲁丝片| 国产高清国产精品国产三级| 亚洲四区av| 免费观看av网站的网址| 亚洲精品国产色婷婷电影| 成人无遮挡网站| 久久久久久久大尺度免费视频| 国产伦精品一区二区三区四那| .国产精品久久| 精品久久久久久久久亚洲| 国产精品一区二区在线不卡| 男人和女人高潮做爰伦理| 日日啪夜夜爽| 五月开心婷婷网| 最黄视频免费看| 久久久久人妻精品一区果冻| 色婷婷久久久亚洲欧美| 久久久久久久大尺度免费视频| 亚洲性久久影院| 天天操日日干夜夜撸| 亚洲一级一片aⅴ在线观看| 精品国产国语对白av| 狠狠精品人妻久久久久久综合| 伦精品一区二区三区| 汤姆久久久久久久影院中文字幕| 国产亚洲欧美精品永久| 久久6这里有精品| 亚洲精品第二区| 国产伦精品一区二区三区四那| 久久久久久久久久久久大奶| 国产高清国产精品国产三级| 狂野欧美激情性bbbbbb| tube8黄色片| 国产精品一区二区三区四区免费观看| 亚洲真实伦在线观看| av天堂久久9| 欧美丝袜亚洲另类| 欧美高清成人免费视频www| 亚洲av中文av极速乱| 男女国产视频网站| 少妇 在线观看| 赤兔流量卡办理| 日韩av不卡免费在线播放| 亚洲精品乱久久久久久| 欧美日韩国产mv在线观看视频| 久久久久久人妻| 国产黄片美女视频| 国内揄拍国产精品人妻在线| 中文字幕av电影在线播放| 狠狠精品人妻久久久久久综合| 十八禁网站网址无遮挡 | 精品少妇久久久久久888优播| 日本黄色片子视频| 人人妻人人看人人澡| 欧美亚洲 丝袜 人妻 在线| a级毛片免费高清观看在线播放| 国产精品久久久久久精品电影小说| 日本欧美视频一区| a级毛片在线看网站| 欧美国产精品一级二级三级 | 夜夜看夜夜爽夜夜摸| 日韩三级伦理在线观看| 高清欧美精品videossex| 99热这里只有是精品50| 永久免费av网站大全| 91久久精品国产一区二区三区| 99热全是精品| 人人妻人人看人人澡| 中文欧美无线码| 人人妻人人添人人爽欧美一区卜| 亚洲精品456在线播放app| 日本黄色日本黄色录像| 精品国产一区二区久久| 国产免费视频播放在线视频| 欧美成人午夜免费资源| 成人亚洲欧美一区二区av| 亚洲精品乱久久久久久| 亚洲精品日本国产第一区| 80岁老熟妇乱子伦牲交| 亚洲精品乱码久久久久久按摩| 亚洲精品国产色婷婷电影| 少妇熟女欧美另类| 亚洲性久久影院| 五月伊人婷婷丁香| 六月丁香七月| 国产在线一区二区三区精| 亚洲av日韩在线播放| 欧美少妇被猛烈插入视频| 国产中年淑女户外野战色| 欧美+日韩+精品| 视频区图区小说| 热99国产精品久久久久久7| 男女无遮挡免费网站观看| 日韩av在线免费看完整版不卡| 久久综合国产亚洲精品| 成人二区视频| 三级经典国产精品| 18禁裸乳无遮挡动漫免费视频| 中文乱码字字幕精品一区二区三区| 久久精品熟女亚洲av麻豆精品| 两个人的视频大全免费| 男女边吃奶边做爰视频| 久久99热6这里只有精品| 国产在线视频一区二区| 亚洲国产精品专区欧美| 成人无遮挡网站| 免费人成在线观看视频色| 日日啪夜夜撸| 精品少妇久久久久久888优播| 乱人伦中国视频| 我的老师免费观看完整版| 久久久久人妻精品一区果冻| 久久久久久久大尺度免费视频| 在线天堂最新版资源| 国产精品欧美亚洲77777| 色吧在线观看| 国产亚洲精品久久久com| 久久狼人影院| 亚洲精品,欧美精品| 日本爱情动作片www.在线观看| 老熟女久久久| 国产成人精品久久久久久| 搡老乐熟女国产| 九九爱精品视频在线观看| 老司机影院毛片| 国产精品免费大片| 成人二区视频| 久久国产精品大桥未久av | 精品亚洲成国产av| 亚洲欧美一区二区三区黑人 | 性色avwww在线观看| 国产午夜精品一二区理论片| 亚洲精品日韩av片在线观看| 王馨瑶露胸无遮挡在线观看| 丝袜喷水一区| 国产免费一级a男人的天堂| 免费大片18禁| 新久久久久国产一级毛片| 亚洲精品aⅴ在线观看| 欧美少妇被猛烈插入视频| 久热久热在线精品观看| 日韩熟女老妇一区二区性免费视频| 亚洲一级一片aⅴ在线观看| 精品久久久精品久久久| 亚洲精品456在线播放app| 国产熟女欧美一区二区| 日韩三级伦理在线观看| 黄色视频在线播放观看不卡| 成年av动漫网址| 熟女电影av网| 最黄视频免费看| 国产综合精华液| 亚洲精品国产av蜜桃| 国产成人精品久久久久久| 9色porny在线观看| 亚洲图色成人| 91久久精品电影网| 精品少妇久久久久久888优播| 精品99又大又爽又粗少妇毛片| 人体艺术视频欧美日本| 一级a做视频免费观看| 特大巨黑吊av在线直播| 寂寞人妻少妇视频99o| 欧美另类一区| 嘟嘟电影网在线观看| 大话2 男鬼变身卡| 2018国产大陆天天弄谢| 亚洲国产成人一精品久久久| 男女免费视频国产| a级毛片免费高清观看在线播放| 欧美 日韩 精品 国产| 熟女电影av网| 亚洲情色 制服丝袜| 内地一区二区视频在线| 91久久精品电影网| 亚洲精品久久午夜乱码| 又黄又爽又刺激的免费视频.| 国产黄色免费在线视频| 亚洲成人手机| 亚洲三级黄色毛片| 国产精品成人在线| 夜夜骑夜夜射夜夜干| 国产伦在线观看视频一区| 国产成人免费观看mmmm| 男人和女人高潮做爰伦理| 亚洲,欧美,日韩| 国产精品人妻久久久影院| 欧美 日韩 精品 国产| 亚洲av欧美aⅴ国产| 91精品国产国语对白视频| 亚洲人成网站在线观看播放| 国产一区二区在线观看日韩| 日韩 亚洲 欧美在线| 老司机影院毛片| 青青草视频在线视频观看| 国产高清三级在线| 日本午夜av视频| av女优亚洲男人天堂| 国产探花极品一区二区| 国产精品国产三级国产av玫瑰| 欧美丝袜亚洲另类| 永久免费av网站大全| 国产av国产精品国产| 日韩精品免费视频一区二区三区 | 亚洲高清免费不卡视频| 美女中出高潮动态图| 国产极品粉嫩免费观看在线 | 噜噜噜噜噜久久久久久91| 国产亚洲av片在线观看秒播厂| 欧美日韩在线观看h| 欧美 亚洲 国产 日韩一| 99国产精品免费福利视频| 777米奇影视久久| 国产爽快片一区二区三区| 汤姆久久久久久久影院中文字幕| 亚洲va在线va天堂va国产| 不卡视频在线观看欧美| 夫妻午夜视频| 又黄又爽又刺激的免费视频.| 欧美日韩视频精品一区| 老司机亚洲免费影院| 免费大片黄手机在线观看| 欧美少妇被猛烈插入视频| 久久久a久久爽久久v久久| 亚洲精品日韩av片在线观看| 国产精品一区www在线观看| 最近中文字幕2019免费版| 三上悠亚av全集在线观看 | 少妇人妻 视频| 国产老妇伦熟女老妇高清| 久久久国产精品麻豆| 黄色怎么调成土黄色| 亚洲人成网站在线观看播放| 日韩大片免费观看网站| 国产无遮挡羞羞视频在线观看| 男女免费视频国产| 在线观看免费高清a一片| 国产又色又爽无遮挡免| 爱豆传媒免费全集在线观看| 美女视频免费永久观看网站| 亚洲av日韩在线播放| 成人国产av品久久久| 国内精品宾馆在线| 99国产精品免费福利视频| 国产男女内射视频| 丝瓜视频免费看黄片| 91aial.com中文字幕在线观看| 春色校园在线视频观看| 欧美日本中文国产一区发布| 亚洲精品国产色婷婷电影| 欧美精品高潮呻吟av久久| 视频区图区小说| 国产无遮挡羞羞视频在线观看| 久久久久久久久久久久大奶| 最后的刺客免费高清国语| 国产一区二区三区综合在线观看 | 亚洲精品一二三| 国产国拍精品亚洲av在线观看| 精品亚洲成国产av| 日韩中字成人| 日本黄色片子视频| 成人毛片a级毛片在线播放| 男人爽女人下面视频在线观看| 永久网站在线| 热99国产精品久久久久久7| 制服丝袜香蕉在线| 性色avwww在线观看| 插逼视频在线观看| 国语对白做爰xxxⅹ性视频网站| 久久久久久久久久人人人人人人| 在线观看美女被高潮喷水网站| 赤兔流量卡办理| 777米奇影视久久| 在线观看人妻少妇| 久久精品国产自在天天线| 在线观看www视频免费| 22中文网久久字幕| 免费观看的影片在线观看| 另类亚洲欧美激情| 国产男女内射视频| av在线老鸭窝| 亚洲av.av天堂| 老熟女久久久| 国产极品粉嫩免费观看在线 | 午夜免费鲁丝| 亚洲精品一二三| 夫妻午夜视频| 国产成人精品一,二区| 极品人妻少妇av视频| 国产精品久久久久久久久免| 亚洲一区二区三区欧美精品| 麻豆成人午夜福利视频| 亚洲欧美精品自产自拍| 亚洲综合精品二区| 精品国产露脸久久av麻豆| 久久97久久精品| 中国美白少妇内射xxxbb| 国产精品三级大全| 国产色婷婷99| 大码成人一级视频| 少妇精品久久久久久久| 亚洲av.av天堂| 亚洲精品一区蜜桃| 看非洲黑人一级黄片| 晚上一个人看的免费电影| 日本与韩国留学比较| 国产乱来视频区| 日韩av不卡免费在线播放| 男人爽女人下面视频在线观看| 亚洲国产精品一区二区三区在线| 中文字幕精品免费在线观看视频 | 大香蕉97超碰在线| 美女脱内裤让男人舔精品视频| 国产精品久久久久久av不卡| 国产成人freesex在线| 亚洲,一卡二卡三卡| 亚洲欧美日韩东京热| 精品少妇久久久久久888优播| 人妻一区二区av| 人人澡人人妻人| 欧美人与善性xxx| 日本爱情动作片www.在线观看| 成年人免费黄色播放视频 | 少妇人妻精品综合一区二区| 久久久久久久久久久久大奶| 水蜜桃什么品种好| 两个人的视频大全免费| 亚洲国产精品成人久久小说| 欧美精品一区二区大全| 亚洲国产精品成人久久小说| 久久热精品热| 欧美xxxx性猛交bbbb| 一级,二级,三级黄色视频| 久久久久国产精品人妻一区二区| 日韩av在线免费看完整版不卡| 亚洲精品成人av观看孕妇| 日日啪夜夜撸| 丰满饥渴人妻一区二区三| 午夜av观看不卡| 国产一级毛片在线| 日本色播在线视频| 特大巨黑吊av在线直播| 少妇被粗大猛烈的视频| 亚洲人成网站在线播| 99热国产这里只有精品6| 高清毛片免费看| 丰满人妻一区二区三区视频av| 又黄又爽又刺激的免费视频.| 中文精品一卡2卡3卡4更新| 久久狼人影院| 97在线视频观看| 欧美高清成人免费视频www| 少妇裸体淫交视频免费看高清| 成年人免费黄色播放视频 | 熟妇人妻不卡中文字幕| 九九久久精品国产亚洲av麻豆| 中文在线观看免费www的网站| 国产精品人妻久久久影院| a级一级毛片免费在线观看| 一级毛片电影观看| 国产成人午夜福利电影在线观看| 国产日韩一区二区三区精品不卡 | 男人添女人高潮全过程视频| 久久久久人妻精品一区果冻| 日韩,欧美,国产一区二区三区| 99久久综合免费| 精品国产露脸久久av麻豆| 日本欧美国产在线视频| 天天操日日干夜夜撸| 少妇 在线观看| 观看av在线不卡| 伊人久久国产一区二区| 免费不卡的大黄色大毛片视频在线观看| 免费少妇av软件| 成人特级av手机在线观看| 亚洲精品一二三| 国产av国产精品国产| 免费大片黄手机在线观看| 一级毛片久久久久久久久女| 欧美老熟妇乱子伦牲交| 亚洲色图综合在线观看| 欧美日韩综合久久久久久|