• <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ī)用材料穿上仿細胞膜外衣
    香芹酚對大腸桿菌和金黃色葡萄球菌細胞膜的影響
    精品久久久精品久久久| 欧美一级毛片孕妇| 男女午夜视频在线观看| 97超级碰碰碰精品色视频在线观看| 丰满饥渴人妻一区二区三| 这个男人来自地球电影免费观看| 亚洲一区中文字幕在线| 欧美成人性av电影在线观看| 国产熟女xx| 中文字幕另类日韩欧美亚洲嫩草| 久久精品国产综合久久久| 久久久国产精品麻豆| 热re99久久国产66热| 美女高潮到喷水免费观看| 可以在线观看毛片的网站| 日本精品一区二区三区蜜桃| 怎么达到女性高潮| 19禁男女啪啪无遮挡网站| 国产精品野战在线观看 | 欧美中文日本在线观看视频| 国产在线观看jvid| 婷婷丁香在线五月| 国产精品香港三级国产av潘金莲| 免费久久久久久久精品成人欧美视频| 国产又色又爽无遮挡免费看| 欧美中文综合在线视频| 国产一卡二卡三卡精品| 怎么达到女性高潮| 国产成人精品在线电影| 国产亚洲精品一区二区www| 午夜亚洲福利在线播放| 亚洲国产精品合色在线| 免费看十八禁软件| 亚洲成人免费电影在线观看| 波多野结衣一区麻豆| 久久欧美精品欧美久久欧美| 十分钟在线观看高清视频www| 亚洲欧美日韩无卡精品| 国产精品免费视频内射| 麻豆av在线久日| 色哟哟哟哟哟哟| 国产麻豆69| 亚洲少妇的诱惑av| 免费不卡黄色视频| 777久久人妻少妇嫩草av网站| 丰满的人妻完整版| xxxhd国产人妻xxx| 日韩视频一区二区在线观看| 手机成人av网站| 桃红色精品国产亚洲av| 国产野战对白在线观看| 午夜福利欧美成人| 国产一区二区在线av高清观看| 精品福利观看| 中文欧美无线码| 成人国产一区最新在线观看| 麻豆一二三区av精品| 在线十欧美十亚洲十日本专区| 亚洲欧美日韩高清在线视频| 国产97色在线日韩免费| 视频在线观看一区二区三区| 日韩三级视频一区二区三区| 一边摸一边抽搐一进一出视频| 91国产中文字幕| 精品久久久久久,| xxx96com| 91av网站免费观看| 午夜两性在线视频| 性少妇av在线| 精品乱码久久久久久99久播| 国产熟女xx| 欧美精品一区二区免费开放| 午夜精品在线福利| 超色免费av| 久久人妻熟女aⅴ| 长腿黑丝高跟| 日韩 欧美 亚洲 中文字幕| 18禁观看日本| 亚洲av成人不卡在线观看播放网| 最好的美女福利视频网| 久久精品91蜜桃| 美女扒开内裤让男人捅视频| 黄色视频不卡| 亚洲五月色婷婷综合| 性少妇av在线| 欧美日韩亚洲国产一区二区在线观看| 成年人免费黄色播放视频| 老司机在亚洲福利影院| 日本wwww免费看| 国产一卡二卡三卡精品| 高清欧美精品videossex| 后天国语完整版免费观看| 亚洲欧美日韩无卡精品| 欧美最黄视频在线播放免费 | 丝袜人妻中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 人人妻人人添人人爽欧美一区卜| 国产伦人伦偷精品视频| 丰满饥渴人妻一区二区三| 精品午夜福利视频在线观看一区| 日本精品一区二区三区蜜桃| 欧美日韩av久久| aaaaa片日本免费| 欧美黄色片欧美黄色片| 欧美精品啪啪一区二区三区| 日韩有码中文字幕| 久久中文看片网| 欧美国产精品va在线观看不卡| 亚洲人成电影免费在线| 国产欧美日韩一区二区三| 日韩国内少妇激情av| 久久久精品国产亚洲av高清涩受| 91麻豆精品激情在线观看国产 | 国产精品一区二区三区四区久久 | 国产精品99久久99久久久不卡| 性少妇av在线| 在线天堂中文资源库| 国产精品98久久久久久宅男小说| 欧美丝袜亚洲另类 | 久久精品亚洲精品国产色婷小说| 亚洲色图 男人天堂 中文字幕| 一a级毛片在线观看| av电影中文网址| 99久久国产精品久久久| 亚洲男人天堂网一区| 成年版毛片免费区| 狂野欧美激情性xxxx| 淫秽高清视频在线观看| 12—13女人毛片做爰片一| 精品高清国产在线一区| 日本黄色视频三级网站网址| 亚洲少妇的诱惑av| 怎么达到女性高潮| 久久国产亚洲av麻豆专区| 亚洲一码二码三码区别大吗| 少妇 在线观看| 色哟哟哟哟哟哟| 亚洲色图av天堂| 男女午夜视频在线观看| 欧美亚洲日本最大视频资源| 亚洲少妇的诱惑av| 人人妻人人爽人人添夜夜欢视频| 国产一区二区三区在线臀色熟女 | 成人三级做爰电影| 看片在线看免费视频| 999久久久精品免费观看国产| 色在线成人网| 悠悠久久av| 欧美日本亚洲视频在线播放| 久久国产精品男人的天堂亚洲| 99re在线观看精品视频| 日韩精品免费视频一区二区三区| 欧美激情极品国产一区二区三区| 久久中文字幕一级| 欧美日韩视频精品一区| 久99久视频精品免费| 男人操女人黄网站| 亚洲一区二区三区欧美精品| 精品国产一区二区久久| 夫妻午夜视频| 91九色精品人成在线观看| 黑人巨大精品欧美一区二区mp4| 久久精品亚洲精品国产色婷小说| 国产av在哪里看| 少妇被粗大的猛进出69影院| 狠狠狠狠99中文字幕| 免费搜索国产男女视频| 又大又爽又粗| 这个男人来自地球电影免费观看| 俄罗斯特黄特色一大片| 一区在线观看完整版| 日本三级黄在线观看| 搡老熟女国产l中国老女人| 欧美日韩黄片免| 亚洲精品成人av观看孕妇| 曰老女人黄片| 亚洲精品中文字幕在线视频| 长腿黑丝高跟| 欧美精品一区二区免费开放| 多毛熟女@视频| 国产精品久久久久成人av| 曰老女人黄片| 黑丝袜美女国产一区| 日本五十路高清| 国产免费男女视频| 久久久久国产一级毛片高清牌| 国产成人精品久久二区二区91| 丝袜美腿诱惑在线| 久久欧美精品欧美久久欧美| 欧美乱色亚洲激情| 美女扒开内裤让男人捅视频| 日韩大码丰满熟妇| 99国产精品免费福利视频| 天天躁夜夜躁狠狠躁躁| 国产熟女xx| 成熟少妇高潮喷水视频| 久久久久国内视频| 国产欧美日韩精品亚洲av| 久久欧美精品欧美久久欧美| 热99国产精品久久久久久7| 熟女少妇亚洲综合色aaa.| 国产黄色免费在线视频| 长腿黑丝高跟| 最好的美女福利视频网| 亚洲色图av天堂| 精品国产乱子伦一区二区三区| 日韩有码中文字幕| 精品免费久久久久久久清纯| 少妇被粗大的猛进出69影院| 老司机靠b影院| 香蕉久久夜色| 久久人人精品亚洲av| 久久精品人人爽人人爽视色| 欧美精品一区二区免费开放| 日本一区二区免费在线视频| aaaaa片日本免费| 婷婷丁香在线五月| 黑人欧美特级aaaaaa片| 欧美人与性动交α欧美软件| xxxhd国产人妻xxx| 成熟少妇高潮喷水视频| 亚洲一卡2卡3卡4卡5卡精品中文| 最近最新中文字幕大全免费视频| 亚洲精品美女久久久久99蜜臀| 亚洲 国产 在线| 欧美不卡视频在线免费观看 | 51午夜福利影视在线观看| 18禁国产床啪视频网站| 免费看a级黄色片| 亚洲成av片中文字幕在线观看| 亚洲av熟女| 日韩高清综合在线| bbb黄色大片| 色老头精品视频在线观看| 日韩中文字幕欧美一区二区| 国产av一区二区精品久久| 村上凉子中文字幕在线| 天堂俺去俺来也www色官网| 一本综合久久免费| 亚洲情色 制服丝袜| 在线av久久热| 欧美色视频一区免费| 啪啪无遮挡十八禁网站| 五月开心婷婷网| 午夜福利一区二区在线看| 国产深夜福利视频在线观看| 国产亚洲av高清不卡| 一本综合久久免费| 久久青草综合色| 国产激情久久老熟女| 18禁黄网站禁片午夜丰满| 一级,二级,三级黄色视频| 成在线人永久免费视频| 中文字幕人妻丝袜一区二区| 天天影视国产精品| 999久久久国产精品视频| 少妇 在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 精品国产超薄肉色丝袜足j| 欧美国产精品va在线观看不卡| 中文字幕av电影在线播放| 在线播放国产精品三级| 日本撒尿小便嘘嘘汇集6| 国产精品偷伦视频观看了| 51午夜福利影视在线观看| 女生性感内裤真人,穿戴方法视频| 精品国产美女av久久久久小说| 久久久久亚洲av毛片大全| 动漫黄色视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 免费av毛片视频| 亚洲第一青青草原| 乱人伦中国视频| 黄色 视频免费看| 精品少妇一区二区三区视频日本电影| 搡老熟女国产l中国老女人| 精品国内亚洲2022精品成人| av网站在线播放免费| 亚洲色图 男人天堂 中文字幕| 亚洲av美国av| 亚洲一区二区三区色噜噜 | 亚洲熟女毛片儿| 激情在线观看视频在线高清| 国产亚洲精品一区二区www| 美女福利国产在线| 亚洲在线自拍视频| 免费av中文字幕在线| 成人三级做爰电影| 久久久久久亚洲精品国产蜜桃av| 国产欧美日韩综合在线一区二区| 中文欧美无线码| 18禁裸乳无遮挡免费网站照片 | 久久天躁狠狠躁夜夜2o2o| 狠狠狠狠99中文字幕| 香蕉丝袜av| 午夜福利免费观看在线| 国产黄a三级三级三级人| 宅男免费午夜| 另类亚洲欧美激情| 一区二区日韩欧美中文字幕| 国产真人三级小视频在线观看| 村上凉子中文字幕在线| 新久久久久国产一级毛片| 999久久久国产精品视频| 少妇粗大呻吟视频| 麻豆久久精品国产亚洲av | 国产亚洲精品久久久久5区| 99国产综合亚洲精品| 亚洲人成77777在线视频| 琪琪午夜伦伦电影理论片6080| 亚洲av成人av| 十八禁网站免费在线| 中文字幕人妻丝袜一区二区| 欧美黑人欧美精品刺激| 亚洲成av片中文字幕在线观看| 国产成年人精品一区二区 | 国产精品亚洲av一区麻豆| av片东京热男人的天堂| 极品人妻少妇av视频| 国产麻豆69| 亚洲欧美日韩另类电影网站| 亚洲熟妇熟女久久| 国产精品亚洲av一区麻豆| 精品少妇一区二区三区视频日本电影| 国产精品日韩av在线免费观看 | 日韩欧美一区视频在线观看| 日韩精品免费视频一区二区三区| 国产不卡一卡二| 亚洲片人在线观看| 精品国产美女av久久久久小说| 亚洲激情在线av| 美女高潮喷水抽搐中文字幕| 亚洲av成人av| 女生性感内裤真人,穿戴方法视频| 国产精品日韩av在线免费观看 | 在线观看日韩欧美| 精品一区二区三区av网在线观看| 在线观看66精品国产| 国产一区在线观看成人免费| 欧美激情高清一区二区三区| 国产成人一区二区三区免费视频网站| 久久久国产一区二区| 天堂影院成人在线观看| 免费观看精品视频网站| 国产成人精品久久二区二区免费| 欧美乱码精品一区二区三区| 一边摸一边抽搐一进一出视频| 在线观看免费视频日本深夜| 亚洲精品国产区一区二| 久久欧美精品欧美久久欧美| 天天躁狠狠躁夜夜躁狠狠躁| 欧美黑人精品巨大| 亚洲av美国av| 90打野战视频偷拍视频| 国产亚洲精品久久久久5区| 男人的好看免费观看在线视频 | 精品第一国产精品| 女同久久另类99精品国产91| 久久狼人影院| 久久99一区二区三区| 中出人妻视频一区二区| 午夜激情av网站| 91成年电影在线观看| 久久久久久久精品吃奶| 99在线视频只有这里精品首页| 国产区一区二久久| 日本三级黄在线观看| 国产精品野战在线观看 | 黑人猛操日本美女一级片| 人人澡人人妻人| 美国免费a级毛片| 麻豆成人av在线观看| 久久人人97超碰香蕉20202| 亚洲国产精品sss在线观看 | 亚洲人成77777在线视频| 精品福利永久在线观看| 欧美日韩视频精品一区| 久久久久久亚洲精品国产蜜桃av| 国产片内射在线| 亚洲中文日韩欧美视频| 老司机亚洲免费影院| 久久精品91蜜桃| 91成人精品电影| 免费在线观看完整版高清| 亚洲狠狠婷婷综合久久图片| 一个人观看的视频www高清免费观看 | 色尼玛亚洲综合影院| 欧美最黄视频在线播放免费 | 欧美日韩黄片免| 亚洲国产精品sss在线观看 | 中文字幕另类日韩欧美亚洲嫩草| 18美女黄网站色大片免费观看| 久久久国产成人免费| 婷婷丁香在线五月| 亚洲精品久久午夜乱码| 免费在线观看视频国产中文字幕亚洲| 国产精品亚洲一级av第二区| 久久欧美精品欧美久久欧美| 国产精品久久久人人做人人爽| 国产亚洲欧美在线一区二区| 一级片免费观看大全| 午夜老司机福利片| 男女高潮啪啪啪动态图| 亚洲av成人一区二区三| 黄网站色视频无遮挡免费观看| 亚洲成人免费av在线播放| 午夜免费观看网址| 少妇裸体淫交视频免费看高清 | 国产无遮挡羞羞视频在线观看| 村上凉子中文字幕在线| 99精品欧美一区二区三区四区| 亚洲精品中文字幕一二三四区| 69精品国产乱码久久久| 成人三级做爰电影| 看片在线看免费视频| 在线观看免费视频日本深夜| 精品久久久久久,| 久久久水蜜桃国产精品网| 色精品久久人妻99蜜桃| 欧美丝袜亚洲另类 | 少妇裸体淫交视频免费看高清 | 国产成人欧美| 在线永久观看黄色视频| 免费观看人在逋| 在线观看免费午夜福利视频| 天堂影院成人在线观看| 国产精品乱码一区二三区的特点 | 搡老熟女国产l中国老女人| 身体一侧抽搐| 少妇被粗大的猛进出69影院| 国产伦人伦偷精品视频| bbb黄色大片| 亚洲欧美日韩高清在线视频| 亚洲欧美精品综合一区二区三区| 国产亚洲精品第一综合不卡| 天堂√8在线中文| 欧美在线一区亚洲| 成人国语在线视频| 国产1区2区3区精品| 亚洲熟女毛片儿| 日本 av在线| 18禁美女被吸乳视频| 日本撒尿小便嘘嘘汇集6| av国产精品久久久久影院| 亚洲国产欧美网| 999精品在线视频| a级片在线免费高清观看视频| 亚洲七黄色美女视频| 亚洲人成伊人成综合网2020| av国产精品久久久久影院| 成人特级黄色片久久久久久久| 日韩欧美一区视频在线观看| 免费不卡黄色视频| 午夜视频精品福利| 欧美精品一区二区免费开放| 久热爱精品视频在线9| 少妇 在线观看| 亚洲黑人精品在线| 两性午夜刺激爽爽歪歪视频在线观看 | 国产97色在线日韩免费| 久久久久久久久中文| 欧美日韩亚洲高清精品| www.999成人在线观看| 大香蕉久久成人网| 级片在线观看| 亚洲全国av大片| 国产欧美日韩一区二区精品| 久久香蕉精品热| 国产精品av久久久久免费| av天堂久久9| 精品人妻1区二区| 又大又爽又粗| 国产成人免费无遮挡视频| 天天影视国产精品| 在线观看日韩欧美| 国产熟女xx| 在线观看一区二区三区激情| 满18在线观看网站| 老汉色∧v一级毛片| 岛国视频午夜一区免费看| 国产精品成人在线| 50天的宝宝边吃奶边哭怎么回事| 色哟哟哟哟哟哟| 国产一区二区在线av高清观看| 久久久国产成人免费| 一级毛片精品| 99国产极品粉嫩在线观看| 亚洲九九香蕉| 精品午夜福利视频在线观看一区| 757午夜福利合集在线观看| 嫩草影院精品99| 国产激情欧美一区二区| 一区二区日韩欧美中文字幕| 国产精品秋霞免费鲁丝片| 亚洲五月色婷婷综合| 亚洲人成77777在线视频| 国产人伦9x9x在线观看| 亚洲成人免费av在线播放| 亚洲色图 男人天堂 中文字幕| 又黄又爽又免费观看的视频| 亚洲人成77777在线视频| 日日爽夜夜爽网站| 婷婷精品国产亚洲av在线| 久久99一区二区三区| 成人影院久久| 久久久久久人人人人人| 极品教师在线免费播放| 亚洲精品国产一区二区精华液| 亚洲五月色婷婷综合| 日本黄色日本黄色录像| 老司机午夜十八禁免费视频| 亚洲自偷自拍图片 自拍| 80岁老熟妇乱子伦牲交| 亚洲色图 男人天堂 中文字幕| 亚洲 欧美一区二区三区| 乱人伦中国视频| 亚洲av熟女| 黄片播放在线免费| 乱人伦中国视频| 久久久久久久精品吃奶| 成人18禁高潮啪啪吃奶动态图| 少妇粗大呻吟视频| 午夜亚洲福利在线播放| 国产高清国产精品国产三级| 欧美一区二区精品小视频在线| 亚洲九九香蕉| 国产高清videossex| 日本免费一区二区三区高清不卡 | 日韩精品免费视频一区二区三区| 亚洲人成电影观看| 色综合婷婷激情| 一区福利在线观看| 波多野结衣一区麻豆| 国产伦一二天堂av在线观看| 精品国产美女av久久久久小说| 亚洲人成电影免费在线| 热99国产精品久久久久久7| 国产一卡二卡三卡精品| 国产成人一区二区三区免费视频网站| 在线十欧美十亚洲十日本专区| 日韩欧美三级三区| 操美女的视频在线观看| 久久国产精品影院| 久久久久久大精品| 欧美中文综合在线视频| 国产精品综合久久久久久久免费 | 一级黄色大片毛片| 女性生殖器流出的白浆| 免费久久久久久久精品成人欧美视频| 麻豆久久精品国产亚洲av | 夜夜躁狠狠躁天天躁| 亚洲欧美日韩无卡精品| 国产色视频综合| 国产免费男女视频| 国产精品一区二区三区四区久久 | av网站在线播放免费| 99久久久亚洲精品蜜臀av| 亚洲性夜色夜夜综合| 久久人妻av系列| 亚洲国产精品999在线| 男人操女人黄网站| 女性被躁到高潮视频| 国产午夜精品久久久久久| 欧美中文综合在线视频| 嫩草影视91久久| 在线观看午夜福利视频| 女人爽到高潮嗷嗷叫在线视频| 在线免费观看的www视频| 最新在线观看一区二区三区| 妹子高潮喷水视频| 男女下面插进去视频免费观看| 啦啦啦在线免费观看视频4| 在线观看一区二区三区| 亚洲人成伊人成综合网2020| 亚洲熟女毛片儿| 亚洲欧美一区二区三区黑人| 99精品在免费线老司机午夜| 两个人免费观看高清视频| 国产国语露脸激情在线看| 黑人巨大精品欧美一区二区mp4| 俄罗斯特黄特色一大片| 免费在线观看日本一区| 国产精品国产av在线观看| 国产成人欧美在线观看| 日韩免费av在线播放| 国产成人影院久久av| 女性被躁到高潮视频| 露出奶头的视频| 淫妇啪啪啪对白视频| 两个人免费观看高清视频| 性欧美人与动物交配| 亚洲全国av大片| 香蕉国产在线看| 久热爱精品视频在线9| 黄频高清免费视频| 在线国产一区二区在线| 一区二区三区激情视频| 久久99一区二区三区| 黑人操中国人逼视频| 欧美日本中文国产一区发布| 欧美大码av| 国产主播在线观看一区二区| 最新在线观看一区二区三区| 18禁观看日本| 欧美日韩视频精品一区| 亚洲人成网站在线播放欧美日韩| 国产精品久久久久久人妻精品电影| 咕卡用的链子| 国产精品爽爽va在线观看网站 | 一级毛片精品| 中文字幕人妻丝袜一区二区| 亚洲欧洲精品一区二区精品久久久| 99re在线观看精品视频| www国产在线视频色| 新久久久久国产一级毛片|