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

    基于能量流網(wǎng)絡(luò)仿真的鋼鐵工業(yè)多能源介質(zhì)優(yōu)化調(diào)配

    2017-03-10 05:20:36孫彥廣梁青艷李文兵賈天云
    自動(dòng)化學(xué)報(bào) 2017年6期
    關(guān)鍵詞:優(yōu)化模型

    孫彥廣 梁青艷 李文兵 賈天云

    鋼鐵工業(yè)是能源密集型產(chǎn)業(yè),我國(guó)鋼鐵工業(yè)能耗約占全國(guó)工業(yè)總能耗的16%[1].目前我國(guó)鋼鐵企業(yè)的噸鋼能耗仍比世界先進(jìn)水平國(guó)家高出10%左右[2].面對(duì)節(jié)能減排的嚴(yán)峻形勢(shì),鋼鐵工業(yè)亟需通過(guò)綠色化、智能化實(shí)現(xiàn)可持續(xù)發(fā)展.

    能源優(yōu)化調(diào)配是鋼鐵企業(yè)系統(tǒng)節(jié)能的關(guān)鍵技術(shù)之一.通過(guò)能源優(yōu)化調(diào)配可以實(shí)現(xiàn)能質(zhì)匹配,提高能源的利用率,降低二次能源放散,在滿足生產(chǎn)對(duì)能源質(zhì)量和數(shù)量要求前提下,降低能源成本.

    鋼鐵工業(yè)能源系統(tǒng)的特點(diǎn)對(duì)能源調(diào)配技術(shù)提出了挑戰(zhàn).首先,鋼鐵企業(yè)能源介質(zhì)種類繁多,包括煤、焦炭、煤氣、電力、蒸汽、技術(shù)氣體、壓縮空氣、水等近30種.其次,各種能源介質(zhì)與鋼鐵生產(chǎn)流程耦合緊密,很多二次能源介質(zhì)直接產(chǎn)生于鋼鐵生產(chǎn)過(guò)程的副產(chǎn)品或余熱余能回收利用,如高爐煤氣、焦?fàn)t煤氣、轉(zhuǎn)爐煤氣、各種蒸汽回收、CDQ(Coke dry quenching)發(fā)電和TRT(Blast furnace top gas recovery turbine unit)發(fā)電等.此外,各種能源介質(zhì)的產(chǎn)生、轉(zhuǎn)換、存儲(chǔ)、輸送和分配使用通過(guò)能源管網(wǎng)實(shí)現(xiàn),構(gòu)成了復(fù)雜的相互制約的能量流網(wǎng)絡(luò).

    對(duì)于鋼鐵工業(yè)單一介質(zhì)優(yōu)化調(diào)配問(wèn)題,很多學(xué)者進(jìn)行了研究.在鋼鐵副產(chǎn)煤氣調(diào)配方面,針對(duì)煤氣產(chǎn)生消耗量預(yù)測(cè)問(wèn)題,劉穎等[3]基于改進(jìn)回聲狀態(tài)網(wǎng)絡(luò)用于高爐煤氣發(fā)生量預(yù)測(cè),張顏顏等[4]采用改進(jìn)的數(shù)據(jù)驅(qū)動(dòng)子空間算法求解鋼鐵企業(yè)能源預(yù)測(cè)問(wèn)題,肖冬峰等[5]基于改進(jìn)BP(Back propagation)網(wǎng)絡(luò)的高爐煤氣發(fā)生量預(yù)測(cè)模型;針對(duì)煤氣優(yōu)化調(diào)配問(wèn)題,Akimoto等[6]采用多周期混合整數(shù)線性規(guī)劃模型研究煤氣柜位控制和自備電廠煤氣分配問(wèn)題,張琦等[7]提出了鋼鐵聯(lián)合企業(yè)煤氣系統(tǒng)優(yōu)化分配模型,Giacomo等[8?9]從煤氣利用效益最大和二氧化碳排放最小角度研究了煤氣系統(tǒng)多目標(biāo)優(yōu)化求解算法,Yang等[10]研究了考慮副產(chǎn)煤氣存儲(chǔ)能力的富裕煤氣最優(yōu)分配問(wèn)題.在氧氣優(yōu)化調(diào)配方面,Han等[11]研究了考慮電費(fèi)的氧氣系統(tǒng)優(yōu)化問(wèn)題.在電力系統(tǒng)調(diào)配方面,王小輝[12]研究了某鋼廠電力負(fù)荷模擬與預(yù)測(cè)研究問(wèn)題,Ashok[13]提出了一種適合工業(yè)熱電聯(lián)產(chǎn)系統(tǒng)的最優(yōu)能源優(yōu)化分配模型,高云龍等[14]研究了高耗能企業(yè)關(guān)口平衡優(yōu)化調(diào)度及其輸出功率控制方式,劉坤等[15]研究了考慮負(fù)荷及煤氣量不確定性的企業(yè)微電網(wǎng)自發(fā)電調(diào)度模型,張玉慶等[16]提出了綜合考慮發(fā)電、配電、用電環(huán)節(jié)的鋼鐵企業(yè)智能微電網(wǎng)系統(tǒng)運(yùn)行框架.鋼鐵工業(yè)單一介質(zhì)優(yōu)化調(diào)配,不能考慮各種能源介質(zhì)的產(chǎn)生、轉(zhuǎn)換間的關(guān)聯(lián)關(guān)系,難以取得多種能源介質(zhì)綜合優(yōu)化的效果.

    對(duì)于鋼鐵工業(yè)多介質(zhì)優(yōu)化調(diào)配問(wèn)題,近年來(lái)也有一些學(xué)者進(jìn)行了研究.李丹[17]提出了基于價(jià)值最大化的能源系統(tǒng)綜合調(diào)整策略,孟華[18]研究了鋼鐵企業(yè)自備電廠機(jī)組配置優(yōu)化及煤氣優(yōu)化調(diào)度問(wèn)題,張琦等[19]建立了鋼鐵企業(yè)副產(chǎn)煤氣–蒸汽–電力耦合模型,通過(guò)ILOG Cplex求解出模型的最優(yōu)解,獲得煤氣、蒸汽和電力的最優(yōu)分配方案,孫彥廣[20]基于能量流網(wǎng)絡(luò)信息模型提出了6大類多種能源介質(zhì)分解–協(xié)調(diào)動(dòng)態(tài)調(diào)控策略,曾玉嬌等[21?22]建立了鋼鐵企業(yè)蒸汽和發(fā)電綜合優(yōu)化調(diào)度模型.何佳毅等[23]研究了鋼鐵企業(yè)能源系統(tǒng)網(wǎng)絡(luò)模型仿真及組態(tài)問(wèn)題,Nawzad[24]提出了仿真和優(yōu)化結(jié)合的方法為提高能效提供決策支持.

    上述鋼鐵工業(yè)多介質(zhì)優(yōu)化調(diào)配方案,并沒(méi)有充分考慮鋼鐵工業(yè)能源系統(tǒng)與鋼材生產(chǎn)系統(tǒng)耦合緊密的特點(diǎn),事實(shí)上,生產(chǎn)系統(tǒng)品種、產(chǎn)量、設(shè)備狀態(tài)和工藝路徑的不同,能源系統(tǒng)各介質(zhì)的產(chǎn)生、轉(zhuǎn)化、分配和使用需求不同,都會(huì)導(dǎo)致各種能源介質(zhì)的平衡關(guān)系、優(yōu)化約束邊界條件發(fā)生變化,使優(yōu)化效果大打折扣.

    論文從鋼鐵制造流程物質(zhì)流能量流耦合特點(diǎn)出發(fā),提出了基于能量流網(wǎng)絡(luò)動(dòng)態(tài)仿真的鋼鐵工業(yè)多能源介質(zhì)綜合優(yōu)化調(diào)配策略.首先,基于主生產(chǎn)工序的能量流模型、分介質(zhì)能量流網(wǎng)絡(luò)模型建立鋼鐵企業(yè)物質(zhì)流能量流集成的能量流網(wǎng)絡(luò)模型,充分表征物質(zhì)流、能量流相互耦合和相互影響;然后,通過(guò)輸入當(dāng)前生產(chǎn)計(jì)劃、工藝路徑、設(shè)備運(yùn)行狀況等信息進(jìn)行鋼鐵制造全流程仿真,來(lái)識(shí)別、調(diào)整生產(chǎn)系統(tǒng)對(duì)能源系統(tǒng)的時(shí)變需求,形成多時(shí)間周期動(dòng)態(tài)變化的優(yōu)化約束邊界條件;最后,針對(duì)不同生產(chǎn)場(chǎng)景進(jìn)行優(yōu)化求解,給出對(duì)應(yīng)的多能源介質(zhì)動(dòng)態(tài)優(yōu)化調(diào)配.

    論文首先探討了能量流網(wǎng)絡(luò)化建模、生產(chǎn)流程與能源系統(tǒng)結(jié)合的能源仿真技術(shù);然后,探討了多能源介質(zhì)綜合優(yōu)化調(diào)配技術(shù);最后,通過(guò)在鋼鐵企業(yè)示范應(yīng)用案例驗(yàn)證了論文所提方案的有效性.

    1 鋼鐵工業(yè)能量流網(wǎng)絡(luò)化模型

    鋼鐵工業(yè)能量流貫穿于能源使用、能源回收和能源轉(zhuǎn)換輸配三個(gè)環(huán)節(jié),在不同環(huán)節(jié)有不同的能源介質(zhì)表現(xiàn)形式.完整描述鋼鐵企業(yè)能量流網(wǎng)絡(luò)的信息模型結(jié)構(gòu)包括主生產(chǎn)工序的能量流模型、分介質(zhì)能量流網(wǎng)絡(luò)模型和鋼鐵企業(yè)能量流網(wǎng)絡(luò)集成模型.

    其中,主生產(chǎn)工序的能量流模型描述各生產(chǎn)工序能源使用和回收情況,分介質(zhì)能量流網(wǎng)絡(luò)模型描述各種能源介質(zhì)產(chǎn)生、轉(zhuǎn)換、輸配情況.鋼鐵企業(yè)能量流網(wǎng)絡(luò)集成模型將能源使用、回收、轉(zhuǎn)換輸配三個(gè)環(huán)節(jié)能源信息關(guān)聯(lián)起來(lái),并形成多種能源介質(zhì)間的調(diào)控[20].圖1為三個(gè)模型之間的信息流和控制流示意圖.

    1.1 主生產(chǎn)工序的能量流模型

    國(guó)內(nèi)外對(duì)鋼鐵企業(yè)主生產(chǎn)工序的能量流模型做了很多研究[25?30],從生產(chǎn)工序物料平衡和能源平衡的角度,分析了各種工藝制度對(duì)鋼鐵企業(yè)能源結(jié)構(gòu)的影響,以及鋼鐵企業(yè)余熱余能循環(huán)利用的潛力.這些模型的能耗計(jì)算是按鋼比進(jìn)行靜態(tài)匹配計(jì)算,不能反映生產(chǎn)節(jié)奏、設(shè)備檢修故障、能源緩沖等動(dòng)態(tài)變化,而且這些模型沒(méi)有考慮與能源網(wǎng)絡(luò)的信息交互,無(wú)法滿足能源動(dòng)態(tài)調(diào)配的需要.

    針對(duì)上述問(wèn)題,采取了以下措施,建立主生產(chǎn)工序的工序能量流模型[20]:

    1)將主生產(chǎn)工序設(shè)備細(xì)分為能源消耗設(shè)備和能源回收設(shè)備,通過(guò)啟動(dòng)/停止等狀態(tài)進(jìn)行標(biāo)識(shí),以反映設(shè)備正常、故障等動(dòng)態(tài)運(yùn)行狀況,設(shè)備狀態(tài)信號(hào)來(lái)源于設(shè)備檢修計(jì)劃和設(shè)備管理系統(tǒng);煤氣管網(wǎng)、蒸汽管網(wǎng)和自發(fā)電機(jī)組等能源調(diào)配響應(yīng)周期為20分鐘到40分鐘,模型循環(huán)計(jì)算的時(shí)間粒度確定為30分鐘.

    2)在能量流模型中引入鋼鐵生產(chǎn)流程動(dòng)態(tài)調(diào)控機(jī)制:根據(jù)生產(chǎn)計(jì)劃,確定工藝路徑、主生產(chǎn)工序單元種類和每單元鋼材產(chǎn)量;根據(jù)作業(yè)計(jì)劃、維護(hù)計(jì)劃決定各單元的運(yùn)行時(shí)間;各生產(chǎn)單元按工藝路徑依次計(jì)算,并考慮前一工序輸出和庫(kù)存對(duì)后工序運(yùn)行的影響.

    3)單元模型的輸入/輸出與各能源介質(zhì)網(wǎng)絡(luò)模型對(duì)應(yīng),便于能源調(diào)控方案的研究.

    主生產(chǎn)工序能源輸入輸出模型由于其與生產(chǎn)品種、產(chǎn)量、工況、原料結(jié)構(gòu)、工藝制度等密切相關(guān),很難用一個(gè)統(tǒng)計(jì)模型描述.論文研究中分析了影響介質(zhì)波動(dòng)的因素,把影響因素分為靜態(tài)因素、動(dòng)態(tài)因素及本身波動(dòng)特性,采用基于工況信息的分段建模方法進(jìn)行模型描述.

    圖1 能量流網(wǎng)絡(luò)模型的信息流與控制Fig.1 Information flow and control of energy flow network model

    以高爐煤氣回收量模型為例.正常生產(chǎn)、休風(fēng)、減風(fēng)等不同工況下煤氣流量模型描述如下:

    1.2 分介質(zhì)能量流網(wǎng)絡(luò)模型

    分介質(zhì)能量流網(wǎng)絡(luò)模型是從燃?xì)?、蒸汽、電、技術(shù)氣體(氧氬氮)、壓縮空氣和水等分介質(zhì)管網(wǎng)角度,將主生產(chǎn)工序分能源介質(zhì)的消耗、回收與能源系統(tǒng)分能源介質(zhì)的產(chǎn)生、儲(chǔ)存、轉(zhuǎn)換分配連接起來(lái),建立分介質(zhì)管網(wǎng)數(shù)學(xué)模型.

    能源管網(wǎng)拓?fù)浣Y(jié)構(gòu)的描述利用圖論的有向圖原理,通過(guò)關(guān)聯(lián)矩陣(樹(shù)枝矩陣、連枝矩陣)和基本回路矩陣將管網(wǎng)圖形(枝狀網(wǎng)與環(huán)狀網(wǎng)混合)信息數(shù)據(jù)化,并與能源節(jié)點(diǎn)相關(guān)聯(lián).能源管網(wǎng)基本方程.根據(jù)流體網(wǎng)絡(luò)的一些基本定律,如質(zhì)量守恒定律、能量守恒定律、阻力定律,確定連續(xù)性方程、能量方程和壓降方程,并對(duì)管段摩阻系數(shù)進(jìn)行辨識(shí).

    1.3 鋼鐵企業(yè)能量流網(wǎng)絡(luò)集成模型

    在主生產(chǎn)工序的能量流模型和分介質(zhì)能量流網(wǎng)絡(luò)模型基礎(chǔ)上,建立鋼鐵企業(yè)能量流網(wǎng)絡(luò)集成模型,將主生產(chǎn)工序的能量流模型的計(jì)算結(jié)果與分介質(zhì)能量流網(wǎng)絡(luò)模型計(jì)算結(jié)果銜接起來(lái),如圖2所示.

    主生產(chǎn)工序作為介質(zhì)能量流網(wǎng)絡(luò)的終端節(jié)點(diǎn)從能源管網(wǎng)獲取所需要的各種能源介質(zhì),產(chǎn)生主產(chǎn)品、副產(chǎn)品,同時(shí)回收能源介質(zhì),回收的能源介質(zhì)又作為介質(zhì)能量流網(wǎng)絡(luò)的始端節(jié)點(diǎn).介質(zhì)能量流網(wǎng)絡(luò)模型始端節(jié)點(diǎn)除了回收能源外包括自產(chǎn)和外購(gòu)能源,中間環(huán)節(jié)包括存儲(chǔ)、轉(zhuǎn)換分配及連接介質(zhì)管網(wǎng),終端節(jié)點(diǎn)包括主產(chǎn)用戶、能源生產(chǎn)單元用戶、公輔及小用戶,多余的能源放散或外送.其中能源生產(chǎn)單元在不同的介質(zhì)能量流網(wǎng)絡(luò)中角色不同,一方面在本身的產(chǎn)出介質(zhì)能量流網(wǎng)絡(luò)中作為始端節(jié)點(diǎn),產(chǎn)生自產(chǎn)能源;另一方面因?yàn)槟茉磫卧蚕哪茉?因此在其他介質(zhì)的能量流網(wǎng)絡(luò)中,作為終端節(jié)點(diǎn),消耗能源.

    鋼鐵企業(yè)能量流網(wǎng)絡(luò)集成模型包括兩部分:1)主生產(chǎn)工序與介質(zhì)能量流網(wǎng)絡(luò)模型集成;2)各介質(zhì)能量流網(wǎng)絡(luò)之間集成.圖3為主生產(chǎn)工序之一的高爐與介質(zhì)能量流網(wǎng)絡(luò)模型集成示意圖,高爐從各介質(zhì)能量流網(wǎng)絡(luò)模型獲取能源作為能源輸入,產(chǎn)生鐵水,同時(shí)回收高爐煤氣,產(chǎn)生的高爐煤氣又作為高爐煤氣能量流網(wǎng)絡(luò)模型的輸入.高爐消耗及回收能源模型在主生產(chǎn)工序能量流網(wǎng)絡(luò)模型中描述.圖4為某鋼鐵企業(yè)電力系統(tǒng)與其他介質(zhì)能量流網(wǎng)絡(luò)集成模型示意圖,電力能量流網(wǎng)絡(luò)始端輸入節(jié)點(diǎn)包括300MW煤氣和燃煤CHP(Combined heat and power)發(fā)電、25MW余熱發(fā)電、CDQ余熱發(fā)電、TRT余能發(fā)電及外購(gòu)能源,終端節(jié)點(diǎn)除給各主工序生產(chǎn)用戶外,還作為其他能量流網(wǎng)絡(luò)(技術(shù)氣體、壓縮空氣、水等)始端節(jié)點(diǎn)的能量輸入,各介質(zhì)能量流網(wǎng)絡(luò)集成模型主要考慮能源介質(zhì)之間的轉(zhuǎn)換.

    圖2 鋼鐵企業(yè)能量流網(wǎng)絡(luò)集成模型結(jié)構(gòu)框圖Fig.2 Structural block diagram of energy flow network integration model for iron and steel enterprises

    圖3 高爐工序與能量流網(wǎng)絡(luò)模型集成Fig.3 Blast furnace process and energy flow network model integration

    2 生產(chǎn)流程與能源系統(tǒng)結(jié)合的能源仿真技術(shù)

    基于上述能量流網(wǎng)絡(luò)化模型,開(kāi)發(fā)了生產(chǎn)流程與能源系統(tǒng)結(jié)合的能源仿真計(jì)算軟件,如圖5所示.

    圖4 電力能量流網(wǎng)絡(luò)與其他介質(zhì)能量流網(wǎng)絡(luò)集成Fig.4 Power energy flow network and other media energy flow network integration

    圖5 基于能量流網(wǎng)絡(luò)模型的能源系統(tǒng)動(dòng)態(tài)仿真Fig.5 Dynamic simulation of energy system based on energy flow network model

    圖6 介質(zhì)之間優(yōu)先級(jí)依賴關(guān)系Fig.6 Priority dependencies between media

    能源系統(tǒng)動(dòng)態(tài)仿真首先確定介質(zhì)之間優(yōu)先級(jí)依賴關(guān)系,如圖6所示.鋼鐵企業(yè)能源介質(zhì)種類繁多,大致可分為燃?xì)夂兔骸⒄羝?、電、技術(shù)氣體、壓縮空氣和水6大類,分析能源介質(zhì)間消耗–轉(zhuǎn)換–產(chǎn)生的“鏈條”關(guān)系,燃?xì)庀到y(tǒng)處于轉(zhuǎn)換“鏈條”上游,電、蒸汽系統(tǒng)處于轉(zhuǎn)換“鏈條”中游,技術(shù)氣體、壓縮空氣、水處于轉(zhuǎn)換“鏈條”下游.據(jù)此,把技術(shù)氣體、風(fēng)、水的需量轉(zhuǎn)化為對(duì)電、蒸汽、燃?xì)獾男枨?重點(diǎn)解決煤氣–蒸汽–電三類能源介質(zhì)的調(diào)配.

    對(duì)于每一種介質(zhì),基于一定的調(diào)配規(guī)則按次序?qū)Ω髫?fù)荷進(jìn)行分配,圖7給出副產(chǎn)煤氣的分配次序.副產(chǎn)煤氣的用戶分為基準(zhǔn)用戶和緩沖用戶.基準(zhǔn)用戶為主生產(chǎn)工序用戶,如高爐熱風(fēng)爐、軋鋼加熱爐等,緩沖用戶為能源系統(tǒng)用戶,如蒸汽鍋爐和自發(fā)電機(jī)組.基準(zhǔn)用戶的優(yōu)先級(jí)高于緩沖用戶.因此仿真計(jì)算時(shí),首先根據(jù)全流程能源需求和可回收計(jì)算結(jié)果,滿足基準(zhǔn)用戶需求,多余煤氣為可供緩沖用戶量.根據(jù)蒸汽和電需求,計(jì)算出緩沖用戶實(shí)際使用量,可供緩沖用戶量與緩沖用戶實(shí)際使用量之差為煤氣多余量.

    動(dòng)態(tài)仿真過(guò)程是事件驅(qū)動(dòng)的.仿真計(jì)算過(guò)程中,在主生產(chǎn)工序能量模型引入了生產(chǎn)設(shè)備啟/停和正常/異常生產(chǎn)工況的事件標(biāo)識(shí),反映生產(chǎn)過(guò)程的動(dòng)態(tài)變化;在分介質(zhì)能量流模型中,設(shè)置了能源設(shè)備啟停標(biāo)識(shí),同時(shí)能源介質(zhì)管網(wǎng)平衡模型考慮了管網(wǎng)壓力、氣柜柜位的計(jì)算,反映能源系統(tǒng)設(shè)備運(yùn)行、能源放散等動(dòng)態(tài)變化.

    通過(guò)鋼鐵企業(yè)能源系統(tǒng)模擬仿真計(jì)算,可以分析各能源介質(zhì)在能源使用、能源回收和轉(zhuǎn)換輸配三個(gè)環(huán)節(jié)動(dòng)態(tài)變化情況,分析評(píng)估各環(huán)節(jié)的效率和流程綜合效率;分析對(duì)比各種鋼鐵制造流程、生產(chǎn)作業(yè)計(jì)劃、中間緩沖能力和二次能源轉(zhuǎn)換輸配方案對(duì)能耗和效率的影響,為鋼鐵企業(yè)能源系統(tǒng)的設(shè)計(jì)方案和運(yùn)行策略的對(duì)比和優(yōu)化提供定量分析手段.圖8為某鋼鐵企業(yè)一段時(shí)間內(nèi)煤氣和電力供需平衡情況.

    3 基于能源仿真的多能源介質(zhì)綜合優(yōu)化調(diào)配

    3.1 基于能源仿真的多能源介質(zhì)綜合優(yōu)化調(diào)配流程

    以上基于調(diào)度規(guī)則的仿真運(yùn)行模式,可以支持以滿足生產(chǎn)用戶需求和放散最小為目標(biāo)的多介質(zhì)能源調(diào)配.在此基礎(chǔ)上,引入優(yōu)化算法實(shí)現(xiàn)以多介質(zhì)能源成本最小為目標(biāo)的多能源介質(zhì)綜合優(yōu)化調(diào)配.

    鋼鐵企業(yè)生產(chǎn)品種、產(chǎn)量、設(shè)備狀態(tài)和工藝路徑的不同,能源系統(tǒng)各介質(zhì)的產(chǎn)生、轉(zhuǎn)化、分配和使用方式不同,都會(huì)導(dǎo)致各種能源介質(zhì)的平衡關(guān)系、優(yōu)化約束邊界條件發(fā)生變化,使靜態(tài)優(yōu)化效果大打折扣.能源仿真的主要目的是輸入當(dāng)前生產(chǎn)計(jì)劃、工藝路徑、設(shè)備運(yùn)行狀況等信息,通過(guò)仿真計(jì)算實(shí)時(shí)識(shí)別、調(diào)整生產(chǎn)工序能量動(dòng)態(tài)平衡約束條件、能源介質(zhì)轉(zhuǎn)換動(dòng)態(tài)平衡約束條件和能源設(shè)備產(chǎn)能約束條件,提高多能源介質(zhì)綜合優(yōu)化調(diào)配對(duì)生產(chǎn)系統(tǒng)、能源系統(tǒng)動(dòng)態(tài)變化的適應(yīng)性.基于能源仿真的多能源介質(zhì)綜合優(yōu)化調(diào)配流程如圖9所示.

    圖7 煤氣的調(diào)配規(guī)則及各類負(fù)荷分配次序Fig.7 Gas distribution rules and various types of load distribution order

    圖8 某鋼鐵企業(yè)一段時(shí)間內(nèi)煤氣和電力供需平衡情況Fig.8 Gas and electricity supply and demand balance in a period of time in a steel enterprise

    圖9 基于能源仿真的多能源介質(zhì)綜合優(yōu)化調(diào)配流程Fig.9 Comprehensive optimization of multi-energy medium based on energy simulation

    1)基于能量流網(wǎng)絡(luò)模型的仿真.首先,根據(jù)生產(chǎn)計(jì)劃確定流程工藝路徑及工序單元,生成作業(yè)計(jì)劃和維修計(jì)劃,調(diào)用主生產(chǎn)工序能量流模型計(jì)算生產(chǎn)工序各種能源介質(zhì)消耗量和回收量,以及生產(chǎn)工序各能源介質(zhì)凈需求,據(jù)此形成主生產(chǎn)工序能量動(dòng)態(tài)平衡約束條件;然后,調(diào)用分介質(zhì)能量流模型進(jìn)行能源系統(tǒng)介質(zhì)轉(zhuǎn)化計(jì)算,將技術(shù)氣體、壓縮空氣、水、氫氣等介質(zhì)需求轉(zhuǎn)換成對(duì)于煤氣、蒸汽和電力需求,形成能源介質(zhì)轉(zhuǎn)換動(dòng)態(tài)平衡約束條件;同時(shí)調(diào)用分介質(zhì)能量流網(wǎng)絡(luò)模型,根據(jù)能源設(shè)備啟/停標(biāo)識(shí)和能源介質(zhì)管網(wǎng)平衡模型生產(chǎn)能源設(shè)備產(chǎn)能動(dòng)態(tài)約束條件.

    2)多能源介質(zhì)動(dòng)態(tài)調(diào)配優(yōu)化.匯總上述仿真信息,生產(chǎn)目標(biāo)函數(shù)和約束條件,進(jìn)行煤氣–蒸汽–電力動(dòng)態(tài)調(diào)配優(yōu)化,進(jìn)而得到技術(shù)氣體、壓縮空氣、水、氫氣、煤氣、蒸汽和電力分介質(zhì)調(diào)配方案.

    3.2 目標(biāo)函數(shù)

    鋼鐵企業(yè)需要充分利用一次能源(如動(dòng)力煤)、副產(chǎn)二次能源(副產(chǎn)煤氣)和余熱余能資源(余熱、余壓),安排各生產(chǎn)設(shè)備進(jìn)行蒸汽和電力的優(yōu)化生產(chǎn),在滿足生產(chǎn)過(guò)程對(duì)各種能源介質(zhì)的需求下,降低整個(gè)系統(tǒng)的運(yùn)行費(fèi)用,減少煤氣放散,降低能源消耗,提高企業(yè)的經(jīng)濟(jì)效益.為此,目標(biāo)函數(shù)選為給定調(diào)度周期內(nèi)整個(gè)鋼鐵企業(yè)煤氣–蒸汽–電力系統(tǒng)運(yùn)行費(fèi)用最低,包括燃料費(fèi)用、鍋爐給水費(fèi)用、設(shè)備維護(hù)費(fèi)用、煤氣放散懲罰費(fèi)用,電能供給不足時(shí)的外購(gòu)電費(fèi)以及電能富余時(shí)的外送電收益[21?22,31].技術(shù)氣體、壓縮空氣和水等能源介質(zhì)的費(fèi)用通過(guò)煤氣、蒸汽、電力三類能源介質(zhì)費(fèi)用來(lái)表征.

    式中,T為一個(gè)調(diào)度周期內(nèi)所包含的時(shí)段數(shù)目,Nse表示蒸汽生產(chǎn)設(shè)備的數(shù)目,Npe表示電力生產(chǎn)設(shè)備的數(shù)目,Ng表示副產(chǎn)煤氣的數(shù)目,τ表示每個(gè)操作時(shí)段的工作時(shí)間,Cg表示副產(chǎn)煤氣g的價(jià)格,表示設(shè)備i在t時(shí)段副產(chǎn)煤氣g消耗量,Ccoal表示外購(gòu)煤的價(jià)格,表示設(shè)備i在t時(shí)段燃料煤的消耗量,CM,i和Mi,t分別表示產(chǎn)汽設(shè)備的制造費(fèi)用(包括設(shè)備折舊、維修費(fèi)用,人工薪酬等)和在t時(shí)段的蒸汽產(chǎn)量,CM,j和Pi,t分別表示發(fā)電設(shè)備的設(shè)備維護(hù)成本(元)以及在t時(shí)段的發(fā)電量,Crel,g表示副產(chǎn)煤氣g的放散懲罰價(jià)格,Rg,t表示副產(chǎn)煤氣g在t時(shí)段的放散量,Pw,t表示在t時(shí)段的關(guān)口交換功率,δt表示外網(wǎng)在t時(shí)段的供電狀態(tài),,Cb,t為t時(shí)段的外購(gòu)電價(jià),Cs,t為時(shí)段t的外送電價(jià).

    3.3 約束條件

    包括生產(chǎn)工序能量動(dòng)態(tài)平衡約束條件、能源介質(zhì)轉(zhuǎn)換動(dòng)態(tài)平衡約束條件和能源設(shè)備產(chǎn)能約束條件,根據(jù)從鋼鐵企業(yè)制造執(zhí)行系統(tǒng)(Manufacturing execution systems,MES)、能源管理系統(tǒng)(Energy management systems,EMS)獲取的當(dāng)前生產(chǎn)計(jì)劃、工藝路徑、設(shè)備運(yùn)行狀況等信息,通過(guò)仿真計(jì)算進(jìn)行動(dòng)態(tài)調(diào)整.

    1)蒸汽、電力、煤氣需求平衡約束

    鋼鐵企業(yè)在不同時(shí)段對(duì)蒸汽、電力、煤氣的需求是不同的,且在每個(gè)時(shí)段系統(tǒng)必須滿足能源供需平衡.

    2)蒸汽鍋爐能量平衡模型

    3)余熱余能回收發(fā)電設(shè)備(如CDQ)能量平衡模型

    4)利用副產(chǎn)煤氣的熱電聯(lián)產(chǎn)(CHP)設(shè)備能量平衡模型

    式中,ηk表示第k臺(tái)熱電聯(lián)產(chǎn)設(shè)備的效率,分別表示自發(fā)電設(shè)備k在t時(shí)段的r等級(jí)蒸汽的產(chǎn)量和發(fā)電功率.

    5)設(shè)備能力約束

    蒸汽鍋爐產(chǎn)汽量能力

    汽輪機(jī)進(jìn)汽量約束

    汽輪機(jī)抽汽能力

    發(fā)電設(shè)備有功出力約束

    6)發(fā)電機(jī)組燃料約束

    發(fā)電機(jī)組燃?xì)忮仩t對(duì)煤氣量使用的上下限范圍

    對(duì)于多種煤氣混燒鍋爐,使用的煤氣不但對(duì)量有要求,而且對(duì)其質(zhì)也有要求,通常要求多種煤氣的混合熱值要大于最低要求.

    式中,βi第i臺(tái)發(fā)電機(jī)對(duì)混合煤氣熱值的最低要求(kJ/kg).

    7)富裕煤氣供應(yīng)約束

    在整個(gè)調(diào)度周期內(nèi)的每一時(shí)間段內(nèi),副產(chǎn)煤氣的使用總量要小于或等于該時(shí)段煤氣供應(yīng)上限.

    8)發(fā)電機(jī)爬坡約束

    式中,URi為發(fā)電機(jī)i瞬間上調(diào)的約束值,即發(fā)電機(jī)在一個(gè)時(shí)段內(nèi)能增加最大功率;DRi為發(fā)電機(jī)i瞬間下調(diào)的約束值,即發(fā)電機(jī)在一個(gè)時(shí)段內(nèi)能減小的最大功率.

    3.4 優(yōu)化求解算法

    上述所建的多周期優(yōu)化模型屬于一類高維、非線性、多約束優(yōu)化問(wèn)題.非線性問(wèn)題由于混合煤氣熱值約束的非線性、鍋爐效率與蒸汽負(fù)荷的非線性、發(fā)電機(jī)組燃料消耗與發(fā)電量和抽汽量的二次曲線關(guān)系等引起.針對(duì)這類問(wèn)題,結(jié)合本文研究的多能源介質(zhì)優(yōu)化調(diào)配的特點(diǎn),采用粒子群算法(Partide swarm optimization,PSO)為基本算法.針對(duì)常規(guī)PSO存在的缺乏有效的約束處理機(jī)制、隨著問(wèn)題維數(shù)的增加容易陷入局部極值和收斂速度變慢等缺點(diǎn),對(duì)常規(guī)PSO算法進(jìn)行了改進(jìn),提出了混沌自適應(yīng)粒子群算法進(jìn)行優(yōu)化求解[21?22,31],該求解算法在文獻(xiàn)[21]中有詳細(xì)介紹.該算法具有以下特點(diǎn):

    1)在參數(shù)選擇上,采用有效動(dòng)態(tài)變化機(jī)制進(jìn)行自適應(yīng)調(diào)整.為了增加種群的多樣性,跳出局部極值,引入了自適應(yīng)變異操作.自適應(yīng)變異機(jī)制采用兩種不同變異算子即差分變異和高斯變異,對(duì)種群中的每個(gè)個(gè)體進(jìn)行變異,增加了種群的多樣性,引出算法跳出局部最優(yōu),并采用貪婪原則選擇適應(yīng)度最優(yōu)個(gè)體作為下一代,提高了算法全局尋優(yōu)性能.

    2)在約束條件的處理上,采用啟發(fā)式規(guī)則而不是常規(guī)的罰函數(shù)方法進(jìn)行約束處理.對(duì)于含有大量約束條件的優(yōu)化問(wèn)題,罰函數(shù)方法不僅帶來(lái)大量懲罰調(diào)整參數(shù)工作量,而且還導(dǎo)致尋優(yōu)效率降低甚至搜索不到最優(yōu)解.本文針對(duì)所建模型的特點(diǎn),采用一種有效的啟發(fā)式約束處理方法,通過(guò)對(duì)不可行個(gè)體進(jìn)行有效的逐次調(diào)整,使所有約束條件逐一得到滿足.

    4 應(yīng)用實(shí)例分析

    4.1 應(yīng)用背景[21]

    某鋼鐵企業(yè)蒸汽電力系統(tǒng)主要由4臺(tái)燃料鍋爐 (B1~B4)、2臺(tái)抽汽背壓式汽輪機(jī) (BT1和 BT2)、2套干熄焦發(fā)電裝置 (CDQ),2臺(tái)熱電聯(lián)產(chǎn)機(jī)組(CHP)組成,為生產(chǎn)工藝過(guò)程提供電力、S1蒸汽 (2.0~3.8MPa,450?C)、S2蒸汽 (0.78~1.27MPa,170~280?C)和 S3蒸汽(0.3~0.78MPa,142~170?C).其系統(tǒng)關(guān)系示意圖如圖10所示.

    其中,4臺(tái)鍋爐(B1~B4)消耗高爐煤氣(Blast furnace gas,BFG)和焦?fàn)t煤氣(Coke oven gas,COG)供S1和S2兩種品質(zhì)蒸汽;鍋爐B3和B4消耗動(dòng)力煤和副產(chǎn)煤氣為汽輪機(jī)(BT1和BT2)提供蒸汽S1發(fā)電,并產(chǎn)生S2抽汽和乏汽;2套CDQ設(shè)備回收紅焦的顯熱發(fā)電,并可以供S1和S2兩種品質(zhì)蒸汽.

    結(jié)合該企業(yè)煤氣–蒸汽–電力系統(tǒng)實(shí)際情況,建立該鋼鐵企業(yè)煤氣–蒸汽–電力集成優(yōu)化調(diào)度模型的目標(biāo)函數(shù)為

    圖10 某鋼鐵企業(yè)煤氣–蒸汽–電力系統(tǒng)示意圖Fig.10 A diagram of a gas-steam-electric power system in a steel enterprise

    各設(shè)備正常運(yùn)行條件約束方程如下:

    該優(yōu)化模型共包括24個(gè)調(diào)度時(shí)段,每個(gè)時(shí)段為1小時(shí).

    4.2 優(yōu)化結(jié)果分析

    采用上述混沌自適應(yīng)粒子群算法進(jìn)行優(yōu)化求解.通過(guò)20次獨(dú)立優(yōu)化計(jì)算可知,單次優(yōu)化計(jì)算時(shí)間小于3分鐘,并且所有的約束條件得到滿足,可以用于解決實(shí)際工程問(wèn)題.

    1)正常工況下優(yōu)化結(jié)果

    表1給出了通過(guò)動(dòng)態(tài)仿真計(jì)算得到的各時(shí)段的富余煤氣供應(yīng)量和蒸汽、電力需求量預(yù)測(cè)結(jié)果.

    圖11給出優(yōu)化前后整個(gè)鋼鐵企業(yè)煤氣–蒸汽–電力系統(tǒng)運(yùn)行費(fèi)用在各個(gè)時(shí)段的變化情況.由圖11分析可知,優(yōu)化后系統(tǒng)總運(yùn)行費(fèi)用在谷段下降幅度較大.

    表2給出了優(yōu)化前后各項(xiàng)費(fèi)用的比較結(jié)果,優(yōu)化后全天系統(tǒng)運(yùn)行總費(fèi)用比實(shí)際費(fèi)用降低了4%,優(yōu)化效果非常明顯.優(yōu)化前的調(diào)配方案是基于煤氣平衡關(guān)系和柜位變化,為避免煤氣放散,有多少富裕煤氣發(fā)多少電,沒(méi)有進(jìn)行煤氣–蒸汽–電力綜合優(yōu)化.優(yōu)化后的調(diào)配方案,在充分利用富余煤氣和蒸汽余熱資源的同時(shí),能有效利用峰平谷時(shí)段電價(jià)不同(見(jiàn)表3),優(yōu)化發(fā)電機(jī)組的出力,減少總的費(fèi)用,進(jìn)而有效地節(jié)約了能源,提高了鋼鐵企業(yè)經(jīng)濟(jì)效益.

    2)異常工況下優(yōu)化結(jié)果

    基于能源仿真的多能源介質(zhì)綜合優(yōu)化調(diào)配不僅可以用于正常工況下能源優(yōu)化,而且可以在異常工況或設(shè)備故障情況下通過(guò)動(dòng)態(tài)改變優(yōu)化的約束條件實(shí)現(xiàn)多能源介質(zhì)綜合優(yōu)化調(diào)配.

    圖11 各周期系統(tǒng)運(yùn)行費(fèi)用優(yōu)化前后對(duì)比Fig.11 Comparison of system operating costs of each cycle before and after optimization

    表4、表5分別為某鋼鐵公司煉鐵煉鋼減產(chǎn)10小時(shí)、2250軋線停產(chǎn)12小時(shí)兩種異常工況下能源優(yōu)化前后各種費(fèi)用對(duì)比.

    表1 各時(shí)段富余煤氣供應(yīng)及蒸汽和電力需求Table 1 Wealthy gas supply and steam and electricity demand for each period

    表2 正常工況優(yōu)化前后各種費(fèi)用比較Table 2 Comparison of various costs in normal working conditions before and after optimization

    表3 電網(wǎng)電價(jià)Table 3 Grid price

    表4 鐵鋼系統(tǒng)減產(chǎn)10小時(shí)工況優(yōu)化前后各種費(fèi)用對(duì)比Table 4 The comparison of various costs for 10 hour reduction of the steel system before and after optimization

    表5 2250軋線停產(chǎn)12小時(shí)工況優(yōu)化前后各種費(fèi)用對(duì)比Table 5 The comparison of various costs for 12 hours cut offof 2250 rolling line before and after the optimization

    5 結(jié)論

    鋼鐵工業(yè)多介質(zhì)優(yōu)化調(diào)配,需要考慮能源系統(tǒng)與鋼材生產(chǎn)系統(tǒng)耦合緊密、多能源介質(zhì)轉(zhuǎn)換的特點(diǎn),基于主生產(chǎn)工序的能量流模型、分介質(zhì)能量流網(wǎng)絡(luò)模型建立鋼鐵企業(yè)物質(zhì)流能量流集成的能量流網(wǎng)絡(luò)模型,是充分反映物質(zhì)流能量流相互耦合和能源介質(zhì)相互影響的有效表達(dá)方式.

    針對(duì)設(shè)備運(yùn)行狀態(tài)變化和生產(chǎn)工況變化,需要提高多能源介質(zhì)優(yōu)化調(diào)配的適應(yīng)能力,基于能量流網(wǎng)絡(luò)動(dòng)態(tài)仿真的多能源介質(zhì)優(yōu)化調(diào)配是一種可行的方案.通過(guò)與鋼鐵企業(yè)目前的信息化系統(tǒng)結(jié)合,及時(shí)獲取生產(chǎn)計(jì)劃、工藝路徑、設(shè)備運(yùn)行狀況等信息,通過(guò)能量流網(wǎng)絡(luò)動(dòng)態(tài)仿真確定不同工況對(duì)能源系統(tǒng)的需求變化,形成動(dòng)態(tài)的約束邊界條件,然后進(jìn)行優(yōu)化求解,得到對(duì)應(yīng)的多能源介質(zhì)動(dòng)態(tài)優(yōu)化調(diào)配方案.

    在某鋼鐵企業(yè)示范應(yīng)用案例表明論文所提方案對(duì)不同工況適應(yīng)性強(qiáng),能取得節(jié)能降成本效果.

    1 Wang Wei-Xing.Steel energy consumption and energy saving margin analysis.China Steel,2011,(4):19?22(王維興.鋼鐵工業(yè)能耗現(xiàn)狀和節(jié)能潛力分析.中國(guó)鋼鐵業(yè),2011,(4):19?22)

    2 Wang Ling,Jiang Fei-Tao.The current situation and prospect of energy saving and emission reduction in China0s steel industry.Industrial Economic Review,2012,(5):81?91(王嶺,江飛濤.中國(guó)鋼鐵工業(yè)節(jié)能減排效果分析與前景.產(chǎn)經(jīng)評(píng)論,2012,(5):81?91)

    3 Liu Ying,Zhao Jun,Wang Wei,Wu Yi-Ping,Chen Wei-Chang.Improved echo state network based on data-driven and its application to prediction of blast furnace gas output.Acta Automatica Sinica,2009,35(6):731?738(劉穎,趙珺,王偉,吳毅平,陳偉昌.基于數(shù)據(jù)的改進(jìn)回聲狀態(tài)網(wǎng)絡(luò)在高爐煤氣發(fā)生量預(yù)測(cè)中的應(yīng)用.自動(dòng)化學(xué)報(bào),2009,35(6):731?738)

    4 Zhang Yan-Yan,Tang Li-Xin.Improved data-driven subspace algorithm for energy prediction in iron and steel industry.Control Theory&Applications,2012,29(12):1616?1622(張顏顏,唐立新.改進(jìn)的數(shù)據(jù)驅(qū)動(dòng)子空間算法求解鋼鐵企業(yè)能源預(yù)測(cè)問(wèn)題.控制理論與應(yīng)用,2012,29(12):1616?1622)

    5 Xiao Dong-Feng,Yang Chun-Jie,Song Zhi-Huan.The forecasting model of blast furnace gas output based on improved BP network.Journal of Zhejiang University(Engineering Science),2012,46(11):2013?2018(肖冬峰,楊春節(jié),宋執(zhí)環(huán).基于改進(jìn)BP網(wǎng)絡(luò)的高爐煤氣發(fā)生量預(yù)測(cè)模型.浙江大學(xué)學(xué)報(bào)(工學(xué)版),2012,46(11):2013?2018)

    6 Akimoto K,Sannomiya N,Nishikawa Y,Tsuda T.An optimal gas supply for a power plant using a mixed integer programming model.Automatica,1991,27(3):513?518

    7 Zhang Qi,Cai Jiu-Ju,Pang Xing-Lu,Jiang Wen-Hao.Optimal distribution of by-product gases in iron and steel complex.Journal of Northeastern University(Natural Science),2011,32(1):98?101(張琦,蔡九菊,龐興露,姜文豪.鋼鐵聯(lián)合企業(yè)煤氣系統(tǒng)優(yōu)化分配模型.東北大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,32(1):98?101)

    8 Porzio G F,Fornai B,Amato A,Matarese N,Vannucci M,Chiappelli L,Colla V.Reducing the energy consumption and CO2emissions of energy intensive industries through decision support systems–an example of application to the steel industry.Applied Energy,2013,112:818?833

    9 Porzio G F,Nastasi G,Colla V,Vannucci M,Branca T A.Comparison of multi-objective optimization techniques applied to off-gas management within an integrated steelwork.Applied Energy,2014,136:1085?1097

    10 Yang J H,Cai J J,Sun W Q,Huang J.Optimal allocation of surplus gas and suitable capacity for buffer users in steel plant.Applied Thermal Engineering,2017,115:586?596

    11 Han Z Y,Zhao J,Wang W.An optimized oxygen system scheduling with electricity cost consideration in steel industry.IEEE/CAA Journal of Automatica Sinica,2017,4(2):216?222

    12 Wang Xiao-Hui.Research on the Power Load Simulation and Forecasting of Baosteel Branch[Master dissertation],Shanghai Jiao Tong University,China,2008.(王小輝.寶鋼分公司電力負(fù)荷模擬與預(yù)測(cè)研究[碩士學(xué)位論文],上海交通大學(xué),中國(guó),2008.)

    13 Ashok S.Peak-load management in steel plants.Applied Energy,2006,83(5):413?424

    14 Gao Yun-Long,Gao Feng,Pan Jin-Yan,Zhai Qiao-Zhu,Guan Xiao-Hong.Self-scheduling for electrical energy balance and output power control of energy-intensive enterprises.Proceedings of the CSEE,2010,30(19):76?83(高云龍,高峰,潘金艷,翟橋柱,管曉宏.高耗能企業(yè)關(guān)口平衡優(yōu)化調(diào)度及其輸出功率控制方式.中國(guó)電機(jī)工程學(xué)報(bào),2010,30(19):76?83)

    15 Liu Kun,Gao Feng,Zhai Qiao-Zhu,Wu Jiang,Guan Xiao-Hong,Wang Zhao-Jie,Zhang Hai-Feng.A self-power generation scheduling model under load demand and uncertainty of a by-product of gas production in enterprises microgrid.Proceedings of the CSEE,2014,34(13):2063?2070(劉坤,高峰,翟橋柱,吳江,管曉宏,王兆杰,張海峰.考慮負(fù)荷及煤氣量不確定性的企業(yè)微電網(wǎng)自發(fā)電調(diào)度模型.中國(guó)電機(jī)工程學(xué)報(bào),2014,34(13):2063?2070)

    16 Zhang Yu-Qing,Sun Yan-Guang.Function and operation frame research of intelligent power system for iron and steel enterprises.Metallurgical Industry Automation,2011,35(3):8?13(張玉慶,孫彥廣.鋼鐵企業(yè)智能電力系統(tǒng)功能與運(yùn)行架構(gòu)研究.冶金自動(dòng)化,2011,35(3):8?13)

    17 Li Dan.Optimization adjustment of large energy systems base on maximizing the integrate value of energy.Energy for Metallurgical Industry,2013,32(1):3?5,58(李丹.基于價(jià)值最大化的能源系統(tǒng)綜合調(diào)整.冶金能源,2013,32(1):3?5,58)

    18 Meng Hua.Study on Optimization of Unit Con figuration for Own Power Plant and Optimization of Gas Dispatching in Iron and Steel Enterprise[Ph.D.dissertation],Kunming University of Science and Technology,China,2013.(孟華.鋼鐵企業(yè)自備電廠機(jī)組配置優(yōu)化及煤氣優(yōu)化調(diào)度研究[博士學(xué)位論文],昆明理工大學(xué),中國(guó),2013.)

    19 Zhang Qi,Ti Wei,Du Tao,Cai Jiu-Ju.Coupling model of gas-steam-electricity and its application in steel works.CIESC Journal,2011,62(3):753?758(張琦,提威,杜濤,蔡九菊.鋼鐵企業(yè)富余煤氣–蒸汽–電力耦合模型及其應(yīng)用.化工學(xué)報(bào),2011,62(3):753?758)

    20 Sun Yang-Guang.Energy flow network information model in iron and steel enterprise and dynamic control of various energy.In:Proceedings of the 2010 National Energy and Environmental Protection Production Technology Conference.Jiujiang,Jiangxi,China:Chinese Society of Metal,2010.(孫彥廣.鋼鐵企業(yè)能量流網(wǎng)絡(luò)信息模型及多種能源介質(zhì)動(dòng)態(tài)調(diào)控.見(jiàn):2010年全國(guó)能源環(huán)保生產(chǎn)技術(shù)會(huì)議論文集.九江,江西,中國(guó):中國(guó)金屬學(xué)會(huì),2010.)

    21 Zeng Yu-Jiao.Research on Active and Reactive Power Optimization Scheduling of the Power System in Iron and Steel Enterprise[Ph.D.dissertation],Central Iron&Steel Research Institute,China,2015.(曾玉嬌.鋼鐵企業(yè)電力系統(tǒng)有功和無(wú)功優(yōu)化調(diào)度問(wèn)題的研究[博士學(xué)位論文],鋼鐵研究總院,中國(guó),2015.)

    22 ZengYu-Jiao,Sun Yan-Guang.Multi-period optimal scheduling of steam power system for iron and steel industry.In:Proceedings of the 25th China Process Control Conference.Dalian,Liaoning,China:China Automation Society Process Control Specialized Committee,2014.(曾玉嬌,孫彥廣.鋼鐵企業(yè)蒸汽–電力系統(tǒng)多時(shí)段優(yōu)化調(diào)度.見(jiàn):第25屆中國(guó)過(guò)程控制會(huì)議論文集.大連,遼寧,中國(guó):中國(guó)自動(dòng)化學(xué)會(huì)過(guò)程控制專業(yè)委員會(huì),2014.)

    23 He Jia-Yi,Ji Yang,Li Wen-Bing,Zhang Yun-Li.Research and software realization of network model and con figuration on energy simulation in iron and steel enterprises.Metallurgical Industry Automation,2012,36(1):7?12(何佳毅,紀(jì)揚(yáng),李文兵,張?jiān)评?鋼鐵企業(yè)能源系統(tǒng)網(wǎng)絡(luò)模型仿真及組態(tài)的研究與實(shí)現(xiàn).冶金自動(dòng)化,2012,36(1):7?12)

    24 Mardan N.Combining Simulation and Optimization for Improved Decision Support on Energy Efficiency in Industry[Ph.D.dissertation],Linkping University,Sweden,2012.

    25 Yamamoto T,Nakagawa T.A vision of energy structure for integrated steel works of future.Transactions of the Iron and Steel Institute of Japan,1983,23(10):862?892

    26 Ohkuma R,Ikegami K,Yasunaga S.Energy problems and energy control system in the Japanese steel industry.A I I E Transactions,1981,13(2):164?174

    27 Gou H,Olynyk S.A corporate mass and energy simulation model for an integrated steel plant.Iron&Steel Technology,2007,4(4):141?150

    28 Liu Liu,Gan Yong,Zhang Jiang-Ling,Li Jing.Research on energy recycling at integrated steel companies.Iron&Steel,2006,41(6):1?4(劉瀏,干勇,張江玲,李菁.鋼鐵聯(lián)合企業(yè)能源循環(huán)利用的分析研究.鋼鐵,2006,41(6):1?4)

    29 Wang Jian-Jun,Cai Jiu-Ju,Zhang Qi,Wu Fu-Zhong,Chen Chun-Xia.Study on energy- flow modelling in iron and steel enterprise.China Metallurgy,2006,16(5):48?52(王建軍,蔡九菊,張琦,吳復(fù)忠,陳春霞.鋼鐵企業(yè)能量流模型化研究.中國(guó)冶金,2006,16(5):48?52)

    30 Qiu Xiao-Lei,Meng Qing-Yu,Hong Xin.Study on mathematical model of process energy consumption of BF-LD process.Energy for Metallurgical Industry,2007,26(3):3?6,53(仇曉磊,孟慶玉,洪新.鋼鐵生產(chǎn)長(zhǎng)流程工序能耗數(shù)學(xué)模型研究.冶金能源,2007,26(3):3?6,53)

    31 Zeng Y J,Sun Y G.An improved particle swarm optimization for the combined heat and power dynamic economic dispatch problem.Electric Power Components and Systems,2014,42(15):1700?1716 and Design Institute of Metallurgical Industry in 2009.Her research interest covers energy management and optimization in iron and steel industry.)

    猜你喜歡
    優(yōu)化模型
    一半模型
    超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
    民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    看十八女毛片水多多多| 亚洲色图综合在线观看| 国产高清国产精品国产三级 | 2018国产大陆天天弄谢| 午夜老司机福利剧场| 韩国高清视频一区二区三区| 精品久久久久久久人妻蜜臀av| 欧美激情久久久久久爽电影| 欧美bdsm另类| 一级二级三级毛片免费看| 国产老妇伦熟女老妇高清| 中文字幕久久专区| 免费看光身美女| 天堂中文最新版在线下载 | 2021少妇久久久久久久久久久| 亚洲最大成人手机在线| av在线app专区| 国产色婷婷99| 九草在线视频观看| 99久久精品热视频| 日韩欧美精品免费久久| 在线亚洲精品国产二区图片欧美 | 国产av国产精品国产| 国产成人免费无遮挡视频| 天天躁日日操中文字幕| 亚洲,一卡二卡三卡| 综合色丁香网| 啦啦啦中文免费视频观看日本| 久久人人爽av亚洲精品天堂 | 69人妻影院| 一区二区三区四区激情视频| 能在线免费看毛片的网站| 最近手机中文字幕大全| 麻豆成人午夜福利视频| 在线亚洲精品国产二区图片欧美 | 亚洲不卡免费看| 九九爱精品视频在线观看| 亚洲va在线va天堂va国产| 午夜免费鲁丝| 国产精品不卡视频一区二区| 国产女主播在线喷水免费视频网站| 人妻夜夜爽99麻豆av| 中文字幕av成人在线电影| 亚洲成色77777| 免费高清在线观看视频在线观看| 亚洲av一区综合| 狂野欧美激情性bbbbbb| 美女高潮的动态| 乱系列少妇在线播放| 人妻少妇偷人精品九色| 免费在线观看成人毛片| 如何舔出高潮| 午夜精品一区二区三区免费看| 日韩一区二区视频免费看| 联通29元200g的流量卡| 国产高清三级在线| 最近最新中文字幕大全电影3| 午夜精品国产一区二区电影 | 欧美激情久久久久久爽电影| 熟女电影av网| 永久免费av网站大全| 久久99蜜桃精品久久| 日本色播在线视频| 五月玫瑰六月丁香| 嫩草影院精品99| 观看美女的网站| 热re99久久精品国产66热6| 日韩在线高清观看一区二区三区| 日本午夜av视频| 啦啦啦中文免费视频观看日本| 国产精品蜜桃在线观看| 成人亚洲精品一区在线观看 | 欧美精品国产亚洲| 99久国产av精品国产电影| 99热这里只有是精品50| 国产永久视频网站| 狂野欧美激情性xxxx在线观看| 久久精品国产亚洲网站| 欧美日本视频| 三级国产精品欧美在线观看| 精品人妻视频免费看| 精华霜和精华液先用哪个| 免费电影在线观看免费观看| 国国产精品蜜臀av免费| 午夜激情久久久久久久| 国产精品三级大全| 99久久中文字幕三级久久日本| av一本久久久久| 视频中文字幕在线观看| 91久久精品国产一区二区成人| 国产精品嫩草影院av在线观看| 成年免费大片在线观看| 新久久久久国产一级毛片| 97超碰精品成人国产| av.在线天堂| 中文精品一卡2卡3卡4更新| 五月天丁香电影| 一级毛片aaaaaa免费看小| 国产乱人视频| 欧美性猛交╳xxx乱大交人| 亚洲精品成人av观看孕妇| kizo精华| 亚洲精品成人久久久久久| 亚洲,一卡二卡三卡| 国产久久久一区二区三区| av福利片在线观看| 欧美精品一区二区大全| 国产精品福利在线免费观看| 好男人视频免费观看在线| 直男gayav资源| 美女脱内裤让男人舔精品视频| 少妇被粗大猛烈的视频| 美女国产视频在线观看| 夫妻午夜视频| 欧美精品国产亚洲| 亚洲aⅴ乱码一区二区在线播放| 国产av国产精品国产| 麻豆成人午夜福利视频| 又爽又黄无遮挡网站| 国产乱人偷精品视频| 人妻 亚洲 视频| 五月开心婷婷网| 中文字幕av成人在线电影| 99久久精品热视频| 听说在线观看完整版免费高清| 国语对白做爰xxxⅹ性视频网站| 国产精品国产三级国产专区5o| 日韩视频在线欧美| 日韩国内少妇激情av| 欧美日韩一区二区视频在线观看视频在线 | 国产综合懂色| 少妇的逼水好多| 看十八女毛片水多多多| 亚洲av国产av综合av卡| 亚洲怡红院男人天堂| 激情 狠狠 欧美| 久久国内精品自在自线图片| 免费看av在线观看网站| 亚洲人成网站在线播| 国产淫语在线视频| 91aial.com中文字幕在线观看| 26uuu在线亚洲综合色| 国产精品三级大全| 一级黄片播放器| 国产伦理片在线播放av一区| 亚洲av.av天堂| 一本一本综合久久| 亚洲精品456在线播放app| 欧美 日韩 精品 国产| 久久久a久久爽久久v久久| h日本视频在线播放| 最后的刺客免费高清国语| 男人爽女人下面视频在线观看| 视频区图区小说| 国产精品99久久久久久久久| 熟妇人妻不卡中文字幕| 日本一二三区视频观看| 国产精品久久久久久av不卡| 国产欧美另类精品又又久久亚洲欧美| 国产免费一级a男人的天堂| 美女脱内裤让男人舔精品视频| 看黄色毛片网站| 97超碰精品成人国产| 亚洲无线观看免费| 我的老师免费观看完整版| 亚洲av中文字字幕乱码综合| 成人一区二区视频在线观看| 国产一区有黄有色的免费视频| 久久精品国产a三级三级三级| 国产精品一二三区在线看| 午夜免费鲁丝| 国产欧美日韩一区二区三区在线 | 五月天丁香电影| 国产在视频线精品| 在线 av 中文字幕| 69av精品久久久久久| 国产av国产精品国产| 精品久久久精品久久久| 亚洲欧美清纯卡通| 日韩 亚洲 欧美在线| 欧美3d第一页| 丝袜喷水一区| 丝袜喷水一区| 国产高清三级在线| 在线观看人妻少妇| 欧美成人精品欧美一级黄| 狂野欧美白嫩少妇大欣赏| 韩国av在线不卡| 亚洲精品国产av蜜桃| 男男h啪啪无遮挡| 大香蕉97超碰在线| 毛片女人毛片| 丰满少妇做爰视频| 狂野欧美激情性bbbbbb| 中文字幕制服av| 国产免费又黄又爽又色| 极品少妇高潮喷水抽搐| 另类亚洲欧美激情| 在线播放无遮挡| 成年版毛片免费区| 亚洲精品国产成人久久av| 免费看日本二区| 99精国产麻豆久久婷婷| 精品一区在线观看国产| 极品少妇高潮喷水抽搐| 日韩欧美精品v在线| 日韩欧美精品免费久久| 日韩成人伦理影院| 亚洲国产欧美在线一区| 亚洲国产成人一精品久久久| 国产高潮美女av| 一级毛片久久久久久久久女| 黄色欧美视频在线观看| 看非洲黑人一级黄片| 国内精品宾馆在线| 九草在线视频观看| 97在线视频观看| 国产精品久久久久久久久免| av在线蜜桃| 午夜亚洲福利在线播放| 亚洲三级黄色毛片| 日韩不卡一区二区三区视频在线| 午夜激情福利司机影院| 一级毛片黄色毛片免费观看视频| 欧美日韩视频高清一区二区三区二| 亚洲精品一区蜜桃| 中文欧美无线码| 亚洲欧美清纯卡通| 大香蕉久久网| 国产成人福利小说| av在线app专区| 国产乱人视频| 18禁在线无遮挡免费观看视频| 秋霞伦理黄片| 亚洲av成人精品一二三区| 丝袜脚勾引网站| 国产毛片a区久久久久| 国产成人精品婷婷| 亚洲不卡免费看| 久久热精品热| 亚洲高清免费不卡视频| 国产伦理片在线播放av一区| 欧美精品国产亚洲| 欧美成人午夜免费资源| 久久99蜜桃精品久久| 国产在视频线精品| 欧美激情久久久久久爽电影| 国产亚洲最大av| av线在线观看网站| 97人妻精品一区二区三区麻豆| 国内精品宾馆在线| 夫妻性生交免费视频一级片| 日本av手机在线免费观看| 又黄又爽又刺激的免费视频.| 国产亚洲5aaaaa淫片| 男插女下体视频免费在线播放| 欧美日本视频| 免费不卡的大黄色大毛片视频在线观看| 99久久精品国产国产毛片| 久久久色成人| 狂野欧美白嫩少妇大欣赏| 在现免费观看毛片| 久久久久久伊人网av| 久久精品熟女亚洲av麻豆精品| 欧美丝袜亚洲另类| 好男人视频免费观看在线| 成人无遮挡网站| 亚洲精品国产av成人精品| 在线 av 中文字幕| 51国产日韩欧美| 色视频www国产| 国产欧美日韩精品一区二区| 中文资源天堂在线| av福利片在线观看| 成人黄色视频免费在线看| 国产成人免费观看mmmm| 少妇裸体淫交视频免费看高清| 国产亚洲av嫩草精品影院| 一级黄片播放器| 男人爽女人下面视频在线观看| 亚洲成人精品中文字幕电影| 精品久久久久久电影网| 在线观看三级黄色| 国产精品人妻久久久影院| 久久久精品欧美日韩精品| 校园人妻丝袜中文字幕| 精品久久久精品久久久| 免费黄网站久久成人精品| 夜夜爽夜夜爽视频| 亚洲性久久影院| 午夜福利在线观看免费完整高清在| av播播在线观看一区| 麻豆国产97在线/欧美| 午夜老司机福利剧场| 国产精品伦人一区二区| 成人午夜精彩视频在线观看| 欧美区成人在线视频| 91aial.com中文字幕在线观看| 国产精品嫩草影院av在线观看| 成人免费观看视频高清| 亚洲av一区综合| 久久这里有精品视频免费| 99热网站在线观看| 日本三级黄在线观看| 免费看a级黄色片| 亚洲欧美中文字幕日韩二区| 肉色欧美久久久久久久蜜桃 | 成人综合一区亚洲| 亚洲av.av天堂| 欧美 日韩 精品 国产| 久久精品人妻少妇| 亚洲国产色片| 日本免费在线观看一区| 一本久久精品| 亚洲av二区三区四区| 久久综合国产亚洲精品| 尾随美女入室| 亚洲精品国产成人久久av| 午夜精品一区二区三区免费看| 少妇人妻 视频| 亚洲成人一二三区av| 成人毛片a级毛片在线播放| 亚洲婷婷狠狠爱综合网| 成人免费观看视频高清| 少妇猛男粗大的猛烈进出视频 | 一区二区三区免费毛片| 网址你懂的国产日韩在线| 色视频在线一区二区三区| 赤兔流量卡办理| 精品久久久久久久久av| 一级毛片aaaaaa免费看小| 国产精品人妻久久久影院| 在线观看人妻少妇| 亚洲精品乱久久久久久| 国产精品麻豆人妻色哟哟久久| 国产精品av视频在线免费观看| 国产精品蜜桃在线观看| 亚洲精品成人久久久久久| 欧美3d第一页| 国产真实伦视频高清在线观看| 亚洲欧美日韩另类电影网站 | 国产午夜精品久久久久久一区二区三区| 国产人妻一区二区三区在| 国产精品久久久久久久电影| 亚洲精品国产av成人精品| 久久久久久久久久人人人人人人| 天天躁日日操中文字幕| 国产老妇女一区| 国产成人a∨麻豆精品| 黄色配什么色好看| 国产免费视频播放在线视频| 少妇高潮的动态图| 色播亚洲综合网| 国产成人freesex在线| 在线观看国产h片| a级一级毛片免费在线观看| 大码成人一级视频| 国产探花在线观看一区二区| 精品久久久久久久久av| 人人妻人人看人人澡| 夫妻午夜视频| 成人午夜精彩视频在线观看| 激情 狠狠 欧美| 老女人水多毛片| 久久久亚洲精品成人影院| 亚洲最大成人手机在线| 黄色配什么色好看| 又粗又硬又长又爽又黄的视频| 男人和女人高潮做爰伦理| 国产免费又黄又爽又色| 免费观看a级毛片全部| 亚洲va在线va天堂va国产| 一级毛片aaaaaa免费看小| 91aial.com中文字幕在线观看| 亚洲婷婷狠狠爱综合网| 国产探花在线观看一区二区| 亚洲精品视频女| 免费大片18禁| 狠狠精品人妻久久久久久综合| 99视频精品全部免费 在线| 激情五月婷婷亚洲| 一级av片app| 十八禁网站网址无遮挡 | 两个人的视频大全免费| 熟女人妻精品中文字幕| 国产黄色免费在线视频| 欧美3d第一页| 日韩中字成人| 国产在线一区二区三区精| 精品久久久久久电影网| 国产一区有黄有色的免费视频| 国产一区亚洲一区在线观看| 成年女人在线观看亚洲视频 | 人妻制服诱惑在线中文字幕| 久久韩国三级中文字幕| 亚洲成人av在线免费| 国产免费一区二区三区四区乱码| 夜夜看夜夜爽夜夜摸| 禁无遮挡网站| 中文天堂在线官网| 国产视频内射| 看非洲黑人一级黄片| 亚洲精品乱码久久久v下载方式| 精品一区二区免费观看| tube8黄色片| 黑人高潮一二区| 天堂网av新在线| 国产av不卡久久| 成年女人在线观看亚洲视频 | 伦理电影大哥的女人| 亚洲国产av新网站| 熟妇人妻不卡中文字幕| 久久久久久伊人网av| 国产高清三级在线| 久久久久久久午夜电影| 国产精品99久久久久久久久| 国产成人免费无遮挡视频| 老司机影院毛片| 免费观看的影片在线观看| 大话2 男鬼变身卡| 国产美女午夜福利| 边亲边吃奶的免费视频| 51国产日韩欧美| 午夜福利在线在线| 成人亚洲精品av一区二区| 久久久久精品久久久久真实原创| 免费黄频网站在线观看国产| 日本猛色少妇xxxxx猛交久久| 小蜜桃在线观看免费完整版高清| 美女xxoo啪啪120秒动态图| 日韩,欧美,国产一区二区三区| 欧美性感艳星| 成人亚洲精品一区在线观看 | 久久久久精品性色| 久久久亚洲精品成人影院| 麻豆精品久久久久久蜜桃| 日韩欧美精品v在线| 91精品国产九色| 国产精品av视频在线免费观看| 久久久a久久爽久久v久久| 高清av免费在线| 在线播放无遮挡| 日本爱情动作片www.在线观看| 一个人看的www免费观看视频| 一区二区三区乱码不卡18| 大话2 男鬼变身卡| 久久精品国产亚洲网站| 黄片wwwwww| 新久久久久国产一级毛片| 亚洲av福利一区| 午夜视频国产福利| 国产综合精华液| 白带黄色成豆腐渣| 亚洲av免费高清在线观看| 中文字幕av成人在线电影| 韩国高清视频一区二区三区| 男女下面进入的视频免费午夜| videos熟女内射| 高清av免费在线| 日韩在线高清观看一区二区三区| 卡戴珊不雅视频在线播放| 插逼视频在线观看| 下体分泌物呈黄色| 国产精品国产三级专区第一集| 国产亚洲91精品色在线| 亚洲激情五月婷婷啪啪| 毛片女人毛片| 在线天堂最新版资源| 一区二区三区乱码不卡18| 国产高清不卡午夜福利| 最近中文字幕高清免费大全6| 人人妻人人看人人澡| 一边亲一边摸免费视频| 伊人久久国产一区二区| 亚洲精品成人久久久久久| 久久99热6这里只有精品| 国内少妇人妻偷人精品xxx网站| 欧美日韩国产mv在线观看视频 | 免费人成在线观看视频色| 国产精品一及| 免费在线观看成人毛片| 天天一区二区日本电影三级| 在线观看国产h片| 亚洲av成人精品一区久久| 午夜日本视频在线| 婷婷色综合www| 在线a可以看的网站| 久久午夜福利片| 少妇高潮的动态图| 69av精品久久久久久| 人妻制服诱惑在线中文字幕| 国产黄片视频在线免费观看| 国产伦理片在线播放av一区| 亚洲av一区综合| 国产高清不卡午夜福利| 51国产日韩欧美| 精品久久久久久久久亚洲| 97在线人人人人妻| 亚洲精品国产色婷婷电影| 九草在线视频观看| 在线 av 中文字幕| 欧美日韩视频精品一区| av线在线观看网站| 精品久久久久久电影网| av线在线观看网站| 欧美三级亚洲精品| 亚洲综合色惰| 国产色爽女视频免费观看| 一本色道久久久久久精品综合| 国产综合懂色| 欧美精品国产亚洲| 国产伦精品一区二区三区四那| 亚洲国产欧美人成| 天堂网av新在线| 国产成年人精品一区二区| 国产日韩欧美亚洲二区| 国产伦精品一区二区三区视频9| 女的被弄到高潮叫床怎么办| 精品久久久久久久久亚洲| 嘟嘟电影网在线观看| 亚洲欧美精品专区久久| 午夜爱爱视频在线播放| 男的添女的下面高潮视频| 成年人午夜在线观看视频| 亚洲欧美精品自产自拍| 精品久久久久久久末码| 精品国产一区二区三区久久久樱花 | 91精品国产九色| 2021天堂中文幕一二区在线观| 久久精品夜色国产| 日韩av不卡免费在线播放| 韩国av在线不卡| 哪个播放器可以免费观看大片| 亚洲国产成人一精品久久久| 日韩不卡一区二区三区视频在线| 亚洲国产精品国产精品| 国产成人免费观看mmmm| 国产亚洲av片在线观看秒播厂| 爱豆传媒免费全集在线观看| 国产欧美日韩精品一区二区| 国产一区二区亚洲精品在线观看| 黄色欧美视频在线观看| 精品一区在线观看国产| 赤兔流量卡办理| 少妇裸体淫交视频免费看高清| 伦精品一区二区三区| 国产视频首页在线观看| 日日摸夜夜添夜夜添av毛片| 在线观看美女被高潮喷水网站| 久久鲁丝午夜福利片| 性插视频无遮挡在线免费观看| 久久精品国产亚洲av涩爱| 亚洲色图综合在线观看| 人人妻人人澡人人爽人人夜夜| 一级毛片我不卡| 人妻夜夜爽99麻豆av| 成年免费大片在线观看| 欧美bdsm另类| 久久久久久久亚洲中文字幕| 国国产精品蜜臀av免费| 日日摸夜夜添夜夜添av毛片| 少妇猛男粗大的猛烈进出视频 | av卡一久久| 全区人妻精品视频| 国产亚洲精品久久久com| 永久免费av网站大全| 最近2019中文字幕mv第一页| 一区二区三区精品91| 国产精品久久久久久久久免| 成人国产麻豆网| 国产高清三级在线| 国产精品三级大全| kizo精华| 久久6这里有精品| 成年女人看的毛片在线观看| 高清欧美精品videossex| 三级国产精品欧美在线观看| 狂野欧美白嫩少妇大欣赏| 日本与韩国留学比较| 晚上一个人看的免费电影| 日韩伦理黄色片| 真实男女啪啪啪动态图| 各种免费的搞黄视频| 成人美女网站在线观看视频| 日本一二三区视频观看| av天堂中文字幕网| 亚洲国产最新在线播放| 亚洲国产高清在线一区二区三| 卡戴珊不雅视频在线播放| 国产精品久久久久久久电影| 国产av国产精品国产| 亚洲av二区三区四区| 禁无遮挡网站| 欧美精品一区二区大全| 99热全是精品| 看非洲黑人一级黄片| 精品亚洲乱码少妇综合久久| 欧美精品国产亚洲| 中国美白少妇内射xxxbb| 国精品久久久久久国模美| 欧美最新免费一区二区三区| 好男人在线观看高清免费视频| 亚洲精品视频女| 极品教师在线视频| 国产伦精品一区二区三区视频9| 亚洲人成网站高清观看| 亚洲精品国产成人久久av| 久久鲁丝午夜福利片| 日韩欧美精品免费久久| 纵有疾风起免费观看全集完整版| 韩国高清视频一区二区三区| 欧美bdsm另类| 精品久久久噜噜| 精品一区在线观看国产| 黄色配什么色好看|