車用燃料電池歐姆內(nèi)阻辨識算法

車用燃料電池歐姆內(nèi)阻辨識算法曾雷;符興鋒;任強;張松林;翟艷霞;王佳騏【摘要】針對車用燃料電池歐姆內(nèi)阻辨識問題，以質(zhì)子交換膜燃料電池工作原理為基礎(chǔ)，建立電池的等效電路模型利用電路相關(guān)知識建立等效電路的數(shù)學關(guān)系模型，在此基礎(chǔ)上推導出模型的差分方程.介紹目前參數(shù)辨識使用比較廣泛的遞推最小二乘法并結(jié)合實車NEDC工況估計出電池的歐姆內(nèi)阻.在Simulink仿真平臺上驗證了該算法的收斂性、有效性以及良好的魯棒性.與現(xiàn)有方法相比,該方法具有在線實時估計、無需額外激勵源以及硬件電路支持、算法較簡單等優(yōu)點.【期刊名稱】《汽車零部件》【年(卷)，期】2018(000)011【總頁數(shù)】5頁(P16-20)【關(guān)鍵詞】質(zhì)子交換膜燃料電池;歐姆內(nèi)阻;遞推最小二乘法【作者】曾雷;符興鋒;任強漲松林;翟艷霞;王佳騏【作者單位】廣州汽車集團股份有限公司汽車工程研究院廣州廣東511434;廣州汽車集團股份有限公司汽車工程研究院廣州廣東511434;廣州汽車集團股份有限公司汽車工程研究院廣州廣東511434;廣州汽車集團股份有限公司汽車工程研究院廣州廣東511434;廣州汽車集團股份有限公司汽車工程研究院廣州廣東511434;廣州汽車集團股份有限公司汽車工程研究院廣州廣東511434【正文語種】中文【中圖分類】TM9110引言燃料電池是一種把燃料所具有的吉布斯自由能部分轉(zhuǎn)換成電能的裝置，它分為堿性燃料電池、磷酸燃料電池、固體氧化物燃料電池、質(zhì)子交換膜燃料電池(ProtonExchangeMembraneFuelCell,PEMFC)。其中PEMFC以氫作為燃料，具有功率密度高、體積小、質(zhì)量輕、反應溫度適中、啟動速度快、比能量高等優(yōu)點，是21世紀最有前途的“綠色能源”裝置［1］，被認為是實現(xiàn)未來汽車工業(yè)可持續(xù)發(fā)展的重要方向之一［2］。對于PEMFC來說，歐姆內(nèi)阻是其運行的一個重要參數(shù)，它反映了燃料電池內(nèi)部的運行狀態(tài)［3］，是衡量質(zhì)子和電子在電池內(nèi)部傳輸難易程度的主要標志，體現(xiàn)發(fā)電效率的大?。?］。另外，通過估計歐姆內(nèi)阻可以間接地預測膜的濕度和電池壽命的衰減。文獻［5-6］中采用交流阻抗法測量燃料電池的內(nèi)阻，該方法對燃料電池的擾動比較小，但也需要額外的激勵源。文獻［7］中采用改進型交流阻抗法測量燃料電池的內(nèi)阻，但仍然需要額外的激勵源。文獻［3］中則采用雙卡爾曼濾波算法估計燃料電池歐姆內(nèi)阻，只需輸入燃料電池電壓電流數(shù)據(jù)即可實時辨識出燃料電池歐姆內(nèi)阻，但雙卡爾曼濾波算法相對較復雜。本文作者采用最小二乘法，根據(jù)燃料電池電壓電流數(shù)據(jù)實時辨識燃料電池歐姆內(nèi)阻。1燃料電池等效電路模型1.1燃料電池工作原理PEMFC由陽極、陰極和質(zhì)子交換膜組成，基本原理如圖1所示。工作時，氫氣和氧氣發(fā)生反應形成離子移動，從而在外部電路形成電流。陽極發(fā)生氧化反應，電子流流向外部電路，氫離子穿過質(zhì)子交換膜到達陰極：H2—2H++2e(1)陰極發(fā)生還原反應，氧氣與氫離子結(jié)合產(chǎn)生水：熱(2)電池總反應：電能+熱⑶圖1PEMFC基本原理1.2燃料電池工作原理文獻［8］中同時考慮燃料電池內(nèi)阻容性和阻性的特點，提出的電池內(nèi)部電路模型如圖2所示，該模型由一個電阻和一個阻容并聯(lián)環(huán)節(jié)組成。圖2燃料電池內(nèi)部電路模型模型中電阻R0即是所求的歐姆內(nèi)阻，由于該電阻的存在，當電流IL發(fā)生突變時也會引起輸出電壓UL的突變，由雙電層分布電容CTh和電阻RTh組成的RC并聯(lián)環(huán)節(jié)，體現(xiàn)了燃料電池電極表面多孔結(jié)構(gòu)的化學特性［3］。燃料電池能斯特電壓與氫氣壓力、氧氣壓力以及溫度等多種因素有關(guān)，實際得到的開路電壓要比理論計算得到的能斯特電壓?。?］。模型中的UOC可理解為燃料電池的開路電壓，在燃料電池正常工作區(qū)域(線性區(qū)域)可視為定值。1.3等效電路模型的差分方程根據(jù)電路相關(guān)知識可列出等效模型的數(shù)學關(guān)系如下：⑷下面推導該模型的差分方程，將公式(4)離散化得式中：&為采樣間隔；T為時間常數(shù)；k為采樣時刻。令Ek為k時刻的動態(tài)電壓:Ek=UL,k-UOC,k(6)由公式(5)和(6)可得Ek+1=-[e-At/TUTh/k+(1-e-At/T)RThIL/k]-IL/k+1R0=-[e-At/T(UOC,k-UL,k-IL,kR0)+(1-e-At/T)RThIL,k]-IL,k+1R0=e-At/TEk+(-R0)IL,k+1+[e-At/TR0-(1-e-At/T)RTh]IL,k⑺因此，差分方程可以寫成：Ek+1=aEk+bIL,k+1+cIL,k(8)式中：a、b、c為待估計參數(shù)，a、b、c分別為⑼由公式(9)可知，只要能辨識出差分方程的參數(shù)b，就可以求出電池等效模型的參數(shù)R0,即燃料電池的歐姆內(nèi)阻。2利用最小二乘法辨識歐姆內(nèi)阻2.1離線最小二乘法最小二乘法是一種古老的數(shù)據(jù)處理方法，由高斯首先提出最小二乘法的概念，并將它用于天文計算。在參數(shù)辨識領(lǐng)域中，最小二乘法是重要的估計方法。它既可以用于動態(tài)系統(tǒng)也可用于靜態(tài)系統(tǒng)，既可用于線性系統(tǒng)也可用于非線性系統(tǒng)，既可用于離線估計也可用于在線估計。在隨機的環(huán)境條件下，利用最小二乘法時，并不要求觀測數(shù)據(jù)提供出有關(guān)噪聲的概率統(tǒng)計特性，而這種方法所得出的結(jié)果卻有相當好的統(tǒng)計特性。下面對最小二乘法基本原理進行簡要介紹，假設被辨識系統(tǒng)的差分方程:z(k)=-a1z(k-1)-......-anz(k-na)+b0u(k)+b1u(k-1)+......+bnu(k-nb)+e(k)(10)式中：z(k)為系統(tǒng)輸出量的第k次觀測值，z(k-1)是系統(tǒng)輸出量的第k-1次觀測值，依次類推；u(k)為系統(tǒng)的第k次輸入值；e(k)為均值為零的噪聲；a1 an，b1 bn為待辨識參數(shù)。令：(11)式中：h為樣本集合；e為待辨識參數(shù)向量，則式(10)可改寫為z(k)=hT(k)0+e(k)(12)對于k=1,2,……，G，構(gòu)成了一個線性方程組，可以寫成：zG(k)=HG(k)0+eG(k)(13)式中：(14)HG=(15)最小二乘法的基本思想是找出一個參數(shù)估計值e,使得判定函數(shù)：(16)取極小值，此時的參數(shù)估計值稱為最小二乘法意義下的估計值，用eGs表示，可以推出：(17)2.2遞推最小二乘法及歐姆內(nèi)阻辨識上一節(jié)討論的最小二乘算法是一種離線算法，適合成批處理觀測數(shù)據(jù)。為了實時跟蹤動態(tài)系統(tǒng)的時變參數(shù)，需轉(zhuǎn)化成遞推算法。遞推算法的基本思想是當被辨識系統(tǒng)在運行時，每獲得一組新的觀測數(shù)據(jù)后，在上次估計結(jié)果的基礎(chǔ)上，利用新獲取的觀測數(shù)據(jù)對上次估計的結(jié)果根據(jù)遞推算法進行修正，從而遞推得到新的參數(shù)估計值。這樣隨著新的觀測數(shù)據(jù)的逐次引入，不斷地進行參數(shù)估計，直到參數(shù)估計值接近于真實值。遞推最小二乘算法的基本思想可以概括：此次新的估計值=上次老的估計值+修正項即新的估計值是在老的估計值的基礎(chǔ)上，利用新的觀測數(shù)據(jù)對老估計值進行修正而得到的。假設上次老的估計值：(18)則此次新的估計值：(19)令：(20)則可以推出：(21)整理得到最小二乘法的遞推公式:(22)最小二乘算法的推導中，沒有做任何近似處理。因此，遞推估計和離線估計在理論上是一致的。結(jié)合公式(8)電池模型的差分方程，可知：(23)式中：Ek=UL,k-UOC,k，開路電壓UOC,k在燃料電池正常工作區(qū)域可視為定值，給定初始值P0和665,0(后面章節(jié)有初始值設定分析)，再輸入燃料電池兩端電壓、電流值即可實時辨識出a、b、c3個參數(shù)。根據(jù)公式(9)可知歐姆內(nèi)阻為-b，從而得到電池的歐姆內(nèi)阻。3算法分析及驗證3.1歐姆內(nèi)阻辨識文中采用廣汽某款燃料電池車上的電堆，額定功率30kW,由120片單體串聯(lián)而成，工作電流0~500A,H作電壓60~120V,開路電壓105V。文中搭載NEDC(NewEuropeanDrivingCycle)!況進行算法的分析及驗證。NED。是新歐洲行駛工況，主要有歐洲、中國、澳大利亞等國家使用，由市區(qū)運轉(zhuǎn)循環(huán)和市郊運轉(zhuǎn)循環(huán)組成，其電壓電流曲線如圖3所示，采樣間隔為0.2s。圖3NEDC工況電壓電流曲線給定初始值如下：0GS,0=[1,1,1]Txm(24)式中：n為1x106;m為0。實時辨識出的歐姆內(nèi)阻曲線如圖4所示?？梢钥闯觯涸诒孀R初期結(jié)果波動比較大，這是因為算法在辨識初期還未收斂，之后曲線逐漸穩(wěn)定，驗證了算法的收斂性，最終的辨識結(jié)果為0.0379Q。圖4歐姆內(nèi)阻辨識結(jié)果3.2算法驗證燃料電池的可測量較少，常用的參數(shù)辨識準確與否的驗證方法是利用真實電壓輸出與辨識得到內(nèi)部參數(shù)值之后，利用辨識參數(shù)值導入仿真模型得到的電壓輸出值，兩者相比較可以得出參數(shù)辨識算法對電池內(nèi)部參數(shù)的辨識能力，兩者之間的誤差值越小，則辨識算法越準確，這是目前較為常用的參數(shù)辨識算法的評價方法。在Simulink仿真平臺上搭建好驗證模型，如圖5所示，輸入電流、真實的電壓以及開路電壓值，利用遞推最小二乘法估計出歐姆內(nèi)阻，然后反推出仿真電壓值，再結(jié)合真實電壓值得出電壓誤差曲線。電壓誤差如圖6所示，可以看出：整個辨識過程中電壓誤差值始終在±1.5V以內(nèi)，相對開路電壓誤差為1.4%，驗證了算法的有效性。3.3算法初始值的設定對于上述算法，在開始計算的第一步必須給定P0和0GS,0的初始值，表1是不同初始值對應的估計結(jié)果，m、n含義見上文所述。圖5驗證模型圖6電壓誤差結(jié)果表1不同初始值的估計結(jié)果Qmn101x1031x1061x101600.03790.03790.03790.0379-1000.03850.03790.03790.03791000.03740.03790.03790.0379從表中可以看出：當n的取值足夠大時(1000以上)，初值的設定對估計結(jié)果無影響。當n初始取值10時，m初始值-100和100對應的相對誤差也僅為1.6%和1.3%，在可接受范圍內(nèi)。因此，初始值的選定對估計結(jié)果影響很小，當n取值足夠大時，初值對估計結(jié)果的影響基本可以忽略。由此可見，算法具有良好的魯棒性。4結(jié)束語文中基于質(zhì)子交換膜燃料電池的等效電路模型，提出采用最小二乘法辨識電池的歐姆內(nèi)阻，所提出的方法與現(xiàn)有方法相比，具有在線實時估計、無需額外激勵源以及硬件電路支持、算法較簡單等優(yōu)點，而且通過搭載整車NEDC工況，在Simulink仿真平臺上驗證了該算法的收斂性、有效性以及良好的魯棒性。當然，該算法僅僅適用于燃料電池在線性工作區(qū)域內(nèi)，結(jié)合實車數(shù)據(jù)，可以發(fā)現(xiàn)燃料電池大部分時間都工作在線性區(qū)域，因此并不影響此方法的推廣。對于PEMFC來說，歐姆內(nèi)阻反映了燃料電池內(nèi)部的運行狀態(tài)，是衡量質(zhì)子和電子在電池內(nèi)部傳輸難易程度的主要標志，體現(xiàn)發(fā)電效率的大小。另外，通過估計歐姆內(nèi)阻可以間接地預測膜的濕度和電池壽命的衰減。由于電堆中各單體本身存在差異性以及電堆裝配工藝的原因，單體電池的輸出特性存在不一致，而單片電池的性能會直接影響到整個電堆的性能。因此，辨識出的歐姆內(nèi)阻可應用于：(1)間接估算膜的濕度，為加濕策略提供參考；(2)間接估算電堆的健康狀態(tài)，從而預測電堆壽命；(3)估算電堆輸出功率、發(fā)電效率等狀態(tài)參數(shù)；⑷辨識出單體的歐姆內(nèi)阻可判斷電堆的不一致性，從而及時采用相應措施?！鞠嚓P(guān)文獻】LARMINIEJ,DICKSA.FuelCellSystemsExplained,SecondEdition[M].Chichster:JohnWiley&SonsLtd,2010.劉宗巍，史天澤，郝瀚，等.中國燃料電池汽車發(fā)展問題研究[J].汽車技術(shù),2018(1):1-9.LIUZW,SHIT乙HAOH,etal.ResearchonMainProblemsAssociatedwithDevelopmentofFuelCellVehicleinChina[J].AutomobileTechnology,2018(1):1-9.魏學哲，陳金干，李佳.應用雙卡爾曼濾波算法估計燃料電池歐姆內(nèi)阻[J].電源技術(shù),2009,33(5):359-362.WEIXZ,CHENJG,LIJ.EstimatingOhmicResistanceofPEMFCBasedonDoubleKalmanFilters[J].ChineseJournalofPowerSources,2009,33(5):359-362.羅良慶.燃料電池內(nèi)阻在線測試軟件系統(tǒng)與健康狀態(tài)監(jiān)測研究[D].武漢：武漢理工大學，2010.戴濤,婁平,李小君.質(zhì)子交換膜燃料電池內(nèi)阻測試系統(tǒng)設計[J].儀表技術(shù),2011(8):20-22.DAIT,LOUP,LIXJ.InternalResistanceofProtonExchangeMembraneFuelCellsTestSystem[J].InstrumentationTechnology,2011(8):20-22.江竑旭，全書海，童亮.新型燃料電池內(nèi)阻測試儀的研究與設計[J].武漢理工大學學報(信息與管理工程版),2012,34(4):467-470.JIANGHX,QUANSH,TONGL.ResearchandDesignonaNew-styleFuelCellResistanceTeste