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

    二氧化碳/甲烷/氮?dú)舛旌衔镌谟行蚪榭滋疾牧螩MK-3中的吸附和分離

    2011-12-11 09:08:08張勤學(xué)曹達(dá)鵬
    物理化學(xué)學(xué)報(bào) 2011年9期
    關(guān)鍵詞:等溫線勢(shì)能混合物

    彭 璇 張勤學(xué) 成 璇 曹達(dá)鵬

    (1北京化工大學(xué)信息科學(xué)與技術(shù)學(xué)院,北京100029;2北京化工大學(xué)理學(xué)院,北京100029; 3北京化工大學(xué)有機(jī)無機(jī)復(fù)合材料國(guó)家重點(diǎn)實(shí)驗(yàn)室材料分子設(shè)計(jì)模擬研究室,北京100029)

    二氧化碳/甲烷/氮?dú)舛旌衔镌谟行蚪榭滋疾牧螩MK-3中的吸附和分離

    彭 璇1,*張勤學(xué)2,3成 璇3曹達(dá)鵬3

    (1北京化工大學(xué)信息科學(xué)與技術(shù)學(xué)院,北京100029;2北京化工大學(xué)理學(xué)院,北京100029;3北京化工大學(xué)有機(jī)無機(jī)復(fù)合材料國(guó)家重點(diǎn)實(shí)驗(yàn)室材料分子設(shè)計(jì)模擬研究室,北京100029)

    采用分子模擬與吸附理論研究了天然氣成分在有序介孔碳材料CMK-3上的吸附和分離.巨正則系綜蒙特卡羅(GCMC)模擬表明,CH4和CO2氣體的較優(yōu)存儲(chǔ)條件分別為208 K、4 MPa和298 K、6 MPa,其最大超額吸附量分別為10.07和14.85 mmol·g-1.基于雙位Langmuir-Freundlich(DSLF)模型,使用理想吸附溶液理論(IAST)預(yù)測(cè)了不同二元混合物在CMK-3中的分離行為,發(fā)現(xiàn)吸附選擇性SCO2/CH4與SCH4/N2比較接近,在298 K和4 MPa下約等于3,而N2-CO2體系中的CO2吸附選擇性較高,可達(dá)到7.5,說明CMK-3是一種適合吸附和分離天然氣組分的碳材料.

    巨正則系綜蒙特卡羅;吸附選擇性;理想吸附溶液理論;甲烷;二氧化碳;氮?dú)?CMK-3

    1 引言

    煤石油天然氣資源的枯竭和溫室效應(yīng)的加劇,迫使各國(guó)研究人員加速研發(fā)高效存儲(chǔ)甲烷氫氣能源和捕捉二氧化碳?xì)怏w的技術(shù).相比于液化和壓縮技術(shù),吸附技術(shù)在操作安全性、低成本和能源效率方面具有一定的優(yōu)勢(shì),因此制備新型多孔材料并應(yīng)用于氣體吸附和分離一直是國(guó)際性的學(xué)術(shù)前沿和研究熱點(diǎn).Jun和Ryoo1基于模板碳化技術(shù),2以蔗糖為碳源和介孔硅材料SBA-15為模板,成功地制備出了一種結(jié)構(gòu)呈六方形排列的新型介孔碳材料CMK-3.該材料的水熱穩(wěn)定性較好、比表面積較高、孔道結(jié)構(gòu)規(guī)則有序、孔徑分布狹窄和孔徑大小連續(xù)可調(diào),有望在吸附、分離、催化和光電磁等領(lǐng)域得到廣泛的應(yīng)用.3

    目前研究者們已對(duì)CMK-3材料上的吸附開展了一些實(shí)驗(yàn)和模擬工作,4-7如Zhou等5制備了CMK-3材料并使用體積法測(cè)量了N2、CH4和CO2等小分子氣體的吸附等溫線;Xia等6在CO2活化的CMK-3材料上研究了H2的吸附存儲(chǔ),發(fā)現(xiàn)活化后儲(chǔ)氫量有明顯改善,在77 K和105Pa下能達(dá)到2.27% (w);Monk等7采用分子動(dòng)力學(xué)方法模擬了CMK-3材料孔中[dmim+][Cl-]型離子液體的結(jié)構(gòu)和動(dòng)力學(xué)性質(zhì).值得一提的是,Zhou等5根據(jù)實(shí)驗(yàn)測(cè)量的純氣體吸附等溫線推測(cè)CMK-3可能是一種適合分離混合物中CO2的材料,但此觀點(diǎn)有待實(shí)驗(yàn)和理論上的支持.更重要的是,CMK-3材料上的微觀吸附機(jī)理還不是很清楚.為此,本文將采用計(jì)算手段去研究該材料上的氣體吸附及分離性質(zhì).

    研究吸附現(xiàn)象的計(jì)算方法主要有分子模擬8和吸附理論9兩種.許多研究者提出了不同的吸附理論模型,如理想吸附溶液理論(IAST),10二維狀態(tài)方程,9,11簡(jiǎn)化局部密度函數(shù)理論,12Ono-Kondo模型.13近年來基于雙位Langmuir-Freundlich模型(DSLF)的IAST在碳材料的吸附分離研究中得到廣泛應(yīng)用,14,15我們將采用該方法來預(yù)測(cè)常溫下CMK-3中的分離效果,為吸附工藝條件的設(shè)計(jì)與優(yōu)化提供理論上的依據(jù).

    本文的工作由以下部分組成:(1)對(duì)CMK-3材料結(jié)構(gòu)進(jìn)行了模型化,建立了有效表征流體與碳層間的相互作用的勢(shì)能模型;(2)利用巨正則蒙特卡羅(GCMC)方法8,16模擬優(yōu)化了CH4和CO2純氣體在CMK-3上的存儲(chǔ)條件;(3)利用IAST評(píng)價(jià)了CH4-N2、CO2-CH4及CO2-N2二元體系的吸附分離性能.

    圖1 CMK-3材料模型示意圖Fig.1 Schematic diagram of CMK-3 material

    2 CMK-3材料結(jié)構(gòu)的模型化

    Ryoo等4通過透射電子顯微鏡(TEM)和X射線衍射(XRD)實(shí)驗(yàn),以及密度泛函理論(DFT)分析,建立了CMK-3的結(jié)構(gòu)模型,其示意圖如圖1所示.本文將借鑒Ryoo等提出的結(jié)構(gòu)模型16進(jìn)行分子模擬研究.由圖1可知,一個(gè)基本的模擬晶胞由四根呈六方排列的碳柱構(gòu)成,其中a為單元晶胞的結(jié)構(gòu)參數(shù),d為模擬碳柱的直徑,孔間距(H)定義為H=(2a/ 31/2)-d.盡管CMK-3材料在結(jié)構(gòu)上長(zhǎng)程有序,但從原子水平上觀察,則近程無序.對(duì)于這種本質(zhì)上為無定形碳的材料,常規(guī)的理論方法難以描述原子水平的微觀結(jié)構(gòu).為對(duì)該結(jié)構(gòu)模型進(jìn)一步簡(jiǎn)化,我們假設(shè)碳柱由一組同軸心的、層間距為0.34 nm的單原子層構(gòu)成.這種處理方式在其他類CMK材料的吸附模擬17-19中取得了較好的效果.

    3 分子模擬

    3.1 理想化勢(shì)能模型

    任一流體分子在CMK-3材料孔中的總勢(shì)能ΦT可表示為兩部分之和,一部分是流體分子與其它流體分子之間的勢(shì)能Φff,另一部分是流體分子與碳層之間的勢(shì)能Φsf,即

    由于模擬中連續(xù)勢(shì)能總是要被截?cái)?因此計(jì)算流體分子之間的勢(shì)能時(shí)采用截?cái)嗥苿?shì)能8

    式中ΦLJ為12-6 Lennard-Jones(L-J)勢(shì)能,rc為截?cái)喟霃?rc=5σff,σff是流體分子LJ勢(shì)的尺寸參數(shù).流體分子與單個(gè)碳層之間的作用勢(shì)由Gordon模型20計(jì)算:

    表1 Lennard-Jones勢(shì)能模型參數(shù)Table 1 Parameters of Lennard-Jones potential models

    其中F[α,β;γ;χ]為超比幾何函數(shù),R表示單個(gè)碳層的半徑,r表示流體分子到碳層的距離,ρsurf為二維碳層密度,ρsurf=38.2 nm-2,εsf與σsf分別表示流體分子與碳原子之間的能量與尺寸相互作用參數(shù),由Lorentz-Berthelot混合規(guī)則計(jì)算.不同流體分子與碳原子的L-J勢(shì)能模型參數(shù)在表117,21,22中列出.

    3.2 有效勢(shì)能模型

    由于上文描述的是一個(gè)高度理想化的模型,該模型對(duì)流體與碳層間的相互作用的表征并不十分準(zhǔn)確,因此為了有效地描述流體與碳層間的相互作用,我們?cè)诟叨壤硐牖P偷幕A(chǔ)上,對(duì)計(jì)算流體與碳層交互作用的Lorentz-Berthelot混合規(guī)則引入可調(diào)參數(shù)ksf21

    通過吻合吸附等溫線的實(shí)驗(yàn)與模擬結(jié)果來確定一個(gè)有效的交互作用參數(shù)ksf,從而使得該模型能有效地表征流體分子與碳層間的相互作用.

    3.3 模擬細(xì)節(jié)

    采用化學(xué)勢(shì)μ、體積V和溫度T固定的GCMC方法研究流體分子在CMK-3材料中的吸附.模擬中,碳柱長(zhǎng)度取20 nm,同時(shí)在x、y、z三維方向上施加周期邊界條件.每個(gè)模擬條件產(chǎn)生2×107個(gè)構(gòu)型,其中最初的1×107個(gè)構(gòu)型將被丟掉,以確保模擬平衡和消除初始構(gòu)型的影響.而在計(jì)算熱力學(xué)量的系綜平均時(shí),只對(duì)后面1×107個(gè)構(gòu)型進(jìn)行統(tǒng)計(jì),其總勢(shì)能和平均數(shù)密度的統(tǒng)計(jì)偏差都在2%以內(nèi).通常實(shí)驗(yàn)中測(cè)量的吸附量為超額性質(zhì),即排除吸附相體積被主體相氣體所填充后的吸附量,而分子模擬計(jì)算則為絕對(duì)性質(zhì).為將分子模擬結(jié)果與實(shí)驗(yàn)結(jié)果比較,根據(jù)以下關(guān)系式將絕對(duì)吸附量Na轉(zhuǎn)換成超額吸附量Nex

    其中ρb為氣體主體密度,由Peng-Robinson狀態(tài)方程計(jì)算得到,V為吸附劑的自由體積.

    4 吸附理論

    4.1 雙位Langmuir-Freundlich模型

    純組分吸附的等溫線由DSLF模型14,15關(guān)聯(lián),回歸的模型參數(shù)用于IAST進(jìn)行混合物的吸附預(yù)測(cè). DSLF模型由如下等式描述

    其中f為主體氣體在吸附平衡時(shí)的主體氣相逸度, Ni、ki和ni為待回歸的模型參數(shù),分別表示i吸附位的最大吸附量、親和常數(shù)以及對(duì)簡(jiǎn)單Langmuir方程的偏離程度.

    4.2 理想吸附溶液理論

    當(dāng)達(dá)到吸附平衡時(shí),氣相與吸附相之間的i組分逸度相等,即

    p為氣相總壓,xi、yi分別為i組分在吸附相、氣相的摩爾分?jǐn)?shù),Φi為由Peng-Robinson狀態(tài)方程計(jì)算的i組分氣相逸度系數(shù),fi0為i組分的吸附相逸度,π為吸附表面的擴(kuò)張壓力.對(duì)于二元混合物,吸附平衡時(shí)兩組分的擴(kuò)張壓力相等,即

    式中Ni0為組分i在逸度fi時(shí)的吸附量.采用數(shù)值積分和方程求根的方法,由上式可解出組分i在吸附相中的摩爾分?jǐn)?shù)xi.組分i對(duì)j的吸附選擇性Sij定義為

    5 結(jié)果與討論

    5.1 純組分吸附的分子模擬

    5.1.1 確定交互作用參數(shù)kfw

    為與實(shí)驗(yàn)吸附等溫線比較,模擬中根據(jù)文獻(xiàn)5的實(shí)驗(yàn)值選取CMK-3材料的結(jié)構(gòu)參數(shù),a=9.8 nm和H=3.8 nm.圖2給出了常溫下N2、CH4和CO2超額吸附等溫線的實(shí)驗(yàn)與模擬結(jié)果.如圖2所示,不進(jìn)行交互作用參數(shù)調(diào)整時(shí)(ksf=1.0),模擬值普遍低于實(shí)驗(yàn)值,表明真實(shí)無定形材料的結(jié)構(gòu)復(fù)雜性難以用一種簡(jiǎn)化模型去精確地描述.另外也說明模擬中流體分子與碳柱間的相互作用偏小.同樣,Nguyen等23發(fā)現(xiàn)如果基于常見的碳-碳作用勢(shì)參數(shù),使用Lorentz-Berthelot混合規(guī)則計(jì)算的碳與氫氣之間作用也存在偏低的現(xiàn)象.因此,我們適當(dāng)?shù)厣险{(diào)ksf,發(fā)現(xiàn)當(dāng)ksf= 1.25時(shí)吸附曲線的實(shí)驗(yàn)值和模擬值均吻合一致,得到的ksf參數(shù)可用于GCMC預(yù)測(cè)中.

    圖2 298 K下超額吸附等溫線的實(shí)驗(yàn)與模擬比較Fig.2 Comparison of experimental and simulated excess adsorption isotherms at 298 K(a)N2;(b)CH4;(c)CO2

    5.1.2 純氣體吸附存儲(chǔ)條件的優(yōu)化

    吸附劑的存儲(chǔ)量主要取決于吸附過程的工藝條件.為對(duì)實(shí)驗(yàn)中的吸附條件進(jìn)行優(yōu)化設(shè)計(jì),我們模擬了不同溫度下CH4和CO2在CMK-3材料中的超額吸附等溫線.如圖3所示,當(dāng)壓力較低時(shí),孔內(nèi)吸附流體密度大于主體相流體,提高壓力可使超額吸附量增加,壓力過高時(shí)則相反.因此,每條吸附等溫線均存在一個(gè)最優(yōu)的壓力使得對(duì)應(yīng)的超額吸附量最大.此外,從圖3還可看出,低壓下溫度的升高會(huì)引起超額吸附量的降低,但是高壓高溫下的超額吸附量反而較大.這可根據(jù)等式(5)來解釋,對(duì)于前者,低壓下的主體相密度很低,故對(duì)超額吸附的影響忽略不計(jì),而溫度的升高使流體分子的動(dòng)能增加,從而削弱了碳材料對(duì)流體分子的吸附能力,導(dǎo)致絕對(duì)吸附量的減少;對(duì)于后者,高溫引起的主體相密度降低對(duì)超額吸附的影響占主導(dǎo)作用,明顯強(qiáng)于高壓引起的主體相密度增加和高溫引起的絕對(duì)吸附量降低這兩種因素.因此,這種主體密度的凈減少使得超額吸附量隨著溫度增加而提高.圖4給出了CH4和CO2在不同溫度下的超額吸附量極值隨著最優(yōu)操作壓力的變化.由圖4可知,最大氣體存儲(chǔ)量隨著溫度的升高快速下降;同時(shí),滿足該極值目標(biāo)需相應(yīng)地提升壓力.對(duì)于CH4和CO2的存儲(chǔ),較優(yōu)的條件為208 K、4 MPa和298 K、6 MPa,其最大吸附量分別為10.07和14.85 mmol·g-1.

    為了從微觀水平上理解氣體在CMK-3中存儲(chǔ)的機(jī)理,我們?cè)趫D5和6中分別給出了常溫條件下CH4和CO2達(dá)到吸附平衡時(shí)的分子構(gòu)型和對(duì)應(yīng)的孔內(nèi)流體局部密度分布.由圖5可知,氣體分子在CMK-3材料的表面出現(xiàn)大量聚集的現(xiàn)象,但基本上仍以單層吸附為主,同時(shí)隨著壓力的升高,吸附量明顯增加;相同壓力條件下,由于CO2具有較大的勢(shì)能作用參數(shù),CMK-3對(duì)CO2的吸附能力要強(qiáng)于CH4.圖6中,兩種物質(zhì)的局部密度在r*=10時(shí)均最高,這對(duì)應(yīng)著材料表面上的吸附單層,而在r*=12,15時(shí)均出現(xiàn)兩個(gè)較弱的峰,其中前者為發(fā)生在吸附單層上面的第二層吸附,后者是發(fā)生在兩相鄰碳柱距離最短處因勢(shì)能疊加引起的吸附層.此外,隨著壓力的提高,峰的強(qiáng)度變得更加明顯.這些觀察與圖5中的流體分子構(gòu)型是相互一致的.

    圖3 GCMC模擬的甲烷與二氧化碳在不同溫度下的超額吸附等溫線Fig.3 Excess adsorption isotherms of methane and carbon dioxide from GCMC simulation at different temperatures(a)CH4;(b)CO2

    圖4 不同溫度下甲烷與二氧化碳的最大超額吸附量與對(duì)應(yīng)的最優(yōu)壓力Fig.4 Maximum adsorption amounts and the corresponding optimum pressures at different temperatures(a)CH4;(b)CO2

    5.2 二元混合物的吸附分離

    5.2.1 IAST方法的驗(yàn)證

    為節(jié)省GCMC方法的計(jì)算耗時(shí),本文選擇更簡(jiǎn)單的IAST預(yù)測(cè)方法來考察二元體系的分離性質(zhì),即以雙位Langmuir-Freundlich模型關(guān)聯(lián)純氣體吸附的實(shí)驗(yàn)數(shù)據(jù),再使用IAST計(jì)算混合物的吸附性質(zhì).我們先使用Sudibandriyo等24在活性炭上測(cè)量的純氣體和混合物吸附數(shù)據(jù)驗(yàn)證IAST方法的可靠性.從圖7可以看出,即使壓力高至13 MPa,Langmuir模型擬合和IAST的預(yù)測(cè)均與實(shí)驗(yàn)吻合,證明這種研究方法是可以接受的.表2給出了擬合CMK-3材料上純氣體的實(shí)驗(yàn)數(shù)據(jù)得到的雙位Langmuir-Freundlich模型參數(shù).由表2看出,對(duì)三種純氣體的吸附,擬合的吸附量與實(shí)驗(yàn)之間的平均相對(duì)偏差(ARD)均小于1%,表明回歸的Langmuir模型參數(shù)可用于混合物吸附的IAST預(yù)測(cè).

    5.2.2 等摩爾比下的吸附選擇性

    圖5 298 K下甲烷與二氧化碳在CMK-3中不同壓力下的分子構(gòu)型圖Fig.5 Molecular configurations of methane and carbon dioxide in CMK-3 at 298 K and different pressures(1,4)1 MPa;(2,5)6 MPa;(3,6)12 MPa

    圖6 CMK-3材料中甲烷與二氧化碳在不同壓力下的局部密度分布圖Fig.6 Local density profiles of methane and carbon dioxide at different pressures in CMK-3(a)CH4;(b)CO2;ρ*:local density,r*:reduced distance from axis

    圖7 318.2 K下活性炭材料上絕對(duì)吸附等溫線計(jì)算值和實(shí)驗(yàn)值比較Fig.7 Comparison of the calculated absolute adsorption amount with experimental data on dry activated carbon at 318.2 K(a)pure gases;(b)CH4-CO2mixture with CH4feed composition of 40.1%

    表2 CMK-3材料中純組分吸附的雙位Langmuir-Freundlich模型參數(shù)Table 2 Parameters of dual-site Langmuir-Freundlich model for pure adsorption in CMK-3

    圖8298 K下IAST預(yù)測(cè)CMK-3材料上等摩爾混合物的選擇性(S)Fig.8 IAST predictions for adsorption selectivity(S)of equimolar mixtures in CMK-3 at 298 K

    圖8為CH4-CO2、N2-CH4及N2-CO2二元混合物在常溫298 K和氣相等摩爾比下的吸附選擇性.從圖8可以看出,低壓(<0.2 MPa)下的吸附選擇性較高,這說明CMK-3材料中存在部分微孔的吸附作用.壓力稍微增加時(shí),微孔容積已被流體分子填充,引起吸附選擇性的陡然下降.當(dāng)壓力高于0.5 MPa時(shí),材料中介孔中的固體孔壁對(duì)流體的多層吸附能力較弱,而單層表面吸附對(duì)分離的貢獻(xiàn)起主導(dǎo)作用,這使得吸附選擇性的變化逐漸平緩.此外,對(duì)于CH4-CO2體系,其CO2吸附選擇性與N2-CH4中的CH4選擇性相差不大,在較高壓力約為3,而CMK-3材料對(duì)N2-CO2中的CO2具有較強(qiáng)的優(yōu)先吸附能力,其選擇性穩(wěn)定后保持在7.5左右.這可以根據(jù)圖9中氣體混合物中組分的絕對(duì)吸附等溫線解釋.在圖9(a)與(b)中,CH4、CO2之間和N2、CH4之間的吸附等溫線差距隨著壓力的變化總是很小,故兩體系分離的選擇性比較接近;而圖9(c)中N2、CO2之間的吸附量差距要明顯大于另外兩種體系,可見CMK-3材料尤其適合分離N2-CO2體系中的CO2.盡管低壓下的氣體吸附選擇性較高,但低壓也不利于組分的存儲(chǔ).綜合權(quán)衡,建議在4 MPa下進(jìn)行吸附操作,以同時(shí)滿足分離和存儲(chǔ)的性能需求.

    圖9298 K下IAST預(yù)測(cè)CMK-3中等摩爾混合物的絕對(duì)吸附等溫線Fig.9 IAST predictions for absolute adsorption isotherms of equimolar mixtures in CMK-3 at 298 K

    圖10 298 K下IAST預(yù)測(cè)CMK-3中二元混合物不同摩爾比下的吸附選擇性Fig.10 IAST predictions for adsorption selectivity of binary mixtures with different molar ratios in CMK-3 at 298 K as a function of bulk pressure

    5.2.3 不同摩爾比下的吸附選擇性

    為考察主體組成對(duì)分離性能的影響,圖10給出了常溫下不同氣相摩爾比的二元混合物吸附選擇性的變化.

    由該圖可知,N2-CH4體系中不同氣相摩爾比的吸附曲線幾乎重合,說明吸附選擇性受主體相組成的影響不大.然而對(duì)于含CO2的體系,CO2選擇性在壓力大于2.0 MPa時(shí)隨著壓力和氣相中CO2濃度而增加,并且壓力越高現(xiàn)象越顯著.這可根據(jù)CO2分子具有較強(qiáng)的四極矩特性解釋.由于孔內(nèi)吸附相CO2分子的密度高于主體相,在四極矩的作用下孔內(nèi)CO2流體會(huì)吸引更多的主體相CO2分子進(jìn)入孔內(nèi),導(dǎo)致孔內(nèi)的吸附量和選擇性提高.

    6 結(jié)論

    通過巨正則系綜蒙特卡羅模擬提出了一種有效表征CMK-3材料結(jié)構(gòu)的理論模型,并對(duì)CH4和CO2純氣體存儲(chǔ)進(jìn)行了吸附條件的優(yōu)化.同時(shí)利用基于雙位Langmuir-Freundlich吸附模型的理想吸附溶液理論考察了N2/CH4/CO2二元混合物在CMK-3材料上的分離性質(zhì).計(jì)算表明,CMK-3材料不僅適合于CH4和CO2氣體的存儲(chǔ),也是一種高效分離天然氣成分的吸附劑.

    (1) Jun,S.;Ryoo,R.J.Am.Chem.Soc.2000,122,10712.

    (2) Kyotani,T.;Sonobe,N.;Tomita,A.Nature 1988,331,331.

    (3) Lee,J.;Han,S.;Hyeon,T.J.Mater.Chem.2004,14,347.

    (4) Joo,S.H.;Ryoo,R.;Kruk,M.J.Phys.Chem.B.2002,106, 4640.

    (5)Zhou,L.;Liu,X.W.;Wang,N.;Wang,Z.;Zhou,Y.P.Chem. Phys.Lett.2005,413,6.

    (6)Xia,K.;Gao,Q.;Wu,C.;Song,S.;Ruan,M.Carbon 2007,45, 1989.

    (7) Monk,J.;Singh,R.;Hung,F.J.Phys.Chem.C 2011,115,3034.

    (8) Allen,M.P.;Tildesley,D.J.Computer Simulation of Liquids; Clarendon Press:Oxford,1987;pp 110-139.

    (9)Peng,X.;Wang,W.C.;Xue,R.S.;Shen,Z.M.AIChE J.2006, 52,994.

    (10) Myers,A.L.;Prausnitz,J.M.AIChE J.1965,11,121.

    (11) Zhou,C.;Hall,F.;Gasem,K.;Robinson,J.Ind.Eng.Chem. Res.1994,33,1280.

    (12) Rangarajan,B.;Lira,C.T.;Subramanian,R.AIChE J.1995,41, 838.

    (13) Ottiger,S.;Pini,R.;Storti,G.;Mazzotti,M.Langmuir 2008,24, 9531.

    (14)Peng,X.;Cao,D.P.;Wang,W.C.Ind.Eng.Chem.Res.2010, 49,8787.

    (15)Peng,X.;Cao,D.P.;Wang,W.C.Chem.Eng.Sci.2011,66, 2266.

    (16)Peng,X.;Zhou,J.;Wang,W.C.;Cao,D.P.Carbon 2010,48, 3760.

    (17) Peng,X.;Cao,D.P.;Wang,W.C.Langmuir 2009,25,10863.

    (18)Peng,X.;Cao,D.P.;Wang,W.C.J.Phys.Chem.C 2008,112, 13024.

    (19) Peng,X.;Cao,D.P.;Zhao,J.S.Sep.Purif.Technol.2009,68, 50.

    (20) Gordon,P.A.;Saeger,R.B.Ind.Eng.Chem.Res.1999,38, 4647.

    (21) Cao,D.P.;Wang,W.C.;Duan,X.J.Colloid Interface Sci. 2002,254,1.

    (22) Peng,X.;Zhao,J.S.;Cao,D.P.J.Colloid Interface Sci.2007, 310,391.

    (23) Nguyen,T.X.;Bae,J.;Wang,Y.;Bhatia,S.K.Langmuir 2009, 25,431.

    (24) Sudibandriyo,M.;Pan,Z.J.;Fitzgerald,J.E.;Robinson,R.L., Jr.;Gasem,K.A.M.Langmuir 2003,19,5323.

    April 20,2011;Revised:June 8,2011;Published on Web:July 12,2011.

    Adsorption and Separation of CO2/CH4/N2Binary Mixtures in an Ordered Mesoporous Carbon Material CMK-3

    PENG Xuan1,*ZHANG Qin-Xue2,3CHENG Xuan3CAO Da-Peng3
    (1College of Information Science and Technology,Beijing University of Chemical Technology,Beijing 100029,P.R.China;2College of Science,Beijing University of Chemical Technology,Beijing 100029,P.R.China;3Division of Molecular and Materials Simulation,State Key Laboratory of Organic-Inorganic Composites, Beijing University of Chemical Technology,Beijing 100029,P.R.China)

    The adsorption and separation of natural gas in the ordered mesoporous carbon material CMK-3 was investigated by molecular simulation and adsorption theory.Grand canonical ensemble Monte Carlo (GCMC)simulations show that a maximum excess uptake of 10.07 and 14.85 mmol·g-1is obtained at the optimum temperature and pressure of 208 K,4 MPa for CH4and 298 K,6 MPa for CO2adsorption, respectively.Based on the dual-site Langmuir-Freundlich(DSLF)model,ideal adsorption solution theory (IAST)was used to predict the adsorption and separation of binary mixtures.The adsorption selectivity of SCO2/CH4is approximately the same as that ofSCH4/N2,with a value of about 3 at 298 K and 4 MPa while the highest CO2selectivity of 7.5 was found in the N2-CO2system under the same conditions.This indicates that theCMK-3 material isapromising candidate for natural gas separation.

    Grand canonical ensemble Monte Carlo;Adsorption selectivity;Ideal adsorption solution theory;Methane;Carbon dioxide;Nitrogen;CMK-3

    ?Corresponding author.Email:pengxuan@mail.buct.edu.cn;Tel:+86-10-64443254.

    The project was supported by the National Natural Science Foundation of China(20806003).

    國(guó)家自然科學(xué)基金(20806003)資助項(xiàng)目

    O647;O641

    猜你喜歡
    等溫線勢(shì)能混合物
    “動(dòng)能和勢(shì)能”知識(shí)鞏固
    作 品:景觀設(shè)計(jì)
    ——《勢(shì)能》
    文化縱橫(2022年3期)2022-09-07 11:43:18
    多組分纖維混合物定量分析通用計(jì)算模型研制
    “動(dòng)能和勢(shì)能”知識(shí)鞏固
    正丁醇和松節(jié)油混合物對(duì)組織脫水不良的補(bǔ)救應(yīng)用
    “動(dòng)能和勢(shì)能”隨堂練
    低溫吸附劑的低溫低壓吸附等溫線分布研究
    混合物按照歐盟CLP進(jìn)行分類標(biāo)簽
    如何在新課改背景下突破等溫線判讀中的難點(diǎn)
    萃取精餾分離甲苯-正庚烷混合物的模擬研究
    亚洲人成网站在线观看播放| 亚洲精华国产精华液的使用体验 | 久久久久久九九精品二区国产| 久久久成人免费电影| 听说在线观看完整版免费高清| 久久国内精品自在自线图片| 十八禁国产超污无遮挡网站| 欧美最新免费一区二区三区| 深夜a级毛片| 成人特级av手机在线观看| 亚洲婷婷狠狠爱综合网| 欧洲精品卡2卡3卡4卡5卡区| 婷婷精品国产亚洲av| 欧美日本亚洲视频在线播放| 国产av不卡久久| 插阴视频在线观看视频| 国产精品一区二区免费欧美| 亚洲中文字幕一区二区三区有码在线看| 亚洲成人中文字幕在线播放| 女人十人毛片免费观看3o分钟| 欧美色欧美亚洲另类二区| 精品久久久噜噜| 国产精品1区2区在线观看.| www日本黄色视频网| 亚洲性夜色夜夜综合| 久99久视频精品免费| 最近最新中文字幕大全电影3| 色综合色国产| 日韩制服骚丝袜av| 日韩av在线大香蕉| 国内久久婷婷六月综合欲色啪| 亚洲最大成人中文| 国产不卡一卡二| 欧美性感艳星| 99久国产av精品| 中文字幕av在线有码专区| 免费在线观看影片大全网站| 亚洲自偷自拍三级| 国产真实乱freesex| 99热这里只有是精品50| av在线亚洲专区| 一区二区三区高清视频在线| 欧美一区二区国产精品久久精品| 热99re8久久精品国产| 中文字幕久久专区| 欧美最新免费一区二区三区| 免费电影在线观看免费观看| 99久久精品国产国产毛片| 精品午夜福利视频在线观看一区| 伦理电影大哥的女人| 搡老岳熟女国产| 中国国产av一级| 美女 人体艺术 gogo| 精品国内亚洲2022精品成人| 大型黄色视频在线免费观看| 18禁在线播放成人免费| 中文字幕免费在线视频6| 特级一级黄色大片| 国内揄拍国产精品人妻在线| 欧美成人a在线观看| 晚上一个人看的免费电影| 亚洲av电影不卡..在线观看| 亚洲最大成人中文| 成人毛片a级毛片在线播放| 波野结衣二区三区在线| 国产精品乱码一区二三区的特点| 国产色婷婷99| 欧美成人免费av一区二区三区| 我要搜黄色片| 亚洲国产色片| 国产精品一区www在线观看| 日本在线视频免费播放| 啦啦啦韩国在线观看视频| 免费看光身美女| 中文资源天堂在线| 国产美女午夜福利| 色av中文字幕| 99热这里只有精品一区| 精品人妻一区二区三区麻豆 | 日韩,欧美,国产一区二区三区 | 国产精品久久久久久av不卡| 97热精品久久久久久| 亚洲精品国产av成人精品 | 成人毛片a级毛片在线播放| 日本五十路高清| 深夜a级毛片| 3wmmmm亚洲av在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产高清视频在线观看网站| 亚洲无线在线观看| 国产高清有码在线观看视频| 一本一本综合久久| 看片在线看免费视频| 国内少妇人妻偷人精品xxx网站| 日韩欧美三级三区| 男插女下体视频免费在线播放| 又黄又爽又刺激的免费视频.| 亚洲电影在线观看av| 色吧在线观看| 欧美xxxx性猛交bbbb| 午夜免费男女啪啪视频观看 | 亚洲人成网站在线播| 婷婷亚洲欧美| АⅤ资源中文在线天堂| 麻豆av噜噜一区二区三区| 国产精品99久久久久久久久| 一a级毛片在线观看| 一个人观看的视频www高清免费观看| 蜜桃久久精品国产亚洲av| h日本视频在线播放| 91麻豆精品激情在线观看国产| 亚洲精品久久国产高清桃花| 男女那种视频在线观看| 精品人妻偷拍中文字幕| 国产精品女同一区二区软件| 极品教师在线视频| 亚洲欧美精品自产自拍| 欧美性猛交黑人性爽| 在线看三级毛片| 国产亚洲精品av在线| 大型黄色视频在线免费观看| 日韩av在线大香蕉| 天堂网av新在线| 一区二区三区高清视频在线| 深爱激情五月婷婷| 久久久a久久爽久久v久久| 国产成人精品久久久久久| 少妇的逼好多水| 51国产日韩欧美| 日日摸夜夜添夜夜添小说| 嫩草影视91久久| 精品乱码久久久久久99久播| 丝袜美腿在线中文| 又黄又爽又刺激的免费视频.| 国内少妇人妻偷人精品xxx网站| 欧美性猛交╳xxx乱大交人| 国产黄色视频一区二区在线观看 | 国产精品一二三区在线看| 精品少妇黑人巨大在线播放 | 午夜福利18| 中文字幕熟女人妻在线| 日本三级黄在线观看| 在线天堂最新版资源| 国国产精品蜜臀av免费| 亚洲av电影不卡..在线观看| 日本免费一区二区三区高清不卡| 国产高清有码在线观看视频| av专区在线播放| 国产一区二区三区av在线 | 久久久久久久久久久丰满| 在线天堂最新版资源| 国产一区二区亚洲精品在线观看| 晚上一个人看的免费电影| 99视频精品全部免费 在线| 精品久久国产蜜桃| 国产91av在线免费观看| 亚洲成人久久性| 亚洲av一区综合| 老司机福利观看| 国产精品一区二区三区四区免费观看 | 美女高潮的动态| 特级一级黄色大片| 国产亚洲91精品色在线| 久久草成人影院| 黄色日韩在线| 狂野欧美激情性xxxx在线观看| 男人狂女人下面高潮的视频| 国产午夜福利久久久久久| 亚洲最大成人av| 99热这里只有精品一区| 99久久久亚洲精品蜜臀av| 亚洲最大成人av| 高清毛片免费看| 国产午夜福利久久久久久| or卡值多少钱| 乱码一卡2卡4卡精品| 国内少妇人妻偷人精品xxx网站| 精品免费久久久久久久清纯| 人妻丰满熟妇av一区二区三区| 99国产极品粉嫩在线观看| 搡老妇女老女人老熟妇| 欧美最黄视频在线播放免费| 精品国产三级普通话版| 一本一本综合久久| 三级经典国产精品| 18禁在线播放成人免费| 国产成人a区在线观看| 免费观看在线日韩| 丝袜美腿在线中文| 尾随美女入室| 精品乱码久久久久久99久播| 美女高潮的动态| 在线天堂最新版资源| 黄色视频,在线免费观看| 在线播放国产精品三级| 免费看a级黄色片| 天天躁日日操中文字幕| 97超碰精品成人国产| 午夜a级毛片| 成人av在线播放网站| 搞女人的毛片| 成年女人永久免费观看视频| 尤物成人国产欧美一区二区三区| 欧美日韩在线观看h| 97热精品久久久久久| 18+在线观看网站| 精品国产三级普通话版| 国产日本99.免费观看| 亚洲av二区三区四区| 色综合站精品国产| 日韩欧美国产在线观看| 黄片wwwwww| 99在线视频只有这里精品首页| 国产伦在线观看视频一区| 亚洲av.av天堂| 午夜激情福利司机影院| 中文亚洲av片在线观看爽| 噜噜噜噜噜久久久久久91| 欧美国产日韩亚洲一区| 色在线成人网| 成人av在线播放网站| 色综合色国产| 18禁在线播放成人免费| 啦啦啦韩国在线观看视频| 少妇丰满av| 亚洲国产精品成人久久小说 | 搡女人真爽免费视频火全软件 | 插逼视频在线观看| 亚洲七黄色美女视频| 可以在线观看毛片的网站| 久久精品久久久久久噜噜老黄 | 高清毛片免费观看视频网站| 人妻制服诱惑在线中文字幕| 综合色av麻豆| 亚洲三级黄色毛片| 亚洲av.av天堂| 啦啦啦啦在线视频资源| 不卡一级毛片| 亚洲婷婷狠狠爱综合网| 日韩,欧美,国产一区二区三区 | 我的女老师完整版在线观看| 如何舔出高潮| 亚洲真实伦在线观看| 最近手机中文字幕大全| 日韩精品中文字幕看吧| 亚洲av一区综合| 联通29元200g的流量卡| 婷婷色综合大香蕉| 麻豆av噜噜一区二区三区| 国内精品美女久久久久久| 欧美激情在线99| 久久久欧美国产精品| 色综合亚洲欧美另类图片| 干丝袜人妻中文字幕| 全区人妻精品视频| 日韩av在线大香蕉| 不卡视频在线观看欧美| 在线免费观看的www视频| 18禁裸乳无遮挡免费网站照片| 国产伦一二天堂av在线观看| 亚洲无线在线观看| 国产av麻豆久久久久久久| 热99在线观看视频| 精品日产1卡2卡| 男女啪啪激烈高潮av片| 精品欧美国产一区二区三| 婷婷亚洲欧美| av在线播放精品| 日韩欧美在线乱码| 国产av不卡久久| 又爽又黄a免费视频| 91av网一区二区| 精品人妻熟女av久视频| 高清午夜精品一区二区三区 | 人妻久久中文字幕网| 在线观看av片永久免费下载| 亚洲国产精品成人久久小说 | 久久久久久久久久黄片| 又爽又黄a免费视频| 精华霜和精华液先用哪个| 午夜激情福利司机影院| 大型黄色视频在线免费观看| 成人av在线播放网站| 国产熟女欧美一区二区| 午夜老司机福利剧场| 国产免费男女视频| 美女xxoo啪啪120秒动态图| 中文字幕久久专区| 国产男人的电影天堂91| 亚洲图色成人| 久久精品久久久久久噜噜老黄 | 国产黄色视频一区二区在线观看 | 免费av不卡在线播放| 日韩成人av中文字幕在线观看 | 女人十人毛片免费观看3o分钟| 18+在线观看网站| 国产国拍精品亚洲av在线观看| 国产精品永久免费网站| 男女下面进入的视频免费午夜| 国产69精品久久久久777片| 12—13女人毛片做爰片一| 一级毛片我不卡| 亚洲av五月六月丁香网| 一个人免费在线观看电影| 看片在线看免费视频| 国产一区亚洲一区在线观看| 亚洲天堂国产精品一区在线| 国产私拍福利视频在线观看| 欧美另类亚洲清纯唯美| 久久精品国产自在天天线| 国产69精品久久久久777片| 日韩精品青青久久久久久| 亚洲最大成人手机在线| 寂寞人妻少妇视频99o| 真人做人爱边吃奶动态| 国产国拍精品亚洲av在线观看| 国产aⅴ精品一区二区三区波| 精品久久久久久久久亚洲| 久久99热6这里只有精品| 日韩av不卡免费在线播放| 午夜福利视频1000在线观看| 高清午夜精品一区二区三区 | 国产 一区 欧美 日韩| 搡女人真爽免费视频火全软件 | 亚洲精品一区av在线观看| 精品一区二区三区视频在线| 成年女人看的毛片在线观看| 男插女下体视频免费在线播放| av福利片在线观看| 天天躁日日操中文字幕| 在线观看66精品国产| 狠狠狠狠99中文字幕| 国产午夜精品论理片| 国产精品日韩av在线免费观看| 亚洲三级黄色毛片| 国产精品三级大全| 成人欧美大片| 亚洲av成人av| 国产精品永久免费网站| 麻豆久久精品国产亚洲av| 狂野欧美激情性xxxx在线观看| 成人美女网站在线观看视频| 免费不卡的大黄色大毛片视频在线观看 | 国产精品1区2区在线观看.| 91在线精品国自产拍蜜月| 成人永久免费在线观看视频| 一夜夜www| 国产精品99久久久久久久久| 国产高清有码在线观看视频| 国产白丝娇喘喷水9色精品| 日本a在线网址| 国产午夜精品论理片| 99九九线精品视频在线观看视频| 日韩成人伦理影院| 偷拍熟女少妇极品色| 尾随美女入室| 一夜夜www| 欧美绝顶高潮抽搐喷水| 丝袜美腿在线中文| 在线观看66精品国产| 精品久久久久久久末码| 在线a可以看的网站| 欧美成人a在线观看| 一个人看的www免费观看视频| 女人十人毛片免费观看3o分钟| 人妻制服诱惑在线中文字幕| 国产高清有码在线观看视频| 黄色日韩在线| 国产精华一区二区三区| 免费在线观看影片大全网站| 最近2019中文字幕mv第一页| av在线观看视频网站免费| 色视频www国产| av中文乱码字幕在线| 国产视频内射| 精华霜和精华液先用哪个| 欧美成人精品欧美一级黄| 18禁黄网站禁片免费观看直播| 亚洲一区高清亚洲精品| 欧美bdsm另类| 国产久久久一区二区三区| 国产不卡一卡二| 蜜桃久久精品国产亚洲av| 草草在线视频免费看| 亚洲美女搞黄在线观看 | 日韩一区二区视频免费看| 亚洲精品粉嫩美女一区| 中文字幕av在线有码专区| 中文字幕熟女人妻在线| 一个人看视频在线观看www免费| 国产精品女同一区二区软件| 国产精品一区二区三区四区免费观看 | 久久久久久久久久久丰满| 日韩欧美一区二区三区在线观看| 日本三级黄在线观看| 久久久久久久久大av| 亚洲第一电影网av| 日韩欧美三级三区| 国产精品野战在线观看| 99九九线精品视频在线观看视频| 黄色配什么色好看| 亚洲高清免费不卡视频| 欧美日韩一区二区视频在线观看视频在线 | av在线天堂中文字幕| 97超碰精品成人国产| 免费大片18禁| 国产精品国产高清国产av| 欧美成人a在线观看| 插逼视频在线观看| 简卡轻食公司| 亚洲成人精品中文字幕电影| 精品日产1卡2卡| 草草在线视频免费看| 亚洲天堂国产精品一区在线| 国产黄a三级三级三级人| av专区在线播放| 激情 狠狠 欧美| 国产精品不卡视频一区二区| 天美传媒精品一区二区| 精品福利观看| 亚洲av免费高清在线观看| 成人毛片a级毛片在线播放| 一级黄片播放器| 国产真实伦视频高清在线观看| 老司机福利观看| 精品一区二区免费观看| 国产成人a区在线观看| 亚洲精品粉嫩美女一区| 成人永久免费在线观看视频| 性插视频无遮挡在线免费观看| 熟女电影av网| 狂野欧美白嫩少妇大欣赏| 午夜福利18| 女的被弄到高潮叫床怎么办| 在线观看一区二区三区| 国产精品一区二区免费欧美| 五月玫瑰六月丁香| 久久久久精品国产欧美久久久| 国产一区亚洲一区在线观看| 国产欧美日韩精品一区二区| 中国国产av一级| 精品少妇黑人巨大在线播放 | 俄罗斯特黄特色一大片| 日韩,欧美,国产一区二区三区 | 自拍偷自拍亚洲精品老妇| 成年女人毛片免费观看观看9| 欧美不卡视频在线免费观看| 欧美在线一区亚洲| 中文亚洲av片在线观看爽| 亚洲国产精品成人久久小说 | 日本免费a在线| 婷婷亚洲欧美| 99久国产av精品国产电影| 亚洲激情五月婷婷啪啪| 老司机福利观看| av免费在线看不卡| 日韩一区二区视频免费看| 超碰av人人做人人爽久久| 欧美激情久久久久久爽电影| videossex国产| 一级黄色大片毛片| 麻豆成人午夜福利视频| 深夜a级毛片| 中文字幕免费在线视频6| 菩萨蛮人人尽说江南好唐韦庄 | 给我免费播放毛片高清在线观看| 欧美又色又爽又黄视频| 欧洲精品卡2卡3卡4卡5卡区| 国产国拍精品亚洲av在线观看| 人人妻人人看人人澡| 日本免费a在线| 欧美日韩综合久久久久久| 午夜免费激情av| 成人毛片a级毛片在线播放| 99热这里只有是精品在线观看| 午夜a级毛片| 精品久久久噜噜| 听说在线观看完整版免费高清| 白带黄色成豆腐渣| 国产淫片久久久久久久久| 亚洲乱码一区二区免费版| 如何舔出高潮| 亚洲国产欧洲综合997久久,| 欧美人与善性xxx| 国模一区二区三区四区视频| 欧美+亚洲+日韩+国产| 人妻久久中文字幕网| 久久久久久久久久黄片| 亚洲成人av在线免费| 在现免费观看毛片| 精品人妻视频免费看| 久久精品国产鲁丝片午夜精品| 国产男人的电影天堂91| 亚洲,欧美,日韩| 性插视频无遮挡在线免费观看| 成人av一区二区三区在线看| 极品教师在线视频| 国产男人的电影天堂91| 日韩国内少妇激情av| 午夜免费激情av| av天堂在线播放| 亚洲精品乱码久久久v下载方式| 精品午夜福利在线看| 97超碰精品成人国产| 嫩草影院入口| 成人av一区二区三区在线看| 老司机福利观看| 欧美精品国产亚洲| 亚洲欧美日韩东京热| 久久国内精品自在自线图片| 亚洲自拍偷在线| 亚洲国产精品国产精品| 日韩强制内射视频| 欧美激情久久久久久爽电影| 亚洲欧美日韩高清专用| 亚洲欧美日韩无卡精品| 精品一区二区三区av网在线观看| АⅤ资源中文在线天堂| 国产在线精品亚洲第一网站| 日本精品一区二区三区蜜桃| 欧美潮喷喷水| 亚洲国产精品sss在线观看| 校园人妻丝袜中文字幕| 99久国产av精品| 亚洲av.av天堂| 美女被艹到高潮喷水动态| 99热网站在线观看| 欧美成人免费av一区二区三区| www日本黄色视频网| 亚洲国产精品久久男人天堂| 国产精品一及| 午夜精品一区二区三区免费看| 国产精品一及| 联通29元200g的流量卡| a级毛片a级免费在线| 国产精品亚洲美女久久久| 免费人成视频x8x8入口观看| av在线亚洲专区| 欧美日韩乱码在线| 成人高潮视频无遮挡免费网站| 欧美成人一区二区免费高清观看| 白带黄色成豆腐渣| 99riav亚洲国产免费| 色哟哟哟哟哟哟| 高清毛片免费观看视频网站| 日本成人三级电影网站| 97热精品久久久久久| 亚洲国产高清在线一区二区三| 男插女下体视频免费在线播放| 国产极品精品免费视频能看的| 精品一区二区三区视频在线| 91午夜精品亚洲一区二区三区| 国产视频内射| 亚洲真实伦在线观看| 91精品国产九色| 日韩中字成人| 国产精品久久电影中文字幕| 日韩欧美一区二区三区在线观看| 99热这里只有精品一区| 久久精品国产自在天天线| 亚洲国产精品久久男人天堂| 国产精品,欧美在线| 最近最新中文字幕大全电影3| 亚洲熟妇熟女久久| 午夜视频国产福利| or卡值多少钱| 亚洲,欧美,日韩| 国产精品久久久久久av不卡| 夜夜看夜夜爽夜夜摸| 最近手机中文字幕大全| 身体一侧抽搐| 偷拍熟女少妇极品色| 欧美极品一区二区三区四区| 精品一区二区三区视频在线观看免费| 免费搜索国产男女视频| 国产精品久久久久久久电影| 国产精品美女特级片免费视频播放器| 欧美日本亚洲视频在线播放| 特大巨黑吊av在线直播| 免费av毛片视频| 一级毛片aaaaaa免费看小| 国产一区二区亚洲精品在线观看| 日韩欧美在线乱码| 日本爱情动作片www.在线观看 | 99久久中文字幕三级久久日本| 国产色爽女视频免费观看| 亚洲精品粉嫩美女一区| av在线亚洲专区| 亚洲精品影视一区二区三区av| 国产午夜福利久久久久久| 美女 人体艺术 gogo| 精华霜和精华液先用哪个| 精品人妻一区二区三区麻豆 | 99久国产av精品国产电影| 男人舔奶头视频| 极品教师在线视频| av免费在线看不卡| 国产伦精品一区二区三区四那| 亚洲av一区综合| 日本成人三级电影网站| 久久午夜亚洲精品久久| 日产精品乱码卡一卡2卡三| 亚洲av电影不卡..在线观看| 精品乱码久久久久久99久播| 久久久久久久久中文| 男女之事视频高清在线观看| 在线a可以看的网站| 国产成人aa在线观看| 国产 一区精品| 亚洲丝袜综合中文字幕| 性色avwww在线观看| 午夜福利18| 久久久久久久亚洲中文字幕| 夜夜爽天天搞| 亚洲真实伦在线观看|