• <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
    基于多級小波域變換的時域擴散熒光層析成像方法
    免费不卡的大黄色大毛片视频在线观看 | 日本免费一区二区三区高清不卡| 97超碰精品成人国产| 亚洲av成人精品一区久久| 国产综合懂色| 亚洲av不卡在线观看| 午夜福利成人在线免费观看| 午夜精品国产一区二区电影 | 免费av不卡在线播放| 国产成人freesex在线 | 变态另类丝袜制服| 麻豆国产97在线/欧美| 女生性感内裤真人,穿戴方法视频| 舔av片在线| 亚洲七黄色美女视频| 一区二区三区四区激情视频 | 国产精品一区二区三区四区久久| a级毛片免费高清观看在线播放| 听说在线观看完整版免费高清| 亚洲七黄色美女视频| 亚洲七黄色美女视频| 国产一区二区三区在线臀色熟女| 国产一区二区三区在线臀色熟女| 国产av不卡久久| 亚洲国产欧美人成| 亚洲最大成人av| 在线观看美女被高潮喷水网站| 日日摸夜夜添夜夜爱| videossex国产| 亚洲国产色片| 男人狂女人下面高潮的视频| 国产精品永久免费网站| 最近最新中文字幕大全电影3| 一级毛片aaaaaa免费看小| 中文字幕av在线有码专区| 亚洲精华国产精华液的使用体验 | 91在线观看av| 欧美绝顶高潮抽搐喷水| 亚洲欧美精品自产自拍| 国产高潮美女av| av在线观看视频网站免费| 深夜a级毛片| 亚洲熟妇熟女久久| 日本爱情动作片www.在线观看 | 成人精品一区二区免费| 午夜日韩欧美国产| 亚洲性久久影院| 人人妻人人澡人人爽人人夜夜 | 日本成人三级电影网站| 亚洲中文字幕日韩| 日韩一本色道免费dvd| 少妇人妻一区二区三区视频| 国产不卡一卡二| 香蕉av资源在线| 久久精品影院6| 亚洲精品粉嫩美女一区| 露出奶头的视频| 亚洲精品影视一区二区三区av| 亚洲国产精品成人久久小说 | 人妻久久中文字幕网| 久久久a久久爽久久v久久| 日韩一区二区视频免费看| 久久精品国产亚洲av涩爱 | 在线a可以看的网站| 久久久久久九九精品二区国产| 国产精品福利在线免费观看| 菩萨蛮人人尽说江南好唐韦庄 | 欧美另类亚洲清纯唯美| 国产精品一区二区免费欧美| 国产精品国产高清国产av| 精品一区二区三区人妻视频| 日韩精品青青久久久久久| 精品久久久久久久久亚洲| 美女黄网站色视频| 噜噜噜噜噜久久久久久91| 欧美在线一区亚洲| 99热这里只有是精品在线观看| 久久久色成人| 啦啦啦啦在线视频资源| 老师上课跳d突然被开到最大视频| 成人一区二区视频在线观看| 国产免费男女视频| 欧美最新免费一区二区三区| 成人av一区二区三区在线看| 国产午夜精品久久久久久一区二区三区 | 永久网站在线| 日韩精品青青久久久久久| 久久精品久久久久久噜噜老黄 | 欧美日本视频| 国产白丝娇喘喷水9色精品| 寂寞人妻少妇视频99o| 日本黄色视频三级网站网址| 亚洲中文字幕日韩| 一进一出抽搐gif免费好疼| 欧美+亚洲+日韩+国产| 3wmmmm亚洲av在线观看| 国产精品国产三级国产av玫瑰| 精品久久国产蜜桃| 91久久精品国产一区二区三区| 一级毛片电影观看 | 校园春色视频在线观看| 大又大粗又爽又黄少妇毛片口| 搞女人的毛片| 国产午夜福利久久久久久| 日本色播在线视频| 亚洲av成人av| 男女视频在线观看网站免费| 久久久久性生活片| 亚洲不卡免费看| 天堂av国产一区二区熟女人妻| 亚洲国产精品合色在线| 亚洲av中文av极速乱| 99久久中文字幕三级久久日本| 综合色av麻豆| 国产成人影院久久av| 亚洲成a人片在线一区二区| 男女那种视频在线观看| 久久鲁丝午夜福利片| 身体一侧抽搐| 午夜a级毛片| 97在线视频观看| 美女内射精品一级片tv| 成人av在线播放网站| 亚洲精品乱码久久久v下载方式| 免费一级毛片在线播放高清视频| 变态另类成人亚洲欧美熟女| 国产av在哪里看| 亚洲av美国av| 欧美性猛交╳xxx乱大交人| 狠狠狠狠99中文字幕| 可以在线观看的亚洲视频| 久久精品国产99精品国产亚洲性色| 国产激情偷乱视频一区二区| 国产精品久久久久久久电影| av天堂在线播放| 哪里可以看免费的av片| 午夜福利在线观看吧| 欧美高清性xxxxhd video| 国产在线男女| 亚洲经典国产精华液单| 18禁黄网站禁片免费观看直播| 久久久久久久午夜电影| 六月丁香七月| 热99re8久久精品国产| 亚洲第一电影网av| 久久亚洲国产成人精品v| 全区人妻精品视频| 亚洲乱码一区二区免费版| 深夜精品福利| 99九九线精品视频在线观看视频| 俺也久久电影网| 中出人妻视频一区二区| 亚洲人成网站高清观看| 色综合色国产| av中文乱码字幕在线| av免费在线看不卡| 国产精品爽爽va在线观看网站| 国产成人a∨麻豆精品| 99热这里只有是精品在线观看| 天堂网av新在线| 黄色一级大片看看| 午夜福利高清视频| 99久久中文字幕三级久久日本| 成人午夜高清在线视频| 欧美xxxx性猛交bbbb| 日韩欧美国产在线观看| 成人永久免费在线观看视频| 日本黄大片高清| 亚洲成人中文字幕在线播放| 亚洲熟妇中文字幕五十中出| 色播亚洲综合网| 可以在线观看的亚洲视频| 国产大屁股一区二区在线视频| 日韩,欧美,国产一区二区三区 | 国模一区二区三区四区视频| 亚洲欧美精品综合久久99| 美女被艹到高潮喷水动态| 在线观看免费视频日本深夜| 少妇丰满av| 欧美性猛交╳xxx乱大交人| 最近在线观看免费完整版| 国产乱人偷精品视频| 免费av不卡在线播放| 少妇丰满av| 99热精品在线国产| 看片在线看免费视频| 全区人妻精品视频| 亚洲成人中文字幕在线播放| 深爱激情五月婷婷| 麻豆乱淫一区二区| 成人美女网站在线观看视频| 在线观看66精品国产| 一级a爱片免费观看的视频| 色综合站精品国产| 亚洲aⅴ乱码一区二区在线播放| 欧美成人精品欧美一级黄| 久久久久久久久久成人| 国产精品无大码| 观看美女的网站| 欧美3d第一页| 久久久久免费精品人妻一区二区| 男人舔奶头视频| 成人漫画全彩无遮挡| 能在线免费观看的黄片| 久久人妻av系列| 免费无遮挡裸体视频| 精华霜和精华液先用哪个| 久久久精品94久久精品| 午夜精品在线福利| 国产又黄又爽又无遮挡在线| 国产色爽女视频免费观看| 免费av不卡在线播放| www日本黄色视频网| 欧美精品国产亚洲| 国产av不卡久久| 嫩草影院精品99| 天美传媒精品一区二区| 久久久久性生活片| 小蜜桃在线观看免费完整版高清| 亚洲中文字幕日韩| 3wmmmm亚洲av在线观看| 91久久精品国产一区二区三区| 国产精品一区www在线观看| 亚洲天堂国产精品一区在线| 亚洲色图av天堂| 最好的美女福利视频网| 色综合站精品国产| 久99久视频精品免费| 淫妇啪啪啪对白视频| 亚洲中文字幕日韩| 色综合亚洲欧美另类图片| 日日摸夜夜添夜夜爱| 老司机午夜福利在线观看视频| 免费看光身美女| 亚洲乱码一区二区免费版| 俄罗斯特黄特色一大片| 女人被狂操c到高潮| 午夜免费男女啪啪视频观看 | 国产麻豆成人av免费视频| 国产亚洲精品综合一区在线观看| 亚洲中文字幕一区二区三区有码在线看| 香蕉av资源在线| 亚洲欧美清纯卡通| 亚洲在线自拍视频| ponron亚洲| 长腿黑丝高跟| 欧美zozozo另类| 午夜久久久久精精品| 国产亚洲精品久久久久久毛片| 亚洲欧美成人综合另类久久久 | 国产免费一级a男人的天堂| 麻豆国产97在线/欧美| 日本黄色片子视频| 99久久九九国产精品国产免费| 看非洲黑人一级黄片| 99国产精品一区二区蜜桃av| 国产aⅴ精品一区二区三区波| 夜夜看夜夜爽夜夜摸| 精品久久国产蜜桃| 99热6这里只有精品| 一个人看视频在线观看www免费| 国产高清激情床上av| 高清毛片免费观看视频网站| 欧美一级a爱片免费观看看| 日韩一本色道免费dvd| 91久久精品国产一区二区三区| 又爽又黄无遮挡网站| 人妻制服诱惑在线中文字幕| 久久这里只有精品中国| 亚洲高清免费不卡视频| 日本一本二区三区精品| 三级国产精品欧美在线观看| 99在线人妻在线中文字幕| aaaaa片日本免费| 国产探花极品一区二区| 3wmmmm亚洲av在线观看| 亚洲精品久久国产高清桃花| 99久久成人亚洲精品观看| 国产亚洲精品久久久久久毛片| 国产亚洲91精品色在线| 久久人人爽人人片av| 亚洲国产精品合色在线| 日本成人三级电影网站| 国产高清视频在线观看网站| 日韩 亚洲 欧美在线| 内射极品少妇av片p| 国产伦精品一区二区三区视频9| av视频在线观看入口| 亚洲不卡免费看| 你懂的网址亚洲精品在线观看 | 白带黄色成豆腐渣| 在线免费观看的www视频| 能在线免费观看的黄片| 午夜精品一区二区三区免费看| 日产精品乱码卡一卡2卡三| 色播亚洲综合网| 一本精品99久久精品77| 寂寞人妻少妇视频99o| 欧美一区二区国产精品久久精品| 99热这里只有是精品在线观看| 国产免费男女视频| 婷婷六月久久综合丁香| 国内精品宾馆在线| a级一级毛片免费在线观看| 日韩精品青青久久久久久| 狠狠狠狠99中文字幕| 亚洲av电影不卡..在线观看| 一个人看的www免费观看视频| 国产av不卡久久| 最新在线观看一区二区三区| 国产熟女欧美一区二区| 久久久久久久久久黄片| aaaaa片日本免费| 99久久九九国产精品国产免费| 性色avwww在线观看| 日韩欧美国产在线观看| 日韩人妻高清精品专区| 日本a在线网址| 精品久久久久久久久久免费视频| 人人妻人人澡欧美一区二区| 日韩av在线大香蕉| 大又大粗又爽又黄少妇毛片口| 国产三级在线视频| 久久天躁狠狠躁夜夜2o2o| 成人综合一区亚洲| 2021天堂中文幕一二区在线观| 最后的刺客免费高清国语| 欧美潮喷喷水| 3wmmmm亚洲av在线观看| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品av在线| 亚洲欧美成人综合另类久久久 | 2021天堂中文幕一二区在线观| 我要看日韩黄色一级片| 国产在视频线在精品| 日日啪夜夜撸| 日本一二三区视频观看| 青春草视频在线免费观看| 日韩欧美精品v在线| 天天躁日日操中文字幕| 亚洲av中文字字幕乱码综合| 婷婷六月久久综合丁香| 国产人妻一区二区三区在| 国产探花极品一区二区| 亚洲精品色激情综合| 老司机午夜福利在线观看视频| 亚洲欧美日韩东京热| 国产精品综合久久久久久久免费| 成人漫画全彩无遮挡| 欧洲精品卡2卡3卡4卡5卡区| 欧美日韩一区二区视频在线观看视频在线 | 99热这里只有是精品在线观看| 欧美日韩一区二区视频在线观看视频在线 | 最新在线观看一区二区三区| 女同久久另类99精品国产91| 亚洲av成人av| 麻豆乱淫一区二区| 在线a可以看的网站| 欧美一区二区亚洲| 国产高清三级在线| 精品人妻偷拍中文字幕| 久久6这里有精品| 大又大粗又爽又黄少妇毛片口| 亚洲成人中文字幕在线播放| av女优亚洲男人天堂| 婷婷亚洲欧美| 色综合站精品国产| 九九爱精品视频在线观看| 精品熟女少妇av免费看| 天美传媒精品一区二区| 国产色爽女视频免费观看| 嫩草影视91久久| 一区二区三区免费毛片| 亚洲精华国产精华液的使用体验 | 久久精品夜夜夜夜夜久久蜜豆| 精品欧美国产一区二区三| 欧美+日韩+精品| 久久精品国产亚洲av香蕉五月| 简卡轻食公司| 亚洲最大成人中文| 变态另类丝袜制服| 精品日产1卡2卡| 日韩av不卡免费在线播放| 夜夜爽天天搞| 国产av麻豆久久久久久久| 亚洲欧美成人综合另类久久久 | 久久这里只有精品中国| 国内精品宾馆在线| 国产探花极品一区二区| 十八禁网站免费在线| 亚洲精品乱码久久久v下载方式| 亚洲精品日韩在线中文字幕 | 国产成人aa在线观看| 亚洲专区国产一区二区| 欧美一区二区亚洲| 国产精品电影一区二区三区| 99久久精品热视频| 国产精品一二三区在线看| 国产极品精品免费视频能看的| av在线天堂中文字幕| 国产色婷婷99| 亚洲欧美成人综合另类久久久 | 美女黄网站色视频| 久久久国产成人精品二区| 亚洲电影在线观看av| 亚洲av第一区精品v没综合| 国产av麻豆久久久久久久| 亚洲av美国av| 免费看日本二区| av天堂在线播放| 毛片女人毛片| 国产精品1区2区在线观看.| 最近视频中文字幕2019在线8| 国产在线精品亚洲第一网站| 亚洲欧美日韩无卡精品| 午夜爱爱视频在线播放| 国产高清不卡午夜福利| 69人妻影院| 赤兔流量卡办理| 91麻豆精品激情在线观看国产| aaaaa片日本免费| 悠悠久久av| 中国美女看黄片| 搞女人的毛片| 色综合站精品国产| 久久久精品大字幕| 色噜噜av男人的天堂激情| 色视频www国产| av在线蜜桃| 岛国在线免费视频观看| 亚洲中文日韩欧美视频| 日本精品一区二区三区蜜桃| 别揉我奶头 嗯啊视频| 国产精品无大码| 非洲黑人性xxxx精品又粗又长| .国产精品久久| 国产高清有码在线观看视频| 成人欧美大片| 日韩欧美三级三区| 级片在线观看| 日韩制服骚丝袜av| 99热这里只有是精品在线观看| 综合色av麻豆| 俺也久久电影网| 国语自产精品视频在线第100页| 免费大片18禁| 看片在线看免费视频| 99久久中文字幕三级久久日本| 又爽又黄无遮挡网站| 婷婷精品国产亚洲av在线| 99在线视频只有这里精品首页| 在线观看66精品国产| 国产黄色小视频在线观看| 国产在线男女| 日韩欧美三级三区| 国产精品久久电影中文字幕| 亚洲精品日韩av片在线观看| 99热6这里只有精品| 国产高清视频在线播放一区| 免费观看的影片在线观看| 国产精品亚洲一级av第二区| 91久久精品国产一区二区成人| 中文资源天堂在线| 色综合站精品国产| 国产色婷婷99| 国产蜜桃级精品一区二区三区| 老司机福利观看| 亚洲成人精品中文字幕电影| 少妇熟女欧美另类| 国产精品一区二区性色av| 高清毛片免费观看视频网站| 欧美激情国产日韩精品一区| 国产成人福利小说| 成人美女网站在线观看视频| 婷婷六月久久综合丁香| 九九爱精品视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 中文字幕熟女人妻在线| 大又大粗又爽又黄少妇毛片口| av天堂在线播放| 亚洲五月天丁香| 国产午夜精品论理片| 九色成人免费人妻av| 丝袜喷水一区| 免费av毛片视频| 国产精品久久电影中文字幕| 欧美日韩精品成人综合77777| 一级av片app| 午夜福利18| 日本黄色片子视频| 村上凉子中文字幕在线| aaaaa片日本免费| 亚洲精品一区av在线观看| 欧美色欧美亚洲另类二区| 人人妻人人澡欧美一区二区| 日本撒尿小便嘘嘘汇集6| av天堂在线播放| 99热全是精品| 国产精品亚洲一级av第二区| 亚洲av熟女| 国产成人freesex在线 | 中文字幕av在线有码专区| 国产一区二区三区在线臀色熟女| 国产在线男女| 日本精品一区二区三区蜜桃| 国产 一区 欧美 日韩| 国产成年人精品一区二区| 亚洲人成网站高清观看| 日韩三级伦理在线观看| 午夜精品一区二区三区免费看| 日本精品一区二区三区蜜桃| 欧美最黄视频在线播放免费| 国产av不卡久久| 波多野结衣巨乳人妻| 亚洲性夜色夜夜综合| 国产精品99久久久久久久久| 嫩草影院新地址| 欧美日韩国产亚洲二区| 观看美女的网站| 国产精品一区二区三区四区免费观看 | 精品无人区乱码1区二区| 亚洲激情五月婷婷啪啪| aaaaa片日本免费| 久久人人精品亚洲av| 欧美成人免费av一区二区三区| 草草在线视频免费看| 直男gayav资源| 亚洲一级一片aⅴ在线观看| 又爽又黄a免费视频| 精品乱码久久久久久99久播| 精品久久久久久久久久免费视频| 免费电影在线观看免费观看| 99久久成人亚洲精品观看| 久久中文看片网| 国产精品一区二区免费欧美| 亚州av有码| 国产色婷婷99| 少妇的逼好多水| 在线国产一区二区在线| 中文字幕熟女人妻在线| 午夜精品国产一区二区电影 | 国产亚洲91精品色在线| 一级毛片aaaaaa免费看小| 男女边吃奶边做爰视频| 18+在线观看网站| 嫩草影院精品99| 国产午夜精品论理片| 亚洲aⅴ乱码一区二区在线播放| 欧美最新免费一区二区三区| 午夜老司机福利剧场| 人妻夜夜爽99麻豆av| 精品欧美国产一区二区三| 国产午夜福利久久久久久| 一本一本综合久久| 国产高清视频在线播放一区| 能在线免费观看的黄片| 大型黄色视频在线免费观看| 亚洲成a人片在线一区二区| 欧美日韩精品成人综合77777| 蜜桃久久精品国产亚洲av| 精品一区二区免费观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲av成人av| 国产精品野战在线观看| 亚洲精品国产成人久久av| 欧美日韩国产亚洲二区| 两个人视频免费观看高清| 色综合色国产| 嫩草影院入口| 国产亚洲精品久久久久久毛片| 男人和女人高潮做爰伦理| 国产成人精品久久久久久| 免费观看的影片在线观看| 国产高潮美女av| 热99re8久久精品国产| 最好的美女福利视频网| 夜夜夜夜夜久久久久| 精华霜和精华液先用哪个| 色尼玛亚洲综合影院| 十八禁网站免费在线| 国产精品电影一区二区三区| 国产一区二区激情短视频| 18+在线观看网站| 一夜夜www| 欧美+日韩+精品| 精品久久久久久久末码| 国产精品久久久久久久电影| 好男人在线观看高清免费视频| 一区二区三区高清视频在线| 午夜激情欧美在线| 亚洲精品一区av在线观看| 校园春色视频在线观看| 97超视频在线观看视频| 天堂av国产一区二区熟女人妻| 日韩强制内射视频| 天天一区二区日本电影三级| 欧美成人a在线观看| 不卡一级毛片| 秋霞在线观看毛片| 搡老岳熟女国产| 亚洲图色成人| 有码 亚洲区| 两性午夜刺激爽爽歪歪视频在线观看| 日韩av在线大香蕉| 国产熟女欧美一区二区| 麻豆国产av国片精品| 国产精品无大码| 国产大屁股一区二区在线视频| 亚洲精品国产成人久久av| 国产精品国产高清国产av| 日本一二三区视频观看| 久久久精品94久久精品| 一个人看视频在线观看www免费| 亚洲高清免费不卡视频| 精品一区二区三区人妻视频|