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

    數(shù)據(jù)中心綜合能源系統(tǒng)優(yōu)化運行研究綜述

    2022-11-05 08:38:04馬曉燕穆云飛李樹榮姜欣陽李華陳長金
    電力建設(shè) 2022年11期
    關(guān)鍵詞:數(shù)據(jù)中心耦合建模

    馬曉燕,穆云飛,李樹榮,姜欣陽,李華,陳長金

    (1. 智能電網(wǎng)教育部重點實驗室(天津大學(xué)),天津市300072;2. 國網(wǎng)河北省電力有限公司雄安新區(qū)供電公司,河北省雄安新區(qū) 071000;3. 國網(wǎng)河北省電力有限公司涉縣供電分公司,河北省邯鄲市 056400;4. 國網(wǎng)河北省電力有限公司電力培訓(xùn)中心,石家莊市 050000)

    0 引 言

    社會經(jīng)濟(jì)的快速發(fā)展和化石能源的大規(guī)模開發(fā)利用導(dǎo)致資源枯竭、環(huán)境污染等問題的日益嚴(yán)重[1-2],同時傳統(tǒng)的能源系統(tǒng)優(yōu)化和運行局限于電/氣/冷/熱等單一能源系統(tǒng)內(nèi)部,且各種能源形式分開且相對獨立,這大大限制了能源的綜合協(xié)調(diào)利用,不利于用能能效的提升和節(jié)能減排等[3]。如何在確保人類社會能源可持續(xù)供應(yīng)的同時盡量減少環(huán)境污染,實現(xiàn)對各類能源的高效利用、清潔轉(zhuǎn)換和互補(bǔ)融合已成為近年來研究的熱點。為此,習(xí)近平總書記提出了推動能源生產(chǎn)和消費革命、加強(qiáng)國際合作等重要論述,明確要求建設(shè)清潔低碳、能源安全、科技創(chuàng)新、擴(kuò)大開放的現(xiàn)代能源體系。針對上述問題和背景,能源互聯(lián)網(wǎng)(energy internet,EI)和綜合能源系統(tǒng)(integrated energy system,IES)[4-6]的概念應(yīng)運而生。

    EI旨在將能源和互聯(lián)網(wǎng)深度融合,強(qiáng)調(diào)信息流和能量流之間的交互作用,注重多能源系統(tǒng)間的泛在互聯(lián)[7-8]。IES是EI的物理載體,特指在“規(guī)劃、建設(shè)和運行等過程中,通過對能源的產(chǎn)生、傳輸、分配、轉(zhuǎn)換、存儲、消費等環(huán)節(jié)進(jìn)行有機(jī)協(xié)調(diào)與優(yōu)化后,逐步形成以電力系統(tǒng)為核心,涵蓋供電/氣/冷/熱等能源產(chǎn)供銷一體化的復(fù)雜系統(tǒng)”[3,6]。IES集電力系統(tǒng)、天然氣系統(tǒng)和熱力系統(tǒng)等于一體,能夠滿足多種能量需求,通過多種能源間的科學(xué)調(diào)度實現(xiàn)對多種能源的綜合管理和經(jīng)濟(jì)調(diào)度,是多能互補(bǔ)、能源梯級利用等理念實現(xiàn)應(yīng)用的關(guān)鍵,同時它能夠提高能源的利用效率、促進(jìn)可再生能源消納、緩解環(huán)境污染、實現(xiàn)重要用戶的不間斷供能等[9]。

    近年來IES已經(jīng)得到各國政府和工業(yè)界的支持,并已列入美國[10]、歐盟[11]、韓日[12]等發(fā)達(dá)國家和地區(qū)能源發(fā)展領(lǐng)域的重點研究方向。國際上的專家學(xué)者于2014年成立能源系統(tǒng)集成國際聯(lián)合研究會以解決IES的協(xié)調(diào)與優(yōu)化問題;而中國近年來已從國家層面對IES進(jìn)行積極探索和研究,2015年7月國務(wù)院印發(fā)《關(guān)于積極推進(jìn)“互聯(lián)網(wǎng)+”行動的指導(dǎo)意見》,并提出“互聯(lián)網(wǎng)+”智慧能源行動的設(shè)想[13]。2016年7月,國家發(fā)展改革委聯(lián)合國家能源局發(fā)布《關(guān)于推進(jìn)多能互補(bǔ)集成優(yōu)化示范工程建設(shè)的實施意見》[14],強(qiáng)調(diào)了創(chuàng)新管理體制和商業(yè)模式,以加快推進(jìn)示范工程建設(shè),提高能源系統(tǒng)效率,帶動有效投資,促進(jìn)經(jīng)濟(jì)穩(wěn)定增長。2017年10月,國家電網(wǎng)有限公司下發(fā)《關(guān)于在各省公司開展綜合能源服務(wù)業(yè)務(wù)的意見》,同時國家電網(wǎng)有限公司與南方電網(wǎng)公司陸續(xù)提出要做強(qiáng)、做優(yōu)、做大綜合能源服務(wù)業(yè)務(wù),推動公司由電能供應(yīng)商向綜合能源服務(wù)商轉(zhuǎn)變。

    隨著互聯(lián)網(wǎng)技術(shù)的飛速發(fā)展,數(shù)據(jù)中心的規(guī)模和數(shù)量高速增加,截至2019年年底,我國數(shù)據(jù)中心機(jī)架總規(guī)模達(dá)到314.5萬架,其中大型及超大型數(shù)據(jù)中心占比達(dá)75%[15]。未來以數(shù)據(jù)中心為應(yīng)用典型的信息能源耦合IES的物質(zhì)流、能量流、信息流三流特性各異但高度耦合的特性[16-18],為IES耦合建模與協(xié)同控制帶來新的困難與挑戰(zhàn)。同時數(shù)據(jù)中心的機(jī)房散熱密度高、空調(diào)設(shè)備不間斷運行、服務(wù)器等IT設(shè)備運行功耗和空調(diào)制冷功耗的時變特性[19],導(dǎo)致機(jī)房內(nèi)部局部過熱成為一個比較突出的問題,影響著整個IES的安全穩(wěn)定運行及能源利用效率[20]。

    目前,國內(nèi)外專家學(xué)者對IES開展了詳細(xì)研究,包括IES的建模、運行、規(guī)劃設(shè)計、優(yōu)化、評估等。其中IES的優(yōu)化運行是通過協(xié)調(diào)控制和優(yōu)化管理各類分布式可再生能源、分布式儲能、分布式冷/熱/電聯(lián)供系統(tǒng)與電網(wǎng)絡(luò)等的運行特性,以滿足用戶的用冷/熱/電等需求,確保以電網(wǎng)為骨干網(wǎng)的IES運行可靠性,以此有效解決系統(tǒng)的經(jīng)濟(jì)、環(huán)保和穩(wěn)定運行問題以及因電/氣/冷/熱等能源種類多樣和數(shù)量眾多造成的耦合管理等難題[21],提高系統(tǒng)內(nèi)部的供電可靠性和電能質(zhì)量,減少能源在傳輸和分配過程中的損耗等。但已有研究大多針對某一特定的或假定的IES進(jìn)行研究,鮮有對基于數(shù)據(jù)中心的IES運行優(yōu)化進(jìn)行系統(tǒng)性的梳理和總結(jié),同時鮮有研究采用IES降低數(shù)據(jù)中心制冷耗電量的策略。在保證數(shù)據(jù)中心正常運行前提下,研究基于數(shù)據(jù)中心的IES優(yōu)化運行,對于降低制冷設(shè)備能耗、提高能源利用效率具有重要意義。

    為此,本文針對數(shù)據(jù)中心IES優(yōu)化運行這一問題對現(xiàn)有研究進(jìn)行綜述并給出相應(yīng)思考。在國內(nèi)外IES相關(guān)研究的基礎(chǔ)上,本文首先引入IES的典型研究框架以及基于數(shù)據(jù)中心的IES傳統(tǒng)研究框架;其次,總結(jié)分析多能源系統(tǒng)中的通用系統(tǒng)建模和含能源集線器的通用建模及擴(kuò)展建模;然后,對IES及基于數(shù)據(jù)中心的IES優(yōu)化運行國內(nèi)外研究現(xiàn)狀進(jìn)行分析總結(jié),并提出一種多能源輸入-多能源耦合-多能源輸出的運行方式;最后,梳理當(dāng)前研究的不足,并對未來基于數(shù)據(jù)中心的IES建模和運行優(yōu)化等問題進(jìn)行展望。

    1 基于數(shù)據(jù)中心的IES研究框架

    1.1 IES典型物理架構(gòu)與設(shè)備

    IES根據(jù)地理因素與電/氣/冷/熱等各類能源的生產(chǎn)、傳輸、存儲和消費環(huán)節(jié)特性,分為跨區(qū)級、區(qū)域級和用戶級IES[22]。根據(jù)已有研究中各類IES的物理架構(gòu)[23-24],現(xiàn)以區(qū)域級IES為例總結(jié)出圖1所示的IES典型物理架構(gòu)和設(shè)備示意。圖中將電力、天然氣、供冷、供熱等多種能源環(huán)節(jié)相互耦合,通過系統(tǒng)間的有機(jī)協(xié)調(diào)和運行調(diào)度,以此實現(xiàn)多種能源的高效利用,提高系統(tǒng)的運行穩(wěn)定性。

    圖1 IES典型物理架構(gòu)與設(shè)備示意Fig.1 Typical physical architecture and equipment diagram of IES

    如圖1所示,根據(jù)電/氣/冷/熱系統(tǒng)間是否存在能流間的相互轉(zhuǎn)化將IES中的設(shè)備分為:獨立型運行設(shè)備(可再生能源和儲能設(shè)備等)和耦合型運行設(shè)備。其中電-氣耦合元件為電轉(zhuǎn)氣(power to gas,P2G)技術(shù),電-熱耦合元件為電鍋爐、熱泵等,電-冷耦合元件為戶用空調(diào)、電制冷機(jī)、冰蓄冷空調(diào)等,電-氣-熱-冷耦合元件為冷/熱/電三聯(lián)供機(jī)組(主要由微型燃?xì)廨啓C(jī)、余熱鍋爐和吸收式制冷機(jī)組成)。通過各類轉(zhuǎn)換元件的相互耦合,各類能源之間可實現(xiàn)能量的相互轉(zhuǎn)換和協(xié)調(diào)運行。

    圖2 數(shù)據(jù)中心用電量分布圖Fig.2 Power consumption distribution of data center

    1.2 基于數(shù)據(jù)中心的IES研究框架

    數(shù)據(jù)中心內(nèi)部含有大量的計算機(jī)及電子設(shè)備進(jìn)行不間斷的運行,由此產(chǎn)生大量熱量。根據(jù)美國采暖制冷與空調(diào)工程師學(xué)會技術(shù)委員會9.9統(tǒng)計報告,數(shù)據(jù)中心用電量分布大致如圖2所示[25],其中空調(diào)制冷耗電約占數(shù)據(jù)中心總耗電量的31%。隨著設(shè)備的增多,IT設(shè)備集成度不斷提高,機(jī)房內(nèi)部局部過熱問題嚴(yán)重,為保證數(shù)據(jù)中心的正常運行,需及時排熱。數(shù)據(jù)中心供電系統(tǒng)如圖3所示,傳統(tǒng)電制冷數(shù)據(jù)中心

    圖3 傳統(tǒng)電制冷數(shù)據(jù)中心結(jié)構(gòu)Fig.3 Structure of data center with traditional electric refrigeration

    結(jié)構(gòu)采用可再生能源(風(fēng)力發(fā)電、光伏發(fā)電)、儲能設(shè)備及市電滿足數(shù)據(jù)中心的電負(fù)荷需求,采用電制冷設(shè)備(如空調(diào)制冷)滿足數(shù)據(jù)中心的冷負(fù)荷需求。

    2 IES的穩(wěn)態(tài)建模及組件分析

    IES中獨立型設(shè)備單元生產(chǎn)、傳輸、存儲電能、熱能或天然氣,相關(guān)設(shè)備主要包括光伏電池[26]、風(fēng)力發(fā)電機(jī)[27]、儲能電池[28]、熱力管網(wǎng)、天然氣管網(wǎng)等,目前對獨立型設(shè)備單元的建模研究較為成熟。隨著多能源供能設(shè)備的不斷組合,其耦合建模難度增加,連接關(guān)系更加復(fù)雜。為解決上述問題,有學(xué)者提出了能源集線器的概念,以研究IES的通用建模。本節(jié)將主要進(jìn)行IES耦合型設(shè)備元件建模和能源樞紐擴(kuò)展建模。

    2.1 IES耦合型設(shè)備元件建模

    2.1.1 電轉(zhuǎn)氣技術(shù)模型

    P2G技術(shù)以一種環(huán)境友好的方式通過電-氣網(wǎng)絡(luò)的深度耦合將電能轉(zhuǎn)換為化學(xué)能,并以H2或天然氣的方式存儲。目前各國的研究機(jī)構(gòu)、學(xué)者對P2G建模、控制和應(yīng)用等方面做了很多研究[29-31]。其中文獻(xiàn)[31]對電轉(zhuǎn)氣進(jìn)行了精細(xì)化建模,模型結(jié)構(gòu)如圖4所示,主要包括以下2個過程:

    圖4 電轉(zhuǎn)氣精細(xì)化模型結(jié)構(gòu)示意圖Fig.4 Structure diagram of refined model of P2G

    1)電轉(zhuǎn)氫氣。

    通過電解水產(chǎn)生H2和O2,生成的H2經(jīng)壓縮機(jī)加壓進(jìn)入儲氫罐儲存,再經(jīng)壓縮機(jī)加壓后用于燃料電池、工業(yè)設(shè)計或者交通運輸?shù)阮I(lǐng)域。

    2)電轉(zhuǎn)天然氣。

    在電轉(zhuǎn)氫氣基礎(chǔ)上,高溫高壓下經(jīng)薩巴蒂催化反應(yīng)將H2和CO2在甲烷反應(yīng)器合成CH4和水,生成的CH4可直接注入天然氣管道供給天然氣負(fù)荷或參與到燃?xì)廨啓C(jī)的運行中。此處CO2可由碳捕捉技術(shù)從大氣、電廠尾氣或沼氣中獲取,同時產(chǎn)生的天然氣由電能轉(zhuǎn)換而來,因而可削弱可再生能源間歇性的影響[32]。

    P2G技術(shù)主要包含電解槽、儲氫罐和甲烷反應(yīng)器三部分,相應(yīng)的通用模型可分別表示為[32-33]:

    EEL,H2,t=uEL,e,tf(PEL,e,t/PEL,rated)πEL,rated

    (1)

    EHT,t=EHT,t-1+EHT,in,t-EHT,out,t

    (2)

    PM,g,t=ηmEM,H2,t×4HL×1/κ

    (3)

    式中:EEL,H2,t、PEL,e,t、PEL,rated分別為t時段電解槽產(chǎn)氫量、輸入電功率和其電功率額定值;πel,rated為電解槽的額定容量;f(PEL,e,t/PEL,rated)為電解槽效率函數(shù);uEL,e,t為t時段電解槽啟動狀態(tài)二進(jìn)制變量,uEL,e,t為0或1分別表示電解槽的停機(jī)和啟動狀態(tài);EHT,t、EHT,in,t、EHT,out,t分別為t時段的有效儲氫罐儲氫量、注入儲氫罐的有效氫量和儲氫罐的輸出氫量;PM,g,t、EM,H2,t分別為t時段甲烷反應(yīng)器生成的天然氣功率和輸入甲烷反應(yīng)器的有效氫量;ηm、HL分別為甲烷反應(yīng)器的運行效率和天然氣低熱值;κ為天然氣管道每立方米對應(yīng)的氣體質(zhì)量;4為H2生成天然氣的摩爾質(zhì)量折算系數(shù)。

    2.1.2 氫燃料電池出力模型

    氫燃料電池(fuel cell,F(xiàn)C)是將燃料的化學(xué)能轉(zhuǎn)變成電能的分布式電源,它也是氣-電耦合設(shè)備的典型代表。FC結(jié)構(gòu)簡單、壽命長、燃燒成本低、燃料利用率高、工作時無噪音、對環(huán)境影響小,其中經(jīng)由P2G技術(shù)產(chǎn)生的H2可供給FC進(jìn)行發(fā)電[34]。氫燃料電池的輸出功率可表示為:

    PFC,t=PH2,FC,tηH2ηFC

    (4)

    式中:PFC,t、PH2,FC,t分別為t時刻FC輸出的有功功率和消耗的H2功率;ηH2和ηFC分別為FC的H2利用率和發(fā)電效率。

    2.1.3 電制冷機(jī)出力模型

    電制冷機(jī)(電-冷耦合設(shè)備)通過消耗電能來進(jìn)行制冷[35],其制冷功率表達(dá)式為:

    Qele,t=Pele,tηele

    (5)

    式中:Qele,t為電制冷機(jī)的制冷功率;Pele,t為電制冷機(jī)消耗的電功率;ηele為電制冷機(jī)的能效比。

    2.1.4 冰蓄冷空調(diào)出力模型

    冰蓄冷空調(diào)主要包括制冷機(jī)和蓄冰槽兩部分,有制冷和蓄冰2種運行方式[36]。冰蓄冷裝置可以利用夜間低谷電價進(jìn)行制冰存儲在蓄冰槽中,利用白天用電高峰時進(jìn)行的融冰制冷量滿足負(fù)荷側(cè)供冷需求。冰蓄冷結(jié)構(gòu)按照制冷機(jī)和蓄冰設(shè)備的連接情況分為并聯(lián)式和串聯(lián)式結(jié)構(gòu),其中并聯(lián)式冰蓄冷結(jié)構(gòu)可聯(lián)合供冷也可單獨供冷,制冷機(jī)可同時制冰和供冷;串聯(lián)式冰蓄冷結(jié)構(gòu)中制冷機(jī)和蓄冰槽的冷量分布則存在一定的比例關(guān)系。以并聯(lián)式冰蓄冷結(jié)構(gòu)為例,其典型的數(shù)學(xué)模型如下所示[37-38]:

    1)制冷機(jī)。

    Qref,i,t=Pref,i,tEref,t

    (6)

    式中:Qref,i,t、Pref,i,t和Eref,t分別為t時段第i個制冷機(jī)的制冷功率、電功率和制冷機(jī)的制冷能效比。

    2)蓄冰槽。

    (7)

    式中:Stank,i,t為第i個蓄冰槽t時段的蓄冰容量;σtank為蓄冰槽的自損耗系數(shù);Ptank,i,t和Qtank,i,t分別為t時段第i個蓄冰槽的電功率和制冷功率;Etank和ηtank分別為蓄冰槽的制冰能效比和融冰效率;Tref、Tmelt分別為冰蓄冷裝置運行在蓄冰和融冰時段。

    3)其他。

    Qref,i,t+Qtank,i,t=Qice,i,t

    (8)

    Pref,i,t+Ptank,i,t=Pice,i,t

    (9)

    (10)

    式中:Qice,i,t、Pice,i,t分別為t時段第i個冰蓄冷裝置的總制冷功率和電功率。在實際運行中,制冷機(jī)和蓄冰槽的供冷量與通過制冷機(jī)和蓄冰槽乙二醇溶液的流量有關(guān)。

    2.1.5 冷熱電三聯(lián)供系統(tǒng)模型

    在IES中,電/氣/冷/熱耦合環(huán)節(jié)通常是通過冷熱電三聯(lián)供系統(tǒng)(combined cooling, heating and power, CCHP)實現(xiàn)的。CCHP主要以微型燃?xì)廨啓C(jī)(micro turbine, MT) (氣-電耦合設(shè)備)進(jìn)行發(fā)電,余熱鍋爐(氣-熱耦合設(shè)備)產(chǎn)生的余熱分別經(jīng)換熱器和吸收式制冷機(jī)(熱-冷耦合設(shè)備)來供熱和制冷。經(jīng)由P2G技術(shù)產(chǎn)生的CH4可供給CCHP系統(tǒng)進(jìn)行發(fā)電、制熱和制冷。CCHP的物理模型[24,39]為:

    1)MT出力模型。

    MT是將燃料的化學(xué)能轉(zhuǎn)化為電能的小型分布式能源系統(tǒng),通常以天然氣為燃料,具有體積小、污染少、可靠性高等性能。其輸出功率可表示為:

    PMT,t=Pgas,MT,tηgasηMT

    (11)

    式中:PMT,t、Pgas,MT,t分別為t時刻MT輸出的有功功率和消耗的天然氣功率;ηgas和ηMT分別為MT的天然氣利用率和發(fā)電效率。

    2)余熱鍋爐出力模型。

    余熱鍋爐輸出的總熱功率表達(dá)式為:

    (12)

    余熱鍋爐供給熱負(fù)荷的制熱功率為:

    Qheat,t=QHE,tηHEηhe

    (13)

    3)吸收式制冷機(jī)出力模型。

    吸收式制冷機(jī)通過利用工業(yè)生產(chǎn)過程中MT的余熱補(bǔ)足冷負(fù)荷需求,所產(chǎn)生的制冷功率為:

    Qcold,t=QHE,tηHEηAC

    (14)

    式中:Qcold,t、ηAC分別為t時刻吸收式制冷機(jī)制冷功率輸出和制冷能效比。

    2.2 能源集線器及其在IES優(yōu)化運行中的擴(kuò)展建模

    2.2.1 能源集線器

    隨著電鍋爐、微型燃?xì)廨啓C(jī)等供能元件的增加和負(fù)荷形式的變化,如何分析IES中的能量轉(zhuǎn)化、分配和存儲性質(zhì)是建模中遇到的主要問題。為此,2007年瑞士蘇黎世聯(lián)邦理工學(xué)院Martin Geidl等人發(fā)起了名為“未來能源愿景”的項目,提出了能源集線器的概念,即一種描述多能源載體中能源、網(wǎng)絡(luò)、負(fù)荷間轉(zhuǎn)換、調(diào)節(jié)和存儲電、熱、冷等多種形式能量的單元[40]。

    圖5為一個含EH的IES端口模型示意圖。無論系統(tǒng)多復(fù)雜,都可以表示成多種能源輸入-多種能源輸出的形式,而EH則用來分析能源轉(zhuǎn)化、傳輸和分配過程,即EH代表不同能源基礎(chǔ)設(shè)施與負(fù)荷之間的接口,其可通過一個能源耦合矩陣C來描述IES中的供能單元和終端需求的輸入輸出特性,以此實現(xiàn)對多種能源系統(tǒng)的協(xié)調(diào)互補(bǔ),并使得IES的規(guī)劃、優(yōu)化等可直接進(jìn)行分析求解。EH數(shù)學(xué)表達(dá)式和簡化表達(dá)式分別為:

    圖5 IES輸入-輸出端口示意圖Fig.5 Schematic of Input-output port for IES

    (15)

    Lout=CPin

    (16)

    式中:α,β, …,ω代表不同能源形式(電力、天然氣、熱能、可再生能源等);Pin和Lout分別代表不同能源的輸入和輸出矩陣;Cαβ為第β種能源輸出與第α種能源的輸入比值。

    2.2.2 能源集線器的擴(kuò)展模型

    IES系統(tǒng)中輸入側(cè)能源種類多樣,耦合關(guān)系復(fù)雜,目前已有通過對輸入側(cè)加入儲能裝置以對EH模型進(jìn)行擴(kuò)展的方法[22],同時文獻(xiàn)[41]對含多種元素的EH建模進(jìn)行了詳細(xì)的推廣,包括充分考慮儲能、需求響應(yīng)、電動汽車、新能源并網(wǎng)等元素。基于此,本文將在基于數(shù)據(jù)中心的IES優(yōu)化中考慮加入儲能,同時將耦合矩陣分解為能源分配和能源傳輸/轉(zhuǎn)化兩個步驟。

    能源分配表示輸入側(cè)的多種能源以一定比例分配到能源傳輸/轉(zhuǎn)化設(shè)備;能源傳輸/轉(zhuǎn)化指通過不同的途徑將輸入該設(shè)備的能源以一定的轉(zhuǎn)化效率進(jìn)行轉(zhuǎn)化,最后傳輸給負(fù)荷側(cè)。相應(yīng)的表達(dá)式為:

    (17)

    (18)

    L′out=ηυPin

    (19)

    式中:vαβ(0≤vαβ≤1)為輸入側(cè)能源α分配到轉(zhuǎn)化第β種能源的比例;ηαβ為能源α轉(zhuǎn)化為能源β的效率;v、η分別為分配矩陣和效率矩陣。

    當(dāng)在IES中加入電儲能和熱儲能后,輸入側(cè)為多種供能設(shè)備P和風(fēng)、光等獨立型運行單元PR;輸出端包括數(shù)據(jù)中心電負(fù)荷Le和冷負(fù)荷Lc;電儲能和熱儲能能量分別用EEB、EHB表示;儲能耦合矩陣分別用SEB、SHB表示。通過耦合矩陣建立的EH擴(kuò)展模型可表示為:

    (20)

    綜上,EH的通用模型及其擴(kuò)展模型通過靜態(tài)轉(zhuǎn)換將多種能源需求和多種能源供給連接,并且可通過優(yōu)化方法得出滿足負(fù)荷需求條件下的能源供給方式,因而適用于基于數(shù)據(jù)中心的IES優(yōu)化運行。

    3 數(shù)據(jù)中心IES優(yōu)化運行國內(nèi)外研究現(xiàn)狀

    3.1 IES優(yōu)化運行國內(nèi)外研究現(xiàn)狀

    IES優(yōu)化運行問題通??沙橄鬄槿缦卤磉_(dá)式:

    (21)

    式中:Ht、Htrans,t分別為t時刻系統(tǒng)中各能源單元的直接供熱功率和轉(zhuǎn)換熱功率;Pt+PR,t=Le,t、Ht+Htrans,t=Lc,t分別為電力和熱力負(fù)荷的等式約束條件;g(Pt,Ht,Htrans,t)≤0為系統(tǒng)的不等式約束條件,通常優(yōu)化運行中目標(biāo)函數(shù)為經(jīng)濟(jì)成本最低、穩(wěn)定性最高或數(shù)據(jù)中心制冷降耗最大等。

    為實現(xiàn)不同能源間的優(yōu)勢互補(bǔ)、可再生能源消納、節(jié)能減排以及實現(xiàn)IES的經(jīng)濟(jì)運行等,針對IES優(yōu)化調(diào)度的相關(guān)研究意義重大。按照系統(tǒng)研究對象來分,IES可分為:電-氣IES、電-熱IES和電-氣-熱IES。按照最優(yōu)化研究方法來說,目前研究熱點包括:單純的最優(yōu)化方法研究、考慮新能源及負(fù)荷隨機(jī)性的概率最優(yōu)能流和隨機(jī)最優(yōu)能流研究。

    3.1.1 電-氣IES優(yōu)化運行研究現(xiàn)狀

    天然氣作為一種重要的化石能源,在電力系統(tǒng)中的占比不斷增加。天然氣資源豐富、易于儲存、清潔高效[42],相較于風(fēng)、光等可再生能源的間歇性和波動性,天然氣可同時為天然氣網(wǎng)和電力網(wǎng)持續(xù)供能。天然氣網(wǎng)絡(luò)與電力網(wǎng)絡(luò)有相似的能量流動形式,且與電力網(wǎng)絡(luò)聯(lián)系最為緊密,目前對電-氣互聯(lián)IES的研究集中在協(xié)調(diào)運行層面。

    文獻(xiàn)[43]研究了電-氣互聯(lián)IES經(jīng)濟(jì)調(diào)度和優(yōu)化運行問題。天然氣網(wǎng)絡(luò)的隨機(jī)性以及運行約束易導(dǎo)致系統(tǒng)優(yōu)化結(jié)果過于樂觀,同時新能源的大規(guī)模接入會增加電力系統(tǒng)的不確定性,進(jìn)而電-氣互聯(lián)系統(tǒng)間的相關(guān)性也會對優(yōu)化調(diào)度方案產(chǎn)生影響,因此考慮隨機(jī)因素影響的概率最優(yōu)潮流受到極大重視[44]。文獻(xiàn)[45]計及負(fù)荷和風(fēng)電的不確定性,考慮電力網(wǎng)與天然氣網(wǎng)的運行約束以及輸入變量間的相關(guān)性,采用基于Nataf變換的點估計法進(jìn)行了電-氣互聯(lián)系統(tǒng)的概率最優(yōu)潮流計算。文獻(xiàn)[46]則考慮風(fēng)電場出力及電力、天然氣負(fù)荷的隨機(jī)性,通過建立基于機(jī)會約束規(guī)劃的電-氣互聯(lián)IES最優(yōu)潮流模型,并采用基于半不變量法和內(nèi)點法的啟發(fā)式算法進(jìn)行求解,最終表明該模型有助于提高電-氣互聯(lián)IES運行的安全性。但上述研究未考慮風(fēng)電等可再生能源并網(wǎng)對電-氣IES的影響。

    針對可再生能源并網(wǎng)問題,目前采用P2G技術(shù)將低谷時段難以消納的風(fēng)電轉(zhuǎn)化為易于大規(guī)模存儲的天然氣,以此實現(xiàn)電力-天然氣網(wǎng)絡(luò)的深度耦合,從而改善系統(tǒng)運行靈活性和存在的棄風(fēng)消納問題[47]。目前文獻(xiàn)[31]已對P2G技術(shù)進(jìn)行了精細(xì)化建模,但天然氣較高的運行成本嚴(yán)重影響了系統(tǒng)運行的經(jīng)濟(jì)性。針對上述問題,文獻(xiàn)[48-51]分別研究了電轉(zhuǎn)氣對電-氣互聯(lián)IES運行的影響、如何提升風(fēng)電調(diào)度的能力、如何解決系統(tǒng)凈負(fù)荷削峰填谷的問題以及如何兼顧系統(tǒng)運行經(jīng)濟(jì)性和風(fēng)電消納能力的問題,以此提高風(fēng)電利用率、降低系統(tǒng)運行成本等,但上述研究均未涉及需求側(cè)響應(yīng)。

    需求響應(yīng)可引導(dǎo)用戶理性用能,合理轉(zhuǎn)移高峰用電以此平移負(fù)荷峰谷,從而有效提升電力系統(tǒng)運行效率。目前對電-氣IES中引入需求響應(yīng)正處于初步探索階段,文獻(xiàn)[52]提出了一種考慮需求側(cè)負(fù)荷響應(yīng)及動態(tài)天然氣潮流的電-氣IES優(yōu)化調(diào)度新模型,并通過仿真驗證了需求響應(yīng)所帶來的經(jīng)濟(jì)性。

    同時考慮負(fù)荷和風(fēng)電的不確定性、風(fēng)電消納以及需求響應(yīng)的電-氣綜合能源系統(tǒng)目前鮮有研究,這是目前亟待解決的問題。

    3.1.2 電-熱IES優(yōu)化運行研究現(xiàn)狀

    除電-氣互聯(lián)外,電-熱互聯(lián)IES中電能和熱能是人類生產(chǎn)生活的基本需求,也是該系統(tǒng)的主要能源形式[53],熱系統(tǒng)的傳輸較慢、慣性較大,其延時效應(yīng)可等效為系統(tǒng)的儲能,特別是當(dāng)傳輸規(guī)模較大時可為系統(tǒng)提供較大的儲能容量,從而提升系統(tǒng)運行的靈活性[32]?,F(xiàn)有文獻(xiàn)已建立了有關(guān)熱網(wǎng)的通用模型[54]。

    在電-熱IES中,傳統(tǒng)電熱聯(lián)合調(diào)度中“以熱定電”的運行模式限制了系統(tǒng)的調(diào)節(jié)能力,導(dǎo)致棄風(fēng)現(xiàn)象高發(fā),除利用電-熱系統(tǒng)的互補(bǔ)特性提升系統(tǒng)消納可再生能源的能力外,常采用配置儲熱或電鍋爐的方式解耦“以熱定電”的運行方式[55-58]。文獻(xiàn)[55]對包含大容量儲熱的電-熱聯(lián)合系統(tǒng)的國內(nèi)外研究現(xiàn)狀、應(yīng)用前景和提高消納能力等關(guān)鍵技術(shù)方面做了綜述,但缺乏對具體實施方式和應(yīng)用效果的分析。文獻(xiàn)[56]構(gòu)建了包含儲熱的熱電聯(lián)產(chǎn)機(jī)組和電供熱系統(tǒng)的調(diào)度模型,將儲熱納入包含風(fēng)電的電力系統(tǒng)有功調(diào)度模型,進(jìn)而以最大化風(fēng)電消納為目標(biāo),通過算例仿真分析和驗證了儲熱提升風(fēng)電消納能力的有效性。文獻(xiàn)[57]提出了在熱電廠中配置大型電鍋爐來解耦熱電廠“以熱定電”約束以參與風(fēng)電調(diào)峰的運行方案,并對方案的國民經(jīng)濟(jì)性和潛在效益進(jìn)行了分析。文獻(xiàn)[58]則將電鍋爐和熱儲相結(jié)合進(jìn)行了風(fēng)電消納能力的分析。

    文獻(xiàn)[59-60]分別考慮熱力系統(tǒng)中供熱管道傳輸時間延遲等熱動態(tài)特性以及用戶供熱需求的柔性、供熱管道動態(tài)特性約束和儲熱裝置的模型,建立了包含儲熱的源-網(wǎng)-荷-儲協(xié)調(diào)優(yōu)化調(diào)度模型和熱電聯(lián)合調(diào)度模型,并分別驗證了相應(yīng)模型提高可再生能源消納能力的水平。除此之外,現(xiàn)有文獻(xiàn)[61]采用能量流法來統(tǒng)一電力和熱力系統(tǒng),提出基于能量流的電熱IES棄風(fēng)消納優(yōu)化調(diào)度模型,并在考慮換熱約束的前提下,協(xié)調(diào)優(yōu)化電-熱IES以實現(xiàn)棄風(fēng)的大規(guī)模消納。文獻(xiàn)[62-63]則從需求側(cè)角度出發(fā),分別考慮了用戶舒適性和多種負(fù)荷需求響應(yīng),建立了電-熱分時間尺度平衡的IES運行優(yōu)化模型和考慮電熱多種復(fù)合綜合需求響應(yīng)的園區(qū)微網(wǎng)IES以提高系統(tǒng)的運行經(jīng)濟(jì)性。不僅如此,文獻(xiàn)[64]在考慮IES經(jīng)濟(jì)性的同時還考慮了傳統(tǒng)發(fā)電機(jī)碳排放對環(huán)境的影響問題。

    目前綜合考慮系統(tǒng)經(jīng)濟(jì)運行、環(huán)境效益和可再生能源消納的電-熱IES優(yōu)化調(diào)度問題少有研究,這將是未來值得深入研究的問題。

    3.1.3 電-氣-熱IES優(yōu)化運行研究現(xiàn)狀

    IES中各類能源轉(zhuǎn)換設(shè)備如電鍋爐、熱泵、冷熱電聯(lián)產(chǎn)裝置、天然氣等使得電力網(wǎng)-天然氣網(wǎng)-熱力網(wǎng)相互耦合,以此實現(xiàn)多能源的互動和轉(zhuǎn)換。

    目前對于電-氣-熱IES建模問題集中在靜態(tài)建模研究。文獻(xiàn)[65]提出了集成電-氣-熱網(wǎng)的能源互聯(lián)網(wǎng)模型。文獻(xiàn)[66]形成了電-氣-熱系統(tǒng)的綜合能量流模型。文獻(xiàn)[67]通過不同能源網(wǎng)絡(luò)間的耦合機(jī)制建立了電-氣-熱綜合能源耦合系統(tǒng)模型,同時考慮階梯碳交易機(jī)制建立了經(jīng)濟(jì)成本最低的區(qū)域IES經(jīng)濟(jì)調(diào)度模型,并提出了包含能源集線器、熱網(wǎng)和天然氣網(wǎng)絡(luò)的優(yōu)化調(diào)度策略,最后采用果蠅優(yōu)化算法對天津某IES“以電定熱”模式、“以熱定電”模式和“綜合最優(yōu)”3種模式下的電力、熱力和天然氣網(wǎng)絡(luò)進(jìn)行算例分析,說明了在“綜合最優(yōu)”模式下的總運行成本更小,碳排放更低,同時也驗證了果蠅優(yōu)化算法的有效性。目前對于同時計及熱網(wǎng)與氣網(wǎng)動態(tài)特性的電-氣-熱互聯(lián)系統(tǒng)運行優(yōu)化的研究尚不多見,大部分側(cè)重于兩種能流的耦合。如文獻(xiàn)[68]提出了一種電-氣-熱互聯(lián)系統(tǒng)日前優(yōu)化調(diào)度模型,并將模型轉(zhuǎn)化為線性問題繼而利用GAMS軟件進(jìn)行求解,最后通過算例仿真分析證明了計及網(wǎng)絡(luò)動態(tài)特性在IES中的可行性與必要性。

    針對風(fēng)電接入背景下的不確定因素及其相關(guān)性,文獻(xiàn)[69]建立了電-氣-熱IES概率最優(yōu)能量流模型,并采用基于Nataf變換的三點估計法求解模型,以此驗證了模型的可行性與有效性。文獻(xiàn)[70]則針對風(fēng)電不確定性提出了一種IES擴(kuò)展規(guī)劃方法,并建立以投資、運行、電能不足與棄風(fēng)成本之和最小為目標(biāo)的混合整數(shù)線性規(guī)劃模型,進(jìn)而采用基于場景法的兩階段規(guī)劃求解策略對模型進(jìn)行求解分析,再利用GAMS對所建電-氣-熱IES進(jìn)行優(yōu)化求解,所得最優(yōu)擴(kuò)展規(guī)劃結(jié)果驗證了所提方法的正確性。

    綜合考慮風(fēng)電不確定性、不同能源網(wǎng)絡(luò)的動態(tài)特性以及經(jīng)濟(jì)、環(huán)境等效益的電-氣-熱IES優(yōu)化目前鮮有研究,這將是未來值得研究的方向。

    3.2 數(shù)據(jù)中心IES優(yōu)化運行國內(nèi)外研究現(xiàn)狀

    目前一般采用水冷、風(fēng)冷、液冷、蒸發(fā)冷卻與機(jī)械制冷聯(lián)合或冰蓄冷空調(diào)制冷等制冷方式來降低數(shù)據(jù)中心內(nèi)部過熱問題,但傳統(tǒng)的電制冷或空調(diào)制冷等方式嚴(yán)重依賴電網(wǎng),導(dǎo)致制冷能耗增加,系統(tǒng)運行效率低下,同時當(dāng)發(fā)生大面積斷電事故時,將嚴(yán)重威脅數(shù)據(jù)中心的安全和可靠運行[71]。

    目前鮮有對數(shù)據(jù)中心IES優(yōu)化運行的研究,但已有采用獨立的單分布式能源運用于數(shù)據(jù)中心,解決數(shù)據(jù)中心冷負(fù)荷和電負(fù)荷需求量大的問題。如文獻(xiàn)[72]從能源、經(jīng)濟(jì)和環(huán)境3個方面評估了冷熱電聯(lián)供的性能,表明數(shù)據(jù)中心可充分發(fā)揮冷熱電聯(lián)供的優(yōu)勢。文獻(xiàn)[73]采用分布式能源如冷熱電聯(lián)供系統(tǒng)為數(shù)據(jù)中心制冷,根據(jù)冷熱電聯(lián)供系統(tǒng)中余熱回收制冷的方式擺脫對電網(wǎng)的依賴性并提高可靠性,同時分布式能源將綜合能源進(jìn)行梯級利用,減少了輸電線路的能源損耗。文獻(xiàn)[74]采用冰蓄冷空調(diào)為數(shù)據(jù)中心供冷,緩解了電網(wǎng)負(fù)荷負(fù)擔(dān)。文獻(xiàn)[75]采用基于FC的冷熱電聯(lián)供機(jī)組以定功率方式為數(shù)據(jù)中心供能,以此降低系統(tǒng)運行成本和碳排放,但該數(shù)據(jù)中心采用的是定負(fù)荷運行,與工程應(yīng)用有一定差距。因目前大多是采用單種類型的分布式電源應(yīng)用于數(shù)據(jù)中心或者是定負(fù)荷數(shù)據(jù)中心的運行,這就導(dǎo)致系統(tǒng)運行效率低下,制冷耗電量增加,運行靈活性差、偏差大。

    若采用多能源系統(tǒng)作為輸入側(cè),以協(xié)調(diào)各類分布式可再生能源、分布式冷/熱/電聯(lián)供系統(tǒng)與電網(wǎng)絡(luò)的運行特性,滿足數(shù)據(jù)中心用冷/熱/電的需求,降低系統(tǒng)制冷耗電量和運行成本,理論上可有效解決數(shù)據(jù)中心內(nèi)部的局部過熱問題。

    如圖6所示,本文提出采用多能源輸入-多能源耦合-多能源輸出的運行方式來滿足數(shù)據(jù)中心的冷、電負(fù)荷需求,其中電負(fù)荷需求由風(fēng)力發(fā)電、光伏發(fā)電、微型燃?xì)廨啓C(jī)、燃料電池、儲能設(shè)備及外界電網(wǎng)來滿足,冷負(fù)荷需求由戶用空調(diào)、電制冷機(jī)、冰蓄冷裝置及熱電聯(lián)產(chǎn)機(jī)構(gòu)中的余熱吸收式制冷設(shè)備來滿足。其中電解槽電解反應(yīng)產(chǎn)生的H2一部分供給燃料電池,燃料電池可進(jìn)行發(fā)電,發(fā)電電量可繼續(xù)參與電解槽的電解;另一部分H2經(jīng)過甲烷反應(yīng)器與CO2反應(yīng)生成天然氣供微型燃?xì)廨啓C(jī)使用,微型燃?xì)廨啓C(jī)再進(jìn)行發(fā)電和制冷制熱,發(fā)電電量可繼續(xù)參與電解槽電解,從而形成“兩回環(huán)”結(jié)構(gòu)。采用上述研究框架原則上可提高電源利用率,因此對于數(shù)據(jù)中心IES調(diào)度優(yōu)化運行的研究顯得尤為重要,這將是未來需要進(jìn)一步研究的問題。

    4 面向數(shù)據(jù)中心IES優(yōu)化運行的研究展望

    4.1 數(shù)據(jù)中心的能源供給側(cè)不確定性及冷、電負(fù)荷預(yù)測研究

    IES是由冷、熱、電、氣等多種能流相互耦合的復(fù)雜多變系統(tǒng),能源供給側(cè)的可再生能源(如風(fēng)、光等)出力受天氣、區(qū)域位置等外界條件的影響有很大的間歇性和不確定性;數(shù)據(jù)中心中的冷、電負(fù)荷性質(zhì)受到機(jī)房布局結(jié)構(gòu)、建筑設(shè)計性質(zhì)、機(jī)房溫度、濕度、照明、建筑材料、通風(fēng)條件、機(jī)器設(shè)備臺數(shù)、用戶行為等因素的影響,造成數(shù)據(jù)中心負(fù)荷建模及求解的困難。因此,可再生能源出力和負(fù)荷需求的預(yù)測精度差。不僅如此,數(shù)據(jù)中心負(fù)荷與供給側(cè)能源的出力具有一定相關(guān)性,因此未來在研究數(shù)據(jù)中心IES優(yōu)化運行過程中,有必要同時研究可再生能源出力和負(fù)荷需求的不確定性以及數(shù)據(jù)中心的負(fù)荷動態(tài)變化特性,以提高數(shù)據(jù)預(yù)測精度。

    目前對于風(fēng)電等可再生能源出力的不確定性問題研究通常采用概率潮流和概率最優(yōu)潮流方法及其改進(jìn)方法;同時通常采用的負(fù)荷預(yù)測建模方法為宏觀類預(yù)測方法(如統(tǒng)計學(xué)方法)或微觀類預(yù)測方法(如物理-統(tǒng)計混合模型)[76-77],但此類方法因影響因素眾多而導(dǎo)致機(jī)理復(fù)雜,無法準(zhǔn)確建立負(fù)荷模型,因此有必要研究新興的深度學(xué)習(xí)、機(jī)器學(xué)習(xí)等人工智能方法對數(shù)據(jù)中心進(jìn)行負(fù)荷預(yù)測,以在不需要人為干預(yù)的情況下,依據(jù)歷史數(shù)據(jù)進(jìn)行自我學(xué)習(xí),從而進(jìn)行有效的負(fù)荷預(yù)測。

    圖6 多能源輸入-多能源耦合-多能源輸出的數(shù)據(jù)中心IES結(jié)構(gòu)Fig.6 Structure of data center IES with multiple energy input, coupling and output

    4.2 IES優(yōu)化運行中多能源動態(tài)模型的構(gòu)建研究

    目前對于IES的建模大多針對特定元件進(jìn)行穩(wěn)態(tài)研究,而對包含電、氣、冷、熱等元素的時間尺度等問題考慮較少,雖然部分學(xué)者已開始對系統(tǒng)進(jìn)行動態(tài)建模,但大都是對有關(guān)CCHP的系統(tǒng)進(jìn)行研究,而對于冰蓄冷裝置、熱泵等一些多時間尺度的耦合型設(shè)備動態(tài)特性和優(yōu)化特性等缺乏深入研究,同時現(xiàn)有的物理模型無法考察含電、氣、冷、熱的多能源系統(tǒng)的動態(tài)特性和多時間尺度特性,導(dǎo)致在進(jìn)行有關(guān)IES的規(guī)劃運行等方面遇到眾多難題。因此為提高研究精度,有必要構(gòu)建有關(guān)熱泵、冰蓄冷裝置等耦合設(shè)備的動態(tài)模型,并開發(fā)相應(yīng)的仿真工具以推進(jìn)有關(guān)IES研究的發(fā)展。

    4.3 采用人工智能方法的數(shù)據(jù)中心IES優(yōu)化調(diào)度策略研究

    目前存在的IES優(yōu)化運行策略因受可再生能源波動、通信延時等擾動問題的影響難以實現(xiàn)實時的自動調(diào)度運行,因此需要研究如何根據(jù)IES中的多能源輸入側(cè)數(shù)據(jù)、多能源耦合模型、數(shù)據(jù)中心冷電負(fù)荷輸出側(cè)數(shù)據(jù),并利用人工智能方法自動生成優(yōu)化運行方案,進(jìn)一步提取相應(yīng)知識和規(guī)則,形成IES優(yōu)化策略庫,以此實現(xiàn)IES在不同時段的實時自動調(diào)度策略,這將是未來的重要研究方向。

    4.4 數(shù)據(jù)中心IES優(yōu)化運行效益評價方法研究

    在運用多能源系統(tǒng)調(diào)配控制策略解決數(shù)據(jù)中心的制冷耗電問題時,如何評價系統(tǒng)的經(jīng)濟(jì)、環(huán)境、社會、投資、能耗效率等效益顯得極其重要,但目前對于IES的研究還未建立起適用于系統(tǒng)的綜合效益評估理論體系,同時一些效益評價方法如:神經(jīng)網(wǎng)絡(luò)算法、概率分析法、蒙特卡洛模擬法等僅限于理論方面的研究,適用于實際落地工程的評價指標(biāo)和方法均未完全成熟地建立。因此進(jìn)行適合中國IES發(fā)展的效益評價及其研究方法的探索將是值得研究的方向。

    5 結(jié) 語

    在能源危機(jī)和環(huán)境污染問題日益突出以及我國數(shù)據(jù)中心規(guī)模日益增大、制冷耗電問題日益嚴(yán)重的背景下,數(shù)據(jù)中心綜合能源系統(tǒng)優(yōu)化運行的研究對于發(fā)展經(jīng)濟(jì)、保護(hù)環(huán)境、提高能源利用效率、降低數(shù)據(jù)中心的制冷耗電等意義重大。本文首先介紹了綜合能源系統(tǒng)以及基于數(shù)據(jù)中心的綜合能源系統(tǒng)典型物理框架;進(jìn)而,根據(jù)給定的物理框架,分析了系統(tǒng)中耦合型運行設(shè)備的物理建模以及含能源集線器的通用建模和擴(kuò)展建模;然后,對國內(nèi)外綜合能源系統(tǒng)及基于數(shù)據(jù)中心的綜合能源系統(tǒng)優(yōu)化運行進(jìn)行了綜述、總結(jié)和探討,提出了多能源輸入-多能源耦合-多能源輸出的數(shù)據(jù)中心運行方式;最后,本文對數(shù)據(jù)中心綜合能源系統(tǒng)優(yōu)化運行研究的重點和難點進(jìn)行了分析,展望了該方向的幾個潛在可研究內(nèi)容,以期為未來相關(guān)的項目落地、仿真研發(fā)等研究提供參考。

    猜你喜歡
    數(shù)據(jù)中心耦合建模
    酒泉云計算大數(shù)據(jù)中心
    非Lipschitz條件下超前帶跳倒向耦合隨機(jī)微分方程的Wong-Zakai逼近
    聯(lián)想等效,拓展建?!浴皫щ娦∏蛟诘刃鲋凶鰣A周運動”為例
    基于PSS/E的風(fēng)電場建模與動態(tài)分析
    電子制作(2018年17期)2018-09-28 01:56:44
    不對稱半橋變換器的建模與仿真
    民航綠色云數(shù)據(jù)中心PUE控制
    電子測試(2018年11期)2018-06-26 05:56:24
    基于“殼-固”耦合方法模擬焊接裝配
    大型鑄鍛件(2015年5期)2015-12-16 11:43:20
    基于云計算的交通運輸數(shù)據(jù)中心實現(xiàn)與應(yīng)用
    求解奇異攝動Volterra積分微分方程的LDG-CFEM耦合方法
    三元組輻射場的建模與仿真
    级片在线观看| 久久国产精品人妻蜜桃| 日本一本二区三区精品| 巨乳人妻的诱惑在线观看| 欧美中文日本在线观看视频| 不卡av一区二区三区| tocl精华| 一级a爱片免费观看的视频| 日韩欧美国产在线观看| 国产欧美日韩精品亚洲av| 欧美不卡视频在线免费观看 | 久热这里只有精品99| 中文资源天堂在线| 两个人免费观看高清视频| 亚洲午夜理论影院| 老鸭窝网址在线观看| av天堂在线播放| 成人国语在线视频| 国产不卡一卡二| 淫妇啪啪啪对白视频| 法律面前人人平等表现在哪些方面| 又黄又粗又硬又大视频| 欧美激情高清一区二区三区| 岛国在线观看网站| 成年女人毛片免费观看观看9| 丝袜美腿诱惑在线| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲三区欧美一区| 色尼玛亚洲综合影院| 午夜亚洲福利在线播放| 一进一出抽搐gif免费好疼| 丁香六月欧美| 色av中文字幕| 亚洲人成电影免费在线| 欧美成人性av电影在线观看| 美女高潮喷水抽搐中文字幕| 91字幕亚洲| 国产精品1区2区在线观看.| 免费一级毛片在线播放高清视频| 色av中文字幕| 午夜福利在线观看吧| 在线观看66精品国产| 亚洲免费av在线视频| 亚洲精品美女久久av网站| 欧美精品啪啪一区二区三区| 婷婷丁香在线五月| 黑丝袜美女国产一区| 日韩欧美国产一区二区入口| 99热这里只有精品一区 | 国产精品国产高清国产av| 成人av一区二区三区在线看| 久久久久久久午夜电影| 中文字幕精品免费在线观看视频| 一区二区三区精品91| 亚洲成人精品中文字幕电影| 九色国产91popny在线| 少妇裸体淫交视频免费看高清 | 亚洲九九香蕉| 国产精品野战在线观看| 久久香蕉国产精品| 久久精品91蜜桃| 日本在线视频免费播放| 少妇裸体淫交视频免费看高清 | 午夜免费激情av| 嫩草影视91久久| 黑人操中国人逼视频| 母亲3免费完整高清在线观看| 人人妻,人人澡人人爽秒播| 国产高清videossex| 国产一区二区激情短视频| 亚洲精品一卡2卡三卡4卡5卡| 老司机午夜十八禁免费视频| 岛国视频午夜一区免费看| 国产精品免费视频内射| 久久精品国产亚洲av高清一级| 又紧又爽又黄一区二区| 精品熟女少妇八av免费久了| 亚洲自拍偷在线| 国产高清videossex| av在线播放免费不卡| netflix在线观看网站| 日韩大码丰满熟妇| 中文字幕人成人乱码亚洲影| 免费在线观看亚洲国产| 久久久国产欧美日韩av| 欧美成人一区二区免费高清观看 | 一个人免费在线观看的高清视频| 777久久人妻少妇嫩草av网站| 欧美丝袜亚洲另类 | 这个男人来自地球电影免费观看| 亚洲美女黄片视频| 精品第一国产精品| 亚洲国产欧美一区二区综合| 久久精品aⅴ一区二区三区四区| 88av欧美| 99国产精品99久久久久| 亚洲精品国产一区二区精华液| 亚洲成国产人片在线观看| 别揉我奶头~嗯~啊~动态视频| a在线观看视频网站| 美女高潮到喷水免费观看| 久久久国产欧美日韩av| 国产亚洲精品久久久久久毛片| 国产熟女xx| 国产99久久九九免费精品| 日本一本二区三区精品| 丰满人妻熟妇乱又伦精品不卡| 欧美色欧美亚洲另类二区| 亚洲人成电影免费在线| 嫁个100分男人电影在线观看| 99久久无色码亚洲精品果冻| 亚洲人成77777在线视频| 亚洲国产欧洲综合997久久, | 99久久无色码亚洲精品果冻| 久久久久久大精品| 哪里可以看免费的av片| 美女高潮喷水抽搐中文字幕| 久久精品夜夜夜夜夜久久蜜豆 | svipshipincom国产片| 一进一出抽搐动态| 午夜精品在线福利| 免费搜索国产男女视频| 欧美性猛交╳xxx乱大交人| АⅤ资源中文在线天堂| 免费电影在线观看免费观看| 日本五十路高清| 熟女电影av网| 日韩国内少妇激情av| 高潮久久久久久久久久久不卡| 亚洲成国产人片在线观看| 午夜精品久久久久久毛片777| 国产精品久久视频播放| 俄罗斯特黄特色一大片| 男女床上黄色一级片免费看| 久久国产精品男人的天堂亚洲| 日韩精品青青久久久久久| 动漫黄色视频在线观看| 欧美性长视频在线观看| 成人精品一区二区免费| 两个人看的免费小视频| 国产精品 欧美亚洲| 老司机午夜十八禁免费视频| 麻豆成人午夜福利视频| 国产成人精品无人区| 天天添夜夜摸| 99在线人妻在线中文字幕| 久久久久久久久久黄片| 国产欧美日韩一区二区三| 日本 av在线| 最好的美女福利视频网| 色哟哟哟哟哟哟| 久久精品91蜜桃| 女同久久另类99精品国产91| 免费在线观看视频国产中文字幕亚洲| 亚洲 欧美一区二区三区| 高潮久久久久久久久久久不卡| 50天的宝宝边吃奶边哭怎么回事| 免费在线观看日本一区| 国产激情久久老熟女| 色av中文字幕| av有码第一页| 午夜福利高清视频| 黄色 视频免费看| 999精品在线视频| 一区二区三区激情视频| 黑人巨大精品欧美一区二区mp4| 国产精品98久久久久久宅男小说| 成人国产综合亚洲| 欧美不卡视频在线免费观看 | 一级毛片高清免费大全| 日韩欧美国产在线观看| 国产片内射在线| 成年免费大片在线观看| 97超级碰碰碰精品色视频在线观看| 日韩一卡2卡3卡4卡2021年| 可以在线观看的亚洲视频| 欧美日韩福利视频一区二区| 18禁裸乳无遮挡免费网站照片 | 久久久久国产一级毛片高清牌| 精品午夜福利视频在线观看一区| 色av中文字幕| 婷婷丁香在线五月| 午夜久久久在线观看| 亚洲第一电影网av| 黄片播放在线免费| 日日摸夜夜添夜夜添小说| 国产激情欧美一区二区| 在线国产一区二区在线| 国产伦在线观看视频一区| 国产精品一区二区精品视频观看| 国产精品亚洲一级av第二区| 动漫黄色视频在线观看| 一进一出抽搐gif免费好疼| 最好的美女福利视频网| 两个人免费观看高清视频| 岛国视频午夜一区免费看| 日本精品一区二区三区蜜桃| 在线永久观看黄色视频| 欧美午夜高清在线| 国产一卡二卡三卡精品| 欧美午夜高清在线| 国产精品免费一区二区三区在线| 亚洲激情在线av| 免费高清视频大片| 日韩视频一区二区在线观看| 国产成人精品无人区| a在线观看视频网站| 亚洲人成伊人成综合网2020| 欧美色欧美亚洲另类二区| 亚洲人成电影免费在线| 中文字幕精品亚洲无线码一区 | 久久这里只有精品19| 99精品久久久久人妻精品| 中文字幕人妻熟女乱码| 欧美色视频一区免费| 中文字幕人妻熟女乱码| 麻豆成人午夜福利视频| 日韩有码中文字幕| 欧美色欧美亚洲另类二区| 精品久久蜜臀av无| 亚洲成国产人片在线观看| 午夜成年电影在线免费观看| 精品第一国产精品| 一级毛片高清免费大全| 国产成+人综合+亚洲专区| 777久久人妻少妇嫩草av网站| 亚洲人成网站在线播放欧美日韩| 亚洲专区国产一区二区| 999精品在线视频| 此物有八面人人有两片| 亚洲一区二区三区色噜噜| 亚洲一区二区三区色噜噜| 国产精品影院久久| 久久人人精品亚洲av| 美女免费视频网站| 亚洲国产毛片av蜜桃av| 亚洲性夜色夜夜综合| 亚洲专区国产一区二区| 超碰成人久久| 一本一本综合久久| av超薄肉色丝袜交足视频| 免费看十八禁软件| 久久这里只有精品19| √禁漫天堂资源中文www| 麻豆一二三区av精品| 午夜激情av网站| 97碰自拍视频| 国产主播在线观看一区二区| 亚洲中文av在线| 欧美zozozo另类| 麻豆av在线久日| 亚洲成a人片在线一区二区| 国产成人啪精品午夜网站| 人人妻人人看人人澡| 身体一侧抽搐| 妹子高潮喷水视频| 99久久久亚洲精品蜜臀av| 久久久久亚洲av毛片大全| 亚洲国产毛片av蜜桃av| 国产亚洲精品第一综合不卡| 午夜a级毛片| 十八禁网站免费在线| 久久精品国产亚洲av高清一级| 久99久视频精品免费| 国产人伦9x9x在线观看| 成人精品一区二区免费| 一级a爱视频在线免费观看| 韩国精品一区二区三区| 岛国视频午夜一区免费看| 欧美zozozo另类| 国产成人av教育| 可以免费在线观看a视频的电影网站| √禁漫天堂资源中文www| 黄色a级毛片大全视频| 欧美日韩精品网址| 无遮挡黄片免费观看| 哪里可以看免费的av片| 日本免费a在线| 亚洲色图 男人天堂 中文字幕| 女警被强在线播放| 一本精品99久久精品77| netflix在线观看网站| 又紧又爽又黄一区二区| 搡老岳熟女国产| 18禁国产床啪视频网站| 俄罗斯特黄特色一大片| av福利片在线| 国内久久婷婷六月综合欲色啪| 好看av亚洲va欧美ⅴa在| 侵犯人妻中文字幕一二三四区| 国产精品久久久久久亚洲av鲁大| 亚洲第一青青草原| 国产一区二区三区在线臀色熟女| 一本久久中文字幕| 精品无人区乱码1区二区| 亚洲国产精品久久男人天堂| 99国产精品一区二区蜜桃av| 久久99热这里只有精品18| 美女高潮到喷水免费观看| 免费搜索国产男女视频| 亚洲人成电影免费在线| 丰满人妻熟妇乱又伦精品不卡| 久久精品亚洲精品国产色婷小说| 欧美一区二区精品小视频在线| 色综合亚洲欧美另类图片| 日本免费a在线| 亚洲人成伊人成综合网2020| 久久青草综合色| 国产精品亚洲一级av第二区| 一a级毛片在线观看| 亚洲精品中文字幕一二三四区| 亚洲五月天丁香| 国产亚洲欧美98| 亚洲国产欧美网| 日韩免费av在线播放| 在线免费观看的www视频| 在线观看午夜福利视频| 精品日产1卡2卡| 中文字幕精品免费在线观看视频| 日韩免费av在线播放| 亚洲精品美女久久久久99蜜臀| 俺也久久电影网| 色尼玛亚洲综合影院| 亚洲av五月六月丁香网| 男人舔女人下体高潮全视频| 国内毛片毛片毛片毛片毛片| 午夜精品久久久久久毛片777| 久久午夜综合久久蜜桃| 香蕉久久夜色| 午夜福利高清视频| 一二三四社区在线视频社区8| 久久九九热精品免费| 91大片在线观看| 久热爱精品视频在线9| 亚洲片人在线观看| 欧美成人性av电影在线观看| 成人亚洲精品av一区二区| 岛国在线观看网站| 一级a爱片免费观看的视频| 黄色视频,在线免费观看| 给我免费播放毛片高清在线观看| 欧美中文综合在线视频| 久久久久免费精品人妻一区二区 | 亚洲精品一区av在线观看| 天堂√8在线中文| 亚洲成av人片免费观看| 色尼玛亚洲综合影院| 亚洲一区中文字幕在线| 亚洲第一av免费看| 国产伦人伦偷精品视频| 无人区码免费观看不卡| 成人欧美大片| 好男人电影高清在线观看| 日韩欧美一区二区三区在线观看| 欧美黑人巨大hd| 欧美成人免费av一区二区三区| 日日夜夜操网爽| 91字幕亚洲| 亚洲精品粉嫩美女一区| 欧美zozozo另类| 免费在线观看亚洲国产| 欧美国产日韩亚洲一区| 熟妇人妻久久中文字幕3abv| 男人舔女人的私密视频| 欧美性猛交黑人性爽| 欧美日韩一级在线毛片| 精品一区二区三区av网在线观看| 中文亚洲av片在线观看爽| 国产精品九九99| 日本免费a在线| 美女扒开内裤让男人捅视频| 日韩精品中文字幕看吧| av在线播放免费不卡| 亚洲av熟女| 在线播放国产精品三级| 黑人巨大精品欧美一区二区mp4| 日本 av在线| 日韩高清综合在线| 日日夜夜操网爽| 91成人精品电影| 国产色视频综合| 国语自产精品视频在线第100页| 免费高清在线观看日韩| 夜夜躁狠狠躁天天躁| 久久亚洲精品不卡| 亚洲成人久久爱视频| 熟女电影av网| 国产精品av久久久久免费| 午夜福利在线观看吧| 国产单亲对白刺激| 色综合婷婷激情| 亚洲国产毛片av蜜桃av| 欧美激情高清一区二区三区| 久久久久国产一级毛片高清牌| 午夜激情av网站| 欧美激情 高清一区二区三区| 丝袜人妻中文字幕| 国产在线精品亚洲第一网站| 老熟妇仑乱视频hdxx| 国产一区二区在线av高清观看| 后天国语完整版免费观看| 精品卡一卡二卡四卡免费| 美女大奶头视频| 久久天躁狠狠躁夜夜2o2o| 午夜免费鲁丝| 男人舔女人下体高潮全视频| 久久香蕉激情| 久久 成人 亚洲| 国产精品久久久av美女十八| 免费女性裸体啪啪无遮挡网站| 人妻久久中文字幕网| 国产精品1区2区在线观看.| 18禁国产床啪视频网站| 在线观看免费午夜福利视频| 欧美日本视频| 美女免费视频网站| 这个男人来自地球电影免费观看| 正在播放国产对白刺激| 亚洲欧美精品综合一区二区三区| 亚洲片人在线观看| 操出白浆在线播放| 国产伦一二天堂av在线观看| 亚洲国产精品999在线| 人成视频在线观看免费观看| 亚洲男人的天堂狠狠| 国产精品亚洲美女久久久| 国产精品久久久久久精品电影 | 看免费av毛片| 国产三级在线视频| 日韩av在线大香蕉| 精品久久久久久久久久久久久 | 久久久久亚洲av毛片大全| 国产黄色小视频在线观看| 国产av一区在线观看免费| ponron亚洲| 国产黄色小视频在线观看| 高清在线国产一区| 18禁观看日本| 丝袜美腿诱惑在线| 丰满人妻熟妇乱又伦精品不卡| 欧美亚洲日本最大视频资源| 欧美激情极品国产一区二区三区| av超薄肉色丝袜交足视频| 国产亚洲欧美精品永久| 妹子高潮喷水视频| 日韩欧美免费精品| 50天的宝宝边吃奶边哭怎么回事| 午夜影院日韩av| 俺也久久电影网| 国产伦一二天堂av在线观看| 亚洲国产看品久久| 在线观看午夜福利视频| 色综合婷婷激情| 亚洲精品久久成人aⅴ小说| 黄片大片在线免费观看| 精品国产国语对白av| 日韩 欧美 亚洲 中文字幕| 99在线视频只有这里精品首页| 欧美绝顶高潮抽搐喷水| 日本精品一区二区三区蜜桃| 欧美黄色淫秽网站| 亚洲专区国产一区二区| 哪里可以看免费的av片| 少妇的丰满在线观看| 国产亚洲精品久久久久久毛片| 哪里可以看免费的av片| 成人国语在线视频| 精品福利观看| 久久久久久国产a免费观看| 国产亚洲欧美精品永久| 一区二区三区高清视频在线| 午夜亚洲福利在线播放| 啦啦啦观看免费观看视频高清| 婷婷精品国产亚洲av在线| 这个男人来自地球电影免费观看| 久久久久久久久中文| 国产私拍福利视频在线观看| 久久这里只有精品19| 国产精品免费视频内射| 美女高潮到喷水免费观看| 色av中文字幕| 午夜视频精品福利| 久久精品91无色码中文字幕| 韩国av一区二区三区四区| 久久天堂一区二区三区四区| 亚洲成a人片在线一区二区| 动漫黄色视频在线观看| 最近最新免费中文字幕在线| 中文资源天堂在线| 日本精品一区二区三区蜜桃| 精品久久蜜臀av无| 999久久久精品免费观看国产| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧美日韩高清在线视频| 日韩欧美一区视频在线观看| 久久久精品国产亚洲av高清涩受| 日本免费一区二区三区高清不卡| 最近最新中文字幕大全免费视频| 看黄色毛片网站| 久久天堂一区二区三区四区| 亚洲成a人片在线一区二区| 女同久久另类99精品国产91| 女生性感内裤真人,穿戴方法视频| 两个人免费观看高清视频| 久99久视频精品免费| 国产精品免费视频内射| 国产精品av久久久久免费| 国产成人av激情在线播放| 亚洲国产看品久久| 亚洲人成伊人成综合网2020| 一区二区三区激情视频| 日本a在线网址| 精品久久久久久久人妻蜜臀av| 久久天堂一区二区三区四区| www国产在线视频色| 波多野结衣高清无吗| 一本久久中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 桃色一区二区三区在线观看| 熟女电影av网| 欧美成人性av电影在线观看| 久久国产亚洲av麻豆专区| 午夜a级毛片| 成年免费大片在线观看| 国产成人欧美| 欧美黑人巨大hd| 中文资源天堂在线| 日韩欧美三级三区| 熟女少妇亚洲综合色aaa.| 97碰自拍视频| 成人手机av| 又黄又爽又免费观看的视频| 嫩草影院精品99| 亚洲 欧美一区二区三区| 久久天堂一区二区三区四区| 国产亚洲欧美98| 黄色 视频免费看| 在线观看午夜福利视频| 一本大道久久a久久精品| 免费av毛片视频| 香蕉国产在线看| 国产激情偷乱视频一区二区| 亚洲自拍偷在线| 亚洲一码二码三码区别大吗| 91大片在线观看| 一进一出抽搐gif免费好疼| 色综合站精品国产| 国内精品久久久久久久电影| 一本大道久久a久久精品| 成人免费观看视频高清| 久久人妻av系列| 久久婷婷人人爽人人干人人爱| 99久久综合精品五月天人人| 好男人电影高清在线观看| 亚洲色图 男人天堂 中文字幕| 韩国av一区二区三区四区| 日韩视频一区二区在线观看| 看片在线看免费视频| av有码第一页| 中文字幕av电影在线播放| 在线天堂中文资源库| 免费在线观看日本一区| 岛国视频午夜一区免费看| 欧洲精品卡2卡3卡4卡5卡区| 男女下面进入的视频免费午夜 | 人人妻人人看人人澡| 欧美日韩中文字幕国产精品一区二区三区| 成人一区二区视频在线观看| 黄色视频不卡| 中文字幕久久专区| 在线天堂中文资源库| 日日爽夜夜爽网站| 黄色视频不卡| 久久伊人香网站| 18禁观看日本| 午夜老司机福利片| 亚洲国产欧美日韩在线播放| cao死你这个sao货| 一本综合久久免费| 91av网站免费观看| 欧美乱色亚洲激情| 亚洲精品中文字幕一二三四区| 91在线观看av| 亚洲av电影在线进入| 亚洲av片天天在线观看| 亚洲性夜色夜夜综合| 国产伦在线观看视频一区| 国产伦人伦偷精品视频| 不卡av一区二区三区| 51午夜福利影视在线观看| 国产亚洲精品一区二区www| 日本免费a在线| 成人国产综合亚洲| 色精品久久人妻99蜜桃| 亚洲五月色婷婷综合| 一区二区三区精品91| 亚洲精品在线美女| 精品欧美国产一区二区三| 丝袜在线中文字幕| 99精品久久久久人妻精品| 伊人久久大香线蕉亚洲五| 97超级碰碰碰精品色视频在线观看| 亚洲欧美日韩无卡精品| 欧美日韩亚洲国产一区二区在线观看| 午夜a级毛片| 精品久久久久久久久久久久久 | 欧美成人性av电影在线观看| 亚洲午夜精品一区,二区,三区| 这个男人来自地球电影免费观看| 国产视频一区二区在线看| 黄片大片在线免费观看| 999精品在线视频| 黄色片一级片一级黄色片| 亚洲真实伦在线观看| 久久久久国产一级毛片高清牌|