• <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è)計
    女性生殖器流出的白浆| 26uuu在线亚洲综合色| 亚洲精品乱码久久久v下载方式| 丝袜脚勾引网站| 国产在视频线精品| 免费大片18禁| 黑人欧美特级aaaaaa片| 视频中文字幕在线观看| 97精品久久久久久久久久精品| 国产爽快片一区二区三区| 久久久欧美国产精品| 免费观看av网站的网址| 97超视频在线观看视频| 成人无遮挡网站| 亚洲av综合色区一区| 精品视频人人做人人爽| 欧美日韩一区二区视频在线观看视频在线| 日韩制服骚丝袜av| 黄色毛片三级朝国网站| 香蕉精品网在线| 成人综合一区亚洲| 在线天堂最新版资源| freevideosex欧美| 激情五月婷婷亚洲| 亚洲欧美中文字幕日韩二区| 天天影视国产精品| 满18在线观看网站| 美女国产视频在线观看| 人妻人人澡人人爽人人| 老女人水多毛片| 亚洲性久久影院| av.在线天堂| 国产av精品麻豆| 欧美成人精品欧美一级黄| 成人二区视频| 亚洲美女搞黄在线观看| 精品人妻一区二区三区麻豆| 狠狠精品人妻久久久久久综合| 亚洲av成人精品一区久久| 伦精品一区二区三区| 国产在线免费精品| 久久久a久久爽久久v久久| 国产精品人妻久久久影院| 丝袜脚勾引网站| 国产精品久久久久久精品古装| 免费黄色在线免费观看| 最近手机中文字幕大全| 丰满少妇做爰视频| 一本大道久久a久久精品| 春色校园在线视频观看| 国产精品一二三区在线看| kizo精华| 亚洲中文av在线| 这个男人来自地球电影免费观看 | 国产亚洲欧美精品永久| 亚洲婷婷狠狠爱综合网| 美女视频免费永久观看网站| 日本黄色日本黄色录像| 视频区图区小说| av在线播放精品| 在线观看www视频免费| av专区在线播放| 你懂的网址亚洲精品在线观看| 成人国产麻豆网| 视频中文字幕在线观看| 免费黄网站久久成人精品| 日韩,欧美,国产一区二区三区| 汤姆久久久久久久影院中文字幕| 欧美国产精品一级二级三级| 久久精品国产亚洲网站| 热99久久久久精品小说推荐| 黄片无遮挡物在线观看| 免费观看性生交大片5| 日本黄色片子视频| www.av在线官网国产| 成人黄色视频免费在线看| 成人漫画全彩无遮挡| 国产乱来视频区| 欧美变态另类bdsm刘玥| 日韩av在线免费看完整版不卡| 有码 亚洲区| 精品人妻熟女av久视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲综合精品二区| av卡一久久| 欧美精品国产亚洲| 国产成人a∨麻豆精品| 最近2019中文字幕mv第一页| 国产精品久久久久久久久免| 内地一区二区视频在线| 考比视频在线观看| 九草在线视频观看| 成年美女黄网站色视频大全免费 | 国产成人精品婷婷| 一级片'在线观看视频| 2018国产大陆天天弄谢| 久久人人爽av亚洲精品天堂| 丰满少妇做爰视频| 老司机亚洲免费影院| 亚洲精品第二区| 国产日韩欧美视频二区| 欧美丝袜亚洲另类| 国产无遮挡羞羞视频在线观看| 在线亚洲精品国产二区图片欧美 | 又大又黄又爽视频免费| 国产精品欧美亚洲77777| 国产精品一区www在线观看| 精品一品国产午夜福利视频| 一区二区三区精品91| 在线观看国产h片| 99热6这里只有精品| 赤兔流量卡办理| 蜜臀久久99精品久久宅男| 制服人妻中文乱码| 国产精品久久久久久精品古装| 日本黄色日本黄色录像| 日本wwww免费看| 91久久精品国产一区二区成人| 国产欧美日韩一区二区三区在线 | 99久国产av精品国产电影| 久久久久久久亚洲中文字幕| 国产极品天堂在线| 午夜福利视频精品| 久久狼人影院| 天天躁夜夜躁狠狠久久av| 男女高潮啪啪啪动态图| 亚洲av不卡在线观看| www.av在线官网国产| 日韩精品有码人妻一区| 亚洲四区av| 久久狼人影院| 最新的欧美精品一区二区| 亚洲av.av天堂| 久久精品熟女亚洲av麻豆精品| 狂野欧美激情性bbbbbb| 一级黄片播放器| 久久久精品区二区三区| 亚洲国产精品国产精品| 精品人妻在线不人妻| 久久久久久久久久久免费av| 国产永久视频网站| 3wmmmm亚洲av在线观看| 我的女老师完整版在线观看| 国产亚洲最大av| 少妇猛男粗大的猛烈进出视频| 国产成人aa在线观看| 精品久久久噜噜| 亚洲av欧美aⅴ国产| 成人影院久久| 亚洲美女搞黄在线观看| 午夜福利网站1000一区二区三区| av一本久久久久| 国产亚洲一区二区精品| 高清视频免费观看一区二区| 久久影院123| 国产成人a∨麻豆精品| 国产精品国产三级国产专区5o| 18禁动态无遮挡网站| av黄色大香蕉| 久久毛片免费看一区二区三区| 成人亚洲精品一区在线观看| 精品国产露脸久久av麻豆| 日韩av不卡免费在线播放| 亚洲美女视频黄频| 男女边吃奶边做爰视频| 久久韩国三级中文字幕| 国产黄色视频一区二区在线观看| 伊人久久国产一区二区| 久久国产精品男人的天堂亚洲 | 成人综合一区亚洲| 五月开心婷婷网| 亚洲精品aⅴ在线观看| 精品久久久久久久久亚洲| 久久综合国产亚洲精品| 亚洲美女搞黄在线观看| 日韩人妻高清精品专区| 91在线精品国自产拍蜜月| 亚洲精品一区蜜桃| 午夜福利网站1000一区二区三区| 国产成人免费无遮挡视频| 黑人巨大精品欧美一区二区蜜桃 | 少妇人妻 视频| 亚洲国产最新在线播放| 日韩在线高清观看一区二区三区| 亚洲av.av天堂| 成年人午夜在线观看视频| 人妻 亚洲 视频| av播播在线观看一区| 日韩成人av中文字幕在线观看| 自线自在国产av| 看十八女毛片水多多多| freevideosex欧美| 蜜臀久久99精品久久宅男| 亚洲伊人久久精品综合| 黑人欧美特级aaaaaa片| 少妇人妻精品综合一区二区| 国产在线一区二区三区精| 国产在线视频一区二区| 99精国产麻豆久久婷婷| 国产国语露脸激情在线看| 色婷婷av一区二区三区视频| 丰满迷人的少妇在线观看| 水蜜桃什么品种好| 久久久国产欧美日韩av| 久久av网站| 亚洲精品视频女| 爱豆传媒免费全集在线观看| 国精品久久久久久国模美| 午夜激情福利司机影院| 一级爰片在线观看| 成人亚洲精品一区在线观看| 久久久午夜欧美精品| 精品久久蜜臀av无| 91在线精品国自产拍蜜月| 人妻少妇偷人精品九色| 亚洲精品国产av蜜桃| 国产成人午夜福利电影在线观看| 亚洲欧美精品自产自拍| 国产极品粉嫩免费观看在线 | 91精品国产国语对白视频| 成年人午夜在线观看视频| 亚洲图色成人| 人妻一区二区av| 22中文网久久字幕| 成人国语在线视频| 日日爽夜夜爽网站| 午夜激情福利司机影院| 精品久久久噜噜| 欧美日韩综合久久久久久| 看非洲黑人一级黄片| 一级二级三级毛片免费看| 久久久久国产精品人妻一区二区| 又黄又爽又刺激的免费视频.| a级毛片在线看网站| 久久99热这里只频精品6学生| 久久久久久久亚洲中文字幕| 亚洲美女黄色视频免费看| 丰满少妇做爰视频| 久久久久国产网址| 成人影院久久| 成人综合一区亚洲| 精品久久久久久电影网| 妹子高潮喷水视频| 高清黄色对白视频在线免费看| 黄色怎么调成土黄色| 久久久国产欧美日韩av| 国产国语露脸激情在线看| 我要看黄色一级片免费的| 久久久国产欧美日韩av| 国产精品一区二区三区四区免费观看| 亚洲精品亚洲一区二区| 777米奇影视久久| av不卡在线播放| 飞空精品影院首页| 考比视频在线观看| 亚洲av中文av极速乱| 精品久久蜜臀av无| 在现免费观看毛片| 日日啪夜夜爽| 人妻 亚洲 视频| 免费观看av网站的网址| 亚洲综合色惰| 能在线免费看毛片的网站| 久热久热在线精品观看| 女人久久www免费人成看片| 观看av在线不卡| 18禁在线播放成人免费| 建设人人有责人人尽责人人享有的| 欧美精品高潮呻吟av久久| 国产在线免费精品| 欧美精品一区二区大全| 女人精品久久久久毛片| 精品少妇久久久久久888优播| 亚洲av综合色区一区| 亚洲国产日韩一区二区| 3wmmmm亚洲av在线观看| 99九九线精品视频在线观看视频| av国产精品久久久久影院| 国产亚洲精品久久久com| 在线亚洲精品国产二区图片欧美 | 亚洲欧美成人综合另类久久久| 女人精品久久久久毛片| 97在线视频观看| 日韩免费高清中文字幕av| 欧美日韩亚洲高清精品| kizo精华| 大香蕉久久成人网| 丰满迷人的少妇在线观看| 久久99热6这里只有精品| 国产不卡av网站在线观看| 亚洲精品自拍成人| 少妇的逼好多水| 午夜免费鲁丝| 三级国产精品片| 亚洲欧美成人精品一区二区| 91国产中文字幕| 免费大片18禁| 黑人巨大精品欧美一区二区蜜桃 | 涩涩av久久男人的天堂| 久久婷婷青草| 女人久久www免费人成看片| 永久网站在线| 亚洲av电影在线观看一区二区三区| 成人二区视频| 九九在线视频观看精品| 亚洲综合色网址| 亚洲精品成人av观看孕妇| 99国产精品免费福利视频| 超色免费av| 一级片'在线观看视频| 久久精品久久精品一区二区三区| 亚洲av免费高清在线观看| 热re99久久精品国产66热6| 极品少妇高潮喷水抽搐| 久久av网站| 国产欧美日韩综合在线一区二区| 午夜老司机福利剧场| 免费黄频网站在线观看国产| 国产精品国产三级国产专区5o| 日日摸夜夜添夜夜爱| 黄色怎么调成土黄色| 男人操女人黄网站| 欧美三级亚洲精品| 在线观看一区二区三区激情| 99久久人妻综合| 国产精品蜜桃在线观看| 一边摸一边做爽爽视频免费| 日韩亚洲欧美综合| 精品卡一卡二卡四卡免费| 日本av手机在线免费观看| 黑人猛操日本美女一级片| 国产男人的电影天堂91| 99久久中文字幕三级久久日本| 国产欧美日韩综合在线一区二区| 中文欧美无线码| 欧美亚洲 丝袜 人妻 在线| 伊人亚洲综合成人网| 国产成人精品久久久久久| 九九在线视频观看精品| 亚洲成色77777| 国产在线免费精品| 看免费成人av毛片| 99re6热这里在线精品视频| 欧美老熟妇乱子伦牲交| 另类亚洲欧美激情| 精品视频人人做人人爽| 午夜影院在线不卡| 国产色爽女视频免费观看| 国产在线一区二区三区精| 亚洲精品国产av蜜桃| 中文字幕制服av| 韩国高清视频一区二区三区| 欧美xxxx性猛交bbbb| 亚洲精品日本国产第一区| 国产在视频线精品| av卡一久久| 毛片一级片免费看久久久久| 一区二区三区四区激情视频| av国产精品久久久久影院| av国产精品久久久久影院| 亚洲色图综合在线观看| 啦啦啦中文免费视频观看日本| 人人妻人人添人人爽欧美一区卜| 人人澡人人妻人| 日韩成人av中文字幕在线观看| 天堂8中文在线网| 免费观看性生交大片5| 亚洲,一卡二卡三卡| 久久国产亚洲av麻豆专区| 美女福利国产在线| 能在线免费看毛片的网站| 国产亚洲精品久久久com| 国产精品熟女久久久久浪| 欧美日韩av久久| 最近的中文字幕免费完整| 午夜激情av网站| a级毛色黄片| 欧美日韩亚洲高清精品| 久久精品国产亚洲av天美| 精品亚洲成a人片在线观看| 国产亚洲精品第一综合不卡 | 久久这里有精品视频免费| 国产精品人妻久久久久久| 高清视频免费观看一区二区| 黄色一级大片看看| 99精国产麻豆久久婷婷| 天天操日日干夜夜撸| av在线观看视频网站免费| 精品久久久久久久久亚洲| 色婷婷av一区二区三区视频| 自拍欧美九色日韩亚洲蝌蚪91| 一级毛片电影观看| xxxhd国产人妻xxx| 亚洲成人手机| 老司机影院成人| 亚洲精品成人av观看孕妇| 飞空精品影院首页| 极品少妇高潮喷水抽搐| 亚洲欧美色中文字幕在线| 80岁老熟妇乱子伦牲交| 午夜免费男女啪啪视频观看| 国产成人freesex在线| 99国产综合亚洲精品| 寂寞人妻少妇视频99o| 国产精品 国内视频| 亚洲美女黄色视频免费看| av福利片在线| 18禁裸乳无遮挡动漫免费视频| 97在线视频观看| 啦啦啦视频在线资源免费观看| 日韩免费高清中文字幕av| 狂野欧美激情性bbbbbb| 晚上一个人看的免费电影| 蜜桃久久精品国产亚洲av| 在线观看一区二区三区激情| 日韩在线高清观看一区二区三区| 亚洲中文av在线| 亚洲国产精品一区二区三区在线| 亚洲精品日韩在线中文字幕| 观看av在线不卡| 五月玫瑰六月丁香| 欧美日韩视频高清一区二区三区二| av卡一久久| 日韩 亚洲 欧美在线| 国产爽快片一区二区三区| 成年女人在线观看亚洲视频| 不卡视频在线观看欧美| 国产免费福利视频在线观看| 国产精品一国产av| 最近最新中文字幕免费大全7| a级毛片免费高清观看在线播放| 91久久精品国产一区二区成人| 我要看黄色一级片免费的| 建设人人有责人人尽责人人享有的| 国产在线一区二区三区精| 国产又色又爽无遮挡免| 久久免费观看电影| 久久国内精品自在自线图片| 大又大粗又爽又黄少妇毛片口| 亚洲精品乱码久久久久久按摩| 亚洲精品国产色婷婷电影| a级毛色黄片| 欧美日本中文国产一区发布| 五月伊人婷婷丁香| 午夜福利视频精品| 午夜免费观看性视频| 永久网站在线| 七月丁香在线播放| 免费观看无遮挡的男女| 秋霞在线观看毛片| 精品一区二区三卡| 99视频精品全部免费 在线| 亚洲色图 男人天堂 中文字幕 | 国产一级毛片在线| 赤兔流量卡办理| 成人免费观看视频高清| 我的老师免费观看完整版| 日韩中字成人| 精品一区二区免费观看| 人人澡人人妻人| 日韩 亚洲 欧美在线| 涩涩av久久男人的天堂| 黄色配什么色好看| 日韩av在线免费看完整版不卡| 国产一区二区在线观看av| 观看美女的网站| 在线免费观看不下载黄p国产| 久久久国产欧美日韩av| 看十八女毛片水多多多| 久久99一区二区三区| 国产成人91sexporn| 啦啦啦中文免费视频观看日本| 麻豆乱淫一区二区| 一级毛片aaaaaa免费看小| 国产精品久久久久久精品古装| 一区二区日韩欧美中文字幕 | 国精品久久久久久国模美| 日韩一区二区三区影片| 国产色婷婷99| 一级毛片黄色毛片免费观看视频| 人人妻人人澡人人看| 在线观看国产h片| 国产成人aa在线观看| 精品久久久噜噜| 99视频精品全部免费 在线| 亚洲怡红院男人天堂| 亚洲,一卡二卡三卡| 亚洲av综合色区一区| 中国美白少妇内射xxxbb| 欧美精品人与动牲交sv欧美| 五月玫瑰六月丁香| 2021少妇久久久久久久久久久| 亚洲美女搞黄在线观看| av有码第一页| 亚洲天堂av无毛| av在线播放精品| 欧美日韩亚洲高清精品| 久久精品国产亚洲av涩爱| 只有这里有精品99| 亚洲国产精品一区三区| 一边摸一边做爽爽视频免费| 久久毛片免费看一区二区三区| 伊人亚洲综合成人网| 一级a做视频免费观看| 国产精品久久久久久精品电影小说| 天天影视国产精品| 国产成人精品久久久久久| 在线观看三级黄色| 国产老妇伦熟女老妇高清| 久久久a久久爽久久v久久| 国产精品久久久久久久电影| 亚洲av男天堂| 一本大道久久a久久精品| 久久久a久久爽久久v久久| 日本免费在线观看一区| 精品久久蜜臀av无| 久久久久精品性色| 欧美变态另类bdsm刘玥| 最新的欧美精品一区二区| 国产日韩一区二区三区精品不卡 | 中文字幕制服av| 国产免费现黄频在线看| 国产成人91sexporn| 一本大道久久a久久精品| 午夜影院在线不卡| 永久免费av网站大全| 久久久精品94久久精品| 国产日韩欧美视频二区| 国产精品99久久久久久久久| 3wmmmm亚洲av在线观看| 91精品伊人久久大香线蕉| 精品人妻一区二区三区麻豆| 天堂中文最新版在线下载| 精品午夜福利在线看| 99热6这里只有精品| 99国产精品免费福利视频| 国产欧美日韩综合在线一区二区| 亚洲婷婷狠狠爱综合网| 国产精品麻豆人妻色哟哟久久| 午夜福利网站1000一区二区三区| 欧美另类一区| 婷婷色av中文字幕| 色94色欧美一区二区| 国产精品人妻久久久影院| 亚洲欧美一区二区三区国产| 久久久久久久久久成人| 国产欧美日韩综合在线一区二区| 丰满少妇做爰视频| 免费人成在线观看视频色| 亚洲av中文av极速乱| 久久精品久久精品一区二区三区| 极品少妇高潮喷水抽搐| 纵有疾风起免费观看全集完整版| 久久午夜综合久久蜜桃| 91精品国产九色| 国产熟女欧美一区二区| 男女高潮啪啪啪动态图| 在线亚洲精品国产二区图片欧美 | 伦理电影免费视频| 中文字幕制服av| 成年美女黄网站色视频大全免费 | 午夜精品国产一区二区电影| 永久网站在线| 国产精品女同一区二区软件| 亚洲av在线观看美女高潮| 边亲边吃奶的免费视频| 色视频在线一区二区三区| av黄色大香蕉| 免费高清在线观看视频在线观看| 91精品国产国语对白视频| 国产高清三级在线| 久久久午夜欧美精品| 嫩草影院入口| 久久久欧美国产精品| 亚洲欧美成人综合另类久久久| 曰老女人黄片| 黄色配什么色好看| 久久人人爽人人爽人人片va| 国产一级毛片在线| 永久免费av网站大全| 久久97久久精品| 成年女人在线观看亚洲视频| 少妇被粗大猛烈的视频| 日韩亚洲欧美综合| 久久鲁丝午夜福利片| 精品久久国产蜜桃| 99热这里只有精品一区| 国产精品 国内视频| 亚洲成人手机| 午夜视频国产福利| 精品久久久久久电影网| 熟女人妻精品中文字幕| av有码第一页| 精品一区在线观看国产| 久久精品国产亚洲av天美| 亚洲av中文av极速乱| 人妻人人澡人人爽人人| 欧美激情国产日韩精品一区| 精品午夜福利在线看| 香蕉精品网在线| 中文字幕人妻丝袜制服| 亚洲精品久久午夜乱码| 日韩一区二区视频免费看| 九九久久精品国产亚洲av麻豆| 青春草国产在线视频| 国产成人精品久久久久久| videossex国产| 夫妻性生交免费视频一级片| 欧美日韩国产mv在线观看视频| 国产在线免费精品| 美女福利国产在线| videosex国产| 国产精品久久久久久av不卡| 亚洲无线观看免费| 黄片播放在线免费|