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

    Huang族校正電容層析成像圖像重建算法

    2018-12-21 09:37:04陳宇李洪宇
    哈爾濱理工大學學報 2018年5期
    關(guān)鍵詞:層析成像流型校正

    陳宇 李洪宇

    摘要:針對電容層析成像(ECT)技術(shù)中的“軟場”效應(yīng)和病態(tài)問題,提出了一種Huang族校正的電容層析成像圖像重建算法。首先依據(jù)ECT系統(tǒng)的基本原理,推導(dǎo)出ECT問題中Huang族校正的校正公式,其次給出校正后用于ECT反問題求解仿真實驗的迭代公式。最后,采用數(shù)字仿真模擬實驗方式,驗證提出方法的有效性。實驗結(jié)果表明,Huang族校正方法對于極低位、低位、核心流而言,圖像誤差分別降到2439%、2581%和4091%,均低于LBP、Landweber、SD和CG方法;對于極低位、低位及柱狀流而言,迭代次數(shù)分別為12、12、27次,比Landweber算法和SD法都要低,綜合分析,可知Huang族校正方法實驗效果良好。

    關(guān)鍵詞:

    電容層析成像; Huang族校正;圖像重建

    DOI:10.15938/j.jhust.2018.05.014

    中圖分類號: TN 91173

    文獻標志碼: A

    文章編號: 1007-2683(2018)05-0080-06

    Abstract:To solve the ‘softfield nature and the illposed problem in electrical capacitance tomography technology, a Huang clan correction image reconstruction algorithm for electrical capacitance tomography is presented Firstly,according to the basic principles of the Electrical Capacitance Tomography system, the formula of Huang clan correction in the problem of capacitance tomography is derived Secondly, the iterative formula for simulation experiment is given after the correctionFinally, the validity of the proposed method is verified by digital simulationThe simulation experiment results show that the error of the image for extremely low layer flow,low layer flow and core flow dropped to 2439%,2581% and 4091% respectively Results were lower than Linear Back Projection method,Landweber method,Steepest Descent method and Conjugate Gradient method In addition,the number of iterations were maintained at 12,12 and 27times,also less than the Landweber method and the Steepest Descent method The results of the analysis show that the effect of the Huang Clan Correction Image Reconstruction Algorithm are good

    Keywords:electrical capacitance tomography; Huang clan correction; image reconstruction

    0引言

    過程層析成像技術(shù)簡稱為(process tomography,PT)技術(shù)。其包括電容層析成像技術(shù)(electrical capacitance tomography, ECT)、電阻層析成像技術(shù)(electrical resistance tomography,ERT)和電磁層析成像技術(shù)(electromagnetic tomography,EMT)。他們分別依據(jù)測量不同的電學特性,來進行相關(guān)信息的重建工作。其中ECT技術(shù)就是利用傳感器測量得到的電容信號,依照某種算法重建出介質(zhì)內(nèi)部的介電常數(shù)分布情況。

    ECT技術(shù)作為PT技術(shù)的一種,相比其他PT技術(shù)而言, 擁有諸多優(yōu)點。如成本低、可靠性高、簡單易實現(xiàn)等[1-2]。另外ECT的非侵入性、高安全性和高速度性使其在解決連續(xù)相非導(dǎo)電物質(zhì)的兩相流、多相流成像問題中產(chǎn)生很好的效果。ECT技術(shù)的適用范圍廣泛,現(xiàn)已在化工、石油和醫(yī)藥等眾多兩相流和多相流檢測領(lǐng)域有相應(yīng)的應(yīng)用[3]。

    自從ECT技術(shù)出現(xiàn)以后,大量的圖像重建算法被國內(nèi)外的眾學者所提出。Isaksen O在1996年就對ECT的各種圖像重建算法進行過一次綜合表述[4]。近些年,Tikhonov正則化方法[5]、正則化的GaussNewton方法[6]、快速模擬退火全局優(yōu)化算法[7]、遺傳算法[8]、支持向量機法等一系列新的ECT圖像重建算法又被提出。

    目前,對于ECT圖像重建最常用幾種方法[9]有直接算法中的LBP(linear back projection)法[10]和Tikhonov正則化法,迭代算法中的共軛梯度(conjugate gradient ,CG)法[11]和Landweber迭代法[12]。

    ECT圖像重建是一個高度非線性問題,因ECT技術(shù)具有的“軟場”特性,該問題屬于病態(tài)問題求解。采用針對與該非線性問題“鄰近”的線性問題的解法不斷逼近要得到的解的這樣一種思想是簡單而有效的。

    ECT技術(shù)是否可以成功應(yīng)用關(guān)鍵取決于ECT反問題求解的精度與速度的優(yōu)劣[13]。LBP算法的成像速度非??欤涑上窬壬蠀s有很大的不足[14]。CG算法在復(fù)雜流型上的成像精度并不特別理想。模擬退火算法、遺傳算法和支持向量機法在某方面的優(yōu)越性上有所強化,但卻以犧牲其它性能作為代價,因此只能在特定領(lǐng)域中有所應(yīng)用。

    電容層析成像技術(shù)要實現(xiàn)滿足工業(yè)實際應(yīng)用需求的目標,緊緊地依賴于ECT圖像重建算法的效果,所以尋找優(yōu)秀的ECT圖像重建算法是值得一直研究的目標。

    本文提出了一種基于Huang族校正的電容層析成像圖像重建算法,目的是獲得更好的電容層析成像圖像重建效果。仿真實驗結(jié)果證明該算法可以達到研究目的。在與采用LBP、經(jīng)典Landweber、CG以及SD(steepest descent)方法的仿真實驗的重建圖像、誤差率和迭代次數(shù)進行對比后,可以發(fā)現(xiàn)本文提出的基于Huang族校正的電容層析成像圖像重建方法顯示出了其自身精度和速度上的優(yōu)勢。

    1電容層析成像系統(tǒng)及ECT技術(shù)原理概要

    在ECT研究中,常用的有6、8、12、16電極的ECT系統(tǒng)。這里來討論典型12電極的ECT系統(tǒng)。12電極ECT系統(tǒng)主要組成部分有:電容傳感器、測量與數(shù)據(jù)采集系統(tǒng)和成像計算機[15-18],如下圖1所示。

    ECT技術(shù)的原理就是,根據(jù)不同相元素介電常數(shù)各不相同的特性,用電容傳感器上的電極板對被測流體施加電壓后,經(jīng)數(shù)據(jù)采集系統(tǒng)測得的流體流動產(chǎn)生的變化的電容值,通過ECT圖像重建算法,得到不同相元素的介電常數(shù)的分布,得到被測流體的截面重建圖像[19-20]。

    通常,一個 N 電極ECT系統(tǒng)中,獨立電極對的總數(shù)M滿足以下公式:

    M=C2N=N(N-1)/2(1)

    首先以這12 個極板的中某個極板為起點,依次對這 12 個極板進行編號。在測量過程中,首先編號為1的電極板被選擇為源極板,即公共電極。然后為源極板加一固定電壓值 U,再依次選取電極板

    2、3、…、12 作為檢測電極板,分別對電極板對 1-2、1-3、…1-12之間的電容值進行測量,且在每次測量時,都將其它閑置電極接地。接下來,繼續(xù)選擇電極 2作為公共電極,同理測出電極對2-3、2-4、…、2-12的電容值。依此類推,直至完成電極對 11-12 的電容值測量。此時,共能得到 66 個獨立的測量電容值。

    現(xiàn)今,大部分的ECT圖像重建算法都是基于介電常數(shù)到電容映射的線性模型。經(jīng)過離散化、線性化和歸一化的模型[21-22]如式(2)所示:

    C=SG(2)

    式中,C∈Rm表示歸一化電容向量,S∈Rm×n表示系數(shù)矩陣(靈敏度矩陣),G∈Rn表示歸一化介質(zhì)分布圖像向量。其中,ECT圖像重建的任務(wù)就是由給定的電容值C最終求解得到介電常數(shù)分布G。

    2基于Huang族校正方法的電容層析成像反問題求解

    有這樣一個目標函數(shù):f(x),設(shè)f∶Rn→R在開集DRn上二次連續(xù)可微,f在xk+1附近的二次近似為:

    f(x)≈f(xk+1)+gTk+1(x-xk+1)+

    12(x-xk+1)Lk+1(x-xk+1)(3)

    對其求導(dǎo)則有:

    g(x)≈gk+1+Lk+1(x-xk+1)(4)

    令x=xk,yk=gk+1-gk,zk=xk+1-xk,那么將有:

    L-1k+1yk≈zk(5)

    現(xiàn)在,設(shè)Tk是第k次迭代的Hesse逆近似。矩陣Tk滿足

    Tk+1yk=γzk(6)

    其中γ是一個參數(shù)。Huang族算法應(yīng)用于二次函數(shù)時,產(chǎn)生共軛方向,具有二次終止性。所有Huang族校正公式都產(chǎn)生相同的迭代點列。對于非二次函數(shù),Huang族校正公式所產(chǎn)生的點列僅依賴于參數(shù)γ[23]。

    設(shè)Huang族校正公式為:

    Tk+1=Tk+ΔM(7)

    (6)式稱為廣義擬牛頓條件。為了滿足我們(6),我們假定式(7)中的ΔM為zk與Tkyk的線性組合。即可以簡單表示為:

    ΔM=zkpTk+TkykqTk(8)

    其中,pk和qk是待定的n維向量。pk和qk滿足下列條件:

    pTkyj=0,j=0,1,…,k-1;ρ,j=k,(9)

    qTkyj=0,j=0,1,…,k-1;-1,j=k,(10)

    選擇pk和qk為zk與Tkyk的線性組合,式(9)式(10)將得到滿足。如下:

    pk=m11zk+m12TTkyk (11)

    qk=m21zk+m22TTkyk(12)

    適當選擇pk和qk滿足

    pTkyk=γ(13)

    qTkyk=-1(14)

    那么,我們得到Huang族校正公式為:

    Tk+1=Tk+zkpTk+TkykqTk(15)

    其中pk和qk分別由式(11)式(12)給出,并且滿足式(13)和式(14)。所以,可知在式(15)中,有五個參數(shù),分別為γ和mrc,這五個參數(shù)中有三個為自由參數(shù),Huang族校正公式依賴于這三個參數(shù)。

    令γ=1,并使Tk對稱,這時m12=m21,我們讓m11為自由參數(shù),這時,

    m12=m21=1-m11sTkykyTkHkyk(16)

    m22=-(1+m12sTkyk)yTkHkyk(17)

    在ECT問題中,有:

    zk=ΔG=Gk+1-Gk(18)

    yk=STS(Gk+1-Gk)(19)

    ρ=m11(ΔGTSTSΔG)2-ΔGTSTSΔG(STSΔG)TTkSTSΔG(20)

    所以此時式(15)將寫成:

    Tk+1=Tk+ΔGΔGTΔGTSTSΔG

    -TkSTSΔG(STSΔG)TT(STSΔG)TTkSTSΔG+ρrkrTk(21)

    其中

    rk=STSΔGTTkSTSΔG12ΔGΔGTSTSΔG

    TkSTSΔGSTSΔGTTkSTSΔG(22)

    下面把ECT問題使用Huang族校正來進行計算:

    1)令ek=-Tkgk;

    2)沿方向ek進行線性搜索,則有:

    Gk+1=Gk+φkek;

    3)對Tk根據(jù)Huang族校正公式進行校正,將得到Tk+1;

    初始的Hesse取為單位矩陣,并且在正定矩陣STS后加上帶系數(shù)的單位矩陣以解決在具有不適定性和病態(tài)性的ECT反問題中當靈敏度矩陣S奇異時,迭代收斂到一個非駐點的問題,所以我們最終得到的ECT反問題求解的迭代公式為下式:

    Gk+1=(STS+kH+Tk)(STC-STSGk)+Gk(23)

    3仿真與實驗結(jié)果

    基于以上理論,進行模擬仿真實驗來驗證算法求解ECT圖像重建這一反問題的效果。針對極低位層流、低位層流、小半徑核心流以及柱狀流這四種流型進行預(yù)設(shè)置并采用MATLAB配合12電極的ECT系統(tǒng)來進行仿真實驗。實驗環(huán)境設(shè)定在配置Intel(R) Core(TM) i76700 @ 34GHz處理器,安裝內(nèi)存40GB,64位操作系統(tǒng)的計算機上。

    依據(jù)以上實驗條件,根據(jù)本文中提出的Huang族校正方法進行ECT圖像重建,得到實驗結(jié)果后,將其與LBP(線性反投影)法、經(jīng)典Landweber法、CG(共軛梯度)法和SD(最速下降)法的圖像重建結(jié)果進行對比。

    ECT反問題求解中良好的精度與合理的速度是ECT技術(shù)能夠成功應(yīng)用的關(guān)鍵。在該仿真實驗中,實驗精度的確定依靠空間圖像誤差來判定。誤差計算公式(24)如下:

    ε=∑ni=1gi(img)-gi(init)∑ni=1gi(init)(24)

    式中:n表示成像區(qū)域單元總數(shù);gimg表示重建圖像向量;ginit表示介質(zhì)分布原型圖像向量;i表示成像區(qū)域剖分單元索引。

    實驗速度則依照實驗過程中的迭代次數(shù)來判定。迭代次數(shù)N越大,表示ECT圖像重建的時間花銷越大。設(shè)定實驗產(chǎn)生的迭代誤差滿足公式(25)即停止迭代。

    ‖SGk-C‖<ξ(25)

    Huang族校正算法(簡寫為Huang)成像結(jié)果與其它算法成像結(jié)果對比以及分析如下:

    從上面的成像結(jié)果對比可以看出,這幾種算法重建的圖像可以說在形狀上均大體接近原始流型,但與標準原始圖像對比上的視覺誤差各不相同。其中,對于極低位層流,Landweber方法和Huang族校正方法的成像效果相對更好,且Huang族校正方法的成像效果更佳。很明顯,對于低位層流和核心流,Huang族校正方法最接近原始圖像。

    表2展示了各種方法重建的圖像與原始圖像的誤差百分比。I 、II 、III 、IV分別代表極低位層流、低位層流、核心流以及柱狀流??梢詮闹锌闯?,對于I 、II 、III這3種流型,采用Huang族校正方法所成圖像的誤差率均為最低,與表1中所顯示結(jié)果基本一致。

    從表3可看出,對于I、II、IV 3種流型,Huang族校正方法的成像迭代次數(shù)均比Landweber方法和SD方法要少,且結(jié)合表2,可以發(fā)現(xiàn),對于這3種流型,Huang族校正方法的成像誤差百分比也是這3種方法中最低的。在流型I、II和III上,Huang族校正方法與CG法相比較,雖然迭代步數(shù)均有增加,但綜合表2我們發(fā)現(xiàn),這3種流型成像的誤差百分比卻明顯要比CG法低,特別對于流型II,CG法和Huang族校正方法的成像迭代次數(shù)明顯低于Landweber方法和SD方法,且成像誤差百分比也均低于其它方法,Huang族校正方法的迭代次數(shù)只比CG法的迭代次數(shù)增加了很少的幾步,但是在成像誤差百分比上卻又有所下降。

    綜合以上分析,可以知道采用該種Huang族校正方法來進行ECT圖像重建是一種值得考慮的方法。其在圖像重建的精度與速度上均有其優(yōu)勢。

    目標算法與其他算法的范數(shù)殘量誤差比較如下:

    圖2分別針對于以上實驗的4種流型,來對經(jīng)典Landweber、 SD、CG及Huang族校正算法的0~20步迭代的2范數(shù)殘量誤差進行比較,其中藍色曲線表示Huang族校正方法的范數(shù)殘量誤差曲線。由圖中信息可得知,對于極低位層流、層流這兩個流型,Huang族校正方法在迭代約第3步的時候就已經(jīng)得到最好的范數(shù)殘量誤差,且對于極低位層流,整體上Huang族校正方法的范數(shù)殘量誤差都要小于經(jīng)典Landweber方法和CG法;對于低位層流,Huang族校正方法的范數(shù)殘量誤差在一定范圍內(nèi)是最小的。對于核心流,Landweber方法、SD法以及Huang族校正方法,范數(shù)殘量誤差均整體在05以下,明顯要穩(wěn)定于CG法。對于柱狀流,Huang族校正方法在迭代一開始,其范數(shù)殘量誤差比其他3種方法的范數(shù)殘量誤差均要小,隨著迭代步數(shù)的增加,逐漸接近于經(jīng)典Landweber方法,并保持在其附近范圍。

    綜合以上分析,可以知道采用該種Huang族校正方法來進行ECT圖像重建是一種值得考慮的方法。其在圖像重建的精度與速度上均有其優(yōu)勢。

    4結(jié)論

    本文依據(jù)ECT圖像重建這一反問題的病態(tài)性提出了一種運用Huang族校正的方法來解決ECT圖像的重建的問題。詳細給出了該種校正方法在ECT圖像重建這一反問題中的校正公式。最后依照利用Huang族校正公式而得到的ECT反問題求解迭代公式而進行仿真實驗。仿真實驗及其結(jié)果證明了此Huang族校正方法是在解決ECT圖像重建問題上是一個有效的方法。在與LBP、經(jīng)典Landweber、CG以及SD方法的重建圖像、誤差率和迭代次數(shù)的對比上,基于Huang族校正的方法顯示出了其自身精度和速度上的優(yōu)勢。該方法高效且簡單,是一種解決ECT圖像重建問題的有效方法。

    參 考 文 獻:

    [1]吳新杰,何在剛,李惠強,等利用多準則Hopfield網(wǎng)絡(luò)對ECT進行圖像重建[J]電機與控制學報,2016,20(8):98-104

    [2]陳德運,高明,李偉,等新型ECT數(shù)據(jù)采集系統(tǒng)設(shè)計與實現(xiàn)[J]電機與控制學報,2013,17(5):87-92

    [3]田海軍,周云龍電容層析成像技術(shù)研究進展[J]化工自動化及儀表,2012,39(11):1387-1392

    [4]ISAKSEN OA Review of Reconstruction Techniques for Capacitance Tomography[J]Measurement Science and Technology,1996,7(3):325-337

    [5]FANG W,CUMBERBATCH EMatrix Properties of Data from Electrical Capacitance Tomography[J]Eng Math, 2005,51(2):127-146

    [6]SOLEIMANI M,LIONHEART W R B Nonlinear Image Reconstruction for Electrical Capacitance Tomography Using Experimental Data[J] Measurement Science and Technology,2005,16(10):1987-1996

    [7]ORTIZALEMAN C,MARTIN R,GAMIO J C,et al Reconstruction of Permittivity Images from Capacitance Tomography Data by Using Very Fast Simulated Annealing[J]Measurement Science and Technology,2004,15(7):1382-1390

    [8]陳德運,鄭貴濱,于曉洋,等基于遺傳算法電容層析成像圖像重建算法的研究[J]電機與控制學報,2003,7(3):207-211

    [9]于金平,陳德運,張云龍,等一種電容層析成像局部方差的圖像融合方法[J]哈爾濱理工大學學報,2015,20(4):116-121

    [10]XI Shaolin,ZHAO FengzhiComputation Methods for Optimization[M]ShAnghai,Shanghai Scientifie & Technical Publishers,1983:113-117

    [11]周云龍,衣得武,高云鵬基于INGA的電容層析成像圖像重建算法的研究[J]制造業(yè)自動化,2012,34(6):81-86

    [12]YANG W Q,PENG L HImage reconstruction algorithms for electrical capacitance tomography[J]Measurement Science and Technology,2003,14(1): R1~R3

    [13]于金平,陳德運,王莉莉一種基于禁忌搜索的電容層析成像圖像重建算法[J]哈爾濱理工大學學報,2016,21(1):51-56

    [14]林甲楠,陳德運,姚玉梅等一種分解型QuasiNewton電容層析成像圖像重建算法[J]哈爾濱理工大學學報,2014,19(6):44-47

    [15]王莉莉,陳德運,于曉洋等電容層析成像系統(tǒng)傳感器優(yōu)化設(shè)計[J]儀器儀表學報,2015,36(3):515-522

    [16]孫犇淵,王化祥,王丕濤基于內(nèi)部陣列電極的電容層析成像系統(tǒng)[J]傳感技術(shù)學報,2013,26(6):820-824

    [17]郭志恒,律德才,邵富群基于差分靈敏度模型的電容層析成像圖像重建方法[J]中國電機工程學報,2012,32(23):75-82

    [18]曹海燕,王化祥電容層析成像系統(tǒng)傳感器仿真研究[J]計算機工程與應(yīng)用,2012,48(15):207-211

    [19]許聞楷,陳德運,高明,等電容層析成像三維傳感器結(jié)構(gòu)設(shè)計及圖像重建[J]哈爾濱理工大學學報,2016,21(5):34-39

    [20]吉洪駐,陳德運,高明.等電容層析成像系統(tǒng)三維傳感器結(jié)構(gòu)分析[J]哈爾濱理工大學學報,2016,21(5):107-112

    [21]陳宇,陳德運基于改進RungeKutta型landweber的電容層析成像圖像重建算法[J]電機與控制學報,2014,18(7):107-112

    [22]陳宇,孫帆,張健基于自適應(yīng)權(quán)重粒子群的電容層析成像邊界灰度補償算法[J]哈爾濱理工大學學報,2010,15(3):44-49

    [23]袁亞湘,孫文瑜最優(yōu)化理論與方法[M]北京:科學出版社,1997:219-259

    [24]張云龍,陳德運,王莉莉一種基于期望最大化條件的電容層析成像圖像重建算法[J]哈爾濱理工大學學報,2016,21(2):13-19

    (編輯:關(guān)毅)

    猜你喜歡
    層析成像流型校正
    水平井油水兩相流型實驗研究
    云南化工(2021年7期)2021-12-21 07:27:32
    基于大數(shù)據(jù)量的初至層析成像算法優(yōu)化
    基于快速行進法地震層析成像研究
    劉光第《南旋記》校正
    國學(2020年1期)2020-06-29 15:15:30
    一類具有校正隔離率隨機SIQS模型的絕滅性與分布
    機內(nèi)校正
    共流型轉(zhuǎn)子的有限元分析
    基于Taitel-Dukler方法的氣液兩相流型邊界計算軟件開發(fā)
    當代化工(2015年8期)2015-02-16 08:38:44
    基于分布式無線網(wǎng)絡(luò)的無線電層析成像方法與實驗研究
    雷達學報(2014年4期)2014-04-23 07:43:22
    基于多級小波域變換的時域擴散熒光層析成像方法
    老司机影院毛片| 亚洲专区字幕在线| 精品国产超薄肉色丝袜足j| 欧美精品高潮呻吟av久久| 国产亚洲欧美在线一区二区| 亚洲熟女精品中文字幕| 久久午夜综合久久蜜桃| 一级a爱视频在线免费观看| 啦啦啦免费观看视频1| 国产免费现黄频在线看| 国产精品成人在线| 欧美国产精品一级二级三级| www.熟女人妻精品国产| 久久久久网色| 国产亚洲精品第一综合不卡| 成在线人永久免费视频| 99久久国产精品久久久| 国产1区2区3区精品| 超碰97精品在线观看| 在线永久观看黄色视频| 满18在线观看网站| 免费不卡黄色视频| 午夜91福利影院| 成年av动漫网址| 久久久水蜜桃国产精品网| 嫁个100分男人电影在线观看| 日本猛色少妇xxxxx猛交久久| 欧美性长视频在线观看| 黑人巨大精品欧美一区二区蜜桃| av有码第一页| 国产极品粉嫩免费观看在线| 久久 成人 亚洲| 美女视频免费永久观看网站| 美女高潮喷水抽搐中文字幕| 久久精品亚洲av国产电影网| 人人澡人人妻人| 日韩中文字幕欧美一区二区| 精品亚洲成国产av| 少妇精品久久久久久久| 男女床上黄色一级片免费看| 欧美成狂野欧美在线观看| 国产男女内射视频| 国产精品一区二区精品视频观看| 国产主播在线观看一区二区| 一级片免费观看大全| av网站在线播放免费| 精品国产国语对白av| 自线自在国产av| 欧美激情极品国产一区二区三区| 操出白浆在线播放| 少妇猛男粗大的猛烈进出视频| 免费在线观看影片大全网站| 啦啦啦 在线观看视频| 久久香蕉激情| 啦啦啦视频在线资源免费观看| 女人被躁到高潮嗷嗷叫费观| 啦啦啦 在线观看视频| 美女主播在线视频| 真人做人爱边吃奶动态| 成人国产一区最新在线观看| 国产免费现黄频在线看| 国产区一区二久久| 亚洲自偷自拍图片 自拍| 精品亚洲成a人片在线观看| 视频区欧美日本亚洲| 国产极品粉嫩免费观看在线| 精品卡一卡二卡四卡免费| 嫩草影视91久久| 日韩熟女老妇一区二区性免费视频| 亚洲天堂av无毛| 脱女人内裤的视频| 久久久久久人人人人人| 女性生殖器流出的白浆| 亚洲伊人久久精品综合| 又黄又粗又硬又大视频| 精品视频人人做人人爽| 大型av网站在线播放| 最黄视频免费看| av在线播放精品| 12—13女人毛片做爰片一| 欧美日本中文国产一区发布| 精品熟女少妇八av免费久了| 免费久久久久久久精品成人欧美视频| 丝袜脚勾引网站| 一区二区日韩欧美中文字幕| 桃红色精品国产亚洲av| 人人澡人人妻人| 国产成人欧美| 人人妻人人澡人人爽人人夜夜| 精品国产一区二区久久| 亚洲自偷自拍图片 自拍| 天天影视国产精品| 亚洲国产欧美一区二区综合| 国产男女内射视频| 国产亚洲一区二区精品| 91精品三级在线观看| 午夜精品久久久久久毛片777| 精品国内亚洲2022精品成人 | 中亚洲国语对白在线视频| svipshipincom国产片| 欧美日韩精品网址| 1024香蕉在线观看| 久9热在线精品视频| 脱女人内裤的视频| 最近中文字幕2019免费版| 亚洲av男天堂| 老司机亚洲免费影院| 精品国产一区二区三区四区第35| 午夜福利在线观看吧| 国产精品一区二区在线不卡| 欧美日韩亚洲高清精品| 亚洲av电影在线进入| 日本撒尿小便嘘嘘汇集6| 国产在视频线精品| 国产精品香港三级国产av潘金莲| 久久免费观看电影| 50天的宝宝边吃奶边哭怎么回事| 亚洲专区中文字幕在线| 免费观看av网站的网址| 夫妻午夜视频| 精品国产超薄肉色丝袜足j| 岛国毛片在线播放| av片东京热男人的天堂| 在线 av 中文字幕| 91麻豆精品激情在线观看国产 | 国产欧美日韩一区二区精品| 99国产精品99久久久久| 亚洲精品一区蜜桃| 日本a在线网址| 午夜成年电影在线免费观看| 另类精品久久| 人人妻人人澡人人爽人人夜夜| 亚洲精品乱久久久久久| 下体分泌物呈黄色| 99精品欧美一区二区三区四区| 一区在线观看完整版| 中文字幕av电影在线播放| 新久久久久国产一级毛片| 51午夜福利影视在线观看| 国产亚洲午夜精品一区二区久久| 啦啦啦啦在线视频资源| 一个人免费在线观看的高清视频 | 亚洲国产成人一精品久久久| 国产精品自产拍在线观看55亚洲 | 两性夫妻黄色片| 亚洲精品久久午夜乱码| 99国产极品粉嫩在线观看| 久久人妻熟女aⅴ| 多毛熟女@视频| 不卡一级毛片| 国产精品.久久久| 老司机午夜福利在线观看视频 | 成人国产一区最新在线观看| 纵有疾风起免费观看全集完整版| 色视频在线一区二区三区| 欧美+亚洲+日韩+国产| 免费女性裸体啪啪无遮挡网站| 两人在一起打扑克的视频| 欧美久久黑人一区二区| 国产精品1区2区在线观看. | 美女福利国产在线| 在线观看一区二区三区激情| 大香蕉久久成人网| 天天躁日日躁夜夜躁夜夜| 青春草视频在线免费观看| 亚洲精品久久久久久婷婷小说| 看免费av毛片| 国产xxxxx性猛交| 视频在线观看一区二区三区| 啦啦啦在线免费观看视频4| 国产亚洲精品第一综合不卡| 又大又爽又粗| 免费一级毛片在线播放高清视频 | 高清黄色对白视频在线免费看| av天堂久久9| 黑人巨大精品欧美一区二区蜜桃| 国产福利在线免费观看视频| 午夜免费鲁丝| 少妇精品久久久久久久| 久久九九热精品免费| 亚洲中文日韩欧美视频| 男女无遮挡免费网站观看| 久久九九热精品免费| 99久久99久久久精品蜜桃| 久久亚洲精品不卡| bbb黄色大片| 99国产精品一区二区三区| 18禁国产床啪视频网站| 少妇人妻久久综合中文| 下体分泌物呈黄色| 成人亚洲精品一区在线观看| √禁漫天堂资源中文www| 老汉色∧v一级毛片| 精品国产一区二区三区久久久樱花| 免费av中文字幕在线| 国产成人欧美| 久久久久久久久久久久大奶| 少妇猛男粗大的猛烈进出视频| 精品国产一区二区三区四区第35| 久久久国产成人免费| 午夜福利影视在线免费观看| 午夜福利,免费看| 菩萨蛮人人尽说江南好唐韦庄| 久久精品亚洲熟妇少妇任你| 久久久久精品人妻al黑| 日韩中文字幕视频在线看片| 女人被躁到高潮嗷嗷叫费观| 亚洲国产欧美一区二区综合| 精品卡一卡二卡四卡免费| 在线永久观看黄色视频| 国产成人av教育| 欧美日韩精品网址| 老司机影院毛片| 欧美精品av麻豆av| 91成人精品电影| 免费看十八禁软件| 亚洲欧美清纯卡通| 少妇被粗大的猛进出69影院| 国产成人一区二区三区免费视频网站| 国产日韩欧美视频二区| 香蕉丝袜av| 精品第一国产精品| 精品视频人人做人人爽| 黄色片一级片一级黄色片| 老熟妇乱子伦视频在线观看 | 欧美大码av| 99久久人妻综合| a级毛片在线看网站| 热re99久久国产66热| 午夜福利一区二区在线看| 18禁黄网站禁片午夜丰满| 亚洲专区国产一区二区| 免费观看人在逋| 国产一区二区 视频在线| 日本91视频免费播放| 日本av免费视频播放| 欧美精品啪啪一区二区三区 | av一本久久久久| 男人舔女人的私密视频| 嫁个100分男人电影在线观看| 无限看片的www在线观看| 亚洲,欧美精品.| 女人精品久久久久毛片| 丝袜在线中文字幕| 亚洲黑人精品在线| 一区二区三区乱码不卡18| 免费在线观看完整版高清| 狂野欧美激情性xxxx| 婷婷成人精品国产| 电影成人av| 王馨瑶露胸无遮挡在线观看| 欧美激情 高清一区二区三区| 少妇人妻久久综合中文| 悠悠久久av| 电影成人av| 在线观看舔阴道视频| 亚洲第一青青草原| 中文精品一卡2卡3卡4更新| 美女视频免费永久观看网站| 丝袜人妻中文字幕| av电影中文网址| 亚洲精品美女久久久久99蜜臀| 精品欧美一区二区三区在线| 极品少妇高潮喷水抽搐| 精品国产乱码久久久久久小说| 欧美日韩精品网址| 国产免费福利视频在线观看| netflix在线观看网站| 免费在线观看完整版高清| av又黄又爽大尺度在线免费看| 久久狼人影院| 菩萨蛮人人尽说江南好唐韦庄| 另类精品久久| 色精品久久人妻99蜜桃| 久久精品亚洲熟妇少妇任你| 亚洲性夜色夜夜综合| 91老司机精品| 另类亚洲欧美激情| 亚洲国产中文字幕在线视频| 国产野战对白在线观看| 女警被强在线播放| 不卡一级毛片| 十八禁人妻一区二区| 99国产精品一区二区三区| 性高湖久久久久久久久免费观看| 亚洲欧美日韩高清在线视频 | 一级,二级,三级黄色视频| 欧美乱码精品一区二区三区| av欧美777| 香蕉国产在线看| 午夜激情av网站| 少妇被粗大的猛进出69影院| 日本wwww免费看| 久久久久久免费高清国产稀缺| 亚洲av日韩在线播放| 黑人操中国人逼视频| 一本一本久久a久久精品综合妖精| 精品欧美一区二区三区在线| 欧美av亚洲av综合av国产av| 亚洲欧美成人综合另类久久久| 国产亚洲一区二区精品| 又大又爽又粗| 亚洲精品成人av观看孕妇| 黄频高清免费视频| 欧美日韩国产mv在线观看视频| 99国产精品免费福利视频| 国产成人啪精品午夜网站| 在线观看人妻少妇| av片东京热男人的天堂| 精品人妻1区二区| 国产一区二区三区av在线| 亚洲欧美精品综合一区二区三区| cao死你这个sao货| 亚洲精品国产av成人精品| 色94色欧美一区二区| 日本av手机在线免费观看| av视频免费观看在线观看| 日韩中文字幕视频在线看片| 波多野结衣av一区二区av| 18禁黄网站禁片午夜丰满| 999久久久精品免费观看国产| 国产精品久久久久久人妻精品电影 | 一本大道久久a久久精品| 国产三级黄色录像| 真人做人爱边吃奶动态| 亚洲中文av在线| 色视频在线一区二区三区| 热99久久久久精品小说推荐| av一本久久久久| 成年动漫av网址| 少妇的丰满在线观看| 免费在线观看影片大全网站| 免费观看a级毛片全部| 亚洲三区欧美一区| 免费高清在线观看视频在线观看| 麻豆乱淫一区二区| 亚洲精品自拍成人| 午夜两性在线视频| 18在线观看网站| 精品福利观看| 丝袜美腿诱惑在线| 黑人操中国人逼视频| 视频区欧美日本亚洲| av国产精品久久久久影院| 国产精品二区激情视频| 色视频在线一区二区三区| 欧美激情 高清一区二区三区| 国产片内射在线| 一区二区日韩欧美中文字幕| 精品一区在线观看国产| 欧美老熟妇乱子伦牲交| 我的亚洲天堂| www.999成人在线观看| 青青草视频在线视频观看| 另类亚洲欧美激情| 国产一区二区 视频在线| 久久天堂一区二区三区四区| 国产精品一区二区免费欧美 | 91精品伊人久久大香线蕉| 十八禁网站网址无遮挡| 久久久久久久精品精品| 十八禁网站免费在线| 人人澡人人妻人| 欧美黑人欧美精品刺激| xxxhd国产人妻xxx| 女性生殖器流出的白浆| 精品卡一卡二卡四卡免费| 国产精品一区二区在线观看99| 淫妇啪啪啪对白视频 | 一级毛片电影观看| 少妇被粗大的猛进出69影院| 久久影院123| 欧美激情 高清一区二区三区| 久久午夜综合久久蜜桃| 国产欧美日韩精品亚洲av| 人人澡人人妻人| 91精品三级在线观看| 久久人人爽av亚洲精品天堂| 欧美激情极品国产一区二区三区| 美女福利国产在线| 日韩欧美免费精品| 色老头精品视频在线观看| 纵有疾风起免费观看全集完整版| 欧美精品亚洲一区二区| 国产精品 国内视频| 9热在线视频观看99| av超薄肉色丝袜交足视频| 精品福利观看| 在线观看免费日韩欧美大片| 老司机靠b影院| 国产精品偷伦视频观看了| 亚洲中文字幕日韩| 99香蕉大伊视频| 国产高清视频在线播放一区 | 亚洲欧美色中文字幕在线| 少妇的丰满在线观看| 欧美av亚洲av综合av国产av| www.精华液| 午夜免费观看性视频| 男女国产视频网站| 国产91精品成人一区二区三区 | 黄色 视频免费看| 成年人黄色毛片网站| 每晚都被弄得嗷嗷叫到高潮| 国产欧美日韩综合在线一区二区| 新久久久久国产一级毛片| 亚洲精品乱久久久久久| 日韩电影二区| 久久久久精品国产欧美久久久 | 国产成人啪精品午夜网站| 80岁老熟妇乱子伦牲交| 9色porny在线观看| 亚洲精品国产精品久久久不卡| 中文字幕高清在线视频| 久久国产精品大桥未久av| 国产成+人综合+亚洲专区| 亚洲欧美一区二区三区黑人| 欧美精品亚洲一区二区| 成人手机av| 亚洲精品中文字幕在线视频| 操出白浆在线播放| 国产精品麻豆人妻色哟哟久久| 婷婷丁香在线五月| 国产又色又爽无遮挡免| 亚洲 欧美一区二区三区| 国产欧美亚洲国产| netflix在线观看网站| 肉色欧美久久久久久久蜜桃| 手机成人av网站| 亚洲av国产av综合av卡| 成人av一区二区三区在线看 | 久久亚洲精品不卡| 天堂8中文在线网| 岛国毛片在线播放| 无限看片的www在线观看| 男女下面插进去视频免费观看| 一二三四在线观看免费中文在| 最新的欧美精品一区二区| 婷婷色av中文字幕| 日韩一卡2卡3卡4卡2021年| 天天影视国产精品| 中文字幕最新亚洲高清| 亚洲国产日韩一区二区| 黑人操中国人逼视频| 在线看a的网站| 欧美精品亚洲一区二区| 一区在线观看完整版| 国产深夜福利视频在线观看| 精品久久蜜臀av无| 国产不卡av网站在线观看| 国产99久久九九免费精品| 国产精品.久久久| 一本久久精品| 色精品久久人妻99蜜桃| 午夜视频精品福利| 亚洲av男天堂| 成年女人毛片免费观看观看9 | 高潮久久久久久久久久久不卡| 五月开心婷婷网| 一本综合久久免费| 亚洲精品自拍成人| 老熟妇仑乱视频hdxx| 精品一区二区三区四区五区乱码| 久久久国产欧美日韩av| 亚洲精品国产色婷婷电影| 男女床上黄色一级片免费看| 在线观看舔阴道视频| 高清黄色对白视频在线免费看| 麻豆乱淫一区二区| 午夜老司机福利片| 欧美97在线视频| 国产精品一区二区在线不卡| 成人亚洲精品一区在线观看| 男人爽女人下面视频在线观看| 18在线观看网站| 在线亚洲精品国产二区图片欧美| 少妇的丰满在线观看| 国产免费av片在线观看野外av| 久久久久久久精品精品| 久久人妻福利社区极品人妻图片| 欧美精品亚洲一区二区| 久久香蕉激情| 欧美精品高潮呻吟av久久| 精品一区在线观看国产| 性高湖久久久久久久久免费观看| 色婷婷久久久亚洲欧美| a在线观看视频网站| 少妇粗大呻吟视频| 亚洲国产欧美在线一区| 美女高潮喷水抽搐中文字幕| 精品久久久精品久久久| 久久精品成人免费网站| 久久久久久久国产电影| 久久久国产一区二区| 免费在线观看完整版高清| 亚洲av片天天在线观看| 老鸭窝网址在线观看| 999久久久精品免费观看国产| 中文欧美无线码| www.自偷自拍.com| 国产成人啪精品午夜网站| 十八禁网站网址无遮挡| 欧美性长视频在线观看| 亚洲美女黄色视频免费看| 精品人妻1区二区| 夫妻午夜视频| 亚洲av日韩精品久久久久久密| 老司机深夜福利视频在线观看 | 亚洲精品国产av蜜桃| 亚洲第一欧美日韩一区二区三区 | 午夜福利在线免费观看网站| 水蜜桃什么品种好| 人人妻人人澡人人看| 男女无遮挡免费网站观看| 精品少妇久久久久久888优播| 国产有黄有色有爽视频| 国产精品一二三区在线看| 免费不卡黄色视频| 18禁国产床啪视频网站| 麻豆乱淫一区二区| 国产一区二区三区av在线| 无遮挡黄片免费观看| 久久午夜综合久久蜜桃| 亚洲九九香蕉| 建设人人有责人人尽责人人享有的| 视频区图区小说| 国产xxxxx性猛交| 久久99热这里只频精品6学生| 天天添夜夜摸| 亚洲中文字幕日韩| 亚洲免费av在线视频| 久久中文看片网| 午夜福利影视在线免费观看| 一区二区三区精品91| 精品亚洲乱码少妇综合久久| 欧美精品av麻豆av| 亚洲九九香蕉| 丰满饥渴人妻一区二区三| 免费不卡黄色视频| 欧美老熟妇乱子伦牲交| 亚洲精品国产av成人精品| 久久精品亚洲av国产电影网| 91成人精品电影| 国产精品一区二区免费欧美 | 亚洲精品第二区| 国产一区二区在线观看av| 久久久欧美国产精品| 精品少妇内射三级| 麻豆av在线久日| 18禁观看日本| 男女高潮啪啪啪动态图| 午夜福利在线免费观看网站| 国产亚洲av高清不卡| 国产97色在线日韩免费| 亚洲av美国av| 亚洲自偷自拍图片 自拍| 香蕉国产在线看| 欧美在线黄色| av国产精品久久久久影院| 久久99一区二区三区| 国产欧美日韩一区二区三区在线| www.精华液| 久久女婷五月综合色啪小说| a级毛片黄视频| 亚洲人成电影免费在线| 久久精品熟女亚洲av麻豆精品| 午夜福利视频精品| 中国美女看黄片| 国产亚洲av片在线观看秒播厂| 我的亚洲天堂| 亚洲黑人精品在线| 日韩欧美国产一区二区入口| 伦理电影免费视频| 美女脱内裤让男人舔精品视频| 一区二区三区激情视频| 婷婷丁香在线五月| av有码第一页| 秋霞在线观看毛片| 首页视频小说图片口味搜索| 欧美日本中文国产一区发布| 亚洲欧美精品综合一区二区三区| 老司机午夜十八禁免费视频| 宅男免费午夜| 不卡av一区二区三区| 免费在线观看黄色视频的| 热99久久久久精品小说推荐| 国产在视频线精品| 别揉我奶头~嗯~啊~动态视频 | 一级毛片女人18水好多| 国产精品久久久久久人妻精品电影 | 国产精品麻豆人妻色哟哟久久| 精品少妇久久久久久888优播| 交换朋友夫妻互换小说| 日韩制服丝袜自拍偷拍| 成人18禁高潮啪啪吃奶动态图| 亚洲精品中文字幕一二三四区 | 视频区欧美日本亚洲| 黑人巨大精品欧美一区二区mp4| 亚洲精品乱久久久久久| 啦啦啦 在线观看视频| 亚洲五月色婷婷综合| 国产精品亚洲av一区麻豆| www.av在线官网国产| 黑人巨大精品欧美一区二区mp4| 18禁裸乳无遮挡动漫免费视频| 黄色视频不卡| 秋霞在线观看毛片| 精品久久蜜臀av无| 色94色欧美一区二区| 丰满饥渴人妻一区二区三| av网站在线播放免费| 天天添夜夜摸| 欧美亚洲日本最大视频资源| 欧美中文综合在线视频| 美女脱内裤让男人舔精品视频|