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

    基于氣流床氣化爐熔渣結(jié)構(gòu)的黏度預(yù)測

    2017-10-16 08:13:24王倩李強張建勝
    化工學(xué)報 2017年10期
    關(guān)鍵詞:熔渣煤灰熔體

    王倩,李強,張建勝

    (清華大學(xué)熱能工程系,熱科學(xué)與動力工程教育部重點實驗室,北京 100084)

    基于氣流床氣化爐熔渣結(jié)構(gòu)的黏度預(yù)測

    王倩,李強,張建勝

    (清華大學(xué)熱能工程系,熱科學(xué)與動力工程教育部重點實驗室,北京 100084)

    在煤氣化工藝中,煤灰高溫黏度特性是影響熔渣流動特性重要參數(shù)之一。以NBO/T、P/M為結(jié)構(gòu)參數(shù),建立熔渣高溫黏度預(yù)測模型,對結(jié)構(gòu)參數(shù)、熔渣黏度的關(guān)系進(jìn)行了研究,并通過文獻(xiàn)的實驗數(shù)據(jù)驗證了模型的準(zhǔn)確性。研究發(fā)現(xiàn)隨著NBO/T的增大,熔體聚合程度降低,熔體流動時內(nèi)摩擦力減小,黏度降低;隨著P/M增大,熔體結(jié)構(gòu)強度減弱,熔體流動所需要克服的阻力增加,熔渣黏度逐漸升高。引入結(jié)構(gòu)參數(shù)后的熔渣黏度預(yù)測模型對一定范圍內(nèi)的熔渣黏度具有良好的預(yù)測效果,在NBO/T < 0.5時,黏度預(yù)測值和實驗值誤差在10% 以內(nèi),采用兩種結(jié)構(gòu)參數(shù)均可;在NBO/T > 0.5時,黏度預(yù)測值和實驗值誤差在25% 以內(nèi),采用結(jié)構(gòu)參數(shù)P/M更合適。

    熔渣結(jié)構(gòu);黏度;氣化;模型;NBO/T;P/M

    Abstract:The viscosity of molten slags is one of the most significant parameters to predict the flow properties of ash in coal gasification process.To estimate the viscosity of molten slags,a model based on the structure of slags was established.NBO/T (number of non-bridging oxygen/tetragonally-bonded oxygen) and P/M(polymer/modifier),as the structural parameters of slags,provide description of the polymerization and bond strengths of molten slags,respectively.The model presented the dependency on both temperature and structural parameters with the viscosity in the SiO2-Al2O3-CaO-MgO-Fe2O3-K2O-Na2O-TiO2slag systems.The accuracy of the model was verified by the experimental data of literatures.The results indicate that the viscosity decreases with the increase of NBO/T due to the reduction of polymerization in the molten slags,whereas the viscosity rises with the increase of P/M because of the strength weaken of molten slags.A good agreement is obtained from the calculated and measured values of the viscosity within a certain range.The deviation between the prediction and experiment results is within 10% when NBO/T < 0.5 and 25% when NBO/T > 0.5.

    Key words:slag structure; viscosity; gasification; model; NBO/T; P/M

    引 言

    煤氣化技術(shù)是煤炭潔凈高效利用的龍頭技術(shù)之一,目前工業(yè)化大型化的工藝以氣流床氣化爐工藝為主。基于氣流床氣化爐工藝高溫、高壓、高熱載荷的特點,操作溫度一般高于所燒煤種的灰熔點,采用液態(tài)排渣工藝。氣化爐運行過程中,由于排渣不暢導(dǎo)致排渣口堵塞甚至迫使氣化爐停爐的現(xiàn)象時常發(fā)生,常采用升溫、減負(fù)荷的方式保證氣化爐的連續(xù)排渣,但是長期高溫運行會造成氧耗、煤耗的增加,加快耐火磚的侵蝕[1]。

    熔渣的流動特性是影響氣化爐穩(wěn)定運行的重要因素,直接關(guān)系到灰渣能否順利排出。熔渣流動特性最重要的參數(shù)之一是煤灰高溫黏度特性,由于氣流床氣化爐在高溫高壓下運行,受到實驗手段的限制,難以直接觀測高溫灰渣的流動現(xiàn)象,建立黏度預(yù)測模型是研究煤灰黏溫特性的有效手段之一。

    煤灰在高溫下是一種硅酸鹽熔體,其黏度特性與熔體結(jié)構(gòu)和溫度密切相關(guān)。黏溫曲線通常用來表征煤灰高溫黏度特性,煤灰熔渣的黏度預(yù)測公式一般是基于組成結(jié)構(gòu)和黏溫數(shù)據(jù),通過數(shù)據(jù)回歸的方式得到。目前擬合得到的經(jīng)驗公式大多數(shù)是以Arrhenius[2-4]、 Weymann[5-7]、 Vogel-Fulcher-Tamman[8-9]模型為基礎(chǔ),建立黏度和溫度的指數(shù)關(guān)系,并采用煤灰化學(xué)組成表示模型參數(shù)。Riboud等[10]將模型參數(shù)表示為網(wǎng)絡(luò)形成組分的線性加和,由于與實際的非線性變化不符而很少采用。Urbain[5]則將模型參數(shù)表示為有關(guān)組分的多項式,擬合過程涉及的參數(shù)較多,但對牛頓流體特性煤灰渣的預(yù)測效果較好[11]。

    煤灰組分對硅酸鹽熔體黏度的影響比較復(fù)雜,基于煤灰組分的黏度預(yù)測模型多采用氧化物的百分含量作為擬合參數(shù)[12-14],來量化組分對黏度的影響。然而,煤灰組分是通過改變?nèi)垠w結(jié)構(gòu)而影響?zhàn)ざ鹊?,為了直接探究熔體結(jié)構(gòu)對黏度的影響,一些研究者提出采用熔渣聚合度(NBO/T)、光學(xué)堿度(Λ)、熔渣強度(P/M)等結(jié)構(gòu)參數(shù)來表征熔體結(jié)構(gòu)。研究者們[5,15]建立了擬合參數(shù)和NBO/T之間的關(guān)系,研究了熔體結(jié)構(gòu)對黏度的影響。Mills等[2]為了考慮陽離子的作用,采用Λ構(gòu)建熔渣黏度預(yù)測的模型參數(shù),發(fā)現(xiàn)對于大多數(shù)組分的灰渣預(yù)測能力強于 Riboud等和 Urbain的模型,但是對煤灰渣的預(yù)測效果較差。而Zhang等[16-18]提出直接采用橋氧、非橋氧和自由氧的數(shù)目來表示模型參數(shù),并較好地預(yù)測了四元和五元模擬渣系,但有待對實際灰渣進(jìn)行預(yù)測。此外,熊友輝等[19]提出了P/M來描述熔體結(jié)構(gòu),用于分析熔體結(jié)構(gòu)對煤灰物理化學(xué)性質(zhì)的影響。

    綜上所述,由于熔渣結(jié)構(gòu)的復(fù)雜性,缺乏對熔體結(jié)構(gòu)仍缺乏全面的認(rèn)識。本研究通過引入熔體結(jié)構(gòu)參數(shù),構(gòu)建黏度預(yù)測模型,探究熔體結(jié)構(gòu)參數(shù)對黏度的影響,更深入地認(rèn)識煤灰黏度的影響機理,為宏觀調(diào)控配煤提供指導(dǎo)性建議。

    1 實驗部分

    1.1 實驗煤樣組分及黏度測量

    實驗選取了11個典型高硅鋁、高硅鋁比煤樣作為研究對象,根據(jù) GB/T 1574的要求制成(815±10)℃的標(biāo)準(zhǔn)灰樣,并采用 X射線熒光光譜分析獲取灰樣中各氧化物的組成成分,分析數(shù)據(jù)見表1。

    表1 煤灰組分Table 1 Components of coal ash

    依據(jù)DL/T 660的標(biāo)準(zhǔn),實驗采用ND-3數(shù)字式高溫旋轉(zhuǎn)黏度計測量灰樣的黏度特性。實驗過程中將煤灰加熱至預(yù)定溫度,使煤灰完全熔融,采用逐個溫度點測量,降溫間隔為10℃,在對應(yīng)溫度點恒溫30 min保證熔渣系統(tǒng)達(dá)到平衡。每個溫度點下測定3次偏轉(zhuǎn)角,取算數(shù)平均值計算對應(yīng)溫度下的黏度值。測試氣氛為弱還原性氣氛(CO/CO2=3/2)。

    1.2 實驗煤樣組分及黏度測量

    1.2.1 熔體結(jié)構(gòu)參數(shù) 熔渣的主要成分是硅酸鹽熔體,其基本結(jié)構(gòu)為 [SiO4]4?四面體構(gòu)成的三維隨機網(wǎng)絡(luò)結(jié)構(gòu)[20-21]。在熔體中氧原子主要以橋氧、非橋氧和自由氧3種結(jié)構(gòu)狀態(tài)存在,3類氧原子的含量比及分布代表熔體聚合程度[22-23]。

    熔渣的聚合度大小直接影響其黏度。不同氧化物在網(wǎng)絡(luò)結(jié)構(gòu)形成過程中作用不同。SiO2是網(wǎng)絡(luò)形成氧化物,是熔渣網(wǎng)絡(luò)結(jié)構(gòu)的主要形成組分,一般與氧原子形成橋氧鍵,SiO2含量越大,熔體結(jié)構(gòu)聚合程度越高[24]。Na2O、K2O、CaO、MgO 等低價氧化物是網(wǎng)絡(luò)改性氧化物,高溫下通過破壞 Si—O共價鍵,形成非橋氧鍵,導(dǎo)致網(wǎng)絡(luò)結(jié)構(gòu)解體,聚合程度降低,黏度減小[25-26]。而Al2O3、Fe2O3為網(wǎng)絡(luò)兩性氧化物,在熔體結(jié)構(gòu)中既可充當(dāng)網(wǎng)絡(luò)形成氧化物,也可成為網(wǎng)絡(luò)改性氧化物,對網(wǎng)絡(luò)結(jié)構(gòu)的影響與自身氧化物和堿(土)金屬氧化物的含量有關(guān)[27]。因此 Mills等[2]采用 NBO/T(熔體中非橋氧數(shù)量與熔體形成四面體聚合物的數(shù)量的比值,即每個四面體中的非橋氧數(shù)量)來衡量熔體聚合程度,NBO/T越大,表明熔渣聚合度越低。

    此外,硅酸鹽熔體的黏度與熔體中的化學(xué)鍵的強度有關(guān),一般堿金屬、堿土金屬按照K>Na>Li,Ba>Sr>Ca>Mg 的順序降低熔體的黏度[28]。因此采用結(jié)構(gòu)參數(shù)P/M(網(wǎng)絡(luò)形成氧化物結(jié)構(gòu)強度與網(wǎng)絡(luò)改性氧化物結(jié)構(gòu)強度之比)來表示熔體化學(xué)鍵強度,其中鍵強用不同金屬陽離子的離子勢表示。從物理意義上講,P/M反映了煤灰熔體的庫侖締合力,這種締合力可能直接反映了熔體的強度。

    基于煤灰特性,在氣流床氣化爐中的高溫熔渣主要以SiO2-Al2O3-CaO-MgO-Fe2O3-K2O-Na2O-TiO2的八元渣系表示。采用NBO/T和P/M兩種結(jié)構(gòu)參數(shù)表征熔體結(jié)構(gòu),進(jìn)行模型構(gòu)建。研究中SiO2、TiO2和 Al2O3為網(wǎng)絡(luò)形成氧化物,Na2O、K2O、CaO、MgO、Fe2O3為網(wǎng)絡(luò)改性氧化物。NBO/T和P/M的計算公式如下

    式中,Xi為氧化物摩爾分?jǐn)?shù);L(i)為陽離子i的離子勢;下角標(biāo)poly指網(wǎng)絡(luò)形成氧化物,modi指網(wǎng)絡(luò)改性氧化物。具體見表2。

    表2 氧化物的離子勢Table 2 Ionic potential of oxides

    1.2.2 黏度預(yù)測模型 煤灰黏度預(yù)測主要是基于溫度和煤灰成分,對于確定組分的熔渣,在給定溫度下對應(yīng)一個確定的黏度值。高硅鋁及高硅鋁比煤灰的黏溫特性表現(xiàn)為牛頓流體特性,采用Urbain模型的預(yù)測結(jié)果強于其他模型[27]。因此,研究主要是在Urbain模型的基礎(chǔ)上,引入熔體結(jié)構(gòu)參數(shù),建立新的黏度預(yù)測模型。

    基于煤灰典型的8種氧化物組分,研究中采用Weymann方程描述熔渣黏度與溫度的關(guān)系

    式中,A、B為模型參數(shù),與煤灰組分密切相關(guān)。對于給定組分的熔渣,根據(jù)其黏溫數(shù)據(jù)利用式(3)擬合得到lnA和B,通過式(1)、式(2)計算結(jié)構(gòu)參數(shù)NBO/T和 P/M,根據(jù)式(4)擬合計算得到m、n、k。模型參數(shù)B的各項系數(shù)可通過式(5)~式(7)擬合得到。其中,S指代結(jié)構(gòu)參數(shù)NBO/T或P/M。

    2 實驗結(jié)果與討論

    2.1 黏度預(yù)測模型結(jié)果

    基于熔渣結(jié)構(gòu)的預(yù)測模型參數(shù)如表3和表4所示。對于給定組分的煤灰,根據(jù)煤灰組分含量可得模型參數(shù)B和結(jié)構(gòu)參數(shù)NBO/T或P/M,再依據(jù)表3中的公式計算lnA,將參數(shù)lnA和B代入Weymann方程[式(3)]即可獲得不同溫度下的黏度預(yù)測值,圖1所示為兩個模型對11個灰樣黏度的預(yù)測值和實驗值對比。由圖1可知,引入結(jié)構(gòu)參數(shù)后,預(yù)測值和實驗值的吻合度較高,模型的預(yù)測效果很好,低黏度下的預(yù)測效果要優(yōu)于高黏度下的預(yù)測結(jié)果。

    表3 黏度預(yù)測模型Table 3 Models of estimating viscosity

    表4 模型參數(shù)B的各項系數(shù)Table 4 Coefficients of model parameterB

    圖1 熔渣黏度模型預(yù)測值和實驗值對比Fig.1 Comparison between model calculation and experimental data of slag viscosity

    2.2 結(jié)構(gòu)參數(shù)與黏度的關(guān)系

    2.2.1 NBO/T與黏度的關(guān)系 圖2所示為不同溫度下黏度和 NBO/T的關(guān)系,從圖中可知,隨著溫度的降低,熔渣逐漸從液態(tài)渣向固態(tài)渣轉(zhuǎn)變,黏度逐漸升高。隨著NBO/T的增加,熔渣黏度降低。這是由于網(wǎng)絡(luò)改性氧化物如堿土金屬增加,破壞Si—O共價鍵,形成非橋氧鍵,導(dǎo)致網(wǎng)絡(luò)結(jié)構(gòu)變松散或者解聚,熔體聚合程度降低,使得熔體流動時內(nèi)摩擦力減小,從而使得黏度降低。

    圖2 不同溫度下熔渣黏度和NBO/T的關(guān)系Fig.2 Relation between slag viscosity and NBO/T at different temperatures

    從圖2可發(fā)現(xiàn)存在NBO/T<0的情況,這是由于網(wǎng)絡(luò)中性氧化物 Al2O3在熔體中作用發(fā)生了變化。圖3所示為不同Al2O3含量下的NBO/T,當(dāng)NBO/T>0 時,[AlO6]9?轉(zhuǎn)化為[AlO4]四面體替代[SiO4]四面體,充當(dāng)網(wǎng)絡(luò)形成氧化物的作用,含量越高,NBO/T越小,黏度越大;當(dāng)NBO/T<0時,部分Al的配位數(shù)從6變?yōu)?,起到網(wǎng)絡(luò)改性氧化物的作用,其黏度與充當(dāng)網(wǎng)絡(luò)改性作用的 Al2O3含量有關(guān)。

    圖3 Al2O3含量和NBO/T的關(guān)系Fig.3 Relation between NBO/T and mass fraction of Al2O3

    2.2.2 P/M與黏度的關(guān)系 圖4所示為不同溫度下黏度和P/M的關(guān)系。由圖4可知,隨著P/M增加,熔體結(jié)構(gòu)強度增大,熔體流動所需要克服的阻力增加,導(dǎo)致熔渣黏度逐漸升高。從圖4中還可發(fā)現(xiàn)當(dāng)溫度為1600℃時,隨著P/M增加到一定程度,熔渣黏度不再增加,這是由于此時熔渣完全熔融,熔渣結(jié)構(gòu)對黏度的影響減小,而溫度對黏度的影響占主導(dǎo)地位。

    圖4 不同溫度下熔渣黏度和P/M的關(guān)系Fig.4 Relation between slag viscosity and P/M at different temperatures

    2.3 NBO/T和P/M的關(guān)系

    通過上述分析可知,NBO/T、P/M這兩個表征熔體結(jié)構(gòu)的結(jié)構(gòu)參數(shù)均考慮了網(wǎng)絡(luò)形成氧化物和網(wǎng)絡(luò)改性氧化物在熔體形成中的作用,但是 P/M 在NBO/T的基礎(chǔ)上進(jìn)一步區(qū)分了同類氧化物中不同陽離子對熔體結(jié)構(gòu)的影響,如網(wǎng)絡(luò)形成氧化物中Si—O和Al—O的作用,網(wǎng)絡(luò)改性氧化物中K、Na、Ca、Mg的作用。圖5建立了NBO/T和P/M兩個結(jié)構(gòu)參數(shù)間的關(guān)系。

    圖5 NBO/T和P/M的關(guān)系Fig.5 Relation between P/M and NBO/T

    從圖中可以發(fā)現(xiàn),當(dāng)NBO/T<0.5時,隨NBO/T的減小,P/M逐漸降低,同類氧化物中不同金屬陽離子對熔體結(jié)構(gòu)的影響較小,此時可以直接用NBO/T代替 P/M 來表征熔體結(jié)構(gòu)。這是由于隨著NBO/T的減小,網(wǎng)絡(luò)改性體氧化物含量較少,不同堿(土)金屬氧化物陽離子對熔體強度的影響可以忽略。當(dāng)NBO/T>0.5時,隨NBO/T的增大,P/M基本不變,這表明金屬陽離子對熔體結(jié)構(gòu)的作用很大,尤其是堿(土)金屬氧化物的含量增加會導(dǎo)致不同金屬陽離子的作用越來越大,用參數(shù) NBO/T表征熔體結(jié)構(gòu)效果不如P/M,此時采用參數(shù)P/M更為合適。

    2.4 預(yù)測模型應(yīng)用

    為了驗證模型的準(zhǔn)確性,從國內(nèi)外文獻(xiàn)搜集了4組相關(guān)黏度實驗數(shù)據(jù),利用該預(yù)測模型進(jìn)行計算并與實驗值進(jìn)行了對比。以NBO/T和P/M為結(jié)構(gòu)參數(shù)的預(yù)測模型分別如圖6所示、圖7所示,4組數(shù)據(jù)計算得到的NBO/T和P/M的值見表5。

    由圖可知,4組黏度數(shù)據(jù)采用兩種預(yù)測模型預(yù)測的結(jié)果均較為理想。兩個模型對Mills等[2]的兩個渣樣數(shù)據(jù)預(yù)測結(jié)果良好,其誤差均在 8%以內(nèi),這是由于該熔渣的結(jié)構(gòu)參數(shù)NBO/T均小于0.5,兩個參數(shù)均能很好地表征熔體結(jié)構(gòu);但是兩個模型對Duchesne等[29]和Folkedahl等[30]兩個熔渣黏度的預(yù)測結(jié)果誤差分別為 23.6%和 16.7%,這是由于NBO/T>0.5,金屬陽離子對熔體結(jié)構(gòu)影響越來越重要,這與上述對比分析兩參數(shù)的結(jié)果一致。

    其中,對 Mills等的兩個渣樣數(shù)據(jù)預(yù)測結(jié)果優(yōu)于 Duchesne等和 Folkedahl等的結(jié)果,由于Duchesne等和Folkedahl等的熔渣為中硅鋁煤灰,在高溫段表現(xiàn)為非牛頓流體,采用以 Urbain模型為基礎(chǔ)的預(yù)測模型誤差較大。

    圖6 以NBO/T為結(jié)構(gòu)參數(shù)的預(yù)測值與實際值Fig.6 Comparison of experimental data and model calculation for NBO/T as structural parameter

    圖7 以P/M為結(jié)構(gòu)參數(shù)的預(yù)測值與實際值Fig.7 Comparison of experimental data and model calculation for P/M as structural parameter

    表5 煤灰結(jié)構(gòu)參數(shù)Fig.5 Structural parameters of coal slags

    3 結(jié) 論

    以NBO/T、P/M為結(jié)構(gòu)參數(shù),建立熔渣高溫黏度預(yù)測模型,對預(yù)測模型進(jìn)行對比分析及應(yīng)用,主要結(jié)論如下。

    隨著 NBO/T的增大,熔體聚合程度降低,熔體流動時內(nèi)摩擦力減小,黏度降低;隨著P/M增大,熔體結(jié)構(gòu)強度增大,熔體流動所需要克服的阻力增加,熔渣黏度逐漸升高。引入結(jié)構(gòu)參數(shù)后的熔渣黏度預(yù)測模型對一定范圍內(nèi)的熔渣黏度具有良好的預(yù)測效果,在NBO/T<0.5時,黏度預(yù)測值和實驗值誤差在 10%以內(nèi),采用兩種結(jié)構(gòu)參數(shù)均可;在NBO/T>0.5時,黏度預(yù)測值和實驗值誤差在25%以內(nèi),此時采用結(jié)構(gòu)參數(shù)P/M更為合適。

    符 號 說 明

    A,B——預(yù)測模型參數(shù)

    L(i) ——氧化物i中陽離子的離子勢

    NBO/T,P/M ——表征熔渣結(jié)構(gòu)的參數(shù)

    Xi——氧化物i的摩爾分?jǐn)?shù),%

    Xmodi——網(wǎng)絡(luò)改性氧化物i的摩爾分?jǐn)?shù),%

    Xpoly——網(wǎng)絡(luò)形成氧化物i的摩爾分?jǐn)?shù),%

    ηexperiment——實驗得到的黏度,Pa·s

    ηmodel——模型計算得到的黏度,Pa·s

    [1]畢大鵬,趙勇,管清亮,等.水冷壁氣化爐內(nèi)熔渣流動特性模型[J].化工學(xué)報,2015,66(3): 888-895.BI D P,ZHAO Y,GUAN Q L,et al.Modeling slag behavior in membrane wall gasifier[J].CIESC Journal,2015,66(3): 888-895.

    [2]MILLS K C,SRIDHAR S.Viscosities of ironmaking and steelmaking slags[J].Ironmaking & Steelmaking,1999,26(4):262-268.

    [3]OGIBAYASHI S,YAMAGUCHI K,MUKAI T,et al.Mold powder technology for continuous casting of low-carbon aluminum-killed steel[J].Nippon Steel Technical Report.Overseas,1987,(34): 1-10.

    [4]IIDA T,SAKAI H,KITA Y,et al.An equation for accurate prediction of the viscosities of blast furnace type slags from chemical composition[J].ISIJ International,2000,40(Suppl.): S110-S114.

    [5]URBAIN G.Viscosity estimation of slags[J].Steel Research International,1987,58(3): 111-116.

    [6]SENIOR C L,SRINIVASACHAR S.Viscosity of ash particles in combustion systems for prediction of particle sticking[J].Energy &Fuels,1995,9(2): 277-283.

    [7]KONADRATIEV A,JAK E.Review of experimental data and modeling of the viscosities of fully liquid slags in the Al2O3-CaO-‘FeO'-SiO2system[J].Metallurgical and Materials Transactions B,2001,32(6): 1015-1025.

    [8]SONG W J,TANG L H,ZHU X D,et al.Prediction of Chinese coal ash fusion temperatures in Ar and H2atmospheres[J].Energy & Fuels,2009,23(4): 1990-1997.

    [9]BROWNING G J,BRYANT G W,HURST H J,et al.An empirical method for the prediction of coal ash slag viscosity[J].Energy &Fuels,2003,17(3): 731-737.

    [10]RIBOUD P V,ROUX Y,LUCAS L D,et al.Improvement of continuous casting powders[J]. Fachberichte Hüttenpraxis Metallweiterverarbeitung,1981,19(8): 859-869.

    [11]VARGAS S,FRANDSEN F J,DAM-JOHANSEN K.Rheological properties of high-temperature melts of coal ashes and other silicates[J].Progress in Energy and Combustion Science,2001,27(3):237-429.

    [12]NEUVILLE D R.Viscosity,structure and mixing in (Ca,Na) silicate melts[J].Chemical Geology,2006,229(1): 28-41.

    [13]NAKAMOTO M,JOONHO L E E,TANAKA T.A model for estimation of viscosity of molten silicate slag[J].ISIJ International,2005,45(5): 651-656.

    [14]GUPTA V K,SINHA S P,RAJ B.A model for describing viscosity-composition relationship for continuous casting mould fluxes[J].Steel India(India),1998,21(1): 22-29.

    [15]唐續(xù)龍,郭敏,王習(xí)東,等.基于修正的 (NBO/T) 比值的黏度模型[J].北京科技大學(xué)學(xué)報,2010,32(12): 1542-1546.TANG X L,GUO M,WANG X D,et al.Estimation model of viscosity based on modified (NBO/T) ratio[J].Journal of University of Science and Technology Beijing,2010,32(12): 1542-1546.

    [16]ZHANG G H,CHOU K C,XUE Q G,et al.Modeling viscosities of CaO-MgO-FeO-MnO-SiO2molten slags[J].Metallurgical and Materials Transactions B,2012,43(1): 64-72.

    [17]ZHANG G H,CHOU K C,MILLS K.Modelling viscosities of CaO-MgO-Al2O3-SiO2molten slags[J].ISIJ International,2012,52(3): 355-362.

    [18]郭靖,程樹森,趙宏博.基于結(jié)構(gòu)理論的 SiO2-CaO-MgO-Al2O3熔渣黏度的預(yù)報模型[J].鋼鐵研究學(xué)報,2013,25(8):6-11.GUO J,CHENG S S,ZHAO H B.Modeling research for estimating viscosity of SiO2-CaO-MgO-Al2O3system molten slag based on slag structure theory[J].Journal of Iron and Steel Research,2013,25(8):6-11.

    [19]熊友輝,孫學(xué)信.動力用煤灰熔體結(jié)構(gòu)特性及其計算研究[J].煤炭轉(zhuǎn)化,1996,19(4):85-91.XIONG Y H,SUN X X.Research on structure characteristic and computational method of coal ash melts[J].Coal Conversion,1996,19(4): 85-91.

    [20]ZACHARIASEN W H.The atomic arrangement in glass[J].Journal of the American Chemical Society,1932,54: 3841-3851.

    [21]MILLS K.The estimation of slag properties[EB/OL].Southern African Pyrometallurgy 2011 International Conference,2011,http://www.pyrometallurgy.co.za/KenMills/KenMills.pdf.

    [22]MYSEN B O.Structure and Properties of Silicate Melts[M].Amsterdam: Elsevier Science Limited,1988.

    [23]MYSEN B O.Relationships between silicate melt structure and petrologic processes[J].Earth-Science Reviews,1990,27(4):281-365.

    [24]MILLS K C.The influence of structure on the physico-chemical properties of slags[J].ISIJ International,1993,33(1): 148-155.

    [25]IIDA T,SAKAI H,KITA Y,et al.Equation for estimating viscosities of industrial mold fluxes[J].High Temperature Materials and Processes,2000,19(3/4): 153-164.

    [26]KUCHARSKI M,STUBINA N M,TOGURI J M.Viscosity measurements of molten Fe-O-SiO2,Fe-O-CaO-SiO2,and Fe-O-MgO-SiO2slags[J].Canadian Metallurgical Quarterly,1989,28(1): 7-11.

    [27]鄭常昊,王倩,張建勝.氣化中煤灰熔點和黏度預(yù)測模型[J].燃料化學(xué)學(xué)報,2016,44(5): 521-527.ZHENG C H,WANG Q,ZHANG J S.Prediction model of ash fusion temperature and viscosity in coal gasification[J].Journal of Fuel Chemistry and Technology,2016,44(5): 521-527.

    [28]KUO Y M,HUANG K L,WANG C T,et al.Effect of Al2O3mole fraction and cooling method on vitrification of an artificial hazardous material(Ⅰ): Variation of crystalline phases and slag tructures[J].Journal of Hazardous Materials,2009,169(1/2/3):626-634.

    [29]DUCHESNE M A,ILYUSHECHKIN A Y,HUGHES R W,et al.Flow behaviour of slags from coal and petroleum coke blends[J].Fuel,2012,97: 321-328.

    [30]FOLKEDAHL B C,SCHOBERT H H.Effects of atmosphere on viscosity of selected bituminous and low-rank coal ash slags[J].Energy & Fuels,2005,19(1): 208-215.

    Model of estimating viscosity based on slag structure in entrained flow gasifier

    WANG Qian,LI Qiang,ZHANG Jiansheng
    (Key Laboratory for Thermal Science and Power Engineering of Ministry of Education,Department of Thermal Engineering,Tsinghua University,Beijing100084,China)

    TQ 546.4

    A

    0438—1157(2017)10—3877—07

    10.11949/j.issn.0438-1157.20170443

    2017-04-21收到初稿,2017-06-29收到修改稿。

    聯(lián)系人:張建勝。

    王倩(1992—),女,博士研究生。

    國家重點基礎(chǔ)研究發(fā)展計劃項目(2014CB744300)。

    Received date:2017-04-21.

    Corresponding author:Prof.ZHANG Jiansheng,zhang-jsh@tsinghua.edu.cn

    Foundation item:supported by the National Basic Research Program of China(2014CB744300).

    猜你喜歡
    熔渣煤灰熔體
    三塘湖煤田煤灰熔融性溫度預(yù)測模型研究
    煤灰分及CaCO3對鶴崗煤灰熔融性和黏溫特性影響研究
    煤化工(2021年3期)2021-07-14 07:25:30
    MgO-B2O3-SiO2三元體系熔渣表面張力計算
    上海金屬(2016年2期)2016-11-23 05:34:45
    銅在熔態(tài)鐵液與熔渣中平衡分配系數(shù)研究
    上海金屬(2016年3期)2016-11-23 05:19:59
    聚合物熔體脈振傳遞過程的協(xié)同學(xué)研究
    中國塑料(2016年4期)2016-06-27 06:33:48
    注射保壓過程中O2/N2分子在PMMA熔體內(nèi)部的擴(kuò)散行為
    中國塑料(2016年3期)2016-06-15 20:30:01
    TG-FTIR在煤灰研究中的應(yīng)用:真實煤灰和模擬灰受熱氣體釋放規(guī)律比較
    含硅芳炔樹脂及其共混物熔體的流變性能
    煤灰黏溫特性的測試條件
    CaO-A12O3-TiO2熔渣表面張力計算模型
    上海金屬(2014年2期)2014-12-18 06:52:45
    国产精品无大码| 国产成人欧美在线观看 | av女优亚洲男人天堂| 国产精品嫩草影院av在线观看| 天天躁夜夜躁狠狠躁躁| 免费黄频网站在线观看国产| 久久久久久久大尺度免费视频| 中文字幕亚洲精品专区| 在现免费观看毛片| 一边摸一边抽搐一进一出视频| 国产亚洲一区二区精品| av在线老鸭窝| 久久精品人人爽人人爽视色| 久久久久久久久久久久大奶| 久久99精品国语久久久| 极品人妻少妇av视频| 成人18禁高潮啪啪吃奶动态图| 丁香六月欧美| 激情五月婷婷亚洲| 欧美日韩国产mv在线观看视频| 中文乱码字字幕精品一区二区三区| 日本欧美国产在线视频| 久热这里只有精品99| 美女高潮到喷水免费观看| 国产免费现黄频在线看| 久久国产精品大桥未久av| 视频在线观看一区二区三区| 大码成人一级视频| 日本av手机在线免费观看| 欧美人与性动交α欧美软件| 亚洲伊人色综图| 巨乳人妻的诱惑在线观看| 久久久国产精品麻豆| 久久久国产精品麻豆| 满18在线观看网站| 制服诱惑二区| 免费久久久久久久精品成人欧美视频| 国产精品秋霞免费鲁丝片| 国产片内射在线| 日韩欧美精品免费久久| 在线亚洲精品国产二区图片欧美| 久久av网站| 欧美激情极品国产一区二区三区| 日本av免费视频播放| 日日撸夜夜添| 久久国产亚洲av麻豆专区| 日韩 亚洲 欧美在线| 女人高潮潮喷娇喘18禁视频| 欧美激情极品国产一区二区三区| 如日韩欧美国产精品一区二区三区| 日韩中文字幕视频在线看片| 9191精品国产免费久久| 麻豆精品久久久久久蜜桃| 两个人免费观看高清视频| 国产精品一区二区在线观看99| 亚洲av国产av综合av卡| 亚洲av男天堂| 制服人妻中文乱码| 大片免费播放器 马上看| 亚洲欧美日韩另类电影网站| 99久久人妻综合| 亚洲欧美成人精品一区二区| av线在线观看网站| 国产成人一区二区在线| 人体艺术视频欧美日本| 精品人妻熟女毛片av久久网站| 一区二区日韩欧美中文字幕| 欧美黑人精品巨大| 欧美激情高清一区二区三区 | 成人手机av| 亚洲伊人久久精品综合| 亚洲欧洲精品一区二区精品久久久 | 免费av中文字幕在线| 精品久久久精品久久久| 天天躁夜夜躁狠狠躁躁| 一级a爱视频在线免费观看| 一边摸一边做爽爽视频免费| 亚洲国产日韩一区二区| 只有这里有精品99| 国产精品亚洲av一区麻豆 | 美女视频免费永久观看网站| 乱人伦中国视频| 十八禁高潮呻吟视频| 久久亚洲国产成人精品v| 中文天堂在线官网| 男女免费视频国产| 欧美在线一区亚洲| 母亲3免费完整高清在线观看| 1024香蕉在线观看| 免费少妇av软件| 韩国精品一区二区三区| 一级黄片播放器| 亚洲欧美一区二区三区久久| 丰满迷人的少妇在线观看| 亚洲,一卡二卡三卡| 一级毛片电影观看| 亚洲情色 制服丝袜| 日日爽夜夜爽网站| 丁香六月欧美| 又黄又粗又硬又大视频| 777久久人妻少妇嫩草av网站| 最近的中文字幕免费完整| 蜜桃国产av成人99| 自线自在国产av| 久久久久久人妻| 国产亚洲精品第一综合不卡| 一区二区日韩欧美中文字幕| 久久久精品免费免费高清| 1024香蕉在线观看| 三上悠亚av全集在线观看| 国产成人免费无遮挡视频| 观看美女的网站| 中文字幕精品免费在线观看视频| 国产黄色视频一区二区在线观看| 日韩中文字幕欧美一区二区 | 国产男女内射视频| 久久综合国产亚洲精品| 色婷婷久久久亚洲欧美| 久久久久久久国产电影| 老司机靠b影院| 婷婷色综合www| 亚洲精品av麻豆狂野| 丰满少妇做爰视频| www日本在线高清视频| 欧美日韩精品网址| 韩国av在线不卡| 男人舔女人的私密视频| 国产成人系列免费观看| 97在线人人人人妻| 亚洲少妇的诱惑av| 男女无遮挡免费网站观看| 久久亚洲国产成人精品v| 女性生殖器流出的白浆| 母亲3免费完整高清在线观看| 国产精品无大码| 少妇人妻久久综合中文| 精品久久久久久电影网| 少妇的丰满在线观看| 久久久国产一区二区| 少妇猛男粗大的猛烈进出视频| 9191精品国产免费久久| 女人爽到高潮嗷嗷叫在线视频| 日韩 亚洲 欧美在线| 女人爽到高潮嗷嗷叫在线视频| 久热爱精品视频在线9| 久久精品久久久久久久性| 免费黄频网站在线观看国产| 国产精品国产av在线观看| 亚洲精品美女久久久久99蜜臀 | 欧美成人精品欧美一级黄| 韩国高清视频一区二区三区| 久久 成人 亚洲| 一本大道久久a久久精品| 18禁观看日本| 亚洲久久久国产精品| 18在线观看网站| 大片电影免费在线观看免费| 国产在视频线精品| 熟妇人妻不卡中文字幕| 久久亚洲国产成人精品v| 久久天躁狠狠躁夜夜2o2o | 桃花免费在线播放| 亚洲成人手机| 自线自在国产av| 国产在线免费精品| 成人三级做爰电影| 亚洲美女搞黄在线观看| 老鸭窝网址在线观看| 欧美日韩视频高清一区二区三区二| 18禁观看日本| 日韩大码丰满熟妇| 久久精品亚洲熟妇少妇任你| 80岁老熟妇乱子伦牲交| av片东京热男人的天堂| 美国免费a级毛片| 久久毛片免费看一区二区三区| 久久免费观看电影| 精品少妇久久久久久888优播| 国产一级毛片在线| 久久国产精品大桥未久av| 婷婷色麻豆天堂久久| 国产视频首页在线观看| 精品国产超薄肉色丝袜足j| 欧美成人精品欧美一级黄| av在线播放精品| 丰满迷人的少妇在线观看| 婷婷色综合www| 国产亚洲一区二区精品| 亚洲精品国产区一区二| 国产精品国产三级国产专区5o| 日本91视频免费播放| 熟女av电影| 90打野战视频偷拍视频| 欧美日韩一区二区视频在线观看视频在线| 国产一区二区三区综合在线观看| 满18在线观看网站| 哪个播放器可以免费观看大片| 日日撸夜夜添| 国产一区亚洲一区在线观看| 在线观看国产h片| 黄色毛片三级朝国网站| 亚洲欧美精品综合一区二区三区| 亚洲av成人不卡在线观看播放网 | 国产男女超爽视频在线观看| 亚洲伊人久久精品综合| av女优亚洲男人天堂| 赤兔流量卡办理| 免费人妻精品一区二区三区视频| 中文字幕高清在线视频| 国产xxxxx性猛交| 伊人亚洲综合成人网| 成人影院久久| 黄色 视频免费看| av视频免费观看在线观看| 午夜日本视频在线| 日韩av免费高清视频| 亚洲精品中文字幕在线视频| 在线天堂最新版资源| 亚洲国产中文字幕在线视频| 亚洲精品中文字幕在线视频| 国产精品.久久久| e午夜精品久久久久久久| 亚洲综合精品二区| 成人午夜精彩视频在线观看| 国产一区二区三区av在线| 在线观看免费高清a一片| 新久久久久国产一级毛片| 亚洲av电影在线进入| 午夜老司机福利片| 女性被躁到高潮视频| 亚洲欧美一区二区三区国产| 国产成人一区二区在线| 亚洲欧洲国产日韩| 最近2019中文字幕mv第一页| av.在线天堂| 亚洲 欧美一区二区三区| 亚洲国产看品久久| 极品人妻少妇av视频| 一区二区三区激情视频| 久久天躁狠狠躁夜夜2o2o | 老熟女久久久| 美女主播在线视频| 亚洲精品在线美女| 在线看a的网站| 亚洲欧洲精品一区二区精品久久久 | 午夜91福利影院| 午夜影院在线不卡| 日本色播在线视频| 久久精品亚洲熟妇少妇任你| 亚洲在久久综合| 精品国产一区二区久久| 女性被躁到高潮视频| 蜜桃国产av成人99| 欧美最新免费一区二区三区| 成人影院久久| 亚洲av在线观看美女高潮| 国产在线免费精品| 亚洲精品国产色婷婷电影| 亚洲熟女毛片儿| 久久精品人人爽人人爽视色| www.熟女人妻精品国产| 欧美成人午夜精品| 岛国毛片在线播放| 久久久国产一区二区| 日韩 亚洲 欧美在线| 19禁男女啪啪无遮挡网站| 日本vs欧美在线观看视频| 午夜激情av网站| 捣出白浆h1v1| 可以免费在线观看a视频的电影网站 | avwww免费| 汤姆久久久久久久影院中文字幕| 王馨瑶露胸无遮挡在线观看| 秋霞伦理黄片| 波多野结衣一区麻豆| 午夜福利视频精品| 亚洲欧美成人精品一区二区| 日韩大片免费观看网站| 成年人午夜在线观看视频| 亚洲精品久久成人aⅴ小说| 久久精品亚洲av国产电影网| 99九九在线精品视频| 国产色婷婷99| 午夜av观看不卡| 咕卡用的链子| 欧美激情 高清一区二区三区| 一区二区av电影网| 国产老妇伦熟女老妇高清| 久久久欧美国产精品| 在线观看免费日韩欧美大片| 99九九在线精品视频| avwww免费| 免费观看a级毛片全部| 亚洲欧美清纯卡通| 人妻一区二区av| 一级毛片 在线播放| 丝袜脚勾引网站| 亚洲精品美女久久av网站| 国产97色在线日韩免费| 777久久人妻少妇嫩草av网站| 国产毛片在线视频| 9色porny在线观看| 天天影视国产精品| 免费观看a级毛片全部| 国产成人a∨麻豆精品| 老司机亚洲免费影院| 美女视频免费永久观看网站| 校园人妻丝袜中文字幕| 欧美乱码精品一区二区三区| 亚洲专区中文字幕在线 | 亚洲精品久久成人aⅴ小说| 看十八女毛片水多多多| 精品国产乱码久久久久久小说| 久久综合国产亚洲精品| 久热爱精品视频在线9| 自线自在国产av| 国产在线视频一区二区| 国产一区二区三区av在线| 中国国产av一级| 麻豆乱淫一区二区| 国产野战对白在线观看| 日韩av在线免费看完整版不卡| av国产精品久久久久影院| h视频一区二区三区| 日韩熟女老妇一区二区性免费视频| 国产精品亚洲av一区麻豆 | 日韩熟女老妇一区二区性免费视频| 国产男女超爽视频在线观看| 亚洲一级一片aⅴ在线观看| 一边亲一边摸免费视频| 欧美成人午夜精品| 国产精品久久久av美女十八| 国产男女超爽视频在线观看| 久久性视频一级片| 亚洲熟女精品中文字幕| 人人妻人人添人人爽欧美一区卜| 午夜福利视频精品| 国产日韩欧美视频二区| 少妇精品久久久久久久| 制服丝袜香蕉在线| 欧美人与性动交α欧美精品济南到| 欧美日韩福利视频一区二区| 中文字幕av电影在线播放| 韩国av在线不卡| 少妇被粗大的猛进出69影院| 国产免费现黄频在线看| 在线观看免费日韩欧美大片| 校园人妻丝袜中文字幕| 晚上一个人看的免费电影| 精品人妻在线不人妻| 久久精品国产亚洲av高清一级| 中文精品一卡2卡3卡4更新| 国产老妇伦熟女老妇高清| 高清在线视频一区二区三区| 国产精品.久久久| 亚洲国产av影院在线观看| 日韩av免费高清视频| 中文字幕精品免费在线观看视频| 少妇精品久久久久久久| 亚洲激情五月婷婷啪啪| 亚洲欧洲日产国产| av在线播放精品| 妹子高潮喷水视频| 久久久久久免费高清国产稀缺| 美女视频免费永久观看网站| 色94色欧美一区二区| 欧美日韩福利视频一区二区| 久久综合国产亚洲精品| 热re99久久精品国产66热6| 老汉色av国产亚洲站长工具| 亚洲熟女精品中文字幕| 性少妇av在线| av卡一久久| 亚洲精品成人av观看孕妇| 另类亚洲欧美激情| 久久精品人人爽人人爽视色| 欧美精品一区二区免费开放| 国产亚洲最大av| 777久久人妻少妇嫩草av网站| 久久狼人影院| 亚洲精品视频女| 精品国产乱码久久久久久小说| 午夜日韩欧美国产| 18禁观看日本| 韩国精品一区二区三区| 丁香六月天网| 免费观看人在逋| 制服人妻中文乱码| 日韩成人av中文字幕在线观看| 王馨瑶露胸无遮挡在线观看| 美女午夜性视频免费| 亚洲一卡2卡3卡4卡5卡精品中文| a 毛片基地| 色94色欧美一区二区| 久久影院123| 国产一区有黄有色的免费视频| 99精品久久久久人妻精品| 国产黄色视频一区二区在线观看| 欧美最新免费一区二区三区| 久久久久久久久免费视频了| 久久婷婷青草| 人人妻人人澡人人看| 国产日韩欧美视频二区| 男的添女的下面高潮视频| 王馨瑶露胸无遮挡在线观看| 十八禁网站网址无遮挡| 欧美国产精品一级二级三级| 久久精品国产a三级三级三级| 青草久久国产| 超碰成人久久| 男女床上黄色一级片免费看| 久久ye,这里只有精品| 丝袜喷水一区| 最新的欧美精品一区二区| 韩国av在线不卡| 纵有疾风起免费观看全集完整版| 综合色丁香网| 伊人久久大香线蕉亚洲五| 国产成人精品福利久久| 中文字幕av电影在线播放| 宅男免费午夜| 亚洲精品乱久久久久久| 亚洲欧美激情在线| 国产伦人伦偷精品视频| 老司机深夜福利视频在线观看 | 国产97色在线日韩免费| 嫩草影视91久久| 午夜日本视频在线| 亚洲婷婷狠狠爱综合网| 人妻 亚洲 视频| 99久久精品国产亚洲精品| 久久精品久久久久久久性| 人体艺术视频欧美日本| 久久久久久久精品精品| 午夜福利一区二区在线看| 啦啦啦在线观看免费高清www| 乱人伦中国视频| 亚洲一卡2卡3卡4卡5卡精品中文| 久久人人爽人人片av| 纵有疾风起免费观看全集完整版| 亚洲精品国产av成人精品| 九草在线视频观看| 我的亚洲天堂| 人妻 亚洲 视频| 色综合欧美亚洲国产小说| 国产野战对白在线观看| 日韩精品有码人妻一区| 老汉色∧v一级毛片| 欧美亚洲 丝袜 人妻 在线| 男人爽女人下面视频在线观看| 韩国精品一区二区三区| 日本午夜av视频| 妹子高潮喷水视频| 欧美日韩成人在线一区二区| 韩国av在线不卡| 这个男人来自地球电影免费观看 | 久久99一区二区三区| 无限看片的www在线观看| 国产成人精品在线电影| 高清不卡的av网站| a级片在线免费高清观看视频| 精品国产国语对白av| 久久精品国产亚洲av涩爱| 成人三级做爰电影| 新久久久久国产一级毛片| 久久久久精品久久久久真实原创| 国产无遮挡羞羞视频在线观看| 亚洲精品视频女| av不卡在线播放| 国产精品av久久久久免费| 久久 成人 亚洲| 久久精品国产亚洲av涩爱| 黑人欧美特级aaaaaa片| 国产成人精品福利久久| 色婷婷久久久亚洲欧美| 香蕉国产在线看| 精品人妻一区二区三区麻豆| 一级毛片电影观看| 黄色一级大片看看| 在线 av 中文字幕| 久久性视频一级片| 男女午夜视频在线观看| 国产亚洲最大av| 伦理电影大哥的女人| 日韩免费高清中文字幕av| 国产女主播在线喷水免费视频网站| 亚洲欧美激情在线| 性少妇av在线| 另类亚洲欧美激情| 丝袜美足系列| a级毛片在线看网站| 婷婷色av中文字幕| 少妇的丰满在线观看| 18禁国产床啪视频网站| 国产1区2区3区精品| 免费女性裸体啪啪无遮挡网站| 韩国av在线不卡| 日本wwww免费看| 国产又爽黄色视频| 制服人妻中文乱码| 欧美日韩综合久久久久久| 久久久久精品性色| 精品国产乱码久久久久久小说| 国产成人欧美| 成人黄色视频免费在线看| 18禁观看日本| 欧美中文综合在线视频| 欧美老熟妇乱子伦牲交| 色视频在线一区二区三区| 亚洲国产日韩一区二区| 国产老妇伦熟女老妇高清| 亚洲色图 男人天堂 中文字幕| 午夜福利视频在线观看免费| 中国国产av一级| 伦理电影大哥的女人| 久久久久精品性色| 日本一区二区免费在线视频| 日韩视频在线欧美| 国产av国产精品国产| 亚洲一码二码三码区别大吗| 免费看av在线观看网站| 精品久久蜜臀av无| 国产伦理片在线播放av一区| 午夜91福利影院| 18在线观看网站| 国产乱人偷精品视频| av一本久久久久| 亚洲av电影在线进入| 一二三四在线观看免费中文在| 成人国产麻豆网| 日韩 欧美 亚洲 中文字幕| 国产99久久九九免费精品| 国产福利在线免费观看视频| 中文字幕亚洲精品专区| 国产成人啪精品午夜网站| 啦啦啦 在线观看视频| av不卡在线播放| 男人操女人黄网站| 精品国产一区二区三区四区第35| 久久97久久精品| 欧美精品一区二区大全| 久久久久久久大尺度免费视频| 大片电影免费在线观看免费| av有码第一页| 少妇人妻精品综合一区二区| 亚洲av男天堂| 亚洲,欧美,日韩| 久久精品亚洲熟妇少妇任你| 中文字幕人妻熟女乱码| 久热这里只有精品99| 亚洲国产欧美网| 一级,二级,三级黄色视频| 免费黄网站久久成人精品| svipshipincom国产片| 久久女婷五月综合色啪小说| 麻豆乱淫一区二区| 各种免费的搞黄视频| 人人妻人人爽人人添夜夜欢视频| 久久精品亚洲熟妇少妇任你| 最黄视频免费看| 久久久精品免费免费高清| 王馨瑶露胸无遮挡在线观看| 久久99精品国语久久久| 天天影视国产精品| 黄色毛片三级朝国网站| 亚洲精品自拍成人| 自线自在国产av| av天堂久久9| 国产有黄有色有爽视频| 老司机靠b影院| √禁漫天堂资源中文www| 国产精品一区二区在线观看99| 国产精品国产三级国产专区5o| 国产成人欧美在线观看 | 久久久久精品久久久久真实原创| 亚洲人成网站在线观看播放| 午夜福利影视在线免费观看| 欧美人与性动交α欧美精品济南到| 日韩一区二区视频免费看| 精品国产国语对白av| 两性夫妻黄色片| 久久精品熟女亚洲av麻豆精品| 丝袜美足系列| 成人手机av| 国产淫语在线视频| 看免费av毛片| 麻豆av在线久日| av在线老鸭窝| 国产1区2区3区精品| 日韩大码丰满熟妇| 久久国产精品男人的天堂亚洲| 国产1区2区3区精品| 国产精品国产三级专区第一集| 国产精品久久久av美女十八| 久久国产亚洲av麻豆专区| 免费高清在线观看视频在线观看| 国产免费又黄又爽又色| 亚洲成人手机| 91aial.com中文字幕在线观看| 久久婷婷青草| 女人爽到高潮嗷嗷叫在线视频| 日本午夜av视频| 久久天堂一区二区三区四区| 99久国产av精品国产电影| 日韩不卡一区二区三区视频在线| 国产伦理片在线播放av一区| 亚洲国产精品成人久久小说| 99热网站在线观看| 国产毛片在线视频| a 毛片基地| 成人国语在线视频| 国产精品.久久久| 在线观看www视频免费| 熟女少妇亚洲综合色aaa.| 亚洲av电影在线进入|