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

    含不飽和金屬位點(diǎn)的金屬有機(jī)骨架材料對(duì)煙道氣中SO2的吸附與分離行為的分子模擬研究

    2018-03-07 01:11:59,,2,
    關(guān)鍵詞:煙道不飽和吸附劑

    ,,2, ,

    (1.浙江理工大學(xué)理學(xué)院,杭州 310018;2.北京師范大學(xué)附屬杭州中學(xué),杭州 310002;3.杭州電子科技大學(xué)材料與環(huán)境工程學(xué)院,杭州 310018)

    0 引 言

    隨著經(jīng)濟(jì)快速發(fā)展,煙氣中日益劇增的SO2排放量對(duì)環(huán)境產(chǎn)生了一系列危害(如形成酸雨、腐蝕建筑物等),此外SO2排放還對(duì)人類(lèi)和動(dòng)植物健康產(chǎn)生了極大的危害。因此有效的捕獲及分離煙道氣中SO2對(duì)改善人類(lèi)生活質(zhì)量具有重要意義。目前,煙氣脫硫是工業(yè)上常用的脫硫方式之一,也是控制SO2排放最有效途徑[1]。按照脫硫方式和產(chǎn)物形態(tài)的不同,煙氣脫硫主要分為濕法脫硫和干法脫硫兩大類(lèi)。濕法脫硫主要包括濕法吸收和濕法氧化,其中濕法吸收具有脫硫反應(yīng)速度快、裝置運(yùn)行可靠性高等優(yōu)點(diǎn),在工業(yè)上應(yīng)用廣泛,但設(shè)備龐大,運(yùn)營(yíng)費(fèi)用高且易造成二次污染;而濕法氧化脫硫效率高,但反應(yīng)緩慢。干法脫硫利用粉狀或粒狀吸收劑、固體吸附劑或催化劑對(duì)硫化物進(jìn)行吸附,以達(dá)到脫硫目的。其作用主要依賴(lài)固體脫硫劑的孔隙度和滲透率等物理特性。干法脫硫主要采用固體吸附劑實(shí)現(xiàn),常用的固體吸附劑有活性炭[2]、沸石[3]、Ca-吸附劑[4-6]、金屬氧化物[7-9]等。干法脫硫主要包括活性炭法、氧化鐵等在內(nèi)的金屬化合物法、分子篩法在內(nèi)的非金屬礦物吸收法[10]、變壓吸附(PSA)法、低溫分離法和膜分離法等[11]。與濕法脫硫相比,該方法因低能耗、操作簡(jiǎn)單、對(duì)設(shè)備的腐蝕程度低等優(yōu)勢(shì)得到了廣泛應(yīng)用,但是也有一些不足之處,主要包括脫硫效率較低、設(shè)備投資大以及吸附劑再生性差等。因此開(kāi)發(fā)更多效率高、成本低、再生性好的固體吸附劑具有重要意義。

    近年來(lái),許多研究者致力于研究含不飽和金屬位點(diǎn)的MOFs材料對(duì)硫化物吸附與分離的影響。例如,Britt等[36]利用動(dòng)態(tài)突破法測(cè)定了室溫下6種MOFs材料對(duì)SO2、Cl2、NH3等有害氣體的捕獲與分離性能。研究結(jié)果表明,SO2在含Zn不飽和位點(diǎn)的MOF-74中動(dòng)態(tài)吸附量(約3.1 mol/kg)最佳,SO2在MOF-74中的吸附量相比于其在BPL型活性炭中增加了6倍;而SO2在含Cu不飽和位點(diǎn)的MOF-199 (也稱(chēng)Cu-BTC)中的吸附量與其在BPL型活性炭中的吸附量相近。此外,其他幾種有害氣體在含不飽和金屬位點(diǎn)的MOFs材料中的吸附量也有所改善,尤其在MOF-199表現(xiàn)更為顯著。該課題組還探討了煙道氣中SO2在這6種MOFs材料對(duì)中的分離行為,研究結(jié)果表明,MOF-74對(duì)煙道氣中SO2具有最優(yōu)的分離效果,其次為MOF-199,即MOFs材料含有的不飽和金屬位點(diǎn)有利于提高SO2的選擇性。這是由于不飽和金屬位點(diǎn)可以充當(dāng)路易斯酸位點(diǎn),有效增強(qiáng)其與SO2之間的相互作用,從而導(dǎo)致SO2在其中的吸附量較大,進(jìn)而提高了該材料對(duì)SO2的分離性能。Glover等[37]利用動(dòng)態(tài)突破實(shí)驗(yàn)探討了SO2、NH3等有害氣體在含不飽和金屬位點(diǎn)的M-MOF-74(M=Zn、Co、Ni以及Mg)系列材料中的吸附/分離行為,研究發(fā)現(xiàn),SO2在M-MOF-74(M=Zn、Co、Ni以及Mg)系列材料中飽和吸附量從大到小依次為:Mg-MOF-74、Co-MOF-74、Zn-MOF-74、Ni-MOF-74,其中在Mg-MOF-74中的最大動(dòng)態(tài)吸附量為1.60 mol/kg。此外,該課題組還利用動(dòng)態(tài)突破曲線探討室溫下SO2/N2混合氣在M-MOF-74系列材料中的分離性能,研究結(jié)果表明,Mg-MOF-74對(duì)混合氣體中SO2具有較好的分離效果。由此可見(jiàn),MOFs材料中不飽和金屬位點(diǎn)有利于提高該材料對(duì)硫化物的吸附與分離性能。

    近年來(lái),隨著計(jì)算機(jī)技術(shù)的迅速發(fā)展,分子模擬方法廣泛用于研究MOFs材料對(duì)氣體的吸附與分離。例如,Sun等[38]和Peng等[39]采用分子模擬方法研究了MOFs和沸石材料對(duì)煙道氣中SO2和NOx的分離性能,研究結(jié)果表明:SO2在含不飽和金屬位點(diǎn)的Mg-MOF-74中的工作容量最大,Mg-MOF-74對(duì)煙道氣中SO2具有較優(yōu)的分離效果。最近,Yu等[40]則進(jìn)一步采用巨正則蒙特卡羅模擬(Grand canonical Monte Carlo,GCMC)方法評(píng)估了煙道氣中的SO2對(duì)Cu-BTC捕獲CO2的影響,研究結(jié)果發(fā)現(xiàn):高壓下煙道氣中的SO2會(huì)減小Cu-BTC對(duì)CO2/N2的選擇性,這是由于SO2與Cu-BTC之間的相互作用更大所致。

    綜上所述,目前主要通過(guò)測(cè)定吸附等溫線及動(dòng)態(tài)突破實(shí)驗(yàn)評(píng)估MOFs材料對(duì)硫化物的吸附分離性能,對(duì)于MOFs材料用于SO2的吸附與分離的研究仍不夠深入,難以從分子尺度上獲得微觀性質(zhì)及相互作用機(jī)制,例如MOFs材料與硫化物之間的相互作用、硫化物在MOFs材料中吸附位點(diǎn)等問(wèn)題目前并不清晰。近年發(fā)展起來(lái)的分子模擬技術(shù)可以有效評(píng)估MOFs材料對(duì)硫化物的吸附與分離性能,還能進(jìn)一步揭示MOFs材料與吸附質(zhì)分子之間的相互作用等微觀信息,獲得MOFs材料的結(jié)構(gòu)與吸附分離性能之間的構(gòu)效關(guān)系。因此本文采用分子模擬來(lái)探討含有不飽和金屬位點(diǎn)的金屬有機(jī)骨架材料對(duì)煙道氣中SO2的吸附與分離行為,揭示金屬離子對(duì)煙道氣中SO2吸附/分離性能的影響,為MOFs材料設(shè)計(jì)提供理論指導(dǎo)。

    1 模型與模擬過(guò)程

    1.1 MOFs材料模型及結(jié)構(gòu)性質(zhì)

    1.1.1 MOFs材料模型

    選擇不同金屬離子的M-DOBDCs(M=Zn、Co、Ni以及Mg)系列材料為吸附劑,探討金屬離子對(duì)SO2吸附分離性能的影響。其中,Zn-DOBDC是最早報(bào)道的含有不飽和金屬位點(diǎn)的M-DOBDCs系列材料,M-DOBDCs系列材料由有機(jī)配體2,5-二羥基-1,4-對(duì)苯二甲酸與二價(jià)金屬離子MⅡ(Mg,Co,Ni以及Zn)自組裝而成,該類(lèi)材料具有一維六邊形的孔道結(jié)構(gòu)。上述MOFs材料的晶體結(jié)構(gòu)從英國(guó)劍橋數(shù)據(jù)庫(kù)(Cambridge crystallographic data centre,CCDC)獲得,除去溶劑分子等多余信息,依據(jù)晶體結(jié)構(gòu)數(shù)據(jù)構(gòu)建MOFs結(jié)構(gòu),并利用分子力學(xué)方法對(duì)該MOFs材料的結(jié)構(gòu)進(jìn)行優(yōu)化,具體優(yōu)化過(guò)程為:優(yōu)化前先固定MOFs材料晶體結(jié)構(gòu)所有原子的位置,然后以MOFs材料的晶胞參數(shù)和原子位置為變量,采用綜合法(Smart minimizer)進(jìn)行能量最小化計(jì)算,從而得到較穩(wěn)定的晶體結(jié)構(gòu)。

    1.1.2 MOFs材料的結(jié)構(gòu)性質(zhì)

    吸附質(zhì)在MOFs材料中的吸附行為與MOFs材料的結(jié)構(gòu)性質(zhì)(孔隙率、比表面積、孔徑分布)密切相關(guān),本文所探討M-DOBDCs系列MOFs材料的結(jié)構(gòu)性質(zhì)(孔隙率、比表面積、孔徑分布)見(jiàn)表1,從表中可以看出Mg-DOBDC的比表面積最大,達(dá)1718 m2/g。

    表1 MOFs材料的結(jié)構(gòu)性質(zhì)

    1.2 勢(shì)能模型與力場(chǎng)參數(shù)

    吸附質(zhì)分子與吸附劑以及吸附質(zhì)分子之間的相互作用能是通過(guò)描述范德華非鍵作用的Lennard-Jones(LJ)相互作用勢(shì)能與靜電作用的庫(kù)侖力勢(shì)能加和獲得的,具體的數(shù)學(xué)表達(dá)式如式(1)所示:

    (1)

    其中:rij代表原子i和原子j之間的距離,εij和σij分別代表原子的LJ勢(shì)能參數(shù),qi和qj分別為原子i和原子j所帶的電荷;ε0為有效介電常數(shù)。不同原子之間的LJ勢(shì)能參數(shù)采用Lorentz-Berthelot(LB)混合規(guī)則計(jì)算[41]。

    MOFs材料中骨架原子的LJ作用參數(shù)來(lái)自DREIDING[42]力場(chǎng),但是對(duì)于DREIDING力場(chǎng)缺失的原子,如Co、Mg、Ni原子,其作用參數(shù)來(lái)自UFF[43]力場(chǎng)。MOFs材料各骨架原子的LJ作用參數(shù)及所屬力場(chǎng)見(jiàn)表2。MOFs材料中骨架各類(lèi)型的原子的電荷基于仲崇立課題組[44]發(fā)展的基于原子連接性貢獻(xiàn)法(Connectivity-based atom contribution,CBAC)獲得, M-DOBDCs材料相同原子所處的環(huán)境不同,電荷不同,本文對(duì)M-DOBDCs材料中不同原子類(lèi)型進(jìn)行分類(lèi)并編號(hào),結(jié)果如圖1所示,M-DOBDCs骨架結(jié)構(gòu)中不同類(lèi)型原子的電荷見(jiàn)表3。對(duì)于吸附質(zhì)分子N2和CO2分子,其電荷和LJ相互作用勢(shì)能參數(shù)均來(lái)自TraPPE力場(chǎng),該力場(chǎng)對(duì)于氣-液共存體系有很好的適用性[45],對(duì)于SO2而言,本文采用Okolic等[46]發(fā)現(xiàn)的三位點(diǎn)的LJ模型來(lái)描述。在模擬過(guò)程中,SO2、N2和CO2采用三位點(diǎn)LJ位點(diǎn)勢(shì)能模型,其LJ相互作用勢(shì)能參數(shù)和原子電荷見(jiàn)表4。

    表2 M-DOBDCs材料中骨架原子的LJ作用參數(shù)及所屬力場(chǎng)

    圖1 M-DOBDCs材料中不同類(lèi)型原子的示意圖注:數(shù)字1-4表示C原子,5表示H原子,6-8表示O原子,9表示金屬原子(M)。

    MOFs材料原子編號(hào)123456789Mg?DOBDC0.800(C)-0.260(C)0.490(C)-0.280(C)0.200(H)-0.710(O)-0.910(O)-0.780(O)1.46(Mg)Co?DOBDC0.830(C)-0.290(C)0.320(C)-0.110(C)0.180(H)-0.680(O)-0.650(O)-0.730(O)1.14(Co)Ni?DOBDC0.710(C)-0.240(C)0.360(C)-0.200(C)0.170(H)-0.680(O)-0.620(O)-0.680(O)1.10(Ni)Zn?DOBDC0.770(C)-0.290(C)0.330(C)-0.150(C)0.170(H)-0.660(O)-0.670(O)-0.700(O)1.21(Zn)

    表4 SO2、N2和CO2 LJ相互作用勢(shì)能參數(shù)和原子電荷

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

    本文采用巨正則蒙特卡羅模擬方法來(lái)研究含不飽和金屬位點(diǎn)的M-DOBDCs系列材料對(duì)SO2的吸附與分離行為,混合物中MSO2∶MN2∶MCO2設(shè)定為0.002∶0.848∶0.150,與工業(yè)中煙道氣一致。模擬過(guò)程中,溫度設(shè)置為313 K。煙道氣中SO2的分離過(guò)程采用真空變壓吸附(Vacuum pressure swing adsorption,VSA),在VSA過(guò)程中吸附壓力為1.0 bar,脫附壓力為0.1 bar。LJ相互作用的截?cái)喟霃綖?.28 nm,模擬盒子的邊長(zhǎng)大于截?cái)喟霃降?倍。對(duì)于每一個(gè)狀態(tài)點(diǎn),GCMC模擬的步長(zhǎng)設(shè)定為1×105步獲得平衡,再運(yùn)行1×105步進(jìn)行結(jié)果統(tǒng)計(jì)。模擬過(guò)程中,為了評(píng)估含有不同金屬離子的M-DOBDCs系列MOF材料對(duì)煙道氣中SO2的分離性能,計(jì)算了該微孔材料對(duì)混合氣中SO2的相對(duì)選擇性系數(shù),該系數(shù)定義如下:

    (2)

    式中:xa和ya分別為物質(zhì)a在吸附相和本體相中的摩爾分?jǐn)?shù)。工作容量(Working capacity,W)是評(píng)估分離性能的另一個(gè)重要參數(shù),為了評(píng)估不飽和金屬位點(diǎn)的M-DOBDCs系列材料對(duì)煙道氣中SO2的分離性能,計(jì)算了VSA過(guò)程中該微孔材料的SO2工作容量(WSO2),用數(shù)學(xué)表達(dá)式表示:

    (3)

    吸附熱可以反映吸附質(zhì)和吸附劑及吸附質(zhì)與吸附質(zhì)之間的相互作用的強(qiáng)弱,在模擬過(guò)程中統(tǒng)計(jì)了吸附熱(Qst),吸附熱的數(shù)學(xué)統(tǒng)計(jì)式如式(4)所示:

    (4)

    式中:V表示體系的總勢(shì)能,N表示蒙特卡羅模擬的分子數(shù),尖括號(hào)表示求系綜平均。

    2 結(jié)果與討論

    2.1 金屬離子對(duì)SO2吸附性能的影響

    為揭示金屬離子對(duì)SO2吸附性能的影響,模擬了313 K時(shí)SO2分子在M-DOBDCs(M=Zn、Co、Ni以及Mg)材料中的吸附等溫線,結(jié)果如圖2所示。從圖2中可發(fā)現(xiàn),吸附等溫線變化趨勢(shì)一致,SO2在M-DOBDCs中的吸附量隨著壓力的增加,首先迅速增加,然后逐步達(dá)到平衡。當(dāng)壓力小于0.1 bar時(shí),SO2在不同金屬離子的M-DOBDCs系列材料中的吸附順序從大到小依次為Zn-DOBDC、Ni-DOBDC、Co-DOBDC、Mg-DOBDC,這可能是因?yàn)榈蛪合耑n-DOBDC與SO2之間的相互作用較大的緣故。

    圖2 313 K時(shí)SO2在不同金屬離子的M-DOBDCs系列材料中的吸附等溫線

    吸附熱大小可以反映吸附質(zhì)與吸附劑間相互作用的強(qiáng)弱,為了揭示吸附質(zhì)與吸附劑間相互作用,統(tǒng)計(jì)了313 K及低壓下SO2在DOBDCs系列材料中的吸附熱,并與該狀態(tài)下的吸附量進(jìn)行關(guān)聯(lián),結(jié)果如圖3所示。從圖3中可發(fā)現(xiàn),在較低壓力時(shí)SO2在M-DOBDCs系列材料中的吸附量隨著其在該材料中的吸附熱的增加呈線性增大的趨勢(shì),吸附熱大小順序?yàn)椋篫n-DOBDC、Ni-DOBDC、Co-DOBDC、Mg-DOBDC,并且SO2在Zn-DOBDC中的吸附熱高達(dá)40 kJ/mol,這是因?yàn)閆nⅡ與DOBDC配體中的1,4-對(duì)苯二甲酸根離子配位時(shí)形成的“夾角”較大(ZnⅡ、NiⅡ、CoⅡ、MgⅡ與上述羧酸根離子配位時(shí)形成的“夾角”分別為153.5°、133.2°、126.4°和91.7°),從而更有利于V型分子SO2接近Zn不飽和位點(diǎn)。該吸附熱變化趨勢(shì)與圖2所示吸附量變化趨勢(shì)相一致,即SO2分子與MOFs材料之間的相互作用越大,吸附熱越大,其在MOFs材料中的吸附量越大。

    圖3 313 K及0.05 bar時(shí)SO2在M-DOBDCs材料中的吸附熱與吸附量之間的關(guān)系曲線

    圖4 313 K及1.0 bar時(shí)不同金屬離子的M-DOBDCs材料的結(jié)構(gòu)性質(zhì)與SO2吸附量的關(guān)系曲線

    高壓下SO2在M-DOBDCs材料中的吸附量接近飽和,飽和吸附量的差異主要與MOFs材料的結(jié)構(gòu)性質(zhì)(比表面積和孔體積)有關(guān)。313 K及1.0 bar時(shí)SO2在M-DOBDCs材料中的吸附量與該類(lèi)材料結(jié)構(gòu)性質(zhì)(比表面積和孔體積)之間的關(guān)系如圖4所示。從圖4不難發(fā)現(xiàn),SO2在M-DOBDCs材料中的吸附量隨著該類(lèi)材料的比表面積及孔體積呈線性增加,即M-DOBDCs材料的比表面積越大、孔體積越高,SO2在其中的吸附性能越好。

    2.2 金屬離子對(duì)煙道氣中SO2分離性能的影響

    MOFs材料除應(yīng)具備較優(yōu)的吸附性能外,還需具備較好的分離性能。為評(píng)估M-DOBDCs系列材料對(duì)煙道氣中SO2分離性能,揭示金屬離子對(duì)分離性能的影響,本文計(jì)算了313 K及0.1 bar時(shí)上述含不同金屬離子的MOFs材料對(duì)煙道氣中SO2的相對(duì)選擇性系數(shù),結(jié)果如圖5所示。由圖5可以發(fā)現(xiàn),M-DOBDCs系列材料對(duì)SO2的吸附選擇性從大到小依次為:Zn-DOBDC、Ni-DOBDC、Co-DOBDC、Mg-DOBDC,其中Zn-DOBDC對(duì)煙道氣中SO2的分離性能最優(yōu),相對(duì)選擇性系數(shù)約為39.3。SO2在M-DOBDCs中的吸附行為受煙道氣中CO2及N2組分的影響,進(jìn)而影響M-DOBDCs材料對(duì)SO2選擇性。M-DOBDCs系列材料對(duì)SO2分離性能的差異主要是由于煙道氣中N2、CO2與SO2存在競(jìng)爭(zhēng)所致,N2和CO2在不同金屬離子的M-DOBDCs材料中吸附行為不同。

    圖5 313 K及0.1 bar下,不同金屬離子的M-DOBDCs材料對(duì)煙道氣中SO2選擇性(MSO2∶MN2∶MCO2=0.002∶0.848∶0.150)

    為解釋不同金屬離子對(duì)煙道氣中SO2分離性能的影響,本文中進(jìn)一步統(tǒng)計(jì)了煙道氣中SO2、N2和CO2各組分在M-DOBDCs材料中的吸附量,結(jié)果如圖6所示。研究結(jié)果表明N2在M-DOBDCs材料中的吸附量相當(dāng)且較低,這是因?yàn)镹2分子為非極性分子,不能有效的與不飽和金屬位點(diǎn)結(jié)合,導(dǎo)致低壓下N2分子在M-DOBDCs材料中吸附量較低。而CO2在M-DOBDCs中的吸附量則存在較大差異,例如CO2在Zn-DOBDC中的吸附量最小,在Mg-DOBDC吸附量最大,即在Zn-DOBDC體系中CO2與SO2之間的競(jìng)爭(zhēng)最小。此外由于煙道氣中的N2分子在M-DOBDCs中吸附量較低且數(shù)值相當(dāng),其對(duì)煙道氣中SO2吸附影響不大,而SO2在Zn-DOBDC體系中的吸附量較大,且對(duì)其它氣體吸附量不高,有效的提高了其對(duì)SO2分離效果。結(jié)合表5中0.1 bar時(shí)煙道氣各組分(SO2、N2和CO2)在M-DOBDCs材料中的吸附熱可知,N2在M-DOBDCs材料中的吸附熱較低(約為11.5~18.6 kJ/mol),CO2在Zn-DOBDC材料中的吸附熱(24.5 kJ/mol)小于CO2在其他M-DOBDCs材料中的吸附熱,而SO2在Zn-DOBDC中的吸附熱(46.7 kJ/mol)最大,即SO2與Zn-DOBDC之間較大的相互作用,增加了其對(duì)SO2的捕獲,有效提高了Zn-DOBDC對(duì)煙道氣中SO2的分離選擇性。該結(jié)果從相互作用角度也進(jìn)一步證明在Zn-DOBDC體系中CO2與SO2之間的競(jìng)爭(zhēng)最小。

    圖6 313 K及0.1 bar下,煙道氣中SO2,N2和CO2在不同金屬離子的M-DOBDCs材料中的吸附量(MSO2∶MN2∶MCO2=0.002∶0.848∶0.150)

    為揭示SO2、N2和CO2混合氣在M-DOBDCs材料的吸附機(jī)制,以Zn-DOBDC為例,在模擬過(guò)程中統(tǒng)計(jì)了SO2、N2和CO2混合氣在該材料中的密度分布,結(jié)果如圖7所示,SO2、N2和CO2主要吸附在不飽和金屬位點(diǎn)周?chē)?,可?jiàn)CO2、N2及SO2在吸附位點(diǎn)處存在競(jìng)爭(zhēng),但是CO2的四極矩較大,其與Zn-DOBDC之間的相互作用較大,所以Zn-DOBDC中煙道氣中存在的CO2對(duì)SO2分離起主要影響。對(duì)比圖7(a)和圖7(c),不難發(fā)現(xiàn)SO2與CO2在Zn-DOBDC中吸附構(gòu)象有所不同,從而造成了其密度分布圖的差異。其中由于SO2分子構(gòu)型為V形分子,Zn-DOBDC材料中金屬離子與其配位的羧酸根離子形成的“夾角”約為153.5°,SO2分子能夠以多種構(gòu)象靠近不飽和金屬原子,增加了SO2與不飽和金屬離子接觸概率,因此SO2分子在Zn-DOBDC材料中的密度分布圖為環(huán)繞狀(圖7(a))。而CO2為直線型分子,在吸附過(guò)程中為保證吸附質(zhì)和吸附劑分子之間具有較大的相互作用,只能以一端“進(jìn)攻”形式靠近不飽和金屬位點(diǎn)(圖7(d)),因此CO2分子在Zn-DOBDC材料中吸附位點(diǎn)較為集中(圖7(c))。隨著M-DOBDCs材料中金屬離子的改變,金屬離子與羧酸根離子形成的“夾角”逐漸變小,SO2越來(lái)越難靠近吸附位點(diǎn),而CO2分子與金屬離子接觸概率則不會(huì)因?yàn)樵揗OF材料中金屬離子和有機(jī)配體之間的夾角改變而改變,反而會(huì)因?yàn)閵A角的逐漸變小,CO2能夠與有機(jī)配體周?chē)嗟脑咏佑|,相互作用增強(qiáng)(例如CO2分子與Mg-DOBDC之間相互作用較強(qiáng),表5)。

    表5 0.1 bar時(shí)煙道氣中SO2、N2和CO2在M-DOBDCs材料中的吸附熱

    在實(shí)際應(yīng)用中,材料除了應(yīng)具有較好的吸附性能和優(yōu)良的選擇性之外,還需要具有良好的再生性能和較大的處理量。工作容量是反映材料處理量大小及可再生性能優(yōu)劣的參數(shù),為此文中計(jì)算了不同M-DOBDCs材料的SO2工作容量,結(jié)果如圖8所示,可以看出VSA過(guò)程中,Mg-DOBDC、Zn-DOBDC具有較高的SO2工作容量,結(jié)合圖4可知Zn-DOBDC對(duì)煙道氣中SO2的選擇性最大,可見(jiàn)Zn-DOBDC適用于VSA過(guò)程中煙道氣中SO2的分離。

    圖7 313K及0.1 bar時(shí)煙道氣中SO2、N2和CO2在Zn-DOBDC材料中的密度分布以及SO2和CO2在Zn-DOBDC中的瞬時(shí)吸附構(gòu)型

    圖8 313 K時(shí)VSA過(guò)程中不同金屬離子的M-DOBDCs材料的SO2工作容量

    3 結(jié) 論

    本文采用GCMC模擬探討了含不飽和金屬位點(diǎn)的M-DOBDCs材料對(duì)SO2的捕獲性能,研究結(jié)果表明,低壓下,SO2在M-DOBDCs材料中的吸附性能主要受吸附質(zhì)和吸附劑之間的相互作用影響,SO2與MOFs材料之間的相互作用越大,SO2在其中的吸附性能越好;高壓下,SO2的吸附性能與M-DOBDCs材料的比表面積、孔體積等結(jié)構(gòu)因素密切相關(guān),比表面積、孔體積越大,SO2在其中的吸附性能越好。此外,論文進(jìn)一步揭示了不同金屬離子的M-DOBDCs材料對(duì)煙道氣中SO2的分離性能,研究結(jié)果表明,煙道氣中N2在不同金屬離子的M-DOBDCs材料中吸附量相當(dāng)且較低,不同金屬離子對(duì)煙道氣分離性能差異主要由CO2與SO2在不同M-DOBDCs中競(jìng)爭(zhēng)效應(yīng)不同所致。M-DOBDCs材料對(duì)煙道氣中SO2分離性能和MOFs材料與吸附質(zhì)之間的相互作用有關(guān)。金屬離子和有機(jī)配體中的羧酸根離子配位時(shí)所形成的夾角越大,SO2與M-DOBDCs中的不飽和金屬位點(diǎn)之間接觸概率越大,相互作用越強(qiáng),而對(duì)于CO2分子,角度越小則越利于CO2分子靠近不飽和金屬離子,最終導(dǎo)致該材料對(duì)煙道氣中SO2的分離性能越好。對(duì)于M-DOBDCs材料,綜合考慮其對(duì)SO2的吸附量和相對(duì)選擇系數(shù)及工作容量,可知Zn-DOBDC是可用于煙道氣中VSA過(guò)程SO2分離的備選材料。

    [1] 薛菲,石勁松.煙氣脫硫技術(shù)現(xiàn)狀與發(fā)展建議[J].江蘇化工,2001,29(2):32-35.

    [2] Davini P. SO2and NOxadsorption properties of activated carbons obtained from a pitch containing iron derivatives[J]. Carbon,2001,39(14):2173-2179.

    [3] Deng H, Yi H, Tang X, et al. Adsorption equilibrium for sulfur dioxide, nitric oxide, carbon dioxide, nitrogen on 13X and 5A zeolites[J]. Chemical Engineering Journal,2012,188:77-85.

    [4] Ridha F N, Manovic V, Macchi A, et al. The effect of SO2on CO2capture by CaO-based pellets prepared with a kaolin derived Al(OH)3binder[J]. Applied Energy,2012,92:415-420.

    [5] Lu H, Smirniotis P G. Calcium oxide doped sorbents for CO2uptake in the presence of SO2at high temperatures[J]. Industrial & Engineering Chemistry Research,2009,48(11):5454-5459.

    [6] Sun P, Grace J R, Lim C J, et al. Removal of CO2by calcium-based sorbents in the presence of SO2[J]. Energy & Fuels,2007,21(1):163-170.

    [7] Peterson G W, Karwacki C J, Feaver W B, et al. Zirconium hydroxide as a reactive substrate for the removal of sulfur dioxide[J]. Industrial & Engineering Chemistry Research,2009,48(4):1694-1698.

    [8] Seredych M, Bandosz T J. Effects of surface features on adsorption of SO2on graphite oxide/Zr(OH)4composites[J]. The Journal of Physical Chemistry C,2010,114(34):14552-14560.

    [9] Seredych M, Mabayoje O, Bandosz T J. Involvement of water and visible light in the enhancement in SO2adsorption at ambient conditions on the surface of zinc(hydr) oxide/graphite oxide composites[J]. Chemical Engineering Journal,2013,223:442-453.

    [10] 舒欣,李春光,任家君,等.國(guó)內(nèi)天然氣干法脫硫劑的比較研究[J].廣東化工,2010(1):89-90.

    [11] 田正山.新型濕法脫硫催化劑的開(kāi)發(fā)與性能研究[D].鄭州:鄭州大學(xué),2006.

    [12] Ou H C, Long J R, Yaghi O M. Introduction to metal-organic frameworks[J]. Chemical Reviews,2012,112(2):673-674.

    [13] Gygi D, Bloch E D, Mason J A, et al. Hydrogen storage in the expanded pore metal-organic frameworks M2(dobpdc)(M=Mg, Mn, Fe, Co, Ni, Zn)[J]. Chemistry of Materials,2016,28(4):1128-1138.

    [14] Savage M, Da Silva I, Johnson M, et al. Observation of binding and rotation of methane and hydrogen within a functional metal-organic framework[J]. Journal of the American Chemical Society,2016,138(29):9119-9127.

    [15] Vieira Soares C, Damasceno Borges D, Wiersum A, et al. Adsorption of small molecules in the porous zirconium-based metal organic framework MIL-140A (Zr): a joint computational-experimental approach[J]. The Journal of Physical Chemistry C,2016,120(13):7192-7200.

    [16] Zhai Q G, Bu X, Mao C, et al. Systematic and dramatic tuning on gas sorption performance in heterometallic metal-organic frameworks[J]. Journal of the American Chemical Society,2016,138(8):2524-2527.

    [17] Chen Y P, Liu Y, Liu D, et al. Direct measurement of adsorbed gas redistribution in metal-organic frameworks[J]. Journal of the American Chemical Society,2015,137(8):2919-2930.

    [18] Foo M L, Matsuda R, Hijikata Y, et al. An adsorbate discriminatory gate effect in a flexible porous coordination polymer for selective adsorption of CO2over C2H2[J]. Journal of the American Chemical Society,2016,138(9):3022-3030.

    [19] Luo F, Yan C, Dang L, et al. UTSA-74: A MOF-74 isomer with two accessible binding sites per metal center for highly selective gas separation[J]. Journal of the American Chemical Society,2016,138(17):5678-5684.

    [20] Fu J, Das S, Xing G, et al. Fabrication of COF-MOF composite membranes and their highly selective separation of H2/CO2[J]. Journal of the American Chemical Society,2016,138(24):7673-7680.

    [21] Li B, Liu J, Nie Z, et al. Metal-organic frameworks as highly active electrocatalysts for high-energy density, aqueous zinc-polyiodide redox flow batteries[J]. Nano Letters,2016,16(7):4335-4340.

    [22] Kim D, Whang D R, Park S Y. Self-healing of molecular catalyst and photosensitizer on metal-organic framework: robust molecular system for photocatalytic H2evolution from water[J]. Journal of the American Chemical Society,2016,138(28):8698-8701.

    [23] Li P Z, Wang X J, Liu J, et al. A Triazole-containing metal-organic framework as a highly effective and substrate size-dependent catalyst for CO2conversion[J]. Journal of the American Chemical Society,2016,138(7):2142-2145.

    [24] He H, Song Y, Sun F, et al. Sorption properties and nitroaromatic explosives sensing based on two isostructural metal-organic frameworks[J]. Crystal Growth & Design,2015,15(4):2033-2038.

    [25] Shustova N B, Cozzolino A F, Reineke S, et al. Selective turn-on ammonia sensing enabled by high-temperature fluorescence in metal-organic frameworks with open metal sites[J]. Journal of the American Chemical Society,2013,135(36):13326-13329.

    [26] Wang Y, Wang L, Chen H, et al. Fabrication of highly sensitive and stable hydroxylamine electrochemical sensor based on gold nanoparticles and metal-metalloporphyrin framework modified electrode[J]. ACS Applied Materials & Interfaces,2016,8(28):18173-18181.

    [27] Levine D J, Runcevski T, Kapelewski M T, et al. Olsalazine-based metal-organic frameworks as biocompatible platforms for H2adsorption and drug delivery[J]. Journal of the American Chemical Society,2016,138(32):10143-10150

    [28] Chowdhuri A R, Singh T, Ghosh S K, et al. Carbon dots embedded magnetic nanoparticles@chitosan@ metal organic framework as a nanoprobe for pH sensitive targeted anticancer drug delivery[J]. ACS Applied Materials & Interfaces,2016,8(26):16573-16583.

    [29] Yang Y, Hu Q, Zhang Q, et al. A large capacity cationic metal-organic framework nanocarrier for physiological pH responsive drug delivery[J]. Molecular Pharmaceutics,2016,13(8):2782-2786.

    [30] Bloch E D, Queen W L, Chavan S, et al. Gradual release of strongly bound nitric oxide from Fe2(NO)2(dobdc)[J]. Journal of the American Chemical Society,2015,137(10):3466-3469.

    [31] Sumida K, Rogow D L, Mason J A, et al. Carbon dioxide capture in metal-organic frameworks[J]. Chemical Reviews,2011,112(2):724-781.

    [32] Suh M P, Park H J, Prasad T K, et al. Hydrogen storage in metal-organic frameworks[J]. Chemical Reviews,2011,112(2):782-835.

    [33] Zhou W, Wu H, Yildirim T. Enhanced H2adsorption in isostructural metal-organic frameworks with open metal sites: strong dependence of the binding strength on metal ions[J]. Journal of the American Chemical Society,2008,130(46):15268-15269.

    [34] Wu H, Zhou W, Yildirim T. High-capacity methane storage in metal-organic frameworks M2(dhtp): the important role of open metal sites[J]. Journal of the American Chemical Society,2009,131(13):4995-5000.

    [35] Caskey S R, Wong-Foy A G, Matzger A J. Dramatic tuning of carbon dioxide uptake via metal substitution in a coordination polymer with cylindrical pores[J]. Journal of the American Chemical Society,2008,130(33):10870-10871.

    [36] Britt D, Tranchemontagne D, Yaghi O M. Metal-organic frameworks with high capacity and selectivity for harmful gases[J]. Proceedings of the National Academy of Sciences,2008,105(33):11623-11627.

    [37] Glover T G, Peterson G W, Schindler B J, et al. MOF-74 building unit has a direct impact on toxic gas adsorption[J]. Chemical Engineering Science,2011,66(2):163-170.

    [38] Sun W, Lin L C, Peng X, et al. Computational screening of porous metal-organic frameworks and zeolites for the removal of SO2and NOxfrom flue gases[J]. AIChE Journal,2014,60(6):2314-2323.

    [39] Peng X, Cao D. Computational screening of porous carbons, zeolites, and metal organic frameworks for desulfurization and decarburization of biogas, natural gas, and flue gas[J]. AIChE Journal,2013,59(8):2928-2942.

    [40] Yu J, Ma Y, Balbuena P B. Evaluation of the impact of H2O, O2, and SO2on postcombustion CO2capture in metal-organic frameworks[J]. Langmuir,2012,28(21):8064-8071.

    [41] Yazaydin A O, Snurr R Q, Park T H, et al. Screening of metal-organic frameworks for carbon dioxide capture from flue gas using a combined experimental and modeling approach[J]. Journal of the American Chemical Society,2009,131(51):18198-18199.

    [42] Mayo S L, Olafson B D, Goddard W A. DREIDING: a generic force field for molecular simulations[J]. Journal of Physical Chemistry,1990,94(26):8897-8909.

    [43] Rappé A K, Casewit C J, Colwell K S, et al. UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations[J]. Journal of the American Chemical Society,1992,114(25):10024-10035.

    [44] Xu Q, Zhong C. A general approach for estimating framework charges in metal-organic frameworks[J]. The Journal of Physical Chemistry C,2010,114(11):5035-5042.

    [45] Potoff J J, Siepmann J I. Vapor-liquid equilibria of mixtures containing alkanes, carbon dioxide, and nitrogen[J]. AIChE Journal,2001,47(7):1676-1682.

    [46] Okolic F, Guissani Y, Guillot B. Molecular dynamics simulations of thermodynamic and structural properties of liquid SO2[J]. Molecular Physics,1985,56:239-253.

    猜你喜歡
    煙道不飽和吸附劑
    固體吸附劑脫除煙氣中SOx/NOx的研究進(jìn)展
    化工管理(2022年13期)2022-12-02 09:21:52
    用于空氣CO2捕集的變濕再生吸附劑的篩選與特性研究
    能源工程(2021年1期)2021-04-13 02:05:50
    昆鋼公司120t轉(zhuǎn)爐固定段煙道漏水分析及治理
    昆鋼科技(2021年6期)2021-03-09 06:10:20
    加熱爐燃燒煙道系統(tǒng)的改造
    n-3多不飽和脂肪酸改善糖脂代謝的研究進(jìn)展
    焦?fàn)t煙道氣治理的難點(diǎn)分析
    超聲波促進(jìn)合成新型吡咯α,β-不飽和酮
    兩個(gè)alpha,beta-不飽和酮藥物的合成及抗腫瘤活性研究
    ω-3多不飽和脂肪酸對(duì)人胃腺癌細(xì)胞系A(chǔ)GS的作用
    茶籽殼吸附劑的制備與表征
    国产在视频线在精品| 欧美在线一区亚洲| 插阴视频在线观看视频| 人妻久久中文字幕网| 国产色爽女视频免费观看| 99久久九九国产精品国产免费| 日本免费一区二区三区高清不卡| 我的女老师完整版在线观看| 99久久精品国产国产毛片| 亚洲av不卡在线观看| 午夜日韩欧美国产| 国产精品久久久久久亚洲av鲁大| 成年女人永久免费观看视频| 欧美一区二区亚洲| 好男人在线观看高清免费视频| 久久久久精品国产欧美久久久| 成人一区二区视频在线观看| 国产在线精品亚洲第一网站| 亚洲真实伦在线观看| av在线蜜桃| ponron亚洲| 99国产极品粉嫩在线观看| 在线观看av片永久免费下载| 国产精品无大码| 亚洲精品日韩av片在线观看| 国产精品久久久久久精品电影| 国产高清有码在线观看视频| 午夜福利成人在线免费观看| 三级经典国产精品| 国产精品,欧美在线| 美女内射精品一级片tv| 又爽又黄a免费视频| 免费av不卡在线播放| 少妇被粗大猛烈的视频| 18禁在线无遮挡免费观看视频 | 亚洲性久久影院| 人妻制服诱惑在线中文字幕| 国产探花极品一区二区| 欧美区成人在线视频| 国产高潮美女av| 亚洲成人精品中文字幕电影| 日日干狠狠操夜夜爽| 美女高潮的动态| 日日摸夜夜添夜夜添av毛片| av国产免费在线观看| 免费观看人在逋| 97超级碰碰碰精品色视频在线观看| 麻豆一二三区av精品| 悠悠久久av| 深夜精品福利| 在线免费观看不下载黄p国产| 成人欧美大片| 午夜福利高清视频| 欧美日韩在线观看h| 成人美女网站在线观看视频| av视频在线观看入口| 日本欧美国产在线视频| 少妇的逼水好多| 大型黄色视频在线免费观看| 国产精品三级大全| 成人av在线播放网站| 免费人成在线观看视频色| 成人av一区二区三区在线看| 男女边吃奶边做爰视频| 久久精品夜色国产| 久久久精品大字幕| 国产成人freesex在线 | 精品一区二区三区视频在线| 禁无遮挡网站| 99热网站在线观看| 我的女老师完整版在线观看| 九九热线精品视视频播放| 免费观看的影片在线观看| 国内久久婷婷六月综合欲色啪| 人妻制服诱惑在线中文字幕| 在线观看66精品国产| 欧美日韩综合久久久久久| 精品国产三级普通话版| 一个人免费在线观看电影| 亚洲欧美中文字幕日韩二区| 精品熟女少妇av免费看| 国产成人精品久久久久久| 国产免费一级a男人的天堂| 波多野结衣高清无吗| 国产一区二区在线观看日韩| eeuss影院久久| 91av网一区二区| 男女下面进入的视频免费午夜| 久久久久久九九精品二区国产| 97超碰精品成人国产| 成人美女网站在线观看视频| 波多野结衣巨乳人妻| 直男gayav资源| 午夜免费男女啪啪视频观看 | 国产精品一区二区性色av| 夜夜看夜夜爽夜夜摸| 网址你懂的国产日韩在线| 久久人人爽人人爽人人片va| 免费不卡的大黄色大毛片视频在线观看 | 日本与韩国留学比较| 亚洲国产高清在线一区二区三| 看免费成人av毛片| 在线观看66精品国产| 一个人看的www免费观看视频| 伊人久久精品亚洲午夜| 日韩,欧美,国产一区二区三区 | 一a级毛片在线观看| 22中文网久久字幕| 在线观看av片永久免费下载| 亚洲一区二区三区色噜噜| 三级国产精品欧美在线观看| 亚洲色图av天堂| 三级男女做爰猛烈吃奶摸视频| 国产黄片美女视频| 看免费成人av毛片| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲,欧美,日韩| 最近最新中文字幕大全电影3| av黄色大香蕉| 日韩欧美国产在线观看| 精品国内亚洲2022精品成人| 日本一本二区三区精品| 别揉我奶头 嗯啊视频| 成年女人毛片免费观看观看9| 欧美色视频一区免费| 欧美zozozo另类| 免费大片18禁| 午夜老司机福利剧场| 欧美xxxx性猛交bbbb| 国产亚洲欧美98| 久久久久精品国产欧美久久久| 97超碰精品成人国产| 国产极品精品免费视频能看的| 久久国产乱子免费精品| 欧美性猛交黑人性爽| 久久国内精品自在自线图片| 亚洲av中文字字幕乱码综合| 精品国产三级普通话版| 国产美女午夜福利| 久久久久国产精品人妻aⅴ院| 久久精品国产亚洲av香蕉五月| 给我免费播放毛片高清在线观看| 搡老熟女国产l中国老女人| 无遮挡黄片免费观看| 色哟哟·www| 18禁裸乳无遮挡免费网站照片| 麻豆乱淫一区二区| 青春草视频在线免费观看| 一个人看视频在线观看www免费| 日本-黄色视频高清免费观看| 男女下面进入的视频免费午夜| 99九九线精品视频在线观看视频| 人妻少妇偷人精品九色| 又爽又黄a免费视频| 波多野结衣巨乳人妻| 女人被狂操c到高潮| 日韩一区二区视频免费看| 日本熟妇午夜| 男女啪啪激烈高潮av片| 色综合亚洲欧美另类图片| 综合色丁香网| 亚洲国产色片| aaaaa片日本免费| 少妇丰满av| 日韩 亚洲 欧美在线| 亚洲av免费高清在线观看| 最近2019中文字幕mv第一页| 国产不卡一卡二| 亚洲自偷自拍三级| 嫩草影院新地址| 精品人妻熟女av久视频| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久久久中文| 成人av一区二区三区在线看| 一进一出好大好爽视频| 日日干狠狠操夜夜爽| 亚洲成a人片在线一区二区| 国产午夜精品论理片| 色综合亚洲欧美另类图片| 一本一本综合久久| 久久久a久久爽久久v久久| 女同久久另类99精品国产91| 99热这里只有精品一区| 国产麻豆成人av免费视频| 国产亚洲欧美98| 亚洲av中文字字幕乱码综合| 99久久九九国产精品国产免费| 少妇人妻一区二区三区视频| 亚洲熟妇熟女久久| 免费观看的影片在线观看| 日本一本二区三区精品| 精品一区二区三区人妻视频| 国产高清激情床上av| 性插视频无遮挡在线免费观看| 免费在线观看成人毛片| 18+在线观看网站| 蜜臀久久99精品久久宅男| 国产精品人妻久久久影院| 欧美成人一区二区免费高清观看| 精品久久久久久成人av| 国产毛片a区久久久久| 日本 av在线| 一边摸一边抽搐一进一小说| 尤物成人国产欧美一区二区三区| 国产男靠女视频免费网站| 嫩草影院入口| 两性午夜刺激爽爽歪歪视频在线观看| 十八禁网站免费在线| 亚洲一区高清亚洲精品| 国内少妇人妻偷人精品xxx网站| а√天堂www在线а√下载| 亚洲av成人精品一区久久| 日韩欧美精品v在线| 亚洲va在线va天堂va国产| 少妇的逼水好多| 久久精品国产自在天天线| 久久国内精品自在自线图片| 老熟妇乱子伦视频在线观看| 免费大片18禁| 午夜免费男女啪啪视频观看 | 亚洲欧美中文字幕日韩二区| 国产精品久久久久久久电影| 女人十人毛片免费观看3o分钟| 中文字幕久久专区| 精品国产三级普通话版| 欧美中文日本在线观看视频| 日韩欧美 国产精品| 欧美日韩乱码在线| 日韩欧美在线乱码| 五月玫瑰六月丁香| 毛片女人毛片| 亚洲三级黄色毛片| 午夜精品在线福利| 欧美xxxx性猛交bbbb| 免费人成在线观看视频色| 国内少妇人妻偷人精品xxx网站| 99在线人妻在线中文字幕| 久久九九热精品免费| 亚洲国产精品久久男人天堂| 国产午夜精品论理片| 国产一级毛片七仙女欲春2| 久久久成人免费电影| 嫩草影院入口| 中文字幕人妻熟人妻熟丝袜美| 国内精品一区二区在线观看| 最新中文字幕久久久久| 午夜福利在线观看吧| 18+在线观看网站| 日日撸夜夜添| 欧美成人精品欧美一级黄| 国产极品精品免费视频能看的| 欧美日韩乱码在线| 看黄色毛片网站| 亚洲美女视频黄频| 久久99热6这里只有精品| 亚洲人成网站在线播放欧美日韩| 99国产精品一区二区蜜桃av| 亚洲无线观看免费| 国产一区二区在线av高清观看| 精品日产1卡2卡| 久久热精品热| 国产亚洲精品久久久com| 国产精品乱码一区二三区的特点| 97在线视频观看| 国产av一区在线观看免费| 最好的美女福利视频网| av女优亚洲男人天堂| 日韩,欧美,国产一区二区三区 | 久久久欧美国产精品| 久久这里只有精品中国| 久久久久性生活片| 国产午夜福利久久久久久| 久久久精品欧美日韩精品| 久久久久久久亚洲中文字幕| 美女大奶头视频| 久99久视频精品免费| 国产高清不卡午夜福利| av在线播放精品| 可以在线观看毛片的网站| 国内揄拍国产精品人妻在线| 日本黄色片子视频| 久久欧美精品欧美久久欧美| 禁无遮挡网站| 精品不卡国产一区二区三区| 亚洲欧美日韩卡通动漫| 特级一级黄色大片| 国产色爽女视频免费观看| 少妇丰满av| 少妇熟女欧美另类| 国产一级毛片七仙女欲春2| 99久久无色码亚洲精品果冻| 欧美+日韩+精品| 国产高清视频在线观看网站| 国产aⅴ精品一区二区三区波| 国产爱豆传媒在线观看| 免费搜索国产男女视频| 国产免费男女视频| 麻豆国产97在线/欧美| 一本精品99久久精品77| 九九久久精品国产亚洲av麻豆| 国产亚洲欧美98| 一进一出抽搐gif免费好疼| 亚洲精品国产av成人精品 | 99riav亚洲国产免费| 日韩,欧美,国产一区二区三区 | 亚洲精品粉嫩美女一区| 亚洲欧美日韩高清在线视频| 久久久久久久久久黄片| 日本一本二区三区精品| 国产伦精品一区二区三区视频9| 免费看日本二区| 人人妻,人人澡人人爽秒播| 欧美色视频一区免费| 你懂的网址亚洲精品在线观看 | 国产精品女同一区二区软件| 97超碰精品成人国产| 国产精品三级大全| 日本一本二区三区精品| 男女啪啪激烈高潮av片| 变态另类丝袜制服| 亚洲成人精品中文字幕电影| 淫秽高清视频在线观看| 国产麻豆成人av免费视频| 精品久久久久久久人妻蜜臀av| 自拍偷自拍亚洲精品老妇| 嫩草影院入口| 成人欧美大片| 精品久久久久久久久亚洲| 免费观看的影片在线观看| 久久亚洲国产成人精品v| 国产精品国产三级国产av玫瑰| 精品久久久久久成人av| 搡女人真爽免费视频火全软件 | 精品不卡国产一区二区三区| 99riav亚洲国产免费| 亚洲专区国产一区二区| videossex国产| 免费无遮挡裸体视频| 日本a在线网址| 日日摸夜夜添夜夜添av毛片| 免费观看人在逋| 91狼人影院| 淫秽高清视频在线观看| 久久久久免费精品人妻一区二区| 欧美一区二区亚洲| 国产免费男女视频| 亚洲av熟女| 欧美激情国产日韩精品一区| 日韩成人av中文字幕在线观看 | 最后的刺客免费高清国语| 亚洲经典国产精华液单| 天天躁日日操中文字幕| 国产白丝娇喘喷水9色精品| av在线播放精品| 黄色配什么色好看| 色播亚洲综合网| 伦理电影大哥的女人| 一区二区三区四区激情视频 | 午夜a级毛片| 99久国产av精品| 国产精品无大码| 免费看av在线观看网站| 18禁黄网站禁片免费观看直播| 老司机午夜福利在线观看视频| 2021天堂中文幕一二区在线观| 国产日本99.免费观看| 国产av在哪里看| 国产视频内射| 婷婷色综合大香蕉| 别揉我奶头 嗯啊视频| 国产伦精品一区二区三区视频9| 亚洲不卡免费看| 婷婷六月久久综合丁香| 久久久久精品国产欧美久久久| 亚洲综合色惰| 天堂网av新在线| 国产一区亚洲一区在线观看| 中文字幕久久专区| 九九热线精品视视频播放| 看黄色毛片网站| 国产一区二区三区在线臀色熟女| 欧美激情国产日韩精品一区| 少妇猛男粗大的猛烈进出视频 | 免费av毛片视频| 国内久久婷婷六月综合欲色啪| 99久久精品热视频| 99热这里只有是精品50| 国产一区二区亚洲精品在线观看| 国产亚洲精品综合一区在线观看| 日韩成人av中文字幕在线观看 | 国产精品嫩草影院av在线观看| av在线蜜桃| 国产精品一区www在线观看| 欧美一区二区精品小视频在线| 精品一区二区三区视频在线| 亚洲精品在线观看二区| 如何舔出高潮| 亚洲综合色惰| 国产探花在线观看一区二区| 国产一区二区三区在线臀色熟女| 不卡视频在线观看欧美| 久久韩国三级中文字幕| 女同久久另类99精品国产91| 久久九九热精品免费| 12—13女人毛片做爰片一| 国产黄色视频一区二区在线观看 | 午夜福利在线观看免费完整高清在 | 久久久精品94久久精品| 丝袜喷水一区| 99久久精品一区二区三区| 欧美+日韩+精品| 97超视频在线观看视频| 精品日产1卡2卡| 国产高清视频在线播放一区| 美女cb高潮喷水在线观看| 一个人看的www免费观看视频| 欧美一区二区亚洲| 午夜视频国产福利| 天堂√8在线中文| 男女那种视频在线观看| 亚洲综合色惰| 免费黄网站久久成人精品| 欧美日韩精品成人综合77777| 韩国av在线不卡| 国产精品综合久久久久久久免费| 免费观看的影片在线观看| 欧美激情久久久久久爽电影| 最近视频中文字幕2019在线8| 大又大粗又爽又黄少妇毛片口| 国产成人91sexporn| 精品无人区乱码1区二区| 亚洲第一区二区三区不卡| 91久久精品国产一区二区成人| 午夜视频国产福利| 国产欧美日韩精品亚洲av| 免费黄网站久久成人精品| 国产精品一区二区性色av| 综合色丁香网| 亚洲av成人精品一区久久| 日日啪夜夜撸| videossex国产| 日韩,欧美,国产一区二区三区 | 色视频www国产| 亚洲国产精品国产精品| 久久人人精品亚洲av| 麻豆成人午夜福利视频| 变态另类成人亚洲欧美熟女| 免费看美女性在线毛片视频| 国产成人a区在线观看| 男女下面进入的视频免费午夜| 国产69精品久久久久777片| 日韩人妻高清精品专区| 一本一本综合久久| 九九在线视频观看精品| 99国产精品一区二区蜜桃av| 亚洲,欧美,日韩| 国产一区二区三区在线臀色熟女| 精品国内亚洲2022精品成人| 我的女老师完整版在线观看| 亚洲国产精品久久男人天堂| 亚洲精品成人久久久久久| 国产日本99.免费观看| 在线观看av片永久免费下载| 欧美色视频一区免费| 成年女人看的毛片在线观看| 亚洲精品久久国产高清桃花| 色综合色国产| 69人妻影院| 日韩成人伦理影院| 亚洲精品粉嫩美女一区| 亚洲av成人精品一区久久| 91麻豆精品激情在线观看国产| 午夜激情欧美在线| 一进一出抽搐gif免费好疼| 最后的刺客免费高清国语| 国产精品一区二区三区四区免费观看 | 成人av在线播放网站| 亚洲中文字幕一区二区三区有码在线看| 黄色视频,在线免费观看| 啦啦啦观看免费观看视频高清| 丰满的人妻完整版| 一级av片app| 亚洲精品日韩av片在线观看| 久久婷婷人人爽人人干人人爱| 91精品国产九色| 国产爱豆传媒在线观看| 亚洲av免费高清在线观看| 亚洲av二区三区四区| 久久6这里有精品| 搡老熟女国产l中国老女人| 久久久久久国产a免费观看| 国产精品一区二区免费欧美| 色综合亚洲欧美另类图片| 九九爱精品视频在线观看| 3wmmmm亚洲av在线观看| 久久久久性生活片| 久久精品人妻少妇| 69av精品久久久久久| 麻豆国产av国片精品| 精品福利观看| 天美传媒精品一区二区| 欧美绝顶高潮抽搐喷水| 两个人视频免费观看高清| 赤兔流量卡办理| 国产一区二区亚洲精品在线观看| 久久久精品大字幕| 99久久精品热视频| av在线老鸭窝| 听说在线观看完整版免费高清| 免费无遮挡裸体视频| av女优亚洲男人天堂| 18禁黄网站禁片免费观看直播| 日本一本二区三区精品| 超碰av人人做人人爽久久| 全区人妻精品视频| 最新在线观看一区二区三区| 亚洲一级一片aⅴ在线观看| 午夜日韩欧美国产| 亚洲精品影视一区二区三区av| 日本免费一区二区三区高清不卡| 亚洲欧美精品综合久久99| 99热只有精品国产| 熟妇人妻久久中文字幕3abv| 午夜免费激情av| 日本与韩国留学比较| 国产aⅴ精品一区二区三区波| 最近视频中文字幕2019在线8| 三级男女做爰猛烈吃奶摸视频| 热99re8久久精品国产| 国产精品av视频在线免费观看| 看片在线看免费视频| 男女视频在线观看网站免费| 亚洲激情五月婷婷啪啪| 国产久久久一区二区三区| 久久久精品94久久精品| 一个人看的www免费观看视频| 亚洲成人久久爱视频| 国内揄拍国产精品人妻在线| 麻豆一二三区av精品| 国产黄色小视频在线观看| 永久网站在线| 99久国产av精品| 国产国拍精品亚洲av在线观看| 不卡视频在线观看欧美| 精品一区二区三区视频在线| 最新中文字幕久久久久| 亚洲人成网站在线观看播放| 97热精品久久久久久| 婷婷精品国产亚洲av在线| 精品久久久久久久久久久久久| 久久久精品94久久精品| 欧美区成人在线视频| 日本黄色片子视频| 国产高清视频在线观看网站| 黄色配什么色好看| 最近2019中文字幕mv第一页| 丰满的人妻完整版| 亚洲性久久影院| 男人和女人高潮做爰伦理| 少妇人妻一区二区三区视频| 大又大粗又爽又黄少妇毛片口| 久久6这里有精品| 国产乱人视频| 亚洲天堂国产精品一区在线| av在线天堂中文字幕| 全区人妻精品视频| 欧美成人精品欧美一级黄| 男女边吃奶边做爰视频| 日本成人三级电影网站| 久久久国产成人精品二区| 天堂√8在线中文| 亚洲色图av天堂| 韩国av在线不卡| 黑人高潮一二区| 欧美激情久久久久久爽电影| 三级毛片av免费| 我的女老师完整版在线观看| 日韩国内少妇激情av| 熟女人妻精品中文字幕| 久久久久久久久大av| 特级一级黄色大片| 乱系列少妇在线播放| 国产蜜桃级精品一区二区三区| 国产精品人妻久久久久久| 99久久精品国产国产毛片| 国产 一区精品| 国产伦精品一区二区三区四那| 六月丁香七月| 亚洲成人精品中文字幕电影| 一级毛片aaaaaa免费看小| 99久久无色码亚洲精品果冻| av中文乱码字幕在线| 国产精品一区二区免费欧美| 女人十人毛片免费观看3o分钟| 波多野结衣高清无吗| 99热这里只有是精品50| 久久韩国三级中文字幕| 欧美日本亚洲视频在线播放| 欧美一级a爱片免费观看看| 欧美成人a在线观看| 国产三级在线视频| 亚洲欧美日韩高清在线视频| 六月丁香七月| 国产 一区精品| 成人精品一区二区免费| 插阴视频在线观看视频| 精品久久久久久久久亚洲| 久久久欧美国产精品| 久久午夜亚洲精品久久| 草草在线视频免费看| 国产精品99久久久久久久久| 麻豆成人午夜福利视频| 亚洲av.av天堂| 久久精品国产亚洲av涩爱 | 九九爱精品视频在线观看|