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

    顯微散焦粒子識(shí)別與追蹤測(cè)速

    2017-01-13 07:23:51趙攀杰王昊利韓威俊
    關(guān)鍵詞:擋板標(biāo)定粒子

    趙攀杰,王昊利,汪 兵,韓威俊

    (中國(guó)計(jì)量學(xué)院計(jì)量測(cè)試工程學(xué)院,杭州 310018)

    顯微散焦粒子識(shí)別與追蹤測(cè)速

    趙攀杰,王昊利,汪 兵,韓威俊

    (中國(guó)計(jì)量學(xué)院計(jì)量測(cè)試工程學(xué)院,杭州 310018)

    根據(jù)三孔擋板散焦成像的原理,使用高速散焦顯微粒子測(cè)速裝置開(kāi)展三維顯微粒子追蹤測(cè)速研究;在粒子散焦像識(shí)別的過(guò)程中利用了其對(duì)應(yīng)三角模型相似的特點(diǎn),通過(guò)散焦像斑灰度呈二維高斯分布特征計(jì)算其質(zhì)心;在粒子三維位置測(cè)量算法分析的過(guò)程中選擇了多達(dá)29顆的樣品粒子,對(duì)29組不同深度的數(shù)據(jù)采用線性擬合的方式進(jìn)行粒子深度方向的標(biāo)定,并用二元多項(xiàng)式對(duì)870組數(shù)據(jù)擬合獲得了關(guān)于粒子平面位置漂移的補(bǔ)償函數(shù);針對(duì)長(zhǎng)直微通道內(nèi)雷諾數(shù)Re分別為0.05和0.1的流場(chǎng)進(jìn)行了驗(yàn)證性測(cè)量,采用20×0.4顯微物鏡,使用高速CMOS相機(jī)對(duì)直徑為2μm的粒子進(jìn)行追蹤;結(jié)果顯示,計(jì)算出的示蹤粒子速度分布與數(shù)值仿真曲線吻合良好。

    散焦成像;模型識(shí)別;粒子測(cè)速;三維追蹤

    0 引言

    在實(shí)驗(yàn)流體力學(xué)研究中,二維粒子圖像/追蹤測(cè)速系統(tǒng)(PIV/PTV)只能對(duì)流體進(jìn)行分層測(cè)量,獲取兩個(gè)共面的速度分量(2D-2C),此類系統(tǒng)難以針對(duì)復(fù)雜三維流場(chǎng)進(jìn)行測(cè)量。為了突破該限制開(kāi)展復(fù)雜流動(dòng)三維三分量(3D-3C)測(cè)量已成為實(shí)驗(yàn)流體力學(xué)研究人員的主要目標(biāo),其關(guān)鍵在于確定粒子在深度方向的位置和速度。

    目前已經(jīng)出現(xiàn)若干三維測(cè)速技術(shù),包括:體視技術(shù)[1]可以從兩個(gè)視角獲得粒子圖像,并通過(guò)兩個(gè)圖像平面計(jì)算出的位移來(lái)重建速度分量;數(shù)字全息顯微鏡[2]是數(shù)字全息術(shù)在顯微物鏡中應(yīng)用的結(jié)果。通過(guò)顯微物鏡放大后使用圖像傳感器記錄散射光與非散射光干涉圖樣作為全息圖,粒子的空間位置可通過(guò)全息圖像進(jìn)行數(shù)值重建;波前傳感器[3]和散光粒子追蹤測(cè)速[4]利用一個(gè)附加的柱面透鏡取代光學(xué)向場(chǎng)透鏡來(lái)創(chuàng)建一個(gè)變形成像系統(tǒng);利用粒子衍射模式的散焦技術(shù)[5]是基于散焦粒子圖像的衍射效應(yīng)。散焦粒子圖像構(gòu)成的同心環(huán)結(jié)構(gòu),其大小隨深度間隔不同而變化,從而確定粒子的深度位置坐標(biāo)。此外,有一類散焦技術(shù)利用多孔擋板使同一顆粒子在成像平面上呈現(xiàn)多個(gè)像斑從而形成一種特定的圖樣。通過(guò)圖樣大小隨其粒子深度變化的相關(guān)關(guān)系便可計(jì)算出粒子的深度位置。這一散焦的概念由Willert和Gharib[6]首次提出,通過(guò)在圖像傳感器與被測(cè)實(shí)驗(yàn)對(duì)象之間嵌入三孔擋板的方式來(lái)成像。2006年,Yoon和Kim[7]首次將這種方法應(yīng)用在了一個(gè)微尺度后臺(tái)階流動(dòng)中。后來(lái),Pereria等人[8]基于此技術(shù)而設(shè)計(jì)了一個(gè)高速三維顯微散焦粒子測(cè)速系統(tǒng),并提出了一種標(biāo)定方法作為復(fù)雜多元化顯微鏡系統(tǒng)的一個(gè)等效光學(xué)模型來(lái)簡(jiǎn)化計(jì)算。而Tien等人[9-10]在擋板的每個(gè)孔上放置了一個(gè)可過(guò)濾不同單色光的濾色片以便容易識(shí)別散焦圖像,并將其從宏觀擴(kuò)展到微觀流速測(cè)量的運(yùn)用中。

    本文利用已經(jīng)搭建的一套高速散焦顯微粒子圖像測(cè)速裝置,在粒子三維空間位置標(biāo)定的過(guò)程中,通過(guò)挑選更多的粒子來(lái)對(duì)其相關(guān)的測(cè)量算法進(jìn)行更深一步的分析與研究。并在此基礎(chǔ)上,針對(duì)長(zhǎng)直微通道內(nèi)Poiseuille流進(jìn)行了測(cè)量,將測(cè)量結(jié)果與數(shù)值仿真的速度廓線進(jìn)行了比較,驗(yàn)證該實(shí)驗(yàn)系統(tǒng)針對(duì)流體不同深度速度測(cè)量的有效性。

    1 散焦追蹤測(cè)速裝置

    搭建的高速散焦Micro-PTV測(cè)速系統(tǒng)裝置,主要包括德國(guó)LaVison顯微系統(tǒng)、MEMRECAMHX-6型CMOS高速相機(jī)(最大幀頻可達(dá)5000幀/秒)、連續(xù)激光器LIF-532-5(功率為5W)、三孔擋板、注射泵、微流控裝置、計(jì)算機(jī)等。

    在實(shí)驗(yàn)中通過(guò)注射泵控制微流控裝置內(nèi)流體的流速,并使用20×0.4倍的顯微物鏡對(duì)微流場(chǎng)中直徑為2μm的示蹤粒子進(jìn)行觀測(cè)。利用連續(xù)激光器發(fā)射的532nm的激光可使這種示蹤粒子激發(fā)出610nm的熒光,這種圖像信號(hào)經(jīng)過(guò)La Vison顯微光學(xué)系統(tǒng)處理后由CMOS高速相機(jī)進(jìn)行采集,在計(jì)算機(jī)顯示屏中顯示示蹤粒子的圖像。為了能夠獲取示蹤粒子的深度信息,在LaVison顯微系統(tǒng)的顯微物鏡后面嵌入了一個(gè)三孔擋板。因此,示蹤粒子發(fā)生散焦現(xiàn)象并形成3個(gè)光斑。設(shè)計(jì)的圓形三孔擋板上3個(gè)圓孔的直徑均為2 mm,三小孔圓心所構(gòu)成的正三角形的幾何中心位于圓形擋板的中心軸線上,其對(duì)應(yīng)的外接圓直徑為4 mm。其中,三小孔的邊緣厚度均只有0.1 mm,且距擋板邊緣間的厚度是成線性分布的。

    三孔擋板利用了散焦成像的原理,示蹤粒子的熒光信號(hào)通過(guò)三孔擋板成像的光路系統(tǒng)原理圖如圖1所示。示蹤粒子在深度方向上的位置不同,在CMOS高速相機(jī)成像平面所成的散焦像的像徑(即正三角形構(gòu)型的外接圓直徑)大小是不同的。當(dāng)示蹤粒子位于參考平面上時(shí),所成像聚焦成一個(gè)像點(diǎn),可將其像徑的大小視為零。而當(dāng)示蹤粒子分別位于z>0和z<0這兩個(gè)方向且距參考平面等距離處,其對(duì)應(yīng)的散焦圖像的像徑大小相同,但方向相反。因此散焦圖像能夠避免粒子在焦平面兩側(cè)的深度位置產(chǎn)生歧義。

    圖1 孔板散焦成像原理簡(jiǎn)圖[11]

    2 三維粒子位置測(cè)量算法

    2.1 粒子散焦圖像識(shí)別

    在粒子深度的標(biāo)定過(guò)程中,當(dāng)熒光粒子聚焦時(shí)所成的圖像如圖2所示;當(dāng)其發(fā)生散焦時(shí),所成的圖像如圖3所示,每個(gè)熒光粒子對(duì)應(yīng)的散焦圖像均為3個(gè)圓斑。為了識(shí)別出熒光粒子對(duì)應(yīng)的散焦圖像,采用了Yoon等人[7]提出散焦圖像所成三角形構(gòu)型相似的理論,其相似度公式如下:

    其中:向量(r1,r2,r3)為參考三角構(gòu)型(以三孔擋板的3個(gè)小孔的圓心作為其頂點(diǎn))三邊的矢量,向量(c1,c2,c3)為熒光粒子散焦圖像所成三角形構(gòu)型三邊的矢量。在熒光粒子散焦圖像的識(shí)別過(guò)程中,如果S≥0.99,則可把正在檢測(cè)中的3個(gè)圓斑當(dāng)作熒光粒子對(duì)應(yīng)的散焦圖像。

    從圖3中可以看出,每個(gè)三角形頂點(diǎn)處的像斑就是熒光粒子對(duì)應(yīng)的散焦圖像。利用像斑灰度值呈二維高斯分布的特點(diǎn),能方便求出像斑的質(zhì)心。這里用(X1,Y1),(X2,Y2)和(X3,Y3)分別表示熒光粒子散焦圖像3個(gè)圓像斑的質(zhì)心,進(jìn)而可計(jì)算出由3個(gè)質(zhì)心所成三角構(gòu)型外接圓的直徑D和圓心(X,Y)(見(jiàn)圖4)。

    圖2 熒光粒子聚焦成像

    圖3 熒光粒子散焦成像及其識(shí)別

    圖4 熒光粒子散焦成像三角構(gòu)型

    2.2 粒子位置標(biāo)定與分析

    在粒子所成散焦圖像像徑與其深度關(guān)系的標(biāo)定過(guò)程中,被測(cè)粒子由之前實(shí)驗(yàn)所采用的11顆粒子[11]增加到了29顆,并將標(biāo)定深度由z=-40μm增加到了z=-70μm。把粒子群聚焦時(shí)的位置定義為參考平面(即粒子深度z=0處),從z=-14μm處開(kāi)始,每隔2μm就用CMOS相機(jī)記錄一次粒子群的散焦圖像,直到z=-70μm為止。為了減少粒子像斑測(cè)量的誤差,取每個(gè)特定深度處粒子群的200張散焦圖片,以其平均像徑作為此深度處粒子的像徑,然后再通過(guò)線性擬合的方法對(duì)粒子的29對(duì)像徑與深度進(jìn)行處理。從圖5(a)中可以看出,粒子像徑與其對(duì)應(yīng)的深度有比較好的線性關(guān)系。雖然在這次的測(cè)量中粒子的像徑與深度關(guān)系中也呈現(xiàn)出了弱非線性的特征(見(jiàn)圖5(b)),但其斜率ΔD/Δz隨粒子深度在參考線附近(圖中水平虛線)上下浮動(dòng),其均方根誤差(RMS)只有0.071,這與Pereira等[12]的研究結(jié)論是相一致的。

    圖5 粒子像經(jīng)與深度圖

    因此,粒子像徑D與深度z之間的關(guān)系式可表示為:

    微通道內(nèi)流體介質(zhì)不同,粒子所對(duì)應(yīng)的深度標(biāo)定函數(shù)的系數(shù)就會(huì)發(fā)生變化。在減少重復(fù)標(biāo)定工作方面,Yoon等人[13]提出了流體介質(zhì)的折射率與其標(biāo)定系數(shù)之間的關(guān)系式,如下:

    其中:nfluid和nref分別為實(shí)驗(yàn)溶液和參考流體的折射率,(?z/?D)fluid和(?z/?D)ref為這兩種溶液對(duì)應(yīng)的標(biāo)定系數(shù)。

    為了消除固定熒光粒子像平面中的位置(X,Y)隨深度變化產(chǎn)生的漂移,獲得補(bǔ)償后粒子在水平面中的實(shí)際位置(x,y),這里,采用二元多項(xiàng)式補(bǔ)償函數(shù)(見(jiàn)式(7)和式(8))[7],通過(guò)最小二乘法對(duì)粒子群的870對(duì)三維坐標(biāo)數(shù)據(jù)進(jìn)行了擬合,獲得了粒子在水平位置的補(bǔ)償函數(shù)。

    通過(guò)對(duì)任意一個(gè)選定的粒子使用水平位置x和y的補(bǔ)償函數(shù)后,從圖6(a)和(b)中可看到,粒子在水平面中的實(shí)際位置(x,y)獲得了較好的改善,其大小隨粒子深度在參考線(過(guò)粒子在z=0處的水平位置(x,y)的虛線)附近輕微浮動(dòng)。

    圖6

    圖7給出的是粒子三維測(cè)量位置x/y/z的誤差分析結(jié)果。從圖中可以看出,29顆熒光粒子在x/y/z3個(gè)方向上的均方根誤差(root mean square error,即RMSE)分別在參考線Rx=0.1468μm、Ry=0.1859μm和Rz=0.4329μm附近上線變化。鑒于深度標(biāo)定中的弱非線性趨勢(shì)和粒子三維空間分辨率的不同等情況,粒子位置測(cè)量在z上的誤差是在x與y的2~3倍。這與Pereira[12]及Yoon[13]等研究人員測(cè)量的實(shí)驗(yàn)結(jié)果是一致的。

    圖7 粒子測(cè)量位置x/y/z誤差分析

    3 粒子追蹤測(cè)速實(shí)驗(yàn)

    在利用散焦粒子圖像測(cè)速技術(shù)對(duì)長(zhǎng)直微通道的速度場(chǎng)進(jìn)行了測(cè)量,通道的截面為矩形和流向的結(jié)構(gòu)見(jiàn)圖8,其中長(zhǎng)、寬及深度分別為5 mm、2 mm和50μm。由于寬深比達(dá)到了40 ∶1,因此可以近似地認(rèn)為是兩平板間流動(dòng)。通道上端為玻璃板、下端由PDMS材料組成,兩種材料具有良好的粘附性,制作簡(jiǎn)單,并在低速流下具有良好的密封效果。

    被測(cè)流體為熒光粒子原溶液與去離子水按照1∶1 000體積比配制的水溶液,采用注射泵(格蘭公司)進(jìn)行驅(qū)動(dòng),流量分別設(shè)定為185.8μl/h和371.6μl/h,對(duì)應(yīng)的雷諾數(shù)Re分別為0.05和0.1。測(cè)量點(diǎn)位于距長(zhǎng)直微通道進(jìn)口3 mm處附近,流場(chǎng)已經(jīng)達(dá)到充分發(fā)展階段,理論速度廓線為平板Poiseuille流。

    圖8 微通道矩形及流向結(jié)構(gòu)截面視圖

    使用CMOS高速相機(jī)以1 000幀/秒的掃描速率對(duì)微通道內(nèi)得流場(chǎng)進(jìn)行拍攝,對(duì)其記錄下來(lái)的序列圖像通過(guò)維納濾波去除噪聲后,可得到較高質(zhì)量的散焦圖像。根據(jù)Yoon等人[7]提出散焦圖像所成三角形構(gòu)型相似的理論公式,能夠識(shí)別出粒子對(duì)應(yīng)的散焦圖像。利用像點(diǎn)追蹤開(kāi)源軟件Vide Spot Tracker v07.02對(duì)粒子3個(gè)散焦像斑進(jìn)行追蹤,根據(jù)像斑灰度值呈二維高斯分布的特點(diǎn),求出像斑的質(zhì)心,進(jìn)而獲得追蹤粒子三像斑所構(gòu)三角形外接圓的直徑與圓心,即追蹤粒子的散焦像徑與質(zhì)心。根據(jù)深度標(biāo)定函數(shù)計(jì)算出粒子的實(shí)際深度位置,通過(guò)粒子水平方向上的補(bǔ)償函數(shù)去對(duì)其質(zhì)心位置進(jìn)行修正,最終根據(jù)所追蹤的示蹤粒子軌跡獲取粒子的速度矢量。

    圖9和圖10分別給出了雷諾數(shù)分別為Re=0.05和Re=0.1的情況下,通道內(nèi)部不同深度處的粒子軌跡。不同顏色代表了不同的深度,而不同深度處粒子在單位時(shí)間內(nèi)移動(dòng)的距離是不同的,且粒子距通道中心越近,其對(duì)應(yīng)的速度越大。利用粒子軌跡分別計(jì)算出兩個(gè)雷諾數(shù)下流場(chǎng)的速度廓線v1和v2,并分別與平板Poiseuille流的數(shù)值仿真結(jié)果進(jìn)行比較,見(jiàn)圖11和圖12??梢钥吹綄?shí)驗(yàn)值和數(shù)值仿真結(jié)果吻合良好,表明實(shí)驗(yàn)值滿足平板Poiseuille流動(dòng)。該實(shí)驗(yàn)結(jié)果驗(yàn)證了高速散焦Micro-PTV測(cè)速系統(tǒng)能夠有效測(cè)量不同深度的流動(dòng)速度,因此能夠進(jìn)一步推廣至三維復(fù)雜流場(chǎng)的測(cè)量中。

    圖9 Re=0.05時(shí)長(zhǎng)直微通道內(nèi)的粒子軌跡

    圖10 Re=0.1時(shí)長(zhǎng)直微通道內(nèi)的粒子軌跡

    圖11 Re=0.05時(shí)粒子速度測(cè)量值與仿真值的比較

    4 結(jié)束語(yǔ)

    在已搭建的一套高速散焦顯微粒子測(cè)速裝置的基礎(chǔ)上,通過(guò)增加粒子數(shù)并增大深度范圍對(duì)粒子深度進(jìn)行了標(biāo)定,并用二元多項(xiàng)式對(duì)870組數(shù)據(jù)擬合計(jì)算出了關(guān)于粒子平面位置漂移的修正函數(shù),完成了粒子三維空間位置測(cè)量算法分析與研究。針對(duì)長(zhǎng)直微通道內(nèi)流速的驗(yàn)證性測(cè)量中,分別對(duì)雷諾數(shù)Re=0.05和Re=0.1這兩種不同的流速進(jìn)行了測(cè)量,獲取的示蹤粒子速度分布與數(shù)值仿真曲線吻合良好。實(shí)驗(yàn)結(jié)果進(jìn)一步驗(yàn)證了該系統(tǒng)裝置在三維測(cè)速方面的可行性和追蹤測(cè)速算法的有效性。

    圖12 Re=0.1時(shí)粒子速度測(cè)量值與仿真值的比較

    [1]Lindken R,Westerweel J,Wieneke B,et al.Stereoscopic microparticle image velocimetry[J].Exp Fluids,2006,41(2):161-171.

    [2]Sheng J,Malkiel E,Katz J,et al.Digital holographic microscope for measuring three-dimensional particle distributions and motions [J].Appl.Opt,2006,45(16):3893-3901.

    [3]Angarita-Jaimes NC,McGhee E,Chennaoui M,et al.Wave front sensing for single view three-component three-dimensional flow velocimetry[J].Exp Fluids.2006,41(6):881-891.

    [4]Cierpka C,Segura R,Hain R,et al.A simple single camera 3C3D velocity measurement technique without errors due to depth of correlation and spatial averaging for micro fluidics[J].Meas.Sci.Technol.,2010,21(4):045401.

    [5]Park J S,Kihm K D.Three-dimensional micro-PTV using deconvolution microscopy[J].Exp Fluids.2006,40(3):491-499.

    [6]Willert C,Gharib M.Three-dimensional particle imaging with a single camera[J].Exp Fluids.1992,12:353-358.

    [7]Yoon S Y,KimKC.3D particle position and 3D velocity field measurement in amicrovolume via the defocusing concept[J].Meas Sci Technol.,2006,17(11):2897-2905.

    [8]Pereira F,Lu J,Castano-Graff E,et al.Microscale 3D flow mapping withμDDPIV[J].Exp.Fluids,2007,42(4):589-599.

    [9]Tien WH,Kartes P,Yamasaki T,et al.A color coded back lighted defocusing digital particle image velocimetry system[J].Exp Fluids,2008,44:1015-1026.

    [10]Wei H,Tien D,Dabiri J,et al.Color-coded three-dimensional micro partilce tracking velocimetry and appl ication to micro backwardfacing step flows[J].Exp.Fluids,2014,55:1684-1697.

    [11]趙攀杰,王昊利.基于散焦的三維顯微粒子追蹤測(cè)速[J].中國(guó)計(jì)量學(xué)院學(xué)報(bào),2015,26(04):428-434.

    [12]Pereira F,Gharib M.Defocusing digital particle image velocimetry and three-dimensional characterization of two phase flows[J].Meas Sci Technol.2002,13:683-694.

    [13]Yoon S Y,Khim K D,Kim K C.Correlation of fluid refractive index with calibration coefficient for micro-defocusing digital particle image velocimetry[J].Meas.Sci.Technology,2011,22:1-5.

    Identification and Tracking About Three-dimensional Microscopic Defocusing Particle

    Zhao Panjie,Wang Haoli,Wang Bing,Han Weijun

    (College of Metrological Technology and Engineering,China Jiliang University,Hangzhou 310018,China)

    According to the principle of a three-pinhole mask’s defocused imaging,a set of high-speed defocusing micro-PTV devices was employed for further research on the three-dimensional micro particle tracking velocimetry.The similar characteristics among particles’corresponding triangle models was used in the process of their defocusing images recognition,and the centroid of particle defocusing image spot was calculated by its gray in two-dimensional Gaussian distribution.As many as 29 sample particles were selected in the process of measurement algorithm analysis about particle’s three-dimensional position,and 29 different depth data were utilized for calibrating in the particles’depth direction with the method of linear fitting,and the drift compensation functions about the in-plane position of particle was acquired by quadratic polynomial fitting with 870 groups of data.A microscope objective lens of 20×0.4 and high speed CMOS camera were adopted to trace the particles of 2μm diameter in the verification measurement of a long straight micro-channel flow field whose Reynolds number were 0.05 and 0.1 respectively.The results showed that the particle velocity distribution calculated was in good agreement with the numerical simulation curve.

    defocusing image;pattern recognition;particle velocity measurement;three-dimensional tracking

    1671-4598(2016)08-0298-03

    10.16526/j.cnki.11-4762/tp.2016.08.081

    :TP317.4;0353.5

    :A

    2016-03-09;

    :2016-03-24。

    國(guó)家自然科學(xué)基金資助項(xiàng)目(11472261;11172287)。

    趙攀杰(1988-),男,河南許昌人,碩士研究生,主要從事顯微顆粒的可視化測(cè)量方向的研究。

    猜你喜歡
    擋板標(biāo)定粒子
    平拋運(yùn)動(dòng)中的“擋板”問(wèn)題探究
    燃燒器二次風(fēng)擋板開(kāi)度對(duì)爐內(nèi)燃燒特性的影響
    使用朗仁H6 Pro標(biāo)定北汽紳寶轉(zhuǎn)向角傳感器
    基于粒子群優(yōu)化的橋式起重機(jī)模糊PID控制
    基于粒子群優(yōu)化極點(diǎn)配置的空燃比輸出反饋控制
    基于勻速率26位置法的iIMU-FSAS光纖陀螺儀標(biāo)定
    船載高精度星敏感器安裝角的標(biāo)定
    基于Harris-張正友平面標(biāo)定法的攝像機(jī)標(biāo)定算法
    基于Matlab的α粒子的散射實(shí)驗(yàn)?zāi)M
    物理與工程(2014年4期)2014-02-27 11:23:08
    基于兩粒子糾纏態(tài)隱形傳送四粒子GHZ態(tài)
    亚洲专区中文字幕在线| 女人被狂操c到高潮| 久久久国产成人免费| 日韩欧美精品免费久久 | 熟女人妻精品中文字幕| 日本a在线网址| 琪琪午夜伦伦电影理论片6080| 欧美黑人欧美精品刺激| av欧美777| 色综合亚洲欧美另类图片| 少妇丰满av| 少妇被粗大猛烈的视频| 欧美在线黄色| 麻豆久久精品国产亚洲av| 身体一侧抽搐| 女同久久另类99精品国产91| 国产伦一二天堂av在线观看| 中文字幕高清在线视频| 日韩大尺度精品在线看网址| 精品久久久久久久久久久久久| 又爽又黄a免费视频| 一本精品99久久精品77| 又爽又黄a免费视频| 久久精品国产亚洲av天美| 久久精品国产亚洲av天美| 在线观看66精品国产| 午夜免费男女啪啪视频观看 | 自拍偷自拍亚洲精品老妇| 两个人的视频大全免费| 动漫黄色视频在线观看| 12—13女人毛片做爰片一| 亚洲成人久久性| 欧美又色又爽又黄视频| 久久热精品热| 亚洲一区二区三区色噜噜| 久久香蕉精品热| 国产午夜福利久久久久久| 国产午夜福利久久久久久| 色哟哟哟哟哟哟| 亚洲中文字幕一区二区三区有码在线看| 久久香蕉精品热| 成人av一区二区三区在线看| 国产免费一级a男人的天堂| 不卡一级毛片| 97碰自拍视频| 欧美高清性xxxxhd video| 亚洲一区二区三区不卡视频| 极品教师在线视频| www.熟女人妻精品国产| 又爽又黄a免费视频| 精品一区二区三区av网在线观看| 国内精品久久久久精免费| 亚洲av中文字字幕乱码综合| 99热精品在线国产| 午夜日韩欧美国产| 在线播放无遮挡| 搡老熟女国产l中国老女人| 国产激情偷乱视频一区二区| 小蜜桃在线观看免费完整版高清| 精品99又大又爽又粗少妇毛片 | 97超视频在线观看视频| 免费高清视频大片| 国产一区二区亚洲精品在线观看| 国产成人福利小说| 九色成人免费人妻av| 琪琪午夜伦伦电影理论片6080| 国产国拍精品亚洲av在线观看| 五月玫瑰六月丁香| 在线免费观看的www视频| 精品午夜福利在线看| 日韩欧美三级三区| 深夜精品福利| 亚洲狠狠婷婷综合久久图片| 久久久精品欧美日韩精品| 偷拍熟女少妇极品色| 欧美日韩乱码在线| 日韩av在线大香蕉| 亚洲精品久久国产高清桃花| 亚洲最大成人手机在线| 日韩 亚洲 欧美在线| 麻豆国产av国片精品| 在线免费观看不下载黄p国产 | 色哟哟·www| 欧美日韩国产亚洲二区| 中文字幕高清在线视频| 日日摸夜夜添夜夜添av毛片 | 欧美成人免费av一区二区三区| 蜜桃亚洲精品一区二区三区| 蜜桃久久精品国产亚洲av| 级片在线观看| 日日摸夜夜添夜夜添av毛片 | 精品国内亚洲2022精品成人| 国产av麻豆久久久久久久| 国产欧美日韩精品一区二区| 特级一级黄色大片| 小说图片视频综合网站| 1000部很黄的大片| 亚洲真实伦在线观看| 亚洲自拍偷在线| 久久国产乱子免费精品| 日韩欧美一区二区三区在线观看| 九九久久精品国产亚洲av麻豆| 热99在线观看视频| 我的女老师完整版在线观看| 久久久久久久久久黄片| 国产三级黄色录像| 69人妻影院| 久久午夜福利片| 国产真实乱freesex| 99热6这里只有精品| 97超级碰碰碰精品色视频在线观看| 性欧美人与动物交配| 亚洲精品粉嫩美女一区| 九九热线精品视视频播放| 久久午夜福利片| 成熟少妇高潮喷水视频| 看片在线看免费视频| 久久国产精品影院| 日本与韩国留学比较| 亚洲熟妇熟女久久| 嫩草影院入口| 性欧美人与动物交配| 真人做人爱边吃奶动态| 99视频精品全部免费 在线| 久久精品国产自在天天线| 噜噜噜噜噜久久久久久91| 麻豆国产av国片精品| 久久久久精品国产欧美久久久| 老司机深夜福利视频在线观看| 久久热精品热| 俄罗斯特黄特色一大片| 精品人妻一区二区三区麻豆 | av在线老鸭窝| 亚洲专区中文字幕在线| 亚洲精华国产精华精| 高潮久久久久久久久久久不卡| 好男人电影高清在线观看| 亚洲av一区综合| 国产精品久久电影中文字幕| 日本免费一区二区三区高清不卡| 老鸭窝网址在线观看| 国产成人a区在线观看| 精品人妻偷拍中文字幕| 国内毛片毛片毛片毛片毛片| 亚洲在线观看片| av天堂在线播放| 成人午夜高清在线视频| 亚洲在线自拍视频| 亚洲成人精品中文字幕电影| 亚洲欧美清纯卡通| 欧美激情国产日韩精品一区| 88av欧美| 国产精品伦人一区二区| 麻豆久久精品国产亚洲av| 深爱激情五月婷婷| 99热这里只有是精品50| 精品人妻视频免费看| 欧美成人a在线观看| 啦啦啦观看免费观看视频高清| 在线免费观看的www视频| 啪啪无遮挡十八禁网站| 欧美潮喷喷水| 久久99热这里只有精品18| 久久久色成人| 特级一级黄色大片| 高清日韩中文字幕在线| 性欧美人与动物交配| 女同久久另类99精品国产91| 成人特级av手机在线观看| 久久草成人影院| 成人特级黄色片久久久久久久| 久久久久久大精品| 日韩欧美国产在线观看| 久久国产乱子伦精品免费另类| 日本 欧美在线| 脱女人内裤的视频| av欧美777| av女优亚洲男人天堂| 伊人久久精品亚洲午夜| 国产成人福利小说| 精品人妻一区二区三区麻豆 | 免费av观看视频| 午夜a级毛片| 国产真实乱freesex| 欧美+日韩+精品| 国产精品,欧美在线| 午夜激情欧美在线| 亚洲 欧美 日韩 在线 免费| 久久6这里有精品| 亚洲欧美日韩高清在线视频| 久久国产乱子免费精品| 国产精品嫩草影院av在线观看 | 欧美日韩综合久久久久久 | 成人午夜高清在线视频| 日韩欧美在线乱码| 99热只有精品国产| 欧美极品一区二区三区四区| 成人国产一区最新在线观看| 桃红色精品国产亚洲av| 熟妇人妻久久中文字幕3abv| 精品久久国产蜜桃| 色尼玛亚洲综合影院| 又粗又爽又猛毛片免费看| 嫩草影院入口| 午夜免费男女啪啪视频观看 | 亚洲在线观看片| 久久精品国产亚洲av天美| 如何舔出高潮| 黄色配什么色好看| av在线观看视频网站免费| 在线观看美女被高潮喷水网站 | 蜜桃亚洲精品一区二区三区| 淫妇啪啪啪对白视频| 久久热精品热| 嫩草影院新地址| 欧美黄色淫秽网站| netflix在线观看网站| 精品久久久久久,| 搡老妇女老女人老熟妇| 亚洲av不卡在线观看| 欧美不卡视频在线免费观看| 波多野结衣高清无吗| 美女xxoo啪啪120秒动态图 | 精品人妻视频免费看| 亚洲av第一区精品v没综合| 简卡轻食公司| 欧美最新免费一区二区三区 | 国产探花在线观看一区二区| 一区福利在线观看| 最新在线观看一区二区三区| 又粗又爽又猛毛片免费看| 97热精品久久久久久| 亚洲午夜理论影院| 美女黄网站色视频| 国产久久久一区二区三区| 亚洲精品亚洲一区二区| 嫩草影院入口| 内射极品少妇av片p| 国产欧美日韩精品亚洲av| 校园春色视频在线观看| 最好的美女福利视频网| 国产亚洲av嫩草精品影院| 国产探花极品一区二区| 国产三级黄色录像| 老女人水多毛片| 无遮挡黄片免费观看| 久久午夜亚洲精品久久| 国语自产精品视频在线第100页| 欧美一区二区国产精品久久精品| 啦啦啦观看免费观看视频高清| 一区福利在线观看| 特大巨黑吊av在线直播| 精品久久久久久,| 亚洲国产色片| 变态另类成人亚洲欧美熟女| 色播亚洲综合网| 国产精品不卡视频一区二区 | 国产激情偷乱视频一区二区| 欧美精品国产亚洲| 国产高清三级在线| 一a级毛片在线观看| 狂野欧美白嫩少妇大欣赏| 国产亚洲精品久久久久久毛片| 免费大片18禁| 人人妻人人看人人澡| 免费电影在线观看免费观看| 欧美高清性xxxxhd video| 免费看a级黄色片| 麻豆国产av国片精品| av福利片在线观看| 欧美黑人巨大hd| 欧美+日韩+精品| 琪琪午夜伦伦电影理论片6080| 中文字幕人妻熟人妻熟丝袜美| 亚洲成人中文字幕在线播放| 国产成人欧美在线观看| 日本五十路高清| 午夜精品在线福利| 中文亚洲av片在线观看爽| 精品免费久久久久久久清纯| 内射极品少妇av片p| a级毛片免费高清观看在线播放| 国产久久久一区二区三区| 国产午夜福利久久久久久| 久久精品久久久久久噜噜老黄 | 在线免费观看不下载黄p国产 | 国语自产精品视频在线第100页| 90打野战视频偷拍视频| 日本在线视频免费播放| 免费在线观看影片大全网站| 搡老妇女老女人老熟妇| 亚洲在线自拍视频| 亚洲精品日韩av片在线观看| 好男人在线观看高清免费视频| 亚洲,欧美,日韩| 五月玫瑰六月丁香| 久久久国产成人精品二区| 看十八女毛片水多多多| 国产亚洲欧美在线一区二区| 又爽又黄无遮挡网站| 国产白丝娇喘喷水9色精品| 久久中文看片网| 欧美黄色淫秽网站| 午夜福利在线观看免费完整高清在 | 成人一区二区视频在线观看| 精品人妻1区二区| 欧美性猛交黑人性爽| 99久久精品国产亚洲精品| 白带黄色成豆腐渣| 亚洲精品粉嫩美女一区| 狂野欧美白嫩少妇大欣赏| 赤兔流量卡办理| 久久精品国产亚洲av香蕉五月| 51国产日韩欧美| h日本视频在线播放| 欧美bdsm另类| 国产亚洲精品久久久com| 91九色精品人成在线观看| avwww免费| 国产精品久久久久久人妻精品电影| 国产高清三级在线| 亚洲人成网站在线播放欧美日韩| av国产免费在线观看| 亚洲av熟女| 人妻丰满熟妇av一区二区三区| 亚洲精品成人久久久久久| 欧美乱妇无乱码| 伦理电影大哥的女人| 欧美一区二区亚洲| 免费在线观看日本一区| 国产一区二区在线观看日韩| 少妇人妻精品综合一区二区 | 91久久精品电影网| 夜夜爽天天搞| 国产亚洲精品综合一区在线观看| 99热只有精品国产| 成熟少妇高潮喷水视频| 国产精品久久久久久亚洲av鲁大| 美女cb高潮喷水在线观看| 亚洲av五月六月丁香网| 深爱激情五月婷婷| 国产精品野战在线观看| 久久精品国产清高在天天线| 免费看日本二区| 久久久久久九九精品二区国产| 男女视频在线观看网站免费| 日韩欧美在线乱码| 天堂√8在线中文| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 嫩草影院新地址| 亚洲午夜理论影院| 国产成人福利小说| 亚洲第一电影网av| 亚洲三级黄色毛片| 免费看a级黄色片| 青草久久国产| 99国产精品一区二区三区| 久久精品影院6| 特大巨黑吊av在线直播| 国产高清视频在线播放一区| 最近最新免费中文字幕在线| 国产亚洲精品av在线| 深夜精品福利| 亚洲久久久久久中文字幕| 此物有八面人人有两片| 俺也久久电影网| 特级一级黄色大片| 变态另类丝袜制服| 欧美高清性xxxxhd video| 国产精品久久久久久亚洲av鲁大| 国产亚洲av嫩草精品影院| 少妇高潮的动态图| 嫩草影院精品99| 精品久久久久久久久久免费视频| 亚洲成人中文字幕在线播放| 欧美xxxx性猛交bbbb| 国产视频一区二区在线看| 偷拍熟女少妇极品色| 99精品在免费线老司机午夜| 国产伦精品一区二区三区视频9| 97超级碰碰碰精品色视频在线观看| 婷婷亚洲欧美| 变态另类丝袜制服| 精品人妻一区二区三区麻豆 | 一本久久中文字幕| 波多野结衣高清无吗| 一边摸一边抽搐一进一小说| 91狼人影院| 一二三四社区在线视频社区8| 校园春色视频在线观看| 大型黄色视频在线免费观看| 日韩高清综合在线| 亚洲,欧美,日韩| 九色成人免费人妻av| 国产在视频线在精品| 亚洲国产精品sss在线观看| 最近在线观看免费完整版| 亚洲av五月六月丁香网| 久久婷婷人人爽人人干人人爱| 欧美激情在线99| 日本一二三区视频观看| 亚洲av第一区精品v没综合| 亚洲av中文字字幕乱码综合| 国产精品,欧美在线| 亚洲精品乱码久久久v下载方式| 国产真实伦视频高清在线观看 | 99riav亚洲国产免费| 日韩国内少妇激情av| 欧美日韩中文字幕国产精品一区二区三区| 日本一本二区三区精品| 色噜噜av男人的天堂激情| 国内精品久久久久久久电影| 亚洲精品一卡2卡三卡4卡5卡| 乱人视频在线观看| 免费人成在线观看视频色| 国产欧美日韩一区二区三| 日本黄色片子视频| 国产一级毛片七仙女欲春2| 欧美成人a在线观看| 韩国av一区二区三区四区| 国产国拍精品亚洲av在线观看| 老司机午夜福利在线观看视频| 2021天堂中文幕一二区在线观| 免费看日本二区| 国产精品av视频在线免费观看| 黄色丝袜av网址大全| 国产私拍福利视频在线观看| 99热精品在线国产| 亚洲狠狠婷婷综合久久图片| 国产精品一区二区三区四区免费观看 | 国产av在哪里看| 欧美色视频一区免费| 午夜免费男女啪啪视频观看 | 免费av毛片视频| 成人av在线播放网站| 国产精品久久久久久久电影| 狠狠狠狠99中文字幕| 最近最新中文字幕大全电影3| 久久亚洲精品不卡| 精品日产1卡2卡| 美女免费视频网站| 免费看日本二区| 日本五十路高清| 午夜激情福利司机影院| 我的老师免费观看完整版| 亚洲av免费高清在线观看| 日本熟妇午夜| 在线免费观看的www视频| 亚洲第一欧美日韩一区二区三区| 国内精品美女久久久久久| 国产三级黄色录像| 欧美高清性xxxxhd video| 午夜免费成人在线视频| 成人性生交大片免费视频hd| 欧美日韩黄片免| 日本 欧美在线| 亚洲欧美日韩东京热| 久久精品综合一区二区三区| 国产真实伦视频高清在线观看 | 国产真实伦视频高清在线观看 | 国产美女午夜福利| 亚洲精品一卡2卡三卡4卡5卡| 欧美一区二区亚洲| 少妇的逼水好多| 国产黄a三级三级三级人| 十八禁国产超污无遮挡网站| 热99在线观看视频| 亚洲欧美精品综合久久99| 亚洲乱码一区二区免费版| 中文字幕精品亚洲无线码一区| 日韩欧美 国产精品| 久久午夜亚洲精品久久| 小说图片视频综合网站| 久久久久久久精品吃奶| 亚洲va日本ⅴa欧美va伊人久久| 舔av片在线| 激情在线观看视频在线高清| 国产亚洲精品久久久com| 99国产极品粉嫩在线观看| 人妻久久中文字幕网| 亚洲av美国av| 9191精品国产免费久久| 嫩草影院入口| 91字幕亚洲| 老司机福利观看| 男女床上黄色一级片免费看| 成人永久免费在线观看视频| 久久久久九九精品影院| 婷婷精品国产亚洲av在线| 亚洲男人的天堂狠狠| 一二三四社区在线视频社区8| 免费观看精品视频网站| 波多野结衣高清作品| 久久久国产成人精品二区| 国内精品久久久久精免费| 亚洲精品色激情综合| 亚洲三级黄色毛片| 成人特级av手机在线观看| 99久久精品一区二区三区| 久久精品综合一区二区三区| 午夜免费男女啪啪视频观看 | 精品久久久久久久人妻蜜臀av| 国产成年人精品一区二区| 国产高清三级在线| 美女黄网站色视频| 国产aⅴ精品一区二区三区波| 国产免费av片在线观看野外av| 亚洲第一电影网av| 最新中文字幕久久久久| 亚洲avbb在线观看| 久久亚洲精品不卡| 亚洲美女搞黄在线观看 | 国产亚洲精品综合一区在线观看| 欧美性猛交黑人性爽| 国产精品野战在线观看| 国产一区二区亚洲精品在线观看| 国产欧美日韩精品一区二区| 亚洲精品粉嫩美女一区| av专区在线播放| a级毛片a级免费在线| 中文字幕av在线有码专区| 男女视频在线观看网站免费| 好男人电影高清在线观看| 国产成人啪精品午夜网站| 精品久久国产蜜桃| 国语自产精品视频在线第100页| 亚洲成人久久爱视频| 一区福利在线观看| 少妇的逼水好多| 色精品久久人妻99蜜桃| av天堂在线播放| 99久久成人亚洲精品观看| 噜噜噜噜噜久久久久久91| 小说图片视频综合网站| 亚洲成a人片在线一区二区| 黄色女人牲交| 国产黄片美女视频| 国产美女午夜福利| 国产精品av视频在线免费观看| 午夜视频国产福利| 嫩草影院入口| 人人妻,人人澡人人爽秒播| 天堂影院成人在线观看| av在线天堂中文字幕| 精品不卡国产一区二区三区| 久久久久性生活片| 高清日韩中文字幕在线| 成人永久免费在线观看视频| 国产精品嫩草影院av在线观看 | 亚洲自偷自拍三级| 少妇熟女aⅴ在线视频| 一级作爱视频免费观看| 男女下面进入的视频免费午夜| 天堂av国产一区二区熟女人妻| 91久久精品国产一区二区成人| 91字幕亚洲| 在线播放国产精品三级| 嫩草影院精品99| 成年女人看的毛片在线观看| 一级av片app| 亚洲av成人av| 国产乱人伦免费视频| 精品无人区乱码1区二区| 色吧在线观看| 欧美激情国产日韩精品一区| av女优亚洲男人天堂| 亚洲美女视频黄频| 国产熟女xx| 毛片女人毛片| 色噜噜av男人的天堂激情| 赤兔流量卡办理| 国产成人影院久久av| 丁香欧美五月| 国产在线男女| 中国美女看黄片| 国产真实伦视频高清在线观看 | 亚洲av成人不卡在线观看播放网| 又爽又黄a免费视频| 欧美日韩乱码在线| 一进一出抽搐gif免费好疼| 色尼玛亚洲综合影院| 国产精品精品国产色婷婷| 日本黄色视频三级网站网址| 人妻久久中文字幕网| 欧美性猛交黑人性爽| 一区二区三区激情视频| 中出人妻视频一区二区| a级毛片免费高清观看在线播放| 午夜两性在线视频| 亚洲精品456在线播放app | 美女 人体艺术 gogo| 亚洲av.av天堂| 亚洲欧美日韩东京热| 午夜激情福利司机影院| 美女高潮喷水抽搐中文字幕| 成人高潮视频无遮挡免费网站| 日韩有码中文字幕| 免费av毛片视频| 亚洲va日本ⅴa欧美va伊人久久| 日韩 亚洲 欧美在线| 亚洲,欧美精品.| 人妻制服诱惑在线中文字幕| 99热6这里只有精品| 国产高清视频在线观看网站| h日本视频在线播放| 成人鲁丝片一二三区免费| 国产精品三级大全| 国产精品久久久久久久电影| 一卡2卡三卡四卡精品乱码亚洲| 噜噜噜噜噜久久久久久91| av在线蜜桃| av黄色大香蕉| 亚洲欧美日韩卡通动漫| 免费看a级黄色片| av专区在线播放|