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

    應用于類球體細胞電旋轉分析中的介電參數(shù)求解算法

    2012-09-17 06:57:06朱曉璐倪中華
    東南大學學報(自然科學版) 2012年2期
    關鍵詞:細胞質介電常數(shù)細胞膜

    張 麗 朱曉璐 倪中華

    (東南大學機械工程學院,南京 211189)

    應用于類球體細胞電旋轉分析中的介電參數(shù)求解算法

    張 麗 朱曉璐 倪中華

    (東南大學機械工程學院,南京 211189)

    根據(jù)薄膜理論建立了一種面向類球體細胞的單殼介電模型,分析了細胞在電旋轉靜流場中的受力,并由旋轉力矩和斯托克斯流體阻力矩平衡原理,推導得到細胞各層介電參數(shù)方程組.利用引入壓縮因子的粒子群優(yōu)化算法求解細胞的介電參數(shù),并討論了壓縮因子的變化對求解結果的影響.以Granulocytes白細胞為例進行算法驗證,得到細胞膜的電導率和相對介電常數(shù)分別為7.25 μS/m和10.3,細胞質的電導率和相對介電常數(shù)分別為0.487 S/m和136,且壓縮因子的變化僅影響各組解的搜索概率,加速系數(shù)為2.015時所需解的搜索概率達到最高值.算法結果與白細胞介電參數(shù)的經(jīng)驗數(shù)據(jù)相吻合,表明該算法具有較高的精確性和穩(wěn)定性.

    介電參數(shù);粒子群優(yōu)化;電旋轉;白細胞

    電動旋轉介電泳簡稱電旋轉[1],其基本原理是粒子在旋轉電場中被極化產(chǎn)生偶極矩,并與電場相互作用而產(chǎn)生旋轉運動.使用電旋轉技術操縱、辨識和表征生物粒子的理論及應用已成為研究熱點[2-4].電旋轉技術易于集成和控制,因而適用于紅細胞、白細胞、酵母細胞等生物粒子的辨識、偵測和分離.已知細胞的介電參數(shù)是識別和操控細胞的前提[5],有利于研究和分析細胞成分[6-7].細胞在旋轉電場中的旋轉角速度與細胞的介電參數(shù)有關,因此,若已知細胞的電旋轉定量運動規(guī)律,則可根據(jù)旋轉角速度與介電參數(shù)的數(shù)值關系,得到細胞的介電參數(shù)值.國內(nèi)外關于生物粒子的電旋轉特性和介電參數(shù)估計的研究很多.Takashima[8]利用內(nèi)部電極法和間隙法直接測量出較大生物細胞的介電性質;Asami等[9]提出應用微小吸管法來測定膜電容和膜電導;Mishima 等[6-7]發(fā)現(xiàn)相對介電常數(shù)增量與細胞內(nèi)部包含的物質成正比,以此估計細胞內(nèi)的生物質;Gascoyne等[5]采用靈敏度分析法和參數(shù)最優(yōu)化法來分析電旋轉數(shù)據(jù),得到介電參數(shù).

    研究電旋轉特性并間接估計細胞介電參數(shù)的工作,目前僅停留在細胞模型的建立[10]、公式推導[5,10]、測量裝置的改進[11-12]等方面,對采用特定的數(shù)學方法進行參數(shù)求解的研究甚少.本文針對類球體生物細胞的電旋轉數(shù)據(jù)特征,將粒子群優(yōu)化方法引入到類球體細胞介電參數(shù)的求解過程中,提出了一種以改進的粒子群優(yōu)化算法原理為基礎的介電參數(shù)求解算法,并以Granulocytes白細胞為例,分析了計算結果.本文算法充分運用了粒子群優(yōu)化的優(yōu)勢,大大提高了求解效率,同時降低了求解介電參數(shù)非線性方程組的難度.

    1 類球體細胞模型與參數(shù)

    生物細胞主要由細胞膜及內(nèi)部介質(簡稱細胞質)構成.假設細胞質內(nèi)部結構均勻,且類球體細胞的2層結構同心(見圖1(a)).根據(jù)Jones的薄膜理論[10],可將其等效為一個連續(xù)球形粒子(見圖1(b)).等效后的細胞復介電常數(shù)為[10]

    圖1 類球體細胞的等效模型

    式中,R1和R2分別為細胞和細胞質的半徑;ε*1,ε*2,ε*3分別為流體介質、細胞膜、細胞質的復介電常數(shù),且 ε*1= ε1- jσ1/ω,ε*2= ε2- jσ2/ω,ε*3=ε3- jσ3/ω,其中 ε1,ε2,ε3分別為流體介質、細胞膜和細胞質的介電常數(shù),σ1,σ2,σ3分別為流體介質、細胞膜和細胞質的電導率,ω為信號的角頻率.

    若細胞為哺乳動物細胞,且細胞膜的厚度遠小于細胞半徑,則其復介電常數(shù)可等效為[5]

    式中,C*m為細胞膜的膜電容,且C*m=ε2/Δ -jσ2/(Δω),其中Δ為細胞膜厚度.

    本文中,介電參數(shù)特指細胞膜電導率σ2、細胞膜介電常數(shù)ε2、細胞質電導率σ3和細胞質介電常數(shù) ε3.

    2 求解原理及方法

    2.1 介電參數(shù)方程組的建立

    將細胞放置于四相旋轉電極之間,通過向4個電極分別施加相位相差90°、幅值和頻率相等的信號,使細胞在偶極矩和旋轉電場的相互作用下發(fā)生旋轉(見圖2).

    圖2 電旋轉電極示意圖

    在介電泳引起的旋轉力矩及流體阻力矩的作用下,細胞勻速旋轉且保持平衡.旋轉力矩可表示為[10]

    式中,K(ω)為克勞修斯-莫索提因子,且;E為電場強度均方根值.

    由流體介質黏度產(chǎn)生的斯托克斯阻力矩為[5]

    式中,η為流體的黏度;Ω為旋轉角速度.

    根據(jù)力矩平衡原理可得TDEP=TStoke.假定細胞介電參數(shù)在靜流場環(huán)境中為定值,在實驗中測得細胞旋轉角速度Ω的頻率譜,已知流體介電參數(shù)及電場規(guī)律,可建立如下的n(n>4)維方程組:

    式中,ω1,ω2,…,ωn為信號角頻率;Ω1,Ω2,…,Ωn為相應的角速度.通過解此方程組,可獲得細胞的介電參數(shù).

    將式(1)代入式(5),可得超定非線性方程組,應用Matlab軟件對其進行直接求解難度高且精度低.本文采用一種基于迭代的進化算法——粒子群優(yōu)化(particle swarm optimization,PSO)算法[13]來求解方程組的數(shù)值解.與通用性較高的遺傳算法相比,PSO算法具有更強的全局優(yōu)化能力,能較快收斂,不需要考慮方程的具體組成形式,且對迭代初始點和目標函數(shù)導數(shù)的依賴性不大,適用于求解范圍未知的細胞介電參數(shù).

    2.2 粒子群優(yōu)化算法

    PSO算法是基于群體的演化算法,源于對鳥群捕食行為的研究[13].用PSO算法求解時,方程組介電參數(shù)的迭代解對應于鳥的位置,鳥被稱為粒子.每個粒子都具有位置和速度.速度決定粒子下一次迭代的方向和前進大小.粒子位置的好壞用適應度來判定,例如對于鳥群而言,可將鳥本身與食物的距離作為捕食行動的適應度,當某一適應度為0或達到一定極限值時則捕食成功.

    假設在一個D維搜索空間中,維數(shù)即為方程組未知數(shù)的個數(shù),隨機產(chǎn)生n個粒子(n也稱為群體規(guī)模),記第i個粒子的位置為 Xi={xi1,xi2,…,xiD},速度為 Vi={vi1,vi2,…,viD},其中i=1,2,…,n.將第i個粒子所經(jīng)歷的最佳位置稱為個體極值,記為Pi={pi1,pi2,…,piD};群體中所有粒子經(jīng)歷的最佳位置稱為全局極值,記為Gi={g1,g2,…,gD}.粒子的速度和位置可根據(jù)以下方程迭代獲得[13]:

    式中為第k次迭代前粒子的速度與位置;為第k次迭代后粒子的速度與位置;c1,c2為加速系數(shù),取值范圍為[0,2];r1,r2為隨機數(shù),取值范圍為[0,1].

    引入壓縮因子[13]可有效地搜索不同區(qū)域,從而解決各介電參數(shù)取值范圍數(shù)量級跨度大的問題,以確保PSO算法能以較快的速度收斂且得到高質量的解.相應的速度及位置迭代公式分別為

    式中,χ為壓縮因子,且

    令欲求解的非線性方程組的一般形式為

    根據(jù)PSO算法的通常處理辦法[14],定義第i個粒子的適應度函數(shù)為

    式中,qj(Xi)為第j個方程代入第i個粒子后得到的實值函數(shù).由此便可將非線性方程組的求解問題轉化為minF(Xi)的最優(yōu)化問題.

    將PSO壓縮因子法應用于求解方程組(5),令搜索粒子的位置 Xi={σi2,εi2,σi3,εi3},飛行速度 Vi={γi1,γi2,γi3,γi4}.根據(jù)式(7)可知,第i個粒子所對應的細胞膜電導率σi2和速度γi1的第k次迭代公式為將方程組(5)代入式(8),得到PSO算法求解類球體細胞介電參數(shù)方程組的適應度函數(shù),即

    采用PSO算法求解類球體細胞介電參數(shù)的非線性方程組的步驟如下:

    ①確定所求介電參數(shù)的個數(shù)(即粒子維數(shù))、粒子個數(shù)及各介電參數(shù)的數(shù)值范圍.

    ②根據(jù)這些變量的范圍,對每個粒子的介電參數(shù)(即位置)及速度進行初始化,并將每個粒子的初始介電參數(shù)作為該個體的個體極值.

    ③計算第i個粒子的適應度F(Xi),并與歷史最佳個體適應度以及全局最佳適應度進行比較.若F(Xi)小于歷史最佳個體適應度,則將第i個粒子的介電參數(shù)作為該粒子的個體極值.同理,與全局最佳適應度比較,若F(Xi)較小,則將當前的介電參數(shù)作為全局極值.

    ④根據(jù)式(9),更新粒子的介電參數(shù)和速度.

    ⑤如果滿足結束要求,則停止迭代;否則,返回步驟③.

    利用PSO算法求解非線性方程組的原理表明,當適應度越趨近于0時,得到的解越精確.為了提高算法的精確性,改進了適應度計算公式,對公式中的每一項進行1/10乘方,以避免計算時的舍入誤差對適應度值產(chǎn)生影響.假設原適應度計算式處理后為F(X).圖 3 為迭代次數(shù)N=1 000時的全局極值適應度曲線.由圖可知,N<50時適應度值下降幅度較大,N>50時下降幅度逐漸減小,且N>300時適應度值基本穩(wěn)定不變.根據(jù)迭代法收斂性的定義,可判斷該算法的收斂性:對于任意ξ>0,該算法存在一個迭代次數(shù)K,使得對于任意大于K的迭代次數(shù)h,滿足故由全局極值適應度構成的序列收斂.綜上可知,采用PSO求解類球體細胞介電參數(shù)方程組的算法收斂,且收斂速度滿足高效率計算的要求.

    圖3 全局極值適應度曲線

    3 計算結果與分析

    以類球體細胞Granulocytes白細胞為例,測試基于PSO的介電參數(shù)求解算法.電極形式為旋轉電極;流體黏度為 0.85 mPa·s,電導率為56 mS/m,介電常數(shù)為80.2;類球體細胞的細胞半徑為4.7 μm,細胞膜厚度為7.5 nm.Granulocytes白細胞的旋轉角速度與頻率的實驗數(shù)據(jù)見圖4.

    設定迭代次數(shù)為1 000,粒子群體規(guī)模為20.運行算法50次,分別統(tǒng)計細胞膜電導率σ2、細胞膜介電常數(shù)ε2、細胞質電導率σ3以及細胞質介電常數(shù)ε3的計算結果.需要說明的是,以下統(tǒng)計的介電常數(shù)為相對介電常數(shù),且各介電參數(shù)為平均值.壓縮因子的取值直接影響了搜索解的能力.試湊法得到的計算結果表明,當加速系數(shù)c1,c2在2.00~2.02內(nèi)取值時,方程組能獲得質量較高的解.表1為算法運算50次后介電參數(shù)的統(tǒng)計數(shù)據(jù).由表可知,運用本文算法可搜索到2組解,且壓縮因子的取值變化僅影響搜到解的概率,不影響解的數(shù)值.

    圖4 Granulocytes白細胞的旋轉角速度-頻率實驗曲線[4]

    表1 介電參數(shù)的統(tǒng)計結果

    對于單殼模型而言,典型的生物粒子介電參數(shù)范圍為[5]:5 mF/m2<Cm<35 mF/m2,20 S/m2<Gm<2 000 S/m2,0.1 S/m <σ3<1.3 S/m,50ε0<ε3<150ε0,其中Cm= ε2/Δ 為膜電容,Gm= σ2/Δ為膜電導,ε0為真空介電系數(shù).表1所列2組解均滿足上述范圍.

    為了進一步驗證本文算法的穩(wěn)定性,考察了c1,c2取值不同時算法運算50次2組解的分布情況(見圖5).由圖5(a)可知,細胞膜電導率主要分布在4~10 μS/m范圍內(nèi),σ2的2組解范圍有交集;細胞膜介電常數(shù)及細胞質電導率分布較聚集,且基本不存在交叉的現(xiàn)象;細胞質介電常數(shù)的第1組解分布集中,第2組解分布跨度較大,但2組數(shù)值沒有交叉.比較4個介電參數(shù)發(fā)現(xiàn),細胞膜電導率和細胞質介電常數(shù)的分布跨度較大,由此可推測,實驗中這2個參數(shù)的取值在一定范圍內(nèi)對類球體細胞的速度影響不大.相比之下,細胞膜介電常數(shù)和細胞質電導率更能體現(xiàn)細胞的介電特性.對比圖5(b)和(c)發(fā)現(xiàn),壓縮因子的取值僅影響各組解的聚集情況及搜索解的概率.

    根據(jù)數(shù)值分析的理論定義,一個算法的輸入數(shù)據(jù)有擾動,但計算過程中誤差不增長,則稱此算法是數(shù)值穩(wěn)定的.在本文算法中,c1,c2取值變化時計算結果基本不變.因此,運用 PSO求解 Granulocytes白細胞介電參數(shù)的算法是數(shù)值穩(wěn)定的.

    圖5 2組解的分布情況

    文獻[4]列出了Granulocytes白細胞介電參數(shù)的經(jīng)驗數(shù)據(jù),即細胞半徑為(4.71±0.23)μm,膜電容為(11.0±3.2)mF/m2,細胞質電導率為(0.60±0.13)S/m,細胞質介電常數(shù)為150.9±39.3.將這些數(shù)據(jù)作為判定本文算法正確性的依據(jù).經(jīng)計算可得ε2為9.32±2.71,比較計算結果可知,第1組解符合經(jīng)驗數(shù)據(jù)的取值范圍,第2組解中細胞質的介電常數(shù)低于經(jīng)驗數(shù)據(jù).第1組解的波動范圍較小,取值較為集中,若以序列來表征50次計算的結果,則第1組解的方差顯然小于第2組解,故第1組解更精密.綜上可知,基于PSO的介電參數(shù)求解算法解得Granulocytes白細胞介電參數(shù)的最優(yōu)解為第1組解.當c1=c2=2.015時,取搜索概率較高的第1組解,則細胞膜電導率為7.25 μS/m,細胞膜相對介電常數(shù)為10.3,細胞質電導率為0.487 S/m,細胞質相對介電常數(shù)為136,符合經(jīng)驗數(shù)據(jù).因此,運用PSO求解Granulocytes白細胞介電參數(shù)的算法是精確的.

    4 結語

    本文描述一種基于PSO的類球體細胞介電參數(shù)求解算法.研究結果表明,引入壓縮因子的PSO介電參數(shù)求解算法具有較好的收斂性,計算效率高,不依賴介電參數(shù)方程組形式且靈活性大,通過選擇合適的壓縮因子可獲得高質量的解.以類球體細胞Granulocytes白細胞為例,在電旋轉電極形式下,采用PSO算法求解介電參數(shù)方程組,可成功得到符合白細胞經(jīng)驗數(shù)據(jù)的計算解,即細胞膜電導率為7.25 μS/m,細胞膜相對介電常數(shù)為10.3,細胞質電導率為0.487 S/m,細胞質相對介電常數(shù)為136.因此,基于PSO的介電參數(shù)求解算法可有效地估算出類球體細胞的介電參數(shù),這對研究未知生物粒子的性能具有重要意義.如何改進算法以及將其應用于橢球體等其他形狀的生物細胞是進一步研究的重點之一.

    [1] Arnold W M,Zimmermann U.Electro-rotation:development of a technique for dielectric measurements on individual cells and particles[J].Journal of Electrostatics,1988,21(2/3):151-191.

    [2] Sancho M,Martínez G,Mu?oz S,et al.Interaction between cells in dielectrophoresis and electrorotation experiments[J].Biomicrofluidics,2010,4(2):1-11.

    [3] Voyer D,F(xiàn)rénéa-Robin M,Buret F,et al.Improvements in the extraction of cell electric properties from their electrorotation spectrum [J].Bioelectrochemistry,2010,79(1):25-30.

    [4] Yang J,Huang Y,Wang X,et al.Dielectric properties of human leukocyte subpopulations determined by electrorotation as a cell separation criterion[J].Biophysical Journal,1999,76(6):3307-3314.

    [5]Gascoyne P R C,Becker F F,Wang X B.Numerical analysis of the influence of experimental conditions on the accuracy of dielectric parameters derived from electrorotation measurements[J].Bioelectrochemistry and Bioenergetics,1995,36(2):115-125.

    [6] Mishima B,Mimura K A,Takahara Y,et al.On-line monitoring of cell concentrations by dielectric measurements[J].Journal of Fermentation and Bioengineer-ing,1991,72(4):291-295.

    [7]Harris C M,Todd R W,Bungard S J,et al.Dielectric permittivity of microbial suspensions at radio frequencies:a novel method for the real-time estimation of microbial biomass[J].Enzyme Microb Technol,1987,9(3):181-186.

    [8] Takashima S.Passive electrical properties and voltage dependent membrane capacitance of single skeletal muscle fibers[J].Pflugers Arch Eur J Physiol,1985,403(2):197-204.

    [9] Asami K,Takashima S.Membrane admittance of cloned muscle cells in culture:use of a micropipette technique[J].Biochim Biophys Acta,1994,1190(1):129-136.

    [10] Jones T B.Basic theory of dielectrophoresis and electrorotation [J].IEEE Engineering in Medicine and Biology Magazine,2003,22(6):33-42.

    [11]Gawad S,Cheung K,Seger U,et al.Dielectric spectroscopy in a micromachined flow cytometer:theoretical and practical considerations[J].Lab Chip,2004,4(3):241-251.

    [12]Holmes D,Green N G,Morgan H.Microdevices for dielectrophoretic flow-through cell separation [J].IEEE Engineering in Medicine and Biology Magazine,2003,22(6):85-89.

    [13]楊維,李歧強.粒子群優(yōu)化算法綜述[J].中國工程科學,2004,6(5):87-94.

    Yang Wei,Li Qiqiang.Survey on particle swarm optimization algorithm [J].Engineering Sciences,2004,6(5):87-94.(in Chinese)

    [14]陳長憶,葉永春.基于粒子群算法的非線性方程組求解[J].計算機應用與軟件,2006,23(5):137-139.

    Chen Changyi,Ye Yongchun.Solving nonlinear systems of equations based on particle swarm optimization[J].Computer Applications and Software,2006,23(5):137-139.(in Chinese)

    Dielectric parameter solving algorithm in electrorotation analysis of spherical cells

    Zhang Li Zhu Xiaolu Ni Zhonghua
    (School of Mechanical Engineering,Southeast University,Nanjing 211189,China)

    A dielectric model with single concentric shell of spherical cells is established based on film theory.The forces of the cells in a static flow field under the action of electrorotation are analyzed,and the force equilibrium equations are derived to obtain dielectric parameters of each shell by the balance of the rotating torque and the Stokes fluid drag torque.Then,the particle swarm optimization algorithm improved with the constriction factor is used to solve the equations.The influence of the change in the constriction factor on the solution results is also discussed.This algorithm is validated with Granulocytes as an example.The calculated membrane conductivity,membrane relative permittivity,cytoplasm conductivity and cytoplasm relative permittivity are 7.25 μS/m,10.3,0.487 S/m and 136,respectively.Besides,the change in the constriction factor just affects the search probability of each solution,and the search probability needed gets the highest value with the acceleration factor of 2.015.The calculated results are consistent with the empirical data,showing that the proposed algorithm can obtain high accuracy and stability.

    dielectric parameter;particle swarm optimization;electrorotation;leukocyte

    TP391

    A

    1001-0505(2012)02-0280-06

    10.3969/j.issn.1001-0505.2012.02.017

    2011-09-25.

    張麗(1987—),女,碩士生;倪中華(聯(lián)系人),男,博士,教授,博士生導師,nzh2003@seu.edu.cn.

    國家自然科學基金重大研究計劃(培育)資助項目(91023024)、教育部博士研究生學術新人獎資助項目、江蘇省普通高校研究生科研創(chuàng)新計劃資助項目(CX10B_062Z)、東南大學優(yōu)秀博士學位論文基金資助項目.

    張麗,朱曉璐,倪中華.應用于類球體細胞電旋轉分析中的介電參數(shù)求解算法[J].東南大學學報:自然科學版,2012,42(2):280-285.[doi:10.3969/j.issn.1001-0505.2012.02.017]

    猜你喜歡
    細胞質介電常數(shù)細胞膜
    作物細胞質雄性不育系實現(xiàn)快速創(chuàng)制
    科學導報(2024年20期)2024-04-22 09:54:13
    無鉛Y5U103高介電常數(shù)瓷料研究
    電子制作(2017年20期)2017-04-26 06:57:40
    低介電常數(shù)聚酰亞胺基多孔復合材料的研究進展
    節(jié)水抗旱細胞質雄性不育系滬旱7A 的選育與利用
    洋蔥細胞質雄性不育基因分子標記研究進展
    中國蔬菜(2015年9期)2015-12-21 13:04:38
    壇紫菜細胞質型果糖1,6-二磷酸酶基因的克隆及表達分析
    低介電常數(shù)聚酰亞胺薄膜研究進展
    中國塑料(2015年8期)2015-10-14 01:10:40
    皮膚磨削術聯(lián)合表皮細胞膜片治療穩(wěn)定期白癜風療效觀察
    宮永寬:給生物醫(yī)用材料穿上仿細胞膜外衣
    香芹酚對大腸桿菌和金黃色葡萄球菌細胞膜的影響
    国产三级中文精品| 搡女人真爽免费视频火全软件| 久久精品人妻少妇| 3wmmmm亚洲av在线观看| 国产成人免费观看mmmm| 能在线免费看毛片的网站| 亚洲国产精品成人久久小说| 成年av动漫网址| 老师上课跳d突然被开到最大视频| 国产欧美日韩精品一区二区| 国产精品国产高清国产av| 三级国产精品欧美在线观看| 又粗又爽又猛毛片免费看| av免费观看日本| 在线免费观看的www视频| 亚洲欧洲日产国产| 精品久久久久久久末码| 免费看日本二区| 最近最新中文字幕大全电影3| 成人一区二区视频在线观看| 国产探花极品一区二区| 日本-黄色视频高清免费观看| 亚洲精品自拍成人| 久久亚洲精品不卡| 久久99热这里只有精品18| 长腿黑丝高跟| 国产亚洲91精品色在线| 亚洲av男天堂| 日本午夜av视频| 亚洲乱码一区二区免费版| 2022亚洲国产成人精品| 神马国产精品三级电影在线观看| 日本与韩国留学比较| 欧美日韩在线观看h| 国产精品福利在线免费观看| 日韩一本色道免费dvd| 午夜福利在线在线| 成人三级黄色视频| 久久精品久久久久久噜噜老黄 | 麻豆成人午夜福利视频| 一区二区三区四区激情视频| 观看美女的网站| 国产黄片视频在线免费观看| 欧美成人a在线观看| 观看美女的网站| 欧美性猛交╳xxx乱大交人| 嫩草影院新地址| 免费观看性生交大片5| 亚洲欧美成人综合另类久久久 | 亚洲国产欧美在线一区| 久久精品影院6| av视频在线观看入口| av在线天堂中文字幕| 亚洲国产欧美在线一区| 国内揄拍国产精品人妻在线| 亚洲精品影视一区二区三区av| 久久久久久大精品| 亚洲18禁久久av| 自拍偷自拍亚洲精品老妇| 久久久久久九九精品二区国产| 熟女电影av网| av福利片在线观看| 欧美最新免费一区二区三区| 国产探花极品一区二区| 看免费成人av毛片| 建设人人有责人人尽责人人享有的 | 97超视频在线观看视频| 亚洲真实伦在线观看| 网址你懂的国产日韩在线| 一级毛片aaaaaa免费看小| 国产 一区 欧美 日韩| 晚上一个人看的免费电影| 热99re8久久精品国产| 日本欧美国产在线视频| 亚洲va在线va天堂va国产| 最后的刺客免费高清国语| 天天躁日日操中文字幕| 男女视频在线观看网站免费| 免费搜索国产男女视频| 免费观看的影片在线观看| 午夜亚洲福利在线播放| 亚洲最大成人av| АⅤ资源中文在线天堂| 一区二区三区高清视频在线| 插逼视频在线观看| 99久久九九国产精品国产免费| 欧美不卡视频在线免费观看| 亚洲色图av天堂| 国产精品一区二区三区四区免费观看| 成年版毛片免费区| 国产精品一区二区三区四区免费观看| 久久99精品国语久久久| 亚洲内射少妇av| 国产精品嫩草影院av在线观看| 乱码一卡2卡4卡精品| 国产探花极品一区二区| 边亲边吃奶的免费视频| 3wmmmm亚洲av在线观看| 亚洲av成人av| 一夜夜www| 国产精华一区二区三区| 免费看av在线观看网站| 天堂影院成人在线观看| 亚洲欧美成人精品一区二区| 色网站视频免费| 国产午夜精品久久久久久一区二区三区| 亚洲自拍偷在线| 亚洲国产高清在线一区二区三| 国产熟女欧美一区二区| 国产三级在线视频| 听说在线观看完整版免费高清| 成人漫画全彩无遮挡| 91精品伊人久久大香线蕉| 毛片女人毛片| 老师上课跳d突然被开到最大视频| 亚洲,欧美,日韩| 一区二区三区高清视频在线| 国内揄拍国产精品人妻在线| 日韩av不卡免费在线播放| 亚洲av二区三区四区| 亚洲精品乱码久久久久久按摩| 寂寞人妻少妇视频99o| 国产午夜精品久久久久久一区二区三区| 欧美精品一区二区大全| 亚洲乱码一区二区免费版| 婷婷色综合大香蕉| 日韩中字成人| 国产一区二区在线观看日韩| www日本黄色视频网| 97人妻精品一区二区三区麻豆| 91aial.com中文字幕在线观看| 久久国产乱子免费精品| 97人妻精品一区二区三区麻豆| 欧美日韩在线观看h| 成人三级黄色视频| 久久婷婷人人爽人人干人人爱| 精华霜和精华液先用哪个| 精品国产三级普通话版| 美女xxoo啪啪120秒动态图| 日本欧美国产在线视频| 亚洲国产精品久久男人天堂| 三级经典国产精品| 亚洲美女视频黄频| 性色avwww在线观看| 简卡轻食公司| 18禁在线无遮挡免费观看视频| 大香蕉97超碰在线| 精品无人区乱码1区二区| 国产v大片淫在线免费观看| 国产久久久一区二区三区| 最近最新中文字幕大全电影3| 少妇的逼水好多| 亚洲最大成人中文| 日本黄色视频三级网站网址| 午夜激情福利司机影院| 麻豆一二三区av精品| 国产精品1区2区在线观看.| 91久久精品电影网| 欧美日本亚洲视频在线播放| 亚州av有码| 欧美一区二区国产精品久久精品| av国产久精品久网站免费入址| 亚洲av免费高清在线观看| 晚上一个人看的免费电影| 中文字幕亚洲精品专区| 国产麻豆成人av免费视频| 亚洲四区av| 美女cb高潮喷水在线观看| 女人十人毛片免费观看3o分钟| 美女xxoo啪啪120秒动态图| 我要看日韩黄色一级片| 日韩av在线大香蕉| 一级毛片我不卡| 99热全是精品| 亚洲av成人精品一区久久| 成人欧美大片| 亚洲经典国产精华液单| 一区二区三区高清视频在线| 亚洲美女视频黄频| 最近视频中文字幕2019在线8| 亚洲国产精品专区欧美| 又粗又硬又长又爽又黄的视频| 丰满乱子伦码专区| 禁无遮挡网站| 麻豆精品久久久久久蜜桃| 乱码一卡2卡4卡精品| 中文字幕久久专区| 免费看光身美女| 日韩精品青青久久久久久| 最近最新中文字幕免费大全7| 成人特级av手机在线观看| 久久午夜福利片| 国产视频内射| 欧美成人免费av一区二区三区| 啦啦啦观看免费观看视频高清| 国产伦一二天堂av在线观看| 欧美人与善性xxx| 日本免费在线观看一区| 2022亚洲国产成人精品| 91av网一区二区| 嫩草影院新地址| 如何舔出高潮| 国产高清不卡午夜福利| 精品久久久久久久久av| 日韩 亚洲 欧美在线| 国产在线男女| 2022亚洲国产成人精品| 成人漫画全彩无遮挡| 免费看av在线观看网站| 一级黄色大片毛片| 国产成人免费观看mmmm| 国产精华一区二区三区| 精品国内亚洲2022精品成人| av在线蜜桃| 精品久久久久久久久av| 非洲黑人性xxxx精品又粗又长| 99久久精品一区二区三区| 午夜福利高清视频| 欧美性猛交╳xxx乱大交人| 欧美3d第一页| 国产黄片美女视频| 丰满人妻一区二区三区视频av| 久久国内精品自在自线图片| 99久久成人亚洲精品观看| 久久精品国产99精品国产亚洲性色| 国产探花极品一区二区| 国产三级在线视频| 欧美zozozo另类| 久热久热在线精品观看| 中文亚洲av片在线观看爽| 美女黄网站色视频| 国产中年淑女户外野战色| 国产又色又爽无遮挡免| 日韩欧美精品v在线| 一个人观看的视频www高清免费观看| 晚上一个人看的免费电影| 看十八女毛片水多多多| 91午夜精品亚洲一区二区三区| 99热网站在线观看| 国产精品久久久久久精品电影小说 | 国模一区二区三区四区视频| av在线老鸭窝| 欧美精品一区二区大全| 日本免费a在线| 精品久久久久久成人av| 亚洲中文字幕日韩| 成人特级av手机在线观看| 韩国高清视频一区二区三区| 天天一区二区日本电影三级| 久久久成人免费电影| 亚洲va在线va天堂va国产| 久久99热6这里只有精品| 尾随美女入室| 婷婷色综合大香蕉| 天天躁日日操中文字幕| 成人欧美大片| 日韩中字成人| 国产片特级美女逼逼视频| 精品99又大又爽又粗少妇毛片| 精品不卡国产一区二区三区| 日韩三级伦理在线观看| 99热全是精品| 老司机影院毛片| 亚洲18禁久久av| 国产一区亚洲一区在线观看| 国产成人免费观看mmmm| 网址你懂的国产日韩在线| 国产精品不卡视频一区二区| 亚洲av成人精品一区久久| 联通29元200g的流量卡| 内射极品少妇av片p| 全区人妻精品视频| 国产精品久久久久久久电影| 蜜臀久久99精品久久宅男| 久久99精品国语久久久| 国产午夜精品久久久久久一区二区三区| 国产精品99久久久久久久久| 黄色配什么色好看| 成人毛片60女人毛片免费| 久久久午夜欧美精品| 热99在线观看视频| 激情 狠狠 欧美| 2021少妇久久久久久久久久久| 精品无人区乱码1区二区| 直男gayav资源| 亚洲av中文字字幕乱码综合| 国内精品一区二区在线观看| 人妻夜夜爽99麻豆av| 亚洲美女视频黄频| 嘟嘟电影网在线观看| 日韩人妻高清精品专区| 国内精品宾馆在线| 色网站视频免费| 欧美性感艳星| 麻豆一二三区av精品| 久久久久久久久大av| 麻豆av噜噜一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 超碰av人人做人人爽久久| 美女被艹到高潮喷水动态| 国产成人91sexporn| 国产精品国产三级国产专区5o | 搞女人的毛片| 欧美一区二区国产精品久久精品| 午夜精品一区二区三区免费看| 亚洲av二区三区四区| 波野结衣二区三区在线| 日韩欧美国产在线观看| 国产精品不卡视频一区二区| 看免费成人av毛片| 夜夜看夜夜爽夜夜摸| 身体一侧抽搐| 成人欧美大片| 看黄色毛片网站| 99久久精品一区二区三区| 久久综合国产亚洲精品| 久久鲁丝午夜福利片| 免费电影在线观看免费观看| 91在线精品国自产拍蜜月| 欧美丝袜亚洲另类| 欧美潮喷喷水| 精品人妻偷拍中文字幕| 97在线视频观看| 国产精品一二三区在线看| 久久人妻av系列| 午夜福利在线观看免费完整高清在| АⅤ资源中文在线天堂| 日韩 亚洲 欧美在线| 久久久精品大字幕| ponron亚洲| av免费在线看不卡| 51国产日韩欧美| 在线播放无遮挡| 日韩欧美精品v在线| 日韩成人av中文字幕在线观看| 干丝袜人妻中文字幕| 日日摸夜夜添夜夜添av毛片| 黄片无遮挡物在线观看| 69人妻影院| 久久精品国产亚洲网站| 中文亚洲av片在线观看爽| 亚洲自偷自拍三级| 国产亚洲精品久久久com| 国产亚洲av片在线观看秒播厂 | 1024手机看黄色片| 国产亚洲av片在线观看秒播厂 | 99热6这里只有精品| 国产精品一区二区三区四区免费观看| 99热网站在线观看| 好男人视频免费观看在线| videos熟女内射| 中文欧美无线码| 搡女人真爽免费视频火全软件| 欧美成人a在线观看| 国产精品久久久久久久久免| 国产免费一级a男人的天堂| 五月玫瑰六月丁香| 午夜爱爱视频在线播放| 国产一区有黄有色的免费视频 | 成人欧美大片| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲色图av天堂| 美女xxoo啪啪120秒动态图| 国产精品久久久久久精品电影小说 | 18禁在线播放成人免费| 99久国产av精品| 亚洲av日韩在线播放| 综合色丁香网| 99热这里只有精品一区| 国内精品美女久久久久久| 成人特级av手机在线观看| 国产精品伦人一区二区| 欧美+日韩+精品| 亚洲自偷自拍三级| 免费看日本二区| 高清午夜精品一区二区三区| 日韩强制内射视频| 极品教师在线视频| 寂寞人妻少妇视频99o| 亚洲国产欧洲综合997久久,| 极品教师在线视频| 欧美三级亚洲精品| 最近视频中文字幕2019在线8| 欧美另类亚洲清纯唯美| 亚洲第一区二区三区不卡| 国产高清国产精品国产三级 | 国产在线男女| 精华霜和精华液先用哪个| 18禁裸乳无遮挡免费网站照片| 青春草视频在线免费观看| 亚洲精品乱码久久久久久按摩| 欧美zozozo另类| 天堂网av新在线| 亚洲精品日韩在线中文字幕| 久久99热这里只频精品6学生 | 欧美最新免费一区二区三区| 国产亚洲91精品色在线| 免费在线观看成人毛片| 少妇高潮的动态图| 亚洲精品乱久久久久久| 非洲黑人性xxxx精品又粗又长| 中文资源天堂在线| 欧美又色又爽又黄视频| 插阴视频在线观看视频| 欧美潮喷喷水| 美女大奶头视频| 一区二区三区乱码不卡18| 欧美变态另类bdsm刘玥| 美女大奶头视频| 男插女下体视频免费在线播放| 美女被艹到高潮喷水动态| av免费观看日本| 日韩欧美精品免费久久| 久久亚洲国产成人精品v| 国产精品久久久久久精品电影| 热99re8久久精品国产| 国内精品一区二区在线观看| 亚洲18禁久久av| 欧美极品一区二区三区四区| 日本黄大片高清| 99九九线精品视频在线观看视频| 亚洲久久久久久中文字幕| 身体一侧抽搐| 日韩高清综合在线| av在线亚洲专区| 丰满人妻一区二区三区视频av| 亚洲成人久久爱视频| 亚洲欧美日韩东京热| 美女黄网站色视频| 51国产日韩欧美| 黄色日韩在线| 色哟哟·www| 中文字幕免费在线视频6| 纵有疾风起免费观看全集完整版 | 美女大奶头视频| 乱系列少妇在线播放| 大又大粗又爽又黄少妇毛片口| 久久久久久九九精品二区国产| 校园人妻丝袜中文字幕| 午夜精品一区二区三区免费看| 男女下面进入的视频免费午夜| 久久久国产成人免费| 国产伦在线观看视频一区| 久久人人爽人人爽人人片va| 国产 一区 欧美 日韩| 免费观看a级毛片全部| 韩国av在线不卡| 午夜日本视频在线| 99久久精品国产国产毛片| 亚洲国产精品合色在线| АⅤ资源中文在线天堂| 亚洲精品一区蜜桃| 国产乱人视频| 一区二区三区乱码不卡18| 国产午夜精品论理片| 亚洲精品色激情综合| 高清在线视频一区二区三区 | 一级毛片久久久久久久久女| 午夜爱爱视频在线播放| 日韩欧美在线乱码| 亚洲av免费在线观看| 99热这里只有精品一区| 亚洲国产精品sss在线观看| 乱系列少妇在线播放| 大话2 男鬼变身卡| 能在线免费看毛片的网站| 国产精品一区二区三区四区免费观看| 偷拍熟女少妇极品色| 赤兔流量卡办理| 国产黄a三级三级三级人| 亚洲在线观看片| 丰满人妻一区二区三区视频av| 久久精品国产亚洲av天美| 国产精品福利在线免费观看| 国产成人aa在线观看| 久久久久久九九精品二区国产| 好男人在线观看高清免费视频| 日本免费在线观看一区| 久久人人爽人人片av| 精华霜和精华液先用哪个| 老女人水多毛片| 亚洲欧洲日产国产| 国产成人freesex在线| 97热精品久久久久久| 国产成人精品久久久久久| 中文字幕精品亚洲无线码一区| 国产真实伦视频高清在线观看| 亚洲精品乱码久久久v下载方式| 97超视频在线观看视频| 欧美高清成人免费视频www| 亚洲国产精品成人综合色| 寂寞人妻少妇视频99o| 韩国高清视频一区二区三区| 日本wwww免费看| 一级毛片aaaaaa免费看小| 青青草视频在线视频观看| 色网站视频免费| 狂野欧美白嫩少妇大欣赏| 欧美bdsm另类| 黄色欧美视频在线观看| 精品免费久久久久久久清纯| av.在线天堂| 三级毛片av免费| 久久亚洲精品不卡| av专区在线播放| 亚洲欧美日韩卡通动漫| 亚洲精品aⅴ在线观看| 国产69精品久久久久777片| 日日撸夜夜添| 国产精品不卡视频一区二区| 丰满少妇做爰视频| 又粗又爽又猛毛片免费看| 亚洲久久久久久中文字幕| 床上黄色一级片| 国产伦理片在线播放av一区| 成人无遮挡网站| 在线观看av片永久免费下载| 亚洲av不卡在线观看| 久久久久久久久久黄片| 久久亚洲国产成人精品v| 一边摸一边抽搐一进一小说| 亚洲人成网站在线观看播放| 欧美性感艳星| 夜夜看夜夜爽夜夜摸| 韩国高清视频一区二区三区| 国产精品国产三级国产av玫瑰| 日韩成人伦理影院| 午夜免费激情av| 波多野结衣高清无吗| 亚洲国产精品国产精品| 舔av片在线| 国产高清不卡午夜福利| 麻豆乱淫一区二区| 精品一区二区免费观看| 波多野结衣巨乳人妻| 丝袜美腿在线中文| 在线观看av片永久免费下载| 亚洲精品乱码久久久v下载方式| 国产一区二区亚洲精品在线观看| 欧美不卡视频在线免费观看| 最近2019中文字幕mv第一页| 最新中文字幕久久久久| 久久精品熟女亚洲av麻豆精品 | 六月丁香七月| 神马国产精品三级电影在线观看| 又爽又黄a免费视频| 最近视频中文字幕2019在线8| 最近中文字幕2019免费版| 日本一本二区三区精品| 美女xxoo啪啪120秒动态图| 日本黄色视频三级网站网址| 免费看日本二区| 长腿黑丝高跟| 小蜜桃在线观看免费完整版高清| 嫩草影院精品99| 欧美激情国产日韩精品一区| 两个人视频免费观看高清| 国产黄片美女视频| 亚洲高清免费不卡视频| 夫妻性生交免费视频一级片| 嫩草影院入口| 最近中文字幕2019免费版| 毛片女人毛片| 一区二区三区高清视频在线| 国产成人aa在线观看| 久久久午夜欧美精品| 桃色一区二区三区在线观看| 久久国产乱子免费精品| 亚洲国产精品久久男人天堂| 高清日韩中文字幕在线| 中文在线观看免费www的网站| 精品酒店卫生间| 美女内射精品一级片tv| 亚洲av一区综合| 国产精品不卡视频一区二区| 免费大片18禁| 国产伦一二天堂av在线观看| 日本免费a在线| 九草在线视频观看| 国模一区二区三区四区视频| 国内少妇人妻偷人精品xxx网站| 麻豆av噜噜一区二区三区| 日本wwww免费看| 欧美潮喷喷水| 观看美女的网站| 你懂的网址亚洲精品在线观看 | 日本黄色视频三级网站网址| 久久精品国产亚洲网站| 97人妻精品一区二区三区麻豆| 一级av片app| 国产成人免费观看mmmm| 国产午夜精品一二区理论片| 最近中文字幕高清免费大全6| 夜夜看夜夜爽夜夜摸| 18+在线观看网站| 一个人免费在线观看电影| 国产黄片美女视频| 亚洲av免费高清在线观看| 精品不卡国产一区二区三区| 2021少妇久久久久久久久久久| 22中文网久久字幕| 超碰97精品在线观看| 精品国产三级普通话版| 精品久久久噜噜| 成人漫画全彩无遮挡| 一个人看的www免费观看视频| 在现免费观看毛片| 国产三级中文精品| 18禁动态无遮挡网站| 丰满人妻一区二区三区视频av| 国产激情偷乱视频一区二区| 免费大片18禁| 美女脱内裤让男人舔精品视频| 亚洲欧美一区二区三区国产| 日韩中字成人|