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

    核孔膜表面重孔分布的研究

    2022-09-21 07:50:18傅楗強
    核技術(shù) 2022年9期
    關(guān)鍵詞:分析程序孔率模擬計算

    傅楗強

    (清華海峽研究院(廈門)廈門 361015)

    利用核反應(yīng)堆或加速器產(chǎn)生的重離子輻照固體膜材,在膜材上造成輻射損傷形成潛徑跡,后經(jīng)化學(xué)蝕刻,損傷區(qū)域可形成納米至微米級的穿透性通道[1?2]。用該方法制成的核孔膜具有表面堅實光滑、孔徑均一、截留率高等特性,已應(yīng)用于醫(yī)藥精密過濾、生化分析、電池隔膜、微量探測等諸多高端領(lǐng)域[3?6]。

    核孔膜的孔隙率及重孔分布直接影響了核孔膜過濾器的濾除率和滲透性[7?8]。當(dāng)孔隙率很小時,孔分布基本為單孔,此時核孔膜對大于蝕刻孔徑的固體微粒的濾除率接近100%,但滲透性較差。增大孔密度可以增加孔隙率,使核孔膜滲透性增強,但孔發(fā)生重疊的概率增大,重孔現(xiàn)象嚴(yán)重影響了孔的平均孔徑,使得核孔膜對接近蝕刻孔徑的固體微粒的濾除率大大降低。核孔膜的許多應(yīng)用中,對濾除率和滲透性的性能指標(biāo)有著嚴(yán)苛的要求,而這兩個參數(shù)通常是負(fù)相關(guān)的。因此,考察核孔膜孔隙率與重孔分布情況顯得尤為重要。

    本文基于蒙特卡羅模擬方法開發(fā)模擬分析程序,模擬計算特定輻照注量、蝕刻條件下的孔隙率與重孔分布之間的關(guān)系,并將計算結(jié)果與理論推導(dǎo)、掃描電子顯微鏡觀測結(jié)果進(jìn)行比照。重孔率的模擬計算結(jié)果與理論公式推導(dǎo)結(jié)果符合,重孔分布情況與掃描電子顯微鏡觀測結(jié)果吻合,驗證了模擬分析程序的準(zhǔn)確性。

    本研究工作對核孔膜的定制化生產(chǎn)具有實用意義,可作為核孔膜生產(chǎn)中優(yōu)化輻照注量和蝕刻孔徑大小等工藝參數(shù)的輔助工具,以滿足特定核孔膜過濾器的最佳參數(shù)配置。

    1 模擬計算方法

    1.1 蒙特卡羅方法

    蒙特卡羅方法是核物理模擬中常用的方法[9]。蒙特卡羅方法的基本思想是:建立一個概率模型或隨機過程,使它的參數(shù)或數(shù)字特征等于問題的解;通過對該模型或過程的抽樣試驗來計算這些參數(shù)或數(shù)字特征;最后給出所求解的近似值。

    本文基于蒙特卡羅思想,模擬加速器重離子入射固體膜材,后經(jīng)化學(xué)蝕刻形成孔道的過程[10],并對核孔膜表面孔道進(jìn)行統(tǒng)計分析,從而得出孔隙率、重孔分布情況等數(shù)據(jù)。

    隨機數(shù)生成是蒙特卡羅模擬的基礎(chǔ),本文中偽隨機數(shù)的產(chǎn)生及后續(xù)的統(tǒng)計分析均通過ROOT編程實現(xiàn)。ROOT是CERN開發(fā)的針對高能物理數(shù)據(jù)處理的面向?qū)ο蟮能浖蚣埽?1],主要用于粒子物理實驗的數(shù)據(jù)分析和統(tǒng)計,具有強大的作圖功能,已被廣泛應(yīng)用在高能物理、核物理等領(lǐng)域中。

    1.2 分析程序簡介

    基于ROOT編寫了核孔膜表面重孔分布分析程序。該程序開發(fā)了可視化界面,操作簡單,可快速計算出結(jié)果。程序界面如圖1所示,在工藝參數(shù)設(shè)置欄,可對孔密度、蝕刻孔徑進(jìn)行設(shè)置。在分析繪制欄,可模擬化學(xué)蝕刻,利用ROOT的作圖功能,繪制核孔膜表面蝕刻微孔的分布;可計算孔隙率和重孔分布情況,并輸出計算結(jié)果到txt文檔中。目前的程序版本中,重離子垂直入射固體膜材,對化學(xué)蝕刻采用了最簡單的圓柱形孔道蝕刻模型。

    圖1 核孔膜表面重孔分布分析程序平臺Fig.1 Analysis platform for hole distribution on the NTEM surface

    1.3 計算方法

    1.3.1 表面蝕刻微孔的繪制

    假設(shè)入射重離子在空間的分布滿足二維的均勻分布。算法步驟如下:

    1)利用ROOT中的gRandom生成器隨機生成N0個 服 從 均 勻 分 布 的 二 維 坐 標(biāo)(xi,yi),i=1,2,3,…,N0。其中:N0=n?S,n為孔密度,S為模擬計算采用的核孔膜面積,S默認(rèn)為1 mm2的正方形。

    2)利用ROOT中的TH2F類新建像素數(shù)目為k×k的二維圖像,圖像對應(yīng)的面積為S。為了保證計算精度,每個像素尺寸設(shè)為孔徑D的1/1 000,即k=

    3)遍歷坐標(biāo)數(shù)組(xi,yi),將(xi,yi)以為圓心,D為直徑區(qū)域內(nèi)的像素值置為1。

    4)繪制黑白二值圖像,即為核孔膜表面蝕刻微孔的分布圖。圖1繪制的圖像,為了顯示效果,已將顯示區(qū)域調(diào)整為孔徑的20倍。

    1.3.2 孔隙率的計算

    若忽略重孔,孔的總面積可認(rèn)為是所有單孔簡單的線性疊加,本文定義標(biāo)稱孔隙率為η,假設(shè)孔密度為n,孔呈圓形、且孔徑為D,則:

    實際的核孔膜中,由于孔之間存在重疊現(xiàn)象,實際孔隙率ηeff<η(下文將實際孔隙率簡稱為孔隙率)。

    孔隙率的計算基于§1.3.1的k×k二維圖像的結(jié)果。對二值圖像而言,用1表示孔道,用0表示膜材,孔的總面積就是統(tǒng)計圖像內(nèi)灰度值為1的像素數(shù)量。對k×k的二維矩陣進(jìn)行求和,得到孔的像素總數(shù),再除以圖像總像素數(shù)k2,即可得到孔隙率ηeff。

    1.3.3 重孔分布的統(tǒng)計

    核孔膜表面的孔可分為單孔、雙重孔、三重孔、直至m重孔。計算重孔分布的算法步驟如下:

    1)利用ROOT中的gRandom生成器隨機生成N0個 服 從 均 勻 分 布 的 二 維 坐 標(biāo)(xi,yi),i=1,2,3,…,N0,其中:N0=n?S。

    2)對于每個二維坐標(biāo)(xi,yi),對應(yīng)一個mi用于標(biāo)記該孔為m重孔,mi初始值為1。

    3)通過循環(huán)函數(shù)遍歷i從1至N0,對每一個孔的圓心(xi,yi),計算(xi,yi)與之前的所有孔的圓心(xj,yj),j

    4)對重孔分布進(jìn)行統(tǒng)計,通過mi值判斷該孔屬于單孔、雙重孔、三重孔或是m重孔,得出重孔概率分布情況(單孔概率P1,雙重孔概率P2,三重孔概率P3,m重孔概率Pm)。

    5)對雙重孔、三重孔或是m重孔,分別計算每個重孔的有效孔道面積(即統(tǒng)計該重孔區(qū)域內(nèi)灰度值為1的像素的數(shù)量),進(jìn)而統(tǒng)計得出重孔面積分布的概率密度函數(shù)。

    2 結(jié)果與討論

    2.1 重孔率的討論

    2.1.1 重孔率的理論推導(dǎo)

    本文定義重孔率為:

    式中:Noverlap和Ntotal分別是核孔膜表面重孔的孔數(shù)和孔的總數(shù)。重孔率的理論推導(dǎo)如下:

    根據(jù)碎片入射分布的統(tǒng)計獨立性,在面積F上發(fā)現(xiàn)k個孔的概率服從泊松分布[12]:

    式中:n為孔密度;nF為分布在面積F上的孔平均數(shù)目。重孔率的表達(dá)式可由式(3)推導(dǎo)得出。假設(shè)直徑為D的孔H0位于坐標(biāo)原點,孔H1是與孔H0相鄰最近的孔,孔H1與原點距離為r。如圖2所示,孔H1出現(xiàn)在r處的概率為:

    圖2 重孔概率的計算示意圖Fig.2 Schematic diagram of calculation for overlap probability

    若孔H1與孔H0發(fā)生重疊,則孔H1必然位于以原點為圓心,以D為半徑的圓形范圍內(nèi),現(xiàn)在計算這一概率。對式(4)進(jìn)行積分,得到重孔率的表達(dá)式為:

    參考式(1),可見重孔率與標(biāo)稱孔隙率η存在簡單的函數(shù)關(guān)系[13]。

    2.1.2 重孔率的模擬計算

    利用模擬分析程序,計算不同孔徑D下,重孔率Woverlap=Pm>1隨孔密度n的變化趨勢。如圖3所示,實點為模擬計算結(jié)果,實線為根據(jù)式(5)繪制的理論曲線??梢娔M計算結(jié)果與理論曲線相吻合,證明了模擬分析程序在重孔率計算上的準(zhǔn)確性。

    根據(jù)式(5),重孔率與標(biāo)稱孔隙率η之間存在的函數(shù)關(guān)系,該函數(shù)關(guān)系不受孔密度n和孔徑D影響。

    下面探討重孔率與實際孔隙率ηeff之間的關(guān)系。對圖3的數(shù)據(jù)實點,分別計算各數(shù)據(jù)點對應(yīng)的孔隙率ηeff。以孔隙率ηeff為橫坐標(biāo),重孔率為縱坐標(biāo),繪制散點圖如圖4所示。從圖4中可見,不同孔密度n、孔徑D下,散點分布也滿足一定趨勢,即重孔率與實際孔隙率ηeff也滿足一定的函數(shù)關(guān)系。

    圖3 不同孔徑下,重孔率隨孔密度的變化曲線Fig.3 The curve of overlap rate with hole density under different hole diameters

    圖4 重孔率隨孔隙率的變化曲線Fig.4 The curve of overlap rate with porosity

    對圖4的散點進(jìn)行擬合,擬合函數(shù)關(guān)系為:

    從式(6)可知,重孔率和孔隙率這兩個產(chǎn)品參數(shù)呈正相關(guān)關(guān)系,在生產(chǎn)中不能獨立控制。當(dāng)孔隙率為10%時,此時的重孔率已接近35%;當(dāng)孔隙率為30%時,重孔率接近80%。在核孔膜的許多應(yīng)用中,都需考量重孔的影響。

    2.2 重孔分布的討論

    2.2.1 重孔分布的模擬計算

    隨著孔隙率的增加,重孔現(xiàn)象越發(fā)顯著。下面對單孔概率P1,雙重孔概率P2,三重孔概率P3,m重孔概率Pm進(jìn)行分析計算。由于重孔形態(tài)的分布比較復(fù)雜,數(shù)學(xué)上難以用較為簡單的公式來描述Pm。通過模擬分析程序計算不同孔隙率下,Pm隨孔隙率ηeff的變化曲線,如圖5所示。單孔出現(xiàn)的概率P1單調(diào)遞減,而重孔出現(xiàn)的概率不斷增加;多重孔出現(xiàn)概率Pm的峰值隨著重孔數(shù)目m的增加而右移。

    圖5 多重孔分布概率隨孔隙率的變化曲線Fig.5 Variation curve of the probability of multiple hole overlap with porosity

    下面對每個重孔的面積開展統(tǒng)計分析。

    雙重孔的形態(tài)分布最為簡單,若孔H1與孔H0相重疊,則孔H1和孔H0兩者圓心的距離r必定小于D。參考圖2,r的分布滿足如下概率密度函數(shù):

    當(dāng)r=0時,孔H1與孔H0完全重疊,此時雙重孔面積最小,為πD2/4。隨著r的增大,該種形態(tài)的雙重孔出現(xiàn)的概率也增大,且對應(yīng)的雙重孔面積也增大。當(dāng)r接近D時,孔H1與孔H0略微重疊,此時雙重孔面積最大,接近πD2/2。三重以上的重孔,形態(tài)分布較為復(fù)雜,難以用解析方法描述。利用模擬計算程序?qū)χ乜椎拿娣e分布情況進(jìn)行統(tǒng)計,分別得到雙重孔、三重孔、m重孔的孔面積分布的概率密度函數(shù)PS,2、PS,3、PS,m。

    計算結(jié)果表明,不同的孔密度n和孔徑D下,重孔重疊的形態(tài)分布是相似的,PS,m概率密度分布函數(shù)也是相似的。圖6的統(tǒng)計中,已將重孔的孔面積除以單孔的面積,對橫坐標(biāo)進(jìn)行歸一化。從圖6中可見,對于m重孔,孔面積分布從1到m。雙重孔面積分布從1到2,面積接近2時幾率最大。隨著重孔數(shù)m的增加,重孔分布的形態(tài)越發(fā)復(fù)雜,孔面積分布也較為分散。對于三重孔,孔面積分布函數(shù)PS,3的極大值對應(yīng)橫坐標(biāo)為2.85。對于四重孔、五重孔,PS,3、PS,4的極大值對應(yīng)橫坐標(biāo)分別為3.25和4.10。

    不同孔隙率下,核孔膜表面孔的面積分布概率函數(shù)PS,total,可通過§1.3.3的算法用蒙特卡羅方法模擬計算得出?;趫D5和圖6的結(jié)論,可推導(dǎo)出PS,total的代數(shù)解析表達(dá)式(8),即對m重孔的出現(xiàn)概率Pm和面積分布概率密度函數(shù)PS,m的進(jìn)行加權(quán)求和。利用該解析表達(dá)式可快速計算PS,total。

    圖6 多重孔面積的概率密度分布函數(shù)Fig.6 The probability density function(PDF)of the area of multiple hole overlap

    核孔膜的應(yīng)用中,經(jīng)常同時要求高滲透性和高濾除率。高滲透性意味著高孔隙率,而高孔隙率下不可避免地出現(xiàn)高重孔率。重孔導(dǎo)致核孔膜實際濾除粒徑大于蝕刻孔徑D,因此需要減小蝕刻孔徑D來保證濾除率。在蝕刻孔徑D較小的情況下保持高孔隙率,需要增大孔密度n。綜上所述,高孔密度且小蝕刻孔徑的核孔膜可同時滿足高滲透性和高濾除率要求。而核孔膜生產(chǎn)中,加速器束流是主要的成本,高孔密度意味著輻照時長變長,成本增加;另一方面,蝕刻過程中對數(shù)十nm以下的孔徑需要更為復(fù)雜的工藝。

    利用該模擬分析程序,在滿足核孔膜性能指標(biāo)的情況下,優(yōu)化生產(chǎn)中輻照注量和蝕刻孔徑等工藝參數(shù),盡量節(jié)約生產(chǎn)成本和簡化生產(chǎn)工藝。

    2.2.2 重孔分布的實驗結(jié)果

    為了驗證模擬程序的可靠性,對兩種核孔膜樣本采用掃描電子顯微鏡(Scanning Electron Microscope,SEM)進(jìn)行觀測[14]。利用電子顯微鏡照片,統(tǒng)計得到視野范圍內(nèi)孔密度和孔徑大小,并對重孔分布情況進(jìn)行統(tǒng)計。

    根據(jù)觀測的孔密度和孔徑大小,利用模擬分析程序繪制相同參數(shù)下的孔分布圖片,并計算其相應(yīng)的重孔分布情況。

    圖7展示了核孔膜的電子顯微鏡照片和模擬分析程序繪制的圖片,圖中核孔膜的孔密度和孔徑等參數(shù)詳見表1。對比模擬分析性程序繪制的圖片和電子顯微鏡的觀測照片,可見兩者的孔分布在宏觀上對照得很好。

    圖7 核孔膜的電子顯微鏡照片(a、b)和模擬分析程序繪制的表面孔分布(c、d)Fig.7 SEM images of nuclear track-etched membranes(a,b)and pictures of hole distribution draw by the analysis platform(c,d)

    表1 核孔膜樣本觀測的多重孔概率和模擬結(jié)果的對比Table 1 The comparison of the probability of multiple hole overlap between experimental observation and simulation

    對掃描電子顯微鏡下多個視野中的核孔膜照片開展重孔統(tǒng)計分析,得到單孔概率P1,雙重孔概率P2,m重孔概率Pm。表1對比了核孔膜樣本的重孔分布觀測結(jié)果和模擬分析程序計算結(jié)果。

    理論上核孔膜表面的微孔在二維空間是均勻分布的,但采用電子顯微鏡觀測時,由于視野范圍很小,不同視野范圍內(nèi)觀測到的微孔數(shù)目及分布均會存在一定的偏差。表1統(tǒng)計的數(shù)據(jù),電子顯微鏡觀測視野面積總計約為0.1 mm2。模擬分析程序中也采用了0.1 mm2的面積,為減小采樣的統(tǒng)計漲落,模擬分析程序重復(fù)采樣1 000次后取平均值。在測量誤差范圍內(nèi),模擬分析程序的計算結(jié)果與實驗觀測結(jié)果匹配、符合得較好,驗證了模擬分析程序的可靠性。

    3 結(jié)語

    核孔膜已被應(yīng)用于在諸多高端領(lǐng)域。核孔膜的孔隙率及重孔分布直接影響了核孔膜過濾器的濾除率和滲透性,因此對核孔膜表面孔隙率與重孔分布的研究具有實際應(yīng)用價值。本文基于蒙特卡羅方法開發(fā)了核孔膜表面重孔分布分析程序,模擬計算了不同孔隙率下單孔、雙重孔、多重孔出現(xiàn)概率的分布,以及重孔面積的分布情況。通過將模擬計算結(jié)果與理論推導(dǎo)值、核孔膜樣本掃描電子顯微鏡觀測結(jié)果進(jìn)行比照,驗證了模擬分析程序的準(zhǔn)確性。

    研究表明,重孔率和孔隙率這兩個產(chǎn)品參數(shù)滿足一定的函數(shù)關(guān)系,且呈正相關(guān)關(guān)系,在生產(chǎn)中不能獨立控制。隨孔隙率ηeff的增加,核孔膜表面單孔出現(xiàn)的概率P1單調(diào)遞減,而重孔出現(xiàn)的概率不斷增加;多重孔出現(xiàn)概率Pm的峰值隨著重孔數(shù)目m的增加而右移。不同的孔密度n和孔徑D下,重孔重疊的形態(tài)分布是相似的,PS,m概率密度分布函數(shù)也是相似的?;趍重孔的出現(xiàn)概率Pm和面積分布概率密度函數(shù)PS,m,得出重孔面積分布的代數(shù)解析表達(dá)式,可用于快速計算不同孔隙率下重孔面積的分布。為滿足高滲透性和高濾除率的性能要求,核孔膜生產(chǎn)中通常需要采用高孔密度和小蝕刻孔徑。利用該模擬分析程序,可優(yōu)化生產(chǎn)中輻照注量和蝕刻孔徑等工藝參數(shù),節(jié)約生產(chǎn)成本和簡化生產(chǎn)工藝。

    致謝感謝歐洲核子中心吳曉偉博士關(guān)于重孔率理論計算的有益探討。

    作者貢獻(xiàn)聲明傅楗強:負(fù)責(zé)研究的提出及設(shè)計、數(shù)據(jù)的收集和整理、文章的起草和最終版本的修訂。

    猜你喜歡
    分析程序孔率模擬計算
    管控經(jīng)營風(fēng)險,以分析程序提升企業(yè)財務(wù)報表審計效能
    中國商人(2022年9期)2022-09-02 08:30:42
    R1234ze PVTx熱物性模擬計算
    能源工程(2022年1期)2022-03-29 01:06:26
    風(fēng)屏障開孔率對高速列車氣動力的影響
    恐龍那么重,為什么沒有壓壞自己的蛋
    恐龍那么重,為什么沒有壓壞自己的蛋?
    大自然探索(2018年7期)2018-09-01 04:36:24
    擠出發(fā)泡片材褶皺分析及模擬計算
    中國塑料(2016年9期)2016-06-13 03:18:54
    基于小波包變換的樂音時—頻綜合分析程序的開發(fā)
    音樂探索(2015年1期)2015-04-27 01:01:52
    實際發(fā)射工況下底排藥柱結(jié)構(gòu)完整性的模擬計算
    包裝袋開孔率對荔枝低溫貨架銷售品質(zhì)的影響
    丙烯酸酯類降凝劑的Monte Carlo模擬計算及分子結(jié)構(gòu)設(shè)計
    亚洲精品乱久久久久久| 天堂影院成人在线观看| 激情五月婷婷亚洲| 午夜福利在线在线| 免费av不卡在线播放| 国产日韩欧美在线精品| 51国产日韩欧美| 精品一区在线观看国产| 91午夜精品亚洲一区二区三区| 波野结衣二区三区在线| 成人欧美大片| 最近最新中文字幕大全电影3| 免费黄色在线免费观看| 国模一区二区三区四区视频| 精品久久国产蜜桃| 看非洲黑人一级黄片| 在线观看av片永久免费下载| 日韩在线高清观看一区二区三区| 蜜桃久久精品国产亚洲av| 99热这里只有是精品在线观看| 色综合站精品国产| 成人毛片60女人毛片免费| 波野结衣二区三区在线| 久久久久免费精品人妻一区二区| 赤兔流量卡办理| 午夜精品国产一区二区电影 | 国产永久视频网站| 别揉我奶头 嗯啊视频| 亚洲精品国产av成人精品| 国产av在哪里看| 国产高潮美女av| 一个人看的www免费观看视频| 不卡视频在线观看欧美| 国产伦精品一区二区三区四那| 又爽又黄无遮挡网站| 久久午夜福利片| 久久综合国产亚洲精品| 中文精品一卡2卡3卡4更新| 免费看a级黄色片| 一区二区三区高清视频在线| 人体艺术视频欧美日本| 人人妻人人看人人澡| 中文字幕av在线有码专区| 亚洲久久久久久中文字幕| 午夜精品在线福利| 欧美成人一区二区免费高清观看| 国产精品女同一区二区软件| 色综合站精品国产| 男女下面进入的视频免费午夜| 免费看a级黄色片| 国产极品天堂在线| 国产乱人视频| 国内少妇人妻偷人精品xxx网站| 国产免费福利视频在线观看| 国产免费福利视频在线观看| 免费观看av网站的网址| 九草在线视频观看| 国产亚洲精品av在线| 久久精品熟女亚洲av麻豆精品 | 国产伦理片在线播放av一区| 日韩电影二区| 精品一区二区三卡| 国产高清三级在线| 极品少妇高潮喷水抽搐| 乱码一卡2卡4卡精品| 日本三级黄在线观看| 天堂俺去俺来也www色官网 | 久99久视频精品免费| 欧美xxxx性猛交bbbb| 亚洲四区av| 国语对白做爰xxxⅹ性视频网站| 午夜福利在线观看吧| 天天躁日日操中文字幕| 午夜福利视频1000在线观看| 亚洲精品影视一区二区三区av| 亚洲国产最新在线播放| 色播亚洲综合网| 久久99热这里只频精品6学生| 国产精品无大码| 国产午夜福利久久久久久| 夜夜看夜夜爽夜夜摸| 精品熟女少妇av免费看| 欧美三级亚洲精品| 日本免费在线观看一区| 欧美区成人在线视频| 大陆偷拍与自拍| 最近中文字幕2019免费版| 老司机影院毛片| 夫妻性生交免费视频一级片| 日本熟妇午夜| 国产伦一二天堂av在线观看| 80岁老熟妇乱子伦牲交| 18+在线观看网站| 美女被艹到高潮喷水动态| 国产又色又爽无遮挡免| 午夜福利视频1000在线观看| 国产v大片淫在线免费观看| 插阴视频在线观看视频| 国产精品一区二区三区四区久久| 久久久精品免费免费高清| 丝瓜视频免费看黄片| 国产一级毛片在线| 青春草国产在线视频| 色视频www国产| 天堂俺去俺来也www色官网 | av黄色大香蕉| 亚洲三级黄色毛片| 久久精品国产亚洲av天美| 久久精品国产亚洲网站| 午夜激情福利司机影院| 精品人妻熟女av久视频| 成年人午夜在线观看视频 | 国产亚洲一区二区精品| 免费黄色在线免费观看| 国产精品一二三区在线看| 久久草成人影院| 青春草国产在线视频| 国产v大片淫在线免费观看| 亚洲精品456在线播放app| 亚洲成人一二三区av| 99热这里只有精品一区| 内地一区二区视频在线| 亚洲久久久久久中文字幕| 直男gayav资源| 国产 一区 欧美 日韩| 精品亚洲乱码少妇综合久久| 亚洲av在线观看美女高潮| 亚洲不卡免费看| 久久久久久久国产电影| 别揉我奶头 嗯啊视频| 中文字幕人妻熟人妻熟丝袜美| 亚洲av成人精品一区久久| av在线老鸭窝| 一级爰片在线观看| 青春草视频在线免费观看| 日韩成人伦理影院| 十八禁网站网址无遮挡 | 天堂av国产一区二区熟女人妻| 青春草视频在线免费观看| 国产成人a区在线观看| 精品国内亚洲2022精品成人| 成年女人在线观看亚洲视频 | 欧美成人精品欧美一级黄| 免费黄色在线免费观看| 韩国av在线不卡| 99热全是精品| 久久精品久久久久久久性| 免费看日本二区| 午夜精品国产一区二区电影 | 美女xxoo啪啪120秒动态图| 久久精品久久精品一区二区三区| 少妇猛男粗大的猛烈进出视频 | 69人妻影院| 日日干狠狠操夜夜爽| 午夜福利视频精品| 国产亚洲精品久久久com| 丰满少妇做爰视频| 2021少妇久久久久久久久久久| av女优亚洲男人天堂| 午夜爱爱视频在线播放| 久热久热在线精品观看| 最近最新中文字幕大全电影3| 亚洲av在线观看美女高潮| 亚洲精品一区蜜桃| 国产精品麻豆人妻色哟哟久久 | 2021天堂中文幕一二区在线观| 久久99热这里只有精品18| 3wmmmm亚洲av在线观看| 精华霜和精华液先用哪个| 热99在线观看视频| 夫妻午夜视频| 91在线精品国自产拍蜜月| 午夜福利视频1000在线观看| 亚洲精品456在线播放app| 国产一区有黄有色的免费视频 | 三级男女做爰猛烈吃奶摸视频| 22中文网久久字幕| 亚洲熟女精品中文字幕| 爱豆传媒免费全集在线观看| 床上黄色一级片| 天美传媒精品一区二区| 一本一本综合久久| 麻豆精品久久久久久蜜桃| 久久久久久久国产电影| 亚洲精品一二三| 肉色欧美久久久久久久蜜桃 | 蜜臀久久99精品久久宅男| 国产亚洲91精品色在线| 18禁在线无遮挡免费观看视频| 男人爽女人下面视频在线观看| 国产精品伦人一区二区| av在线老鸭窝| 国产精品一区二区三区四区免费观看| freevideosex欧美| xxx大片免费视频| 欧美成人a在线观看| 好男人在线观看高清免费视频| 亚洲久久久久久中文字幕| 亚洲欧美成人精品一区二区| 免费看光身美女| 国产白丝娇喘喷水9色精品| 麻豆av噜噜一区二区三区| 18禁在线无遮挡免费观看视频| 日日干狠狠操夜夜爽| av在线播放精品| 亚洲四区av| 久久99热6这里只有精品| 干丝袜人妻中文字幕| 午夜福利成人在线免费观看| av在线老鸭窝| 波多野结衣巨乳人妻| 黄色一级大片看看| 国产精品一区二区在线观看99 | 国产精品一区二区性色av| 日韩人妻高清精品专区| 国产精品1区2区在线观看.| 亚洲国产高清在线一区二区三| 永久网站在线| 亚洲成人中文字幕在线播放| 亚洲av免费在线观看| 亚洲av电影在线观看一区二区三区 | 一区二区三区乱码不卡18| 成年人午夜在线观看视频 | 欧美97在线视频| 看黄色毛片网站| 亚洲人成网站在线观看播放| 天堂av国产一区二区熟女人妻| 欧美一级a爱片免费观看看| 国产精品人妻久久久久久| 国产精品99久久久久久久久| 女人久久www免费人成看片| 熟女电影av网| 久久久久久久久久人人人人人人| 夫妻性生交免费视频一级片| 男的添女的下面高潮视频| 黄片无遮挡物在线观看| eeuss影院久久| 亚洲国产欧美在线一区| 丰满乱子伦码专区| 青春草亚洲视频在线观看| 亚洲国产最新在线播放| 国产精品久久久久久精品电影小说 | 色吧在线观看| 日本wwww免费看| 97超视频在线观看视频| 久久久精品94久久精品| 亚洲国产精品成人综合色| 亚洲欧美一区二区三区国产| 久久99热6这里只有精品| 国产91av在线免费观看| 免费少妇av软件| 人人妻人人看人人澡| 亚洲成色77777| 亚洲av国产av综合av卡| 国产视频首页在线观看| 亚洲欧美一区二区三区国产| 联通29元200g的流量卡| 国产成人精品福利久久| 边亲边吃奶的免费视频| 日韩中字成人| 国产精品av视频在线免费观看| 亚洲av福利一区| 一级毛片久久久久久久久女| 一本久久精品| 超碰97精品在线观看| 十八禁国产超污无遮挡网站| 中文字幕av在线有码专区| 日日摸夜夜添夜夜添av毛片| 国产有黄有色有爽视频| 亚洲国产精品sss在线观看| 一级毛片 在线播放| 午夜亚洲福利在线播放| 亚洲国产精品成人久久小说| 青春草国产在线视频| 国产亚洲av片在线观看秒播厂 | 日韩成人伦理影院| 精品熟女少妇av免费看| 亚洲av福利一区| 免费观看的影片在线观看| 婷婷色av中文字幕| 99热这里只有是精品50| 国产精品.久久久| 婷婷色综合大香蕉| 欧美日韩视频高清一区二区三区二| 乱系列少妇在线播放| 国产欧美日韩精品一区二区| 高清午夜精品一区二区三区| 亚洲成人av在线免费| 久久99精品国语久久久| 久久久久久久久大av| 偷拍熟女少妇极品色| 麻豆成人午夜福利视频| 如何舔出高潮| 国产又色又爽无遮挡免| 亚洲人成网站在线观看播放| 三级男女做爰猛烈吃奶摸视频| 天天躁夜夜躁狠狠久久av| 搡女人真爽免费视频火全软件| 99热网站在线观看| av在线播放精品| 午夜激情福利司机影院| 亚洲在线观看片| 十八禁网站网址无遮挡 | 亚洲国产av新网站| 久久久久久久久久黄片| 天堂√8在线中文| 最后的刺客免费高清国语| 麻豆av噜噜一区二区三区| 午夜视频国产福利| 亚洲成色77777| 亚洲av中文av极速乱| 国产成人免费观看mmmm| 国产精品国产三级国产av玫瑰| 美女xxoo啪啪120秒动态图| 欧美另类一区| 久久久久九九精品影院| 六月丁香七月| 一本一本综合久久| 亚洲精品国产av蜜桃| 亚洲人成网站在线观看播放| 午夜福利在线观看免费完整高清在| 男人舔女人下体高潮全视频| 夫妻午夜视频| 午夜老司机福利剧场| 免费在线观看成人毛片| 国产日韩欧美在线精品| 黄色一级大片看看| 亚洲成人精品中文字幕电影| 九色成人免费人妻av| 欧美成人a在线观看| 少妇人妻精品综合一区二区| 天堂网av新在线| 三级经典国产精品| 亚洲精品久久久久久婷婷小说| 亚洲精品乱久久久久久| 欧美另类一区| 18禁在线无遮挡免费观看视频| 麻豆精品久久久久久蜜桃| 黄色欧美视频在线观看| 黄色日韩在线| 欧美xxⅹ黑人| 丰满人妻一区二区三区视频av| 国产成人a∨麻豆精品| 亚洲欧美一区二区三区国产| 乱码一卡2卡4卡精品| www.av在线官网国产| 蜜臀久久99精品久久宅男| 美女黄网站色视频| 在线观看美女被高潮喷水网站| 日韩亚洲欧美综合| 丝瓜视频免费看黄片| 91狼人影院| 成人毛片60女人毛片免费| 汤姆久久久久久久影院中文字幕 | 国产永久视频网站| 久久久成人免费电影| 一级毛片久久久久久久久女| videossex国产| 免费在线观看成人毛片| 欧美一级a爱片免费观看看| 高清午夜精品一区二区三区| 少妇人妻一区二区三区视频| 久久人人爽人人爽人人片va| 欧美成人一区二区免费高清观看| 色尼玛亚洲综合影院| 国产av码专区亚洲av| 2018国产大陆天天弄谢| 简卡轻食公司| 丝袜喷水一区| 亚洲色图av天堂| 日韩欧美精品v在线| 最近视频中文字幕2019在线8| 亚洲精品久久午夜乱码| 水蜜桃什么品种好| 午夜久久久久精精品| 国产成人一区二区在线| 亚洲精品视频女| 亚洲精品国产成人久久av| 99热全是精品| 日韩一区二区视频免费看| 精品不卡国产一区二区三区| 亚洲伊人久久精品综合| 天天躁日日操中文字幕| 国产乱人偷精品视频| 亚洲性久久影院| 国产高潮美女av| 国产精品一及| 丝袜美腿在线中文| 99久国产av精品| 国产精品不卡视频一区二区| 天天躁夜夜躁狠狠久久av| 精华霜和精华液先用哪个| 爱豆传媒免费全集在线观看| 亚洲最大成人手机在线| 亚洲在线观看片| 亚洲成人av在线免费| 欧美三级亚洲精品| av国产免费在线观看| 天美传媒精品一区二区| 99久国产av精品| 少妇高潮的动态图| 日本wwww免费看| 欧美成人一区二区免费高清观看| 久久久久精品性色| 欧美 日韩 精品 国产| 91狼人影院| 亚洲成人精品中文字幕电影| 天堂√8在线中文| 97热精品久久久久久| 亚洲精品乱码久久久v下载方式| 国产成人精品福利久久| 97人妻精品一区二区三区麻豆| 日日干狠狠操夜夜爽| 亚洲人成网站在线播| 久久久久国产网址| 男女下面进入的视频免费午夜| 青青草视频在线视频观看| 男人狂女人下面高潮的视频| 久久久久久久久久成人| 久久99精品国语久久久| a级毛片免费高清观看在线播放| 三级国产精品片| av在线亚洲专区| 街头女战士在线观看网站| 少妇裸体淫交视频免费看高清| 国产精品蜜桃在线观看| 天天一区二区日本电影三级| 少妇被粗大猛烈的视频| 国产黄片视频在线免费观看| av线在线观看网站| 全区人妻精品视频| 搡女人真爽免费视频火全软件| 亚洲婷婷狠狠爱综合网| 精品亚洲乱码少妇综合久久| 在线天堂最新版资源| 亚洲一级一片aⅴ在线观看| 男人爽女人下面视频在线观看| 青春草视频在线免费观看| 中文字幕av成人在线电影| 国产av不卡久久| 狂野欧美白嫩少妇大欣赏| 久久久久免费精品人妻一区二区| 日本色播在线视频| 日日啪夜夜撸| 国产爱豆传媒在线观看| 有码 亚洲区| 建设人人有责人人尽责人人享有的 | av免费在线看不卡| 国产精品久久视频播放| 如何舔出高潮| 免费av不卡在线播放| 美女黄网站色视频| 亚洲精品日本国产第一区| 欧美变态另类bdsm刘玥| 亚洲乱码一区二区免费版| 卡戴珊不雅视频在线播放| 亚洲av在线观看美女高潮| 夜夜看夜夜爽夜夜摸| 一级毛片 在线播放| 精品人妻熟女av久视频| 国产探花在线观看一区二区| 国产乱人视频| 91久久精品国产一区二区成人| 99热这里只有是精品50| 蜜臀久久99精品久久宅男| 91久久精品国产一区二区成人| 高清在线视频一区二区三区| 成人毛片a级毛片在线播放| 噜噜噜噜噜久久久久久91| 99热全是精品| 国产精品.久久久| 久久久久久久久大av| 熟女人妻精品中文字幕| 欧美变态另类bdsm刘玥| 我的老师免费观看完整版| 纵有疾风起免费观看全集完整版 | 欧美日韩国产mv在线观看视频 | 777米奇影视久久| 激情 狠狠 欧美| 狠狠精品人妻久久久久久综合| 最近中文字幕2019免费版| 久久久精品94久久精品| 亚洲国产成人一精品久久久| 日本-黄色视频高清免费观看| 一边亲一边摸免费视频| 国产精品一区www在线观看| 国产成人一区二区在线| 亚洲欧洲国产日韩| 欧美xxxx性猛交bbbb| .国产精品久久| 简卡轻食公司| 亚洲性久久影院| 九色成人免费人妻av| 国产91av在线免费观看| 日本黄大片高清| 搡老乐熟女国产| 亚洲人成网站在线播| 欧美日韩一区二区视频在线观看视频在线 | 一夜夜www| 老师上课跳d突然被开到最大视频| 国产免费一级a男人的天堂| 久久久精品欧美日韩精品| 特大巨黑吊av在线直播| 人妻夜夜爽99麻豆av| 日本av手机在线免费观看| 丝袜喷水一区| 久久久午夜欧美精品| 午夜老司机福利剧场| 亚洲国产日韩欧美精品在线观看| 国产午夜精品久久久久久一区二区三区| 插阴视频在线观看视频| 80岁老熟妇乱子伦牲交| 少妇丰满av| 五月天丁香电影| 岛国毛片在线播放| 有码 亚洲区| 免费黄网站久久成人精品| 看非洲黑人一级黄片| 免费看日本二区| 亚洲乱码一区二区免费版| 亚洲精品中文字幕在线视频 | 成年版毛片免费区| 最近视频中文字幕2019在线8| 国产伦理片在线播放av一区| 99热网站在线观看| 婷婷色av中文字幕| 国产精品无大码| 国产单亲对白刺激| 精华霜和精华液先用哪个| 亚洲成人一二三区av| 免费高清在线观看视频在线观看| 日韩 亚洲 欧美在线| 国产高潮美女av| 久久久久网色| 精品一区二区三区人妻视频| 日韩亚洲欧美综合| www.av在线官网国产| 免费大片18禁| 高清视频免费观看一区二区 | 观看免费一级毛片| av播播在线观看一区| 国产精品美女特级片免费视频播放器| 夫妻性生交免费视频一级片| 色综合站精品国产| 中文字幕免费在线视频6| 91久久精品电影网| 精品久久久精品久久久| 欧美高清成人免费视频www| 深夜a级毛片| 精品人妻熟女av久视频| 黄色日韩在线| 舔av片在线| 午夜久久久久精精品| 亚洲精品成人久久久久久| 亚洲人与动物交配视频| 欧美区成人在线视频| 91在线精品国自产拍蜜月| 国产欧美日韩精品一区二区| 七月丁香在线播放| 成人亚洲精品av一区二区| av网站免费在线观看视频 | 日本一二三区视频观看| 中文资源天堂在线| 久久久久久国产a免费观看| 日本三级黄在线观看| 精华霜和精华液先用哪个| 久久精品久久久久久久性| 国产永久视频网站| 天天躁日日操中文字幕| av在线亚洲专区| 亚洲内射少妇av| 中文字幕久久专区| 五月玫瑰六月丁香| 尤物成人国产欧美一区二区三区| 久久99蜜桃精品久久| 亚洲最大成人手机在线| 色网站视频免费| 日本爱情动作片www.在线观看| 国产午夜精品论理片| 国产成年人精品一区二区| 免费黄网站久久成人精品| 亚洲欧美精品专区久久| av女优亚洲男人天堂| 欧美区成人在线视频| 嫩草影院入口| 欧美激情久久久久久爽电影| 18禁在线播放成人免费| 免费观看精品视频网站| 一级毛片 在线播放| 特大巨黑吊av在线直播| 自拍偷自拍亚洲精品老妇| 日韩 亚洲 欧美在线| 夜夜看夜夜爽夜夜摸| 亚洲欧洲日产国产| 国产伦精品一区二区三区四那| 我的女老师完整版在线观看| 免费黄网站久久成人精品| 一级毛片 在线播放| 黄色欧美视频在线观看| 精品一区在线观看国产| 亚洲,欧美,日韩| 毛片一级片免费看久久久久| 亚洲国产高清在线一区二区三| 在线免费十八禁| 国产黄频视频在线观看| 国产69精品久久久久777片| 婷婷色综合大香蕉| 午夜久久久久精精品| 少妇的逼好多水| 尾随美女入室| 少妇猛男粗大的猛烈进出视频 | 伦精品一区二区三区| 国产免费视频播放在线视频 | 最近视频中文字幕2019在线8| 午夜精品一区二区三区免费看|