兆瓦級堿性電解槽多變量數(shù)字孿生仿真系統(tǒng)構(gòu)建

摘 要：堿性電解水制氫系統(tǒng)耦合性強、運行工況復(fù)雜，屬于危險化工生產(chǎn)環(huán)境，因此難以開展可再生能源制氫寬功率波動測試等極端實驗。以張家口崇禮制氫廠的1 MW堿性電解槽為研究對象，通過生產(chǎn)實驗數(shù)據(jù)，分析溫度、壓力和堿液流量等對氣體純度和小室電壓等參數(shù)的影響。首先對運行機理進行分析，對數(shù)據(jù)進行多元非線性擬合出經(jīng)驗公式，建立兆瓦級堿性電解槽的電壓模型、法拉第效率模型、氧中氫濃度模型，并對氣液分離等后處理系統(tǒng)建模；其次對模型結(jié)果進行理論分析與驗證，通過模型分析得出制氫系統(tǒng)安全運行邊界；最后將制氫系統(tǒng)的數(shù)學(xué)模型與三維模型相結(jié)合，實現(xiàn)制氫系統(tǒng)在虛擬空間中的映射，為研究可再生能源制氫寬功率波動提供先驗平臺。

關(guān)鍵詞：堿性電解槽；數(shù)字孿生；制氫；寬功率波動；多變量；安全邊界

中圖分類號：TK911 文獻(xiàn)標(biāo)志碼：A

0 引 言

被視為“人類終極能源”的氫能具有熱值高、可循環(huán)、清潔無污染等優(yōu)點，“十四五”發(fā)展規(guī)劃進一步鞏固了氫能在中國未來能源體系中的地位［1］。中國早在2020年就明確提出雙碳愿景，而氫氣在加快能源轉(zhuǎn)型中發(fā)揮著重要作用［2-3］?！熬G氫”是通過風(fēng)、光等可再生能源給電解槽供電制得，因其成本低、零碳排放而備受世界各國學(xué)者的青睞［4-5］。目前具有容量大、成本低等特點的堿性溶液電解水制氫方式被廣泛采用［6］，但由于風(fēng)、光等具有隨機性和波動性［7-8］，直接作用于堿性電解槽之后會使系統(tǒng)安全參數(shù)超過臨界值，同時頻繁的啟停會使產(chǎn)氫效率大大降低，因此堿性電解槽目前還不能很好的適應(yīng)可再生能源的寬功率波動。為更好的研究適應(yīng)寬功率波動的控制策略，需構(gòu)建大功率堿性電解水制氫數(shù)字孿生系統(tǒng)，實現(xiàn)制氫系統(tǒng)的虛擬映射，以測試堿性電解槽的寬功率安全穩(wěn)定運行條件。

數(shù)字孿生在智能制造領(lǐng)域發(fā)展迅速［9-10］，可再生能源和電力系統(tǒng)也相繼推出相關(guān)數(shù)字孿生模型［11-12］。構(gòu)建堿性電解水制氫數(shù)字孿生系統(tǒng)需要探索小室電壓、氧中氫等參數(shù)的特性，而氧中氫作為影響系統(tǒng)安全的重要參數(shù)成為國內(nèi)外研究的熱點，文獻(xiàn)［13］通過對250 kW的電解槽進行實驗，發(fā)現(xiàn)系統(tǒng)壓力可影響氧中氫的含量，通過調(diào)節(jié)系統(tǒng)壓力可使電解槽在30%～100%的寬功率范圍下穩(wěn)定運行，但并未考慮其他工藝參數(shù)的影響；文獻(xiàn)［14］從安全性角度對電解槽展開寬功率實驗研究，研究表明適當(dāng)調(diào)節(jié)堿液流量可使氧中氫含量穩(wěn)定在安全范圍之內(nèi)，但未考慮系統(tǒng)壓力；文獻(xiàn)［15］從工藝條件展開實驗，通過改變溫度、堿液的流速和濃度觀察氣體純度的變化，但并未調(diào)節(jié)系統(tǒng)壓力。上述文獻(xiàn)從實驗的角度給出了氧中氫等特性的影響因素。文獻(xiàn)［16］基于熱力學(xué)、傳熱學(xué)和電化學(xué)理論建立了熱模型和堿性電解槽的數(shù)學(xué)模型，該模型適用性強，被廣泛運用于電解槽數(shù)學(xué)模型的研究，然而文章［16］僅考慮了溫度的影響，忽略了其他因素對電解槽的影響；文獻(xiàn)［17］在文獻(xiàn)［16］的基礎(chǔ)上增加了系統(tǒng)壓力的影響，建立了15 kW電解槽的半經(jīng)驗數(shù)學(xué)模型，但未考慮堿液流量，研究表明堿液流量對氧中氫含量有影響作用；文獻(xiàn)［18］在文獻(xiàn)［17］所建模型的基礎(chǔ)上利用Aspen對工藝進行了仿真，但并未用控制器進行順序過程控制；文獻(xiàn)［19］完全基于機理和物理基礎(chǔ)建立堿性電解槽的小室電壓模型，并未利用運行數(shù)據(jù)來減弱繁雜公式推導(dǎo)，也未建立其他模型；文獻(xiàn)［20］給出堿性電解槽的數(shù)字孿生框架，并詳細(xì)闡述了層級結(jié)構(gòu)；文獻(xiàn)［21］基于實驗與原理，以溫度為變量，建立堿性電解槽的數(shù)字孿生模型，但并未建立氧中氫濃度和法拉第效率模型。

綜合上述文獻(xiàn)，目前對堿性電解水制氫數(shù)字孿生系統(tǒng)的研究停留在千瓦級，而且未綜合考慮系統(tǒng)壓力和堿液流量的影響。本文探索溫度、系統(tǒng)壓力、堿液流量等對小室電壓、氧中氫等特性的影響，并考慮氣液分離等后處理系統(tǒng)，采用機理分析和數(shù)據(jù)驅(qū)動的方式建立兆瓦級堿性電解水制氫數(shù)字孿生系統(tǒng)，該系統(tǒng)利用用于過程控制的OLE（OLE for process centrol，OPC）等技術(shù)可與控制器交互，對數(shù)字孿生系統(tǒng)進行了測試驗證并得到系統(tǒng)的安全運行邊界，以期為研究可再生能源寬功率波動的控制策略提供實驗平臺和參考。

1 堿性電解水制氫數(shù)字孿生系統(tǒng)框架

堿性電解水制氫系統(tǒng)的數(shù)字孿生模型分為實體物理對象層、數(shù)據(jù)采集與控制實體層、數(shù)字孿生核心實體層和數(shù)字孿生模型驗證與應(yīng)用層。層與層之間相互聯(lián)系，逐層向電解水制氫系統(tǒng)數(shù)字孿生模型遞進。制氫系統(tǒng)的數(shù)字孿生模型框架如圖1所示。

1.1 制氫系統(tǒng)數(shù)字孿生模型層級結(jié)構(gòu)

實體物理對象層即為堿性電解水制氫系統(tǒng)的實體對象，本文研究的電解水制氫系統(tǒng)主要包括堿性電解槽、氣液分離和后處理系統(tǒng)，電解水制氫系統(tǒng)的工藝圖如圖2所示，電源供給電解槽直流電，堿性電解槽產(chǎn)生的氣體經(jīng)由堿液換熱器流向氣體分離器與洗滌器，經(jīng)分離器與洗滌器處理后走向下一工藝過程。由氣體帶出的堿液在堿液循環(huán)泵的作用下經(jīng)由堿液過濾器重新回到電解槽內(nèi)。數(shù)據(jù)采集與控制實體層即為電解水制氫系統(tǒng)中各種傳感器與控制流程。本文重點論述數(shù)字孿生核心實體層，本層是堿性電解水制氫系統(tǒng)數(shù)字孿生模型的核心。將所建數(shù)學(xué)模型在虛擬仿真平臺進行測試，利用可編程邏輯控制器（programmable logic controller，PLC）對模型進行控制與驗證，最后將運行數(shù)據(jù)實時傳送給三維模型，實現(xiàn)制氫系統(tǒng)數(shù)字孿生應(yīng)用層的構(gòu)建。

1.2 堿性電解槽及后處理系統(tǒng)的數(shù)學(xué)模型結(jié)構(gòu)

堿性電解槽的數(shù)學(xué)模型建立過程綜合機理建模與數(shù)據(jù)建模的方法，首先根據(jù)堿性電解槽的運行機理，利用電化學(xué)領(lǐng)域知識進行公式的推導(dǎo)與論證；再根據(jù)堿性電解槽的型號對機理模型的參數(shù)進行計算；為避免機理建模中繁雜的數(shù)學(xué)推導(dǎo)，利用數(shù)據(jù)建模的方式對機理模型進行簡化與優(yōu)化，形成經(jīng)驗公式，利用制氫系統(tǒng)的運行數(shù)據(jù)對經(jīng)驗公式進行多元非線性擬合，最終建立堿性電解槽的數(shù)學(xué)模型。此外在氣體純度的研究中，壓力和堿液流量是兩個重要的影響因素，因此將堿液流量考慮到堿性電解槽的數(shù)學(xué)模型中并繪制出安全邊界是必要的。

除了對堿性電解槽的數(shù)學(xué)模型進行重點構(gòu)建，本文還對后處理系統(tǒng)的數(shù)學(xué)特性進行模擬，后處理系統(tǒng)主要包括氣液處理系統(tǒng)和堿液循環(huán)系統(tǒng)。根據(jù)實驗數(shù)據(jù)，后處理系統(tǒng)的數(shù)學(xué)模型采用二階大慣性加時滯環(huán)節(jié)。為了能配合控制端進行順序流程控制，本文設(shè)置多種閥門特性并建立虛擬仿真平臺，后處理系統(tǒng)的參數(shù)可根據(jù)實驗要求自行設(shè)置。

2 堿性電解水制氫數(shù)字孿生系統(tǒng)構(gòu)建

2.1 堿性電解水制氫原理概述

純水中存在極少量的H+和OH-，因此將兩根惰性電極插入水中并通電會分別從陰陽兩極析出氫氣和氧氣，體積比為2∶1，但由于純水中的離子含量較少，常溫下水的離子積為[1×10-14]，H+和OH-濃度均為[1×10-7] mol/L，因此需加入強電解質(zhì)增強導(dǎo)電能力，加速電解過程。水電解制氫原理如圖3所示。

電解水制氫過程的電化學(xué)方程式如下：

陰極： [2H++2e-H2↑]

陽極： [4OH--4e-O2↑+2H2O]

總電解方程式：[ 2H2O通電 2H2↑+O2↑]

2.2 堿性電解水制氫系統(tǒng)數(shù)學(xué)模型構(gòu)建

2.2.1 電壓特性建模

電壓模型描述了堿性電解槽的電化學(xué)行為，由于本文實驗中電流為給定，根據(jù)歐姆定律，理論上即可得到電解槽的伏安特性曲線，但事實并非如此，由于堿性電解槽的強非線性，為了使電壓模型足夠準(zhǔn)確，必須考慮熱力學(xué)、動力學(xué)、流體力學(xué)和電阻效應(yīng)。無論是純水還是堿水，只有當(dāng)電壓達(dá)到一個限度才會開始制氫，稱為可逆電壓（[Urev]），它受溫度和壓力的影響，可用能斯特方程計算［22］。

最常用的電壓-電流模型由文獻(xiàn)［16］提出，該文基于熱力學(xué)和堿性電解槽的電化學(xué)特性，推導(dǎo)出堿性電解槽電壓模型的經(jīng)驗公式，即：

[Ucell=Urev+r1+r2TA I+s lnt1+t2/T+t3/T2AI+1] （1）

式中：[Ucell]——堿性電解槽電壓（小室電壓），V；[Urev]——可逆電壓，V；r1、r2、t1、t2、t3——待定系數(shù)；[T]——槽溫，℃；[I]——給定電流，A；[A]——電極的面積，cm2；[s]——特定的系數(shù)，與堿性電解槽的型號和功率有關(guān)。

文獻(xiàn)［17］基于15 kW堿性電解槽的實驗，提出電壓模型應(yīng)加入系統(tǒng)壓力的影響，因此將堿性電解槽的電壓公式優(yōu)化為：

[Ucell=Urev+r1+d1+r2T+d2pIs+s lnt1+t2T+t3T2Is+1]

（2）

式中：[d1]、[d2]——待定系數(shù)；p——系統(tǒng)壓力，MPa；[Is]——電流密度，A/cm2。

式（2）中加入了系統(tǒng)壓力的影響作用，但仍未考慮堿液流量對電壓的影響。流速的大小直接影響堿液的導(dǎo)電能力，從而影響堿性電解槽的電壓，因此堿性電解槽的電壓模型應(yīng)加入堿液流量的影響作用，如式（3）所示。

[Ucell=Urev+a+b+cIs+s lnt1+t2T+t3T2Is+1] （3）

其中：

[a=a1+a2T] （4）

[b=b1+b2p] （5）

[c=c1+c2v] （6）

式中：[v]——堿液流量，m3/h，隨著堿液流量的加入，堿性電解槽的電壓模型更加準(zhǔn)確；a、b、c、[a1]、[b1]、[c1]、a2、b2、c2、——待定系數(shù)。

最終堿性電解槽的電壓模型為：

[Ucell=Urev+a1+b1+c1+a2T+b2p+c2vIs+ " " " " " "s lnt1+t2T+t3T2Is+1] （7）

式（7）為每個小室的電壓，堿性電解槽的總電壓為小室數(shù)量乘以小室電壓即：

[Uelec=ncell·Ucell] （8）

式中：[Uelec]——電解槽總電壓，V；[ncell]——小室數(shù)量。

2.2.2 法拉第效率特性建模

在電化學(xué)領(lǐng)域，通常把實際生成物占理論生成物的百分比稱為法拉第效率。在電解水制氫系統(tǒng)中，單位時間內(nèi)的產(chǎn)氫量與堿性電解槽在單位時間內(nèi)消耗的電量緊密相關(guān)，法拉第效率被定義為單位時間內(nèi)系統(tǒng)實際生成氫氣與單位時間內(nèi)理論生成氫氣的物質(zhì)的量之比，可用經(jīng)驗公式來表示，文獻(xiàn)［16］推導(dǎo)出了法拉第效率的經(jīng)驗公式，即：

[ηF=I2sf1+Is2·f2] （9）

式中：[f1]、[f2]——待定系數(shù)。

理論上[ηF=1]，但在制氫系統(tǒng)中法拉第效率往往比1小，這是因為在電解過程中容易產(chǎn)生氣體雜質(zhì)，電解槽消耗的電量不能完全用于生產(chǎn)氫氣。而且由于堿性電解槽的密封性不好，在電解過程中有少量的空氣進入電解槽，因此與堿液發(fā)生如下反應(yīng)：CO2+2OH-→[CO2-3]+H2O。在實際產(chǎn)氫過程中，制氫系統(tǒng)的法拉第效率約為95%。

式（9）只考慮了法拉第效率與電流密度的關(guān)系，并未考慮溫度、壓力和堿液流量的影響。之前的研究證實壓力對法拉第效率的影響很小，但溫度直接影響了化學(xué)反應(yīng)的快慢，隨著溫度的升高，化學(xué)反應(yīng)的活化能降低，更利于反應(yīng)的進行，使電解水反應(yīng)平衡正向移動，加快氫氣產(chǎn)生的速率。堿液流動的速度直接決定反應(yīng)物的多少，如果堿液流速過慢則會導(dǎo)致堿液在分離器中堆積，堿性電解槽內(nèi)無足夠的堿液，從而降低產(chǎn)氫速率。

為了將溫度和流量的作用考慮在內(nèi)，對式（9）進行修訂，將[f1]、[f2]更新為以溫度和流量為變量的函數(shù)：

[f1=f11+ft1T+fv1v] （10）

[f2=f22+ft2T+fv2v] （11）

式中：[f11]、[f22]、[ft1]、[ft2]、[fv1]、[fv2]——待定系數(shù)。

利用法拉第效率可計算堿性電解槽實際產(chǎn)氫率，首先計算理論產(chǎn)氫率，根據(jù)電解水制氫的電化學(xué)方程式：

[H2O通電H2↑+12O2↑～2e_] （12）

堿水通電發(fā)生歧化反應(yīng)，生成1 mol氫氣需轉(zhuǎn)移2 mol電子，因此只需計算堿性電解槽單位時間內(nèi)消耗的電荷量，再根據(jù)法拉第系數(shù)計算出單位時間內(nèi)轉(zhuǎn)移的電荷的物質(zhì)的量，即可求得產(chǎn)氫速率，公式為：

[nH2=It2F·ncell] （13）

式中：[F]——法拉第常數(shù)，C/mol；[t]——時間，s；在標(biāo)準(zhǔn)大氣壓下，25 ℃時氣體摩爾體積為24.5 L /mol，即可求得單位時間內(nèi)產(chǎn)生氫氣的體積為24.5[nH2]。

2.2.3 氧中氫含量特性建模

制氫系統(tǒng)的氧中氫含量是一個重要參數(shù)，為使堿性電解槽在寬功率波動下長時間穩(wěn)定安全的運行，必須保證系統(tǒng)的氧中氫含量在安全范圍之內(nèi)，本文實驗取2%（50%LFL），如果高于此濃度則會存在爆炸的風(fēng)險。有很多原因可導(dǎo)致產(chǎn)品氣體含有其他氣體雜質(zhì)，特別是當(dāng)堿性電解槽以低功率運行時氣體純度將會大大下降，之前的研究表明壓力和溫度會影響氣體純度并給出了經(jīng)驗公式，但沒有經(jīng)驗公式將堿液流量考慮在內(nèi)，文獻(xiàn)［14］研究表明可通過調(diào)節(jié)堿液流量從而改變氧中氫的含量。由此可見，溫度、壓力、堿液流量均可影響氧中氫的含量。

氫中氧的含量一般在0.5%，氫氣的體積分?jǐn)?shù)遠(yuǎn)遠(yuǎn)高于爆炸極限，因此本文只針對氧中氫的濃度進行建模。根據(jù)前文的論述，氧中氫含量可根據(jù)溫度、壓力和堿液流量等從理論上進行建模。文獻(xiàn)［23］給出計算氧中氫濃度的經(jīng)驗公式：

[cHIO=C1+C2T+C3T2+C4+C5T+C6T2eC7+C8T+C9T2Is] （14）

式（14）準(zhǔn)確描述了在不同溫度下的氧中氫濃度，但此公式的變量較單一，并未考慮系統(tǒng)壓力的影響，在文獻(xiàn)［13］的實驗中已經(jīng)明確表明系統(tǒng)壓力會對氧中氫濃度產(chǎn)生影響，因此文獻(xiàn)［17］在式（14）的基礎(chǔ)上做了更新，如式（15）所示。

[cHIO=C1+E1+C2T+E2p+ " " " " "C3T2+C4+C5T+C6T2eC7+C8T+C9T2Is+ " " " " "E3p2+E4+E5p+E6p2eE7+E8p+E9p2Is] （15）

式（15）雖然加入了系統(tǒng)壓力的影響，但仍未考慮堿液流量的影響，堿液流量已經(jīng)被證實對氧中氫含量有影響，因此構(gòu)建的模型還是不夠精確，本文在式（15）的基礎(chǔ)上增加了堿液流量的影響，式（16）使氧中氫濃度的計算更加準(zhǔn)確，將堿性電解槽的安全運行邊界拓寬。

[cHIO=C1+C2T+E2p+F2v+C3T2+E3p2+F3v2+ " " " " "C4+C5T+C6T2eC7+C8T+C9T2Is+ E4+E5p+E6p2eE7+E8p+E9p2Is+ " " " " "F4+F5v+F6v2eF7+F8v+F9v2Is] （16）

式中：[C1]——常數(shù)；[C2～C9]——氧中氫濃度關(guān)于溫度的系數(shù)；[E2～E9]——氧中氫濃度關(guān)于壓力的系數(shù)；[F2～F9]——氧中氫關(guān)于堿液流量的系數(shù)。首先固定壓力和堿液流量，對方程進行多元非線性回歸得到[C2～C9]的值，然后在恒定的溫度和堿液流量下調(diào)整[F2～F9]的值，最后在固定的溫度和壓力下校正[F2～F9]的值，最終完成25個參數(shù)值的確定。

2.2.4 系統(tǒng)壓力特性建模

在制氫系統(tǒng)中，系統(tǒng)壓力是一個很重要的參數(shù)，壓力過高或過低都會對系統(tǒng)產(chǎn)生負(fù)面影響，而且前文中電壓和氧中氫的模型中都有壓力這個參數(shù)。因此，建立制氫系統(tǒng)壓力的數(shù)學(xué)模型十分必要。本文以氧分離器中的壓力為系統(tǒng)壓力，氫側(cè)壓力同理可推導(dǎo)。氣液分離器結(jié)構(gòu)示意如圖4所示，電解槽產(chǎn)生氫氣和氧氣分別經(jīng)堿液換熱器進入氫氧分離器，分離器中被氣體帶入的堿液在循環(huán)泵的作用下經(jīng)過堿液過濾器回到電解槽中，形成堿液循環(huán)。

氣液分離器、堿液換熱器與電解槽可近似看做連通器，通過調(diào)節(jié)閥調(diào)節(jié)氫側(cè)壓力使兩側(cè)液位相平?？蓪ρ醴蛛x器中的氣體列理想氣體狀態(tài)方程：

[pVO2=nO2RT] （17）

式中：[p]——系統(tǒng)壓力，Pa；[VO2]——氧分離器中氣體的體積，m3，理想情況下均為氧氣；[nO2]——生成的氧氣的物質(zhì)的量，mol；[R]——摩爾氣體常數(shù)，J/（mol·K）；[T]——溫度，K。

1） 當(dāng)[LO2lt;r]時，氧分離器透視圖如圖5所示。本文實驗中氧分離器為臥式分離器，形狀近似為圓柱形。

為了公式表示的統(tǒng)一性，由于arctan函數(shù)為奇函數(shù)，因此采用反正切函數(shù)來表示角度：

[θ=arctanr2-r-LO22r-LO2] （18）

式中：[LO2]——液面高度，mm；[h]——分離器長度，mm；[r]——側(cè)面半徑，mm。

由此算出扇形的面積為：

[S=12θ·r2] （19）

根據(jù)幾何關(guān)系算出堿液體積和氧氣的體積分別為：

[VKOH=[S-r·（r-LO2）sinθ]·hVO2=πr2h-VKOH] （20）

2） 當(dāng)[LO2≥r]時，分離器透視圖如圖6所示。

分析過程同1），這里直接給出扇形的面積和氧氣的體積：

[S=θ·r2VO2=S-r·LO2-rsinθ] （21）

將[L=LO2-r]表示為分離器中液位與分離器側(cè)面半徑[r]的差值，令[r2-L2=m]，這樣式（21）得到進一步簡化，氧氣分離器中氣體體積為：

[VO2=-r2arctanmL+La·h， "Llt;0πr2h-r2arctanmL-Lm·h ， "L≥0] （22）

將式（13）轉(zhuǎn)換成氧氣的物質(zhì)的量后代入式（17），再將式（22）帶入式（17），得到最終堿性電解水制氫系統(tǒng)的壓力模型為：

[p=14F·ncellRT·0tIdt-r2arctanmL+Lm·h， "Llt;014F·ncellRT·0tIdtπr2h-r2arctanmL-Lm·h， "L≥0] （23）

式（23）中的電荷量為電流對時間的積分，將液位高度與半徑[r]作差后代入式（23）即可得到系統(tǒng)壓力。

3 電解水制氫數(shù)字孿生系統(tǒng)仿真研究

本文根據(jù)張家口崇禮1 MW堿性電解槽的多日穩(wěn)定運行數(shù)據(jù)和通過在不同功率下調(diào)節(jié)壓力和堿液流量的運行數(shù)據(jù)，利用Matlab對第2節(jié)建立的數(shù)學(xué)模型進行多元非線性回歸，并對擬合結(jié)果進行闡述和分析。1 MW堿性電解槽的參數(shù)如表1所示，制氫系統(tǒng)的堿性電解槽如圖7a所示，單日實驗運行數(shù)據(jù)如圖7b所示。

圖8為實驗平臺的工藝流程，由于空氣中的CO2會與堿液發(fā)生反應(yīng)：CO2+2OH-→[CO2-3]+H2O，會削弱導(dǎo)電能力并產(chǎn)生雜質(zhì)堵塞管道，并且堿液長時間不使用會析出結(jié)晶，因此需在實驗結(jié)束后利用氮氣排空管道中的堿液并使氮氣充滿管道以起到保護的作用。

3.1 電壓模型仿真研究

根據(jù)式（7）進行電壓模型的仿真，經(jīng)過更新的電壓模型中有4個參數(shù)，因此分別固定其中的兩個參數(shù)，以另外兩個參數(shù)為變量繪制電壓模型的曲面圖。如圖9所示，本文分別以溫度-電流密度、壓力-電流密度、堿液流量-電流密度為變量分析影響堿性電解槽電壓的因素。

從圖9可看出，隨著電流的增大小室電壓增大，符合基本電學(xué)規(guī)律。如圖9a所示，取溫度為30～80 ℃，可看出隨著溫度的升高電解槽小室電壓降低。原因之一是溫度升高提高了反應(yīng)物的能量，使活化分子的體積分?jǐn)?shù)增大，因此需要更少的電能；二是因為電解水是吸熱反應(yīng)，焓變大于零，升高溫度會使化學(xué)平衡正向移動，因此也會導(dǎo)致所需的電能減少。如圖9b所示，由于本文的壓力變化范圍較寬，所以可看出電壓隨壓力變化的趨勢，一方面可逆電壓會隨壓力的增大而有微小的增大，但系統(tǒng)壓力的增大會使附著在電極表面以及溶解在堿液中的微小氣泡的體積減小，從而降低堿液的電阻，可見后者對電壓的影響較大，因此才會使小室電壓下降。如圖9c所示，過低的堿液流量會導(dǎo)致電解槽中缺少堿液從而減少壽命，因此將其下限定為10 m3/h，隨著堿液流量的增加，小室電壓呈下降趨勢，堿液流量的增加會及時帶走附著在電極上的氣泡，使電解液的電阻降低，因此會使小室電壓降低。

圖10為電解槽小室電壓在不同溫度、壓力和堿液流量下的擬合曲線，并將實測的數(shù)據(jù)點在圖中加以標(biāo)記。從圖10可看出，溫度對壓力的影響最為劇烈，這也是文獻(xiàn)［16］只將溫度作為影響小室電壓的唯一變量的原因。圖10結(jié)果與上文分析一致，并且擬合數(shù)值與實測數(shù)值的誤差在0.1 V之內(nèi)，因此可將其作為制氫系統(tǒng)的電壓模型，以供后續(xù)寬功率波動等試驗研究。

3.2 法拉第效率模型仿真研究

根據(jù)式（9）～式（11）進行法拉第模型的結(jié)果分析，分別以溫度-電流密度、堿液流量-電流密度為變量繪制堿性電解槽法拉第效率的模型，如圖11所示，從圖11可看出，法拉第效率隨電流密度的增大而升高，當(dāng)電流為額定電流時，最終的法拉第效率約為95%，與工程實際相吻合。如圖11a所示，取溫度為30～80 ℃，隨著溫度的升高法拉第效率略微降低，由于溫度升高會導(dǎo)致堿液中陰陽離子的運動加快，因此電解液的電阻降低，這會導(dǎo)致少量的電流不能用于產(chǎn)生氫氣，從而導(dǎo)致法拉第效率的降低；如圖11b所示，取堿液流量為10～19 m3/h，可看出堿液流量對法拉第效率的影響并不大。

圖12為不同參數(shù)下的法拉第效率。如圖12b所示，分別取14、16、18 m3/h這3個堿液流量值，觀察其對法拉第效

率的影響，對曲線進行局部放大之后可更明顯的看出，隨著堿液流量的增大法拉第效率略微下降，第2.2節(jié)中指出堿液流量的增大直接導(dǎo)致電阻的降低，因此堿液流量對效率的影響與溫度相似。另外，曲線的擬合結(jié)果與實測實驗數(shù)據(jù)的差值在0.1%之內(nèi)，可作為制氫系統(tǒng)的產(chǎn)氫效率模型。

3.3 氧中氫濃度模型仿真研究

根據(jù)本文建立的氧中氫濃度的數(shù)學(xué)模型，分別以溫度-電流密度、壓力-電流密度、堿液流量-電流密度為變量繪制氧中氫濃度的變化規(guī)律，如圖13所示，從圖13中可看出，當(dāng)給定電流小于額定電流的20%時，氧中氫的濃度呈指數(shù)上升，位于氫氣的爆炸極限范圍之內(nèi)，整個系統(tǒng)處于非常危險的狀態(tài)。因此，必須控制電流的給定值大于額定電流的20%，在本實驗中給定電流大于400 A。從圖13還可看出，溫度對氧中氫的影響不大，而壓力和堿液流量對氧中氫濃度的影響較大，壓力的增大會加快氣體分子的擴散速率，并且會增加氧氣的溶解度，最終呈現(xiàn)的結(jié)果是壓力增大會使氧中氫的濃度升高；如果堿液流量增加，回流到氫氧分離器中的堿液也會增加，而且氫氣的滲透率較高，會使堿液中溶有微量氫氣，從而導(dǎo)致氧中氫的濃度升高，這一點在文獻(xiàn)［24］中也有說明。

圖14為不同溫度、壓力和堿液流量下氧中氫濃度隨電流密度的變化曲線，可看出，隨著電流密度的增大，氧中氫濃度都在下降。如圖14a所示，隨著溫度的升高，氧中氫濃度略微升高，溫度升高會導(dǎo)致溶解在堿液中的氣體析出并且氣體分子的熱運動加快，從而擴散速率增加。圖14b和圖14c更加明顯的展示了壓力和堿液流量對氧中氫的影響，并且在圖中繪制了（50%LFL）和（75%LFL）氧中氫的安全范圍，曲線實驗數(shù)據(jù)和實際可測得的工程運行數(shù)據(jù)的偏差在0.1%內(nèi)，部分?jǐn)?shù)據(jù)由于工藝要求無法測得。

圖15為根據(jù)圖14的仿真結(jié)果繪制的不同給定電流下電解水制氫系統(tǒng)的安全運行邊界。由實驗結(jié)果和模型預(yù)測可看出溫度對氧中氫影響較小，而且由于換熱器的大慣性純遲延特性，溫度調(diào)節(jié)比較困難，因此本文繪制壓力和堿液流量的安全邊界。圖15a為系統(tǒng)壓力的安全運行邊界，圖15b為堿液流量的運行邊界。在圖15中將區(qū)域分為5部分，1區(qū)域代表氧中氫濃度在2%（50%LFL）以下，在此區(qū)域內(nèi)系統(tǒng)可長期穩(wěn)定運行；2區(qū)域代表氧中氫濃度在2%以上3%（75%LFL）以下，此區(qū)域內(nèi)系統(tǒng)也可穩(wěn)定運行，但如果工藝對氧中氫含量有嚴(yán)格的要求，那么制氫系統(tǒng)不應(yīng)長期運行在此區(qū)域；如圖15a所示，3區(qū)域代表氧中氫模型雖然在此區(qū)域內(nèi)是2%以下，但在實際工程中此區(qū)域較難測得，因為在接近額定電流時調(diào)節(jié)閥即使達(dá)到飽和也無法將系統(tǒng)壓力降的很低，而且也不利于工藝的進行；4區(qū)域是氧中氫濃度高于4%的區(qū)域，屬于極端危險區(qū)域，制氫系統(tǒng)進入此區(qū)域后應(yīng)立即停機，然而在數(shù)字孿生模型中可允許進入此區(qū)域以開展寬功率波動實驗；5區(qū)域是氧中氫濃度高于3%低于4%的區(qū)域，制氫系統(tǒng)在此區(qū)域內(nèi)不應(yīng)長時間運行，可短時進入。通過前文的氧中氫濃度模型繪制的安全邊界可在安全區(qū)域內(nèi)進行尋優(yōu)，以達(dá)到最大的運行效率，為寬功率波動的研究提供了準(zhǔn)確有效的實驗平臺。

3.4 制氫系統(tǒng)虛擬仿真平臺搭建與三維監(jiān)控

根據(jù)第2節(jié)建立的數(shù)學(xué)模型，本文利用易控天地軟件將數(shù)學(xué)模型加以集成，并加上閥門、分離器、換熱器和傳感器等特性，最終建立電解水制氫系統(tǒng)的虛擬仿真平臺。其中調(diào)節(jié)閥和傳感器均近似為一階慣性環(huán)節(jié)，時間常數(shù)可自行設(shè)置，另外調(diào)節(jié)閥有飽和和死區(qū)特性可供用戶選擇，仿真平臺可根據(jù)需求設(shè)置仿真步長，由于在實驗過程中發(fā)現(xiàn)換熱器遲延比較嚴(yán)重，因此將其近似為二階慣性加純遲延環(huán)節(jié)，傳遞函數(shù)為：

[G（S）=KT1s+1T2s+1e-τs] （24）

式（24）中[K]、[T1]、[T2]和[τ]可根據(jù)實驗需求自行設(shè)置，也可根據(jù)實際工程的管路和工藝進行數(shù)學(xué)推理計算。

利用SketchuUp等3D建模軟件根據(jù)崇禮制氫現(xiàn)場建立電解水制氫系統(tǒng)的三維模型，包括堿性電解槽、氣液處理與洗滌系統(tǒng)、堿液循環(huán)系統(tǒng)、換熱系統(tǒng)、補水系統(tǒng)和各種閥門，并對管道進行貼圖處理，便于后期設(shè)計液體和氣體的流動效果，盡最大可能的還原崇禮電解水制氫系統(tǒng)的物理實體。

本文首先將堿性電解槽的數(shù)學(xué)模型和閥門傳感器等器件的特性集成到易控天地軟件，進行制氫畫面和動畫的設(shè)計，將系統(tǒng)壓力等模擬量與建立的數(shù)學(xué)模型和器件特性進行連接，形成虛擬仿真平臺；然后進行順序控制和過程控制的程序設(shè)計，控制端對虛擬仿真平臺進行實時控制，輸出對應(yīng)的閥門開度等信號，虛擬仿真平臺按照器件特性和數(shù)學(xué)模型輸出相應(yīng)的信號返還給控制端，實現(xiàn)控制端與虛擬仿真平臺的實時交互。虛擬仿真平臺以MySQL為后臺數(shù)據(jù)庫，將仿真結(jié)果實時傳入數(shù)據(jù)庫，如圖16所示，前端界面利用WebGL、js等技術(shù)進行設(shè)計，將建立的三維模型導(dǎo)入進來，三維監(jiān)控系統(tǒng)實時讀取數(shù)據(jù)庫中的內(nèi)容，將閥門狀態(tài)、模擬量的數(shù)值以可視化形式進行輸出，并且管道中的流速會根據(jù)實際數(shù)學(xué)模型的輸出結(jié)果而變化。三維監(jiān)控系統(tǒng)實時顯示虛擬仿真平臺的仿真結(jié)果，最終完成了1 MW堿性電解水制氫數(shù)字孿生仿真系統(tǒng)的建立。

4 結(jié) 論

1）本文基于機理建模和數(shù)據(jù)建模的方式構(gòu)建了1 MW堿性電解水制氫系統(tǒng)的數(shù)學(xué)模型，包括電壓模型、法拉第效率模型、氣體純度模型、系統(tǒng)壓力模型以及換熱器、分離器和閥門等特性。模型的構(gòu)建以數(shù)據(jù)和數(shù)學(xué)公式為基礎(chǔ)，曲線的預(yù)測結(jié)果與工程實測結(jié)果的誤差在工程允許范圍之內(nèi)。

2）本文根據(jù)建立的數(shù)學(xué)模型繪制了1 MW堿性電解槽在典型參數(shù)（系統(tǒng)壓力和堿液流量）下的安全運行邊界，并對邊界進行了分類與說明，為后續(xù)進行寬功率波動實驗提供了參考。

3）本文根據(jù)數(shù)學(xué)模型和各種特性搭建虛擬仿真平臺，然后進行控制端程序的設(shè)計，其次利用3D制圖軟件，基于WebGL技術(shù)，建立堿性電解水制氫系統(tǒng)的三維模型，最后將三者有效的連接起來并測試運行特性，完成1 MW堿性電解水制氫數(shù)字孿生系統(tǒng)的搭建。

本文只針對制氫系統(tǒng)做了相應(yīng)的工作，而后續(xù)的純化干燥和儲氫系統(tǒng)涉及的機理有待研究，相應(yīng)的數(shù)字孿生模型亟待建立。

［參考文獻(xiàn)］

［1］ 程文姬， 趙磊， 郗航， 等. “十四五” 規(guī)劃下氫能政策與電解水制氫研究［J］. 熱力發(fā)電， 2022， 51（11）： 181-188.

CHENG W J， ZHAO L， XI H， et al. Research on hydrogen energy policy and water-electrolytic hydrogen under "the "14th "Five-Year "Plan［J］. "Thermal "power generation， 2022， 51（11）： 181-188.

［2］ 徐鋼， 薛小軍， 張鐘， 等. 一種基于電解水制氫及甲醇合成的碳中和能源技術(shù)路線［J］. 中國電機工程學(xué)報， 2023， 43（1）： 191-200.

XU G， XUE X J， ZHANG Z， et al. A new carbon neutral energy technology route based on electrolytic water to hydrogen and methanol synthesis［J］. Proceedings of the CSEE， 2023， 43（1）： 191-200.

［3］ ACHARYA A. Scaling-up green hydrogen development with " "effective " "policy " "interventions［J］. " "Journal " "of sustainable development， 2022， 15（5）： 135.

［4］ 李亮榮， 彭建， 付兵， 等. 碳中和愿景下綠色制氫技術(shù)發(fā)展趨勢及應(yīng)用前景分析［J］. 太陽能學(xué)報， 2022， 43（6）： 508-520.

LI L R， PENG J， FU B， et al. Development trend and application prospect of green hydrogen production technologies "under "carbon "neutrality "vision［J］. "Acta energiae solaris sinica， 2022， 43（6）： 508-520.

［5］ 陳夢萍， 任建興， 李芳芹. 風(fēng)光互補與電解水制氫系統(tǒng)負(fù)荷的協(xié)調(diào)穩(wěn)定運行［J］. 太陽能學(xué)報， 2023， 44（3）： 344-350.

CHEN M P， REN J X， LI F Q. Coordinated and stable operation of wind solar complementarity and load of electrolytic water hydrogen production system［J］. Acta energiae solaris sinica， 2023， 44（3）： 344-350.

［6］ 王培燦， 萬磊， 徐子昂， 等. 堿性膜電解水制氫技術(shù)現(xiàn)狀與展望［J］. 化工學(xué)報， 2021， 72（12）： 6161-6175.

WANG P C， WAN L， XU Z A， et al. Hydrogen production based-on anion exchange membrane water electrolysis： a critical review and perspective［J］. CIESC journal， 2021， 72（12）： 6161-6175.

［7］ ZHENG Y， YOU S， BINDNER H W， et al. Optimal day-ahead dispatch of an alkaline electrolyser system concerning thermal-electric properties and state-transitional dynamics［J］. Applied energy， 2022， 307： 118091.

［8］ HONG Z P， WEI Z X， Han X J. Optimization schedulingcontrol strategy of wind-hydrogen system considering hydrogen production efficiency［J］. Journal of energy storage， 2022， 47： 103609.

［9］ HAW J， SING S L， LIU Z H. Digital twins in design for additive manufacturing［J］. Materials today： proceedings， 2022， 70： 352-357.

［10］ HUANG W H， ZHANG Y J， ZENG W. Development and application of digital twin technology for integrated regional energy systems in smart cities［J］. Sustainable computing： informatics and systems， 2022， 36： 100781.

［11］ 臧釗. 基于BIM+GIS的京張高速鐵路空地一體 “數(shù)字孿生” 智能化運維技術(shù)研究［J］. 鐵道運輸與經(jīng)濟， 2022， 44（9）： 139-145.

ZANG Z. Operation and maintenance technology of air-ground integrated＼ “digital twins＼” based on BIM and GIS for intelligent Beijing-Zhangjiakou high speed railway［J］. Railway transport and economy， 2022， 44（9）： 139-145.

［12］ YANG D G， KARIMI H R， KAYNAK O， et al. Developments of digital twin technologies in industrial， smart city and healthcare sectors： a survey［J］. Complex engineering systems， 2021， 1（1）： 1-21.

［13］ ZHANG C， WANG J Y， REN Z B， et al. Wind-powered 250 kW electrolyzer for dynamic hydrogen production： a pilot study［J］. International journal of hydrogen energy， 2021， 46（70）： 34550-34564.

［14］ 寧楠. 水電解制氫裝置寬功率波動適應(yīng)性研究［J］. 艦船科學(xué)技術(shù)， 2017， 39（6）： 133-136.

NING N. Research on hydrogen generation system by water electrolysis under wide power fluctuation［J］. Ship science and technology， 2017， 39（6）： 133-136.

［15］ HAUG P， KOJ M， TUREK T. Influence of process conditions on gas purity in alkaline water electrolysis［J］. International journal of hydrogen energy， 2017， 42（15）： 9406-9418.

［16］ ULLEBERG ?. Modeling of advanced alkaline electrolyzers： a system simulation approach［J］. International journal of hydrogen energy， 2003， 28（1）： 21-33.

［17］ SáNCHEZ M， AMORES E， RODRíGUEZ L， et al. Semi-empirical model and experimental validation for the performance evaluation of a 15 kW alkaline water electrolyzer［J］. International journal of hydrogen energy， 2018， 43（45）： 20332-20345.

［18］ SáNCHEZ M， AMORES E， ABAD D， et al. Aspen Plus model of an alkaline electrolysis system for hydrogen production［J］. International journal of hydrogen energy， 2020， 45（7）： 3916-3929.

［19］ ABDIN Z， WEBB C J， GRAY E M. Modelling and simulation of an alkaline electrolyser cell［J］. Energy， 2017， 138： 316-331.

［20］ 江悅， 沈小軍. 堿性電解槽制氫設(shè)備數(shù)字孿生體構(gòu)建及應(yīng)用［J］. 高電壓技術(shù)， 2022， 48（5）： 1673-1683.

JIANG Y， SHEN X J. Construction and application of digital twin in hydrogen production system of alkaline water electrolyzer［J］. High voltage engineering， 2022， 48（5）： 1673-1683.

［21］ 江悅， 沈小軍， 呂洪， 等. 堿性電解槽運行特性數(shù)字孿生模型構(gòu)建及仿真［J］. 電工技術(shù)學(xué)報， 2022， 37（11）： 2897-2908.

JIANG Y， SHEN X J， LYU H， et al. Construction and simulation of operation digital twin model for alkaline water electrolyzer［J］. Transactions of China Electrotechnical Society， 2022， 37（11）： 2897-2908.

［22］ 劉承錫， 曾冠維， 廖敏芳， 等. 大容量電解槽動態(tài)仿真建模及其快速頻率響應(yīng)分析［J］. 電網(wǎng)技術(shù)， 2023， 47（11）： 4638-4646.

LIU C X， ZENG G W， LIAO M F， et al. Modeling of large-scale electrolyzers for real-time simulation and analysis of its fast frequency response［J］. Power system technology， 2023， 47（11）： 4638-4646.

［23］ HUG W， BUSSMANN H， BRINNER A. Intermittent operation and operation modeling of an alkaline electrolyzer［J］. International journal of hydrogen energy， 1993， 18（12）： 973-977.

［24］ 黃登高， 尹玉國， 胡石林， 等. 堿液流量對水電解槽運行的影響［J］. 河南化工， 2016， 33（6）： 26-28.

HUANG D G， YIN Y G， HU S L， et al. Influence of alkali flow on operation of water electrolysis cell［J］. Henan chemical industry， 2016， 33（6）： 26-28.

CONSTRUCTION OF MULTIVARIATE DIGITAL TWIN SIMULATION SYSTEM FOR MEGAWATT-SCALE ALKALINE ELECTROLYZER

Liang Tao1，Liu Zicong1，Tan Jianxin2，Jing Yanwei2，Lyu Liangnian3

（1. School of Artificial Intelligence， Hebei University of Technology， Tianjin 300401， China；

2. Hebei Jiantou New Energy Co.， Ltd，， Shijiazhuang 050051， China；

3. Goldwind Science amp; Technology Coltd， Beijing 102600， China）

Abstract：Alkaline electrolytic water hydrogen production system with strong coupling and complex operating conditions is a hazardous chemical production environment， so it is difficult to carry out extreme experiments such as wide power fluctuation test for renewable energy hydrogen production. In this paper， the 1 MW alkaline electrolyzer of Zhangjiakou Chongli hydrogen production plant is used as the research object， and the effects of variables such as temperature， pressure and alkali flow rate on parameters such as gas purity and cell voltage are analyzed by production experimental data. Firstly， the operating mechanism is analyzed， the data are fitted with multivariate nonlinearities to produce empirical equations， and the voltage model， Faraday efficiency model， and hydrogen concentration in oxygen model of the MW alkaline electrolyzer were established， and the post-treatment systems such as gas-liquid separation are modeled. Next， the model results are analyzed and verified theoretically， and the safe operation boundary of the hydrogen production system is derived from the model analysis. Finally， the mathematical model of the hydrogen production system is combined with the 3D model to realize the mapping of the hydrogen production system in the virtual space， which provides an a priori platform to study the wide power fluctuation of renewable energy hydrogen production.

Keywords：alkaline electrolyzer； digital twins； hydrogen production； wide power fluctuation； multivariate； security boundaries