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

    秸稈纖維素在亞/超臨界乙醇中液化反應(yīng)路徑及5集總反應(yīng) 動力學(xué)模型的建立

    2015-07-25 09:12:32戰(zhàn)曉青解新安李璐王婭莉陶紅秀
    化工進(jìn)展 2015年11期
    關(guān)鍵詞:集總超臨界液化

    戰(zhàn)曉青,解新安,李璐,王婭莉,陶紅秀

    (華南農(nóng)業(yè)大學(xué)食品學(xué)院,廣東 廣州 510641)

    隨著人們?nèi)找嬖鲩L的能源需求以及化石燃料的枯竭,尋找可替代新型能源迫在眉睫。生物質(zhì)作為一種可再生資源因其分布廣、價格低、產(chǎn)量高而備受關(guān)注。其中,秸稈生物質(zhì)的熱化學(xué)轉(zhuǎn)化制備燃料和平臺化合物也已成為研究熱點。而生物質(zhì)熱化學(xué)轉(zhuǎn)化方法中,超臨界液化技術(shù)因為其產(chǎn)品品質(zhì)較好、產(chǎn)率高等優(yōu)點備受研究者青睞。生物質(zhì)液化產(chǎn)物種類繁多,且隨原料種類不同變化明顯,導(dǎo)致其液化機理非常復(fù)雜。Durak[1]研究了光果甘草在超臨界乙醇和丙酮中液化過程,發(fā)現(xiàn)在290℃、無催化劑、丙酮為溶劑時,生物油收率最高達(dá)到79%。Koriakin等[2]對微晶纖維素在超臨界有機溶劑中液化進(jìn)行了研究,發(fā)現(xiàn)轉(zhuǎn)化率最高可達(dá)92%(含氣體產(chǎn)物),但其溫度和壓力要求較高。廖益強等[3]對竹粉在亞臨界無水乙醇中液化研究發(fā)現(xiàn)影響液化率的反應(yīng)條件順序為:反應(yīng)溫度>反應(yīng)時間>液固質(zhì)量比>催化劑用量,在最佳條件下液化率可達(dá)86.44%。本文作者課題組研究玉米秸稈生物質(zhì)在亞/超臨界乙醇中液化過程,并建立了包含氣體、輕油、重油、易揮發(fā)組分和C、H、N、O 的8 集總反應(yīng)網(wǎng)絡(luò)[4],同時亦研究了玉米秸稈生物質(zhì)液化過程,根據(jù)溫度和乙醇用量的不同、反應(yīng)階段和升溫階段的不同分析建立了含有氣體、輕油、重油、揮發(fā)分和殘渣的5集總反應(yīng)網(wǎng)絡(luò)[5]。但以上研究均局限于玉米秸稈生物質(zhì)液化的反應(yīng)路徑和反應(yīng)網(wǎng)絡(luò),但對秸稈纖維素的液化反應(yīng)路徑研究,尤其是動力學(xué)模型的建立和動力學(xué)參數(shù)的求取則是少有研究。而王連山等[6]建立了催化重整反應(yīng)38 集總動力學(xué)模型,并在模擬計算后用該模型優(yōu)化提高了芳烴收率。王建平等[7]建立了加氫渣油催化裂化14 集總動力學(xué)模型,求取了動力學(xué)參數(shù),驗證表明模型可以較為準(zhǔn)確地預(yù)測催化裂化產(chǎn)品分布和主要產(chǎn)物性質(zhì)。本文以集總方法研究玉米秸稈纖維素在亞/超臨界乙醇中液化過程,建立其液化反應(yīng)路徑,以此建立反應(yīng)網(wǎng)絡(luò),并進(jìn)行動力學(xué)模擬,簡化和優(yōu)化動力學(xué)模型,求取動力學(xué)參數(shù),為將來生物質(zhì)的應(yīng)用和液化反應(yīng)器的設(shè)計提供理論依據(jù)。

    1 反應(yīng)路徑與反應(yīng)網(wǎng)絡(luò)

    1.1 集總劃分

    秸稈纖維素在亞/超臨界乙醇中液化過程非常復(fù)雜,包含大量的平行反應(yīng)和交叉反應(yīng),且產(chǎn)物種類繁多,包含大量酮類、酯類、醇類和酸類等化合物[8-10]。因此,為研究液化過程動力學(xué),必須在保證研究的準(zhǔn)確性和可行性的基礎(chǔ)上,盡可能地使之簡化,所以采用集總的方法,根據(jù)本文作者課題組秸稈纖維素在亞/超臨界乙醇中液化實驗和機理分析研究結(jié)果,并考慮液化產(chǎn)物結(jié)構(gòu)、分子量、產(chǎn)物特性等[11],將液化產(chǎn)物劃分為5 個集總,分別為秸稈纖維素在反應(yīng)釜液化后,打開反應(yīng)釜通氣閥,收集氣體(GAS),打開反應(yīng)釜,用無水乙醇將液化產(chǎn)物轉(zhuǎn)移至燒杯,置于室溫條件待乙醇自然揮發(fā)(揮發(fā)部分為乙醇+VO),加入蒸餾水提取水溶性有機物即輕油(WSO);剩余固體物再用丙酮提取,提取物為重油(HO),不溶物即殘渣(RE),數(shù)據(jù)如表1[12-13](原料量15g,乙醇160mL),并分析液化過程,建立反應(yīng)路徑,形成反應(yīng)網(wǎng)絡(luò),進(jìn)而研究其反應(yīng)動力學(xué)[5,11,14-15]。

    表1 不同溫度及反應(yīng)停留時間下各集總收率

    1.2 反應(yīng)路徑

    纖維素液化過程中,輕油收率隨溫度升高而提高,且隨反應(yīng)停留時間增加而減少。這可能是由于纖維素在較低溫度時發(fā)生脫水、裂解等形成小分子片段,而隨著溫度提高,小分子片段活性增強,導(dǎo)致其相互之間發(fā)生縮合、聚合等反應(yīng)形成氣體和重油[4-5,16]。而乙醇自由基的存在,穩(wěn)定了中間產(chǎn)物,抑制了小分子片段之間的相互作用,減少了其發(fā)生縮合形成殘渣,從而使得輕油收率隨反應(yīng)停留時間增加而減少得較為緩慢,見式(1)~式(3)。

    重油收率隨溫度升高而逐漸增多,且隨反應(yīng)停留時間增加而提高。這可能是由于溫度升高,乙醇產(chǎn)生的自由基濃度和活性提高,使其能夠促進(jìn)纖維素裂解進(jìn)程,從而形成重油。與此同時,乙醇自由基可能作用于液化中間產(chǎn)物,使液化中間產(chǎn)物穩(wěn)定,抑制了中間產(chǎn)物之間重聚形成殘渣。同時,乙醇自由基也可能與小分子中間產(chǎn)物轉(zhuǎn)化為分子量較大,碳鏈較長的重油成分,見式(4)[4,15]。

    揮發(fā)分與輕油類似,收率均是隨反應(yīng)溫度升高和反應(yīng)停留時間增加而提高。這可能是由于揮發(fā)分中多為分子量較小、沸點較低化合物組成,其隨著溫度升高,易于發(fā)生相互反應(yīng)生成碳鏈相對較長的輕油。同時,一些分子量相對較大的化合物也可能發(fā)生加氫裂化生成分子量較小的揮發(fā)分,見式(5)、式(6)[17]。

    氣體收率隨反應(yīng)停留時間增加而逐漸提高,且氣體收率隨著溫度升高而增加,這可能是由于溫度升高,乙醇自由基活性提高,首先,其能夠促進(jìn)纖維素裂解形成小分子化合物導(dǎo)致氣體增多,繼而乙醇自由基亦可以作用于液化中間產(chǎn)物,促進(jìn)其加氫裂化形成小分子氣體,同時,纖維素裂解產(chǎn)生的小分子氣體可能會穩(wěn)定性較低,在乙醇自由基的作用下發(fā)生聚合反應(yīng)形成輕油,見式(7)[4]。

    殘渣收率逐漸降低,表明纖維素斷鍵裂解形成各種產(chǎn)物,在280℃時,殘渣收率降低較快,當(dāng)溫度達(dá)到300℃時,反應(yīng)停留時間對殘渣收率的影響變得較為平緩,這可能是因為當(dāng)溫度較高時,纖維素已轉(zhuǎn)化較為完全,而隨著反應(yīng)停留時間的提高,發(fā)生反應(yīng)更多的是產(chǎn)物之間的相互轉(zhuǎn)化。而殘渣的形成,可能是由于乙醇自由基的存在,使液化中間產(chǎn)物穩(wěn)定,抑制了中間產(chǎn)物之間聚合形成殘渣,見式(8)~式(11)[5,11]。

    1.3 反應(yīng)網(wǎng)絡(luò)

    綜合上述反應(yīng)路徑,可以得到秸稈纖維素超臨界液化的反應(yīng)網(wǎng)絡(luò)如圖1。

    由圖1 可以看出,秸稈纖維素超臨界液化過程中存在諸多交叉反應(yīng),產(chǎn)物之間相互轉(zhuǎn)化較多,這也是纖維素液化機理非常復(fù)雜的重要原因。輕油、重油(即生物油)是纖維素液化的主要產(chǎn)物,將其作為兩大集總來研究,既大幅降低了研究的難度, 又對纖維素實際利用具有良好的指導(dǎo)意義,基于此網(wǎng)絡(luò)秸稈纖維素液化動力學(xué)研究也可以為將來秸稈生物質(zhì)液化的大規(guī)模應(yīng)用提供理論依據(jù)。

    圖1 秸稈纖維素液化5 集總反應(yīng)網(wǎng)絡(luò)

    2 動力學(xué)

    2.1 動力學(xué)方程

    根據(jù)前文建立的反應(yīng)路徑和反應(yīng)網(wǎng)絡(luò),結(jié)合秸稈纖維素超臨界液化實際來建立動力學(xué)方程。同時,由于纖維素液化過程復(fù)雜,平行和交叉反應(yīng)眾多,所以應(yīng)盡可能簡化模型,降低計算難度,故動力學(xué)方程的建立不考慮環(huán)境因素和操作因素。本實驗為質(zhì)量守恒體系,且為封閉系統(tǒng),故僅考慮反應(yīng)本身的影響,動力學(xué)方程按式(12)建立。

    式中,cM為集總M 的濃度;t 為反應(yīng)停留時間; Ni為與M 相關(guān)的第i 個集總;ki為反應(yīng)i 的反應(yīng)速率常數(shù);nmi、nni為對應(yīng)的反應(yīng)級數(shù);j 為與M 相關(guān)的集總數(shù)。

    因此,初始動力學(xué)方程建立如式(13)~式(17)。

    當(dāng)上述方程中各參數(shù)求解出后,將其代入阿倫尼烏斯公式[式(18)]。

    式中,k 為反應(yīng)速率常數(shù);A 為指前因子;Ea為活化能;R 為通用氣體常數(shù);T 為反應(yīng)溫度。

    然后根據(jù)不同溫度下反應(yīng)速率常數(shù)k 來求取各反應(yīng)活化能Ea和指前因子A,最終用實驗數(shù)據(jù)驗證模型的準(zhǔn)確性。

    反應(yīng)級數(shù)、指前因子和活化能是動力學(xué)模型的主要參數(shù),反應(yīng)級數(shù)能夠直觀地反映各個集總在反應(yīng)過程中的權(quán)重,即反應(yīng)級數(shù)與集總收率對反應(yīng)速率影響成正比?;罨芎椭盖耙蜃觿t反映反應(yīng)體系中各組分濃度和溫度與反應(yīng)速率之間的關(guān)系,以滿足反應(yīng)過程開發(fā)和反應(yīng)器設(shè)計的需要。

    2.2 參數(shù)求解和優(yōu)化

    考慮到反應(yīng)路徑復(fù)雜、各集總之間相互作用較多,本工作以遺傳算法為基礎(chǔ)求取和優(yōu)化動力學(xué)參數(shù)[18]。遺傳算法是模擬自然界遺傳機制和進(jìn)化論來搜索最優(yōu)解的算法,由于其多次迭代的方式和良好的魯棒性,以及多種群平行搜索的特點,使其一般能夠獲得較優(yōu)解。但遺傳算法全局搜索能力較強而局部搜索能力較弱,這導(dǎo)致遺傳算法往往會得到全局較優(yōu)解而非最優(yōu)解,非線性規(guī)劃恰恰可以彌補這一點,而最小二乘法操作簡單,適用性強,在已有模型基礎(chǔ)上,非常適合用來進(jìn)行參數(shù)估計和局部優(yōu)化,將最小二乘法與遺傳算法結(jié)合相對容易實現(xiàn),且經(jīng)過測試發(fā)現(xiàn)與單獨使用遺傳算法相比,其結(jié)果較優(yōu)。故本文采用遺傳算法結(jié)合最小二乘法進(jìn)行非線性規(guī)劃來求取動力學(xué)參數(shù)和優(yōu)化模型[19-22]。

    2.2.1 遺傳算法

    遺傳算法的基本單位是個體,根據(jù)要解決的問題需要將參數(shù)編碼而成,然后根據(jù)適應(yīng)度函數(shù)來進(jìn)化。之后經(jīng)過選擇、交叉、變異來獲得適應(yīng)度好的新的種群,如此循環(huán)。高適應(yīng)度的染色體更容易被選擇進(jìn)入下一代,常用的選擇方法有輪盤賭法、隨機排序等。交叉方法也很多,可以根據(jù)問題需要來設(shè)置,比如固定位置基因的交叉、均勻交叉等。變異則可以維持種群的多樣性,防止因早熟而無法獲得最優(yōu)解[19,23]。經(jīng)過以上幾步,新的個體已經(jīng)產(chǎn)生,然后投入新的循環(huán),直到達(dá)到預(yù)先設(shè)定的遺傳代數(shù)。

    2.2.2 遺傳算法應(yīng)用

    遺傳算法參數(shù)設(shè)置如下:遺傳代數(shù)為200,初始種群規(guī)模為20,每個染色體基因數(shù)為33,分別對應(yīng)k1~k11和n1~n22,交叉概率為0.6,指定染色體段交叉,即k 與k 交叉、n 與n 交叉而k 與n 不交叉。變異概率為0.01。適應(yīng)度函數(shù)為各模擬值與實驗點差值的絕對值之和。為了提高解的精度,將非線性規(guī)劃融入循環(huán)中,每10 代進(jìn)行一次非線性規(guī)劃,以優(yōu)化局部解。具體流程如圖2[24]。

    圖2 遺傳算法結(jié)合非線性規(guī)劃流程圖

    計算結(jié)果得出后,在保持準(zhǔn)確度前提下,盡量簡化方程形式,得出優(yōu)化結(jié)果后代入阿倫尼烏斯公式求取Ea和A,整理如表2。

    表2 遺傳算法優(yōu)化后各動力學(xué)參數(shù)

    在反應(yīng)級數(shù)選擇過程中,考慮到模型復(fù)雜、交叉反應(yīng)多,為了提高模型準(zhǔn)確性,不考慮其按照一級反應(yīng)進(jìn)行,而將其范圍設(shè)置符合常規(guī),并通過調(diào)節(jié)使其更加合理。由表2 可以看出,求得的反應(yīng)級數(shù)最大為1.6,對應(yīng)GAS 集總,表明在反應(yīng)過程中GAS 收率對反應(yīng)速率影響較大,而相對應(yīng)的最小反應(yīng)級數(shù)為0.2,則表明RE 收率對反應(yīng)速率影響較小。

    將模擬優(yōu)化結(jié)果代入動力學(xué)方程中,整理后可得優(yōu)化后動力學(xué)方程,如式(19)~式(23)。

    從上述方程可以看出,方程已經(jīng)大幅簡化,從而降低了方程的復(fù)雜程度。反應(yīng)級數(shù)范圍在0.2~1.6 之間,活化能在88.55~168.53kJ/mol 之間,同時從活化能可以看出,纖維素裂解所需活化能較大,而各集總之間交叉反應(yīng)的活化能相對較小,也表明隨著反應(yīng)的進(jìn)行,纖維素液化過程中更多的是各中間產(chǎn)物之間的相互轉(zhuǎn)化。

    3 模型驗證

    3.1 模擬效果

    模型已經(jīng)計算和優(yōu)化完成,然后用所建立的動力學(xué)對液化過程進(jìn)行模擬計算,將模擬值與實驗值作對比,以驗證模型準(zhǔn)確性,結(jié)果如圖3 所示。

    圖3 280℃下各集總收率模型計算值與實驗值對比

    由圖3 可以看出,280℃下動力學(xué)模型準(zhǔn)確性較好,但RE 值與實驗值相差較大,分析認(rèn)為,這是由于280℃時,纖維素剩余量較大,纖維素液化仍處于較高速度狀態(tài),所以反應(yīng)停留時間越長,RE實驗值降低越快,而模擬值則是呈現(xiàn)出RE 收率降低越來越慢,這可能是由于選取的實驗數(shù)據(jù)溫度范圍為280~320℃,當(dāng)溫度高于280℃后,RE 降速放緩,而使得整個模型在RE 收率變化呈現(xiàn)同高溫時相同的趨勢。

    圖4 顯示了300℃下的模型求解值與實驗值對比,由圖4 可以看出模型準(zhǔn)確性好,且RE 變化與280℃下相比跟實驗數(shù)據(jù)更為接近,尤其是主要產(chǎn)物輕油和重油,模擬值與實驗值擬合較好。

    從圖5 可以看出,主要產(chǎn)品輕油和重油模擬結(jié)果與實驗值非常接近,且能看出RE 變化起初較快,反應(yīng)停留時間增長,其變化也變得平緩,也說明了此時纖維素已經(jīng)液化較為徹底。氣體和殘渣的模擬結(jié)果則略有出入,這可能是由于秸稈纖維素液化過程屬于一個封閉系統(tǒng),其模擬的可調(diào)性較差,從而使得氣體的相對誤差略高。

    3.2 誤差分析

    從圖3~圖5 可以粗略看出,動力學(xué)模型較為可信,現(xiàn)對模擬結(jié)果進(jìn)行誤差分析,以更精確地反映建立的動力學(xué)模型的準(zhǔn)確度,各集總實驗值與模擬值相關(guān)系數(shù)見表3。

    由表3 可以看出,相關(guān)系數(shù)范圍在0.9345~0.9995 之間,且大部分集總實驗值與模擬計算值相關(guān)系數(shù)在0.99 以上,由此可見,文中建立的動力學(xué)模型能夠較好地反映秸稈纖維素在亞/超臨界乙醇中液化的過程,所獲得的各個動力學(xué)參數(shù)較為 可靠。

    圖4 300℃下各集總收率模型計算值與實驗值對比

    圖5 320℃下各集總收率模型計算值與實驗值對比

    表3 實驗值與模擬計算值相關(guān)系數(shù)表

    4 結(jié) 論

    綜上所述,本文系統(tǒng)研究了秸稈纖維素在亞/超臨界乙醇中液化過程,建立了含輕油(WSO)、重油(HO)、揮發(fā)分(VO)、氣體(GAS)和殘渣(RE)的5 集總反應(yīng)網(wǎng)絡(luò),并建立了動力學(xué)模型。結(jié)論如下。

    (1)求出了反應(yīng)動力學(xué)常數(shù)包括活化能(Ea)范圍在88.55~168.53kJ/mol、指前因子(A)范圍在6.27×105~7.28×1013min-1以及反應(yīng)級數(shù)(n)范圍在0.2~1.6。

    (2)利用建立的動力學(xué)模型對反應(yīng)過程進(jìn)行模擬,并用實驗數(shù)據(jù)進(jìn)行驗證,其相關(guān)系數(shù)在0.9345~0.9995 之間,證明所建立的動力學(xué)模型是可行的。

    (3)采用遺傳算法結(jié)合非線性規(guī)劃對秸稈纖維素在亞/超臨界乙醇中液化過程的動力學(xué)參數(shù)進(jìn)行估值,其求解策略是可行的。

    [1] Durak H. Bio-oil production from Glycyrrhiza glabra through supercritical fluid extraction[J]. The Journal of Supercritical Fluids,2014,95:373-386.

    [2] Koriakin A,Van Nguyen H,Kim D W,et al. Direct thermochemical liquefaction of microcrystalline cellulose by sub- and supercritical organic solvents[J]. The Journal of Supercritical Fluids,2014,95:175-186.

    [3] 廖益強,郭銀清,盧澤湘,等. 竹粉乙醇液化及其產(chǎn)物表征[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報,2014,19(2):43-50.

    [4] Liu H,Xie X,Ren J,et al. 8-Lump reaction pathways of cornstalk liquefaction in sub- and super-critical ethanol[J]. Industrial Crops and Products,2012,35(1):250-256.

    [5] 劉華敏,解新安,丁年平,等. 亞/超臨界乙醇液化玉米秸稈反應(yīng)路徑與機理[J]. 農(nóng)業(yè)工程學(xué)報,2010,26(6):277-282.

    [6] 王連山,張泉靈,梁超. 催化重整反應(yīng)38 集總動力學(xué)模型及其在連續(xù)催化重整中的應(yīng)用[J]. 化工學(xué)報,2012,63(4):1076-1082.

    [7] 王建平,許先焜,翁惠新,等. 加氫渣油催化裂化14 集總動力學(xué)模型的建立[J]. 化工學(xué)報,2007,58(1):86-94.

    [8] Hashimoto K,Hasegawa I,Hayashi J,et al. Correlations of kinetic parameters in biomass pyrolysis with solid residue yield and lignin content[J]. Fuel,2011,90(1):104-112.

    [9] 王昶,李丹,郝慶蘭,等. 粉粒流化床中松木生物質(zhì)熱解特性的研究[J]. 燃料化學(xué)學(xué)報,2012,40(2):156-163.

    [10] 武宏香,李海濱,趙增立. 煤與生物質(zhì)熱重分析及動力學(xué)研究[J]. 燃料化學(xué)學(xué)報,2009,37(5):538-545.

    [11] Zheng C,Xie X,Tao H. Distribution and characterizations of liquefaction of celluloses in sub- and super-critical ethanol[J]. BioResources,2013,8(1):648-662.

    [12] 鄭朝陽. 秸稈纖維素在亞-超臨界乙醇中液化行為與機理分析[D]. 廣州:華南農(nóng)業(yè)大學(xué),2012.

    [13] 陶紅秀. 亞/超臨界乙醇對纖維素液化及酮類、酯類化合物形成機理的研究[D]. 廣州:華南農(nóng)業(yè)大學(xué),2013.

    [14] 鄭朝陽,解新安,陶紅秀,等. 亞/超臨界乙醇液化秸稈纖維素解聚反應(yīng)研究與機理初探[J]. 燃料化學(xué)學(xué)報,2012,40(5):526-532.

    [15] 陶紅秀,解新安,鄭朝陽,等. 玉米秸稈纖維素在亞/超臨界乙醇中的液化行為研究[J]. 西北農(nóng)林科技大學(xué)學(xué)報:自然科學(xué)版,2014(1):196-204.

    [16] 向波濤,王濤,沈忠耀. 含乙醇廢水的超臨界水氧化反應(yīng)動力學(xué)及反應(yīng)機理[J]. 化工學(xué)報,2003,64(1):80-85.

    [17] Wang G,Li W,Li B,et al. Direct liquefaction of sawdust under syngas[J]. Fuel,2007,86(10/11):1587-1593.

    [18] Holland J H. Adaptation in Natural and Artificial Systems[M]. Detroit MI:University of Michigan Press,1975.

    [19] Homayouni S M,Tang S H,Motlagh O. A genetic algorithm for optimization of integrated scheduling of cranes,vehicles,and storage platforms at automated container terminals[J]. Journal of Computational and Applied Mathematics,2014,270:545-556.

    [20] Lu H,Wen X,Lan L,et al. A self-adaptive genetic algorithm to estimate JA model parameters considering minor loops [J]. Journal of Magnetism and Magnetic Materials,2015,374:502-507.

    [21] Camacho-Rodríguez J,Cerón-García M C,F(xiàn)ernández-Sevilla J M,et al. Genetic algorithm for the medium optimization of the microalga Nannochloropsis gaditana cultured to aquaculture [J]. Bioresource Technology C,2015,177:102-109.

    [22] 覃柏英. 非線性規(guī)劃的遺傳算法在多峰函數(shù)優(yōu)化中的應(yīng)用[J]. 廣西工學(xué)院學(xué)報,2013,24(2):25-31.

    [23] 涂啟玉,梅亞東. 遺傳算法在水庫(群)優(yōu)化調(diào)度研究中的應(yīng)用綜述[J]. 水電自動化與大壩監(jiān)測,2008,32(1):16-18.

    [24] 史峰,王輝,胡斐,等. MATLAB 智能算法30 個案例分析[M]. 北京:北京航空航天大學(xué)出版社,2011.

    猜你喜歡
    集總超臨界液化
    基于撕裂法的變壓器繞組集總參數(shù)等效電路頻率響應(yīng)計算方法
    超臨界CO2在頁巖氣開發(fā)中的應(yīng)用研究進(jìn)展
    云南化工(2021年5期)2021-12-21 07:41:20
    蠟油加氫脫硫集總動力學(xué)模型研究
    辨析汽化和液化
    面部液化隨意改變表情
    煤的液化原理及應(yīng)用現(xiàn)狀
    600MW超臨界機組熱經(jīng)濟性定量分析
    1200MW等級超超臨界機組可行性研究
    天然氣液化廠不達(dá)產(chǎn)的改進(jìn)方案
    一種加載集總器件的可調(diào)三維周期結(jié)構(gòu)
    中文资源天堂在线| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品国产一区二区精华液| 亚洲精品av麻豆狂野| 国产亚洲av高清不卡| 亚洲乱码一区二区免费版| 国产黄a三级三级三级人| 桃红色精品国产亚洲av| 91老司机精品| 精品不卡国产一区二区三区| 精品乱码久久久久久99久播| 欧美zozozo另类| 国产成人av教育| 国产99久久九九免费精品| 国产精品久久久久久久电影 | 亚洲真实伦在线观看| 中文字幕人妻丝袜一区二区| 别揉我奶头~嗯~啊~动态视频| 99热这里只有精品一区 | 亚洲精品在线美女| 老汉色av国产亚洲站长工具| 制服诱惑二区| 蜜桃久久精品国产亚洲av| 欧美日本视频| 五月伊人婷婷丁香| av超薄肉色丝袜交足视频| 午夜福利成人在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区在线观看日韩 | 美女午夜性视频免费| 91av网站免费观看| 欧美极品一区二区三区四区| 国产精品1区2区在线观看.| 国产又色又爽无遮挡免费看| 日韩av在线大香蕉| 伦理电影免费视频| 又紧又爽又黄一区二区| 日韩欧美三级三区| 国产亚洲av嫩草精品影院| 最新美女视频免费是黄的| 国产精品1区2区在线观看.| 欧美黑人精品巨大| 麻豆av在线久日| www.精华液| 久久久久国产精品人妻aⅴ院| 成人三级黄色视频| 亚洲五月天丁香| 亚洲男人的天堂狠狠| 无限看片的www在线观看| 欧美日本亚洲视频在线播放| 大型av网站在线播放| 99久久综合精品五月天人人| 国产精品1区2区在线观看.| 男女做爰动态图高潮gif福利片| 亚洲av熟女| 老司机福利观看| 无限看片的www在线观看| 午夜福利18| 俄罗斯特黄特色一大片| 欧美一级毛片孕妇| 熟女电影av网| 亚洲精品粉嫩美女一区| 欧美中文日本在线观看视频| 少妇熟女aⅴ在线视频| 亚洲国产精品sss在线观看| 精品国产乱码久久久久久男人| 亚洲一区高清亚洲精品| 亚洲av五月六月丁香网| 国产蜜桃级精品一区二区三区| 老鸭窝网址在线观看| 精品久久久久久成人av| 一边摸一边做爽爽视频免费| 变态另类丝袜制服| 国产区一区二久久| 国产精品1区2区在线观看.| 国产精品久久久久久精品电影| 久久亚洲真实| 高潮久久久久久久久久久不卡| 色哟哟哟哟哟哟| 十八禁网站免费在线| 精品欧美国产一区二区三| 怎么达到女性高潮| 亚洲自拍偷在线| 一夜夜www| 精品午夜福利视频在线观看一区| 变态另类成人亚洲欧美熟女| 久久香蕉国产精品| 免费看日本二区| 色老头精品视频在线观看| av在线播放免费不卡| 色精品久久人妻99蜜桃| 欧美一级a爱片免费观看看 | 在线播放国产精品三级| 久久久久久人人人人人| 免费看美女性在线毛片视频| 亚洲第一电影网av| 亚洲成人久久性| 午夜福利免费观看在线| 午夜福利成人在线免费观看| 精品不卡国产一区二区三区| 麻豆国产97在线/欧美 | 亚洲精品一区av在线观看| 叶爱在线成人免费视频播放| 日本三级黄在线观看| 午夜两性在线视频| 舔av片在线| 久久精品国产99精品国产亚洲性色| 欧美精品啪啪一区二区三区| 深夜精品福利| 亚洲免费av在线视频| 欧美性长视频在线观看| 深夜精品福利| 亚洲成a人片在线一区二区| 欧美 亚洲 国产 日韩一| 欧美在线黄色| 欧美人与性动交α欧美精品济南到| 啦啦啦观看免费观看视频高清| ponron亚洲| 精品久久久久久久毛片微露脸| 国产亚洲欧美98| 成人高潮视频无遮挡免费网站| 黄色女人牲交| 热99re8久久精品国产| 国产亚洲精品一区二区www| 精华霜和精华液先用哪个| cao死你这个sao货| 成人一区二区视频在线观看| 一级a爱片免费观看的视频| 国产精品久久久久久人妻精品电影| 丰满的人妻完整版| 99久久国产精品久久久| 亚洲av五月六月丁香网| 国产精品九九99| 99久久精品国产亚洲精品| 超碰成人久久| 欧美黑人精品巨大| 美女 人体艺术 gogo| 国产伦一二天堂av在线观看| 久久中文看片网| 国产成人aa在线观看| 亚洲一区二区三区不卡视频| 国产私拍福利视频在线观看| 99热6这里只有精品| 久久精品亚洲精品国产色婷小说| 亚洲成人国产一区在线观看| 又爽又黄无遮挡网站| 午夜成年电影在线免费观看| 91av网站免费观看| 老司机午夜福利在线观看视频| 国内精品一区二区在线观看| 欧美黑人精品巨大| 亚洲七黄色美女视频| 精品欧美国产一区二区三| 丰满人妻熟妇乱又伦精品不卡| 亚洲av五月六月丁香网| 日本 av在线| 九九热线精品视视频播放| 亚洲国产欧美一区二区综合| 丁香六月欧美| av在线播放免费不卡| 成人手机av| 麻豆一二三区av精品| 啦啦啦免费观看视频1| 这个男人来自地球电影免费观看| 国产成人精品久久二区二区免费| 国产精品香港三级国产av潘金莲| 国产成人av教育| 又黄又爽又免费观看的视频| 国产一区二区在线av高清观看| e午夜精品久久久久久久| 亚洲av五月六月丁香网| 亚洲人成77777在线视频| 亚洲精品国产一区二区精华液| 亚洲中文字幕一区二区三区有码在线看 | 一区福利在线观看| 狂野欧美白嫩少妇大欣赏| www日本在线高清视频| 宅男免费午夜| 99久久综合精品五月天人人| 天天躁夜夜躁狠狠躁躁| 搡老妇女老女人老熟妇| 国产激情偷乱视频一区二区| 少妇被粗大的猛进出69影院| 欧美zozozo另类| 亚洲色图 男人天堂 中文字幕| 丰满人妻熟妇乱又伦精品不卡| 久久香蕉精品热| 美女黄网站色视频| 99久久无色码亚洲精品果冻| 亚洲欧美激情综合另类| 男男h啪啪无遮挡| 亚洲 国产 在线| 淫秽高清视频在线观看| 一级作爱视频免费观看| 国产成人精品久久二区二区免费| 日韩欧美国产一区二区入口| 一个人免费在线观看的高清视频| 女人被狂操c到高潮| 成人国产综合亚洲| 两个人免费观看高清视频| 极品教师在线免费播放| 亚洲精华国产精华精| 蜜桃久久精品国产亚洲av| 国语自产精品视频在线第100页| 床上黄色一级片| 校园春色视频在线观看| 男人舔奶头视频| 1024视频免费在线观看| 不卡av一区二区三区| 日本精品一区二区三区蜜桃| 欧洲精品卡2卡3卡4卡5卡区| 村上凉子中文字幕在线| 最近在线观看免费完整版| 又黄又爽又免费观看的视频| 国产精品精品国产色婷婷| 欧美大码av| 婷婷精品国产亚洲av在线| www.www免费av| 91九色精品人成在线观看| 91麻豆精品激情在线观看国产| 国产成人一区二区三区免费视频网站| 国产乱人伦免费视频| 欧美性猛交黑人性爽| 黄色丝袜av网址大全| 国产精品久久久久久亚洲av鲁大| 中国美女看黄片| 欧美高清成人免费视频www| 国产av一区二区精品久久| 精品无人区乱码1区二区| 女同久久另类99精品国产91| 制服丝袜大香蕉在线| 男女那种视频在线观看| 亚洲成人中文字幕在线播放| 美女 人体艺术 gogo| 黑人巨大精品欧美一区二区mp4| 国产一区二区三区在线臀色熟女| 日韩大尺度精品在线看网址| 免费av毛片视频| 非洲黑人性xxxx精品又粗又长| 精品不卡国产一区二区三区| 人成视频在线观看免费观看| 好看av亚洲va欧美ⅴa在| 亚洲av美国av| 欧美一级毛片孕妇| 国产亚洲精品第一综合不卡| 18禁黄网站禁片午夜丰满| 三级男女做爰猛烈吃奶摸视频| 欧美性猛交╳xxx乱大交人| 97超级碰碰碰精品色视频在线观看| 欧美绝顶高潮抽搐喷水| 亚洲性夜色夜夜综合| 精品久久久久久久久久免费视频| 久久99热这里只有精品18| 给我免费播放毛片高清在线观看| 激情在线观看视频在线高清| 性色av乱码一区二区三区2| 日韩大尺度精品在线看网址| 黄片小视频在线播放| 亚洲第一欧美日韩一区二区三区| 99国产极品粉嫩在线观看| 欧美乱妇无乱码| 欧美激情久久久久久爽电影| 性色av乱码一区二区三区2| 久久 成人 亚洲| 一进一出好大好爽视频| 最近视频中文字幕2019在线8| 精品久久久久久成人av| 亚洲乱码一区二区免费版| 久久香蕉激情| 久久婷婷成人综合色麻豆| 久久精品国产亚洲av高清一级| √禁漫天堂资源中文www| 欧洲精品卡2卡3卡4卡5卡区| 无人区码免费观看不卡| 国产高清videossex| 亚洲国产精品合色在线| 色精品久久人妻99蜜桃| 欧美国产日韩亚洲一区| 熟女电影av网| 国产欧美日韩精品亚洲av| 亚洲人成电影免费在线| 成人欧美大片| 日韩欧美国产一区二区入口| 欧美日韩乱码在线| 法律面前人人平等表现在哪些方面| 在线a可以看的网站| 免费高清视频大片| 精品无人区乱码1区二区| 亚洲午夜精品一区,二区,三区| 99久久精品热视频| 琪琪午夜伦伦电影理论片6080| 中文字幕人成人乱码亚洲影| 俄罗斯特黄特色一大片| 狂野欧美激情性xxxx| 黑人操中国人逼视频| 午夜免费成人在线视频| а√天堂www在线а√下载| 午夜成年电影在线免费观看| 在线观看www视频免费| 欧美日韩中文字幕国产精品一区二区三区| 成年免费大片在线观看| 手机成人av网站| 99热只有精品国产| 99热6这里只有精品| 欧美精品啪啪一区二区三区| 18禁国产床啪视频网站| 午夜老司机福利片| 国产精品乱码一区二三区的特点| 国产91精品成人一区二区三区| 日韩高清综合在线| 亚洲av五月六月丁香网| 性色av乱码一区二区三区2| 国产亚洲精品综合一区在线观看 | 可以免费在线观看a视频的电影网站| 亚洲全国av大片| 在线永久观看黄色视频| 91麻豆精品激情在线观看国产| 欧美大码av| 免费高清视频大片| 免费高清视频大片| 欧美在线黄色| 免费av毛片视频| 国产精品综合久久久久久久免费| 中文字幕久久专区| 久久婷婷人人爽人人干人人爱| 国产精品九九99| 欧美日韩精品网址| 国产亚洲av高清不卡| 日韩精品青青久久久久久| 国产av在哪里看| 美女午夜性视频免费| 亚洲九九香蕉| 嫁个100分男人电影在线观看| 久久香蕉激情| 亚洲第一欧美日韩一区二区三区| 欧美一级a爱片免费观看看 | 国产欧美日韩一区二区精品| 国产1区2区3区精品| 久久亚洲精品不卡| 婷婷丁香在线五月| 婷婷精品国产亚洲av在线| 在线观看免费日韩欧美大片| 精品少妇一区二区三区视频日本电影| 搡老妇女老女人老熟妇| 亚洲专区字幕在线| 最近视频中文字幕2019在线8| 啦啦啦观看免费观看视频高清| 国产高清激情床上av| 可以免费在线观看a视频的电影网站| 国产1区2区3区精品| 成人欧美大片| 欧美人与性动交α欧美精品济南到| 亚洲黑人精品在线| 操出白浆在线播放| 可以免费在线观看a视频的电影网站| av片东京热男人的天堂| 欧美最黄视频在线播放免费| 草草在线视频免费看| 欧美黄色淫秽网站| 村上凉子中文字幕在线| 国产区一区二久久| 国产成人系列免费观看| 少妇的丰满在线观看| 午夜免费成人在线视频| 亚洲人与动物交配视频| 国产爱豆传媒在线观看 | 成年免费大片在线观看| 成人国产一区最新在线观看| 国产成人av教育| 日韩欧美 国产精品| 亚洲一卡2卡3卡4卡5卡精品中文| 国产99白浆流出| 精华霜和精华液先用哪个| 长腿黑丝高跟| 岛国视频午夜一区免费看| 亚洲最大成人中文| 欧美最黄视频在线播放免费| 欧美极品一区二区三区四区| 国产av一区在线观看免费| 欧洲精品卡2卡3卡4卡5卡区| 久久这里只有精品中国| 日韩欧美一区二区三区在线观看| 香蕉av资源在线| 免费av毛片视频| 国产精品av视频在线免费观看| 亚洲av美国av| 999久久久国产精品视频| 亚洲av成人一区二区三| 少妇熟女aⅴ在线视频| 成年版毛片免费区| 小说图片视频综合网站| 国产精品一区二区精品视频观看| 每晚都被弄得嗷嗷叫到高潮| 五月玫瑰六月丁香| 国内精品久久久久久久电影| 亚洲最大成人中文| 波多野结衣高清作品| 久久久国产成人精品二区| 三级男女做爰猛烈吃奶摸视频| 亚洲精品粉嫩美女一区| 狂野欧美激情性xxxx| 亚洲电影在线观看av| 少妇裸体淫交视频免费看高清 | 麻豆一二三区av精品| 日韩欧美精品v在线| 国产日本99.免费观看| av中文乱码字幕在线| 亚洲成a人片在线一区二区| 国产三级黄色录像| netflix在线观看网站| 亚洲无线在线观看| 亚洲一区二区三区不卡视频| 两个人视频免费观看高清| 国内精品一区二区在线观看| 久久中文字幕一级| 亚洲avbb在线观看| 欧美一级毛片孕妇| av视频在线观看入口| 99久久国产精品久久久| 一级毛片精品| 午夜福利欧美成人| 在线观看免费日韩欧美大片| 国产伦在线观看视频一区| 国产免费av片在线观看野外av| 女警被强在线播放| 国产乱人伦免费视频| 久久精品综合一区二区三区| 国产精品久久电影中文字幕| 国产高清视频在线播放一区| 最近最新中文字幕大全免费视频| 成人一区二区视频在线观看| 国产真人三级小视频在线观看| 99在线视频只有这里精品首页| 国产精品国产高清国产av| 久久久久久久午夜电影| 日本三级黄在线观看| 亚洲最大成人中文| 一级毛片高清免费大全| 亚洲无线在线观看| 亚洲最大成人中文| 高清毛片免费观看视频网站| 国内少妇人妻偷人精品xxx网站 | 天堂影院成人在线观看| 亚洲国产欧洲综合997久久,| 中文字幕高清在线视频| 久久中文字幕人妻熟女| 日韩精品中文字幕看吧| 久久性视频一级片| 人人妻,人人澡人人爽秒播| 成人国语在线视频| 久久精品aⅴ一区二区三区四区| 伦理电影免费视频| 久久香蕉国产精品| 久久精品亚洲精品国产色婷小说| 99久久无色码亚洲精品果冻| 亚洲 欧美 日韩 在线 免费| 日本一区二区免费在线视频| 无遮挡黄片免费观看| 色av中文字幕| 久久中文看片网| 免费在线观看日本一区| 我要搜黄色片| 欧美成人午夜精品| 在线观看日韩欧美| 村上凉子中文字幕在线| 伦理电影免费视频| 亚洲中文av在线| 搡老熟女国产l中国老女人| а√天堂www在线а√下载| 国产精品一区二区精品视频观看| 欧美人与性动交α欧美精品济南到| 亚洲精品一区av在线观看| 婷婷丁香在线五月| 中文字幕人妻丝袜一区二区| 欧美高清成人免费视频www| 黄频高清免费视频| 黑人巨大精品欧美一区二区mp4| 一级作爱视频免费观看| 无遮挡黄片免费观看| 国产av一区在线观看免费| 九色成人免费人妻av| 国产亚洲精品av在线| 亚洲av熟女| 亚洲人成网站在线播放欧美日韩| 午夜激情福利司机影院| 国产高清激情床上av| 免费在线观看日本一区| 亚洲av第一区精品v没综合| 中文亚洲av片在线观看爽| 成人av在线播放网站| 国产三级在线视频| 色尼玛亚洲综合影院| 波多野结衣高清作品| 欧美日韩国产亚洲二区| 精品久久久久久久末码| 午夜视频精品福利| 国产成人欧美在线观看| 香蕉av资源在线| 男男h啪啪无遮挡| 狠狠狠狠99中文字幕| 精品一区二区三区视频在线观看免费| e午夜精品久久久久久久| 日韩欧美一区二区三区在线观看| 精品高清国产在线一区| 久久久久久久精品吃奶| 一区二区三区高清视频在线| 狂野欧美激情性xxxx| 黄片小视频在线播放| 一区二区三区国产精品乱码| 大型av网站在线播放| 老汉色∧v一级毛片| 国产区一区二久久| 国产精品免费一区二区三区在线| 中文字幕av在线有码专区| av视频在线观看入口| 色哟哟哟哟哟哟| 亚洲免费av在线视频| 国产免费av片在线观看野外av| 麻豆一二三区av精品| 日韩av在线大香蕉| 每晚都被弄得嗷嗷叫到高潮| 成熟少妇高潮喷水视频| 黄色视频不卡| 国产伦在线观看视频一区| 一本久久中文字幕| 村上凉子中文字幕在线| 亚洲av成人精品一区久久| 黄色视频,在线免费观看| 精品不卡国产一区二区三区| 热99re8久久精品国产| 在线观看一区二区三区| 啦啦啦韩国在线观看视频| 午夜亚洲福利在线播放| 色哟哟哟哟哟哟| 91麻豆av在线| www.999成人在线观看| 亚洲va日本ⅴa欧美va伊人久久| a在线观看视频网站| 亚洲专区中文字幕在线| www国产在线视频色| 香蕉av资源在线| 麻豆一二三区av精品| 69av精品久久久久久| 日韩免费av在线播放| 亚洲精品国产一区二区精华液| 国产一区二区在线av高清观看| 亚洲午夜精品一区,二区,三区| 麻豆国产av国片精品| 久久人人精品亚洲av| 制服人妻中文乱码| 国产精品爽爽va在线观看网站| 国产精品99久久99久久久不卡| 亚洲欧美日韩高清专用| 久久国产乱子伦精品免费另类| 淫妇啪啪啪对白视频| 欧美成人性av电影在线观看| 亚洲熟女毛片儿| 在线观看免费午夜福利视频| 悠悠久久av| 国产欧美日韩一区二区精品| 中文亚洲av片在线观看爽| 亚洲av熟女| 天堂影院成人在线观看| 999久久久国产精品视频| 中亚洲国语对白在线视频| 啦啦啦观看免费观看视频高清| 91麻豆精品激情在线观看国产| 香蕉久久夜色| 成人手机av| 亚洲国产精品999在线| 精品久久久久久久久久免费视频| 日本免费a在线| 国产一区二区三区视频了| 婷婷丁香在线五月| 国产亚洲欧美在线一区二区| 婷婷精品国产亚洲av| 美女 人体艺术 gogo| 美女黄网站色视频| 搞女人的毛片| 在线观看日韩欧美| 丰满的人妻完整版| 91九色精品人成在线观看| 制服诱惑二区| 亚洲精品久久成人aⅴ小说| 一级毛片高清免费大全| 高清毛片免费观看视频网站| 丁香欧美五月| 男男h啪啪无遮挡| 男人的好看免费观看在线视频 | 高清在线国产一区| 国产精品精品国产色婷婷| 身体一侧抽搐| 一卡2卡三卡四卡精品乱码亚洲| 变态另类丝袜制服| 在线a可以看的网站| 少妇的丰满在线观看| 麻豆国产97在线/欧美 | 日本免费a在线| 欧美一区二区精品小视频在线| 国产高清videossex| 国内毛片毛片毛片毛片毛片| 日韩欧美 国产精品| 伦理电影免费视频| 性色av乱码一区二区三区2| 国产高清视频在线播放一区| 正在播放国产对白刺激| 我要搜黄色片| 999久久久国产精品视频| 久久国产精品人妻蜜桃| 亚洲欧美一区二区三区黑人| 丰满人妻熟妇乱又伦精品不卡| 免费看美女性在线毛片视频| 欧美一区二区精品小视频在线| 中文字幕最新亚洲高清| 亚洲国产欧美人成|