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

    微型流化床內(nèi)松木屑和煤泥等溫混合熱解特性

    2017-10-16 08:13:32李天濤郭飛強(qiáng)王巖郭成龍董玉平
    化工學(xué)報(bào) 2017年10期
    關(guān)鍵詞:木屑煤泥流化床

    李天濤,郭飛強(qiáng),王巖,郭成龍,董玉平

    (1中國(guó)礦業(yè)大學(xué)電氣與動(dòng)力工程學(xué)院,江蘇 徐州 221116;2山東大學(xué)機(jī)械工程學(xué)院,山東 濟(jì)南 250061)

    微型流化床內(nèi)松木屑和煤泥等溫混合熱解特性

    李天濤1,郭飛強(qiáng)1,王巖1,郭成龍1,董玉平2

    (1中國(guó)礦業(yè)大學(xué)電氣與動(dòng)力工程學(xué)院,江蘇 徐州 221116;2山東大學(xué)機(jī)械工程學(xué)院,山東 濟(jì)南 250061)

    采用微型流化床反應(yīng)器對(duì)松木屑和煤泥的等溫混合熱解氣體釋放行為進(jìn)行實(shí)驗(yàn)研究。探討不同溫度和摻混比例對(duì)CH4、CO、CO2與H2釋放特性的影響,并通過(guò)模型配合法求解其動(dòng)力學(xué)參數(shù),研究松木屑和煤泥混合熱解過(guò)程的相互作用。通過(guò)FT-IR檢測(cè)發(fā)現(xiàn)煤泥的主要成分為含有C—O和CO鍵的芳香化合物,松木屑則以帶—OH鍵的長(zhǎng)鏈脂肪烴為主。在等溫穩(wěn)定反應(yīng)階段,松木屑熱解氣體生成速率高于煤泥,隨著生物質(zhì)摻入比例的不斷提高,混合原料氣體生成反應(yīng)速率亦呈現(xiàn)不同程度的增加。利用模型積分法求解了松木屑、煤泥及其混合物熱解氣體生成動(dòng)力學(xué)參數(shù),并通過(guò)實(shí)驗(yàn)值和計(jì)算值對(duì)比篩選出了最概然機(jī)理函數(shù)。通過(guò)活化能對(duì)比發(fā)現(xiàn),混合熱解對(duì)4種氣體組分生成具有不同的影響作用,其中CO實(shí)驗(yàn)活化能明顯低于計(jì)算值,表現(xiàn)為二者協(xié)同作用利于CO的生成釋放;對(duì) H2而言,在 75%混合比例條件下,混合反應(yīng)導(dǎo)致其生成活化能呈現(xiàn)協(xié)同負(fù)效應(yīng),使得活化能實(shí)驗(yàn)值明顯高于計(jì)算值;相較而言,CH4在混合熱解過(guò)程影響相對(duì)較弱,并呈現(xiàn)小幅度的協(xié)同負(fù)效應(yīng),而CO2的生成特性則受混合比例的影響較為明顯。

    微型流化床;生物質(zhì);煤泥;熱解;等溫反應(yīng);反應(yīng)動(dòng)力學(xué);相互作用

    Abstract:The releasing characteristics of gas species during isothermal co-pyrolysis of pine sawdust (PS) and coal slime (CS) were studied in a micro fluidized bed reactor.Four main pyrolysis gaseous products,CH4,CO,CO2and H2,were investigated under different temperature and blending ratios.The kinetic parameters were calculated employing the universal integral method to examine the interaction of PS and CS during co-pyrolysis.It was observed by FT-IR that CS mainly contained aromatic compounds carrying C—O and CO,while the PS was mainly consist of long-chain aliphatic hydrocarbon with —OH.The gas forming reaction rate of PS was higher than that of CS,and its value increased in a certain extent with the increase of biomass blending ratio.The kinetic parameters of pyrolysis gas formation of PS,CS and their blends were obtained by the mechanism function models,and the most probable mechanism function was determined by comparing the experimental and calculated values.The effect of blending PS and CS on the forming activations of the four gas species varied obviously.The activation energies of CO were significantly lower than the calculated ones,indicating that a positive synergisticeffect existed between PS and CS.The co-pyrolysis reaction led to a negative effect on the formation of H2at a PS blending ratio of 75% as the activation energy was obviously higher than the calculated one.In comparison,biomass blending showed relatively small positive synergistic effect on the formation of CH4,while that of CO2depended upon the blending ratio significantly.

    Key words:micro fluidized bed reactor; biomass,coal slime; pyrolysis; isothermal reaction; reaction kinetics;interaction

    引 言

    煤泥是洗煤過(guò)程中產(chǎn)生的一種副產(chǎn)品,顆粒較細(xì),水分與灰分含量較高,黏度較大,運(yùn)輸十分不便,利用不當(dāng)將造成能量的浪費(fèi)并存在潛在的環(huán)境污染風(fēng)險(xiǎn)[1-5]。目前常見的煤泥燃燒利用容易產(chǎn)生大量的細(xì)灰以及氮和硫的化合物,不利于環(huán)境的保護(hù)[6]。熱解氣化作為一種重要的潔凈煤利用技術(shù),可以實(shí)現(xiàn)煤泥在相對(duì)較低的溫度條件下轉(zhuǎn)化制備清潔燃?xì)?,但煤泥中氫元素含量較低,若在熱解氣化過(guò)程摻入富氫生物質(zhì)則能一定程度上對(duì)煤泥補(bǔ)充氫元素,改善氣體產(chǎn)物分布[7]。同時(shí),生物質(zhì)儲(chǔ)量巨大,來(lái)源廣泛,但密度和熱值較低,若與煤泥混合可提高整體能量密度,實(shí)現(xiàn)與煤泥的規(guī)?;旌限D(zhuǎn)化。再者,生物質(zhì)屬于二氧化碳零排放的新型能源,與煤泥混合熱解氣化轉(zhuǎn)化實(shí)現(xiàn)優(yōu)勢(shì)互補(bǔ),是新能源與傳統(tǒng)能源協(xié)同利用的有效途徑[8]。

    目前國(guó)內(nèi)外對(duì)生物質(zhì)和各類煤樣的混合熱解研究主要采用固定床反應(yīng)器、下降管反應(yīng)器和熱重分析儀等,尤其是熱重分析儀被廣泛用于混合熱解過(guò)程研究,是當(dāng)前研究煤和生物質(zhì)混合熱解過(guò)程相互作用的主要設(shè)備。部分研究結(jié)果證實(shí)煤與生物質(zhì)混合熱解過(guò)程可以降低焦炭產(chǎn)率,并改變氣體生成特性[9-10],視二者混合熱解存在相互作用或協(xié)同效應(yīng)。然而,部分研究結(jié)果認(rèn)為煤和生物質(zhì)混合熱解過(guò)程相互作用關(guān)系微弱,并不能說(shuō)明二者之間存在協(xié)同效應(yīng)[7,11-13]。由此可以看出,煤和生物質(zhì)混合熱解反應(yīng)過(guò)程復(fù)雜,受到多種因素的限制,包括原料種類、H/C比、升溫速率及無(wú)機(jī)金屬離子成分等。因此,研究二者之間混合熱解特性不僅要從整體轉(zhuǎn)化率的角度,還要注意主要產(chǎn)物組分生成的特性,并立足于二者處于同一反應(yīng)條件。

    近期,微型流化床反應(yīng)器研究氣固反應(yīng)動(dòng)力學(xué)特性被研究開發(fā),其可實(shí)現(xiàn)原料等溫快速流化反應(yīng),降低內(nèi)擴(kuò)散影響,在線快速檢測(cè)氣體產(chǎn)物釋放特性,進(jìn)而推測(cè)氣體生成動(dòng)力學(xué)特性[14-15]?;谀P头e分法,微型流化床分析儀已成功應(yīng)用于熱解反應(yīng)等研究中。Liu等[16]采用微型流化床對(duì)稻殼進(jìn)行催化熱解實(shí)驗(yàn),通過(guò)模型積分法求解其動(dòng)力學(xué)參數(shù),進(jìn)而研究了鉀離子對(duì)稻殼熱解的影響。Guo等[17]采用微型流化床研究了焦油催化裂解生成 H2和 CO的特性。Yu等[18]在微型流化床中進(jìn)行石墨粉燃燒反應(yīng),證實(shí)微型流化床是一種簡(jiǎn)單且可靠的動(dòng)力學(xué)研究?jī)x器。

    基于此,本文采用微型流化床分析儀研究煤泥與生物質(zhì)混合熱解動(dòng)力學(xué)特性,分析主要?dú)怏w組分的生成和轉(zhuǎn)化規(guī)律,并通過(guò)模型積分法求算動(dòng)力學(xué)參數(shù),進(jìn)而從氣體生成的角度闡明煤泥和生物質(zhì)混合熱解過(guò)程相互作用特性。

    1 實(shí)驗(yàn)材料和方法

    1.1 實(shí)驗(yàn)原料

    本文實(shí)驗(yàn)原料煤泥(CS)和松木屑(PS)取自徐州本地,其工業(yè)分析及元素分析如表1所示,煤泥有較高的固定碳和灰分含量,松木則是有較高的揮發(fā)分含量。兩種原料分別粉碎、篩分至150~250 μm,并按照生物質(zhì)質(zhì)量分?jǐn)?shù)分別為0%、25%、50%、75%、100% 5種比例進(jìn)行混合,在105℃的鼓風(fēng)干燥箱中干燥24 h,取出密封保存?zhèn)溆谩?/p>

    表1 原料工業(yè)分析及元素分析Table 1 Proximate and ultimate analysis of samples

    1.2 實(shí)驗(yàn)裝置與操作

    實(shí)驗(yàn)裝置和流程如圖1所示。

    實(shí)驗(yàn)在流化床石英反應(yīng)器上進(jìn)行,石英反應(yīng)器出口連接質(zhì)譜儀對(duì)出口氣體成分進(jìn)行連續(xù)監(jiān)測(cè)記錄。反應(yīng)器內(nèi)徑為20 mm,高150 mm,由兩層布風(fēng)板分成3段,兩層布風(fēng)板之間40 mm的區(qū)域?yàn)橹饕磻?yīng)區(qū)。詳盡的實(shí)驗(yàn)裝置介紹可以參考文獻(xiàn)[14-15]。每次實(shí)驗(yàn)前,將3 g石英砂(74~125 μm)放在石英反應(yīng)器底層作為流化介質(zhì)。前期研究表明[19],載氣流速在300 ml·min?1即可實(shí)現(xiàn)良好的流化性,并可以忽略外擴(kuò)散的影響,因而本實(shí)驗(yàn)氬氣的流量設(shè)定為 300 ml·min?1。首先,在流態(tài)化狀態(tài)下將石英反應(yīng)器加熱到設(shè)定溫度,然后用壓縮氣將10 mg樣品迅速注入石英反應(yīng)器,4種主要熱解產(chǎn)氣 H2、CO、CH4、CO2的釋放強(qiáng)度通過(guò)質(zhì)譜儀進(jìn)行連續(xù)檢測(cè)。為了保證實(shí)驗(yàn)結(jié)果的可靠性,每次實(shí)驗(yàn)重復(fù)3次。

    圖1 實(shí)驗(yàn)系統(tǒng)流程圖Fig.1 Schematic diagram of experimental system

    1.3 動(dòng)力學(xué)模型

    快速熱解實(shí)驗(yàn)中,各氣體的釋放特性由在線質(zhì)譜儀實(shí)時(shí)連續(xù)測(cè)得,產(chǎn)氣轉(zhuǎn)化率計(jì)算如下

    式中,t0,t和te分別代表開始反應(yīng)、反應(yīng)進(jìn)行到時(shí)間t和結(jié)束反應(yīng)的時(shí)間,s;分別代表t時(shí)刻和反應(yīng)結(jié)束時(shí)刻的氣體產(chǎn)量,g;φi代表氣體i的體積濃度分?jǐn)?shù),%;qv代表氣體的流量,L·min?1;MR表示氣體摩爾質(zhì)量;22.4 為標(biāo)準(zhǔn)狀態(tài)下氣體摩爾常數(shù);x代表氣體的轉(zhuǎn)化率;dx/dt代表氣體轉(zhuǎn)化速率,s?1。

    其中,f(x)代表反應(yīng)的微分機(jī)理函數(shù);k(T)代表Arrhenius速率常數(shù),由Arrhenius方程所定義,在等溫過(guò)程中是一個(gè)常量。

    本文采用模型積分法對(duì)各氣體生成過(guò)程進(jìn)行動(dòng)力學(xué)求解,計(jì)算方法如下

    表2 用于氣固反應(yīng)的典型機(jī)理函數(shù)Table 2 Typical mechanism functions using in gas-solid reactions

    式中,G(x)是機(jī)理函數(shù)的積分形式;Ea代表表觀活化能,kJ·mol?1;A表示指前因子,s?1;T代表溫度,K;R代表氣體常數(shù),8.314 J·mol?1·K?1。由式(7)可以擬合出一條機(jī)理函數(shù)G(x)關(guān)于時(shí)間t的直線,直線斜率即為k(T)。再根據(jù)式(7),擬合得到速率常數(shù)k(T)與溫度T的直線,即可根據(jù)直線的斜率與截距求得表觀活化能Ea和指前因子A。將求得的動(dòng)力學(xué)參數(shù)進(jìn)行逆向運(yùn)算可以獲得轉(zhuǎn)化率計(jì)算值,進(jìn)而通過(guò)對(duì)實(shí)驗(yàn)值和計(jì)算值的誤差分析,篩選最概然機(jī)理函數(shù)。對(duì)于氣相與固相之間的化學(xué)反應(yīng),常見的非勻相反應(yīng)模型有20種,如表2所示[20],這些函數(shù)已廣泛應(yīng)用于描述熱解反應(yīng)過(guò)程[18,20-23]。

    2 實(shí)驗(yàn)結(jié)果與討論

    2.1 生物質(zhì)和煤的FT-IR分析

    采用傅里葉變換紅外光譜儀(VERTEX 80v)分析松木屑和煤泥的結(jié)構(gòu)特征,結(jié)果如圖2所示。煤泥和松木屑在結(jié)構(gòu)上既有相似之處也有明顯不同,第一個(gè)主要區(qū)別出現(xiàn)在波數(shù) 3568~3730 cm?1范圍內(nèi),煤泥有3個(gè)連續(xù)峰峰值分別在3622 cm?1、3658 cm?1和3695 cm?1,表示煤泥中帶有小分子側(cè)鏈的亞氨基N—H結(jié)構(gòu),而松木屑在這一范圍內(nèi)沒有波峰。吸收峰 3000~3568 cm?1主要是由羥基的振動(dòng)引起,松木屑在此范圍出現(xiàn)一個(gè)較強(qiáng)的吸收峰,而煤泥在此區(qū)域的峰值較弱,其相應(yīng)基團(tuán)含量明顯低于松木屑。在 2765~3000 cm?1波數(shù)范圍內(nèi),煤泥分別在2854 cm?1和2919 cm?1出現(xiàn)兩個(gè)微弱吸收峰,松木屑則只有一個(gè)強(qiáng)峰在2903 cm?1處,表明煤泥中存在少量環(huán)烷或脂肪族的亞甲基—CH2和次甲基—CH,而松木屑中有較多的亞甲基—CH2,即松木屑中有較長(zhǎng)的脂肪鏈。由此可以看出,松木屑中含氫基團(tuán)明顯多于煤泥,在熱解過(guò)程中可以作為主要的氫源,釋放含氫氣體組分。1378~1802 cm?1范圍的吸收峰主要是由芳香結(jié)構(gòu)CC和CO伸縮振動(dòng)產(chǎn)生的,煤泥有1個(gè)主峰和1個(gè)肩峰,峰值分別在1442 cm?1和1602 cm?1處,而松木屑在這一波數(shù)范圍內(nèi)波峰多且雜,表明松木屑中有著豐富的含氧基團(tuán)。在833~1378 cm?1范圍內(nèi),煤泥出現(xiàn)了較強(qiáng)的吸收峰,在 1032 cm?1處,主要為有機(jī)質(zhì)中醇和烷基醚的C—O鍵,可能會(huì)導(dǎo)致較高的CO和CO2產(chǎn)率;松木屑中波峰仍然繁雜,主峰出現(xiàn)在1037 cm?1處,對(duì)應(yīng)于醇和烷基醚的C—O鍵和羧酸中的醚鍵。由此可見,松木屑中有豐富的脂肪族結(jié)構(gòu),且脂肪鏈較長(zhǎng),羧酸基團(tuán)及醚鍵的連接方式多樣,而煤泥中更多的是芳香族結(jié)構(gòu),相較而言更為穩(wěn)定,熱解溫度更高[24]。反應(yīng)過(guò)程中若兩種基團(tuán)發(fā)生交互作用實(shí)現(xiàn)元素互補(bǔ),可以促進(jìn)二者熱解反應(yīng)。

    圖2 松木屑和煤泥的紅外光譜圖Fig.2 FT-IR spectra of pine sawdust and coal slime samples

    2.2 松木屑和煤泥單獨(dú)熱解特性

    松木屑和煤泥在注射入石英反應(yīng)器后迅速進(jìn)入熱解反應(yīng)階段。實(shí)驗(yàn)由在線質(zhì)譜儀實(shí)時(shí)連續(xù)檢測(cè)樣本在750~900℃熱解條件下的 H2、CO、CO2和CH4釋放強(qiáng)度。圖3展示了松木屑和煤泥在750℃和 900℃下的氣體釋放強(qiáng)度隨時(shí)間的變化。可以看到在微型流化床等溫條件下,松木屑和煤泥產(chǎn)氣熱解釋放基本都在反應(yīng)開始后的15~20 s內(nèi)完成且不同氣體釋放的起止時(shí)間并不相同,不同溫度條件下各氣體的釋放強(qiáng)度也有明顯差異。

    在750~900℃熱解條件下,比較煤泥和松木屑熱解氣體釋放強(qiáng)度可以發(fā)現(xiàn),松木屑的CO釋放強(qiáng)度約為煤泥的3倍。同時(shí),松木屑熱解釋放CO2、H2和CH4的釋放強(qiáng)度同樣高于煤泥,這也與常規(guī)煤與生物質(zhì)熱解特性的預(yù)期相符合。生物質(zhì)的主要組分包括纖維素,半纖維素和木質(zhì)素,相互之間以醚鍵(R—O—R)連接,醚鍵的鍵能為 380~420 kJ·mol?1[25],熱解條件下能夠釋放大量的氣體。相比之下,煤泥中固定碳和灰分含量較高,有機(jī)物主要組成為稠多環(huán)芳烴,連接鍵為芳香環(huán)鍵(CC),鍵能達(dá)到1000 kJ·mol?1[26],熱解溫度較高。另外,在相同溫度條件下生物質(zhì)含有更多的氫元素,使得含氫氣體組分(H2和CH4)產(chǎn)率更高,因而如果混合熱解,可為煤泥的熱解提供氫源,使煤熱解產(chǎn)生的大量不穩(wěn)定分子自由基得到配對(duì)電子穩(wěn)定成為低分子揮發(fā)性物質(zhì)[27]。

    在 750℃熱解條件下,松木屑和煤泥最先釋放的組分均為CO2,其次是CO、CH4和H2。隨著溫度的升高,氣體釋放順序發(fā)生變化。900℃下H2的釋放明顯提前,H2的釋放強(qiáng)度也明顯增強(qiáng)。另外CH4和CO的釋放強(qiáng)度也隨著溫度的升高而有所增強(qiáng),而CO2的釋放強(qiáng)度沒有太大變化。各氣體不同的釋放行為意味著其不同的生成路徑和反應(yīng)機(jī)理。在生物質(zhì)熱解中,CO通常由醚鍵(C—O—C)和羰基(CO)的斷裂和揮發(fā)分的二次裂解產(chǎn)生,在高溫下更容易發(fā)生[28]。H2通常由自由基縮聚反應(yīng)和脫氫反應(yīng)產(chǎn)生,CH4的形成通常伴隨—OCH3和—CH2—的斷裂[29],同樣與溫度密切相關(guān)。相較之下,CO2主要由 CO鍵和羧基(—COOH)的裂解和重整產(chǎn)生[28],鍵能相對(duì)較弱,因此CO2的產(chǎn)生受熱解溫度的影響較小。

    圖3 熱解氣體釋放特性Fig.3 Gas releasing characteristics during pyrolysis process

    圖4 不同溫度下的各氣體轉(zhuǎn)化率隨時(shí)間的變化Fig.Conversion profiles of gas species related to time and temperature

    圖4為不同溫度條件下4種氣體相對(duì)生成率隨時(shí)間的變化。對(duì)于松木屑和煤泥,隨著溫度的升高,4種氣體完成轉(zhuǎn)化的時(shí)間都相應(yīng)變短,表明產(chǎn)氣反應(yīng)與溫度高度相關(guān)。除此之外,不同氣體完成轉(zhuǎn)化的時(shí)間也有明顯差別,同樣暗示著其不同的化學(xué)反應(yīng)生成路徑?;谑?5),將圖4中x對(duì)t求微分,可以進(jìn)一步得到轉(zhuǎn)化速率dx/dt與轉(zhuǎn)化率x之間的關(guān)系,如圖5所示。從圖5可以看出在同一轉(zhuǎn)化率下,各氣體的轉(zhuǎn)化速率都隨著溫度的升高而加快。在反應(yīng)初始階段,原料快速升溫并迅速釋放氣體,使得氣體轉(zhuǎn)化速率迅速增長(zhǎng)并在轉(zhuǎn)化率0.25附近達(dá)到最大值。隨后,氣體轉(zhuǎn)化速率隨著轉(zhuǎn)化率的增加逐漸降低,表明熱解反應(yīng)已經(jīng)進(jìn)行到等溫?zé)峤怆A段。反應(yīng)進(jìn)行到后期,氣體組分的生成主要依賴于芳香核的縮合反應(yīng),導(dǎo)致反應(yīng)速率逐漸下降[30]。

    圖5 松木屑和煤泥熱解氣生成反應(yīng)速率隨轉(zhuǎn)化率的變化Fig.5 Reaction rate of gas releasingversusconversion of CS and PS

    圖6 混合樣本產(chǎn)氣釋放速率與轉(zhuǎn)化率的關(guān)系Fig.6 Reaction rate of gas releasingversusconversion of blends

    2.3 混合物的混合熱解特性

    當(dāng)在煤泥中添加松木屑后,以850℃為例,不同混合比例的混合樣品在熱解過(guò)程中各氣體轉(zhuǎn)化速率隨轉(zhuǎn)化率的變化如圖6所示。整體而言,4種氣體在穩(wěn)定的等溫?zé)峤夥磻?yīng)區(qū)(x>0.25),隨著添加松木屑添加比例的提高,反應(yīng)速率均呈現(xiàn)不同程度的升高,表明松木屑熱解反應(yīng)活性高于煤泥。但松木屑的添加對(duì)4種氣體的影響也存在明顯區(qū)別,其中H2受生物質(zhì)摻入比例的影響最為明顯,反應(yīng)速率隨松木屑比例增加呈現(xiàn)明顯增加的趨勢(shì),當(dāng)木屑摻入比例達(dá)到75%時(shí),H2的釋放反應(yīng)速率接近于純松木屑。對(duì)CO而言,松木屑的混合明顯增加了CO的轉(zhuǎn)化速率,特別是當(dāng)生物質(zhì)摻入比例為50%時(shí),其反應(yīng)速率高于75%比例的混合樣品。然而煤泥和松木屑熱解生成CH4和CO2的反應(yīng)速率較為接近,混合熱解對(duì)于兩種氣體轉(zhuǎn)化速率影響較為微弱。由此可見,混合原料熱解過(guò)程中,氣體生成速率一定程度上取決于煤泥和松木屑自身的反應(yīng)速率,另外在部分混合比例條件下反應(yīng)速率出現(xiàn)一定的變化,可能出現(xiàn)了一定的相互作用。

    2.4 模型配合法求解動(dòng)力學(xué)參數(shù)

    本文采用模型積分法對(duì)各氣體進(jìn)行動(dòng)力學(xué)求解。將各機(jī)理模型函數(shù)分別代入式(8)進(jìn)行擬合,得到線性擬合程度較好的機(jī)理模型函數(shù)為G3(三維球形擴(kuò)散模型)、G7(1級(jí)成核生長(zhǎng)模型)和G17(1.5級(jí)化學(xué)反應(yīng)模型)。在松木屑與煤泥 1:1混合比例下,機(jī)理模型函數(shù)G3、G7、G17與反應(yīng)時(shí)間t的關(guān)系如圖7所示,其擬合相關(guān)系數(shù)列于表3,4種氣體相關(guān)系數(shù)均在0.95以上,其他比例下的擬合結(jié)果與之類似,表明3種模型能夠較好地描述4種氣體在微型流化床中的等溫?zé)峤忉尫盘匦浴?種機(jī)理函數(shù)模型獲得的動(dòng)力學(xué)參數(shù)如表4和表5所示。為篩選出最概然機(jī)理函數(shù),本文將模型配合法計(jì)算得到的表觀活化能依次代入式(8)和式(7)求算轉(zhuǎn)化率的計(jì)算值,并與實(shí)驗(yàn)值進(jìn)行比較,采用式(9)計(jì)算實(shí)驗(yàn)值與計(jì)算值的平均相對(duì)誤差[31],式中定義 OF為平均相對(duì)誤差值,OF最小的機(jī)理函數(shù)即被認(rèn)定為最概然機(jī)理函數(shù)。

    圖7 煤泥與松木屑1:1混合樣本的機(jī)理模型函數(shù)Fig.7 Linear goodness of fit mechanism function model determination for 50% blending ratio of PS

    以松木屑與煤泥1:1摻混為例比較實(shí)驗(yàn)值與計(jì)算值的差異如圖8所示。計(jì)算結(jié)果顯示G3模型擬合的OF值范圍為0.03~0.07,G7的OF值集中在0.03~0.05,G17的OF值大部分在0.04~0.08之間,G7可以實(shí)現(xiàn)對(duì)于氣體轉(zhuǎn)化更好的計(jì)算預(yù)測(cè)。由表5可以看出,通過(guò)G7機(jī)理模型函數(shù)進(jìn)行計(jì)算,煤泥熱解生成H2、CO、CO2和CH4的活化能分別為75.0、50.3、37.6、50.8 kJ·mol?1;相應(yīng)的松木屑 4 種氣體生成活化能分別為 70.3、37.9、54.9、68.0 kJ·mol?1。Guo等[15]和Yu等[14]采用類似微型流化床分別研究中藥渣和酒糟熱解動(dòng)力學(xué)特性,獲得的4種氣體熱解活化能數(shù)值與本實(shí)驗(yàn)結(jié)果相當(dāng)。煤泥熱解反應(yīng)過(guò)程 CO2生成活化能最低,H2最高,說(shuō)明熱解反應(yīng)中 CO2最易生成,H2釋放最為困難;而松木屑熱解反應(yīng) CO生成活化能最低,H2最高,因此所獲得的活化能值在一定程度上也與氣體的釋放順序相一致。

    表3 煤泥松木屑1:1混合樣本的機(jī)理模型函數(shù)的線性相關(guān)系數(shù)Table 3 Linear correlation coefficient of mechanism function model determination for 50% blending ratio of PS

    表4 不同混合比例下各氣體指前因子Table 4 Pre-exponential factor of gas species related to blend ratios in different mechanism functions

    表5 不同混合比例下各氣體表觀活化能Table 5 Apparent activation energy of gas species related to blend ratios in different mechanism functions

    2.5 煤泥與生物質(zhì)相互作用特性分析

    在混合熱解反應(yīng)中,煤泥和生物質(zhì)的相互作用可以通過(guò)預(yù)測(cè)其混合樣本的活化能進(jìn)行討論。Goldfarb等[32]認(rèn)為,如果混合樣本在反應(yīng)動(dòng)力學(xué)方面沒有協(xié)同效應(yīng),則混合樣本的活化能應(yīng)為獨(dú)立樣本活化能的加權(quán)平均值?;旌蠘颖镜挠?jì)算活化能可由式(10)得到

    其中,EC代表計(jì)算活化能,kJ·mol?1;EPS和ECS分別表示松木屑和煤泥的實(shí)驗(yàn)活化能;kJ·mol?1;vPS代表松木屑的混合質(zhì)量比,%。

    圖9為采用G7模型計(jì)算得到的各氣體生成活化能的實(shí)驗(yàn)值和計(jì)算值對(duì)比。可以看出混合樣本與單獨(dú)樣本之間并不遵從線性關(guān)系,4種產(chǎn)氣中,混合樣本的H2和CH4實(shí)驗(yàn)活化能較計(jì)算活化能偏高,CO的實(shí)驗(yàn)活化能較計(jì)算活化能偏低。

    為了進(jìn)一步表征松木屑和煤泥間的相互作用,定義計(jì)算活化能與實(shí)驗(yàn)活化能的相對(duì)偏差δ如式(11)所示

    式中,Emeasured代表實(shí)驗(yàn)活化能,kJ·mol?1。不同比例下松木屑與煤泥混合樣本活化能相對(duì)偏差如圖10所示。

    對(duì)于CH4而言,在松木屑混合比例為25%、50%和75%下,相對(duì)偏差均為正值,分別為8.0%、8.6%、7.4%,表明實(shí)驗(yàn)活化能偏高,兩者的混合熱解對(duì)CH4的產(chǎn)生具有一定的不利作用。對(duì)于CO而言,當(dāng)松木屑混合比例為 25%、50%、75%時(shí),活化能實(shí)驗(yàn)值和計(jì)算值呈現(xiàn)明顯偏差,分別為?26.0%、?30.7%、?29.9%,說(shuō)明混合原料熱解過(guò)程CO釋放所需能量壘勢(shì)降低,對(duì)CO的生成有明顯的促進(jìn)作用。在CO2生成過(guò)程,在松木屑與煤泥比例為1:3時(shí),相對(duì)偏差幾乎為 0,兩者的混合熱解沒有協(xié)同效應(yīng)發(fā)生;當(dāng)松木屑與煤泥比例為1:1時(shí),相關(guān)系數(shù)為7.6%,混合熱解不利于CO2的生成;當(dāng)松木屑摻混比例提高到75%時(shí),相關(guān)系數(shù)降到?10.6%,反而促進(jìn)了 CO2的產(chǎn)生。對(duì) H2而言,當(dāng)松木屑混合比例為 50%和 75%時(shí),相關(guān)系數(shù)分別為 10.1%、20.5%,說(shuō)明在較高的松木屑混合比例條件下,H2生成所需能量增加,可能由于在生物質(zhì)比例過(guò)高時(shí),低密度的生物質(zhì)顆粒在熱解中軟化、變形,包裹黏附在煤泥顆粒的表面,堵塞煤泥的孔隙,抑止煤泥中揮發(fā)分的逸出和擴(kuò)散[33]。通過(guò)活化能的變化可以發(fā)現(xiàn),混合熱解過(guò)程對(duì)于不同氣體的活化能影響不同,整體而言混合熱解對(duì)于降低CO生成活化能影響較為明顯,表現(xiàn)為CO生成過(guò)程發(fā)生協(xié)同效應(yīng),降低了其生成的表觀活化能;在高松木屑摻入比例條件下H2的活化能高于理論計(jì)算值,表現(xiàn)出一定的負(fù)協(xié)同作用。

    圖8 煤泥與松木屑1:1混合樣本的機(jī)理模型函數(shù)值與實(shí)驗(yàn)值比較Fig.8 Comparison between mechanism function model value and experimental value of 50% blending ratio of PS

    圖9 G7模型下不同生物質(zhì)比例的各氣體表觀活化能Fig.9 Apparent activation energy of gas species in different biomass blending ratio in G7

    3 結(jié) 論

    (1)煤泥和松木屑在官能團(tuán)結(jié)構(gòu)存在明顯差異,主要體現(xiàn)在含氫官能團(tuán)和含氧官能團(tuán),在熱解過(guò)程中有利于官能團(tuán)相互作用和元素互補(bǔ)。

    (2)松木屑和煤泥等溫?zé)峤膺^(guò)程,松木屑?xì)怏w釋放反應(yīng)速率高于煤泥,混合比例的增加一定程度可以提高氣體生成速率。

    (3)松木屑和煤泥的混合熱解過(guò)程相互作用對(duì)于不同氣體影響存在差異,其中對(duì)于CO生成影響最為明顯,混合反應(yīng)降低了CO生成的活化能;在高生物質(zhì)摻入比例條件下,二者相互作用導(dǎo)致 H2活化能升高,表現(xiàn)出協(xié)同負(fù)效應(yīng);混合反應(yīng)對(duì)于CH4和CO2生成影響相對(duì)較弱,表現(xiàn)為混合反應(yīng)會(huì)導(dǎo)致CH4活化能小幅度增加,而CO2活化能變化受原料混合比例影響明顯。

    [1]SONG Z L,JING C M,YAO L S,et al.Coal slime hot air/microwave combined drying characteristics and energy analysis[J].Fuel Processing Technology,2017,156: 491-499.

    [2]ZHOU K,LIN Q Z,HU H W,et al.The ignition characteristics and combustion processes of the single coal slime particle under different hot-coflow conditions in N2/O2atmosphere[J].Energy,2017,136:173-184.

    [3]張坤.火電廠煤泥型煤制備及工程應(yīng)用方案初步設(shè)計(jì)[D].北京:華北電力大學(xué),2013.ZHANG K.Preliminary design on preparation and engineering applications of slime briquette in thermal power plant[D].Beijing:North China Electric Power University,2013.

    [4]劉明強(qiáng),劉建忠,王睿坤,等.煤泥燃燒動(dòng)力學(xué)機(jī)理研究的新方法探討[J].煤炭學(xué)報(bào),2012,37(a02): 444-448.LIU M Q,LIU J Z,WANG R K,et al.Discussion of a new approach to study kinetic mechanism of coal slurry combustion[J].Journal of China Coal Society,2012,37(a02): 444-448.

    [5]徐曉虹,李法香,吳建鋒,等.利用煤泥制備蓄熱陶瓷及抗熱震性能[J].武漢理工大學(xué)學(xué)報(bào),2016,36(11): 6-10.XU X H,LI F X,WU J F,et al.Preparation and thermal shock resistance of thermal storage ceramic from coal slime[J].Journal of Wuhan University of Technology,2016,36(11): 6-10.

    [6]陳恒寶,周敏,賀國(guó)章.煤泥的二氧化碳?xì)饣瘎?dòng)力學(xué)[J].煤炭轉(zhuǎn)化,2012,35(3): 19-22.CHEN H B,ZHOU M,HE G Z.Carbon dioxide gasification kinetics of coal slime[J].Coal Conversion,2012,35(3): 19-22.

    [7]WEILAND N T,MEANS N C,MORREALE B D.Product distributions from isothermal co-pyrolysis of coal and biomass[J].Fuel,2011,94(1): 563-570.

    [8]WU Z Q,WANG S Z,ZHAO J,et al.Product distribution during co-pyrolysis of bituminous coal and lignocellulosic biomass major components in a drop-tube furnace[J].Energy & Fuels,2015,29(7):4168-4180.

    [9]ZHANG L,XU S P,ZHAO W,et al.Co-pyrolysis of biomass and coal in a free fall reactor[J].Fuel,2007,86(3): 353-359.

    [10]DU Y Y,JIANG X G,LV G,et al.Thermal behavior and kinetics of bio-ferment residue/coal blends during co-pyrolysis[J].Energy Conversion & Management,2014,88: 459-463.

    [11]GIL M V,CASAL D,PEVIDA C,et al.Thermal behaviour and kinetics of coal/biomass blends during co-combustion[J].Bioresource Technology,2010,101(14): 5601-5608.

    [12]ZHOU C C,LIU G J,FANG T,et al.Investigation on thermal and trace element characteristics during co-combustion biomass with coal gangue.[J].Bioresource Technology,2015,175(175C):454-462.

    [13]LU K M,LEE W J,CHEN W H,et al.Thermogravimetric analysis and kinetics of co-pyrolysis of raw/torrefied wood and coal blends[J].Applied Energy,2013,105(2): 57-65.

    [14]YU J,YAO C B,ZENG X,et al.Biomass pyrolysis in a micro-fluidized bed reactor: characterization and kinetics[J].Chemical Engineering Journal,2011,168(2): 839-847.

    [15]GUO F Q,DONG Y P,LV Z C,et al.Pyrolysis kinetics of biomass(herb residue) under isothermal condition in a micro fluidized bed[J].Energy Conversion & Management,2015,93: 367-376.

    [16]LIU Y,WANG Y,GUO F Q,et al.Characterization of the gas releasing behaviors of catalytic pyrolysis of rice husk using potassium over a micro-fluidized bed reactor[J].Energy Conversion and Management,2017,136: 395-403.

    [17]GUO F Q,DONG Y P,FAN P F,et al.Catalytic decomposition of biomass tar compound by calcined coal gangue: a kinetic study[J].International Journal of Hydrogen Energy,2016,41(31):13380-13389.

    [18]YU J,ZENG X,ZHANG J,et al.Isothermal differential characteristics of gas-solid reaction in micro-fluidized bed reactor[J].Fuel,2013,103(1): 29-36.

    [19]GUO F Q,DONG Y P,FAN P F,et al.Detailed kinetic study of phenol decomposition under isothermal conditions to understand tar catalytic cracking process[J].Journal of Analytical & Applied Pyrolysis,2016,118: 155-163.

    [20]MAO Y B,DONG L,DONG Y P,et al.Fast co-pyrolysis of biomass and lignite in a micro fluidized bed reactor analyzer[J].Bioresource Technology,2015,181: 155-162.

    [21]朱劍虹.微型流化床反應(yīng)分析儀的應(yīng)用[D].湘潭: 湘潭大學(xué),2010.ZHU J H.Application of micro fluidized bed reaction analyzer[D].Xiangtan: Xiangtan University,2010.

    [22]余劍,李強(qiáng),段正康,等.微型流化床中的等溫微分反應(yīng)特性[J].中國(guó)科學(xué):化學(xué),2011,41(1): 152-160.YU J,LI Q,DUAN Z K,et al.Isothermal differential characteristics of the reaction in micro fluidized bed[J].Scientia Sinica Chimica,2011,41(1):152-162.

    [23]GAI C,DONG Y P,LV Z C,et al.Pyrolysis behavior and kinetic study of phenol as tar model compound in micro fluidized bed reactor[J].International Journal of Hydrogen Energy,2015,40(25):7956-7964.

    [24]車德勇.基于物化結(jié)構(gòu)特征的生物質(zhì)與煤共氣化特性研究[D].北京: 華北電力大學(xué),2013.CHE D Y.Research on co-gasification mechanism of biomass and coal based on surface physicochemical structure characteristics[D].Beijing: North China Electric Power University,2013.

    [25]BLA?EJ A,KO?íK M.Phytomass : A Raw Material for Chemistry and Biotechnology[M].Ellis Horwood,1993.

    [26]MOGHTADERI B,MEESRI C,WALL T F.Pyrolytic characteristics of blended coal and woody biomass[J].Fuel,2004,83(6): 745-750.

    [27]王鵬,文芳,邊文,等.煤與生物質(zhì)共熱解特性初步研究[J].煤炭轉(zhuǎn)化,2008,31(4): 40-44.WANG P,WEN F,BIAN W,et al.Study on the co-pyrolysis characteristics of coal and biomass[J].Coal Conversion,2008,31(4):40-44.

    [28]YANG H P,YAN R,CHEN H P,et al.Characteristics of hemicellulose,cellulose and lignin pyrolysis[J].Fuel,2007,86(12/13):1781-1788.

    [29]LIU Q,WANG S R,ZHENG Y,et al.Mechanism study of wood lignin pyrolysis by using TG-FTIR analysis[J].Journal of Analytical& Applied Pyrolysis,2008,82(1): 170-177

    [30]STANIS?AW P.The reactions of formation of selected gas products during coal pyrolysis[J].Fuel,2004,83(9): 1191-1196.

    [31]CAPART R,KHEZAMI L,BURNHAM A K.Assessment of various kinetic models for the pyrolysis of a microgranular cellulose[J].Thermochimica Acta,2004,417(1): 79-89.

    [32]GOLDFARB J L,CEYLAN S.Second-generation sustainability:application of the distributed activation energy model to the pyrolysis of locally sourced biomass-coal blends for use in co-firing scenarios[J].Fuel,2015,160: 297-308.

    [33]閻維平,陳吟穎.生物質(zhì)混合物與褐煤共熱解特性的試驗(yàn)研究[J].動(dòng)力工程學(xué)報(bào),2006,26(6): 865-870.YAN W P,CHEN Y Y.Experimental study on co-pyrolysis characteristics of lignite mixed with biomass mixture[J].Journal of Power Engineering,2006,26(6): 865-870.

    Characterization of co-pyrolysis of pine sawdust and coal slime under isothermal conditions in micro fluidized bed reactor

    LI Tiantao1,GUO Feiqiang1,WANG Yan1,GUO Chenglong1,DONG Yuping2
    (1School of Electrical and Power Engineering,China University of Mining and Technology,Xuzhou221116,Jiangsu,China;2School of Mechanical Engineering,Shandong University,Jinan250061,Shandong,China)

    X 752

    A

    0438—1157(2017)10—3923—11

    10.11949/j.issn.0438-1157.20170393

    2017-04-13收到初稿,2017-05-25收到修改稿。

    聯(lián)系人:郭飛強(qiáng)。

    李天濤(1992—),男,碩士研究生。

    國(guó)家自然科學(xué)基金項(xiàng)目(51406226)。

    Received date:2017-04-13.

    Corresponding author:GUO Feiqiang,fqguo@cumt.edu.cn

    Foundation item:supported by the National Natural Science Foundation of China (51406226).

    猜你喜歡
    木屑煤泥流化床
    木屑和鐵屑的分離實(shí)驗(yàn)
    煤泥重介質(zhì)旋流器在動(dòng)力煤選煤廠粗煤泥分選中的應(yīng)用
    MG-9#捕收劑在極細(xì)粒煤泥浮選中的應(yīng)用
    煤泥水中煤與不同礦物相互作用的模擬研究
    木屑
    飛天(2020年4期)2020-04-17 14:48:24
    流化床丙烷脫氫反應(yīng)段的模擬及優(yōu)化
    關(guān)于循環(huán)流化床鍋爐集控運(yùn)行研究
    煤泥分選超凈煤的藥劑研究
    油頁(yè)巖與木屑混合熱解特性研究
    單沉浸管流化床內(nèi)離散顆粒數(shù)值模擬
    2021天堂中文幕一二区在线观| 亚洲国产成人一精品久久久| 我的女老师完整版在线观看| 成人亚洲欧美一区二区av| 日本三级黄在线观看| 美女cb高潮喷水在线观看| 床上黄色一级片| 成人毛片a级毛片在线播放| 欧美高清成人免费视频www| 建设人人有责人人尽责人人享有的 | 亚洲人成网站高清观看| 午夜老司机福利剧场| 国产成年人精品一区二区| 亚洲av男天堂| 国产一区二区亚洲精品在线观看| 日韩强制内射视频| 毛片一级片免费看久久久久| 亚洲婷婷狠狠爱综合网| 18禁裸乳无遮挡免费网站照片| 禁无遮挡网站| 亚洲精品日本国产第一区| 欧美xxxx性猛交bbbb| 午夜免费激情av| 亚洲av免费在线观看| 亚洲熟女精品中文字幕| 91aial.com中文字幕在线观看| 久久99精品国语久久久| freevideosex欧美| 国产成人精品久久久久久| 大又大粗又爽又黄少妇毛片口| 97精品久久久久久久久久精品| 国产片特级美女逼逼视频| 国产麻豆成人av免费视频| 亚洲精品成人av观看孕妇| 国产91av在线免费观看| 国产高清三级在线| 久久这里有精品视频免费| 天天一区二区日本电影三级| 国产在视频线在精品| 国产精品久久久久久精品电影小说 | 美女xxoo啪啪120秒动态图| 亚洲国产欧美在线一区| 美女黄网站色视频| 一级a做视频免费观看| 综合色丁香网| 成人无遮挡网站| 99久国产av精品国产电影| 亚洲18禁久久av| 国产精品99久久久久久久久| 成人综合一区亚洲| 99久久精品一区二区三区| 赤兔流量卡办理| 国产精品伦人一区二区| 一夜夜www| 国产亚洲最大av| 亚洲欧美精品专区久久| 最近中文字幕2019免费版| 久久精品熟女亚洲av麻豆精品 | 亚洲三级黄色毛片| 亚洲天堂国产精品一区在线| 一级爰片在线观看| 日韩欧美一区视频在线观看 | 亚洲成人久久爱视频| 日日摸夜夜添夜夜爱| 麻豆国产97在线/欧美| 亚洲av成人精品一二三区| 日韩精品青青久久久久久| 人人妻人人澡人人爽人人夜夜 | 人体艺术视频欧美日本| 尾随美女入室| 黄色欧美视频在线观看| 蜜桃亚洲精品一区二区三区| 久久久色成人| 免费观看无遮挡的男女| 亚洲欧洲日产国产| 亚洲伊人久久精品综合| 国内精品美女久久久久久| 女人久久www免费人成看片| 人人妻人人看人人澡| 久久精品国产鲁丝片午夜精品| 99热6这里只有精品| 亚洲欧美清纯卡通| 国产高清有码在线观看视频| 国产av不卡久久| 亚洲怡红院男人天堂| 国产精品综合久久久久久久免费| 在线天堂最新版资源| 国产色婷婷99| 亚洲图色成人| 国产精品蜜桃在线观看| 国产久久久一区二区三区| 日日啪夜夜爽| 亚洲不卡免费看| 国产在线一区二区三区精| 亚洲美女搞黄在线观看| 国产精品国产三级国产专区5o| 欧美xxxx黑人xx丫x性爽| 一级二级三级毛片免费看| 校园人妻丝袜中文字幕| 黄片无遮挡物在线观看| 国产不卡一卡二| 爱豆传媒免费全集在线观看| 男女国产视频网站| 黑人高潮一二区| 午夜日本视频在线| 亚洲欧美日韩东京热| 别揉我奶头 嗯啊视频| 最近最新中文字幕大全电影3| 久久久欧美国产精品| 日日摸夜夜添夜夜添av毛片| 九草在线视频观看| 高清午夜精品一区二区三区| 亚洲欧美一区二区三区黑人 | 国产一级毛片在线| 国产麻豆成人av免费视频| 99热6这里只有精品| 亚洲av在线观看美女高潮| 亚洲天堂国产精品一区在线| 午夜精品一区二区三区免费看| 国产色爽女视频免费观看| 99热这里只有是精品在线观看| 午夜亚洲福利在线播放| 爱豆传媒免费全集在线观看| 男女国产视频网站| 69av精品久久久久久| 亚洲欧美清纯卡通| 亚洲精品乱久久久久久| 亚洲精品乱码久久久久久按摩| 伦理电影大哥的女人| 免费av观看视频| 乱码一卡2卡4卡精品| 内地一区二区视频在线| 欧美日韩亚洲高清精品| 内地一区二区视频在线| 国产综合懂色| 久久99热6这里只有精品| 国产麻豆成人av免费视频| 中文乱码字字幕精品一区二区三区 | 在线观看免费高清a一片| 精品久久久精品久久久| 偷拍熟女少妇极品色| 一级片'在线观看视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲在线观看片| 最近最新中文字幕大全电影3| 尤物成人国产欧美一区二区三区| 成年av动漫网址| 免费黄网站久久成人精品| 亚洲国产成人一精品久久久| 嘟嘟电影网在线观看| 国产单亲对白刺激| 韩国高清视频一区二区三区| 日韩一区二区三区影片| 51国产日韩欧美| 人妻少妇偷人精品九色| 欧美日韩精品成人综合77777| 国产高清不卡午夜福利| 久久精品人妻少妇| 99热网站在线观看| 国产精品女同一区二区软件| 老司机影院成人| 国产一级毛片在线| av专区在线播放| 亚洲色图av天堂| 亚洲乱码一区二区免费版| 久久这里只有精品中国| av播播在线观看一区| 日韩制服骚丝袜av| 久久综合国产亚洲精品| 秋霞在线观看毛片| 80岁老熟妇乱子伦牲交| 成年女人看的毛片在线观看| www.色视频.com| 超碰av人人做人人爽久久| 国产精品嫩草影院av在线观看| 日本欧美国产在线视频| 一本一本综合久久| 国产麻豆成人av免费视频| 国产乱来视频区| 伊人久久国产一区二区| 伦精品一区二区三区| 2022亚洲国产成人精品| 又粗又硬又长又爽又黄的视频| 国产成人a∨麻豆精品| 麻豆av噜噜一区二区三区| 亚洲天堂国产精品一区在线| 日韩视频在线欧美| 国产午夜精品一二区理论片| 国内精品美女久久久久久| 国产精品久久久久久精品电影| 国产 亚洲一区二区三区 | 一级毛片 在线播放| 日本免费a在线| 国产中年淑女户外野战色| 激情五月婷婷亚洲| 男女边吃奶边做爰视频| 男插女下体视频免费在线播放| 老女人水多毛片| 看免费成人av毛片| 成人亚洲欧美一区二区av| 国产老妇女一区| 九九在线视频观看精品| av国产免费在线观看| 国产亚洲最大av| 熟女人妻精品中文字幕| 亚洲三级黄色毛片| 国内精品美女久久久久久| 亚洲av福利一区| 五月天丁香电影| 久久精品久久久久久久性| 青春草国产在线视频| 男人狂女人下面高潮的视频| 精品99又大又爽又粗少妇毛片| 又爽又黄a免费视频| 久久精品熟女亚洲av麻豆精品 | 国产精品久久久久久久电影| 午夜福利网站1000一区二区三区| 日本色播在线视频| 我的女老师完整版在线观看| 卡戴珊不雅视频在线播放| 麻豆精品久久久久久蜜桃| 18禁动态无遮挡网站| 亚洲人成网站在线播| 在线观看人妻少妇| 日韩一区二区视频免费看| 高清视频免费观看一区二区 | 久久久精品欧美日韩精品| 国产黄片视频在线免费观看| 国产免费视频播放在线视频 | 99久久九九国产精品国产免费| 免费观看av网站的网址| 日韩av在线免费看完整版不卡| 深夜a级毛片| 国产精品久久久久久精品电影| 高清av免费在线| 综合色丁香网| 丰满少妇做爰视频| 能在线免费观看的黄片| 亚洲一区高清亚洲精品| 直男gayav资源| 午夜免费观看性视频| 国产成人freesex在线| 国产精品一区二区性色av| 亚洲精品aⅴ在线观看| av免费在线看不卡| 青春草亚洲视频在线观看| 在现免费观看毛片| 男人狂女人下面高潮的视频| 免费看不卡的av| 男的添女的下面高潮视频| av在线天堂中文字幕| 青春草国产在线视频| 麻豆成人av视频| 2021天堂中文幕一二区在线观| 国产精品女同一区二区软件| 国产单亲对白刺激| 网址你懂的国产日韩在线| www.色视频.com| 性色avwww在线观看| 又黄又爽又刺激的免费视频.| 精品国产露脸久久av麻豆 | 精品一区二区免费观看| 丰满乱子伦码专区| 成年人午夜在线观看视频 | 欧美xxⅹ黑人| 国产成人免费观看mmmm| 校园人妻丝袜中文字幕| 精品一区二区三区视频在线| 午夜精品国产一区二区电影 | 91久久精品国产一区二区成人| 亚洲av电影不卡..在线观看| 白带黄色成豆腐渣| 国产一区二区在线观看日韩| 久久精品久久精品一区二区三区| 成人二区视频| 成年版毛片免费区| 日本与韩国留学比较| 国产一区二区在线观看日韩| 日日干狠狠操夜夜爽| 晚上一个人看的免费电影| 久久精品国产亚洲av涩爱| 免费观看av网站的网址| 日韩制服骚丝袜av| 国产精品熟女久久久久浪| 成人综合一区亚洲| 免费看a级黄色片| 亚洲国产精品成人综合色| 又黄又爽又刺激的免费视频.| 黄色一级大片看看| 亚洲av中文av极速乱| a级毛片免费高清观看在线播放| 午夜激情欧美在线| 亚洲av中文字字幕乱码综合| 高清在线视频一区二区三区| 丝袜喷水一区| 天堂影院成人在线观看| 不卡视频在线观看欧美| 精品久久久久久久久亚洲| 欧美成人a在线观看| 国产高清三级在线| 成人毛片a级毛片在线播放| 亚洲最大成人av| 天堂俺去俺来也www色官网 | 国产在线一区二区三区精| 国产亚洲精品av在线| 亚洲第一区二区三区不卡| 最新中文字幕久久久久| 亚洲av一区综合| 春色校园在线视频观看| 狠狠精品人妻久久久久久综合| 国产成人午夜福利电影在线观看| 欧美日韩一区二区视频在线观看视频在线 | av在线老鸭窝| 日本-黄色视频高清免费观看| 街头女战士在线观看网站| 亚洲熟妇中文字幕五十中出| 色播亚洲综合网| 国产69精品久久久久777片| 国产日韩欧美在线精品| 99久久九九国产精品国产免费| 2022亚洲国产成人精品| 一级片'在线观看视频| 国模一区二区三区四区视频| 超碰97精品在线观看| 国产伦精品一区二区三区四那| 国产老妇女一区| 日韩伦理黄色片| or卡值多少钱| 简卡轻食公司| 国产大屁股一区二区在线视频| 久久久久国产网址| 免费看a级黄色片| 国国产精品蜜臀av免费| 毛片女人毛片| 亚洲精品一区蜜桃| 人妻制服诱惑在线中文字幕| 国产黄色免费在线视频| 一边亲一边摸免费视频| 春色校园在线视频观看| 色尼玛亚洲综合影院| www.av在线官网国产| 只有这里有精品99| 国产精品日韩av在线免费观看| 亚洲四区av| 久久精品国产鲁丝片午夜精品| 国内精品一区二区在线观看| 少妇人妻一区二区三区视频| 欧美zozozo另类| 晚上一个人看的免费电影| 国产精品.久久久| 黑人高潮一二区| 久久久久性生活片| 国产av码专区亚洲av| 国产精品三级大全| 男人狂女人下面高潮的视频| 色网站视频免费| 国产精品日韩av在线免费观看| 欧美高清性xxxxhd video| 纵有疾风起免费观看全集完整版 | 久久精品夜色国产| 亚洲在久久综合| 在线免费观看的www视频| 非洲黑人性xxxx精品又粗又长| 国产精品.久久久| 国产成人精品婷婷| 少妇熟女欧美另类| 国产成人一区二区在线| eeuss影院久久| 久久久精品94久久精品| 午夜激情欧美在线| kizo精华| 99久久精品热视频| 成年人午夜在线观看视频 | av.在线天堂| 亚洲最大成人av| 天堂俺去俺来也www色官网 | h日本视频在线播放| 91av网一区二区| 日日干狠狠操夜夜爽| 亚洲av在线观看美女高潮| 中文字幕av成人在线电影| 青青草视频在线视频观看| 99久久人妻综合| 亚洲激情五月婷婷啪啪| 久久6这里有精品| av在线老鸭窝| av网站免费在线观看视频 | 午夜免费激情av| 亚洲精华国产精华液的使用体验| 精品国产三级普通话版| 久久97久久精品| 亚洲国产成人一精品久久久| 日本-黄色视频高清免费观看| 亚洲av电影在线观看一区二区三区 | 亚洲内射少妇av| 99热这里只有精品一区| 91久久精品电影网| 日本爱情动作片www.在线观看| 国产久久久一区二区三区| 三级经典国产精品| 少妇被粗大猛烈的视频| 丝袜美腿在线中文| 国产乱人偷精品视频| 亚洲电影在线观看av| 少妇裸体淫交视频免费看高清| 亚洲不卡免费看| 91午夜精品亚洲一区二区三区| av一本久久久久| 国产69精品久久久久777片| 亚洲在久久综合| 亚洲欧洲国产日韩| 性插视频无遮挡在线免费观看| av黄色大香蕉| 天堂中文最新版在线下载 | 黄色日韩在线| 我的女老师完整版在线观看| 成人av在线播放网站| 国产伦理片在线播放av一区| 看非洲黑人一级黄片| 熟女人妻精品中文字幕| 亚洲成人av在线免费| 久久综合国产亚洲精品| 夜夜爽夜夜爽视频| 国产午夜精品久久久久久一区二区三区| 国产黄色小视频在线观看| 99热这里只有精品一区| 熟女电影av网| 国产黄a三级三级三级人| 99热全是精品| 伊人久久精品亚洲午夜| 亚洲欧美成人综合另类久久久| 最近手机中文字幕大全| 日韩一区二区视频免费看| 人人妻人人澡欧美一区二区| 亚洲婷婷狠狠爱综合网| 中文字幕免费在线视频6| 国产av码专区亚洲av| 日本猛色少妇xxxxx猛交久久| 99热这里只有是精品在线观看| 国内精品一区二区在线观看| 亚洲综合色惰| 成人性生交大片免费视频hd| 亚洲国产日韩欧美精品在线观看| 精华霜和精华液先用哪个| 性色avwww在线观看| 日韩一区二区视频免费看| 精品国内亚洲2022精品成人| 国产成人精品福利久久| 亚洲国产欧美在线一区| 中文字幕亚洲精品专区| 亚洲精品日韩在线中文字幕| 性色avwww在线观看| 国产成人一区二区在线| 日韩三级伦理在线观看| 亚洲精品中文字幕在线视频 | 观看免费一级毛片| 婷婷色麻豆天堂久久| 插阴视频在线观看视频| 亚洲在久久综合| 日韩大片免费观看网站| 国产黄色视频一区二区在线观看| 欧美激情久久久久久爽电影| 欧美日韩精品成人综合77777| 在现免费观看毛片| 夫妻性生交免费视频一级片| 国产高清国产精品国产三级 | 久久这里只有精品中国| 国产精品嫩草影院av在线观看| 欧美日本视频| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲精品久久久com| 伦精品一区二区三区| 少妇的逼好多水| 精品人妻一区二区三区麻豆| 亚洲自拍偷在线| 亚洲精品第二区| 午夜福利在线观看免费完整高清在| 美女内射精品一级片tv| 国产探花极品一区二区| 天堂√8在线中文| 一级av片app| 91午夜精品亚洲一区二区三区| 国语对白做爰xxxⅹ性视频网站| 人妻制服诱惑在线中文字幕| 深爱激情五月婷婷| 中文字幕久久专区| 亚洲综合色惰| 国产精品久久久久久久久免| 亚洲精品亚洲一区二区| 久久久色成人| 网址你懂的国产日韩在线| 我要看日韩黄色一级片| 国产亚洲av片在线观看秒播厂 | 爱豆传媒免费全集在线观看| 亚洲精品中文字幕在线视频 | 老女人水多毛片| 国产 一区精品| 国内精品宾馆在线| 日韩,欧美,国产一区二区三区| 国产成人精品婷婷| 尤物成人国产欧美一区二区三区| 亚洲国产精品专区欧美| 九九爱精品视频在线观看| 婷婷色综合大香蕉| 麻豆av噜噜一区二区三区| 亚洲精品自拍成人| 久久久久精品性色| 免费观看在线日韩| 欧美日韩国产mv在线观看视频 | 国产亚洲午夜精品一区二区久久 | 国产单亲对白刺激| 寂寞人妻少妇视频99o| 最近中文字幕高清免费大全6| 美女被艹到高潮喷水动态| 菩萨蛮人人尽说江南好唐韦庄| 麻豆精品久久久久久蜜桃| 午夜福利在线观看免费完整高清在| 久久精品久久精品一区二区三区| 天堂av国产一区二区熟女人妻| 国产精品国产三级国产av玫瑰| 中文字幕亚洲精品专区| 日产精品乱码卡一卡2卡三| 国内精品宾馆在线| 免费观看av网站的网址| av国产久精品久网站免费入址| 看免费成人av毛片| 色综合站精品国产| 成人美女网站在线观看视频| 久久久久精品久久久久真实原创| 国产成人aa在线观看| 全区人妻精品视频| 女人十人毛片免费观看3o分钟| 久久久久久久久中文| 最近最新中文字幕免费大全7| 亚洲一区高清亚洲精品| 亚洲av日韩在线播放| 精品久久久久久久末码| 国产探花极品一区二区| 国产高清国产精品国产三级 | 精品久久久久久电影网| 欧美极品一区二区三区四区| 欧美bdsm另类| 亚洲av在线观看美女高潮| 免费av毛片视频| 丝袜喷水一区| 菩萨蛮人人尽说江南好唐韦庄| 日日啪夜夜爽| 99热网站在线观看| 老女人水多毛片| 又粗又硬又长又爽又黄的视频| 日韩电影二区| 亚洲精品国产成人久久av| 伊人久久精品亚洲午夜| 国产免费视频播放在线视频 | 日韩成人av中文字幕在线观看| 中国美白少妇内射xxxbb| 青春草亚洲视频在线观看| 日日摸夜夜添夜夜添av毛片| 一级二级三级毛片免费看| 欧美日本视频| 18禁在线播放成人免费| 久久久久久久久久成人| 免费大片18禁| 日日摸夜夜添夜夜添av毛片| 成人特级av手机在线观看| 日日摸夜夜添夜夜添av毛片| 九九爱精品视频在线观看| 亚洲欧洲国产日韩| 丝瓜视频免费看黄片| 国产精品久久久久久精品电影| 国产在视频线在精品| 五月玫瑰六月丁香| 可以在线观看毛片的网站| 国产精品三级大全| 亚洲国产日韩欧美精品在线观看| 青青草视频在线视频观看| 不卡视频在线观看欧美| 日日干狠狠操夜夜爽| 久久久久久久久久人人人人人人| 亚洲精品国产av成人精品| 日本免费在线观看一区| 大片免费播放器 马上看| 老司机影院毛片| 国模一区二区三区四区视频| 不卡视频在线观看欧美| 夫妻午夜视频| 中文字幕久久专区| 国国产精品蜜臀av免费| 国产伦精品一区二区三区四那| 一级毛片电影观看| 亚洲欧美成人精品一区二区| 欧美一区二区亚洲| 天天一区二区日本电影三级| 亚洲aⅴ乱码一区二区在线播放| 99久久九九国产精品国产免费| 久久午夜福利片| 身体一侧抽搐| 久久久久精品性色| 久久热精品热| 97精品久久久久久久久久精品| 91午夜精品亚洲一区二区三区| 少妇熟女欧美另类| 日产精品乱码卡一卡2卡三| 美女高潮的动态| eeuss影院久久| 欧美极品一区二区三区四区| 亚洲欧美成人综合另类久久久| 九九爱精品视频在线观看| 最近最新中文字幕免费大全7| 联通29元200g的流量卡| 狂野欧美激情性xxxx在线观看| 免费大片黄手机在线观看| 插阴视频在线观看视频| 边亲边吃奶的免费视频| 久久久久性生活片|