• <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ù)的局部估計
    亚洲美女搞黄在线观看| 你懂的网址亚洲精品在线观看| 欧美激情 高清一区二区三区| 国产成人精品福利久久| 亚洲色图 男人天堂 中文字幕 | 国产免费视频播放在线视频| 亚洲国产精品成人久久小说| 中文字幕免费在线视频6| 80岁老熟妇乱子伦牲交| 国产av码专区亚洲av| 人人妻人人澡人人爽人人夜夜| 亚洲精品456在线播放app| 精品午夜福利在线看| 在线天堂最新版资源| 90打野战视频偷拍视频| 大话2 男鬼变身卡| 亚洲欧美一区二区三区黑人 | 少妇猛男粗大的猛烈进出视频| 少妇人妻精品综合一区二区| 国产亚洲av片在线观看秒播厂| 欧美激情极品国产一区二区三区 | 狠狠精品人妻久久久久久综合| 两个人免费观看高清视频| 香蕉丝袜av| 免费观看a级毛片全部| 大片电影免费在线观看免费| 成人无遮挡网站| 超碰97精品在线观看| 80岁老熟妇乱子伦牲交| 午夜福利网站1000一区二区三区| 日本vs欧美在线观看视频| av国产精品久久久久影院| 国产精品一二三区在线看| 97超碰精品成人国产| 精品人妻偷拍中文字幕| 中文精品一卡2卡3卡4更新| 王馨瑶露胸无遮挡在线观看| 黄色 视频免费看| av视频免费观看在线观看| 熟女电影av网| 精品一区二区免费观看| 综合色丁香网| 一级毛片我不卡| 大香蕉久久成人网| 18+在线观看网站| 国产日韩欧美在线精品| 看十八女毛片水多多多| av片东京热男人的天堂| av电影中文网址| 亚洲伊人色综图| 亚洲国产日韩一区二区| 最新中文字幕久久久久| 久久午夜福利片| 女人被躁到高潮嗷嗷叫费观| 亚洲国产精品999| 亚洲精品视频女| 亚洲av福利一区| 一本大道久久a久久精品| tube8黄色片| 久久国产精品大桥未久av| 亚洲美女视频黄频| 色婷婷久久久亚洲欧美| 男人添女人高潮全过程视频| 在线看a的网站| 国产一区亚洲一区在线观看| 大香蕉97超碰在线| 久久精品aⅴ一区二区三区四区 | 午夜福利在线观看免费完整高清在| 亚洲一级一片aⅴ在线观看| 一区二区三区乱码不卡18| 男人舔女人的私密视频| 中文字幕另类日韩欧美亚洲嫩草| 十分钟在线观看高清视频www| 最后的刺客免费高清国语| 伊人久久国产一区二区| 男女免费视频国产| 夫妻午夜视频| 成人二区视频| 久久国产精品大桥未久av| av免费在线看不卡| 嫩草影院入口| 国产精品麻豆人妻色哟哟久久| 99国产综合亚洲精品| 免费播放大片免费观看视频在线观看| 亚洲欧美清纯卡通| 国产日韩欧美在线精品| 欧美精品人与动牲交sv欧美| 日韩av免费高清视频| 五月天丁香电影| 国产精品久久久av美女十八| 中文字幕人妻丝袜制服| 色哟哟·www| 亚洲国产精品成人久久小说| 久热这里只有精品99| 亚洲国产精品一区二区三区在线| 色婷婷av一区二区三区视频| 国语对白做爰xxxⅹ性视频网站| 色婷婷av一区二区三区视频| 成人影院久久| 久久99热6这里只有精品| a级片在线免费高清观看视频| 日韩一本色道免费dvd| 国产成人免费观看mmmm| 午夜福利视频精品| 涩涩av久久男人的天堂| 你懂的网址亚洲精品在线观看| 一级毛片 在线播放| 免费观看性生交大片5| 日韩一区二区三区影片| 久久亚洲国产成人精品v| 欧美精品一区二区免费开放| 久久青草综合色| 熟妇人妻不卡中文字幕| 色网站视频免费| 精品一区二区三区四区五区乱码 | av卡一久久| 国产精品.久久久| 性高湖久久久久久久久免费观看| av网站免费在线观看视频| 亚洲av成人精品一二三区| 免费av中文字幕在线| 啦啦啦中文免费视频观看日本| 免费观看性生交大片5| 18禁在线无遮挡免费观看视频| 久久精品国产亚洲av天美| 久久精品人人爽人人爽视色| 26uuu在线亚洲综合色| 免费观看无遮挡的男女| 精品一品国产午夜福利视频| 亚洲精品视频女| 久久影院123| 激情视频va一区二区三区| 亚洲,欧美精品.| 两性夫妻黄色片 | 婷婷色av中文字幕| 欧美日本中文国产一区发布| 观看av在线不卡| 国产1区2区3区精品| 天天躁夜夜躁狠狠躁躁| 天天操日日干夜夜撸| 自拍欧美九色日韩亚洲蝌蚪91| 9热在线视频观看99| 18+在线观看网站| 亚洲熟女精品中文字幕| 久久韩国三级中文字幕| 国产乱人偷精品视频| 最近中文字幕2019免费版| 各种免费的搞黄视频| a级毛片黄视频| 精品国产一区二区三区四区第35| 国产一区二区三区综合在线观看 | 亚洲精品国产av蜜桃| 国产免费又黄又爽又色| 18在线观看网站| 久久青草综合色| 亚洲少妇的诱惑av| 久久久久久久久久久久大奶| 免费久久久久久久精品成人欧美视频 | 一区二区三区乱码不卡18| 少妇精品久久久久久久| 日本黄大片高清| 国产精品三级大全| 日韩制服丝袜自拍偷拍| 天堂中文最新版在线下载| 91久久精品国产一区二区三区| 另类精品久久| 亚洲成av片中文字幕在线观看 | 久久久久精品人妻al黑| 少妇的逼水好多| www.av在线官网国产| 91午夜精品亚洲一区二区三区| 亚洲伊人久久精品综合| 精品少妇久久久久久888优播| 欧美日韩av久久| 欧美精品高潮呻吟av久久| 成人亚洲精品一区在线观看| 18在线观看网站| 狂野欧美激情性xxxx在线观看| 亚洲四区av| 自线自在国产av| 国产成人精品福利久久| 亚洲精品自拍成人| 九色成人免费人妻av| 高清欧美精品videossex| 午夜精品国产一区二区电影| 日本-黄色视频高清免费观看| 亚洲欧洲精品一区二区精品久久久 | 26uuu在线亚洲综合色| 亚洲欧美清纯卡通| 国产在线免费精品| 日韩av免费高清视频| 亚洲精品第二区| 亚洲中文av在线| 国产精品人妻久久久影院| 麻豆乱淫一区二区| 日本黄大片高清| 黄色怎么调成土黄色| 国产亚洲精品第一综合不卡 | 天天躁夜夜躁狠狠躁躁| av网站免费在线观看视频| 黑人欧美特级aaaaaa片| 天天影视国产精品| 日韩伦理黄色片| 色婷婷久久久亚洲欧美| videosex国产| 亚洲av电影在线观看一区二区三区| av线在线观看网站| 黄片无遮挡物在线观看| 日韩欧美精品免费久久| 亚洲精品美女久久av网站| 啦啦啦在线观看免费高清www| 热re99久久精品国产66热6| 国产av码专区亚洲av| 久久热在线av| 精品一区二区三卡| 日韩精品有码人妻一区| 最近2019中文字幕mv第一页| 久久久久国产网址| www.熟女人妻精品国产 | 国产国拍精品亚洲av在线观看| 国产成人欧美| 欧美xxⅹ黑人| 精品国产国语对白av| 成人午夜精彩视频在线观看| 街头女战士在线观看网站| 99久久综合免费| 欧美变态另类bdsm刘玥| 爱豆传媒免费全集在线观看| 97精品久久久久久久久久精品| 免费高清在线观看日韩| 色视频在线一区二区三区| 国产一区有黄有色的免费视频| 在线观看免费日韩欧美大片| 99久久精品国产国产毛片| 深夜精品福利| 国产男人的电影天堂91| 精品亚洲成a人片在线观看| 大香蕉久久成人网| 最近的中文字幕免费完整| 久久精品人人爽人人爽视色| 日本猛色少妇xxxxx猛交久久| 成年av动漫网址| 中文天堂在线官网| 欧美 亚洲 国产 日韩一| 深夜精品福利| 狂野欧美激情性xxxx在线观看| 国产亚洲av片在线观看秒播厂| 中文欧美无线码| 中文乱码字字幕精品一区二区三区| 国产在线免费精品| 人人妻人人添人人爽欧美一区卜| 日本欧美视频一区| 久久这里只有精品19| 日韩一区二区三区影片| 伊人亚洲综合成人网| 精品少妇内射三级| 十八禁网站网址无遮挡| 久久热在线av| 18+在线观看网站| 黑人高潮一二区| 飞空精品影院首页| 日韩av不卡免费在线播放| 精品久久久精品久久久| 人人妻人人澡人人爽人人夜夜| 9色porny在线观看| kizo精华| 九色成人免费人妻av| 久久精品国产鲁丝片午夜精品| 中文欧美无线码| 国产 精品1| 久久久久久久久久人人人人人人| 中国国产av一级| 999精品在线视频| 国产精品人妻久久久久久| 国产成人免费观看mmmm| 麻豆精品久久久久久蜜桃| 日韩在线高清观看一区二区三区| 国产一区二区三区av在线| 色哟哟·www| 日本与韩国留学比较| 欧美人与善性xxx| 亚洲国产看品久久| av在线播放精品| 亚洲内射少妇av| 最近最新中文字幕免费大全7| 久久狼人影院| 亚洲综合色惰| 久久99精品国语久久久| 亚洲精品自拍成人| 在线看a的网站| 寂寞人妻少妇视频99o| 成人亚洲精品一区在线观看| 人妻一区二区av| 午夜免费观看性视频| av在线app专区| 老女人水多毛片| 香蕉精品网在线| 亚洲精品中文字幕在线视频| 一二三四在线观看免费中文在 | 亚洲一级一片aⅴ在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 香蕉国产在线看| 久久99蜜桃精品久久| 美女大奶头黄色视频| 亚洲精品aⅴ在线观看| 精品国产一区二区三区久久久樱花| 18禁国产床啪视频网站| 午夜免费男女啪啪视频观看| 亚洲中文av在线| 亚洲精品乱码久久久久久按摩| 爱豆传媒免费全集在线观看| 久久青草综合色| 亚洲欧美日韩另类电影网站| 精品亚洲乱码少妇综合久久| 亚洲在久久综合| 久久免费观看电影| 国产成人91sexporn| 亚洲精品久久午夜乱码| 欧美国产精品va在线观看不卡| 激情视频va一区二区三区| 国产xxxxx性猛交| 国产一区二区三区综合在线观看 | 久久国产亚洲av麻豆专区| 黄色 视频免费看| 精品一品国产午夜福利视频| 国产有黄有色有爽视频| 日本av免费视频播放| √禁漫天堂资源中文www| 咕卡用的链子| 日韩制服丝袜自拍偷拍| 伦精品一区二区三区| 亚洲一区二区三区欧美精品| 精品久久蜜臀av无| 性色avwww在线观看| 久久精品国产鲁丝片午夜精品| 在线天堂中文资源库| 日本色播在线视频| 两个人看的免费小视频| 国产老妇伦熟女老妇高清| 午夜影院在线不卡| 精品人妻偷拍中文字幕| 午夜久久久在线观看| 国产亚洲精品久久久com| 丝袜喷水一区| 日本免费在线观看一区| 最近最新中文字幕免费大全7| 国产精品熟女久久久久浪| 99久久人妻综合| 国产一区二区激情短视频 | 757午夜福利合集在线观看| 精品人妻1区二区| 亚洲人成77777在线视频| 一级作爱视频免费观看| 亚洲欧美日韩高清在线视频| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久久av美女十八| 成在线人永久免费视频| 女人被狂操c到高潮| 天天躁日日躁夜夜躁夜夜| 国产日韩欧美亚洲二区| 国内毛片毛片毛片毛片毛片| 18禁裸乳无遮挡动漫免费视频| 99国产精品99久久久久| 亚洲av成人不卡在线观看播放网| 亚洲熟妇熟女久久| 欧美色视频一区免费| 欧美最黄视频在线播放免费 | 精品国内亚洲2022精品成人 | 国产蜜桃级精品一区二区三区 | 757午夜福利合集在线观看| 久久亚洲精品不卡| av中文乱码字幕在线| 欧美另类亚洲清纯唯美| 天天影视国产精品| 亚洲av美国av| 脱女人内裤的视频| 一级毛片高清免费大全| 建设人人有责人人尽责人人享有的| 午夜福利欧美成人| 亚洲一码二码三码区别大吗| 18禁裸乳无遮挡动漫免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 香蕉国产在线看| 成年女人毛片免费观看观看9 | 久久精品国产a三级三级三级| 亚洲九九香蕉| 精品午夜福利视频在线观看一区| 亚洲va日本ⅴa欧美va伊人久久| 男人舔女人的私密视频| 成人免费观看视频高清| 制服诱惑二区| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av成人一区二区三| 91字幕亚洲| 国产成人精品在线电影| 可以免费在线观看a视频的电影网站| 亚洲av片天天在线观看| 亚洲av成人不卡在线观看播放网| 成人亚洲精品一区在线观看| 女人被狂操c到高潮| 欧美国产精品va在线观看不卡| 国产高清激情床上av| 制服诱惑二区| 午夜成年电影在线免费观看| 麻豆成人av在线观看| 欧美黑人欧美精品刺激| 久9热在线精品视频| 国产精品电影一区二区三区 | 交换朋友夫妻互换小说| 欧美精品av麻豆av| 国产无遮挡羞羞视频在线观看| 国产一区二区三区综合在线观看| 叶爱在线成人免费视频播放| 狂野欧美激情性xxxx| 午夜福利影视在线免费观看| 久久亚洲精品不卡| 日韩大码丰满熟妇| 黑人猛操日本美女一级片| 欧美日韩精品网址| 在线免费观看的www视频| 人人妻人人澡人人看| 一本综合久久免费| 午夜福利欧美成人| a在线观看视频网站| av天堂久久9| 久久久精品区二区三区| 国产精品自产拍在线观看55亚洲 | 黑人巨大精品欧美一区二区蜜桃| 亚洲av电影在线进入| 色综合欧美亚洲国产小说| 无限看片的www在线观看| 在线视频色国产色| 夜夜爽天天搞| 成人18禁高潮啪啪吃奶动态图| 久久久久久久国产电影| 国产精品自产拍在线观看55亚洲 | 亚洲成人免费av在线播放| 1024视频免费在线观看| 午夜免费鲁丝| 国产欧美日韩一区二区三区在线| 久久精品国产亚洲av高清一级| 国产精品久久电影中文字幕 | 久久国产精品男人的天堂亚洲| 国产精品久久视频播放| √禁漫天堂资源中文www| 亚洲综合色网址| 日本wwww免费看| 成人永久免费在线观看视频| 精品一品国产午夜福利视频| 欧美色视频一区免费| 又黄又粗又硬又大视频| 精品国产美女av久久久久小说| 国产国语露脸激情在线看| 岛国在线观看网站| 99精品欧美一区二区三区四区| 国产成人一区二区三区免费视频网站| 欧美精品亚洲一区二区| 亚洲熟妇熟女久久| 大码成人一级视频| 999精品在线视频| 国内毛片毛片毛片毛片毛片| 丰满人妻熟妇乱又伦精品不卡| 超碰97精品在线观看| 久久国产精品大桥未久av| 欧美大码av| 成人国产一区最新在线观看| 亚洲欧美一区二区三区黑人| 久久国产精品人妻蜜桃| 国产在线一区二区三区精| 最近最新免费中文字幕在线| 黄色毛片三级朝国网站| 久久亚洲精品不卡| 欧美亚洲 丝袜 人妻 在线| 黄色怎么调成土黄色| 欧美+亚洲+日韩+国产| 亚洲一码二码三码区别大吗| 亚洲欧美日韩高清在线视频| 亚洲欧美色中文字幕在线| 在线观看免费午夜福利视频| 一边摸一边抽搐一进一出视频| 人人妻人人澡人人爽人人夜夜| 久久人人爽av亚洲精品天堂| 国产蜜桃级精品一区二区三区 | 十八禁高潮呻吟视频| 在线观看日韩欧美| 天天操日日干夜夜撸| 欧美最黄视频在线播放免费 | 黄色成人免费大全| 国产有黄有色有爽视频| 黄色毛片三级朝国网站| 丰满的人妻完整版| 久久久久久亚洲精品国产蜜桃av| 国产一区二区三区视频了| 视频在线观看一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利在线免费观看网站| 亚洲精品久久午夜乱码| 欧美激情 高清一区二区三区| 大片电影免费在线观看免费| 亚洲少妇的诱惑av| 精品国产一区二区久久| 飞空精品影院首页| 免费看十八禁软件| netflix在线观看网站| 色94色欧美一区二区| 色综合欧美亚洲国产小说| 国产成人系列免费观看| 久久久精品国产亚洲av高清涩受| 老汉色∧v一级毛片| 十八禁人妻一区二区| 国产精品免费视频内射| 交换朋友夫妻互换小说| 99国产精品一区二区三区| 国产一区二区激情短视频| 国产激情欧美一区二区| 国产蜜桃级精品一区二区三区 | 自线自在国产av| 99精品久久久久人妻精品| 精品国产超薄肉色丝袜足j| 亚洲一码二码三码区别大吗| 搡老熟女国产l中国老女人| 人人妻人人澡人人爽人人夜夜| 亚洲精品av麻豆狂野| 宅男免费午夜| √禁漫天堂资源中文www| 巨乳人妻的诱惑在线观看| 国产精品亚洲av一区麻豆| 嫩草影视91久久| 日韩有码中文字幕| 大型黄色视频在线免费观看| 69精品国产乱码久久久| 婷婷成人精品国产| 国产精华一区二区三区| 麻豆乱淫一区二区| 99riav亚洲国产免费| 欧美日韩精品网址| 老鸭窝网址在线观看| 美女国产高潮福利片在线看| 亚洲国产精品一区二区三区在线| 美国免费a级毛片| 久久人妻福利社区极品人妻图片| 国产成人精品久久二区二区免费| 人人妻人人爽人人添夜夜欢视频| 黄色丝袜av网址大全| 在线永久观看黄色视频| 国产精品久久电影中文字幕 | a级毛片在线看网站| 一级a爱视频在线免费观看| 色精品久久人妻99蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 免费人成视频x8x8入口观看| 色婷婷久久久亚洲欧美| 亚洲av欧美aⅴ国产| 成年版毛片免费区| 亚洲av熟女| 女性被躁到高潮视频| 久久亚洲真实| 一级黄色大片毛片| 中亚洲国语对白在线视频| 欧美日韩精品网址| 亚洲av欧美aⅴ国产| 欧美精品啪啪一区二区三区| 亚洲国产精品一区二区三区在线| 老司机亚洲免费影院| 男女免费视频国产| 色婷婷久久久亚洲欧美| 亚洲黑人精品在线| 久久久久国产精品人妻aⅴ院 | 午夜福利影视在线免费观看| 欧美激情高清一区二区三区| 超碰成人久久| 亚洲精品中文字幕一二三四区| 激情在线观看视频在线高清 | 757午夜福利合集在线观看| 欧美精品人与动牲交sv欧美| 久热爱精品视频在线9| 久久久久久久久久久久大奶| 亚洲精品一卡2卡三卡4卡5卡| 一边摸一边抽搐一进一小说 | 亚洲av第一区精品v没综合| 亚洲自偷自拍图片 自拍| 久久精品国产亚洲av高清一级| 亚洲国产精品合色在线| 亚洲aⅴ乱码一区二区在线播放 | 国产极品粉嫩免费观看在线| 久久久久久人人人人人| 国产精品美女特级片免费视频播放器 | 亚洲午夜理论影院| 国产一区二区三区综合在线观看| 啦啦啦视频在线资源免费观看| 老司机影院毛片| av不卡在线播放| 国产成人精品久久二区二区91| 亚洲精品一卡2卡三卡4卡5卡| 超色免费av| 美女视频免费永久观看网站| 午夜老司机福利片| 久久久国产欧美日韩av| 狂野欧美激情性xxxx| 国产又色又爽无遮挡免费看| 国产高清国产精品国产三级| 老汉色av国产亚洲站长工具| 黄色毛片三级朝国网站| 国产av精品麻豆| 老司机深夜福利视频在线观看| 夜夜爽天天搞| 精品一区二区三区视频在线观看免费 | 女同久久另类99精品国产91| 欧美亚洲 丝袜 人妻 在线| 18禁观看日本| 色婷婷av一区二区三区视频| avwww免费| 免费观看精品视频网站| 久久中文字幕人妻熟女| 九色亚洲精品在线播放|