• <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ò)散熒光層析成像方法
    国产高清视频在线观看网站| av黄色大香蕉| 热99re8久久精品国产| 三级毛片av免费| 久久久国产成人免费| 欧美色视频一区免费| 欧美zozozo另类| 色精品久久人妻99蜜桃| 在线观看午夜福利视频| 国产黄色小视频在线观看| 少妇裸体淫交视频免费看高清| 亚洲欧美日韩卡通动漫| 精品久久久久久久末码| 亚洲精品久久国产高清桃花| 欧美xxxx黑人xx丫x性爽| 午夜福利在线观看吧| 99国产精品一区二区蜜桃av| 天天躁日日操中文字幕| 日本免费a在线| cao死你这个sao货| 中亚洲国语对白在线视频| 欧美日韩乱码在线| 性色av乱码一区二区三区2| 亚洲狠狠婷婷综合久久图片| 亚洲国产欧美人成| 99国产极品粉嫩在线观看| 很黄的视频免费| 久久国产乱子伦精品免费另类| 免费看日本二区| 久久精品91蜜桃| 国产精华一区二区三区| av国产免费在线观看| 又黄又粗又硬又大视频| xxx96com| aaaaa片日本免费| 国产精品久久电影中文字幕| 级片在线观看| 久久久精品大字幕| 久久久久久久久中文| 久久精品国产综合久久久| 观看免费一级毛片| 欧美中文日本在线观看视频| 天堂影院成人在线观看| 日本一本二区三区精品| 精品一区二区三区视频在线 | 美女午夜性视频免费| 美女午夜性视频免费| 一进一出好大好爽视频| 88av欧美| 国产亚洲欧美在线一区二区| 免费大片18禁| av视频在线观看入口| 99riav亚洲国产免费| 国产成人一区二区三区免费视频网站| 欧美中文日本在线观看视频| 国产精品久久电影中文字幕| 日本一本二区三区精品| 丁香六月欧美| 又黄又粗又硬又大视频| 亚洲av五月六月丁香网| 变态另类成人亚洲欧美熟女| 国产精品99久久99久久久不卡| 免费av不卡在线播放| 99久久久亚洲精品蜜臀av| 村上凉子中文字幕在线| 国产真人三级小视频在线观看| 日韩欧美一区二区三区在线观看| 91麻豆精品激情在线观看国产| 国产黄a三级三级三级人| 亚洲无线在线观看| 少妇裸体淫交视频免费看高清| 亚洲五月天丁香| 后天国语完整版免费观看| x7x7x7水蜜桃| 特大巨黑吊av在线直播| 亚洲男人的天堂狠狠| 国产 一区 欧美 日韩| 亚洲 欧美一区二区三区| 看黄色毛片网站| 一个人看视频在线观看www免费 | 中文字幕av在线有码专区| 成人高潮视频无遮挡免费网站| 亚洲色图av天堂| 亚洲色图av天堂| 身体一侧抽搐| 国产野战对白在线观看| 少妇丰满av| 麻豆国产97在线/欧美| 可以在线观看的亚洲视频| 亚洲成人久久爱视频| 在线播放国产精品三级| 亚洲av片天天在线观看| h日本视频在线播放| 亚洲 国产 在线| 美女大奶头视频| 亚洲片人在线观看| 波多野结衣高清作品| 久久中文字幕一级| 99久久无色码亚洲精品果冻| 亚洲成av人片免费观看| 精品无人区乱码1区二区| 国产激情久久老熟女| 高清在线国产一区| 亚洲 欧美一区二区三区| 亚洲午夜精品一区,二区,三区| 久久精品91蜜桃| 亚洲精品一区av在线观看| 亚洲,欧美精品.| 极品教师在线免费播放| 欧美zozozo另类| 国产亚洲精品一区二区www| 亚洲,欧美精品.| 精品久久久久久久毛片微露脸| av视频在线观看入口| 男人和女人高潮做爰伦理| 久久久久性生活片| 精品99又大又爽又粗少妇毛片 | 午夜激情欧美在线| 亚洲欧美精品综合久久99| 91老司机精品| 欧美乱色亚洲激情| 一级黄色大片毛片| 亚洲成人久久性| 国内揄拍国产精品人妻在线| 一二三四社区在线视频社区8| 观看美女的网站| 久9热在线精品视频| www.精华液| 国产亚洲精品av在线| 999久久久精品免费观看国产| 日韩欧美精品v在线| 亚洲精品色激情综合| 久久精品综合一区二区三区| 给我免费播放毛片高清在线观看| 国产v大片淫在线免费观看| 啦啦啦韩国在线观看视频| 老司机在亚洲福利影院| 一进一出抽搐gif免费好疼| 美女大奶头视频| 日韩三级视频一区二区三区| 噜噜噜噜噜久久久久久91| 后天国语完整版免费观看| 成人三级做爰电影| 免费一级毛片在线播放高清视频| 国产熟女xx| 999精品在线视频| 国产不卡一卡二| 天天添夜夜摸| 欧美色视频一区免费| 色av中文字幕| 日韩欧美在线乱码| 日韩中文字幕欧美一区二区| 宅男免费午夜| 国产91精品成人一区二区三区| 国产亚洲欧美在线一区二区| 成人午夜高清在线视频| 国产一区二区激情短视频| 日本a在线网址| 久久久国产成人免费| 在线播放国产精品三级| 日韩精品中文字幕看吧| 日韩人妻高清精品专区| 久久精品国产亚洲av香蕉五月| 国产精品精品国产色婷婷| 后天国语完整版免费观看| 欧美一级毛片孕妇| 啦啦啦韩国在线观看视频| 国产精品99久久99久久久不卡| 亚洲美女黄片视频| 免费高清视频大片| 日本 欧美在线| 成年女人毛片免费观看观看9| 亚洲国产欧美网| 999精品在线视频| 国产不卡一卡二| 午夜亚洲福利在线播放| 非洲黑人性xxxx精品又粗又长| 久久热在线av| 国产精品99久久99久久久不卡| 国产精品久久久av美女十八| 麻豆国产97在线/欧美| 精品欧美国产一区二区三| 久久伊人香网站| 免费一级毛片在线播放高清视频| 精品国产亚洲在线| 99re在线观看精品视频| 噜噜噜噜噜久久久久久91| 亚洲成人久久爱视频| 一级毛片女人18水好多| 日本成人三级电影网站| 国产一区二区在线观看日韩 | 美女免费视频网站| 嫩草影院精品99| 狂野欧美白嫩少妇大欣赏| 夜夜躁狠狠躁天天躁| 网址你懂的国产日韩在线| 欧美av亚洲av综合av国产av| 欧美性猛交╳xxx乱大交人| 欧美日韩福利视频一区二区| 国产 一区 欧美 日韩| 国产毛片a区久久久久| 免费看a级黄色片| 成年女人永久免费观看视频| 这个男人来自地球电影免费观看| 男插女下体视频免费在线播放| 成人鲁丝片一二三区免费| 国产伦精品一区二区三区四那| 亚洲中文字幕日韩| 草草在线视频免费看| 色播亚洲综合网| 中出人妻视频一区二区| 中文字幕高清在线视频| 日韩欧美 国产精品| 美女高潮喷水抽搐中文字幕| 中文字幕熟女人妻在线| 岛国在线免费视频观看| 成人鲁丝片一二三区免费| 精品久久久久久久人妻蜜臀av| 欧美乱色亚洲激情| 小说图片视频综合网站| 亚洲国产精品sss在线观看| 中文字幕人成人乱码亚洲影| 精品久久久久久久毛片微露脸| 99久久国产精品久久久| 最近最新中文字幕大全电影3| 国产视频一区二区在线看| av欧美777| 91在线精品国自产拍蜜月 | x7x7x7水蜜桃| 中文字幕最新亚洲高清| 国产精品一区二区精品视频观看| 国产主播在线观看一区二区| 亚洲国产精品成人综合色| 狂野欧美白嫩少妇大欣赏| 免费一级毛片在线播放高清视频| 欧美日韩黄片免| 久久久水蜜桃国产精品网| 1024手机看黄色片| 天堂影院成人在线观看| 亚洲人成网站在线播放欧美日韩| 国产精品一及| 午夜福利视频1000在线观看| 男女那种视频在线观看| 麻豆成人av在线观看| 一区福利在线观看| 69av精品久久久久久| 人妻久久中文字幕网| 又黄又爽又免费观看的视频| 久久久国产欧美日韩av| 制服丝袜大香蕉在线| 欧美成人免费av一区二区三区| 国内精品一区二区在线观看| 日韩 欧美 亚洲 中文字幕| 国产成人精品久久二区二区91| 精品久久蜜臀av无| 国产精品99久久久久久久久| 女人高潮潮喷娇喘18禁视频| 欧美日韩乱码在线| 欧美成人免费av一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 久久国产乱子伦精品免费另类| 搡老妇女老女人老熟妇| 亚洲成a人片在线一区二区| 天堂av国产一区二区熟女人妻| 精品人妻1区二区| 国产精品av久久久久免费| 色视频www国产| 99国产精品一区二区三区| 久久这里只有精品中国| 天天躁日日操中文字幕| 国产精品永久免费网站| 亚洲成人免费电影在线观看| 哪里可以看免费的av片| 亚洲色图 男人天堂 中文字幕| 免费观看精品视频网站| 男人舔奶头视频| 香蕉丝袜av| 蜜桃久久精品国产亚洲av| 亚洲精品一卡2卡三卡4卡5卡| 久久久成人免费电影| 俺也久久电影网| 国产av麻豆久久久久久久| 国产精品久久久久久亚洲av鲁大| 极品教师在线免费播放| 亚洲欧美精品综合久久99| 欧美日韩精品网址| 91麻豆精品激情在线观看国产| 桃红色精品国产亚洲av| 男女那种视频在线观看| 1000部很黄的大片| 色综合站精品国产| 精品一区二区三区视频在线 | 在线国产一区二区在线| 亚洲av成人不卡在线观看播放网| 亚洲狠狠婷婷综合久久图片| 成年免费大片在线观看| 成人高潮视频无遮挡免费网站| 日本 av在线| 日本一二三区视频观看| 国产精品野战在线观看| 久久久国产欧美日韩av| 一个人观看的视频www高清免费观看 | 两个人视频免费观看高清| 欧美另类亚洲清纯唯美| 黑人巨大精品欧美一区二区mp4| 亚洲av成人av| 很黄的视频免费| 日韩国内少妇激情av| 久久天堂一区二区三区四区| 一级作爱视频免费观看| 欧美高清成人免费视频www| 女人高潮潮喷娇喘18禁视频| 国产成人精品久久二区二区91| 国产人伦9x9x在线观看| 99久久久亚洲精品蜜臀av| 国产精品av久久久久免费| 亚洲精品乱码久久久v下载方式 | 亚洲午夜精品一区,二区,三区| 99久久精品热视频| 午夜精品久久久久久毛片777| 真实男女啪啪啪动态图| 一进一出抽搐gif免费好疼| 久久中文字幕人妻熟女| 精品99又大又爽又粗少妇毛片 | 动漫黄色视频在线观看| 亚洲真实伦在线观看| 黄频高清免费视频| 成熟少妇高潮喷水视频| 国内精品一区二区在线观看| 国产91精品成人一区二区三区| 天天添夜夜摸| 婷婷丁香在线五月| 久久国产精品影院| 麻豆国产97在线/欧美| 免费看日本二区| 亚洲专区国产一区二区| 十八禁网站免费在线| 桃色一区二区三区在线观看| 黄色日韩在线| 99在线人妻在线中文字幕| 久久这里只有精品19| 99国产精品一区二区三区| 91麻豆精品激情在线观看国产| 熟女人妻精品中文字幕| netflix在线观看网站| 日韩欧美一区二区三区在线观看| 国产成人av教育| 757午夜福利合集在线观看| 亚洲,欧美精品.| 久久这里只有精品19| 成人国产综合亚洲| 欧美黑人欧美精品刺激| 在线免费观看不下载黄p国产 | 国产午夜福利久久久久久| 九九热线精品视视频播放| 99久久综合精品五月天人人| 日本成人三级电影网站| 久久香蕉国产精品| 精品日产1卡2卡| 欧美日韩黄片免| 一个人免费在线观看电影 | 一区福利在线观看| 桃红色精品国产亚洲av| 波多野结衣巨乳人妻| 美女大奶头视频| 亚洲成av人片免费观看| 免费在线观看视频国产中文字幕亚洲| 国产精品久久久久久精品电影| 久9热在线精品视频| 国产伦在线观看视频一区| 久久精品夜夜夜夜夜久久蜜豆| 国产精品美女特级片免费视频播放器 | 国产视频一区二区在线看| 欧美性猛交黑人性爽| www.精华液| 伊人久久大香线蕉亚洲五| 免费人成视频x8x8入口观看| 可以在线观看毛片的网站| 无遮挡黄片免费观看| 18禁裸乳无遮挡免费网站照片| 久久精品国产综合久久久| 国产精品久久久久久精品电影| 好男人在线观看高清免费视频| 天天一区二区日本电影三级| 天堂动漫精品| 老司机午夜福利在线观看视频| 动漫黄色视频在线观看| 精品免费久久久久久久清纯| 真实男女啪啪啪动态图| 国产野战对白在线观看| 中亚洲国语对白在线视频| 国产亚洲精品综合一区在线观看| 免费av毛片视频| 久久这里只有精品19| 日本成人三级电影网站| 黄色视频,在线免费观看| 美女扒开内裤让男人捅视频| 精品国产超薄肉色丝袜足j| 欧美性猛交黑人性爽| 色尼玛亚洲综合影院| 日韩免费av在线播放| 亚洲五月天丁香| 99久久精品一区二区三区| 国产精品av视频在线免费观看| 成年人黄色毛片网站| 亚洲最大成人中文| 深夜精品福利| 一卡2卡三卡四卡精品乱码亚洲| 久久欧美精品欧美久久欧美| 啪啪无遮挡十八禁网站| 人人妻人人看人人澡| 99热只有精品国产| 久久久久性生活片| 国产亚洲精品久久久久久毛片| 亚洲性夜色夜夜综合| 一本精品99久久精品77| 久久久国产欧美日韩av| 好男人在线观看高清免费视频| 国内精品一区二区在线观看| 久久精品综合一区二区三区| www.熟女人妻精品国产| 欧美又色又爽又黄视频| 午夜福利视频1000在线观看| 亚洲一区二区三区不卡视频| 18禁国产床啪视频网站| 亚洲一区高清亚洲精品| 精品一区二区三区视频在线观看免费| av在线天堂中文字幕| 亚洲欧洲精品一区二区精品久久久| 一级a爱片免费观看的视频| 在线观看免费午夜福利视频| 在线a可以看的网站| 国内毛片毛片毛片毛片毛片| 精品久久久久久成人av| 国产高潮美女av| 麻豆av在线久日| 好看av亚洲va欧美ⅴa在| АⅤ资源中文在线天堂| 亚洲精品乱码久久久v下载方式 | 亚洲真实伦在线观看| 久久天躁狠狠躁夜夜2o2o| 黄片小视频在线播放| 色综合站精品国产| 观看免费一级毛片| 五月伊人婷婷丁香| 很黄的视频免费| 99久久久亚洲精品蜜臀av| 桃红色精品国产亚洲av| 国产精品98久久久久久宅男小说| 偷拍熟女少妇极品色| 日日摸夜夜添夜夜添小说| 黄片小视频在线播放| 亚洲男人的天堂狠狠| 国产av在哪里看| 一本一本综合久久| 国产毛片a区久久久久| 国产精品自产拍在线观看55亚洲| 精品国产乱子伦一区二区三区| 午夜免费观看网址| 99热这里只有精品一区 | 丰满的人妻完整版| 特大巨黑吊av在线直播| 色哟哟哟哟哟哟| 亚洲国产欧美一区二区综合| 国产亚洲精品av在线| a级毛片a级免费在线| 窝窝影院91人妻| 亚洲精品色激情综合| 在线播放国产精品三级| 国产av在哪里看| 狂野欧美白嫩少妇大欣赏| 身体一侧抽搐| 俺也久久电影网| 手机成人av网站| 国产伦人伦偷精品视频| 9191精品国产免费久久| 99在线人妻在线中文字幕| 国产aⅴ精品一区二区三区波| 国产毛片a区久久久久| 99热精品在线国产| 亚洲av免费在线观看| www.熟女人妻精品国产| 无遮挡黄片免费观看| 精品一区二区三区四区五区乱码| 麻豆久久精品国产亚洲av| 美女cb高潮喷水在线观看 | 久久精品国产清高在天天线| 久久草成人影院| 99久久综合精品五月天人人| 身体一侧抽搐| 一本一本综合久久| 国产99白浆流出| 久久性视频一级片| 99视频精品全部免费 在线 | 一级毛片高清免费大全| 国产一区二区激情短视频| 最近最新中文字幕大全电影3| 男女做爰动态图高潮gif福利片| 国产精品亚洲av一区麻豆| av天堂在线播放| 亚洲人成伊人成综合网2020| 麻豆久久精品国产亚洲av| 成人无遮挡网站| 搞女人的毛片| 亚洲九九香蕉| 国产精品亚洲美女久久久| 欧美性猛交╳xxx乱大交人| 色综合婷婷激情| 亚洲av免费在线观看| 日韩大尺度精品在线看网址| 久久精品影院6| 99热精品在线国产| 99精品欧美一区二区三区四区| www.精华液| 午夜福利18| 欧美丝袜亚洲另类 | 国产野战对白在线观看| 欧美中文日本在线观看视频| 天天躁狠狠躁夜夜躁狠狠躁| e午夜精品久久久久久久| 国产精品影院久久| 国产日本99.免费观看| 国产av在哪里看| 最新在线观看一区二区三区| 村上凉子中文字幕在线| 亚洲一区高清亚洲精品| 无遮挡黄片免费观看| 精品久久蜜臀av无| 成人三级黄色视频| 久久久久国产一级毛片高清牌| 亚洲乱码一区二区免费版| 国产日本99.免费观看| 男女午夜视频在线观看| 国产精品亚洲av一区麻豆| 五月伊人婷婷丁香| 国产精品久久电影中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 亚洲男人的天堂狠狠| 日本 av在线| 91麻豆av在线| 在线观看一区二区三区| 亚洲国产精品久久男人天堂| 久久精品aⅴ一区二区三区四区| 国产激情久久老熟女| 黄色女人牲交| 国产av一区在线观看免费| 黄片小视频在线播放| 欧美一级a爱片免费观看看| 国产午夜精品久久久久久| 小蜜桃在线观看免费完整版高清| 久久精品综合一区二区三区| 久久久久久久久中文| 99国产综合亚洲精品| 97超视频在线观看视频| 国产高清videossex| 国产成年人精品一区二区| 亚洲 欧美一区二区三区| 麻豆国产97在线/欧美| 亚洲成av人片在线播放无| 成在线人永久免费视频| 欧美绝顶高潮抽搐喷水| 91在线精品国自产拍蜜月 | 欧美激情久久久久久爽电影| 日本 欧美在线| 亚洲中文字幕日韩| av女优亚洲男人天堂 | 麻豆国产av国片精品| 男人舔奶头视频| 久久久国产成人精品二区| 亚洲成人中文字幕在线播放| 亚洲精品乱码久久久v下载方式 | 欧美一区二区国产精品久久精品| 99久久99久久久精品蜜桃| 日本免费a在线| 老熟妇乱子伦视频在线观看| 嫩草影院入口| 亚洲中文字幕一区二区三区有码在线看 | 日韩欧美在线二视频| 18禁国产床啪视频网站| 亚洲 欧美一区二区三区| 成人特级av手机在线观看| 无限看片的www在线观看| 国产高清有码在线观看视频| 女同久久另类99精品国产91| 欧美色视频一区免费| 欧美又色又爽又黄视频| 精品日产1卡2卡| 99久久精品热视频| 高潮久久久久久久久久久不卡| 88av欧美| 色综合婷婷激情| 亚洲av熟女| 亚洲中文字幕日韩| 99久久国产精品久久久| 亚洲国产精品999在线| 亚洲欧洲精品一区二区精品久久久| 国产精品免费一区二区三区在线| 久久国产精品人妻蜜桃| 国产成人福利小说| 国产亚洲精品一区二区www| 国产精品久久久久久人妻精品电影| 国产亚洲av高清不卡| 成年人黄色毛片网站| 91字幕亚洲| 亚洲精品456在线播放app | 国内毛片毛片毛片毛片毛片| 国产精品乱码一区二三区的特点| 亚洲国产欧美一区二区综合| 国产99白浆流出| 国产成人一区二区三区免费视频网站| 午夜福利18| 国产精品永久免费网站| 亚洲精品国产精品久久久不卡| 9191精品国产免费久久| 日日干狠狠操夜夜爽|