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

    基于COMSOL的超聲層析成像仿真及圖像重建

    2016-09-07 06:26:50顧建飛蘇明旭蔡小舒
    聲學(xué)技術(shù) 2016年3期
    關(guān)鍵詞:層析成像二值聲場(chǎng)

    顧建飛,蘇明旭,蔡小舒

    ?

    基于COMSOL的超聲層析成像仿真及圖像重建

    顧建飛,蘇明旭,蔡小舒

    (上海理工大學(xué)顆粒與兩相流測(cè)量研究所,上海市動(dòng)力工程多相流動(dòng)與傳熱重點(diǎn)實(shí)驗(yàn)室,上海 200093)

    為了驗(yàn)證三角形和四邊形準(zhǔn)則的二值邏輯反投影(Binary Logic Back Projection,BLBP)算法在超聲層析成像重建過程中的可行性,以超聲層析成像的重建區(qū)域中同時(shí)含有圓柱形和橢圓柱形障礙物為例,利用COMSOL對(duì)成像區(qū)域內(nèi)存在障礙物的聲場(chǎng)進(jìn)行仿真,得到聲場(chǎng)聲壓分布。然后,通過提取和處理所獲取的聲壓分布,獲得信號(hào)矩陣。最后,將信號(hào)矩陣代入二值邏輯反投影算法進(jìn)行圖像重建。仿真結(jié)果表明:當(dāng)入射波的波長(zhǎng)λ遠(yuǎn)小于障礙物的尺寸時(shí),聲衍射的現(xiàn)象并不明顯。重建結(jié)果表明:兩種準(zhǔn)則在幾何近似理論下可行,且四邊形準(zhǔn)則的效果優(yōu)于三角形準(zhǔn)則。

    超聲波;有限元;圖像重建;超聲過程層析成像

    0 引言

    在冶金、化工、電力等工業(yè)過程中經(jīng)常出現(xiàn)兩相流或多相流介質(zhì),例如在冶金行業(yè)的選礦廠、冶煉廠和礦山充填工藝生產(chǎn)過程中,常常要求對(duì)漿狀的混合液中的固體(即金屬)量進(jìn)行測(cè)量,以利于成本核算、質(zhì)量控制和金屬回收。作為一種用以測(cè)量多相流流動(dòng)的過程參數(shù)的技術(shù),超聲過程層析成像技術(shù)(Ultrasonic Process Tomography,UPT)在20世紀(jì)80年代中后期逐漸形成和發(fā)展起來,并受到各界的廣泛關(guān)注[1-3]。層析成像的理論要求只有在投影角無限多,且每個(gè)投影角度下的投影線足夠多的情況下才能獲得對(duì)象的“精確重建”。目前,扇形束掃描方式由于有效投影數(shù)據(jù)量大,已經(jīng)得到了廣泛的應(yīng)用。其主要的實(shí)現(xiàn)方式是采用一種半柱面式的壓電晶片,換能器的結(jié)構(gòu)如圖1所示,最大特點(diǎn)是其聲輻射/接收的角度范圍為180°,其聲場(chǎng)可以覆蓋檢測(cè)區(qū)域。但是,超聲層析成像技術(shù)仍有如下方面的制約:超聲波發(fā)射/接收同步和掃描切換技術(shù)不夠成熟,實(shí)現(xiàn)超聲波快速成像較為困難。此外,若超聲波波長(zhǎng)與障礙物尺寸相當(dāng), 會(huì)產(chǎn)生聲散射、繞射等現(xiàn)象,且隨機(jī)發(fā)生,給測(cè)量分析帶來難度。

    COMSOL Multiphysics是一款大型數(shù)值仿真軟件,通過有限元方法模擬,在科研和工程中能用偏微分方程(Partial Differential Equation,PDE)描述各種問題。它是集前處理、求解器和后處理于一體的有限元軟件,提供友好的圖形用戶界面和建模工具,簡(jiǎn)化了建模工作。所以本文的工作將利用COMSOL軟件模擬柱面壓電換能器在超聲層析成像技術(shù)中的聲場(chǎng)。在層析成像的信號(hào)提取階段,利用COMSOL多物理場(chǎng)耦合軟件聲學(xué)模塊模擬聲場(chǎng),得到理論上陣列傳感器的聲壓分布,這對(duì)UPT實(shí)驗(yàn)系統(tǒng)的搭建具有一定的指導(dǎo)意義,為信號(hào)矩陣的精確度提供參照作用。同時(shí)也可以對(duì)改進(jìn)的二值邏輯反投影算法做進(jìn)一步說明和驗(yàn)證。

    1 模型建立

    1.1 幾何模型

    本文采用如圖2所示的二維幾何模型,建立與實(shí)驗(yàn)系統(tǒng)相同的測(cè)量區(qū)域。其中,區(qū)域C0安裝了16個(gè)換能器的陣列裝置;區(qū)域C1是半徑為55 mm的成像區(qū)域。區(qū)域C2是一個(gè)半徑為10 mm的圓形區(qū)域,圓心位置在(-0.01,0.02) m,表示的是聚四氟乙烯圓柱障礙物。區(qū)域C3是一個(gè)逆時(shí)針旋轉(zhuǎn)45°的長(zhǎng)半軸為18 mm、短半軸為12 mm、中點(diǎn)在(0,-0.02) m的聚四氟乙烯橢圓柱障礙物。障礙物聚四氟乙烯的聲速0=5200 m/s,密度0= 7865 kg/m3。障礙物外部為介質(zhì)水,聲速s=1473 m/s,密度s= 997 kg/m3。假定某時(shí)刻,根據(jù)實(shí)際測(cè)量過程中采用的扇形束掃描模式,位于最左邊的換能器為發(fā)射探頭,剩下其余的換能器為接收探頭。由于換能器晶片的發(fā)射/接收面是半徑為5 mm的半圓柱面,因此設(shè)定成像區(qū)域外邊界5 mm位置,即(-0.06,0) m處有一柱面超聲波源,頻率=2 MHz,則聲波在障礙物區(qū)域的波長(zhǎng)0=0/=2.6 mm,在介質(zhì)水的區(qū)域波長(zhǎng)s=s/= 0.736 mm。這里,計(jì)算出的波長(zhǎng)對(duì)后續(xù)設(shè)定該區(qū)域的網(wǎng)格單元尺寸尤為重要。

    1.2 超弱變分公式(Ultra Weak Variational Formulation,UWVF)模塊

    在固體和流體中,線性時(shí)諧波場(chǎng)可以由Helmholtz和Navier方程來描述。但是,由于波的振蕩特性難以對(duì)波動(dòng)方程作數(shù)值逼近。在能接受的精確度范圍內(nèi),對(duì)Helmholtz或Navier方程問題作近似求解,需要有一個(gè)相對(duì)密集的空間描述。對(duì)于標(biāo)準(zhǔn)的數(shù)值計(jì)算方法,比如低階有限差分法(Finite Difference,F(xiàn)D)和有限元法(Finite Element,F(xiàn)E),根據(jù)經(jīng)驗(yàn),每個(gè)波長(zhǎng)取10個(gè)點(diǎn)是合適的。然而高頻超聲,由于聲波波長(zhǎng)非常短,因此要求將網(wǎng)格劃分得更細(xì),這將帶來計(jì)算復(fù)雜性問題。針對(duì)高頻聲波,文中利用COMSOL獨(dú)有的超弱變分公式可以解決上述問題,和低階有限元分析相比,UWVF能顯著減少計(jì)算時(shí)間,具有明顯的優(yōu)勢(shì)。這種優(yōu)勢(shì)來自UWVF的單元屬性,UWVF單元屬性包含了求解自由空間Helmholtz波方程的信息。對(duì)于二維、三維Helmholtz問題以及二維Navier問題,該方法被證明是有效的[4]。

    1.3 網(wǎng)格劃分

    網(wǎng)格的劃分對(duì)于模擬計(jì)算的收斂性和準(zhǔn)確性有較大的影響。要給出可接受的準(zhǔn)確性,聲場(chǎng)模塊的超弱變分公式需要一個(gè)相對(duì)均勻的網(wǎng)格,并通過最大單元尺寸參數(shù)max控制,而不是通過明顯的幾何細(xì)節(jié)。COMSOL軟件的超聲散射模型指出,當(dāng)使用默認(rèn)的UWVF設(shè)置時(shí),定義的max等于或稍微大于2λ(為正整數(shù))時(shí),會(huì)在收斂性和準(zhǔn)確性之間給出較好的平衡,可以獲得較好的仿真結(jié)果。由于本文建立的模型中介質(zhì)水和障礙物區(qū)域中的波長(zhǎng)并不一致,需要在COMSOL中對(duì)這兩個(gè)不同區(qū)域設(shè)置并定義不同大小的max。根據(jù)上述原則,分別對(duì)介質(zhì)水和障礙物兩個(gè)求解域的max設(shè)置為0.0015 m和0.0052 m。網(wǎng)格劃分結(jié)果如圖3所示。

    1.4 求解過程

    模擬采用了連續(xù)波(continuous wave)形式,聲源發(fā)出的柱面波與第二類零階Hankel函數(shù)成正比:

    其中:=[(-0)2+(-0)2]1/2指成像區(qū)域內(nèi)任意一點(diǎn)坐標(biāo)(,)離聲源位置(0,0)的距離;s=2/s,表示介質(zhì)中的波數(shù);函數(shù)J和Y分別指第一、二類Bessel函數(shù),i為虛部單位。成像區(qū)域的外邊界施加輻射邊界條件,吸收柱面波。此外,要取得可接受的準(zhǔn)確性,聲場(chǎng)模塊的超弱變分公式采用默認(rèn)20階超弱Helmholtz單元。對(duì)求解域進(jìn)行求解即可獲得聲場(chǎng)分布。

    2 仿真結(jié)果分析

    本文算例的模擬結(jié)果如圖4所示,結(jié)果以聲壓級(jí)云圖表示,顯示了超聲層析成像系統(tǒng)中液態(tài)介質(zhì)和固態(tài)障礙物的聲壓級(jí)分布情況。從圖中可以看出,超聲波在水中的擴(kuò)散致使聲壓沿傳播路徑逐漸減弱。另外,每一個(gè)柱體對(duì)超聲波存在明顯的遮擋作用,柱體后向聲壓非常弱,形成陰影區(qū)。考慮到扇形束掃描的幾何近似理論:只有當(dāng)入射波的波長(zhǎng)遠(yuǎn)小于目標(biāo)物的尺寸時(shí),才能忽略由目標(biāo)引起的衍射[5],即>> 1(=2π /,表示聲波在介質(zhì)中的波數(shù);是障礙物的等效半徑)。所以,滿足該條件的2 MHz超聲波所測(cè)障礙物的等效半徑須遠(yuǎn)大于0.12 mm。由于仿真幾何模型中設(shè)定的障礙物尺寸均遠(yuǎn)大于該頻率下所測(cè)障礙物等效半徑的最小值,所以從該圖中還可以看到,聲波存在的散射現(xiàn)象,非常微弱,符合幾何近似理論。作者在實(shí)驗(yàn)中也注意到了類似現(xiàn)象,其表現(xiàn)為障礙物后方對(duì)應(yīng)的傳感器仍可以接收到極其微弱的超聲信號(hào),可以忽略不計(jì)。然后,提取和處理所獲取的該聲壓分布結(jié)果,得到理論上陣列傳感器的聲壓分布。通過對(duì)代入二值邏輯反投影算法的信號(hào)矩陣進(jìn)行閾值控制,得到一個(gè)信號(hào)矩陣(,),其中為所使用的換能器個(gè)數(shù)。那么,第個(gè)發(fā)射位置對(duì)應(yīng)的第個(gè)接收位置的信號(hào)為:

    3 圖像重建算法

    目前,采用二值邏輯反投影反演算法采用了上述的幾何聲學(xué)特性,例如基于三角形判斷的反投影反演算法[6],其主要思想如圖5所示。當(dāng)1號(hào)換能器發(fā)射時(shí),4號(hào)換能器被其傳播路徑中物體遮擋,則該探頭的信號(hào)經(jīng)過二值化處理后置0。這樣,由換能器間隔的中點(diǎn)a、b、c三點(diǎn)組成的區(qū)域表示目標(biāo)物影響區(qū)間,圖像矩陣對(duì)應(yīng)區(qū)域增加權(quán)重“1”。在此基礎(chǔ)上,本文提出了一種基于四邊形判斷的二值邏輯反投影反演算法,如圖6所示,當(dāng)1號(hào)換能器發(fā)射時(shí),4號(hào)換能器被其傳播路徑中物體遮擋,則該探頭的信號(hào)經(jīng)過二值化處理后置0。這樣,由換能器間隔的中點(diǎn)a、b、c、d四個(gè)點(diǎn)組成的四邊形區(qū)域表示目標(biāo)物影響區(qū)間,圖像矩陣對(duì)應(yīng)區(qū)域增加權(quán)重“1”,以此類推,待所有掃描信號(hào)處理完畢后通過閾值對(duì)圖像進(jìn)行目標(biāo)和背景分割,獲得重建圖像。分別按三角形判斷準(zhǔn)則和四邊形判斷準(zhǔn)則對(duì)其進(jìn)行反演,得到的結(jié)果如圖7和圖8所示。

    圖6 四邊形判斷法

    將兩種方式的重建結(jié)果與仿真對(duì)象進(jìn)行對(duì)比發(fā)現(xiàn),聲波的散射和邊緣效應(yīng)沒有對(duì)兩種算法的圖像重建帶來較大影響;障礙物的位置信息、形狀信息和尺寸信息,都較為準(zhǔn)確。說明了二值邏輯反投影算法所用的幾何聲學(xué)理論在入射波的波長(zhǎng)遠(yuǎn)小于目標(biāo)物的尺寸時(shí)具備可行性。圖7和圖8雖然可以看出兩種判斷準(zhǔn)則的重建效果大致一樣,但仔細(xì)對(duì)比上述兩種算法,可以發(fā)現(xiàn)細(xì)微的差別。主要表現(xiàn)在:兩種準(zhǔn)則都能較好地得到圓柱障礙物的截面形狀,但是四邊形準(zhǔn)則得到的圓柱障礙物的尺寸更接近真實(shí)模型,而三角形準(zhǔn)則測(cè)得的偏大;此外,四邊形準(zhǔn)則得到的橢圓柱障礙物的截面信息也稍微接近真實(shí)模型。

    4 結(jié)論

    在超聲層析成像領(lǐng)域,采用COMSOL軟件建立不同形狀、尺寸的柱體障礙物的有限元模型,并利用超弱變分公式模擬仿真成像區(qū)域中聲場(chǎng)的各點(diǎn)聲壓級(jí),提取各個(gè)換能器位置的聲壓級(jí)數(shù)值大小可獲得超聲信號(hào)矩陣。

    結(jié)合三角形判斷和四邊形判斷準(zhǔn)則的二值邏輯反投影算法的反演結(jié)果表明,兩種判斷準(zhǔn)則的重建效果大致一樣,同時(shí)也驗(yàn)證了當(dāng)入射波的波長(zhǎng)遠(yuǎn)小于目標(biāo)物時(shí)該反演算法的可行性。

    [1] Ibrahim S, Yunus M A M, Khairi M T M, et al. A review on ultrasonic process tomography system[J]. Jurnal Teknologi, 2014, 70(3): 1-5.

    [2] Abbaszadeh J, Rahim H A, Rahim R A, et al. Ultrasonic tomography system: optimizing the frequency in a metal pipe conveyor[C]//Applied Mechanics & Materials, 2013, 284: 572-576.

    [3] Hafiz M, Rahiman F, Ruzairi, et al. Ultrasonic tomography - image reconstruction algorithms[J]. International Journal of Innovative Computing Information & Control Ijicic, 2012, 8(1): 527-538.

    [4] Huttunen T. The ultra weak variational formulation for ultrasound transmission problems[M]. Kuopion yliopisto, 2004.

    [5] Xu L J, Xu L A. Ultrasound tomography system used for monitoring bubbly gas/liquid two-phase flow[J]. Ultrasonics Ferroelectrics & Frequency Control IEEE Transactions on, 1997, 44(1): 67-76.

    [6] 陳孝震, 蘇明旭, 蔡小舒. 液固兩相系中的超聲波層析成像實(shí)驗(yàn)研究[C]//中國(guó)工程熱物理學(xué)會(huì). 2012.

    CHEN Xiaozheng, SU Mingxu, CAI Xiaosu. Liquid-solid two-phase system of ultrasonic tomographic imaging experiments[C]//Chinese Society of Engineering Thermophysics. 2012.

    TheCOMSOL basedsound filed simulation and image reconstruction of ultrasonic tomography

    GU Jian-fei, SU Ming-xu, CAI Xiao-shu

    (Institute of Particle and Two-phase Flow Measurement, University of Shanghai for Science and Technology,Shanghai Key Laboratory of Multiphase Flow and Heat Transfer in Power Engineering, Shanghai 200093, China)

    In order to verify the feasibility of the Binary Logic Back Projection (BLBP) algorithm of triangular and quadrilateral criteria in ultrasonic tomography reconstruction, a typical UPT reconstruction area, which simultaneously contains a cylindrical obstacle and an elliptic cylindrical obstacle,is taken as an example. COMSOL is used to make simulation of sound field when the obstacles exist, and an acoustic pressure distribution is obtained. Then, by extracting and processing the obtained sound pressure distribution, a signal matrix is formed. Finally, the signal matrix is substituted into the binary logic back projection algorithm for image reconstruction. It could be clearly observed that the effect of diffraction is insignificant from the computation results whenthe wavelength of the incident wave is smaller than the size of the obstacle. The reconstruction results show that the two criteria are feasible under the geometric acoustic, and the reconstructed quality of the quadrilateral criterion is better than that of the triangular one.

    ultrasound; finite element; image reconstruction; Ultrasonic Process Tomography(UPT)

    TB551

    A

    1000-3630(2016)-03-0231-04

    10.16300/j.cnki.1000-3630.2016.03.009

    2015-08-31;

    2015-11-10

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51176128, 51206113)

    顧建飛(1990-), 男, 江蘇海門人, 碩士研究生, 研究方向?yàn)槌曔^程層析成像。

    蘇明旭, E-mail: sumx@usst.edu.cn

    猜你喜歡
    層析成像二值聲場(chǎng)
    混沌偽隨機(jī)二值序列的性能分析方法研究綜述
    基于大數(shù)據(jù)量的初至層析成像算法優(yōu)化
    支持CNN與LSTM的二值權(quán)重神經(jīng)網(wǎng)絡(luò)芯片
    基于快速行進(jìn)法地震層析成像研究
    基于BIM的鐵路車站聲場(chǎng)仿真分析研究
    探尋360°全聲場(chǎng)發(fā)聲門道
    基于二值形態(tài)學(xué)算子的軌道圖像分割新算法
    視頻圖像文字的二值化
    基于分布式無線網(wǎng)絡(luò)的無線電層析成像方法與實(shí)驗(yàn)研究
    基于多級(jí)小波域變換的時(shí)域擴(kuò)散熒光層析成像方法
    999久久久国产精品视频| 欧美人与性动交α欧美软件| 丝袜脚勾引网站| 三上悠亚av全集在线观看| av视频免费观看在线观看| 你懂的网址亚洲精品在线观看| 只有这里有精品99| 久久久久人妻精品一区果冻| 日韩一区二区三区影片| a级片在线免费高清观看视频| 亚洲av日韩精品久久久久久密 | 国产成人精品福利久久| 中文欧美无线码| 美女扒开内裤让男人捅视频| 国产精品久久久久久精品电影小说| videos熟女内射| 美女高潮到喷水免费观看| 在线观看免费日韩欧美大片| 黄片播放在线免费| 国产无遮挡羞羞视频在线观看| 亚洲人成网站在线观看播放| 亚洲精品日韩在线中文字幕| 久久影院123| 丝袜脚勾引网站| 国产成人系列免费观看| 黄频高清免费视频| 大片电影免费在线观看免费| 97人妻天天添夜夜摸| 巨乳人妻的诱惑在线观看| 叶爱在线成人免费视频播放| 一本—道久久a久久精品蜜桃钙片| 亚洲国产欧美网| 欧美成人精品欧美一级黄| 午夜福利,免费看| 男的添女的下面高潮视频| 国产一区亚洲一区在线观看| 男人爽女人下面视频在线观看| 免费黄频网站在线观看国产| 亚洲一码二码三码区别大吗| 大片电影免费在线观看免费| 观看美女的网站| 大香蕉久久网| 91成人精品电影| av线在线观看网站| 伊人久久国产一区二区| 伊人久久国产一区二区| 国产亚洲精品第一综合不卡| 天天影视国产精品| 亚洲国产日韩一区二区| 成人毛片60女人毛片免费| 在线天堂最新版资源| 久久天躁狠狠躁夜夜2o2o | 多毛熟女@视频| 国产麻豆69| 国产一区有黄有色的免费视频| 国产在线一区二区三区精| 91成人精品电影| 亚洲,欧美精品.| 亚洲成av片中文字幕在线观看| 国产精品久久久人人做人人爽| 少妇被粗大的猛进出69影院| 国产精品久久久人人做人人爽| 赤兔流量卡办理| 亚洲中文av在线| 一区在线观看完整版| 亚洲国产精品一区三区| 亚洲精品国产av蜜桃| 精品人妻熟女毛片av久久网站| 街头女战士在线观看网站| 两个人看的免费小视频| 女的被弄到高潮叫床怎么办| 国产男女超爽视频在线观看| av网站在线播放免费| 在线 av 中文字幕| 搡老岳熟女国产| 肉色欧美久久久久久久蜜桃| 国产激情久久老熟女| 老司机影院成人| 在线观看一区二区三区激情| 男女边摸边吃奶| 激情五月婷婷亚洲| 欧美成人精品欧美一级黄| 黄色视频不卡| 99久久精品国产亚洲精品| 久久 成人 亚洲| 久久人人97超碰香蕉20202| 亚洲中文av在线| 日韩一卡2卡3卡4卡2021年| 久久精品国产综合久久久| 女人高潮潮喷娇喘18禁视频| 天天操日日干夜夜撸| 国产激情久久老熟女| 精品少妇一区二区三区视频日本电影 | 精品一区二区三卡| 超色免费av| 激情五月婷婷亚洲| 久久久久久人人人人人| 美女高潮到喷水免费观看| 99re6热这里在线精品视频| 欧美日本中文国产一区发布| 国产精品麻豆人妻色哟哟久久| 日本欧美视频一区| 男女床上黄色一级片免费看| 亚洲成色77777| 最近最新中文字幕大全免费视频 | 国产极品粉嫩免费观看在线| 国产亚洲午夜精品一区二区久久| 国产精品 欧美亚洲| 午夜福利影视在线免费观看| 国产精品免费视频内射| 十八禁高潮呻吟视频| 国产淫语在线视频| 成年人免费黄色播放视频| 91老司机精品| 精品第一国产精品| 午夜91福利影院| 国产精品嫩草影院av在线观看| 成年av动漫网址| 亚洲精华国产精华液的使用体验| 欧美最新免费一区二区三区| 只有这里有精品99| 色吧在线观看| 久久ye,这里只有精品| 亚洲婷婷狠狠爱综合网| 欧美日韩av久久| 三上悠亚av全集在线观看| 一区二区日韩欧美中文字幕| 亚洲第一区二区三区不卡| 少妇被粗大猛烈的视频| 国产精品 欧美亚洲| 免费在线观看视频国产中文字幕亚洲 | 精品福利永久在线观看| 亚洲精品久久午夜乱码| 国产色婷婷99| av在线观看视频网站免费| 久久97久久精品| 亚洲国产精品国产精品| 欧美在线一区亚洲| 丁香六月欧美| 亚洲精品久久成人aⅴ小说| av在线观看视频网站免费| 高清不卡的av网站| 秋霞伦理黄片| 日韩不卡一区二区三区视频在线| 在线 av 中文字幕| 欧美激情极品国产一区二区三区| 男女边吃奶边做爰视频| 亚洲国产成人一精品久久久| 久久人妻熟女aⅴ| 国产精品欧美亚洲77777| 热re99久久精品国产66热6| av有码第一页| 色综合欧美亚洲国产小说| 麻豆精品久久久久久蜜桃| 成人免费观看视频高清| 日本wwww免费看| 美女主播在线视频| 青春草视频在线免费观看| 国产日韩欧美在线精品| 国产精品国产av在线观看| 99精国产麻豆久久婷婷| 亚洲色图 男人天堂 中文字幕| 国产精品欧美亚洲77777| 夜夜骑夜夜射夜夜干| 亚洲第一av免费看| 欧美人与性动交α欧美精品济南到| 精品国产国语对白av| 国产又爽黄色视频| 中文字幕人妻熟女乱码| 婷婷色综合www| 国产一卡二卡三卡精品 | 久久 成人 亚洲| 一级爰片在线观看| 国产精品国产三级专区第一集| 久久韩国三级中文字幕| 久久久欧美国产精品| 最黄视频免费看| 老司机深夜福利视频在线观看 | 美女中出高潮动态图| 两性夫妻黄色片| 日韩欧美精品免费久久| 久热这里只有精品99| 日本午夜av视频| 婷婷色综合www| 亚洲成人国产一区在线观看 | 日韩免费高清中文字幕av| 另类亚洲欧美激情| 极品人妻少妇av视频| 国产精品一区二区精品视频观看| 成人亚洲欧美一区二区av| 日本猛色少妇xxxxx猛交久久| 欧美日韩亚洲高清精品| 日韩人妻精品一区2区三区| 国产成人av激情在线播放| 中文字幕制服av| 19禁男女啪啪无遮挡网站| 成人毛片60女人毛片免费| 九色亚洲精品在线播放| 一个人免费看片子| 久久久国产精品麻豆| 国产成人91sexporn| netflix在线观看网站| 久久狼人影院| 老鸭窝网址在线观看| 男女高潮啪啪啪动态图| 少妇 在线观看| 亚洲成av片中文字幕在线观看| av免费观看日本| www.精华液| 最近手机中文字幕大全| 99久久人妻综合| 国产成人精品在线电影| 欧美人与性动交α欧美精品济南到| 一本色道久久久久久精品综合| 天天躁日日躁夜夜躁夜夜| 丝袜喷水一区| 日韩一本色道免费dvd| 中文字幕色久视频| 午夜福利视频精品| 在线亚洲精品国产二区图片欧美| 青春草视频在线免费观看| 日本欧美视频一区| 肉色欧美久久久久久久蜜桃| 亚洲精品aⅴ在线观看| 免费人妻精品一区二区三区视频| √禁漫天堂资源中文www| 国产精品久久久久久精品电影小说| 国产又爽黄色视频| 天天躁夜夜躁狠狠久久av| 国产激情久久老熟女| 亚洲精品国产av成人精品| 国产精品麻豆人妻色哟哟久久| 老汉色av国产亚洲站长工具| 中文欧美无线码| 国产激情久久老熟女| 亚洲精品第二区| 日韩 欧美 亚洲 中文字幕| 丝袜美腿诱惑在线| 日韩精品有码人妻一区| 国产精品久久久久久人妻精品电影 | 男人舔女人的私密视频| 18禁国产床啪视频网站| 久久久久久久久久久久大奶| 国产精品一区二区在线不卡| 精品一品国产午夜福利视频| 国产精品亚洲av一区麻豆 | 久久免费观看电影| 亚洲欧洲国产日韩| 国产99久久九九免费精品| 自线自在国产av| 最新的欧美精品一区二区| 夜夜骑夜夜射夜夜干| 人人妻人人爽人人添夜夜欢视频| 亚洲精品乱久久久久久| 亚洲美女搞黄在线观看| 美女主播在线视频| 国产精品国产三级国产专区5o| 叶爱在线成人免费视频播放| 午夜福利视频精品| 成年人免费黄色播放视频| 国产男女内射视频| 日韩一区二区三区影片| 蜜桃在线观看..| 欧美日韩成人在线一区二区| 亚洲 欧美一区二区三区| 女的被弄到高潮叫床怎么办| 久久久国产精品麻豆| 久久久久人妻精品一区果冻| 精品人妻在线不人妻| 日本vs欧美在线观看视频| 婷婷色av中文字幕| 最新的欧美精品一区二区| 国产精品熟女久久久久浪| 亚洲成av片中文字幕在线观看| 午夜激情久久久久久久| 丰满饥渴人妻一区二区三| 国产伦理片在线播放av一区| 亚洲av欧美aⅴ国产| 欧美日韩亚洲国产一区二区在线观看 | 女人被躁到高潮嗷嗷叫费观| 国产精品av久久久久免费| 久久精品久久精品一区二区三区| 亚洲专区中文字幕在线 | 欧美97在线视频| 日本黄色日本黄色录像| 成人国产麻豆网| 又黄又粗又硬又大视频| 不卡视频在线观看欧美| 免费av中文字幕在线| 久久久国产欧美日韩av| 一级毛片 在线播放| 中文字幕制服av| av在线观看视频网站免费| 无限看片的www在线观看| 一区福利在线观看| 精品福利永久在线观看| 国产男女内射视频| 久热爱精品视频在线9| 日本欧美视频一区| 亚洲一区二区三区欧美精品| 精品少妇内射三级| 欧美黑人精品巨大| 日本av免费视频播放| 99久久综合免费| 巨乳人妻的诱惑在线观看| 日韩av在线免费看完整版不卡| 久久久久久久精品精品| 三上悠亚av全集在线观看| 最新的欧美精品一区二区| 丰满少妇做爰视频| 老司机影院毛片| 啦啦啦中文免费视频观看日本| kizo精华| 免费av中文字幕在线| 欧美在线黄色| 亚洲av福利一区| 亚洲少妇的诱惑av| e午夜精品久久久久久久| 欧美日韩亚洲高清精品| 国产 一区精品| 国产精品一二三区在线看| 国产精品久久久久久人妻精品电影 | 国产淫语在线视频| 久久久国产一区二区| 少妇精品久久久久久久| 国产精品亚洲av一区麻豆 | svipshipincom国产片| tube8黄色片| 亚洲美女黄色视频免费看| 午夜日韩欧美国产| 在线观看人妻少妇| 亚洲av国产av综合av卡| 一边摸一边抽搐一进一出视频| 亚洲av电影在线进入| 亚洲,一卡二卡三卡| 免费在线观看完整版高清| 精品一区在线观看国产| 日韩电影二区| 国产亚洲午夜精品一区二区久久| 少妇被粗大的猛进出69影院| 国产高清国产精品国产三级| 亚洲成人国产一区在线观看 | 99热全是精品| 国产极品天堂在线| 一区二区三区激情视频| 制服诱惑二区| 中文字幕亚洲精品专区| 精品国产一区二区久久| 亚洲伊人色综图| 亚洲国产av新网站| 十八禁高潮呻吟视频| 精品久久久久久电影网| 午夜福利视频在线观看免费| 成年女人毛片免费观看观看9 | 亚洲,欧美精品.| 欧美少妇被猛烈插入视频| 久热这里只有精品99| 欧美中文综合在线视频| 精品视频人人做人人爽| 日韩精品有码人妻一区| 五月天丁香电影| 国产成人一区二区在线| 最黄视频免费看| 午夜日韩欧美国产| h视频一区二区三区| 国产一区亚洲一区在线观看| 巨乳人妻的诱惑在线观看| 在线观看三级黄色| 欧美黑人精品巨大| 制服人妻中文乱码| 精品人妻熟女毛片av久久网站| 色视频在线一区二区三区| 精品少妇黑人巨大在线播放| 国产1区2区3区精品| 久久精品国产a三级三级三级| 纵有疾风起免费观看全集完整版| 亚洲欧美中文字幕日韩二区| 人成视频在线观看免费观看| 亚洲精品成人av观看孕妇| 日韩大码丰满熟妇| 美女视频免费永久观看网站| 日本一区二区免费在线视频| 新久久久久国产一级毛片| 十分钟在线观看高清视频www| 男女午夜视频在线观看| 可以免费在线观看a视频的电影网站 | 国产麻豆69| 午夜精品国产一区二区电影| 精品久久久久久电影网| 波多野结衣av一区二区av| 中文字幕色久视频| 久热这里只有精品99| 伦理电影免费视频| 欧美人与善性xxx| 日韩成人av中文字幕在线观看| 精品少妇黑人巨大在线播放| 少妇猛男粗大的猛烈进出视频| 午夜福利免费观看在线| 精品久久久精品久久久| 国产成人精品无人区| 9热在线视频观看99| 亚洲激情五月婷婷啪啪| 国产精品国产三级国产专区5o| 91国产中文字幕| 老司机靠b影院| 69精品国产乱码久久久| 日日啪夜夜爽| 日韩中文字幕视频在线看片| 午夜福利一区二区在线看| 成人影院久久| 捣出白浆h1v1| 久久精品国产综合久久久| 中文字幕人妻丝袜制服| 日韩 亚洲 欧美在线| 亚洲精品视频女| 国产亚洲一区二区精品| 亚洲欧洲国产日韩| 一级爰片在线观看| 国产视频首页在线观看| 成年人免费黄色播放视频| 精品少妇内射三级| 大话2 男鬼变身卡| 国产亚洲av高清不卡| 青草久久国产| 99久久人妻综合| 一级毛片 在线播放| 精品一区二区三卡| 成人18禁高潮啪啪吃奶动态图| 不卡视频在线观看欧美| 丰满少妇做爰视频| 伊人亚洲综合成人网| 亚洲精品成人av观看孕妇| 国产精品一国产av| 久久精品国产a三级三级三级| 操美女的视频在线观看| 午夜影院在线不卡| 亚洲欧洲精品一区二区精品久久久 | 丝袜人妻中文字幕| 亚洲国产成人一精品久久久| 一本—道久久a久久精品蜜桃钙片| 国产一区有黄有色的免费视频| 国产激情久久老熟女| 亚洲国产精品999| 黄色毛片三级朝国网站| 天堂中文最新版在线下载| 亚洲欧美一区二区三区国产| 日本色播在线视频| 国产在线一区二区三区精| 国产成人欧美在线观看 | 国产精品偷伦视频观看了| 亚洲国产欧美日韩在线播放| 国产精品免费视频内射| 女性生殖器流出的白浆| 免费日韩欧美在线观看| 超色免费av| 亚洲精品在线美女| 日韩av在线免费看完整版不卡| 一边摸一边抽搐一进一出视频| 制服丝袜香蕉在线| 中文字幕精品免费在线观看视频| 老司机影院毛片| 国产精品久久久久成人av| 中文字幕人妻丝袜制服| 99九九在线精品视频| 99国产综合亚洲精品| 下体分泌物呈黄色| 哪个播放器可以免费观看大片| 十八禁高潮呻吟视频| 日韩av免费高清视频| 宅男免费午夜| 王馨瑶露胸无遮挡在线观看| 母亲3免费完整高清在线观看| 免费看av在线观看网站| 国产成人免费观看mmmm| 久久久欧美国产精品| 亚洲,一卡二卡三卡| 欧美日韩一区二区视频在线观看视频在线| 免费黄色在线免费观看| 看非洲黑人一级黄片| 国产精品麻豆人妻色哟哟久久| 精品一区在线观看国产| 亚洲av日韩精品久久久久久密 | 亚洲五月色婷婷综合| 69精品国产乱码久久久| 欧美精品高潮呻吟av久久| 国产精品三级大全| 岛国毛片在线播放| 国产片特级美女逼逼视频| 中文欧美无线码| 97在线人人人人妻| 欧美精品av麻豆av| 成人亚洲精品一区在线观看| 国产老妇伦熟女老妇高清| 国产精品无大码| 天天躁夜夜躁狠狠躁躁| 国产欧美日韩综合在线一区二区| 18禁动态无遮挡网站| 久久精品久久久久久噜噜老黄| 久久精品国产a三级三级三级| 亚洲精品国产一区二区精华液| 深夜精品福利| 无限看片的www在线观看| 这个男人来自地球电影免费观看 | 在线观看三级黄色| 99re6热这里在线精品视频| 美国免费a级毛片| 久久免费观看电影| www.自偷自拍.com| 欧美少妇被猛烈插入视频| 亚洲中文av在线| 曰老女人黄片| 操出白浆在线播放| 秋霞伦理黄片| 91精品伊人久久大香线蕉| 中文天堂在线官网| 欧美日韩亚洲高清精品| 欧美日韩亚洲国产一区二区在线观看 | 日韩 欧美 亚洲 中文字幕| 国产成人av激情在线播放| 精品国产乱码久久久久久男人| 麻豆乱淫一区二区| 亚洲伊人久久精品综合| 亚洲欧美成人精品一区二区| 精品久久久精品久久久| 大香蕉久久成人网| 天天躁狠狠躁夜夜躁狠狠躁| 久久人人爽av亚洲精品天堂| 国产伦人伦偷精品视频| 国产精品久久久久久精品古装| 一级毛片 在线播放| 久久亚洲国产成人精品v| 欧美乱码精品一区二区三区| 久久99精品国语久久久| 午夜福利视频在线观看免费| 91老司机精品| 大话2 男鬼变身卡| 亚洲国产成人一精品久久久| 亚洲欧洲日产国产| 午夜福利在线免费观看网站| 国产日韩欧美视频二区| 菩萨蛮人人尽说江南好唐韦庄| 两个人免费观看高清视频| 久久久久久久精品精品| 最近最新中文字幕大全免费视频 | 女性生殖器流出的白浆| 国产黄色视频一区二区在线观看| 久久这里只有精品19| 免费看不卡的av| 女人爽到高潮嗷嗷叫在线视频| 在线观看免费日韩欧美大片| 一级爰片在线观看| 亚洲精品av麻豆狂野| 老鸭窝网址在线观看| 老熟女久久久| 久久影院123| 中文字幕人妻丝袜制服| 丝袜在线中文字幕| www.精华液| 人人妻人人澡人人爽人人夜夜| 国产 精品1| 人人妻人人澡人人看| 80岁老熟妇乱子伦牲交| 日韩欧美精品免费久久| 欧美精品一区二区免费开放| 欧美97在线视频| 少妇猛男粗大的猛烈进出视频| 日韩av在线免费看完整版不卡| 极品少妇高潮喷水抽搐| 观看美女的网站| 亚洲精品aⅴ在线观看| 九九爱精品视频在线观看| 成人影院久久| 久久人妻熟女aⅴ| 亚洲国产成人一精品久久久| 丰满乱子伦码专区| 国产有黄有色有爽视频| 国产精品久久久久久久久免| 欧美中文综合在线视频| 国产高清不卡午夜福利| 在线观看www视频免费| 午夜激情久久久久久久| 亚洲成色77777| 国产亚洲一区二区精品| 一级片'在线观看视频| 亚洲成色77777| 国产一区亚洲一区在线观看| 黄色一级大片看看| 韩国av在线不卡| 少妇的丰满在线观看| www.av在线官网国产| 国产 精品1| 日日撸夜夜添| 亚洲男人天堂网一区| 老司机亚洲免费影院| 看非洲黑人一级黄片| 国产又色又爽无遮挡免| 人妻人人澡人人爽人人| 欧美人与性动交α欧美精品济南到| 午夜免费观看性视频| 美国免费a级毛片| 欧美人与性动交α欧美精品济南到| 久久久精品免费免费高清| av电影中文网址| 波多野结衣一区麻豆| 午夜福利一区二区在线看| 中文字幕精品免费在线观看视频| 欧美日韩一级在线毛片| 赤兔流量卡办理| 欧美 亚洲 国产 日韩一| 日韩精品有码人妻一区| 欧美激情 高清一区二区三区| 亚洲四区av| 亚洲国产av新网站| 国产探花极品一区二区| 免费黄色在线免费观看|