• <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ù)值模擬
    精品酒店卫生间| 久久久欧美国产精品| 国产精品免费视频内射| 热re99久久精品国产66热6| 老司机深夜福利视频在线观看 | 一级黄片播放器| 国产一区二区 视频在线| 如日韩欧美国产精品一区二区三区| 亚洲在久久综合| 久久国产精品男人的天堂亚洲| 9热在线视频观看99| 综合色丁香网| av片东京热男人的天堂| 精品免费久久久久久久清纯 | 国产精品无大码| 一区福利在线观看| 欧美精品亚洲一区二区| 亚洲精品av麻豆狂野| 国产欧美日韩综合在线一区二区| 国产男女超爽视频在线观看| tube8黄色片| 国产成人欧美| 久久韩国三级中文字幕| √禁漫天堂资源中文www| 免费观看a级毛片全部| 人体艺术视频欧美日本| 久久亚洲国产成人精品v| 亚洲成人免费av在线播放| 天堂俺去俺来也www色官网| 美女福利国产在线| 国产精品一区二区在线不卡| 久久99一区二区三区| 考比视频在线观看| 久久精品国产a三级三级三级| 亚洲自偷自拍图片 自拍| 亚洲成人手机| 好男人视频免费观看在线| 日韩一卡2卡3卡4卡2021年| 天天操日日干夜夜撸| 人人妻人人澡人人爽人人夜夜| 国产又爽黄色视频| 亚洲精品中文字幕在线视频| 国语对白做爰xxxⅹ性视频网站| 一级毛片电影观看| 午夜激情久久久久久久| 成人国产av品久久久| 久久久精品免费免费高清| 丁香六月天网| 国产在视频线精品| 亚洲国产中文字幕在线视频| 国产精品偷伦视频观看了| 亚洲人成77777在线视频| 深夜精品福利| 国产野战对白在线观看| 另类精品久久| 咕卡用的链子| 午夜激情久久久久久久| 一级毛片 在线播放| 亚洲国产欧美一区二区综合| 国产精品.久久久| 日韩 欧美 亚洲 中文字幕| 伊人久久国产一区二区| 9191精品国产免费久久| 亚洲av电影在线观看一区二区三区| 国产成人欧美在线观看 | 欧美日韩亚洲高清精品| 自线自在国产av| 亚洲四区av| 亚洲 欧美一区二区三区| 国产成人免费观看mmmm| 91国产中文字幕| 亚洲国产毛片av蜜桃av| 婷婷色综合www| 精品一区在线观看国产| 又大又黄又爽视频免费| 99香蕉大伊视频| 男的添女的下面高潮视频| 女人久久www免费人成看片| 另类亚洲欧美激情| 在线观看免费视频网站a站| av在线老鸭窝| 韩国高清视频一区二区三区| 侵犯人妻中文字幕一二三四区| av一本久久久久| 欧美日韩亚洲国产一区二区在线观看 | 一二三四在线观看免费中文在| 国产av国产精品国产| 美女扒开内裤让男人捅视频| 午夜av观看不卡| 电影成人av| 黄色视频不卡| 天天躁狠狠躁夜夜躁狠狠躁| 狂野欧美激情性bbbbbb| 国产精品久久久久成人av| 一区二区三区激情视频| 亚洲欧美日韩另类电影网站| 国产精品免费视频内射| 黄色毛片三级朝国网站| 亚洲成国产人片在线观看| 如日韩欧美国产精品一区二区三区| 超色免费av| 丝袜人妻中文字幕| 欧美精品人与动牲交sv欧美| 精品国产露脸久久av麻豆| 大香蕉久久网| 久久天堂一区二区三区四区| 深夜精品福利| 欧美变态另类bdsm刘玥| 两个人免费观看高清视频| 极品少妇高潮喷水抽搐| 妹子高潮喷水视频| 国产亚洲av高清不卡| 亚洲精品日本国产第一区| 国产精品久久久久成人av| 91精品国产国语对白视频| 精品一区二区三区av网在线观看 | 亚洲欧美成人精品一区二区| 亚洲图色成人| 赤兔流量卡办理| 一级爰片在线观看| 18禁动态无遮挡网站| 麻豆精品久久久久久蜜桃| 亚洲欧美一区二区三区国产| 一边摸一边做爽爽视频免费| 国产一级毛片在线| 韩国高清视频一区二区三区| 一边摸一边抽搐一进一出视频| 美国免费a级毛片| 亚洲欧美日韩另类电影网站| 中文字幕制服av| 青春草国产在线视频| 不卡av一区二区三区| 天堂8中文在线网| 精品第一国产精品| 欧美少妇被猛烈插入视频| 777米奇影视久久| 久久久精品国产亚洲av高清涩受| 精品国产乱码久久久久久小说| 香蕉国产在线看| 国产精品成人在线| 别揉我奶头~嗯~啊~动态视频 | 毛片一级片免费看久久久久| 少妇猛男粗大的猛烈进出视频| 我的亚洲天堂| 亚洲欧洲日产国产| 中文精品一卡2卡3卡4更新| 欧美日韩亚洲国产一区二区在线观看 | av片东京热男人的天堂| 日韩制服骚丝袜av| 精品国产乱码久久久久久小说| 免费在线观看黄色视频的| 久久性视频一级片| 一区二区三区精品91| 亚洲精品自拍成人| 视频在线观看一区二区三区| 只有这里有精品99| 久久久亚洲精品成人影院| 天天操日日干夜夜撸| 国产福利在线免费观看视频| 久久久久久久久久久免费av| 制服丝袜香蕉在线| 国产极品粉嫩免费观看在线| 国产免费福利视频在线观看| 久久久亚洲精品成人影院| 日韩一本色道免费dvd| 欧美日韩成人在线一区二区| 国产极品天堂在线| 人人妻人人澡人人看| 日日摸夜夜添夜夜爱| 国产成人免费无遮挡视频| 看免费av毛片| 免费观看性生交大片5| 日本爱情动作片www.在线观看| 观看美女的网站| 夫妻午夜视频| 国产精品免费大片| 中文字幕av电影在线播放| 丁香六月欧美| av国产久精品久网站免费入址| 大片免费播放器 马上看| 欧美乱码精品一区二区三区| 午夜福利视频精品| 欧美成人午夜精品| 性色av一级| 丁香六月天网| 蜜桃国产av成人99| 日韩一区二区三区影片| 999精品在线视频| 18禁动态无遮挡网站| 老司机在亚洲福利影院| 日本欧美国产在线视频| 丝袜在线中文字幕| 青春草国产在线视频| 午夜精品国产一区二区电影| 美女中出高潮动态图| 亚洲精华国产精华液的使用体验| 国产成人午夜福利电影在线观看| 欧美日本中文国产一区发布| 中文字幕人妻丝袜制服| 制服丝袜香蕉在线| 69精品国产乱码久久久| 国产国语露脸激情在线看| 日韩大片免费观看网站| 国产精品免费大片| 亚洲国产精品一区三区| 亚洲情色 制服丝袜| 亚洲精品av麻豆狂野| 中文字幕av电影在线播放| 高清欧美精品videossex| 咕卡用的链子| 婷婷成人精品国产| 亚洲少妇的诱惑av| 无遮挡黄片免费观看| 另类亚洲欧美激情| 国产男女超爽视频在线观看| 日韩av在线免费看完整版不卡| 国产黄色免费在线视频| 亚洲免费av在线视频| 一本—道久久a久久精品蜜桃钙片| 国产成人啪精品午夜网站| 男女国产视频网站| 啦啦啦视频在线资源免费观看| 欧美人与性动交α欧美软件| 亚洲三区欧美一区| 悠悠久久av| 日日啪夜夜爽| tube8黄色片| 午夜福利乱码中文字幕| 尾随美女入室| 亚洲欧美一区二区三区黑人| 久久久久久久久久久久大奶| 观看av在线不卡| 亚洲欧美成人综合另类久久久| 免费久久久久久久精品成人欧美视频| 精品一区在线观看国产| 99久久综合免费| 欧美 亚洲 国产 日韩一| 中文字幕亚洲精品专区| 女人高潮潮喷娇喘18禁视频| 深夜精品福利| 丝袜人妻中文字幕| 亚洲美女视频黄频| 色视频在线一区二区三区| 一级黄片播放器| 亚洲精品久久成人aⅴ小说| 色网站视频免费| 少妇人妻 视频| 人人妻人人添人人爽欧美一区卜| 丝袜脚勾引网站| 欧美精品高潮呻吟av久久| 在线天堂最新版资源| 又大又黄又爽视频免费| 欧美日韩视频高清一区二区三区二| 欧美精品亚洲一区二区| 精品一品国产午夜福利视频| 午夜福利网站1000一区二区三区| 18禁裸乳无遮挡动漫免费视频| 两个人看的免费小视频| 精品国产超薄肉色丝袜足j| 97在线人人人人妻| 黄色一级大片看看| 两性夫妻黄色片| 一边亲一边摸免费视频| a 毛片基地| 国产精品麻豆人妻色哟哟久久| 高清在线视频一区二区三区| 国产成人一区二区在线| 欧美最新免费一区二区三区| 制服丝袜香蕉在线| 久久久亚洲精品成人影院| bbb黄色大片| 亚洲av福利一区| 看免费av毛片| 水蜜桃什么品种好| 丝瓜视频免费看黄片| 欧美老熟妇乱子伦牲交| 熟女av电影| 一级毛片黄色毛片免费观看视频| 国产黄色免费在线视频| av视频免费观看在线观看| 99国产精品免费福利视频| 欧美日韩亚洲综合一区二区三区_| 免费女性裸体啪啪无遮挡网站| 久久天躁狠狠躁夜夜2o2o | 视频区图区小说| 一本色道久久久久久精品综合| 日韩成人av中文字幕在线观看| 亚洲综合精品二区| 如何舔出高潮| 男人添女人高潮全过程视频| 精品人妻熟女毛片av久久网站| 国产精品偷伦视频观看了| 亚洲一级一片aⅴ在线观看| 久久99热这里只频精品6学生| 2021少妇久久久久久久久久久| 欧美亚洲日本最大视频资源| 天天躁狠狠躁夜夜躁狠狠躁| 91精品三级在线观看| 午夜老司机福利片| xxxhd国产人妻xxx| 老司机影院成人| 波多野结衣av一区二区av| 久久精品aⅴ一区二区三区四区| 国产免费现黄频在线看| 精品福利永久在线观看| 18禁裸乳无遮挡动漫免费视频| 热re99久久国产66热| 欧美日韩亚洲高清精品| 少妇被粗大的猛进出69影院| 亚洲一码二码三码区别大吗| 国产麻豆69| 综合色丁香网| 韩国高清视频一区二区三区| 免费在线观看完整版高清| 天天躁夜夜躁狠狠久久av| 爱豆传媒免费全集在线观看| 午夜日韩欧美国产| 欧美 日韩 精品 国产| 免费在线观看完整版高清| 国产精品国产三级专区第一集| 男女高潮啪啪啪动态图| 日本欧美国产在线视频| 无遮挡黄片免费观看| 香蕉丝袜av| 99热国产这里只有精品6| 久久97久久精品| 在线亚洲精品国产二区图片欧美| 亚洲第一青青草原| 亚洲一区二区三区欧美精品| 国产一区亚洲一区在线观看| 我要看黄色一级片免费的| 国产又爽黄色视频| 肉色欧美久久久久久久蜜桃| 国产成人免费观看mmmm| 99久久综合免费| 亚洲国产精品国产精品| 欧美激情极品国产一区二区三区| 久久性视频一级片| 久久热在线av| 丰满饥渴人妻一区二区三| 18禁国产床啪视频网站| 久久久亚洲精品成人影院| 大码成人一级视频| 天天躁日日躁夜夜躁夜夜| 亚洲第一av免费看| 久久久久久久久久久久大奶| 久久久久精品国产欧美久久久 | 成人午夜精彩视频在线观看| 满18在线观看网站| 欧美国产精品一级二级三级| a级毛片黄视频| 亚洲在久久综合| 99热网站在线观看| 亚洲成人手机| 男人舔女人的私密视频| 欧美激情高清一区二区三区 | 国产亚洲一区二区精品| av视频免费观看在线观看| 亚洲美女黄色视频免费看| 午夜日韩欧美国产| 国产精品久久久av美女十八| av有码第一页| 伊人亚洲综合成人网| 一边亲一边摸免费视频| av不卡在线播放| 久久久国产欧美日韩av| 国产av码专区亚洲av| 久久精品久久久久久噜噜老黄| 国产1区2区3区精品| 18禁动态无遮挡网站| 亚洲人成网站在线观看播放| 国产精品.久久久| 免费观看a级毛片全部| 午夜福利视频精品| 91老司机精品| 亚洲精品中文字幕在线视频| 国产一区有黄有色的免费视频| 91精品三级在线观看| 看非洲黑人一级黄片| 亚洲精品久久午夜乱码| 丰满饥渴人妻一区二区三| 在线观看免费视频网站a站| 欧美精品一区二区大全| 国产高清国产精品国产三级| 搡老乐熟女国产| 777久久人妻少妇嫩草av网站| 嫩草影视91久久| 99国产精品免费福利视频| 国产精品秋霞免费鲁丝片| 久久午夜综合久久蜜桃| 欧美亚洲日本最大视频资源| 亚洲精品美女久久久久99蜜臀 | 国产伦理片在线播放av一区| 欧美日韩福利视频一区二区| 国产成人av激情在线播放| 亚洲精品一区蜜桃| 啦啦啦 在线观看视频| 如日韩欧美国产精品一区二区三区| 91老司机精品| 婷婷成人精品国产| 免费不卡黄色视频| av在线播放精品| 狠狠精品人妻久久久久久综合| 国产一区二区激情短视频 | 丝袜美足系列| 天天躁夜夜躁狠狠躁躁| 国产av一区二区精品久久| 国产福利在线免费观看视频| 一本色道久久久久久精品综合| 国产精品嫩草影院av在线观看| 一个人免费看片子| 91国产中文字幕| 黄色怎么调成土黄色| 日韩中文字幕欧美一区二区 | 久久精品国产a三级三级三级| 91精品三级在线观看| 久久这里只有精品19| 国产老妇伦熟女老妇高清| 美女中出高潮动态图| av网站免费在线观看视频| 天堂8中文在线网| 亚洲五月色婷婷综合| 欧美激情高清一区二区三区 | 亚洲欧美一区二区三区国产| videos熟女内射| 国产成人精品久久久久久| 色综合欧美亚洲国产小说| 超碰成人久久| 在线观看一区二区三区激情| av在线观看视频网站免费| 一边亲一边摸免费视频| 建设人人有责人人尽责人人享有的| 国产一区二区激情短视频 | 建设人人有责人人尽责人人享有的| 老熟女久久久| 五月开心婷婷网| 免费在线观看黄色视频的| 夜夜骑夜夜射夜夜干| 尾随美女入室| 欧美亚洲 丝袜 人妻 在线| 亚洲精品av麻豆狂野| 伊人久久国产一区二区| 国产av码专区亚洲av| 大片电影免费在线观看免费| 色婷婷久久久亚洲欧美| 午夜免费观看性视频| 国产 一区精品| 久久青草综合色| 人妻 亚洲 视频| 久久久久久久久久久久大奶| 亚洲国产欧美在线一区| 亚洲成av片中文字幕在线观看| av在线观看视频网站免费| 啦啦啦在线免费观看视频4| 精品视频人人做人人爽| 亚洲国产日韩一区二区| 欧美日韩亚洲国产一区二区在线观看 | 久久精品久久久久久久性| 赤兔流量卡办理| 久久精品久久久久久噜噜老黄| 亚洲婷婷狠狠爱综合网| 99国产精品免费福利视频| 免费在线观看黄色视频的| 亚洲三区欧美一区| kizo精华| 国产精品一国产av| 欧美日韩精品网址| 久久久欧美国产精品| 亚洲欧美成人综合另类久久久| 观看美女的网站| 亚洲欧美色中文字幕在线| 国产精品久久久久久人妻精品电影 | 成年人午夜在线观看视频| 欧美黄色片欧美黄色片| 十分钟在线观看高清视频www| 狂野欧美激情性bbbbbb| 久久亚洲国产成人精品v| 成人黄色视频免费在线看| 女人爽到高潮嗷嗷叫在线视频| 丝袜喷水一区| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美清纯卡通| 宅男免费午夜| 制服诱惑二区| 99久久99久久久精品蜜桃| 一区福利在线观看| 午夜av观看不卡| av女优亚洲男人天堂| 久久精品aⅴ一区二区三区四区| 免费人妻精品一区二区三区视频| 国产1区2区3区精品| av在线app专区| 亚洲精品自拍成人| 99久久99久久久精品蜜桃| 人妻 亚洲 视频| 国产成人欧美在线观看 | av在线老鸭窝| 成人18禁高潮啪啪吃奶动态图| 久久久久国产精品人妻一区二区| 国产不卡av网站在线观看| 最近中文字幕2019免费版| 亚洲欧洲国产日韩| 99九九在线精品视频| 亚洲美女搞黄在线观看| 亚洲精品aⅴ在线观看| 国产成人免费无遮挡视频| 国产淫语在线视频| 日韩不卡一区二区三区视频在线| 欧美 亚洲 国产 日韩一| 在线精品无人区一区二区三| 99香蕉大伊视频| 欧美人与性动交α欧美软件| 欧美人与善性xxx| 国产精品免费大片| 啦啦啦啦在线视频资源| 色网站视频免费| 久久国产亚洲av麻豆专区| 精品亚洲成a人片在线观看| 国产福利在线免费观看视频| 亚洲成人手机| 两个人免费观看高清视频| 国产熟女午夜一区二区三区| 咕卡用的链子| 99国产精品免费福利视频| 中国三级夫妇交换| 欧美日韩亚洲高清精品| 男女免费视频国产| netflix在线观看网站| 另类亚洲欧美激情| av在线观看视频网站免费| 老司机亚洲免费影院| 亚洲,欧美,日韩| 日本猛色少妇xxxxx猛交久久| 天天躁夜夜躁狠狠躁躁| 日韩大码丰满熟妇| 欧美日韩国产mv在线观看视频| 欧美黑人欧美精品刺激| 一本色道久久久久久精品综合| 建设人人有责人人尽责人人享有的| 色网站视频免费| 丝袜美足系列| 人妻 亚洲 视频| 中文字幕色久视频| 婷婷色av中文字幕| 国产福利在线免费观看视频| 成人漫画全彩无遮挡| 极品少妇高潮喷水抽搐| 操出白浆在线播放| 美女国产高潮福利片在线看| 国产高清不卡午夜福利| 我的亚洲天堂| 热99久久久久精品小说推荐| 亚洲欧美成人综合另类久久久| 超碰97精品在线观看| 欧美日韩一区二区视频在线观看视频在线| 日本爱情动作片www.在线观看| 美女中出高潮动态图| 日日撸夜夜添| 丰满饥渴人妻一区二区三| 欧美日韩一区二区视频在线观看视频在线| 制服人妻中文乱码| 狠狠精品人妻久久久久久综合| 国产精品一国产av| 久久久久精品性色| 香蕉丝袜av| 久久青草综合色| 欧美黑人欧美精品刺激| 啦啦啦在线观看免费高清www| 十八禁高潮呻吟视频| 国产极品天堂在线| 一边摸一边抽搐一进一出视频| 国产在线视频一区二区| 青春草亚洲视频在线观看| a级片在线免费高清观看视频| 一边亲一边摸免费视频| 欧美在线黄色| 久久久久视频综合| 熟妇人妻不卡中文字幕| 欧美日韩一级在线毛片| 亚洲欧美一区二区三区国产| 久久久亚洲精品成人影院| 热re99久久国产66热| av有码第一页| 两个人看的免费小视频| 成人午夜精彩视频在线观看| 亚洲一级一片aⅴ在线观看| 久久人人97超碰香蕉20202| 日本av免费视频播放| 精品国产超薄肉色丝袜足j| 在线看a的网站| 在现免费观看毛片| 大陆偷拍与自拍| 丰满少妇做爰视频| 岛国毛片在线播放| av国产精品久久久久影院| 久久久精品免费免费高清| 欧美中文综合在线视频| 韩国精品一区二区三区| 精品国产乱码久久久久久小说| a 毛片基地| 亚洲国产欧美一区二区综合| 热re99久久国产66热| 午夜免费男女啪啪视频观看| 一级片'在线观看视频| 免费少妇av软件| 如日韩欧美国产精品一区二区三区| 亚洲精品成人av观看孕妇| 亚洲国产欧美一区二区综合| 中文乱码字字幕精品一区二区三区| 久久人妻熟女aⅴ| av国产久精品久网站免费入址| 免费高清在线观看日韩| 黄色怎么调成土黄色| 亚洲,欧美,日韩| svipshipincom国产片|