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

    基于滲透率修正因子的氣體有效擴散系數(shù)分形模型

    2021-06-06 22:44:52王昌進徐靜磊
    巖性油氣藏 2021年3期
    關(guān)鍵詞:擴散系數(shù)維數(shù)分形

    王昌進,張 賽,徐靜磊

    (昆明理工大學機電工程學院,昆明 650500)

    0 引言

    多孔介質(zhì)中的氣體擴散問題是能源開采、燃料電池、化學工程等領(lǐng)域研究的重點,但是氣體在多孔介質(zhì)中的擴散包含了很多復(fù)雜的傳質(zhì)機理,如何建立多孔介質(zhì)中氣體有效擴散系數(shù)模型成為研究熱點[1]。

    傳統(tǒng)的物質(zhì)擴散定律一般由菲克定律來描述,該定律描述了物質(zhì)擴散通量與濃度梯度之間的定量關(guān)系。Cao 等[2]通過建立多孔介質(zhì)三維模型,證明在穩(wěn)態(tài)情況下,菲克定律能夠在多孔介質(zhì)氣體擴散中適用,然而該模型使用的是100×100×100 的立方網(wǎng)絡(luò)模型,對于多孔介質(zhì)的不規(guī)則性和復(fù)雜性方面的描述不夠準確。氣體在多孔介質(zhì)中復(fù)雜的輸運機理與多孔介質(zhì)孔隙大小和結(jié)構(gòu)差異有很大關(guān)系,根據(jù)多孔介質(zhì)孔隙大小和結(jié)構(gòu)的不同,多孔介質(zhì)中的氣體流動與擴散常常劃分為4 個流動區(qū)域:連續(xù)介質(zhì)區(qū)、滑移區(qū)、過渡區(qū)和自由分子區(qū)[3]。Beskok等[4]通過大量實驗給出了一維圓管中氣體體積流量的表達式,該表達式引入克努森數(shù)kn和滲透率修正因子f(kn),可以適用于所有流動區(qū)域。符東宇等[5]運用格子Boltzmann 方法研究了頁巖氣藏的滲流規(guī)律,發(fā)現(xiàn)克努森數(shù)kn是影響頁巖氣滲流的主要參數(shù)。Roy 等[6]使用Ar,N2,O2和Al2O3作為實驗材料進行了擴散實驗,引入了克努森擴散系數(shù),建立了微納米孔道中的氣體擴散方程,但是該方程只能適用于自由分子區(qū)。張烈輝等[7]考慮了多重運移機制,建立了頁巖氣藏的綜合滲流模型。李亞雄等[8]通過研究頁巖氣的滑脫和擴散機理,提出了“壁聯(lián)擴散”的新概念來表征克努森擴散和表面擴散的總效應(yīng),但該方法的數(shù)學模型與論證過程還有待完善。在一些科學研究和實際情況中,例如川南地區(qū)頁巖孔隙結(jié)構(gòu)和川西坳陷致密砂巖孔隙結(jié)構(gòu),它們的孔道直徑很小,達到了微納米級別,此時,孔隙結(jié)構(gòu)顯示出分形的特性。陳居凱等[9],朱漢卿等[10]和鄧浩陽等[11]均使用分形理論研究了頁巖與致密砂巖等致密多孔介質(zhì)的孔隙結(jié)構(gòu)。目前,也有很多文獻將分形理論運用到氣體滲透和擴散中[12-14]。Zheng等[15]基于分形理論和毛細管束模型,得到了氣體的有效擴散系數(shù)模型。該模型能與已有實驗數(shù)據(jù)相吻合,并且能適用于不同的多孔介質(zhì)。Shi 等[16]在考慮了2 種運輸機制的情況下,運用分形理論建立了H2在燃料電池氣體擴散層中的有效擴散系數(shù)模型,并且分析了分形理論相關(guān)參數(shù)(孔隙度、分形維數(shù)等)的影響。Zhang[17]考慮了3 種多孔介質(zhì)中的氣體流動機理,得到了滲透率和氣體有效擴散分形模型,但忽略了氣體的滑脫效應(yīng)。Woignier 等[18]研究了分形多孔介質(zhì)中氣體在滑移區(qū)的擴散機理并建立了滑移模型。張賽等[19]和牟新竹等[20]在研究多尺度分形多孔介質(zhì)氣體擴散時加入了孔隙連通性參數(shù),得到了新的數(shù)學模型。以上研究大多只考慮了1 種或2 種氣體擴散機理,不能完整地描述氣體在多孔介質(zhì)中的擴散問題。

    基于Beskok 等[4]給出的表達式、迂曲毛細管束模型、分形理論和菲克定律,建立氣體在多孔介質(zhì)中擴散的分形模型;該模型包含面積分形維數(shù)、迂曲分形維數(shù)、孔隙度、最大最小孔隙直徑、氣體密度、黏度、摩爾質(zhì)量等參數(shù),然后分析這些參數(shù)與氣體有效擴散系數(shù)之間的函數(shù)關(guān)系與變化趨勢,以期能夠正確地描述氣體在多孔介質(zhì)中的擴散規(guī)律,并應(yīng)用于油氣田的開發(fā)。

    1 多孔介質(zhì)的分形特征

    1.1 孔隙數(shù)目與橫截面積描述

    在具有分形特性的多孔介質(zhì)中,孔隙直徑大于等于λ的孔隙數(shù)量的總和與孔隙直徑大小滿足下面的分形標度關(guān)系[21]

    式中:N為孔隙長度大于孔隙直徑的多孔介質(zhì)的孔隙總數(shù)目;L為毛細管長度,m;λmax與λ分別為最大孔隙直徑和孔隙直徑,m;Df為無因式多孔介質(zhì)的面積分形維數(shù)。

    在真實多孔介質(zhì)中,孔隙數(shù)目十分巨大,根據(jù)統(tǒng)計理論和方法,可以將式(1)看成是連續(xù)和可微的函數(shù),對式(1)進行微分可以得到

    式(2)得到了在(λ,λ+dλ) 的孔隙數(shù)目,表明了孔隙數(shù)目隨著孔隙直徑的增大而減少,-dN>0。由于截面上的孔隙可以看作是直徑不同的圓,因此通過式(2),用微元法可得到分形集中截面上總的孔隙面積AP為

    式中:AP為截面上的總的孔隙面積,m2;λmin為最小孔隙直徑,m。

    則總的截面積At可以計算得到

    式中:At為總的截面積,m2;ε為面孔隙度。

    一般情況下,體孔隙度和面孔隙度不相等,兩者之間相差一個迂曲度因子

    式中:Φ為體孔隙度;為迂曲度因子[1]。在計算多孔介質(zhì)相關(guān)參數(shù)時通常使用的是體孔隙度,通過式(5)可以實現(xiàn)面孔隙度與體孔隙度的轉(zhuǎn)化。

    上述式(1)—(5)描述了多孔介質(zhì)內(nèi)部結(jié)構(gòu)和孔隙分布情況,是描述多孔介質(zhì)內(nèi)部輸運物理現(xiàn)象的重要公式。

    1.2 孔隙迂曲特征描述

    假設(shè)多孔介質(zhì)是由一束橫截面積不同的迂曲毛細管道組成,在該模型中,單根迂曲毛細管道的長度L(λ) 與管道直徑λ的關(guān)系滿足分形冪規(guī)律,表示為

    式中:L0為毛細管流動方向的特征長度,m。根據(jù)分形理論,L0可以通過L0=計算得到。須要注意的是L0不是多孔介質(zhì)的實際邊長[21]。Dt是迂曲分形維數(shù),Dt越高表示毛細管道越彎曲。在二維平面內(nèi),當Dt=1 時,表示毛細管道是直管道,當Dt=2 時,表示毛細管道非常彎曲,已經(jīng)填充了整個二維平面。根據(jù)式(6),迂曲度可以定義為

    2 多孔介質(zhì)中氣體擴散的分形模型

    氣體在多孔介質(zhì)中的擴散問題影響因素眾多,其中孔隙結(jié)構(gòu)為主要的影響因素。當孔隙直徑不同時,氣體擴散的規(guī)律也不同。Beskok 等[4]通過實驗給出了描述氣體在一維圓管中氣體體積流量的模型,該模型根據(jù)Hagen-Poiseuille 方程修正得到,其表達式為

    式中:kn為克努森數(shù),kn=l/λ,定義為氣體分子的平均自由程l與孔隙直徑λ的比值,l的單位為m;Δp為毛細管中的氣體壓力差,Pa;μ為氣體黏度,Pa·s。f(kn)為滲透率修正因子,此因子只與kn有關(guān),其表達式為

    根據(jù)kn的不同,氣體在致密多孔介質(zhì)中可以分為連續(xù)介質(zhì)區(qū)、滑移區(qū)、過渡區(qū)和自由分子區(qū)等4 個流動區(qū)域,但是,將式(9)和式(10)代入式(8)后得到的公式具有很強的非線性特征,不易求解,應(yīng)用價值不大[22]。為了得到廣泛適用的滲透率修正因子,使用matlab 多項式曲線擬合工具,對式(9)進行分段擬合。

    當kn≤0.001和0.001<kn≤0.100時,氣體流動處于連續(xù)介質(zhì)區(qū)和滑移區(qū)。運用泰勒公式將式(10)化簡后再代入式(9),然后忽略滲透率修正因子的二階和高階修正,保留一級修正。此時,f(kn)=1+4kn,據(jù)文獻[23]報道,該公式正是Klinkengerg 于1941 年給出的考慮滑脫效應(yīng)的氣測滲透率的數(shù)學表達式,這里的4kn即為滑移因子,但是忽略了二階和高級修正會導(dǎo)致簡化模型值小于Beskok模型值,而通過曲線擬合的方式得到的擬合公式f(kn)=1+4.8 kn其簡化模型值與Beskok 模型值比較吻合。相對誤差為0.060%~0.272%,度量擬合優(yōu)度的統(tǒng)計量r2可達0.999 4。

    當0.1 <kn≤10.0 和kn>10 時,氣體處于過渡區(qū)和自由分子區(qū)。過渡區(qū)的擬合公式為f(kn)=0.8+6.2kn,相對誤差為1.3%~8.5%,統(tǒng)計量r2為0.992。自由分子區(qū)的擬合公式為f(kn)=-4.1 +6.65kn,相對誤差為0.11%~1.93%,統(tǒng)計量r2為0.999。

    綜上,通過對滲透率修正因子曲線擬合而簡化得到的模型如下

    通過對式(8)的擬合得到分段函數(shù)與原函數(shù)的對比(圖1),可以看出簡化后的模型與原模型擬合誤差較小,精確度較高。通過對滲透率修正因子的簡化,不僅為后續(xù)理論計算和求解提供了基礎(chǔ),同時也有利于Beskok 模型的推廣。

    圖1 Beskok 模型與簡化模型的對比Fig.1 Comparison between Beskok model and simplified model

    所以一維圓管中氣體體積流量公式也可以用一個分段函數(shù)來描述

    由于多孔介質(zhì)具有不規(guī)則性和復(fù)雜性,在同一多孔介質(zhì)中孔隙大小差別很大,這使得多種流動區(qū)域在同一多孔介質(zhì)中同時存在,根據(jù)氣體流動區(qū)域的劃分,計算氣體通過橫截面的流量公式如下[17]

    式中:λk=0.1l為滑移區(qū)和過渡區(qū)的臨界直徑;λp=10l為過渡區(qū)和自由分子區(qū)的臨界直徑。

    聯(lián)合式(2)、式(5)、式(11)和式(12)可得

    根據(jù)菲克定律,氣體在多孔介質(zhì)中擴散的總的摩爾流率為[15]

    式中:QF為氣體摩爾流率,mol/s;At為多孔介質(zhì)總的截面積,m2;Deff為有效擴散系數(shù),m2/s;ΔC為多孔介質(zhì)孔道兩端的濃度差,ΔC=Δp/(RT),其中,R為氣體常數(shù)8.314J/(mol·K) ;T為溫度,K。計算時須要將氣體摩爾流率轉(zhuǎn)化為氣體體積流量。

    聯(lián)合式(5)、式(6)、式(13)和式(14),可以求出

    式(16)是基于滲透率修正因子、迂曲毛細管束模型、分形理論和菲克定律求出來的多孔介質(zhì)氣體擴散模型。此模型的每一個參數(shù)都具有物理意義,不包含經(jīng)驗常數(shù)。式(16)表明,氣體有效擴散系數(shù)與多孔介質(zhì)面積分形維數(shù)、迂曲分形維數(shù)、孔隙度、溫度、最大與最小孔隙直徑、氣體分子平均自由程、氣體摩爾質(zhì)量、氣體密度和黏度等均有關(guān),并且描述了在不同流動區(qū)域共同作用下多孔介質(zhì)中氣體擴散的詳細機理。

    通過式(16)可得知這些參數(shù)對氣體有效擴散系數(shù)的影響,并應(yīng)用于工程實踐中的相關(guān)參數(shù)分析。

    3 分析和討論

    本文建立的氣體有效擴散系數(shù)包含了2 類參數(shù)。一類是關(guān)于多孔介質(zhì)孔隙數(shù)目和結(jié)構(gòu)的參數(shù),如多孔介質(zhì)面積分形維數(shù)、迂曲分形維數(shù)、特征長度、孔隙度等。另一類是關(guān)于氣體的性質(zhì)參數(shù),如密度、摩爾質(zhì)量、黏度、氣體分子平均自由程等。

    多孔介質(zhì)的結(jié)構(gòu)參數(shù)通過分形理論的描述和計算得出,面積分形維數(shù)由以下公式來確定

    式中:DE為歐幾里德維數(shù),在二維空間中,DE=2;在三維空間中,DE=3。該文是在二維平面上建立的模型,所以DE=2。Feng 等[24]歸納了孔隙度和面積分形維數(shù)的實驗數(shù)據(jù),結(jié)果顯示當λmin/λmax≈10-3時,理論模型值與實驗結(jié)果比較吻合。

    對于迂曲分形維數(shù),尹帥等[25]分析了常規(guī)及非常規(guī)儲層巖石的分形特征。運用壓汞法和吸附法研究得到了具有分形特性的多孔介質(zhì)的迂曲分形維數(shù),其中砂巖多孔介質(zhì)的迂曲分形維數(shù)為1.207~1.704,頁巖多孔介質(zhì)的迂曲分形維數(shù)為1.668~2.187。

    然后將已有實驗數(shù)據(jù)與本文提出的有效擴散系數(shù)的理論模型進行對比。Currie[26]將H2作為擴散氣體,使用了多種不同的多孔介質(zhì)來進行擴散實驗。須要注意的是,實驗中一般使用體孔隙度,而理論模型中是用面孔隙度。須要借助關(guān)系式Φ=,把面孔隙度轉(zhuǎn)化為體孔隙度。其中=1+0.63ln(1/Φ),該公式由實驗數(shù)據(jù)擬合得到[1]。

    表1 為建立的模型的部分參數(shù):

    表1 分形模型中的部分參數(shù)與數(shù)值Table 1 Some parameters and values in the fractal model

    將以上數(shù)據(jù)和公式(17)代入式(16)可以得到圖2 中的函數(shù)關(guān)系,圖2 顯示氣體有效擴散系數(shù)的預(yù)測值隨著孔隙度的增加而增加,因為隨著孔隙度的增加,氣體擴散的區(qū)域和路徑會增多,有效擴散系數(shù)相應(yīng)增加。當孔隙度超過0.65 時,模型預(yù)測值會低于實際值,因為該模型中的最大孔隙直徑使用的是固定值。通過對比模型預(yù)測值與實驗數(shù)據(jù),可知該模型能夠準確地計算氣體在分形多孔介質(zhì)中的有效擴散系數(shù)。本次建立的理論模型可以用于多孔介質(zhì)微納米級孔隙情況下的單相氣體擴散計算,可以為非常規(guī)天然氣資源開采提供關(guān)鍵參數(shù),也可以為氣液兩相或者多相滲透與擴散提供參考。

    圖2 分形模型的氣體有效擴散系數(shù)與文獻[26]實驗數(shù)據(jù)對比Fig.2 Comparison of effective gas diffusion coefficient of fractal model with experimental data in reference

    圖3 為孔道迂曲分形維數(shù)分別取1.5,1.6,1.7時,氣體有效擴散系數(shù)隨著多孔介質(zhì)面積分形維數(shù)的變化趨勢。從圖3 可以得出,面積分形維數(shù)越大,氣體有效擴散系數(shù)也越大。因為面積分形維數(shù)增加時,孔隙數(shù)目也會增加,孔隙度也會增大,導(dǎo)致擴散相對更加容易。當面積分形維數(shù)一定時,迂曲分形維數(shù)增加,氣體有效擴散系數(shù)會減小。因為迂曲分形維數(shù)越大,表明孔道毛細管越彎曲,擴散阻力越大,導(dǎo)致擴散性能減小。

    圖3 面積分形維數(shù)和迂曲分形維數(shù)對氣體有效擴散系數(shù)的影響Fig.3 Effect of area fractal dimension and tortuous fractal dimension on effective gas diffusion coefficient

    圖4 為考慮滑移因子的模型與不考慮滑移因子模型的對比。據(jù)文獻[23]報道,Klinkengerg 于1941 年發(fā)現(xiàn)了滑移現(xiàn)象,具體表現(xiàn)為氣體的表觀滲透率明顯大于絕對滲透率。從圖4 可以看出當考慮滑移因子時,氣體的有效擴散系數(shù)大于未考慮滑移因子的氣體有效擴散系數(shù),這與Klinkengerg 描述的現(xiàn)象一致。由圖4 分析可知,在孔隙直徑較大的情況下,滑移因子的影響越來越小,當孔隙直徑較小時,滑移因子對有效擴散系數(shù)的影響較大,這與滑脫效應(yīng)一般發(fā)生在低滲透率和低孔隙度相一致。從圖4 還可以看出,考慮了滑移因子的分形模型與不同多孔介質(zhì)的實驗數(shù)據(jù)更加吻合,這表明考慮了滑移因子的分形模型能夠適用于更多種類的多孔介質(zhì)。

    圖4 滑移因子對氣體有效擴散系數(shù)的影響Fig.4 Effect of slip factor on effective gas diffusion coefficient

    圖5 為在不同最大孔隙直徑下,有效擴散系數(shù)隨孔隙度變化的情況。

    圖5 最大孔隙直徑對氣體有效擴散系數(shù)的影響Fig.5 Effect of maximum pore diameter on effective gas diffusion coefficient

    從圖5 可以得出,當最大孔隙直徑減小時,氣體有效擴散系數(shù)也會減小。因為根據(jù)分形理論,分形多孔介質(zhì)的面積分形維數(shù)、迂曲分形維數(shù)、最小最大孔徑比、擴散截面面積等都與最大孔隙直徑有關(guān),所以最大孔隙直徑是氣體有效擴散系數(shù)的主要影響因素。

    圖6 為使用不同的擴散氣體時,有效擴散系數(shù)隨孔隙度的變化。由圖6 可知,不同的氣體在多孔介質(zhì)中的有效擴散系數(shù)不同,當擴散的氣體的摩爾質(zhì)量,或者氣體密度增加時,有效擴散系數(shù)會減小。因為氣體質(zhì)量越大,擴散速率就會越小,有效擴散系數(shù)也會越小,這與氣體擴散定律描述的現(xiàn)象一致。

    圖6 不同氣體對氣體有效擴散系數(shù)的影響Fig.6 Effect of different kinds of gases on effective gas diffusion coefficient

    4 結(jié)論

    (1)基于多孔介質(zhì)的分形特征、迂曲毛細管束模型、菲克定律、滲透率修正因子建立了多孔介質(zhì)中氣體擴散的分形模型。該模型包含了多孔介質(zhì)面積分形維數(shù)、迂曲分形維數(shù)、孔隙度、最大孔隙直徑、氣體密度和摩爾質(zhì)量等參數(shù)。通過參數(shù)分析可以發(fā)現(xiàn):多孔介質(zhì)孔隙度、面積分形維數(shù)、最大孔隙直徑、氣體的物理性質(zhì)等對有效擴散系數(shù)均有顯著影響。

    (2)對比分形多孔介質(zhì)氣體擴散模型與Currie關(guān)于氣體擴散率的數(shù)據(jù)發(fā)現(xiàn):體孔隙度為0.10~0.65 時,模型與多種不同的多孔介質(zhì)實驗數(shù)據(jù)均具有較好的一致性,說明該模型可以用于計算具有分形特性的多孔介質(zhì)的氣體有效擴散系數(shù),對非常規(guī)天然氣在實際多孔介質(zhì)中的有效擴散系數(shù)的計算具有一定的參考意義。同時也能為后續(xù)研究多孔介質(zhì)氣液兩相或多相滲透與擴散提供參考。

    (3)考慮滑移因子的氣體有效擴散系數(shù)大于未考慮滑移因子的氣體有效擴散系數(shù),在較小的孔隙直徑下,兩者的差距比較明顯,隨著孔隙直徑的增大,兩者的差異逐漸減小,直至相等。考慮了滑移因子的氣體有效擴散系數(shù)與Currie 關(guān)于氣體擴散率的實驗數(shù)據(jù)更加吻合。

    猜你喜歡
    擴散系數(shù)維數(shù)分形
    β-變換中一致丟番圖逼近問題的維數(shù)理論
    感受分形
    一類齊次Moran集的上盒維數(shù)
    分形之美
    分形空間上廣義凸函數(shù)的新Simpson型不等式及應(yīng)用
    關(guān)于齊次Moran集的packing維數(shù)結(jié)果
    涉及相變問題Julia集的Hausdorff維數(shù)
    基于Sauer-Freise 方法的Co- Mn 體系fcc 相互擴散系數(shù)的研究
    上海金屬(2015年5期)2015-11-29 01:13:59
    FCC Ni-Cu 及Ni-Mn 合金互擴散系數(shù)測定
    上海金屬(2015年6期)2015-11-29 01:09:09
    非時齊擴散模型中擴散系數(shù)的局部估計
    亚洲欧美日韩无卡精品| 晚上一个人看的免费电影| 国产大屁股一区二区在线视频| 亚洲无线观看免费| 国产乱人视频| 国产亚洲午夜精品一区二区久久| 国产精品久久久久成人av| 久久久久久人妻| 久久精品国产亚洲av涩爱| 日韩一本色道免费dvd| 一级片'在线观看视频| 99久久综合免费| 欧美成人精品欧美一级黄| 久久亚洲国产成人精品v| 人人妻人人澡人人爽人人夜夜| 这个男人来自地球电影免费观看 | 你懂的网址亚洲精品在线观看| 日韩伦理黄色片| 欧美性感艳星| 国产精品国产三级国产专区5o| 日韩一区二区视频免费看| 免费观看av网站的网址| 日韩中文字幕视频在线看片 | 欧美三级亚洲精品| 国产欧美另类精品又又久久亚洲欧美| 成年人午夜在线观看视频| 日日啪夜夜爽| 亚洲电影在线观看av| 亚洲av成人精品一二三区| 国产黄片视频在线免费观看| 在线亚洲精品国产二区图片欧美 | 欧美变态另类bdsm刘玥| 我的老师免费观看完整版| 色吧在线观看| 国产精品嫩草影院av在线观看| 精品亚洲成a人片在线观看 | 亚洲成人一二三区av| 久热这里只有精品99| 国产免费福利视频在线观看| 精品久久国产蜜桃| 男人爽女人下面视频在线观看| 免费黄频网站在线观看国产| 两个人的视频大全免费| 国产视频首页在线观看| 亚洲av综合色区一区| kizo精华| 久久国产乱子免费精品| 久久韩国三级中文字幕| 老女人水多毛片| 黑丝袜美女国产一区| 日韩不卡一区二区三区视频在线| 国产精品久久久久久av不卡| 中文在线观看免费www的网站| 黄色日韩在线| 国产欧美另类精品又又久久亚洲欧美| 麻豆成人午夜福利视频| 国产午夜精品一二区理论片| 国产精品福利在线免费观看| 国产av国产精品国产| 欧美 日韩 精品 国产| 十八禁网站网址无遮挡 | 观看美女的网站| 建设人人有责人人尽责人人享有的 | 国产精品99久久久久久久久| 日韩不卡一区二区三区视频在线| av国产久精品久网站免费入址| 蜜桃久久精品国产亚洲av| 亚洲最大成人中文| 亚洲美女视频黄频| 久久影院123| 成人影院久久| 男人狂女人下面高潮的视频| 精品亚洲乱码少妇综合久久| 人妻制服诱惑在线中文字幕| 男人爽女人下面视频在线观看| 久久精品国产亚洲av天美| 91精品一卡2卡3卡4卡| 女的被弄到高潮叫床怎么办| 久久婷婷青草| 亚洲,欧美,日韩| 看十八女毛片水多多多| 国产精品.久久久| 1000部很黄的大片| 成人黄色视频免费在线看| 26uuu在线亚洲综合色| 99热网站在线观看| 男女国产视频网站| 插阴视频在线观看视频| 水蜜桃什么品种好| 国产欧美日韩一区二区三区在线 | 我的女老师完整版在线观看| 99久国产av精品国产电影| 日本欧美视频一区| 亚洲精品日韩在线中文字幕| 国产高清有码在线观看视频| 最近最新中文字幕大全电影3| 国产精品久久久久久久电影| 爱豆传媒免费全集在线观看| 99热6这里只有精品| 97超碰精品成人国产| 亚洲精品,欧美精品| 久久久久久久精品精品| 国产爽快片一区二区三区| 老司机影院毛片| 久久国产乱子免费精品| 午夜精品国产一区二区电影| 在线观看免费高清a一片| 亚洲高清免费不卡视频| 国产精品嫩草影院av在线观看| 毛片女人毛片| 欧美3d第一页| 一级爰片在线观看| 十八禁网站网址无遮挡 | 91精品一卡2卡3卡4卡| 日韩成人伦理影院| h日本视频在线播放| 久久久久久九九精品二区国产| 欧美bdsm另类| 大片电影免费在线观看免费| 日韩成人av中文字幕在线观看| 高清欧美精品videossex| 久久久久久久亚洲中文字幕| 国产精品偷伦视频观看了| 国产av国产精品国产| 97在线人人人人妻| 日韩一本色道免费dvd| 在线观看美女被高潮喷水网站| 国产精品伦人一区二区| 欧美三级亚洲精品| 综合色丁香网| 久久久a久久爽久久v久久| 亚洲精品aⅴ在线观看| 91久久精品国产一区二区三区| 成年免费大片在线观看| 美女脱内裤让男人舔精品视频| 99久久中文字幕三级久久日本| 日产精品乱码卡一卡2卡三| 麻豆精品久久久久久蜜桃| av黄色大香蕉| av福利片在线观看| 一区二区av电影网| 亚洲欧美日韩卡通动漫| 毛片一级片免费看久久久久| 蜜桃亚洲精品一区二区三区| 欧美变态另类bdsm刘玥| 久久国产精品大桥未久av | 国产高清不卡午夜福利| 一级a做视频免费观看| 麻豆乱淫一区二区| 成人亚洲欧美一区二区av| 777米奇影视久久| 日韩一区二区三区影片| 五月开心婷婷网| 菩萨蛮人人尽说江南好唐韦庄| 大香蕉97超碰在线| 又爽又黄a免费视频| 夫妻性生交免费视频一级片| 国产精品一及| 国产在视频线精品| 18+在线观看网站| 26uuu在线亚洲综合色| 国产精品久久久久成人av| 欧美高清成人免费视频www| 亚洲熟女精品中文字幕| 久久久久网色| 少妇高潮的动态图| 一级片'在线观看视频| 色视频www国产| 三级经典国产精品| 97在线视频观看| 久热久热在线精品观看| 欧美xxⅹ黑人| 综合色丁香网| 国产 一区 欧美 日韩| 亚洲精品亚洲一区二区| 亚洲高清免费不卡视频| 91狼人影院| 九九久久精品国产亚洲av麻豆| 国产精品久久久久久精品电影小说 | 国产成人免费观看mmmm| 国产视频内射| 亚洲电影在线观看av| 成人18禁高潮啪啪吃奶动态图 | 最后的刺客免费高清国语| 一级毛片aaaaaa免费看小| 黄片无遮挡物在线观看| 国产在视频线精品| 尾随美女入室| 精品午夜福利在线看| 三级国产精品片| 免费av不卡在线播放| 欧美少妇被猛烈插入视频| 欧美 日韩 精品 国产| 亚洲欧美一区二区三区黑人 | 高清日韩中文字幕在线| 一级二级三级毛片免费看| 国产久久久一区二区三区| 身体一侧抽搐| 精品久久久久久久末码| 久久精品国产亚洲av涩爱| 黑人高潮一二区| 青春草国产在线视频| 色网站视频免费| 97在线视频观看| 亚洲精品乱码久久久v下载方式| 人妻夜夜爽99麻豆av| 嫩草影院新地址| 亚洲激情五月婷婷啪啪| 激情 狠狠 欧美| 97在线人人人人妻| 精品久久久久久久久亚洲| 久久国内精品自在自线图片| 九草在线视频观看| 亚洲国产欧美人成| 亚洲国产毛片av蜜桃av| 视频区图区小说| 国产 精品1| 国产毛片在线视频| 在线观看一区二区三区| 中国国产av一级| 我的老师免费观看完整版| 国产欧美另类精品又又久久亚洲欧美| 国产片特级美女逼逼视频| 久久这里有精品视频免费| 男女边摸边吃奶| 成人18禁高潮啪啪吃奶动态图 | 一级毛片aaaaaa免费看小| 成年美女黄网站色视频大全免费 | 国产伦在线观看视频一区| 边亲边吃奶的免费视频| 精品一区在线观看国产| 51国产日韩欧美| av在线老鸭窝| 久久国产精品大桥未久av | 观看av在线不卡| 大片电影免费在线观看免费| 夜夜骑夜夜射夜夜干| av在线app专区| 五月开心婷婷网| 少妇被粗大猛烈的视频| 精品少妇久久久久久888优播| 亚洲成人一二三区av| 伦精品一区二区三区| 2022亚洲国产成人精品| 黑人猛操日本美女一级片| 久久国内精品自在自线图片| 国产探花极品一区二区| 国产精品三级大全| 涩涩av久久男人的天堂| 日本与韩国留学比较| 国产精品嫩草影院av在线观看| 少妇高潮的动态图| 国产亚洲最大av| 精品亚洲成国产av| 草草在线视频免费看| 最近的中文字幕免费完整| 最近中文字幕2019免费版| 深爱激情五月婷婷| 国产伦精品一区二区三区四那| 久久久成人免费电影| 性色av一级| 永久免费av网站大全| 在线观看免费高清a一片| 国产精品人妻久久久久久| 日日撸夜夜添| 亚洲精品一区蜜桃| 免费播放大片免费观看视频在线观看| 男人狂女人下面高潮的视频| 欧美亚洲 丝袜 人妻 在线| 男人和女人高潮做爰伦理| 九草在线视频观看| 校园人妻丝袜中文字幕| 国产伦在线观看视频一区| 日本wwww免费看| 搡老乐熟女国产| 男人舔奶头视频| 日韩大片免费观看网站| 黄色视频在线播放观看不卡| 亚洲性久久影院| 欧美日韩精品成人综合77777| 男男h啪啪无遮挡| 高清毛片免费看| 免费大片黄手机在线观看| 精品国产乱码久久久久久小说| 精品久久久久久久末码| 久久久a久久爽久久v久久| 熟女电影av网| 亚洲一级一片aⅴ在线观看| 精品久久久精品久久久| 亚洲图色成人| 大陆偷拍与自拍| 成人国产麻豆网| 久久国产亚洲av麻豆专区| 国产精品蜜桃在线观看| 国产在线视频一区二区| 高清日韩中文字幕在线| 国产在线一区二区三区精| 国产精品福利在线免费观看| 国产视频内射| 一区二区三区精品91| 国产人妻一区二区三区在| 大陆偷拍与自拍| 亚洲aⅴ乱码一区二区在线播放| 色5月婷婷丁香| 少妇的逼好多水| av在线app专区| 最新中文字幕久久久久| 极品教师在线视频| 日韩中字成人| 欧美日韩亚洲高清精品| 少妇人妻一区二区三区视频| 免费看av在线观看网站| 亚洲av国产av综合av卡| 国产人妻一区二区三区在| 精品熟女少妇av免费看| 亚洲av免费高清在线观看| 最近中文字幕2019免费版| 久久韩国三级中文字幕| 国产精品一区二区三区四区免费观看| 久久青草综合色| 午夜视频国产福利| 26uuu在线亚洲综合色| 在线观看国产h片| 一本色道久久久久久精品综合| 久久久久精品性色| 日日摸夜夜添夜夜添av毛片| 精品人妻偷拍中文字幕| 91久久精品国产一区二区成人| 人妻 亚洲 视频| 国产精品.久久久| 精品国产三级普通话版| 免费在线观看成人毛片| 久久久精品免费免费高清| av天堂中文字幕网| 国产精品偷伦视频观看了| 亚洲精品日韩av片在线观看| 色吧在线观看| 国产淫片久久久久久久久| 在线精品无人区一区二区三 | 三级国产精品片| 成人免费观看视频高清| 国产视频首页在线观看| 欧美极品一区二区三区四区| 久久久久久久久久人人人人人人| 免费久久久久久久精品成人欧美视频 | 少妇丰满av| 熟女人妻精品中文字幕| 亚洲最大成人中文| 国产一区二区三区av在线| 国精品久久久久久国模美| 欧美一级a爱片免费观看看| 久久国产乱子免费精品| 简卡轻食公司| 久久鲁丝午夜福利片| 黑人猛操日本美女一级片| 日韩免费高清中文字幕av| 国产成人一区二区在线| 亚洲欧洲国产日韩| 乱系列少妇在线播放| 欧美人与善性xxx| 色婷婷av一区二区三区视频| 99久国产av精品国产电影| 男人添女人高潮全过程视频| 久久99蜜桃精品久久| 日韩欧美精品免费久久| 一级爰片在线观看| 国产精品偷伦视频观看了| 六月丁香七月| 久久久久久人妻| 日韩av在线免费看完整版不卡| 啦啦啦中文免费视频观看日本| 国产欧美日韩精品一区二区| 免费高清在线观看视频在线观看| 美女xxoo啪啪120秒动态图| 在线观看人妻少妇| 黄片无遮挡物在线观看| 国内精品宾馆在线| 国产探花极品一区二区| 国产精品免费大片| 久久久久久人妻| 午夜福利高清视频| 国语对白做爰xxxⅹ性视频网站| 亚洲三级黄色毛片| 五月伊人婷婷丁香| 熟女人妻精品中文字幕| 岛国毛片在线播放| 少妇的逼好多水| 狠狠精品人妻久久久久久综合| 亚洲一区二区三区欧美精品| 国产精品嫩草影院av在线观看| 久久女婷五月综合色啪小说| 超碰av人人做人人爽久久| 成人午夜精彩视频在线观看| 精品一区二区免费观看| 久久精品国产亚洲网站| 欧美激情极品国产一区二区三区 | 啦啦啦视频在线资源免费观看| 国产欧美日韩精品一区二区| 国产一区亚洲一区在线观看| 18禁动态无遮挡网站| 大话2 男鬼变身卡| 好男人视频免费观看在线| 黄色一级大片看看| 我要看日韩黄色一级片| 在线观看人妻少妇| 日本色播在线视频| 色视频www国产| 亚洲精品自拍成人| 97热精品久久久久久| 一个人看视频在线观看www免费| 国产一区亚洲一区在线观看| 中文字幕制服av| 亚洲精品乱码久久久久久按摩| 精品久久久精品久久久| 久久精品熟女亚洲av麻豆精品| 99热国产这里只有精品6| 国产淫片久久久久久久久| 欧美激情国产日韩精品一区| 18禁裸乳无遮挡动漫免费视频| 欧美3d第一页| 亚洲欧美一区二区三区黑人 | 精华霜和精华液先用哪个| 97热精品久久久久久| 欧美日本视频| 菩萨蛮人人尽说江南好唐韦庄| 国产视频内射| .国产精品久久| 一级毛片久久久久久久久女| 特大巨黑吊av在线直播| 午夜福利网站1000一区二区三区| 少妇猛男粗大的猛烈进出视频| 亚洲欧美成人精品一区二区| a级一级毛片免费在线观看| 国产精品偷伦视频观看了| 国产精品嫩草影院av在线观看| 女人久久www免费人成看片| 亚洲欧洲日产国产| 高清在线视频一区二区三区| 哪个播放器可以免费观看大片| 久久韩国三级中文字幕| 秋霞伦理黄片| 精品国产乱码久久久久久小说| 久久久久国产精品人妻一区二区| 国产精品.久久久| 久久影院123| 在线观看人妻少妇| 精品久久久久久久久亚洲| 久久久久国产精品人妻一区二区| 网址你懂的国产日韩在线| 国产亚洲午夜精品一区二区久久| 新久久久久国产一级毛片| 黄片无遮挡物在线观看| 2021少妇久久久久久久久久久| 国产探花极品一区二区| 成人影院久久| 亚洲在久久综合| 亚洲高清免费不卡视频| 不卡视频在线观看欧美| 亚洲综合色惰| 少妇裸体淫交视频免费看高清| av在线app专区| 国产精品福利在线免费观看| 亚洲人与动物交配视频| 亚洲欧美一区二区三区黑人 | 国产白丝娇喘喷水9色精品| 亚洲国产欧美人成| 免费少妇av软件| 99久久精品热视频| 中文字幕免费在线视频6| 国产真实伦视频高清在线观看| 国产黄色视频一区二区在线观看| 一级a做视频免费观看| 老师上课跳d突然被开到最大视频| 搡女人真爽免费视频火全软件| 91精品国产九色| 亚洲欧美精品专区久久| 免费在线观看成人毛片| 久久精品人妻少妇| 国产成人aa在线观看| 男的添女的下面高潮视频| 亚洲欧美清纯卡通| .国产精品久久| 少妇猛男粗大的猛烈进出视频| 国产成人freesex在线| 一级二级三级毛片免费看| 成人午夜精彩视频在线观看| 国产精品免费大片| 另类亚洲欧美激情| 在线观看免费高清a一片| 日本黄大片高清| 99热全是精品| 亚洲伊人久久精品综合| 久久婷婷青草| 久久99热这里只频精品6学生| 精品国产露脸久久av麻豆| 中文乱码字字幕精品一区二区三区| 国产精品福利在线免费观看| 国产日韩欧美亚洲二区| 国产高清国产精品国产三级 | 国产精品熟女久久久久浪| 亚洲精品久久久久久婷婷小说| 欧美日韩亚洲高清精品| 欧美精品国产亚洲| 免费看光身美女| 大香蕉97超碰在线| 欧美日韩在线观看h| 精品亚洲成国产av| 成人漫画全彩无遮挡| 久久国产精品大桥未久av | 日韩在线高清观看一区二区三区| 久久国产亚洲av麻豆专区| 成人毛片60女人毛片免费| 亚洲av综合色区一区| 日本vs欧美在线观看视频 | 亚洲av国产av综合av卡| 丝瓜视频免费看黄片| 亚洲国产欧美人成| 99re6热这里在线精品视频| 校园人妻丝袜中文字幕| 久久久欧美国产精品| 午夜激情久久久久久久| 又黄又爽又刺激的免费视频.| 视频区图区小说| 乱系列少妇在线播放| 亚洲欧美精品自产自拍| 在线观看免费高清a一片| 18+在线观看网站| 免费不卡的大黄色大毛片视频在线观看| 午夜精品国产一区二区电影| 中文欧美无线码| a 毛片基地| 免费av不卡在线播放| 国产视频首页在线观看| 久久国产精品男人的天堂亚洲 | 久久99热这里只有精品18| 狂野欧美白嫩少妇大欣赏| 亚洲精品乱码久久久久久按摩| 亚洲高清免费不卡视频| 大陆偷拍与自拍| 男人和女人高潮做爰伦理| 国产永久视频网站| 亚洲人与动物交配视频| 男女边摸边吃奶| 欧美3d第一页| 亚洲怡红院男人天堂| 国产色爽女视频免费观看| 九色成人免费人妻av| 亚洲精品一二三| 天美传媒精品一区二区| 波野结衣二区三区在线| 国产免费视频播放在线视频| 看免费成人av毛片| 久久6这里有精品| 波野结衣二区三区在线| 插阴视频在线观看视频| 免费av不卡在线播放| 夫妻午夜视频| 中文字幕人妻熟人妻熟丝袜美| 国产一区亚洲一区在线观看| 一级爰片在线观看| 一级毛片 在线播放| 日日啪夜夜爽| www.色视频.com| 久久久久久人妻| 国产精品av视频在线免费观看| 18禁动态无遮挡网站| 在线看a的网站| 春色校园在线视频观看| 人体艺术视频欧美日本| 天天躁日日操中文字幕| 国产精品.久久久| 亚洲精品久久午夜乱码| 午夜免费男女啪啪视频观看| 日本猛色少妇xxxxx猛交久久| 久久久久久久久久久丰满| 亚洲天堂av无毛| 国产精品伦人一区二区| 精品一区二区免费观看| 乱码一卡2卡4卡精品| 一级片'在线观看视频| av专区在线播放| 色哟哟·www| 国产日韩欧美在线精品| 亚州av有码| 亚洲精华国产精华液的使用体验| 多毛熟女@视频| 日韩欧美一区视频在线观看 | 亚洲婷婷狠狠爱综合网| 久久99精品国语久久久| 欧美97在线视频| 亚洲成人手机| 国产极品天堂在线| 一个人免费看片子| 国产高清有码在线观看视频| 亚洲激情五月婷婷啪啪| 国产色爽女视频免费观看| 久久97久久精品| 亚洲av国产av综合av卡| 国产色爽女视频免费观看| 久久6这里有精品| 午夜福利在线在线| 国产免费又黄又爽又色| 国产精品久久久久久av不卡| 午夜福利在线在线| 久久国产精品男人的天堂亚洲 | 久久热精品热| 在线亚洲精品国产二区图片欧美 | 国产男人的电影天堂91| 国产高清三级在线| 各种免费的搞黄视频| 国产精品人妻久久久久久| 亚洲高清免费不卡视频| av网站免费在线观看视频| 久久人人爽人人片av| 国产精品国产三级专区第一集|