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

    準(zhǔn)東煤中鈉高溫形態(tài)的化學(xué)動(dòng)力學(xué)模擬

    2017-10-17 09:32:34季杰強(qiáng)程樂鳴劉炎泉王勤輝駱仲泱
    動(dòng)力工程學(xué)報(bào) 2017年10期
    關(guān)鍵詞:準(zhǔn)東堿金屬氣相

    季杰強(qiáng), 程樂鳴, 劉炎泉, 王勤輝, 駱仲泱

    (浙江大學(xué) 熱能工程研究所,能源清潔利用國家重點(diǎn)實(shí)驗(yàn)室,杭州 310027)

    準(zhǔn)東煤中鈉高溫形態(tài)的化學(xué)動(dòng)力學(xué)模擬

    季杰強(qiáng), 程樂鳴, 劉炎泉, 王勤輝, 駱仲泱

    (浙江大學(xué) 熱能工程研究所,能源清潔利用國家重點(diǎn)實(shí)驗(yàn)室,杭州 310027)

    采用化學(xué)動(dòng)力學(xué)軟件Chemkin對(duì)堿金屬含量高的準(zhǔn)東煤在燃燒過程中氣相鈉的高溫形態(tài)分布進(jìn)行模擬計(jì)算,將模型計(jì)算結(jié)果與文獻(xiàn)中的試驗(yàn)結(jié)果進(jìn)行對(duì)比,研究了溫度、過量空氣系數(shù)α和壓力對(duì)堿金屬鈉在高溫下形態(tài)分布的影響,并對(duì)含鈉氣相產(chǎn)物進(jìn)行了敏感性分析. 結(jié)果表明:計(jì)算結(jié)果與試驗(yàn)結(jié)果一致,驗(yàn)證了該模型的正確性;溫度對(duì)含鈉氣相產(chǎn)物的轉(zhuǎn)化有很大影響;當(dāng)α<1時(shí),含鈉氣相產(chǎn)物的組分分布受α的影響較大;壓力對(duì)含鈉氣相產(chǎn)物物質(zhì)的量分?jǐn)?shù)的影響較小.

    準(zhǔn)東煤; 鈉形態(tài); 化學(xué)動(dòng)力學(xué); 敏感性; 模擬

    Abstract: Chemical kinetics software Chemkin was applied to perform analog calculation on morphology distribution of gas-phase sodium at high temperatures in Zhundong coal with high alkali content during combustion process, and the calculation results were subsequently compared with reference data, so as to study the effects of following factors on the distribution of high-temperature morphology of alkali sodium, such as the temperature, excess air coefficient (α) and pressure, etc., and to analyze the sensitivity of relevant products. Results show that the calculation results agree well with actual measurements, proving the model to be correct; temperature significantly influences the conversion of sodium-containing gas-phase products; the component distribution of sodium-containing gas-phase products is remarkably influenced byαin the case ofα<1, which is slightly influenced by pressure.

    Key words: Zhundong coal; sodium morphology; chemical kinetics; sensitivity; simulation

    準(zhǔn)東煤田發(fā)現(xiàn)于2005年,其蘊(yùn)含的煤炭資源非常豐富,預(yù)測儲(chǔ)量可達(dá)3.9×1011t[1],占新疆煤炭儲(chǔ)量(2.19×1012t)的17%~18%,占全國煤炭儲(chǔ)量(5×1012~6×1012t)的7%~8%. 準(zhǔn)東煤的煤質(zhì)較好,水分和揮發(fā)分含量較高,灰分和硫含量較低,應(yīng)用前景非常廣闊. 然而研究表明,準(zhǔn)東煤中鈉、鉀等堿金屬含量偏高,國內(nèi)傳統(tǒng)用煤中鈉的質(zhì)量分?jǐn)?shù)一般在1%~2%(以灰分計(jì)),而準(zhǔn)東煤中鈉的質(zhì)量分?jǐn)?shù)一般在5%以上(以灰分計(jì)).燃用含高堿金屬和堿土金屬的煤會(huì)出現(xiàn)結(jié)渣、沾污、積灰和腐蝕等各種問題,嚴(yán)重影響鍋爐的正常運(yùn)行. Wang等[2]認(rèn)為準(zhǔn)東煤在燃燒過程中析出的堿金屬和堿土金屬化合物(NaCl、Na2SO4、CaSO4和K2SO4)[3]先在受熱面附近沉積,隨后部分硅鋁酸鹽被初始沉積層捕獲,這是沾污形成的主要原因. 因此,研究準(zhǔn)東煤這類典型的高鈉煤中氣相堿金屬的析出規(guī)律及高溫形態(tài)分布具有重要意義.

    由于氣相堿金屬在線測量比較困難,因此針對(duì)氣相堿金屬在高溫下形態(tài)分布的試驗(yàn)研究很少. Bl?sing等[4]采用分子束質(zhì)譜儀對(duì)煙氣進(jìn)行在線測量后發(fā)現(xiàn),58NaCl+的析出量與硫、氯元素的質(zhì)量比呈負(fù)相關(guān),同時(shí)也與鋁和硅的總含量呈負(fù)相關(guān). Bl?sing等[5]研究了水蒸氣對(duì)煤氣化過程中堿金屬等無機(jī)物析出的影響,發(fā)現(xiàn)加入水蒸氣將導(dǎo)致23Na+、39K+/39NaO+、56KOH+、58NaCl+和74KCl+的釋放量增加. Oleschko等[6]采用分子束質(zhì)譜儀分析了6種無煙煤在燃燒過程中氣相堿金屬的釋放情況,結(jié)果表明隨著溫度從800 ℃升至1 200 ℃,NaCl的釋放量減少,KCl的釋放量增加;增大壓力導(dǎo)致NaCl的釋放量減少. 然而,上述試驗(yàn)研究還存在一些不足:(1)具有相同質(zhì)荷比的物質(zhì)無法區(qū)分,如40NaOH+與40Ar+、39K+與39NaO+;(2)含量較低的堿金屬化合物信號(hào)易受環(huán)境噪音影響而難以測定,如142Na2SO4+等.這些都會(huì)影響測量結(jié)果的準(zhǔn)確性,因此通常需要結(jié)合模擬結(jié)果對(duì)堿金屬高溫形態(tài)的變化規(guī)律進(jìn)行完善.

    譚力[7]采用FactSage研究了稻草秸稈燃燒過程中燃燒溫度、爐內(nèi)壓力、過量空氣系數(shù)、燃料中的硫元素含量和氯元素含量等因素對(duì)堿金屬遷移的影響. 陳安和等[8]采用FactSage計(jì)算了秸稈、樹皮、木屑、廢木和橄欖渣5種生物質(zhì)在熱解條件下氯及堿金屬元素的分布和遷移特性. Wei等[9]采用FactSage研究了過量空氣系數(shù)和壓力對(duì)秸稈燃燒過程中氯及堿金屬釋放行為的影響. 目前,熱力學(xué)計(jì)算均基于以下前提:(1)不研究過程的速率,即忽略停留時(shí)間這個(gè)因素;(2)僅研究系統(tǒng)的平衡狀態(tài),忽略反應(yīng)機(jī)理. 對(duì)于爐內(nèi)燃燒或氣化過程,停留時(shí)間對(duì)堿金屬析出的影響很大. 另外,由于爐內(nèi)復(fù)雜的物理化學(xué)反應(yīng),系統(tǒng)難以達(dá)到或者維持在平衡狀態(tài),因此熱力學(xué)計(jì)算存在一定的局限性. 化學(xué)動(dòng)力學(xué)研究的是化學(xué)過程的速率及化學(xué)反應(yīng)機(jī)理,可有效彌補(bǔ)熱力學(xué)計(jì)算的不足,更接近爐內(nèi)真實(shí)狀況. 但目前的資料顯示,幾乎沒有應(yīng)用化學(xué)動(dòng)力學(xué)研究氣相堿金屬在高溫下分布的相關(guān)報(bào)道.雖然Takuwa等[10]基于Chemkin軟件模擬了在H2、O2和NaOH混合系統(tǒng)中燃燒時(shí)含鈉氣相產(chǎn)物的轉(zhuǎn)化過程,但并未考慮HCl及SO2的影響. 化學(xué)動(dòng)力學(xué)是研究化學(xué)反應(yīng)機(jī)理的有效途徑,針對(duì)準(zhǔn)東煤燃燒過程中氣相堿金屬的高溫形態(tài)分布,筆者采用化學(xué)動(dòng)力學(xué)軟件Chemkin Pro計(jì)算并分析了溫度、壓力及過量空氣系數(shù)對(duì)含鈉氣相產(chǎn)物的影響,通過分析各基元反應(yīng)了解反應(yīng)產(chǎn)物隨各影響因素變化的原因.

    1 氣相化學(xué)組分的化學(xué)動(dòng)力學(xué)模型

    1.1 煤樣參數(shù)

    以含高堿金屬的準(zhǔn)東煤為計(jì)算對(duì)象,其工業(yè)分析、元素分析及灰分分析如表1和表2所示,其中氯元素質(zhì)量分?jǐn)?shù)的測定采用GB/T 3558—1996 《煤中氯的測定方法》.

    表1 準(zhǔn)東煤的工業(yè)分析和元素分析

    表2 準(zhǔn)東煤灰分分析

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

    基元反應(yīng)的速率常數(shù)k為依賴于溫度的Arrehenius形式:

    k=A·Tβexp[-E/(RT)]

    (1)

    式中:k為速率常數(shù),s-1;A為指前因子;T為溫度,K;β為指數(shù);E為活化能,J/mol;R為氣體常數(shù),取8.314 J/(mol·K).

    計(jì)算模型選用PSR(Perfect stirred reactor),即物料以穩(wěn)定的流量進(jìn)入反應(yīng)器后,瞬間就在整個(gè)反應(yīng)器內(nèi)分散均勻并與反應(yīng)器內(nèi)原存留的物料完全混合,因此反應(yīng)物轉(zhuǎn)化為生成物的速率由化學(xué)反應(yīng)速率控制,而不是由混合過程控制.

    基于涉及C、H、O、N元素的GADM98詳細(xì)化學(xué)反應(yīng)機(jī)理[11],同時(shí)考慮了Glarborg等[12-13]提出的S、Cl、K、Na元素的詳細(xì)化學(xué)反應(yīng)機(jī)理,建立了Na高溫形態(tài)分布的化學(xué)動(dòng)力學(xué)模型.該模型涉及的主要元素和化合物如表3所示,模型中考慮了139步反應(yīng). 計(jì)算中使用的熱力學(xué)數(shù)據(jù)來自于Chemkin自帶的therm.dat文件、NIST(National Institute of Standards and Technology)數(shù)據(jù)庫和NASA(National Aeronautics and Space Administration)發(fā)表的熱力學(xué)數(shù)據(jù)[14].

    表3 化學(xué)動(dòng)力學(xué)模型中包含的組分

    1.3 參數(shù)設(shè)置

    根據(jù)表1和表2中準(zhǔn)東煤的分析數(shù)據(jù),初始組分分2部分輸入計(jì)算,C、H、O、N、S元素折算為理論煙氣含量輸入,Cl、K、Na元素按以下形式輸入計(jì)算.

    煤中的鈉以4種形態(tài)存在:無機(jī)鈉、羧酸鹽形式的有機(jī)鈉、配位形式的有機(jī)鈉和不可溶鈉[15]. 在該煤樣中,上述4種鈉的質(zhì)量分?jǐn)?shù)依次為61%、11%、4%和24%[16]. 劉大海等[15]認(rèn)為以NaCl晶體和水合離子形式存在的無機(jī)鈉都以NaCl的形式釋放,有機(jī)鈉以鈉原子的形式釋放,而不可溶鈉在加熱過程中不揮發(fā),故鈉元素以NaCl(61%w(Natotal))和鈉原子(15%w(Natotal))的形式輸入. 煤中的鉀也以4種形態(tài)存在,其中水溶性鉀、醋酸銨溶鉀、鹽酸溶鉀和不可溶鉀的質(zhì)量分?jǐn)?shù)依次為29%、8%、5%和58%[16],鉀的析出形式與鈉相似[17],故鉀元素以KCl(29%w(Ktotal))和鉀原子(13%w(Ktotal))的形式輸入. 對(duì)于氯,除了以堿金屬氯化物的形式析出外,還會(huì)以HCl的形式析出[18].可通過NaCl和KCl的含量來確定與堿金屬結(jié)合的氯含量(NaCl:43.0%w(Cltotal);KCl:1.6%w(Cltotal)),剩余的即為以HCl(55.4%w(Cltotal))形態(tài)存在的氯,故氯有NaCl、KCl和HCl 3部分輸入.

    此外,該模型中未將硅、鋁、鈣等包括在內(nèi).雖然這些元素對(duì)氣相堿金屬的析出有一定影響,但該化學(xué)動(dòng)力學(xué)計(jì)算主要研究的是氣相堿金屬與不同化合物間的相互轉(zhuǎn)化,且含硅、鋁、鈣礦物對(duì)堿金屬吸附的化學(xué)動(dòng)力學(xué)機(jī)理尚不明確,因此該模型暫時(shí)不考慮以上幾種元素. 計(jì)算中變量的控制如表4所示.

    表4 計(jì)算中的參數(shù)設(shè)置

    1.4 敏感性分析

    敏感性分析是一種能夠深入直觀地分析組分、基元反應(yīng)和反應(yīng)條件等因素對(duì)系統(tǒng)反應(yīng)參數(shù)敏感程度的方法. 總體敏感性分析用于機(jī)理簡化,局部敏感性分析用于研究詳細(xì)機(jī)理的化學(xué)動(dòng)力學(xué)特性. 局部敏感性分析的正交表達(dá)式[19]如下:

    (2)

    2 結(jié)果及討論

    2.1 Chemkin計(jì)算結(jié)果與試驗(yàn)結(jié)果對(duì)比

    由于氣相堿金屬在線測量比較困難,故研究者大多給出了氣相堿金屬的定性或半定量分析.在與Dleschko等[6]試驗(yàn)條件相同工況下(T=1 473 K,壓力p=100 kPa、300 kPa和900 kPa)進(jìn)行Chemkin計(jì)算.為了驗(yàn)證該模型的可靠性,將模型計(jì)算結(jié)果與試驗(yàn)結(jié)果進(jìn)行對(duì)比,如圖1所示,其中左邊縱坐標(biāo)為計(jì)算所得產(chǎn)物的物質(zhì)的量分?jǐn)?shù),右邊縱坐標(biāo)為質(zhì)譜儀得出的峰面積數(shù)據(jù),其數(shù)值大小與產(chǎn)物的物質(zhì)的量分?jǐn)?shù)成正比. 由圖1可知,隨著壓力升高,NaCl、HCl和SO2的物質(zhì)的量分?jǐn)?shù)都減少,計(jì)算結(jié)果與試驗(yàn)結(jié)果相符,驗(yàn)證了該模型的正確性.

    (a) NaCl

    (b) HCl

    (c) SO2

    2.2 含鈉氣相產(chǎn)物的物質(zhì)的量分?jǐn)?shù)與停留時(shí)間的關(guān)系

    圖2給出了含鈉氣相產(chǎn)物的物質(zhì)的量分?jǐn)?shù)隨停留時(shí)間(t)的變化,圖中列出了主要的反應(yīng)產(chǎn)物,忽略了低物質(zhì)的量分?jǐn)?shù)的產(chǎn)物,如NaOH和Na2SO4等. 由圖 2可知,反應(yīng)初期NaCl(g)的物質(zhì)的量分?jǐn)?shù)約為60%,其他幾種鈉化合物的物質(zhì)的量分?jǐn)?shù)很少,這與初始加入的組分有關(guān);隨著停留時(shí)間延長,NaCl(g)的物質(zhì)的量分?jǐn)?shù)繼續(xù)增加,成為主要的產(chǎn)物;當(dāng)停留時(shí)間達(dá)到10-3s時(shí),各產(chǎn)物的物質(zhì)的量分?jǐn)?shù)基本保持不變,在此狀態(tài)下鈉主要以NaCl(g)形態(tài)存在,其物質(zhì)的量分?jǐn)?shù)約為70%.

    圖2 含鈉氣相產(chǎn)物的物質(zhì)的量分?jǐn)?shù)隨停留時(shí)間的變化

    由此可見,在給定的煙氣流量下,反應(yīng)進(jìn)行10-3s以后,各產(chǎn)物物質(zhì)的量分?jǐn)?shù)基本處于穩(wěn)定值,可以認(rèn)為停留時(shí)間達(dá)到10-3s后,反應(yīng)達(dá)到了平衡.

    2.3 含鈉氣相產(chǎn)物的物質(zhì)的量分?jǐn)?shù)與溫度的關(guān)系

    圖3給出了含鈉氣相產(chǎn)物的物質(zhì)的量分?jǐn)?shù)隨溫度的變化.由圖3可知,各產(chǎn)物的物質(zhì)的量分?jǐn)?shù)隨溫度的變化主要經(jīng)歷以下過程:

    (1)溫度低于1 200 K時(shí),產(chǎn)物以NaCl(g)和(NaCl)2(g)為主;在溫度升至1 300 K的過程中,(NaCl)2(g)的物質(zhì)的量分?jǐn)?shù)逐漸減少,而NaCl(g)的物質(zhì)的量分?jǐn)?shù)則逐漸增加,并在1 300 K左右達(dá)到最大值75%,此過程中(NaCl)2(g)逐漸向NaCl(g)轉(zhuǎn)化. 對(duì)于該過程,劉敬等[16]認(rèn)為,從1 000 K左右開始,準(zhǔn)東煤中的鈉主要以NaCl(g)的形式析出,部分以中間體(NaCl)2(g)的形式析出,且(NaCl)2(g)在1 273 K后消失,這與圖3的結(jié)果相近.

    (2)溫度高于1 400 K后,NaCl(g)的物質(zhì)的量分?jǐn)?shù)逐漸減少,在1 900 K時(shí)減至60%左右;與此同時(shí),NaOH(g)和Na(g)的物質(zhì)的量分?jǐn)?shù)逐漸增加,在1 900 K時(shí)分別達(dá)到12%和3%,此時(shí)部分NaCl(g)轉(zhuǎn)化成NaOH(g)和Na(g). 對(duì)于此過程,Bl?sing等[5]認(rèn)為在高溫下特別是還原性氣氛下,NaCl(g)會(huì)分解生成Na(g),但是Na(g)不穩(wěn)定,一部分Na(g)會(huì)與H2O(g)反應(yīng)生成NaOH(g);此外,堿金屬氯化物也會(huì)直接轉(zhuǎn)化成氫氧化物[16].

    (3)由圖3(b)可知, 產(chǎn)物中以硫酸鹽形式存在的鈉的物質(zhì)的量分?jǐn)?shù)在0.3%以下.這主要是由于該煤樣中氯的物質(zhì)的量分?jǐn)?shù)比鈉的物質(zhì)的量分?jǐn)?shù)高(Cl、Na摩爾比為1.4∶1),根據(jù)Kim等[20]的研究結(jié)果可知,氯的物質(zhì)的量分?jǐn)?shù)直接影響NaCl的析出,NaCl的析出量隨氯質(zhì)量分?jǐn)?shù)的增加而增多. 另外,在1 200 K左右,NaHSO4(g)和Na2SO4(g)的物質(zhì)的量分?jǐn)?shù)都達(dá)到了峰值.

    (a) 高物質(zhì)的量分?jǐn)?shù)產(chǎn)物

    (b) 低物質(zhì)的量分?jǐn)?shù)產(chǎn)物

    選取圖3中的NaCl和Na2SO4作為分析對(duì)象,對(duì)NaCl在1 100 K和1 900 K溫度下進(jìn)行敏感性分析,對(duì)Na2SO4在1 200 K溫度下進(jìn)行敏感性分析,結(jié)果如圖4所示.在圖4(a)中,反應(yīng)(3)的敏感性系數(shù)的絕對(duì)值較大,且為正值,這將促進(jìn)NaCl的生成;反應(yīng)(4)對(duì)NaCl的生成起抑制作用,但其敏感性系數(shù)的絕對(duì)值較小,故在1 100 K下,NaCl(g)的生成反應(yīng)占據(jù)主導(dǎo),NaCl(g)的物質(zhì)的量分?jǐn)?shù)增加,這與圖3(a)中NaCl(g)在1 100 K左右的變化規(guī)律相符.

    NaO2+HCl=NaCl+HO2

    (3)

    NaCl+NaCl=(NaCl)2

    (4)

    在圖4(b)中,對(duì)NaCl生成影響較大的是反應(yīng)(4)和反應(yīng)(5). 與反應(yīng)(4)相比,反應(yīng)(6)的敏感性系數(shù)的絕對(duì)值更大,即對(duì)NaCl的抑制作用更強(qiáng).因此,在1 900 K時(shí),NaCl(g)的物質(zhì)的量分?jǐn)?shù)比1 100 K時(shí)少. 根據(jù)反應(yīng)(6),NaCl(g)將轉(zhuǎn)化為NaOH(g),NaCl(g)的物質(zhì)的量分?jǐn)?shù)減少,而NaOH(g)的物質(zhì)的量分?jǐn)?shù)增加,這與圖3(a)的轉(zhuǎn)化規(guī)律相符,也與劉敬等[16]的結(jié)論一致.

    (a) 1 100 K

    (b) 1 900 K

    (c) 1 200 K

    Na+HCl=NaCl+H

    (5)

    NaCl+H2O=NaOH+HCl

    (6)

    在圖4(c)中,主要是反應(yīng)(7)和反應(yīng)(8)促進(jìn)Na2SO4物質(zhì)的量分?jǐn)?shù)增加.在815 ℃下燃燒時(shí),大部分有機(jī)鈉以鈉原子形式釋放[15],此時(shí)生成的鈉接觸到SO2便會(huì)發(fā)生反應(yīng)生成Na2SO4. 這一結(jié)論與Bl?sing等[4]的結(jié)論有所不同,Bl?sing[4]認(rèn)為Na2SO4是由NaCl直接轉(zhuǎn)化而來的,發(fā)生的反應(yīng)見式(9).

    Na+SO2(+M)=NaSO2(+M)

    (7)

    NaSO2+NaO2=Na2SO4

    (8)

    4NaCl+2H2O+2SO2+O2=2Na2SO4+4HCl

    (9)

    2.4 含鈉氣相產(chǎn)物的物質(zhì)的量分?jǐn)?shù)與過量空氣系數(shù)的關(guān)系

    圖5給出了NaCl(g)、NaOH(g)和Na(g)的物質(zhì)的量分?jǐn)?shù)與過量空氣系數(shù)α和溫度的關(guān)系. 由圖5可知,溫度高于1 100 K時(shí),α對(duì)NaCl(g)、NaOH(g)和Na(g)的物質(zhì)的量分?jǐn)?shù)的影響才有所體現(xiàn);當(dāng)α≥1時(shí),即在氧化性氣氛下,Na(g)、NaCl(g)和NaOH(g)的物質(zhì)的量分?jǐn)?shù)的差異不大;當(dāng)α<1時(shí),隨著α增大,NaCl(g)和NaOH(g)的物質(zhì)的量分?jǐn)?shù)逐漸增加,Na(g)的物質(zhì)的量分?jǐn)?shù)逐漸減少. 這與Wei等[9]的熱力學(xué)計(jì)算結(jié)果相近,Wei等[9]認(rèn)為與燃燒過程相比,氣化過程(即α<1)對(duì)堿金屬析出的影響更大,而隨過量空氣系數(shù)α增大,Na(g)的物質(zhì)的量分?jǐn)?shù)減少及NaOH(g)的物質(zhì)的量分?jǐn)?shù)增加的主要原因是發(fā)生了如下反應(yīng):

    (b) NaOH

    (c) Na

    (10)

    由圖4(b)可知,溫度較高時(shí),反應(yīng)(5)會(huì)促進(jìn)NaCl(g)的生成;在氧化性氣氛下,氫原子被快速氧化,致使該反應(yīng)的平衡右移,NaCl(g)的物質(zhì)的量分?jǐn)?shù)增加;而在還原性氣氛下,隨著α減小,該反應(yīng)的作用減弱,NaCl(g)的物質(zhì)的量分?jǐn)?shù)減少,Na(g)的物質(zhì)的量分?jǐn)?shù)增加,這與圖5(a)和圖5(c)的變化規(guī)律相符.

    在1 900 K下NaOH的敏感性分析結(jié)果如圖6所示. NaOH的熔點(diǎn)為323 ℃,沸點(diǎn)為1 388 ℃,加熱至分解時(shí)產(chǎn)生Na2O. 由圖6可知,在1 900 K下,影響NaOH生成的主要是反應(yīng)(6)、反應(yīng)(11)和反應(yīng)(12).高溫下,反應(yīng)(6)促進(jìn)了NaCl的消耗以及NaOH物質(zhì)的量分?jǐn)?shù)的增加,即高溫下NaCl向NaOH轉(zhuǎn)化. 在氧化性氣氛下,O2的存在以及氫的消耗會(huì)促進(jìn)反應(yīng)(11)和反應(yīng)(12)中NaOH和NaO的生成;在還原性氣氛下,NaOH的生成將受到抑制,隨著α減小,NaOH(g)的物質(zhì)的量分?jǐn)?shù)逐漸減少,Na(g)的物質(zhì)的量分?jǐn)?shù)逐漸增加,這與圖5(b)的變化規(guī)律相符.

    Na+H2O=NaOH+H

    (11)

    Na+O2=NaO+O

    (12)

    圖6 1 900 K下NaOH的敏感性分析結(jié)果

    2.5 含鈉氣相產(chǎn)物物質(zhì)的量分?jǐn)?shù)與壓力的關(guān)系

    圖7給出了Na(g)、NaCl(g)和NaOH(g)的物質(zhì)的量分?jǐn)?shù)與壓力的關(guān)系. 由圖7可知,當(dāng)溫度高于1 500 K時(shí),Na(g)和NaOH(g)的物質(zhì)的量分?jǐn)?shù)隨壓力的升高而略有減少;當(dāng)溫度低于1 300 K時(shí),NaCl(g)的物質(zhì)的量分?jǐn)?shù)受壓力的影響較大,并隨壓力的升高而減少;當(dāng)溫度高于1 300 K時(shí),NaCl(g)的物質(zhì)的量分?jǐn)?shù)隨壓力升高而略有增加. Bl?sing等[21]和Oleschko等[6]的試驗(yàn)結(jié)果也表明升高壓力會(huì)減少NaCl(g)的物質(zhì)的量分?jǐn)?shù),Wei等[9]的熱力學(xué)計(jì)算結(jié)果表明壓力升高會(huì)導(dǎo)致NaCl(g)和NaOH(g)的物質(zhì)的量分?jǐn)?shù)減少.

    (a) NaCl

    (b) NaOH

    (c) Na

    總體來看,壓力主要影響較低溫度下NaCl(g)的物質(zhì)的量分?jǐn)?shù),其主要原因在于:溫度低于1 300 K時(shí),NaCl(g)與(NaCl)2(g)相互轉(zhuǎn)化,從化學(xué)反應(yīng)平衡角度來看,升高壓力使反應(yīng)(4)的平衡右移,即NaCl(g)的物質(zhì)的量分?jǐn)?shù)減少.

    氣相堿金屬析出后會(huì)與Al2O3、SiO2等氧化物反應(yīng).研究表明,NaCl和Na2SO4是參與反應(yīng)的主要化合物,主要發(fā)生反應(yīng)(13)~反應(yīng)(16),生成的二硅酸鈉(Na2Si2O5)及鈉長石(Na2Al2Si6O16)的熔點(diǎn)分別為874 ℃和1 100 ℃,這些低溫共熔體易加劇爐內(nèi)沾污、結(jié)渣,因此可根據(jù)上述結(jié)論選取合適的工況,控制氣相堿金屬化合物的含量,從而達(dá)到減少鍋爐沾污、結(jié)渣的目的.

    2NaCl+2SiO2+H2O=Na2Si2O5+2HCl

    (13)

    2NaCl+Al2O3+6SiO2+H2O=Na2Al2Si6O16+2HCl

    (14)

    2Na2SO4+6SiO2=2Na2O·6SiO2+2SO2+O2

    (15)

    Na2SO4+Al2O3·2SiO2=Na2O·Al2O3·2SiO2+SO3

    (16)

    3 結(jié) 論

    (1)對(duì)比試驗(yàn)結(jié)果與化學(xué)動(dòng)力學(xué)模型計(jì)算結(jié)果可知,兩者的變化趨勢相符,驗(yàn)證了該模型的準(zhǔn)確性.

    (2)計(jì)算結(jié)果表明,含鈉氣相產(chǎn)物隨溫度的變化主要經(jīng)歷以下過程:溫度低于1 200 K時(shí),產(chǎn)物以NaCl(g)和(NaCl)2(g)為主;在溫度升至1 300 K過程中,(NaCl)2(g)逐漸轉(zhuǎn)化為NaCl(g);當(dāng)溫度高于1 400 K后,NaCl(g)逐漸轉(zhuǎn)化為NaOH(g)和Na(g). 硫酸鹽形式的含鈉氣相產(chǎn)物的物質(zhì)的量分?jǐn)?shù)很少,NaHSO4(g)和Na2SO4(g)的物質(zhì)的量分?jǐn)?shù)均在1 200 K達(dá)到峰值.

    (3)溫度高于1 100 K時(shí),過量空氣系數(shù)α對(duì)含鈉氣相產(chǎn)物的影響較明顯:當(dāng)α≥1時(shí),Na(g)、NaCl(g)和NaOH(g)的物質(zhì)的量分?jǐn)?shù)基本不受α的影響;當(dāng)α<1時(shí),隨著α增大,Na(g)的物質(zhì)的量分?jǐn)?shù)逐漸減少,NaCl(g)和NaOH(g)的物質(zhì)的量分?jǐn)?shù)逐漸增加.

    (4)壓力對(duì)含鈉氣相產(chǎn)物的影響不大,當(dāng)溫度高于1 500 K時(shí),隨著壓力生高,Na(g)和NaOH(g)的物質(zhì)的量分?jǐn)?shù)有所減少;當(dāng)溫度低于1 300 K時(shí),NaCl(g)的物質(zhì)的量分?jǐn)?shù)隨著壓力升高而減少;當(dāng)溫度高于1 300 K時(shí),NaCl(g)的物質(zhì)的量分?jǐn)?shù)隨著壓力升高而略有增加.

    (5)不同堿金屬化合物的物理性質(zhì)不同,與其他氧化物生成低溫共熔體的能力也不同,故可根據(jù)以上結(jié)論選取合適的工況,控制氣相堿金屬化合物的組成,從而達(dá)到減少鍋爐沾污、結(jié)渣的目的.

    [1] 崔育奎, 張翔, 烏曉江. 配煤對(duì)新疆準(zhǔn)東高堿煤沾污結(jié)渣特性的影響[J].動(dòng)力工程學(xué)報(bào), 2105, 35(5): 361-365.

    CUI Yukui, ZHANG Xiang, WU Xiaojiang. The effect of coal blending ratio on slagging/fouling characteristics of Zhundong Xinjiang high-alkali coal[J].JournalofChineseSocietyofPowerEngineering, 2015, 35(5): 361-365.

    [2] WANG X B, XU Z X, WEI B, et al. The ash deposition mechanism in boilers burning Zhundong coal with high contents of sodium and calcium: A study from ash evaporating to condensing[J].AppliedThermalEngineering, 2015, 80: 150-159.

    [3] 代百乾, 烏曉江, 張忠孝. 高堿煤燃燒過程中灰中主要元素的遷移規(guī)律[J].動(dòng)力工程學(xué)報(bào), 2014, 34(6): 438-442.

    DAI Baiqian, WU Xiaojiang, ZHANG Zhongxiao. Transition behavior of main elements in fly ash during high alkali coal combustion[J].JournalofChineseSocietyofPowerEngineering, 2014, 34(6): 438-442.

    [6] OLESCHKO H, MüLLER M. Influence of coal composition and operating conditions on the release of alkali species during combustion of hard coal[J].Energy&Fuels, 2007, 21(6): 3240-3248.

    [7] 譚力. 生物質(zhì)燃燒過程中堿金屬遷移及結(jié)渣特性的研究[D]. 北京: 華北電力大學(xué), 2014.

    [8] 陳安和, 楊學(xué)民, 林偉剛. 生物質(zhì)熱解和氣化過程中Cl及堿金屬逸出行為的化學(xué)熱力學(xué)平衡分析[J].燃料化學(xué)學(xué)報(bào), 2007, 35(5): 539-547.

    CHEN Anhe, YANG Xuemin, LIN Weigang. Release characteristics of chlorine and alkali metals during pyrolysis and gasification of biomass by thermodynamical equilibrium analysis[J].JournalofFuelChemistryandTechnology, 2007, 35(5): 539-547.

    [9] WEI X L, SCHNELL U, HEIN K R G. Behaviour of gaseous chlorine and alkali metals during biomass thermal utilisation[J].Fuel, 2004, 84(7/8): 841-848.

    [10] TAKUWA T, NARUSE I. Detailed kinetic and control of alkali metal compounds during coal combustion[J].FuelProcessTechnology, 2007, 88(11/12): 1029-1034.

    [11] GLARBORG P, ALZUETA M U, JOHANSEN K D, et al. Kinetic modeling of hydrocarbon/nitric oxide interactions in a flow reactor[J].CombustionandFlame, 1998, 115(1/2): 1-27.

    [12] GLARBORG P, MARSHALL P. Mechanism and modeling of the formation of gaseous alkali sulfates[J].CombustionandFlame, 2005, 141(1/2): 22-39.

    [13] HINDIYARTI L, FRANDSEN F, LIVBJERG H, et al. An exploratory study of alkali sulfate aerosol formation during biomass combustion[J].Fuel, 2008, 87(8/9): 1591-1600.

    [14] MCBRIDE B J, ZEHE M J, GORDON S. NASA glenn coefficients for calculating thermodynamic properties of individual species[R]. Cleveland, Ohio: Glenn Research Center, 2002.

    [15] 劉大海, 張守玉, 涂圣康, 等. 五彩灣煤中鈉在燃燒過程中的遷移釋放規(guī)律[J].化工進(jìn)展, 2015, 24(3): 705-709.

    LIU Dahai, ZHANG Shouyu, TU Shengkang, et al. Modes of occurrence and transformation of alkali metals in Zhundong coal during combustion[J].ChemicalIndustryandEngineeringProgress, 2015, 24(3): 705-709.

    [16] 劉敬, 王智化, 項(xiàng)飛鵬, 等. 準(zhǔn)東煤中堿金屬的賦存形式及其在燃燒過程中的遷移規(guī)律實(shí)驗(yàn)研究[J].燃料化學(xué)學(xué)報(bào), 2014, 42(3): 316-322.

    LIU Jing, WANG Zhihua, XIANG Feipeng, et al. Experimental study on the occurrence of alkali metal of Zhundong coal and its migration during combustion[J].JournalofFuelChemistryandTechnology, 2014, 42(3): 316-322.

    [17] JOHANSEN J M, JAKOBSEN J G, FRANDSEN F J, et al. Release of K, Cl, and S during pyrolysis and combustion of high-chlorine biomass[J].Energy&Fuels, 2011, 25(11): 4961-4971.

    [18] LI G Y, WANG C A, YAN Y, et al. Release and transformation of sodium during combustion of Zhundong coals[J].JournaloftheEnergyInstitute, 2016, 89(1): 48-56.

    [19] 徐瑩. 氣體再燃化學(xué)動(dòng)力學(xué)分析及機(jī)理簡化研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2010.

    [20] KIM S S, KANG Y S, LEE H D, et al. Release of potassium and sodium species during combustion of various rank coals, biomass, sludge and peats[J].JournalofIndustrialandEngineeringChemistry, 2012, 18(6): 2199-2203.

    Chemical Kinetics Simulation on High-temperature Morphology of Alkali Sodium in Zhundong Coal

    JIJieqiang,CHENGLeming,LIUYanquan,WANGQinhui,LUOZhongyang

    (State Key Laboratory of Clean Energy Utilization, Institute for Thermal Power Engineering, Zhejiang University, Hangzhou 310027, China)

    2016-10-17

    2016-11-03

    國家科技支撐計(jì)劃資助項(xiàng)目(2015BAA04B02-09)

    季杰強(qiáng)(1991-),男,浙江金華人,博士研究生,主要從事準(zhǔn)東煤沾污結(jié)渣方面的研究. 程樂鳴(通信作者),男,教授,博士,電話(Tel.):0571-87953208;E-mail:lemingc@zju.edu.cn.

    1674-7607(2017)10-0780-08

    TF536

    A

    470.30

    猜你喜歡
    準(zhǔn)東堿金屬氣相
    氣相過渡金屬鈦-碳鏈團(tuán)簇的研究
    新疆首個(gè)煤炭資源主題能源館開館
    中國煤炭(2019年8期)2019-01-17 10:16:53
    大型工業(yè)園需水預(yù)測方法探討——以新疆準(zhǔn)東經(jīng)濟(jì)技術(shù)開發(fā)區(qū)為例
    新型釩基催化劑催化降解氣相二噁英
    新疆準(zhǔn)東煤田西部八道灣組層序地層及聚煤規(guī)律分析
    預(yù)縮聚反應(yīng)器氣相管“鼓泡”的成因探討
    準(zhǔn)東地區(qū)電網(wǎng)短路電流面臨的形勢及對(duì)策
    《堿金屬》疑難問題解析
    有關(guān)堿金屬的計(jì)算舉隅
    氣相防銹技術(shù)在電器設(shè)備防腐中的應(yīng)用
    一二三四中文在线观看免费高清| 国产成人精品婷婷| 最近中文字幕高清免费大全6| 在现免费观看毛片| 少妇精品久久久久久久| 永久网站在线| 亚州av有码| 亚洲欧美日韩卡通动漫| 人妻系列 视频| 日韩精品免费视频一区二区三区 | 97超碰精品成人国产| 全区人妻精品视频| 极品人妻少妇av视频| 久久97久久精品| 日本黄色片子视频| kizo精华| 大香蕉久久成人网| 亚洲精品久久午夜乱码| 中文字幕免费在线视频6| 精品少妇黑人巨大在线播放| 日本午夜av视频| 黄片播放在线免费| 啦啦啦啦在线视频资源| 国产成人一区二区在线| av免费观看日本| 黑人欧美特级aaaaaa片| 美女视频免费永久观看网站| 国产精品一区二区三区四区免费观看| 一区二区av电影网| 春色校园在线视频观看| 亚洲三级黄色毛片| 日本黄大片高清| 成人国产av品久久久| 十八禁网站网址无遮挡| 亚洲熟女精品中文字幕| 少妇猛男粗大的猛烈进出视频| 国产精品一二三区在线看| 国产av一区二区精品久久| 交换朋友夫妻互换小说| 丰满乱子伦码专区| 日本黄色片子视频| 欧美+日韩+精品| 91国产中文字幕| 97在线视频观看| 国产伦理片在线播放av一区| 午夜日本视频在线| 久久久久人妻精品一区果冻| 亚洲精品第二区| 欧美变态另类bdsm刘玥| 日韩人妻高清精品专区| 国产精品女同一区二区软件| 国产无遮挡羞羞视频在线观看| 免费看不卡的av| 亚洲精品亚洲一区二区| 欧美另类一区| 欧美3d第一页| 亚洲经典国产精华液单| 搡女人真爽免费视频火全软件| 中文字幕av电影在线播放| 如何舔出高潮| 亚洲成人av在线免费| 精品久久久久久电影网| 2021少妇久久久久久久久久久| 乱码一卡2卡4卡精品| 九九在线视频观看精品| 国产精品人妻久久久影院| 插阴视频在线观看视频| 国产淫语在线视频| 汤姆久久久久久久影院中文字幕| 人人妻人人添人人爽欧美一区卜| 蜜桃国产av成人99| 精品酒店卫生间| 69精品国产乱码久久久| 亚洲高清免费不卡视频| 一个人免费看片子| 日韩人妻高清精品专区| 日本爱情动作片www.在线观看| 999精品在线视频| 色5月婷婷丁香| 欧美日韩成人在线一区二区| 免费观看无遮挡的男女| 日本91视频免费播放| 你懂的网址亚洲精品在线观看| 十分钟在线观看高清视频www| 午夜91福利影院| 国产亚洲av片在线观看秒播厂| 熟女电影av网| 亚洲精品国产色婷婷电影| 51国产日韩欧美| 国产午夜精品一二区理论片| 欧美日韩精品成人综合77777| 成人二区视频| 免费高清在线观看日韩| 国产国拍精品亚洲av在线观看| 男人爽女人下面视频在线观看| 久热这里只有精品99| 国产黄色免费在线视频| 黄色毛片三级朝国网站| 夜夜爽夜夜爽视频| 亚洲熟女精品中文字幕| 日本爱情动作片www.在线观看| 亚洲精品国产av蜜桃| 免费高清在线观看日韩| 精品久久蜜臀av无| 久久久精品区二区三区| 成人亚洲欧美一区二区av| 爱豆传媒免费全集在线观看| 欧美日韩av久久| 国产精品国产av在线观看| 国产精品偷伦视频观看了| 亚洲精品久久午夜乱码| av线在线观看网站| 中文天堂在线官网| 亚洲精品av麻豆狂野| 久久精品国产a三级三级三级| 欧美三级亚洲精品| 人人妻人人添人人爽欧美一区卜| av在线app专区| 三级国产精品欧美在线观看| 精品久久国产蜜桃| av卡一久久| 在线免费观看不下载黄p国产| 欧美日韩亚洲高清精品| 亚洲国产毛片av蜜桃av| 日韩视频在线欧美| 国产成人午夜福利电影在线观看| 黄色一级大片看看| 亚洲欧美精品自产自拍| 黄片无遮挡物在线观看| 高清av免费在线| 国产精品熟女久久久久浪| 欧美bdsm另类| 最后的刺客免费高清国语| 亚洲av成人精品一区久久| 午夜福利网站1000一区二区三区| 91在线精品国自产拍蜜月| 亚洲美女搞黄在线观看| 午夜日本视频在线| 一区二区三区免费毛片| 赤兔流量卡办理| 亚洲精品日韩av片在线观看| 婷婷色综合www| av视频免费观看在线观看| 91国产中文字幕| 一区二区三区四区激情视频| 精品亚洲乱码少妇综合久久| 成年美女黄网站色视频大全免费 | 国产乱来视频区| 国产成人精品福利久久| 免费久久久久久久精品成人欧美视频 | 国产探花极品一区二区| 亚洲精品色激情综合| 国产黄色视频一区二区在线观看| .国产精品久久| 2018国产大陆天天弄谢| 精品99又大又爽又粗少妇毛片| 天堂中文最新版在线下载| 国产高清国产精品国产三级| 精品亚洲成国产av| 九九在线视频观看精品| 精品一区二区免费观看| 在线观看美女被高潮喷水网站| 精品一区二区三区视频在线| 三上悠亚av全集在线观看| 麻豆乱淫一区二区| 一级黄片播放器| 久久婷婷青草| av网站免费在线观看视频| 国产亚洲欧美精品永久| 高清欧美精品videossex| 男人添女人高潮全过程视频| 国产精品久久久久久久电影| av免费在线看不卡| 久久久久网色| 丝袜喷水一区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 看非洲黑人一级黄片| 大片免费播放器 马上看| 国产精品蜜桃在线观看| 国产探花极品一区二区| 91aial.com中文字幕在线观看| 校园人妻丝袜中文字幕| 国产乱来视频区| 亚洲人成网站在线播| 中国美白少妇内射xxxbb| .国产精品久久| 国产午夜精品久久久久久一区二区三区| 日韩中文字幕视频在线看片| 在线观看国产h片| 亚洲婷婷狠狠爱综合网| 久久97久久精品| 亚洲综合精品二区| 丰满饥渴人妻一区二区三| 高清视频免费观看一区二区| 国产精品一国产av| xxxhd国产人妻xxx| 精品99又大又爽又粗少妇毛片| 久久久国产精品麻豆| 日本wwww免费看| 王馨瑶露胸无遮挡在线观看| 特大巨黑吊av在线直播| 亚洲人成网站在线观看播放| 久久99一区二区三区| 99热网站在线观看| 少妇猛男粗大的猛烈进出视频| 国产片特级美女逼逼视频| 97在线人人人人妻| videosex国产| 欧美日韩成人在线一区二区| 有码 亚洲区| 日本欧美国产在线视频| 一二三四中文在线观看免费高清| 下体分泌物呈黄色| 精品亚洲成a人片在线观看| 国产av精品麻豆| av福利片在线| 久久国产亚洲av麻豆专区| 熟妇人妻不卡中文字幕| 人妻制服诱惑在线中文字幕| 一区在线观看完整版| 插阴视频在线观看视频| 亚洲中文av在线| 一级片'在线观看视频| 国产精品久久久久久久久免| 18在线观看网站| av国产久精品久网站免费入址| 国产精品国产av在线观看| 国产成人91sexporn| 国产日韩欧美视频二区| 久久久久久久久久久久大奶| 婷婷色综合www| 欧美激情国产日韩精品一区| 97在线视频观看| 亚洲国产精品国产精品| av免费观看日本| 精品人妻熟女av久视频| 免费观看av网站的网址| 一区二区三区免费毛片| 国产成人一区二区在线| 亚洲欧洲精品一区二区精品久久久 | 欧美一级a爱片免费观看看| 免费观看a级毛片全部| 五月天丁香电影| 国产在线免费精品| 久久综合国产亚洲精品| 80岁老熟妇乱子伦牲交| xxx大片免费视频| 狠狠精品人妻久久久久久综合| 久久免费观看电影| 日本欧美国产在线视频| 国产精品一区二区三区四区免费观看| 人人澡人人妻人| 久久精品熟女亚洲av麻豆精品| 26uuu在线亚洲综合色| 99热6这里只有精品| 久久久久久久久久久免费av| av在线播放精品| 成年女人在线观看亚洲视频| 日本vs欧美在线观看视频| 国产日韩欧美视频二区| 成人毛片a级毛片在线播放| 99久久综合免费| 人妻一区二区av| 亚洲精品色激情综合| 丝袜在线中文字幕| 成人手机av| 欧美日韩av久久| 亚洲中文av在线| 免费久久久久久久精品成人欧美视频 | 色婷婷久久久亚洲欧美| 人体艺术视频欧美日本| 热re99久久国产66热| 亚洲无线观看免费| 91久久精品国产一区二区三区| 欧美精品高潮呻吟av久久| 成人国产麻豆网| 丰满饥渴人妻一区二区三| 亚洲精品aⅴ在线观看| 亚洲成人手机| 91精品三级在线观看| 国产精品欧美亚洲77777| 一个人免费看片子| 十分钟在线观看高清视频www| 亚洲av不卡在线观看| 亚洲精品久久久久久婷婷小说| 高清欧美精品videossex| 看十八女毛片水多多多| 最近中文字幕高清免费大全6| 另类亚洲欧美激情| 亚洲欧洲日产国产| 中文字幕久久专区| 自拍欧美九色日韩亚洲蝌蚪91| 伊人久久国产一区二区| 成人国产麻豆网| 欧美亚洲日本最大视频资源| 午夜av观看不卡| 久久久久国产精品人妻一区二区| 男的添女的下面高潮视频| 一区二区三区乱码不卡18| 欧美日本中文国产一区发布| 黄色怎么调成土黄色| 成人午夜精彩视频在线观看| 一级a做视频免费观看| 精品亚洲成国产av| 色吧在线观看| 国产亚洲精品第一综合不卡 | 美女福利国产在线| 在线精品无人区一区二区三| 在线观看美女被高潮喷水网站| 少妇被粗大的猛进出69影院 | 天堂8中文在线网| 国产av一区二区精品久久| 国产男人的电影天堂91| 丁香六月天网| 18禁裸乳无遮挡动漫免费视频| 亚洲精品乱码久久久v下载方式| 欧美精品一区二区免费开放| 热re99久久精品国产66热6| 精品人妻偷拍中文字幕| 欧美日韩av久久| 国产成人精品久久久久久| 国产成人免费无遮挡视频| 伊人久久精品亚洲午夜| 亚洲精品久久午夜乱码| 久久久久国产网址| 亚洲精品av麻豆狂野| 99国产精品免费福利视频| 国产精品一区二区在线不卡| 国产精品不卡视频一区二区| 久久久久人妻精品一区果冻| 久久亚洲国产成人精品v| 国产精品三级大全| 亚洲人成网站在线观看播放| 中文字幕最新亚洲高清| 亚洲国产精品999| 97在线视频观看| 日韩精品免费视频一区二区三区 | 天堂8中文在线网| 在线播放无遮挡| 国国产精品蜜臀av免费| 精品视频人人做人人爽| 国产男女内射视频| av视频免费观看在线观看| 免费久久久久久久精品成人欧美视频 | 99热这里只有精品一区| 国产精品秋霞免费鲁丝片| 三上悠亚av全集在线观看| 观看av在线不卡| 另类精品久久| 少妇被粗大猛烈的视频| 中文字幕制服av| 免费观看a级毛片全部| 另类精品久久| 欧美bdsm另类| 乱人伦中国视频| 国产一区二区在线观看av| 亚洲图色成人| 美女国产视频在线观看| 久久久a久久爽久久v久久| av免费在线看不卡| 久久国产精品男人的天堂亚洲 | 大香蕉97超碰在线| 国产免费一级a男人的天堂| 在线 av 中文字幕| 久久热精品热| 国产精品一二三区在线看| 大香蕉久久成人网| 嘟嘟电影网在线观看| 国产午夜精品久久久久久一区二区三区| 久久久精品区二区三区| 美女主播在线视频| 成年av动漫网址| 少妇人妻 视频| 九九爱精品视频在线观看| 日韩伦理黄色片| 成年av动漫网址| 亚洲av成人精品一二三区| 自线自在国产av| 国产伦精品一区二区三区视频9| 亚洲美女搞黄在线观看| 99久久中文字幕三级久久日本| 观看美女的网站| 国产毛片在线视频| 亚洲欧美日韩卡通动漫| 夫妻性生交免费视频一级片| 国产在视频线精品| 男女免费视频国产| 水蜜桃什么品种好| 蜜臀久久99精品久久宅男| 亚洲五月色婷婷综合| 汤姆久久久久久久影院中文字幕| 婷婷成人精品国产| 搡女人真爽免费视频火全软件| av线在线观看网站| 日韩中文字幕视频在线看片| 精品人妻在线不人妻| 狂野欧美激情性xxxx在线观看| 久久亚洲国产成人精品v| 激情五月婷婷亚洲| 欧美一级a爱片免费观看看| 欧美人与性动交α欧美精品济南到 | 看免费成人av毛片| 国产一区二区三区av在线| 日本欧美国产在线视频| 尾随美女入室| 国产精品人妻久久久影院| 老熟女久久久| 国产精品偷伦视频观看了| 国语对白做爰xxxⅹ性视频网站| 蜜桃在线观看..| 777米奇影视久久| 乱码一卡2卡4卡精品| 草草在线视频免费看| 妹子高潮喷水视频| 九九在线视频观看精品| 在线观看美女被高潮喷水网站| 精品一品国产午夜福利视频| 桃花免费在线播放| 女的被弄到高潮叫床怎么办| 国产免费福利视频在线观看| 日日摸夜夜添夜夜爱| 青春草视频在线免费观看| 五月伊人婷婷丁香| 99热6这里只有精品| 夫妻性生交免费视频一级片| 人人澡人人妻人| 大又大粗又爽又黄少妇毛片口| 国产精品一国产av| a级毛片黄视频| 精品国产一区二区三区久久久樱花| 久久人人爽av亚洲精品天堂| h视频一区二区三区| 日韩三级伦理在线观看| 乱人伦中国视频| 少妇被粗大猛烈的视频| 天堂俺去俺来也www色官网| 国产成人精品无人区| 午夜91福利影院| 欧美激情极品国产一区二区三区 | 一本色道久久久久久精品综合| 成年av动漫网址| 91aial.com中文字幕在线观看| 亚洲人与动物交配视频| 寂寞人妻少妇视频99o| 日韩欧美精品免费久久| 亚洲av欧美aⅴ国产| 青春草国产在线视频| 3wmmmm亚洲av在线观看| 午夜福利视频在线观看免费| 精品久久国产蜜桃| 国产成人freesex在线| av国产久精品久网站免费入址| 寂寞人妻少妇视频99o| 亚洲精品第二区| 婷婷色综合大香蕉| 国产亚洲av片在线观看秒播厂| 一边亲一边摸免费视频| 欧美激情国产日韩精品一区| 午夜老司机福利剧场| 亚洲精品aⅴ在线观看| 国产精品一二三区在线看| 亚洲精品一区蜜桃| 亚洲伊人久久精品综合| 亚洲丝袜综合中文字幕| 国产综合精华液| 成人亚洲精品一区在线观看| 亚洲av免费高清在线观看| a级毛片在线看网站| 999精品在线视频| 成人国产麻豆网| 国产日韩欧美在线精品| 美女xxoo啪啪120秒动态图| 久久久久久久久久久免费av| a级毛片黄视频| av播播在线观看一区| 99久久人妻综合| 久久人人爽人人爽人人片va| 在线精品无人区一区二区三| 插阴视频在线观看视频| freevideosex欧美| 亚洲精品乱码久久久v下载方式| 自拍欧美九色日韩亚洲蝌蚪91| 日韩在线高清观看一区二区三区| 嫩草影院入口| 久久热精品热| 中文字幕免费在线视频6| av又黄又爽大尺度在线免费看| 热99久久久久精品小说推荐| 在线观看免费高清a一片| 满18在线观看网站| 欧美亚洲日本最大视频资源| 各种免费的搞黄视频| 国产午夜精品一二区理论片| 男男h啪啪无遮挡| 午夜福利网站1000一区二区三区| 国产精品久久久久久精品古装| 极品人妻少妇av视频| 秋霞伦理黄片| 亚洲精品av麻豆狂野| 国产成人精品久久久久久| 天堂俺去俺来也www色官网| 日韩制服骚丝袜av| 成人午夜精彩视频在线观看| 久久韩国三级中文字幕| 天天操日日干夜夜撸| 亚洲国产欧美日韩在线播放| 亚洲av福利一区| 亚洲久久久国产精品| 国语对白做爰xxxⅹ性视频网站| 亚洲一区二区三区欧美精品| 成年av动漫网址| 国产精品免费大片| 五月开心婷婷网| 亚洲国产精品成人久久小说| 日本色播在线视频| 亚洲国产毛片av蜜桃av| 欧美日韩视频高清一区二区三区二| 国产有黄有色有爽视频| 国产亚洲一区二区精品| 18禁裸乳无遮挡动漫免费视频| 日韩成人伦理影院| 亚洲国产最新在线播放| 视频在线观看一区二区三区| 九九爱精品视频在线观看| 亚洲欧洲日产国产| 久久久国产精品麻豆| 国产精品99久久99久久久不卡 | 亚洲av成人精品一区久久| 黑人欧美特级aaaaaa片| 中文字幕人妻丝袜制服| 超碰97精品在线观看| 久久久久视频综合| 精品少妇黑人巨大在线播放| 亚洲av免费高清在线观看| 亚洲色图 男人天堂 中文字幕 | 亚洲,欧美,日韩| 我要看黄色一级片免费的| 在线观看免费日韩欧美大片 | 亚洲av电影在线观看一区二区三区| 久久国产精品大桥未久av| 99热6这里只有精品| 美女国产视频在线观看| 午夜福利,免费看| 免费观看的影片在线观看| 欧美日韩在线观看h| 亚洲第一av免费看| 亚洲欧美日韩卡通动漫| 久久久亚洲精品成人影院| 久久精品国产a三级三级三级| 日韩强制内射视频| 插逼视频在线观看| 亚洲av综合色区一区| 精品午夜福利在线看| 人妻制服诱惑在线中文字幕| 国产精品一区www在线观看| 观看av在线不卡| 一级毛片aaaaaa免费看小| 日日爽夜夜爽网站| 亚洲激情五月婷婷啪啪| av在线app专区| 黄片播放在线免费| 高清黄色对白视频在线免费看| 国产欧美日韩一区二区三区在线 | 九草在线视频观看| 午夜激情av网站| www.色视频.com| 午夜福利,免费看| 国产精品嫩草影院av在线观看| 69精品国产乱码久久久| 久热这里只有精品99| 一区二区三区免费毛片| 日本与韩国留学比较| 日本欧美视频一区| 一二三四中文在线观看免费高清| 亚洲精品视频女| 午夜福利视频在线观看免费| 国产黄片视频在线免费观看| 精品人妻偷拍中文字幕| 在线天堂最新版资源| 久久精品久久久久久久性| 欧美精品国产亚洲| 国产伦理片在线播放av一区| 卡戴珊不雅视频在线播放| 成年美女黄网站色视频大全免费 | 国产有黄有色有爽视频| 一本—道久久a久久精品蜜桃钙片| 国产色婷婷99| 有码 亚洲区| 免费大片黄手机在线观看| 国内精品宾馆在线| 久久99精品国语久久久| 精品亚洲乱码少妇综合久久| 一级a做视频免费观看| www.av在线官网国产| 久久久久网色| 国产熟女午夜一区二区三区 | 日韩制服骚丝袜av| 精品久久久噜噜| 久久ye,这里只有精品| 黄色配什么色好看| 亚洲四区av| 久久99蜜桃精品久久| 久久精品国产鲁丝片午夜精品| 我的女老师完整版在线观看| 91精品一卡2卡3卡4卡| 91久久精品国产一区二区成人| 三级国产精品欧美在线观看| 蜜桃国产av成人99| 久久久久视频综合| 91精品三级在线观看| 国产精品嫩草影院av在线观看| 一区二区三区精品91| 在线观看人妻少妇| 一级二级三级毛片免费看| 午夜免费男女啪啪视频观看|