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

    低階煤干燥過(guò)程水分析出動(dòng)力學(xué)行為分析

    2017-12-22 05:47:57張守玉張一帆江鋒浩王才威呂俊復(fù)
    煤炭學(xué)報(bào) 2017年11期
    關(guān)鍵詞:低階級(jí)數(shù)等溫

    李 昊,張守玉,李 尤,慕 晨,張一帆,江鋒浩,王才威,呂俊復(fù)

    (1.上海理工大學(xué) 能源與動(dòng)力工程學(xué)院,上海 200093; 2.清華大學(xué) 熱能工程系,北京 100084)

    低階煤干燥過(guò)程水分析出動(dòng)力學(xué)行為分析

    李 昊1,張守玉1,李 尤1,慕 晨1,張一帆1,江鋒浩1,王才威1,呂俊復(fù)2

    (1.上海理工大學(xué) 能源與動(dòng)力工程學(xué)院,上海 200093; 2.清華大學(xué) 熱能工程系,北京 100084)

    使用熱重分析儀對(duì)伊敏褐煤(YM煤)和天池能源次煙煤(TC煤)進(jìn)行干燥過(guò)程實(shí)驗(yàn)研究,采用等轉(zhuǎn)化率等方法對(duì)低階煤熱重?cái)?shù)據(jù)進(jìn)行分析,得到煤中水分析出活化能E。根據(jù)求得的活化能,由Achar方程和Coats-Redfern方程確定最概然機(jī)理函數(shù),并確定干燥過(guò)程級(jí)數(shù)n和指前因子A。同時(shí),采用差式掃描量熱法(DSC)對(duì)煤中不同賦存形式的水分進(jìn)行定量分析,以此驗(yàn)證熱重分析結(jié)果。結(jié)果表明:低階煤水分析出的活化能E為45~55 kJ/mol,干燥過(guò)程遵循簡(jiǎn)單級(jí)數(shù)模型f(α)=(1-α)n,n≠1。YM煤水分析出活化能為53.28 kJ/mol,干燥過(guò)程級(jí)數(shù)n=1.82~1.95,指前因子A=(3.52~3.96)×108s-1;TC煤水分析出活化能為47.44 kJ/mol,干燥過(guò)程級(jí)數(shù)n=1.39~1.44,指前因子A=(6.62~8.17)×107s-1。DSC的定量分析結(jié)果表明,低階煤中以自由水和弱束縛水為主的不穩(wěn)定組分含量達(dá)到70%~80%。干燥過(guò)程中,自由水先析出,其析出活化能略高于純水蒸發(fā)活化能;強(qiáng)束縛水在干燥末期析出,其析出活化能隨干燥過(guò)程的進(jìn)行急劇升高,且遠(yuǎn)大于純水蒸發(fā)活化能。

    低階煤;干燥動(dòng)力學(xué);活化能;機(jī)理函數(shù)

    我國(guó)低階煤儲(chǔ)量豐富,占煤炭總儲(chǔ)量的50%以上。然而,豐富的孔隙結(jié)構(gòu)和親水含氧基團(tuán)使低階煤含有較多水分,增加了運(yùn)輸成本,降低了燃燒效率,嚴(yán)重限制了低階煤的利用[1-2]。WANG,LIU等[3-4]研究發(fā)現(xiàn)低階煤進(jìn)行干燥預(yù)處理不僅可節(jié)省輔機(jī)能耗,而且有助于提高其整體利用率;劉長(zhǎng)勝等[5]對(duì)非蒸發(fā)式脫水干燥進(jìn)行總結(jié)得到:干燥過(guò)程,煤質(zhì)提高,熱值增加,而且有利于煤工藝性的增強(qiáng)。此外,降低低階煤含水量(≤10%)還有利于氣化、液化等利用過(guò)程[6]。綜上所述,干燥脫水提質(zhì)不僅提高了煤階煤質(zhì),而且有利于拓寬低階煤的利用途徑,對(duì)于提高低階煤的利用價(jià)值和市場(chǎng)競(jìng)爭(zhēng)力,具有十分重要的意義。

    干燥動(dòng)力學(xué)是低階煤干燥技術(shù)開發(fā)的理論基礎(chǔ),其中動(dòng)力學(xué)模型及參數(shù)的確定是描述干燥過(guò)程的核心。常用的研究方法包括等溫法和非等溫法。等溫法理論成熟,在低階煤干燥方面已有較多應(yīng)用。ZHENG等[7]對(duì)褐煤進(jìn)行高溫干燥,通過(guò)等溫法得出水分有效擴(kuò)散系數(shù)及水分析出活化能;KOMATSU[8]、楊亞利等[9]通過(guò)實(shí)驗(yàn)分別建立了褐煤干燥動(dòng)力學(xué)模型,并求得褐煤干燥的界面蒸發(fā)活化能E和指前因子A。相比于等溫法,非等溫法理論體系形成較晚,但是其實(shí)驗(yàn)樣品用量少,檢測(cè)誤差小且分析速度快。因此,逐漸成為熱分析動(dòng)力學(xué)的核心。滕英躍[10]、趙鋒鋒[6]等采用非等溫法研究了褐煤低溫干燥過(guò)程水分的變化,得出水分析出活化能隨干燥過(guò)程變化的大致規(guī)律。綜上所述,國(guó)內(nèi)外主要使用等溫法研究低階煤的干燥過(guò)程,非等溫法使用較少。而且,現(xiàn)有的研究大多集中于活化能等動(dòng)力學(xué)參數(shù)的求解,缺乏機(jī)理函數(shù)等較深層次的理論研究。

    本文使用熱重分析儀分析水分析出過(guò)程,采用等轉(zhuǎn)化率法求得干燥活化能。根據(jù)求得的活化能,采用微分法(Achar方程)和積分法(Coats-Redfern方程)相結(jié)合的方式,確定最概然機(jī)理函數(shù),進(jìn)而求解干燥過(guò)程級(jí)數(shù)n和指前因子A。最后,利用差示掃描量熱法(DSC)進(jìn)行定量分析并驗(yàn)證熱重分析結(jié)果。

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

    1.1 實(shí)驗(yàn)樣品制備

    選取內(nèi)蒙伊敏褐煤(YM)和新疆準(zhǔn)東天池能源次煙煤(TC)作為研究對(duì)象,工業(yè)分析和元素分析結(jié)果見(jiàn)表1。按照GB474—2008煤樣的制備方法,將原始樣品破碎并研磨,用標(biāo)準(zhǔn)篩篩分出100目的合格煤樣,裝入樣品袋中,密封保存。

    表1樣品的工業(yè)分析和元素分析
    Table1Proximateandultimateanalysesofthecoalsamples

    樣品工業(yè)分析/%MarMadAadVadFCad元素分析/%CadHadNadSadOadYM36.707.0813.5036.1943.2356.963.610.950.5917.31TC24.005.145.3327.6661.8766.813.881.100.6917.05

    參考GB/T 4632—2008對(duì)兩種煤樣進(jìn)行飽水實(shí)驗(yàn),具體步驟如下:取上述樣品2 g,放入敞口稱量瓶中鋪平。將稱量瓶置于恒溫恒濕室中(溫度30 ℃,相對(duì)濕度97%),通入高純氮?dú)猓3? d左右。待樣品水分恒定,取出樣品,稱量并裝入樣品袋中密封保存,命名為YM煤和TC煤,作為后續(xù)研究待分析煤樣。

    1.2 實(shí)驗(yàn)裝置和步驟

    使用美國(guó)TA公司生產(chǎn)的Q500熱重儀對(duì)樣品進(jìn)行實(shí)驗(yàn)。稱取10 mg的樣品置于Al2O3坩堝內(nèi),放入熱重儀中。根據(jù)實(shí)驗(yàn)要求,利用程序控溫系統(tǒng)設(shè)定4種升溫速率(5,10,20,30 ℃/min)進(jìn)行實(shí)驗(yàn)。熱重系統(tǒng)會(huì)實(shí)時(shí)記錄樣品質(zhì)量變化,并可直接導(dǎo)出相關(guān)數(shù)據(jù),以備后期使用。

    1.3 分析方法

    使用德國(guó)耐馳公司生產(chǎn)的DSC200型差示掃描量熱儀,配置液氮制冷裝置,對(duì)煤中水分進(jìn)行分析。實(shí)驗(yàn)過(guò)程:取5 mg的飽水樣品裝入鋁坩堝中鋪平,置于氮?dú)?10 mL/min)保護(hù)的量熱儀中,以8 ℃/min的速率將樣品從20 ℃降至-150 ℃。利用凝固放熱特性,DSC能在線測(cè)量出煤樣中水分冷凝放出熱量,并繪制DSC熱分析圖。根據(jù)放熱峰,對(duì)煤樣中水分進(jìn)行定性和定量分析。為方便對(duì)比,實(shí)驗(yàn)中使用純水作為參照樣品。

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

    2.1 低階煤水分析出過(guò)程

    利用熱重分析儀對(duì)不同升溫速率(5,10,20,30 ℃/min)下兩種低階煤水分析出過(guò)程進(jìn)行動(dòng)力學(xué)分析,結(jié)果如圖1所示。由圖1可知,不同升溫速率下的TG,DTG曲線變化規(guī)律相似。低階煤中水分的析出主要分為兩個(gè)階段:初始階段(加速干燥段)失重率與溫度近似呈線性關(guān)系,隨溫度的升高失重率急劇增加。該階段主要是煤中自由水析出階段;第2階段(降速干燥段)失重率隨溫度的升高而下降。這一階段煤中束縛水開始析出,而束縛水與孔隙結(jié)構(gòu)和表面含氧官能團(tuán)之間有較強(qiáng)的煤水相互作用,因此析出速率降低。

    圖1 不同升溫速率下兩種煤的TG與DTG曲線Fig.1 TG and DTG curves of YM and TC coal during drying process at different heating rates

    2.2 低階煤水分析出動(dòng)力學(xué)分析

    2.2.1水分析出動(dòng)力學(xué)模型

    研究發(fā)現(xiàn),如果將干燥看成是濕物質(zhì)變?yōu)楦晌镔|(zhì)和水分蒸發(fā)的過(guò)程,那么干燥過(guò)程與熱分解反應(yīng)相似,可定性為一種固體物質(zhì)分成兩種物質(zhì)的過(guò)程[11]。因此,低階煤的干燥過(guò)程可表示如下:

    式中,Asolid表示含水煤;Bsolid表示干燥煤;Cgas表示水蒸氣。

    前人研究表明相變動(dòng)力學(xué)理論、熱分析動(dòng)力學(xué)理論與薄層干燥模型具有相似性,可將其中的一些方法應(yīng)用于干燥動(dòng)力學(xué)中。因此,類比分解反應(yīng),本文對(duì)低階煤干燥過(guò)程動(dòng)力學(xué)行為進(jìn)行分析。以文獻(xiàn)[12-13]研究為基礎(chǔ),得到干燥過(guò)程動(dòng)力學(xué)方程的微分和積分形式:

    式中,α為煤中水分的轉(zhuǎn)化率,α=(m0-mt)/(m0-m∞)(m0,mt,m∞分別為煤樣初始、t時(shí)刻以及干燥結(jié)束時(shí)的質(zhì)量);f(α)為機(jī)理函數(shù);A為指前因子;E為活化能,kJ/mol;R為理想氣體常數(shù);T為溫度,℃;β為升溫速率,℃/min;p(u)為溫度積分。

    以動(dòng)力學(xué)微分和積分方程為基礎(chǔ),采用非等溫法,可計(jì)算煤干燥過(guò)程中水分析出的活化能E。利用求得的活化能檢測(cè)不同機(jī)理函數(shù)的單升溫速率法,從而確定過(guò)程機(jī)理函數(shù)f(α),并通過(guò)迭代法確定級(jí)數(shù)n及指前因子A的范圍。

    2.2.2水分析出動(dòng)力學(xué)參數(shù)

    干燥動(dòng)力學(xué)研究分為等溫法和非等溫法,非等溫法又包括單升溫速率法和多重升溫速率法(等轉(zhuǎn)化率法)。相比單升溫速率法,多重升溫速率法不涉及動(dòng)力學(xué)模型函數(shù)即可獲得較可靠的活化能值E。因此,本文選擇多重升溫速率法進(jìn)行研究。采用經(jīng)典的FWO 法、KSA 法進(jìn)行求解[14-15]。在式(2)和(3)的基礎(chǔ)上,引入p(u)的近似解,得到FWO 法和KSA 法的計(jì)算公式。

    FWO公式:

    KAS公式:

    由于FWO和KAS公式在推導(dǎo)過(guò)程中采用了積分近似解,所以結(jié)果與真實(shí)值仍存在一定的差距。為了獲得更加精確的動(dòng)力學(xué)參數(shù),以FWO和KAS公式為基礎(chǔ),參考前人的研究,采用迭代法進(jìn)行求解[14-15],基本公式如下:

    其中,h(u)為p(u)真實(shí)值和近似值之比。

    迭代步驟如下:

    (1)先假設(shè)h(u)=1,可以求出活化能E1;

    (2)將某一轉(zhuǎn)化率α下不同升溫速率β的各T值以及E1代入u=E/RT,再由Seman-Yang[15]公式求出h(u),將ln[β/h(u)T2]對(duì)1/T作圖并線性擬合,由斜率可得該轉(zhuǎn)化率α?xí)r的活化能E2;

    (3)以E2代替E1,重復(fù)步驟(2),直至Ei-Ei-1<0.01 kJ/mol,此時(shí)的Ei即可認(rèn)為是E的真實(shí)值。采用同樣方法可求出不同轉(zhuǎn)化率時(shí)E的真實(shí)值。

    表2為3種方法求得的不同轉(zhuǎn)化率下水分析出活化能。隨轉(zhuǎn)化率增大,水分析出活化能先降低后增加,其中α在0.2~0.8時(shí),為水分析出的主要階段,且析出活化能較穩(wěn)定。因此,將α為0.2~0.8對(duì)應(yīng)的平均值作為煤中水分析出活化能。

    表2不同轉(zhuǎn)化率下水分析出活化能
    Table2Activationenergyofmoisturetransferenceduringthecoaldryingprocessatdifferentconversionrates

    αE/(kJ·mol-1)YM褐煤FWO法KAS法迭代法TC次煙煤FWO法KAS法迭代法0.161.9260.3160.0164.9162.3262.040.258.5255.6355.9956.5253.4953.120.357.0154.2753.8753.4050.2849.850.455.7152.7952.3551.3948.0947.610.555.2052.1451.6850.0946.5746.040.654.8751.6751.1849.2045.5244.940.755.5752.2851.7748.8245.3044.690.859.8256.5856.1050.3146.4445.810.986.3884.1983.8760.7957.3556.85平均56.6753.6253.2851.3947.9647.44

    結(jié)合上文對(duì)水分析出過(guò)程的分析和活化能的求解,借鑒前人的研究方法[15],選用Achar法和Coats-Redfern法,采用微分法和積分法結(jié)合的方式確定低階煤干燥過(guò)程中水分析出機(jī)理。

    Achar方程:

    Coats-Redfern方程:

    根據(jù)前人的研究[15],選取以下幾種典型的機(jī)理模型:隨機(jī)成核和生長(zhǎng)(A1)、相邊界反應(yīng)(B2)、二維擴(kuò)散(D2)、三維擴(kuò)散(D3)、簡(jiǎn)單級(jí)數(shù)反應(yīng)1和2(O1和O2)、冪函數(shù)法則(P2)。以熱重結(jié)果為基礎(chǔ),將式(7),(8)等號(hào)左邊的值對(duì)1/T線性回歸,由斜率和截距分別求出E和A。研究發(fā)現(xiàn),不同升溫速率下低階煤的干燥實(shí)驗(yàn)結(jié)果相似。所以僅以β= 30 ℃/min為例,表3和4分別列出了由不同機(jī)理函數(shù)求得YM與TC煤干燥過(guò)程中的活化能E、指前因子A以及線性相關(guān)系數(shù)。

    表3根據(jù)常用機(jī)理函數(shù)求得YM煤干燥動(dòng)力學(xué)參數(shù)(β=30℃/min)
    Table3DryingkineticparametersofYMcoalobtainedbythecommonmechanismfunction(β=30℃/min)

    代號(hào)Achar方程E/(kJ·mol-1)A/s-1R2Coats-Redfern方程E/(kJ·mol-1)A/s-1R2A170.063.02×10-90.98969.475.60×1040.1025B213.956.66×1020.732931.931.30×1040.9761D243.628.28×1060.889765.282.94×1080.9730D359.674.83×1080.965474.711.99×1090.9860O129.753.21×1050.978839.272.06×1050.9927O261.361.86×10100.997257.301.53×1080.9992P217.812.41×1020.78849.743.920.9063

    表4根據(jù)常用機(jī)理函數(shù)求得TC煤干燥動(dòng)力學(xué)參數(shù)(β=30℃/min)
    Table4DryingkineticparametersofTCcoalobtainedbythecommonmechanismfunction(β=30℃/min)

    代號(hào)Achar方程E/(kJ·mol-1)A/s-1R2Coats-Redfern方程E/(kJ·mol-1)A/s-1R2A167.186.07×10-90.980161.27-1.28×10-90.9841B218.333.06×1030.828532.802.07×1040.9854D248.695.37×1070.917566.876.75×1080.9822D364.943.55×1090.976576.424.94×1090.9923O134.321.67×1060.987140.223.46×1050.9972O266.311.24×10110.994658.462.95×1080.9966P213.968.17×1020.642510.185.030.9337

    最概然機(jī)理函數(shù)的確定需要滿足以下幾個(gè)條件:① 較大的線性回歸相關(guān)系數(shù);② 積分法與微分法所求出的動(dòng)力學(xué)參數(shù)相近;③ 求得的E值與等轉(zhuǎn)化率法求出的E值相近。由表3可知,D3,O1和O2的微分、積分方程式線性擬合度較高,但是僅O2的兩種方法求得的活化能相近且與多升溫速率法值吻合。表4中的情況與表3相似。由此推斷低階煤干燥過(guò)程中水分析出遵循簡(jiǎn)單級(jí)數(shù)模型,其一般表達(dá)式如下。

    微分形式:

    積分形式:

    整理式(3)和(10)得

    由迭代法求得的活化能介于O1,O2模型之間,推斷實(shí)際的n值在1~2。在1~2之間以0.001為步長(zhǎng)取n值,對(duì)式(11)等號(hào)左端[(1-α)1-n-1]/(n-1)和右端T2exp(-E/RT)做線性回歸,相關(guān)系數(shù)最接近1時(shí)的n值即為所求,此時(shí)的斜率即為A值。借助Matlab中的循環(huán)及線性回歸算法,求得不同β下兩種煤的干燥過(guò)程級(jí)數(shù)n及指前因子A,見(jiàn)表5。由表5可得,YM煤的干燥過(guò)程級(jí)數(shù)n=1.82~1.95,指前因子A=(3.52~3.96)×108s-1;TC煤的干燥過(guò)程級(jí)數(shù)n=1.39~1.44,指前因子A=(6.62~8.17)×107s-1。

    表5YM煤和TC煤的干燥過(guò)程級(jí)數(shù)n及指前因子A
    Table5Processordernandpre-exponentialfactorAofYMcoalandTCcoaldryingprocess

    升溫速率/(℃·min-1)YM煤nA/s-1R2TC煤nA/s-1R251.953.40×1080.99961.397.15×1070.9956101.923.80×1080.99871.407.48×1070.9914201.863.96×1080.99621.448.17×1070.9900301.823.52×1080.99311.436.62×1070.9878

    2.3 低階煤水分析出特性分析

    2.3.1水分析出特性

    由2.2節(jié)表2可初步得出干燥過(guò)程水分析出活化能變化規(guī)律,為精確反映水分析出率和活化能的關(guān)系,增加實(shí)驗(yàn)數(shù)據(jù)點(diǎn)作圖2(a),并添加空白對(duì)照組——純水蒸發(fā)活化能,如圖2(b)所示。

    圖2 干燥過(guò)程活化能隨轉(zhuǎn)化率的變化Fig.2 Changes in the activation energy with conversion during the coal drying process

    由圖2(a)可知,兩種煤干燥過(guò)程活化能變化分為3個(gè)階段:干燥初期,活化能稍高,隨反應(yīng)的進(jìn)行逐漸降低;平穩(wěn)階段,活化能基本維持不變;干燥末期,活化能急劇升高。而圖2(b)顯示,純水蒸發(fā)過(guò)程僅分為初期和平穩(wěn)兩階段。對(duì)比圖2(a),(b),做出以下推測(cè)及分析[16-17]:

    (1)干燥初期,水分蒸發(fā)活化能下降與煤的理化結(jié)構(gòu)無(wú)關(guān),可能由于溫度測(cè)點(diǎn)在坩堝處,具有熱延遲性,不能反映干燥初期樣品真實(shí)溫度所致。

    (2)平穩(wěn)階段,煤中水分析出活化能略高于純水蒸發(fā)活化能,主要因?yàn)榈碗A煤內(nèi)部孔隙結(jié)構(gòu)。水分蒸發(fā)時(shí),低階煤孔隙中的水分子需克服較大表面張力才能擴(kuò)散到氣相中。此外,孔隙結(jié)構(gòu)給水蒸氣分子的擴(kuò)散帶來(lái)阻力,因此,相比純水的蒸發(fā),低階煤中的水分析出活化能亦升高。

    (3)干燥末期,煤中水分析出活化能急劇升高,而純水的析出活化能幾乎不變,說(shuō)明此階段與低階煤理化結(jié)構(gòu)相關(guān)。究其原因,干燥末期的水分多屬于單層水及少量多層水,極易與含氧官能團(tuán)形成氫鍵,水分子析出要克服氫鍵的束縛[16]。此外,此時(shí)的水分多存在于毛細(xì)孔中,隨干燥的進(jìn)行水分賦存的孔隙孔徑減小,水分子通過(guò)兩相界面需要克服的表面張力極大。綜上兩點(diǎn),末期階段為煤水相互作用較強(qiáng)的水分析出階段。

    2.3.2DSC分析

    以上利用水分蒸發(fā)特性采用熱重法從轉(zhuǎn)化率與活化能的關(guān)系角度對(duì)煤中水分析出特性進(jìn)行了研究。本部分采用差式掃描量熱法,從水分凝固放熱特性的角度對(duì)兩種煤進(jìn)行分析,并對(duì)上述研究結(jié)果進(jìn)行驗(yàn)證。

    參考前人的研究,DSC測(cè)得煤中可凝固水主要為自由水和部分束縛水,通過(guò)定量分析可得到兩種可凝固水含量,但是兩種水分含量之和與干燥法測(cè)得水分總量不同,因此,煤中存在不能凝固的水,即不凍水,可由差減法測(cè)得其含量[18-20]。表6為計(jì)算得到的低階煤中3種水分含量。

    表6YM原煤和TC原煤不同水分含量
    Table6ContentofthewaterwithdifferentexistenceforminYMcoalandTCcoal

    煤樣自由水/%可凝固束縛水/%不凍水/%YM原煤0.4770.1970.326TC原煤0.4860.2910.223

    由表6可知,YM煤中自由水和可凝固束縛水的總含量為67.4%,TC煤中為77.7%。結(jié)合表2和圖2可知,當(dāng)水分析出率達(dá)到70%~80%,水分析出活化能也急劇增大。DSC的定量分析很好地驗(yàn)證了熱重分析結(jié)果,表明無(wú)論以何種方式對(duì)煤中水分進(jìn)行檢測(cè)和分類,從水分的穩(wěn)定性這一本質(zhì)特性出發(fā),低階煤中自由水和弱束縛水穩(wěn)定性差,含量達(dá)到70%~80%,其余20%~30%的強(qiáng)束縛水主要以單分子層及多分子層狀態(tài)存在,在孔隙和分子力的作用下存在煤中孔隙內(nèi)難以脫除。

    3 結(jié) 論

    (1)由多升溫速率法求得低階煤中水分析出活化能為45~55 kJ/mol,其中以迭代法求得的精度最高。

    (2)低階煤中水分析出過(guò)程符合簡(jiǎn)單級(jí)數(shù)反應(yīng)模型f(α)=(1-α)n,n≠1。求得YM褐煤的干燥過(guò)程級(jí)數(shù)n=1.82~1.95,指前因子A=(3.52~3.96)×108s-1;TC準(zhǔn)東煤的干燥過(guò)程級(jí)數(shù)n=1.39~1.44,指前因子A=(6.62~8.17)×107s-1。

    (3)低階煤中自由水和弱束縛水屬于不穩(wěn)定組分,含量達(dá)到70%~80%,剩余20%~30%的水分主要以單分子層及多分子層狀態(tài)存在。干燥過(guò)程中,自由水先析出,其析出活化能略高于純水蒸發(fā)活化能;強(qiáng)束縛水在干燥末期析出,其析出活化能遠(yuǎn)大于純水蒸發(fā)活化能,并隨干燥過(guò)程的進(jìn)行,水分析出活化能急劇升高。

    [1] 趙孟浩,張守玉,鄭紅俊,等.低階煤中含氧官能團(tuán)干燥前后的演變規(guī)律[J].煤炭學(xué)報(bào),2016,41(2):483-489.

    ZHAO Menghao,ZHANG Shouyu,ZHENG Hongjun,et al.Transition of the oxygen-containing functional groups of low rank coal during drying process[J].Journal of China Coal Society,2016,41(2):483-489.

    [2] 王勁草,陽(yáng)金,許鋒.中低溫?zé)峤獾碗A煙煤水焦?jié){的制備[J].黑龍江科技大學(xué)學(xué)報(bào),2014,24(2):115-118.

    WANG Jincao,YANG Jin,XU Feng.Preparation of char water slurry from low rank bituminous coal pyrolyzed at medium-low temperature[J].Journal of Heilongjiang University of Science and Technology,2014,24(2):115-118.

    [3] WANG W C.Laboratory investigation of drying process of Illinois coals[J].Powder Technology,2012,225:72-85.

    [4] LIU M,YAN J,CHONG D,et al.Thermodynamic analysis of pre-drying methods for pre-dried lignite-fired power plant[J].Energy,2013,49(1):107-118.

    [5] 劉長(zhǎng)勝,劉輝.煤非蒸發(fā)干燥脫水技術(shù)述評(píng)[J].煤炭加工與綜合利用,2016(6):41-44.

    [6] 趙鋒鋒.錫林浩特褐煤干燥動(dòng)力學(xué)及水分復(fù)吸特性的研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2014.

    ZHAO Fengfeng.Experimental study Xilinhot lignite drying kinetics and moisture readsorption characteristics[D].Harbin:Harbin Institute of Technology,2014.

    [7] ZHENG H,ZHANG S,GUO X,et al.An experimental study on the drying kinetics of lignite in high temperature nitrogen atmosphere[J].Fuel Processing Technology,2014,126(10):259-265.

    [8] KOMATSU Y,SCIAZKO A,ZAKRZEWSKI M,et al.An experimental investigation on the drying kinetics of a single coarse particle of Belchatow lignite in an atmospheric superheated steam condition[J].Fuel Processing Technology,2015,131:356-369.

    [9] 楊亞利,王志超,楊忠燦,等.褐煤等溫干燥過(guò)程及動(dòng)力學(xué)研究[J].潔凈煤技術(shù),2017,23(1):19-23.

    YANG Yali,WANG Zhichao,YANG Zhongcan,et al.Study on isothermal drying process and kinetics of lignite[J].Clean Coal Technology,2017,23(1):19-23.

    [10] 滕英躍,廉士俊,宋銀敏,等.基于1H-NMR的勝利褐煤原位低溫干燥過(guò)程中弛豫時(shí)間及孔結(jié)構(gòu)變化[J].煤炭學(xué)報(bào),2014,39(12):2525-2530.

    TENG Yingyue,LIAN Shijun,SONG Yinmin,et al.Relaxation time and pore structure evolution of Shengli lignite during in situ low temperature drying using1H-NMR[J].Journal of China Coal Society,2014,39(12):2525-2530.

    [11] 王寶和.相變動(dòng)力學(xué)和熱分析動(dòng)力學(xué)在干燥動(dòng)力學(xué)研究中應(yīng)用的可能性[J].干燥技術(shù)與設(shè)備,2009(3):103-110.

    WANG Baohe.Applicability probability of phase transformation kinetics and thermal analysis kinetics in the drying kinetics[J].Drying Technology & Equipment,2009(3):103-110.

    [12] KANG Taejin,HUEON Namkung,XU Lihua,et al.The drying kinetics of Indonesian low rank coal(IBC) using a lab scale fixed-bed reactor and thermobalance to apply catalytic gasification process[J].Renewable Energy,2013,54(6):138-143.

    [13] JIN Lijun,LI Yang,LIN Lin,et al.Drying characteristic and kinetics of Huolinhe lignite in nitrogen and methane atmospheres[J].Fuel,2015,152:80-87.

    [14] MEHMOOD Muhammad Aamer,YE Guangbin,LUO Huibo,et al.Pyrolysis and kinetic analyses of Camel grass(Cymbopogon schoenanthus) for bioenergy[J].Bioresource Technology,2017,228:18-24.

    [15] 趙學(xué)偉,魏益民,張波.熱重分析法研究小米擠壓膨化產(chǎn)品的干燥動(dòng)力學(xué)[J].鄭州輕工業(yè)學(xué)院學(xué)報(bào)(自然科學(xué)版),2011,26(4):17-21.

    ZHAO Xuewei,WEI Yimin,ZHANG Bo.Study on drying kinetics of millet extrudate by thermogravimetry[J].Journal of Zhengzhou University of Light Industry(Natural Science),2011,26(4):17-21.

    [16] MRAW S C,NAAS-O′ROURKE D F.Water in coal pores:low-temperature heat capacity behavior of the moisture in wyodak coal[J].Science,1979,205:901-902.

    [17] LIN Baiquan,LI He,YUAN Desheng,et al.Development and application of an efficient gas extraction model for low-rank high-gas coal beds[J].International Journal of Coal Science & Technology,2015,2(1):76-83.

    [18] NORINAGA K,KUMAGAI H,HAYASHI J,et al.Classification of water sorbed in coal on the basis of congelation characteristics[J].Energy & Fuels,1998,12(3):574-579.

    [19] FEI Y,CHAFFEE A L,MARSHALL M,et al.Lignite-water interactions studied by phase transition-differential scanning calorimetry[J].Fuel,2005,84(12-13):1557-1562.

    [20] TAHMASEBI A,YU J,SU H,et al.A differential scanning calorimetric(DSC) study on the characteristics and behavior of water in low-rank coals[J].Fuel,2014,135(11):243-252.

    Analysisofthekineticcharacteristicofwateremissionduringlow-rankcoaldryingprocess

    LI Hao1,ZHANG Shouyu1,LI You1,MU Chen1,ZHANG Yifan1,JIANG Fenghao1,WANG Caiwei1,Lü Junfu2

    (1.SchoolofEnergyandPowerEngineering,UniversityofShanghaiforScienceandTechnology,Shanghai200093,China; 2.DepartmentofThermalEngineering,TsinghuaUniversity,Beijing100084,China)

    The drying processes of Yimin lignite (YM) and Tianchi subbituminous (TC) were studied using thermogravimetric analysis.Isoconversion method was used to analyze the TG data to obtain the drying activation energyE.Based on the activation energy obtained,the most probable mechanism function was determined by Achar equation and Coats-Redfern equation,and then the reaction ordernand the pre-exponential factorAwere obtained.At the same time,the differential scanning calorimetry (DSC) method was used to quantitatively analyze the moisture contained in the coals and verified the thermogravimetric results.The experimental results show that the activation energyEof the water emission process is 45-55 kJ/mol and the drying process follows the simple series model for the low-rank coals.For YM coal,E=53.28 kJ/mol,n=1.82-1.95,A=(3.52-3.96)×108s-1and for TC coal,E=47.44 kJ/mol,n=1.39-1.44,A=(6.62-8.17)×107s-1.The quantitative analysis on the water existence form using DSC method indicates that the content of unstable water mainly in the form of free water and weak bound water is 70%-80%.During the drying process,the free water first releases from the coal particle and the strong bound water releases at the end of the period.The activation energy of the strong bound water emission increases abruptly with the drying process and is much larger than that of pure water.

    low-rank coal;drying dynamics;activation energy;mechanism function

    李昊,張守玉,李尤,等.低階煤干燥過(guò)程水分析出動(dòng)力學(xué)行為分析[J].煤炭學(xué)報(bào),2017,42(11):3014-3020.

    10.13225/j.cnki.jccs.2017.0611

    LI Hao,ZHANG Shouyu,LI You,et al.Analysis of the kinetic characteristic of water emission during low-rank coal drying process[J].Journal of China Coal Society,2017,42(11):3014-3020.doi:10.13225/j.cnki.jccs.2017.0611

    TQ536.1

    A

    0253-9993(2017)11-3014-07

    2017-05-05

    2017-07-11責(zé)任編輯許書閣

    “十二五”國(guó)家科技支撐計(jì)劃資助項(xiàng)目(2012BAA04B01);上海理工大學(xué)科技發(fā)展資助項(xiàng)目(16KJFZ094)

    李 昊(1994—),男,江蘇南通人,碩士研究生。E-mail:haoli_usst@163.com。

    張守玉(1971—),男,吉林集安人,教授,博士生導(dǎo)師。E-mail:zhangsy-guo@163.com

    猜你喜歡
    低階級(jí)數(shù)等溫
    EPDM/PP基TPV非等溫結(jié)晶行為的研究
    山西低階煤分布特征分析和開發(fā)利用前景
    一類具低階項(xiàng)和退化強(qiáng)制的橢圓方程的有界弱解
    Dirichlet級(jí)數(shù)及其Dirichlet-Hadamard乘積的增長(zhǎng)性
    Extended Fisher-Kolmogorov方程的一類低階非協(xié)調(diào)混合有限元方法
    幾個(gè)常數(shù)項(xiàng)級(jí)數(shù)的和
    p級(jí)數(shù)求和的兩種方法
    快速檢測(cè)豬鏈球菌的環(huán)介導(dǎo)等溫?cái)U(kuò)增方法
    納米CaCO3對(duì)FEP非等溫結(jié)晶動(dòng)力學(xué)的影響
    國(guó)內(nèi)外低階煤煤層氣開發(fā)現(xiàn)狀和我國(guó)開發(fā)潛力研究
    久久香蕉精品热| 久久人妻av系列| 色av中文字幕| 欧美乱码精品一区二区三区| 亚洲中文字幕日韩| 欧美在线黄色| 露出奶头的视频| 亚洲国产欧美一区二区综合| 法律面前人人平等表现在哪些方面| 日日摸夜夜添夜夜添小说| 婷婷精品国产亚洲av| 国产精品久久久久久亚洲av鲁大| 长腿黑丝高跟| 1024手机看黄色片| 久久国产精品人妻蜜桃| 亚洲三区欧美一区| 美女高潮喷水抽搐中文字幕| 欧美人与性动交α欧美精品济南到| 亚洲第一电影网av| 亚洲五月色婷婷综合| 国产三级在线视频| 99国产综合亚洲精品| cao死你这个sao货| 国产高清有码在线观看视频 | 又黄又爽又免费观看的视频| 黑人欧美特级aaaaaa片| 久久伊人香网站| 成人国语在线视频| 国产亚洲欧美98| 人人妻人人澡人人看| 男人的好看免费观看在线视频 | 国产国语露脸激情在线看| 国产色视频综合| 亚洲国产精品成人综合色| 1024视频免费在线观看| 老鸭窝网址在线观看| 真人一进一出gif抽搐免费| 久久久国产成人精品二区| 中出人妻视频一区二区| 久久人妻av系列| 看免费av毛片| 午夜久久久久精精品| 大香蕉久久成人网| 亚洲黑人精品在线| 女性生殖器流出的白浆| 一二三四社区在线视频社区8| 国产精品免费一区二区三区在线| 丁香欧美五月| 日韩 欧美 亚洲 中文字幕| 91成人精品电影| 亚洲国产欧美日韩在线播放| 亚洲精品一卡2卡三卡4卡5卡| 亚洲色图av天堂| 又黄又爽又免费观看的视频| 亚洲人成网站在线播放欧美日韩| 国产一区在线观看成人免费| 久久人妻福利社区极品人妻图片| 最新在线观看一区二区三区| 在线十欧美十亚洲十日本专区| 丁香欧美五月| 日韩高清综合在线| 亚洲欧美日韩高清在线视频| 悠悠久久av| 国产精品免费一区二区三区在线| 午夜福利视频1000在线观看| 欧美国产日韩亚洲一区| 在线观看www视频免费| 丰满的人妻完整版| 俺也久久电影网| 高清在线国产一区| 国产一级毛片七仙女欲春2 | 中出人妻视频一区二区| 精品国产乱子伦一区二区三区| 香蕉丝袜av| 日韩成人在线观看一区二区三区| 不卡一级毛片| 中文亚洲av片在线观看爽| 欧美绝顶高潮抽搐喷水| 欧美 亚洲 国产 日韩一| 色婷婷久久久亚洲欧美| 99国产精品一区二区三区| 搡老熟女国产l中国老女人| 欧美绝顶高潮抽搐喷水| 国产亚洲精品一区二区www| 丁香六月欧美| 亚洲av成人一区二区三| 日日爽夜夜爽网站| 日日夜夜操网爽| 久久久久久久久中文| 亚洲 欧美一区二区三区| 韩国av一区二区三区四区| 亚洲中文日韩欧美视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产黄色小视频在线观看| xxxwww97欧美| 国产伦在线观看视频一区| 久久久久亚洲av毛片大全| a级毛片a级免费在线| 最近在线观看免费完整版| 可以在线观看毛片的网站| 在线观看免费视频日本深夜| 18禁观看日本| 日本五十路高清| 99久久国产精品久久久| 在线av久久热| 久久人妻av系列| 精品久久久久久成人av| 精品久久久久久久人妻蜜臀av| 欧美乱色亚洲激情| 无限看片的www在线观看| 色老头精品视频在线观看| 国产伦在线观看视频一区| 免费在线观看视频国产中文字幕亚洲| 曰老女人黄片| 国产一级毛片七仙女欲春2 | a级毛片在线看网站| 禁无遮挡网站| 亚洲人成伊人成综合网2020| 国产v大片淫在线免费观看| 久久人妻av系列| 在线国产一区二区在线| 日本成人三级电影网站| 18禁黄网站禁片免费观看直播| 国产高清有码在线观看视频 | 国产精品一区二区免费欧美| 国产精品久久久久久亚洲av鲁大| 69av精品久久久久久| 亚洲国产中文字幕在线视频| 国产精品1区2区在线观看.| 亚洲第一青青草原| 免费在线观看黄色视频的| 色综合婷婷激情| 色在线成人网| 这个男人来自地球电影免费观看| 久久国产亚洲av麻豆专区| 2021天堂中文幕一二区在线观 | 国产亚洲av高清不卡| 欧美激情 高清一区二区三区| 禁无遮挡网站| 国产1区2区3区精品| 午夜免费成人在线视频| 18美女黄网站色大片免费观看| 欧美大码av| 亚洲精品国产区一区二| 午夜亚洲福利在线播放| 欧美大码av| 午夜福利高清视频| 欧美zozozo另类| 国产人伦9x9x在线观看| 久久久久久九九精品二区国产 | 日本撒尿小便嘘嘘汇集6| 搞女人的毛片| 婷婷六月久久综合丁香| 无人区码免费观看不卡| 久久精品国产综合久久久| 亚洲国产精品sss在线观看| 亚洲av熟女| 欧美中文综合在线视频| 国产成人欧美| 免费观看人在逋| 午夜久久久在线观看| 一级a爱视频在线免费观看| 天堂动漫精品| e午夜精品久久久久久久| 欧美日本视频| 国产蜜桃级精品一区二区三区| 国产精华一区二区三区| 国产成人精品久久二区二区免费| 久久精品国产亚洲av香蕉五月| www.www免费av| 免费在线观看日本一区| 国产成人系列免费观看| 青草久久国产| 一本综合久久免费| 久久久久久久久久黄片| 国产亚洲欧美98| 亚洲国产精品合色在线| 久久香蕉国产精品| 精品日产1卡2卡| 女性被躁到高潮视频| 最好的美女福利视频网| 日本 av在线| 岛国在线观看网站| 激情在线观看视频在线高清| 91成人精品电影| 成在线人永久免费视频| 可以在线观看的亚洲视频| 午夜久久久在线观看| 手机成人av网站| 亚洲精品一区av在线观看| 久久久久久久久久黄片| 国内毛片毛片毛片毛片毛片| 91成年电影在线观看| 国产熟女午夜一区二区三区| 91大片在线观看| 久久精品91蜜桃| 老汉色av国产亚洲站长工具| 国产真实乱freesex| 欧美一级a爱片免费观看看 | 国产国语露脸激情在线看| e午夜精品久久久久久久| 亚洲天堂国产精品一区在线| 欧美绝顶高潮抽搐喷水| 国产男靠女视频免费网站| 精品乱码久久久久久99久播| 视频区欧美日本亚洲| 制服诱惑二区| 99热这里只有精品一区 | 天天添夜夜摸| 国产av又大| 亚洲欧洲精品一区二区精品久久久| 亚洲国产精品合色在线| 亚洲国产精品999在线| 国内毛片毛片毛片毛片毛片| 非洲黑人性xxxx精品又粗又长| 亚洲va日本ⅴa欧美va伊人久久| 国产一卡二卡三卡精品| 欧美黄色淫秽网站| av免费在线观看网站| 欧美黑人精品巨大| 久久久水蜜桃国产精品网| 亚洲av成人一区二区三| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美黑人欧美精品刺激| 99久久无色码亚洲精品果冻| 欧美色欧美亚洲另类二区| 怎么达到女性高潮| 成人特级黄色片久久久久久久| 亚洲片人在线观看| 天天添夜夜摸| xxxwww97欧美| 精品高清国产在线一区| 亚洲人成网站在线播放欧美日韩| 国产又色又爽无遮挡免费看| 亚洲精品粉嫩美女一区| 精品欧美国产一区二区三| 怎么达到女性高潮| 亚洲成人久久性| 宅男免费午夜| 国产一区二区激情短视频| 啪啪无遮挡十八禁网站| 在线观看免费日韩欧美大片| 啦啦啦观看免费观看视频高清| 老司机靠b影院| 亚洲中文av在线| 非洲黑人性xxxx精品又粗又长| 免费一级毛片在线播放高清视频| 久久人妻av系列| 中文字幕高清在线视频| 一进一出抽搐gif免费好疼| 国产一区二区在线av高清观看| 久99久视频精品免费| 怎么达到女性高潮| 亚洲精品国产区一区二| 亚洲久久久国产精品| 视频在线观看一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 757午夜福利合集在线观看| 国产在线观看jvid| 国产成人精品无人区| 在线av久久热| 欧美性长视频在线观看| 久久狼人影院| 亚洲黑人精品在线| 日本一本二区三区精品| 亚洲中文字幕一区二区三区有码在线看 | 午夜福利成人在线免费观看| 少妇熟女aⅴ在线视频| 一区二区日韩欧美中文字幕| av超薄肉色丝袜交足视频| 欧美大码av| ponron亚洲| 日韩欧美免费精品| 老司机午夜福利在线观看视频| 亚洲一区中文字幕在线| 久久久国产欧美日韩av| 中文字幕另类日韩欧美亚洲嫩草| 日本一区二区免费在线视频| 国产伦在线观看视频一区| 精品久久久久久,| 日日干狠狠操夜夜爽| 久久草成人影院| 亚洲自偷自拍图片 自拍| 亚洲精品一区av在线观看| 18禁美女被吸乳视频| 午夜精品在线福利| 亚洲欧美日韩高清在线视频| 黄网站色视频无遮挡免费观看| 国产av不卡久久| 最近最新免费中文字幕在线| 欧美黄色淫秽网站| 国产激情偷乱视频一区二区| 此物有八面人人有两片| av有码第一页| 成年免费大片在线观看| 天堂影院成人在线观看| 色老头精品视频在线观看| 日韩 欧美 亚洲 中文字幕| 国产黄色小视频在线观看| 日韩一卡2卡3卡4卡2021年| 啦啦啦韩国在线观看视频| 精品国产超薄肉色丝袜足j| 日本 av在线| 天天添夜夜摸| 欧美黄色淫秽网站| 美女国产高潮福利片在线看| 午夜福利在线在线| 天天添夜夜摸| tocl精华| 高清在线国产一区| 黄网站色视频无遮挡免费观看| 一二三四社区在线视频社区8| 美女国产高潮福利片在线看| 大型黄色视频在线免费观看| 日韩视频一区二区在线观看| av福利片在线| 日日爽夜夜爽网站| 久久久久久人人人人人| 亚洲片人在线观看| 日本 欧美在线| 亚洲男人的天堂狠狠| 久久国产精品影院| 黄色女人牲交| 国产黄色小视频在线观看| 国产成人欧美在线观看| 少妇熟女aⅴ在线视频| 国产又爽黄色视频| 精品久久久久久久末码| 熟女少妇亚洲综合色aaa.| 一边摸一边做爽爽视频免费| 草草在线视频免费看| 久久亚洲真实| 亚洲第一电影网av| 亚洲中文字幕一区二区三区有码在线看 | 国产真实乱freesex| 丰满人妻熟妇乱又伦精品不卡| 久久国产精品影院| 色综合站精品国产| 国产精品亚洲av一区麻豆| 侵犯人妻中文字幕一二三四区| 国产一区二区激情短视频| 成人亚洲精品av一区二区| 老司机靠b影院| 亚洲av第一区精品v没综合| 精品第一国产精品| 亚洲自拍偷在线| 日韩高清综合在线| 亚洲国产欧美网| 久久亚洲精品不卡| 免费在线观看亚洲国产| 麻豆成人av在线观看| 国产精品精品国产色婷婷| 亚洲第一电影网av| 天天躁夜夜躁狠狠躁躁| 国内毛片毛片毛片毛片毛片| 熟女少妇亚洲综合色aaa.| 国产一级毛片七仙女欲春2 | 日本一本二区三区精品| ponron亚洲| 免费观看人在逋| 色尼玛亚洲综合影院| 国产成人欧美| 日本免费a在线| 中文字幕久久专区| 日本五十路高清| 国产高清有码在线观看视频 | xxxwww97欧美| 999久久久精品免费观看国产| 狠狠狠狠99中文字幕| 国产精品野战在线观看| 在线观看一区二区三区| 两性夫妻黄色片| 中文字幕精品免费在线观看视频| 国产高清激情床上av| 中文字幕人妻熟女乱码| 99国产精品一区二区三区| 男女床上黄色一级片免费看| 久久草成人影院| 午夜视频精品福利| 精品久久久久久久末码| 人人妻人人看人人澡| av欧美777| 国产熟女午夜一区二区三区| 好看av亚洲va欧美ⅴa在| 两个人免费观看高清视频| 黄色片一级片一级黄色片| 国产亚洲av高清不卡| 久久中文字幕人妻熟女| 99热这里只有精品一区 | 狠狠狠狠99中文字幕| 日韩欧美在线二视频| 久久精品人妻少妇| 色在线成人网| 日韩大码丰满熟妇| 久久精品国产99精品国产亚洲性色| 侵犯人妻中文字幕一二三四区| 在线永久观看黄色视频| 亚洲精品久久国产高清桃花| 丁香六月欧美| 香蕉av资源在线| a在线观看视频网站| 成人18禁高潮啪啪吃奶动态图| 一区二区三区激情视频| 国产一卡二卡三卡精品| e午夜精品久久久久久久| 国产亚洲精品综合一区在线观看 | 久久香蕉国产精品| 1024视频免费在线观看| 一级毛片精品| 少妇粗大呻吟视频| 日韩欧美一区二区三区在线观看| 久久精品成人免费网站| 精品久久久久久成人av| 无限看片的www在线观看| 亚洲精品在线观看二区| 露出奶头的视频| 99久久国产精品久久久| 久久久精品国产亚洲av高清涩受| 成人18禁在线播放| 成人国产综合亚洲| 久久国产乱子伦精品免费另类| 免费高清在线观看日韩| 国产麻豆成人av免费视频| 淫妇啪啪啪对白视频| 午夜免费成人在线视频| 久久草成人影院| 自线自在国产av| 国产精品久久久久久人妻精品电影| 亚洲av五月六月丁香网| a级毛片a级免费在线| 亚洲avbb在线观看| 91麻豆精品激情在线观看国产| 久久精品国产99精品国产亚洲性色| 一级a爱片免费观看的视频| 国产精品美女特级片免费视频播放器 | 亚洲久久久国产精品| 亚洲avbb在线观看| 国产在线观看jvid| 欧洲精品卡2卡3卡4卡5卡区| 2021天堂中文幕一二区在线观 | 亚洲精品中文字幕在线视频| 国产精品 国内视频| 国产视频一区二区在线看| 女性生殖器流出的白浆| 美女扒开内裤让男人捅视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产精品久久男人天堂| 精品人妻1区二区| 成年女人毛片免费观看观看9| 国产一区二区三区视频了| 欧美日韩乱码在线| 曰老女人黄片| 听说在线观看完整版免费高清| 国产成人av激情在线播放| 国产午夜精品久久久久久| 丰满人妻熟妇乱又伦精品不卡| 99精品欧美一区二区三区四区| 久久久久久亚洲精品国产蜜桃av| 国产亚洲精品一区二区www| 99热这里只有精品一区 | 精品第一国产精品| 色在线成人网| 亚洲五月天丁香| 国产精品一区二区免费欧美| av欧美777| 国产国语露脸激情在线看| 亚洲精品国产区一区二| 久久人妻福利社区极品人妻图片| 人妻丰满熟妇av一区二区三区| 欧美+亚洲+日韩+国产| 视频在线观看一区二区三区| 午夜福利高清视频| 少妇裸体淫交视频免费看高清 | 老司机深夜福利视频在线观看| 久久亚洲真实| 白带黄色成豆腐渣| 99在线视频只有这里精品首页| 欧美三级亚洲精品| 非洲黑人性xxxx精品又粗又长| 国产一区二区三区视频了| 国内精品久久久久久久电影| 色哟哟哟哟哟哟| 久久这里只有精品19| 18禁黄网站禁片免费观看直播| 可以免费在线观看a视频的电影网站| 在线永久观看黄色视频| 成人18禁在线播放| 琪琪午夜伦伦电影理论片6080| 国产99久久九九免费精品| 精品国产乱码久久久久久男人| 欧美日本视频| 丝袜人妻中文字幕| 国产精品99久久99久久久不卡| 欧美人与性动交α欧美精品济南到| 亚洲精品一卡2卡三卡4卡5卡| 久久久久九九精品影院| 大型黄色视频在线免费观看| 18禁黄网站禁片午夜丰满| 中文字幕人成人乱码亚洲影| 亚洲人成伊人成综合网2020| 亚洲欧美日韩无卡精品| 国内久久婷婷六月综合欲色啪| 亚洲国产精品sss在线观看| www国产在线视频色| 精品一区二区三区视频在线观看免费| 免费高清视频大片| 97超级碰碰碰精品色视频在线观看| 亚洲精品在线美女| 久久人妻福利社区极品人妻图片| 色综合婷婷激情| 久热这里只有精品99| x7x7x7水蜜桃| 欧美日本视频| 亚洲国产欧洲综合997久久, | 男人舔奶头视频| 黄频高清免费视频| 久久精品亚洲精品国产色婷小说| 亚洲av美国av| 欧美一级a爱片免费观看看 | 亚洲欧洲精品一区二区精品久久久| 亚洲成人免费电影在线观看| 19禁男女啪啪无遮挡网站| 午夜成年电影在线免费观看| 麻豆久久精品国产亚洲av| 在线免费观看的www视频| bbb黄色大片| 成人国产综合亚洲| ponron亚洲| 欧美性猛交╳xxx乱大交人| 可以在线观看毛片的网站| 久久久久亚洲av毛片大全| 香蕉丝袜av| 亚洲第一欧美日韩一区二区三区| 久热爱精品视频在线9| 国内揄拍国产精品人妻在线 | 国产精品九九99| 国产人伦9x9x在线观看| 伊人久久大香线蕉亚洲五| 精品久久久久久久末码| 丁香六月欧美| 久久香蕉精品热| 欧美成人一区二区免费高清观看 | 日韩大尺度精品在线看网址| 精品国产美女av久久久久小说| 热99re8久久精品国产| 亚洲色图av天堂| 久久久久久大精品| 看黄色毛片网站| 亚洲激情在线av| 性欧美人与动物交配| 一级黄色大片毛片| 女性生殖器流出的白浆| 在线播放国产精品三级| 妹子高潮喷水视频| 丝袜美腿诱惑在线| 午夜成年电影在线免费观看| 国语自产精品视频在线第100页| 在线永久观看黄色视频| av福利片在线| 91字幕亚洲| 一边摸一边抽搐一进一小说| 午夜两性在线视频| 制服丝袜大香蕉在线| 日本熟妇午夜| 19禁男女啪啪无遮挡网站| 免费观看精品视频网站| 巨乳人妻的诱惑在线观看| 一级片免费观看大全| 又黄又爽又免费观看的视频| 男女床上黄色一级片免费看| 夜夜夜夜夜久久久久| 欧美午夜高清在线| 国产高清videossex| 99re在线观看精品视频| 欧洲精品卡2卡3卡4卡5卡区| 久久欧美精品欧美久久欧美| 国产97色在线日韩免费| 国产亚洲精品久久久久久毛片| 亚洲性夜色夜夜综合| 欧美性长视频在线观看| 国产精品自产拍在线观看55亚洲| 久久久久免费精品人妻一区二区 | 每晚都被弄得嗷嗷叫到高潮| 午夜亚洲福利在线播放| 日本熟妇午夜| 国产精品精品国产色婷婷| 国产乱人伦免费视频| a在线观看视频网站| 97人妻精品一区二区三区麻豆 | 国产精品国产高清国产av| 日韩欧美免费精品| 看黄色毛片网站| 国产亚洲av嫩草精品影院| 亚洲一区中文字幕在线| 国产精华一区二区三区| 久久狼人影院| 天堂影院成人在线观看| а√天堂www在线а√下载| 久久久久久人人人人人| 99热6这里只有精品| 深夜精品福利| 国产黄片美女视频| aaaaa片日本免费| 成人三级黄色视频| 成人一区二区视频在线观看| 男女视频在线观看网站免费 | 亚洲 欧美一区二区三区| 久久久水蜜桃国产精品网| 老司机福利观看| 伊人久久大香线蕉亚洲五| 首页视频小说图片口味搜索| 又紧又爽又黄一区二区| 欧美人与性动交α欧美精品济南到| 亚洲国产高清在线一区二区三 | 久久中文字幕一级|