999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于多次滑動均值濾波的混合儲能功率分配與定容研究

2024-03-14 07:20:42田博文張志禹楊夢飛
電工技術(shù)學(xué)報 2024年5期
關(guān)鍵詞:指令分配模態(tài)

田博文 張志禹 楊夢飛

基于多次滑動均值濾波的混合儲能功率分配與定容研究

田博文 張志禹 楊夢飛

(西安理工大學(xué)電氣工程學(xué)院 西安 710048)

為了解決混合儲能系統(tǒng)功率分配時因模態(tài)混疊導(dǎo)致功率分配不精確、儲能系統(tǒng)成本過高的問題,提出一種多次滑動均值濾波(MMAF)的功率分配方法用于削弱模態(tài)混疊現(xiàn)象,降低混合儲能成本。首先,獲取滿足平抑要求的混合儲能最小總功率指令,采用MMAF算法對其進(jìn)行濾波,獲得蓄電池與超級電容各自的功率指令,引入Pearson相關(guān)系數(shù)量化模態(tài)混疊現(xiàn)象,作為判斷濾波次數(shù)和每次濾波滑窗大小的指標(biāo),將蓄電池和超級電容各自的功率指令作為儲能定容的可行域,考慮電池荷電狀態(tài)約束求取儲能適配的最小額定功率和額定容量;然后,基于等效運行時間建立蓄電池全壽命周期量化模型,為經(jīng)濟(jì)性分析提供依據(jù);最后,以改進(jìn)的混合儲能全壽命周期成本模型驗證了該文方法可以有效地限制模態(tài)混疊,降低混合儲能系統(tǒng)成本。

混合儲能 功率分配 模態(tài)混疊 多次滑動均值濾波(MMAF)算法 蓄電池全壽命周期量化模型

0 引言

高比例風(fēng)電接入電網(wǎng),由于其自身的波動性和隨機(jī)性會對電力系統(tǒng)造成沖擊[1]。儲能設(shè)備對電能具有時空轉(zhuǎn)移和靈活調(diào)度的能力,在風(fēng)電場側(cè)配置儲能設(shè)備是解決風(fēng)電波動的一個可行方案[2-3]。蓄電池作為能量型儲能的代表具有儲能能量密度大、存儲能力強(qiáng)等特性,但由于其存在功率密度小、使用壽命短等缺陷,不宜高頻率地切換充放電狀態(tài)[4-5]。超級電容作為功率型儲能,具有功率密度大、響應(yīng)時間短、使用壽命長等優(yōu)點[6-7]。將二者結(jié)合組成的混合儲能系統(tǒng)可以形成優(yōu)勢互補(bǔ)[8-9],克服單一儲能的劣勢,提高風(fēng)電平抑效率[10]。但是如何將混合儲能總功率精確地分配給蓄電池和超級電容卻是一個亟須解決的問題。混合儲能總功率分配不準(zhǔn)確引起混合儲能系統(tǒng)定容不合理從而導(dǎo)致混合儲能系統(tǒng)成本高的問題成為阻礙儲能行業(yè)發(fā)展的重要因素[11-13]。如何通過精準(zhǔn)分配混合儲能總功率來合理定容以實現(xiàn)經(jīng)濟(jì)效益最大化已成為近年來眾多學(xué)者研究的重大課題[14-16]。

國內(nèi)外研究人員針對上述問題開展了許多研究,得到了一些結(jié)論。文獻(xiàn)[17]采用經(jīng)驗?zāi)B(tài)分解(Empirical Mode Decomposition, EMD)對混合儲能總功率進(jìn)行分析,按照不同的截止頻率將總功率分配給蓄電池和超級電容。但EMD分解后各模態(tài)易發(fā)生混疊且分解過程中易存在噪聲,這對后續(xù)混合儲能的經(jīng)濟(jì)性分析十分不利。為了限制模態(tài)混疊,文獻(xiàn)[18-19]提出了基于變分模態(tài)分解(Variational Modal Decomposition, VMD)的頻域分析方法用于對混合儲能總功率進(jìn)行分解。雖然VMD在一定程度上限制了模態(tài)混疊,但其并不能在所有非平穩(wěn)信號上直接使用,且分解層數(shù)的選擇有較大的主觀性,因此文獻(xiàn)[20]提出了一種采用智能算法對VMD參數(shù)尋優(yōu)的混合儲能功率分配策略。這種方法在一定程度上克服了VMD參數(shù)選取的人為性,但尋優(yōu)算法無法保證優(yōu)化所得的參數(shù)是全局最優(yōu)解,進(jìn)而也無法判斷是否最大化地避免了模態(tài)混疊。文獻(xiàn)[21-23]則直接從儲能定容入手,建立風(fēng)光儲系統(tǒng),以經(jīng)濟(jì)成本為目標(biāo)探究了不同場景下混合儲能定容配置問題。文獻(xiàn)[24]基于合作博弈建立了共享儲能的運行機(jī)制,分析了不同合作模式下的收益,以此規(guī)劃儲能的容量配置。文獻(xiàn)[25]通過模型預(yù)測控制(Model Predictive Control, MPC)建立了混合儲能動態(tài)功率分配策略并獲得了較好的功率分配效果。文獻(xiàn)[26]利用動態(tài)小波對混合儲能系統(tǒng)功率進(jìn)行分配,克服了傳統(tǒng)小波不能根據(jù)時序自適應(yīng)地調(diào)整小波參數(shù)的局限性,但并沒有對動態(tài)小波分解后各模態(tài)是否發(fā)生混疊做進(jìn)一步分析。文獻(xiàn)[27]改進(jìn)了一階低通濾波器,在功率分配之前加入判斷指令以減少因模態(tài)混疊使混合儲能功率誤分配的情況。

盡管已經(jīng)有眾多研究人員針對此類問題進(jìn)行探究,但現(xiàn)階段仍存在如下問題:①在對混合儲能總功率分配時缺少可以量化模態(tài)混疊的指標(biāo);②由于風(fēng)電出力的不確定性,現(xiàn)有功率分配方法并不能對所有非平穩(wěn)信號都有較好的分解效果。

滑動均值濾波(Moving Average Filtering, MAF)因其原理簡單、適用性強(qiáng)的優(yōu)勢近年來被廣泛應(yīng)用于平滑風(fēng)功率波動中[28]。本文提出一種基于多次滑動均值濾波(Multiple Moving Average Filtering, MMAF)的混合儲能系統(tǒng)功率分配方法及定容策略。首先,獲取混合儲能總功率指令,通過MMAF對混合儲能總功率進(jìn)行分解,同時引入Pearson相關(guān)系數(shù)判斷濾波次數(shù)與滑窗大小,根據(jù)功率分配結(jié)果對混合儲能系統(tǒng)進(jìn)行定容;然后,基于等效運行時間建立了蓄電池全壽命周期量化模型,為混合儲能經(jīng)濟(jì)性分析提供理論依據(jù);最后,通過改進(jìn)的混合儲能全壽命周期成本(Life Cycle Cost, LCC)模型驗證了本文方法能夠更好地量化并限制模態(tài)混疊,降低混合儲能系統(tǒng)的綜合成本。本文的整體框架如圖1所示。

1 混合儲能系統(tǒng)總平抑功率指令的獲取

1.1 平抑風(fēng)電波動的風(fēng)儲發(fā)電系統(tǒng)模型

用于平滑風(fēng)功率波動的混合儲能系統(tǒng)結(jié)構(gòu)拓?fù)淙鐖D2所示,其功率流關(guān)系為

式中,Bat為蓄電池的有功出力;Cap為超級電容的有功出力;H為混合儲能的有功出力;W為風(fēng)電原始功率;G為風(fēng)電并網(wǎng)功率。

圖2 平滑風(fēng)功率的混合儲系統(tǒng)結(jié)構(gòu)拓?fù)?/p>

1.2 混合儲能系統(tǒng)總功率指令的獲取

為了平抑風(fēng)電并網(wǎng)功率波動,以混合儲能系統(tǒng)作為研究對象。在求解混合儲能總功率指令H時將混合儲能系統(tǒng)(Hybrid Energy Storage System, HESS)看作整體,不考慮蓄電池與超級電容各自的運行狀態(tài),通過定義式(3)~式(9)所示的多目標(biāo)函數(shù)獲取滿足平抑要求的HESS最小出力指令。

1)目標(biāo)Ⅰ:功率波動絕對值均值

2)目標(biāo)Ⅱ:儲能充放電總能量

式中,H為混合儲能在整個控制周期內(nèi)吞吐能量的總值。

3)目標(biāo)Ⅲ:儲能出力能力評價系數(shù)

式中,SOC()為時刻混合儲能的荷電狀態(tài)。

4)目標(biāo)Ⅳ: 儲能進(jìn)入死區(qū)時間

式中,D為控制周期內(nèi)儲能裝置由于荷電狀態(tài)(State of Charge, SOC)越限而無法吞吐能量的時間,即死區(qū)時間;SOCmax和SOCmin分別為儲能電池荷電狀態(tài)的上限和下限。

上述四個目標(biāo)函數(shù)均會影響HESS的運行狀態(tài)。通過多目標(biāo)哈里斯鷹算法(Multi Objective Harris Hawks Optimizer, MOHHO)[29]求取在四個目標(biāo)值最小時的HESS總功率指令H,為后續(xù)分析如何將HESS總功率指令分配給蓄電池與超級電容以實現(xiàn)混合儲能系統(tǒng)成本最低提供數(shù)據(jù)基礎(chǔ)。

2 HESS的功率分配策略

近些年常見的儲能功率分配方法有EMD、VMD等。這些方法普遍存在的問題是信號分解層數(shù)需要依靠經(jīng)驗判斷,主觀性較強(qiáng),分解后所得各模態(tài)的頻譜時常發(fā)生混疊,對功率分配及儲能設(shè)備的健康運行極其不利。針對上述問題,本文借鑒模態(tài)分解中“多層分解”的思想,采用MMAF對混合儲能總功率指令進(jìn)行多次濾波,同時引入Pearson相關(guān)系數(shù)作為判斷滑窗大小與濾波次數(shù)的一項重要指標(biāo),避免主觀地選擇濾波參數(shù),以達(dá)到最大程度抑制模態(tài)混疊的目的。

2.1 MMAF算法功率分配原理

MAF本質(zhì)屬于低通濾波,因其原理簡單,易于實現(xiàn)等優(yōu)點,被廣泛應(yīng)用于信號處理領(lǐng)域。將混合儲能總功率指令作為輸入信號,通過MMAF可以獲得相應(yīng)的低頻分量和高頻分量。本文重點對每次濾波后獲取的低頻分量進(jìn)行處理,最大限度地剔除混雜在低頻分量中的高頻分量。取最后一次濾波所得低頻分量作為蓄電池的功率指令,總功率指令與最終所得的低頻分量做差即可得到高頻分量,以此作為超級電容的功率指令。MMAF算法功率分配原理如圖3所示。

圖3 MMAF算法功率分配功率原理

2.2 MMAF算法參數(shù)的確定

2.2.1 滑窗大小的確定

2.2.2 濾波次數(shù)的確定

3 HESS定容策略

儲能定容包括對儲能電池額定功率與額定容量的配置。儲能定容方案的合理性及精確性直接影響其經(jīng)濟(jì)性。本文以獲取的蓄電池與超級電容各自分配的功率指令為基礎(chǔ),綜合考慮儲能電池荷電狀態(tài)的限制,分別對蓄電池與超級電容進(jìn)行定容規(guī)劃。定容策略框架如圖4所示。

圖4 混合儲能定容策略框圖

3.1 HESS額定功率的確定

同理超級電容的最小額定功率配置方式為

3.2 HESS額定容量的確定

儲能容量反映的是儲能設(shè)備吞吐電量多少的能力。儲能設(shè)備某一時刻剩余容量占額定容量的百分比用荷電狀態(tài)表征,有

基于式(19)可分別建立儲能電池充放電時SOC的動態(tài)遞推關(guān)系式。

在儲能電池正常運行的過程中,其額定容量的最小值應(yīng)滿足任意時刻的能量吞吐需求,同時還要滿足SOC的約束關(guān)系,即

將式(20)代入式(21)中可得

通過式(22)可求解出充電過程與放電過程中儲能額定容量的取值范圍,即

在儲能電池的充電(放電)過程中,其SOC逐漸上升(下降),如若要滿足式(21)的約束,則充、放電過程中儲能的額定容量均應(yīng)滿足不等式(23)中的最小值,即

在實際過程中,儲能電池的額定容量既要滿足充電過程的容量需求,又要滿足放電過程中的容量需求,因此取式(24)的最大值作為最終確定的儲能電池額定容量,即

基于上述分析和理論推導(dǎo)可知,蓄電池和超級電容的最小適配額定容量分別為

4 蓄電池全壽命周期量化模型與混合儲能LCC模型

4.1 蓄電池儲能全壽命周期量化模型

本文以鉛酸電池作為電化學(xué)儲能的代表,建立電池全壽命周期量化模型。由于超級電容的循環(huán)次數(shù)一般為幾十萬次到幾百萬次,遠(yuǎn)遠(yuǎn)超過蓄電池的循環(huán)次數(shù),因此本文將超級電容的使用壽命設(shè)為固定值。一般的蓄電池實際循環(huán)次數(shù)由式(27)的Gaussian壽命函數(shù)進(jìn)行擬合。

傳統(tǒng)的蓄電池循環(huán)壽命模型由放電深度的冪函數(shù)擬合,即

傳統(tǒng)模型以最大放電深度下對應(yīng)的循環(huán)次數(shù)代替蓄電池運行周期內(nèi)的全部循環(huán)次數(shù),模型的精確性有待提高[30]。為此,本文提出一種新的蓄電池儲能全壽命周期量化模型。由面向事件的蓄電池壽命模型(式(30))可知[31],電池實際放電量在使用中實時變化且不易準(zhǔn)確測量,而電池循環(huán)次數(shù)一般由生產(chǎn)廠家直接提供,相對較為準(zhǔn)確,因此可采用循環(huán)次數(shù)代替其中的電池放電量,得到在每一時刻實際放電深度下對應(yīng)的循環(huán)次數(shù)(式(31))。

折算后的蓄電池壽命模型不僅考慮了循環(huán)次數(shù)對電池壽命折損的累加效應(yīng),而且克服了傳統(tǒng)模型中用最大放電深度下的循環(huán)次數(shù)代替全壽命周期循環(huán)次數(shù)的局限性。

故在蓄電池的運行周期內(nèi),其實際循環(huán)壽命可通過對等效循環(huán)次數(shù)求和,以等效運行時間表示為

為了提高所建蓄電池儲能全壽命周期量化模型的魯棒性,現(xiàn)對模型做如下假設(shè):

(1)假設(shè)蓄電池儲能所處密閉空間內(nèi)的溫度和濕度均恒定。

(3)假設(shè)蓄電池儲能的充放電效率在運行過程中保持不變。

(4)假設(shè)蓄電池儲能不發(fā)生自放電現(xiàn)象。

4.2 混合儲能LCC模型

傳統(tǒng)的混合儲能LCC模型中電池維護(hù)成本通常為

在電池年維護(hù)成本中添加“電池維護(hù)成本增量”這一因子,用來描述每隔十年由于電池不斷老化的累積效應(yīng)所帶來維護(hù)成本的非線性增加現(xiàn)象,以便更加符合儲能電池在實際運行過程中的情況,提高LCC模型的可實用性。

改進(jìn)后的混合儲能LCC模型各部分表達(dá)式如下。

1)蓄電池年綜合成本

2)蓄電池年維護(hù)成本

3)超級電容年綜合成本

4)超級電容年維護(hù)成本

結(jié)合式(36)~式(39)可得混合儲能系統(tǒng)全壽命周期成本為

5 算例分析

本文以某100 MW級風(fēng)電場一天的連續(xù)風(fēng)電功率數(shù)據(jù)為實驗數(shù)據(jù),采樣時間間隔為1 min,共1 440個采樣點。獲得滿足在風(fēng)電平抑要求的混合儲能總功率指令序列如圖6a所示,平抑后風(fēng)電并網(wǎng)功率波動如圖6b所示,滿足我國GB規(guī)定的100 MW風(fēng)電場1 min有功功率波動最大值不得超過10 MW的限制。儲能系統(tǒng)相關(guān)參數(shù)見表1[32]。

圖6 混合儲能總出力與平抑效果

圖6a展示了不同算法對混合儲能總功率指令的求取結(jié)果。可以看出,MOHHO算法相較于MOPSO算法可以獲得更小的儲能出力,因此本文后續(xù)均以MOHHO算法獲得的混合儲能總功率指令進(jìn)行分析與討論。

圖6b展示了風(fēng)電并網(wǎng)功率波動情況。從中可以看出,未經(jīng)平抑的風(fēng)電并網(wǎng)功率部分時刻的波動值已經(jīng)超出10 MW,經(jīng)過HESS平抑后的風(fēng)電并網(wǎng)功率波動有了明顯減小,且1 min內(nèi)波動均未超過 10 MW,符合我國相應(yīng)風(fēng)電場裝機(jī)容量GB規(guī)定的1 min風(fēng)電有功功率并網(wǎng)波動要求。

表1 儲能相關(guān)參數(shù)

Tab.1 Relevant parameters of energy storage

5.1 最優(yōu)滑窗大小與濾波次數(shù)的確定

MMAF滑動窗口大小及濾波次數(shù)的不同對模態(tài)混疊程度均有直接影響,繼而蓄電池和超級電容所分配的功率指令就會有相應(yīng)的變化,影響儲能電池的使用壽命及最終混合儲能額定功率和容量的配置。

為了盡可能地避免模態(tài)混疊,本文通過引入Pearson相關(guān)系數(shù)作為判斷最優(yōu)滑窗大小與最優(yōu)濾波次數(shù)的一項指標(biāo)。在不斷增大滑窗大小時,取相關(guān)系數(shù)最小時對應(yīng)的滑動窗口作為當(dāng)前濾波的窗口大小,直至將遍歷完全。獲取了最優(yōu)滑窗大小之后,再對濾波次數(shù)進(jìn)行遍歷,求出每次濾波后的低頻分量與高頻分量的Pearson相關(guān)系數(shù)。

圖7 濾波次數(shù)和滑窗大小與不同類相關(guān)系數(shù)間的關(guān)系

圖8a展示了混合儲能總功率指令的頻譜圖。從中可以看出未,經(jīng)分解的HESS總功率頻譜較為雜亂,高低頻混雜在一起,這對混合儲能系統(tǒng)的使用壽命和經(jīng)濟(jì)性十分不利。

圖8 混合儲能系統(tǒng)功率指令頻譜關(guān)系

5.2 HESS的功率分配與定容

5.2.1 HESS總功率指令的分配結(jié)果

圖9 MMAF蓄電池與超級電容的功率指令

圖10展示了在本文所提MMAF方法下蓄電池在一個調(diào)度日內(nèi)的實際放電深度。

圖10 MMAF蓄電池放電深度變化

從圖10可以看出,蓄電池在整個調(diào)度日內(nèi)的最大放電深度僅為15.2%,所有時段的放電深度均在25%以下,處于淺循環(huán)放電狀態(tài)。電池的淺循環(huán)放電狀態(tài)可大大增加蓄電池的使用壽命,從而降低了因壽命問題而更換儲能電池所帶來的額外成本。

5.2.2 混合儲能的定容優(yōu)化結(jié)果

我國多個省市發(fā)布最新相關(guān)文件,明確新能源整體配儲標(biāo)準(zhǔn):每100MW風(fēng)電場按不低于其裝機(jī)容量的10%~20%配建儲能系統(tǒng)。但由于地理位置差異等因素,導(dǎo)致有些風(fēng)能資源匱乏的地區(qū)并不需要建設(shè)與標(biāo)準(zhǔn)等規(guī)模大小的儲能設(shè)備,否則會使得混合儲能系統(tǒng)整體的經(jīng)濟(jì)性較差。

對上述問題本文采取了針對性的混合儲能定容優(yōu)化方案。為了檢驗本文所提方案的有效性及建立蓄電池儲能全壽命周期量化模型的準(zhǔn)確性,分不同場景對比了不同儲能組合方案下儲能設(shè)備的使用壽命與儲能的LCC;對比了本文定容方法與標(biāo)準(zhǔn)配置方案下的經(jīng)濟(jì)性差異,結(jié)果見表2與表3。

表2 MMAF實驗結(jié)果

Tab.2 Experimental results of MMAF

表3 單一儲能實驗結(jié)果

Tab.3 Single energy storage experiment results

按照規(guī)定,本文算例中100MW風(fēng)電場應(yīng)配置10 ~20MW?h的儲能。表2的實驗結(jié)果得出滿足當(dāng)前混合儲能總功率指令要求的最小定容參數(shù),表明蓄電池和超級電容的額定功率和額定容量并不需要按標(biāo)準(zhǔn)配置即可滿足當(dāng)前的總平抑功率指令需求,避免了因儲能配置與實際情況不適配而導(dǎo)致的經(jīng)濟(jì)性差的現(xiàn)象。驗證了本文提出的混合儲能總功率指令分配方法與定容策略高度匹配。

為了更好地說明混合儲能系統(tǒng)在響應(yīng)功率指令時較單一儲能有更顯著的經(jīng)濟(jì)性,本文使蓄電池與超級電容分別各自響應(yīng)混合儲能總功率指令。在“單一儲能實驗”這組對比實驗中,儲能按照風(fēng)電場裝機(jī)容量的20%進(jìn)行定容。

通過表2和表3兩組實驗結(jié)果對比可知相較于單一儲能,由蓄電池和超級電容組成的混合儲能系統(tǒng)全壽命周期成本最優(yōu),LCC相比于單蓄電池和單超級電容分別降低了60%和90%;使用單一超級電容儲能進(jìn)行功率指令的響應(yīng)時經(jīng)濟(jì)性最差,這是由于超級電容的單位容量成本過高導(dǎo)致;使用單一蓄電池儲能進(jìn)行功率指令的響應(yīng)時經(jīng)濟(jì)性也較差,而且相比較于混合儲能系統(tǒng),蓄電池的等效運行時間從19.9降至6.5。這是由于在使用單一蓄電池儲能響應(yīng)功率指令時蓄電池充放電深度增加,頻繁地進(jìn)行充放電狀態(tài)切換,大大影響了其使用壽命,最終導(dǎo)致了成本的升高。

為了進(jìn)一步說明本文MMAF方法多次濾波對降低混合儲能經(jīng)濟(jì)性的顯著性,現(xiàn)對混合儲能總功率指令只進(jìn)行一次MAF分析混合儲能的LCC,以此驗證本文多次濾波的必要性,結(jié)果如圖11、圖12、表4所示。

結(jié)合圖11和表4的實驗數(shù)據(jù)可以看出,只經(jīng)過一次MAF獲取的蓄電池和超級電容功率指令由于存在模態(tài)混疊,使得原本應(yīng)分配給超級電容的功率指令分配給了蓄電池。如圖11中480 min和1 000 min時,蓄電池分得的功率指令明顯高于圖9中對應(yīng)時刻的功率指令,而圖11中超級電容此刻對應(yīng)的功率指令要比圖9中小很多。與表2的實驗數(shù)據(jù)相比,由于存在模態(tài)混疊現(xiàn)象,蓄電池響應(yīng)了部分高頻分量指令從而影響了其使用壽命,致使其等效運行時間從19.9減少至19.2。圖11展示了單次MAF蓄電池放電深度變化情況。需要說明的是由于混合儲能總功率指令在1 440min的信號不含高頻分量的原因才導(dǎo)致了通過MAF方法與MMAF方法得到的蓄電池最大放電深度相同。另外通過圖10與圖12對比分析可以看出,其余采樣點時單次MAF時蓄電池的放電深度普遍要高于MMAF時蓄電池的放電深度,這也是導(dǎo)致單次MAF蓄電池壽命減少的一個原因。此外由于MMAF模態(tài)混疊程度較MAF時有了大幅削弱,因此LCC也從566.2萬元減少至565.2萬元,驗證了本文MMAF方法降低混合儲能成本的必要性。

圖11 單次MAF蓄電池與超級電容的功率指令

圖12 單次MAF蓄電池放電深度變化

表4 單次MAF實驗結(jié)果

Tab.4 Single MAF experiment results

為了說明MMAF方法較其他方法的優(yōu)勢,將本文混合儲能功率分配方法與EMD、VMD、小波變換方法進(jìn)行比較,得到如圖13所示結(jié)果。

不同功率分配方法下HESS成本見表5。通過所有實驗結(jié)果橫向比較可以看出,本文MMAF方法比EMD方法獲取的混合儲能LCC降低了87.6%、比VMD方法獲取的LCC降低了3.3%、比WT方法獲取的LCC降低了56.8%。證實了本文提出的MMAF方法在功率分配時能夠更好地削弱模態(tài)混疊現(xiàn)象,更加精確地分離高低頻分量,從而定容時可以依據(jù)功率分配結(jié)果最小化混合儲能的額定功率和容量。將儲能功率分配問題與定容問題聯(lián)合考慮,求得了混合儲能LCC的最優(yōu)解。

5.3 故障工況時HESS功率指令的分配

為了更好地契合混合儲能系統(tǒng)在實際運行工況中的某些特殊情況,以蓄電池為例,對其發(fā)生放電越限情況時HESS總功率指令如何分配做了進(jìn)一步研究。定義混合儲能中電池發(fā)生放電越限現(xiàn)象時的工況為“故障工況”。需要指出的是:不同于正常工況,故障工況下HESS的經(jīng)濟(jì)效益將不再是首要目標(biāo),如何將蓄電池因放電越限而無法響應(yīng)的這部分功率指令進(jìn)行再分配才是重中之重。因為一旦蓄電池因放電越限而無法響應(yīng)本身的功率指令時就會造成風(fēng)電平抑效果變差,風(fēng)電波動的增大嚴(yán)重時有可能導(dǎo)致風(fēng)機(jī)脫網(wǎng),造成更大的經(jīng)濟(jì)損失。故本節(jié)不再對HESS的經(jīng)濟(jì)性做分析,而是更關(guān)注于故障工況下的HESS功率指令的分配情況。

表5 不同功率分配方法下HESS成本的對比

Tab.5 Comparison of HESS cost under different power allocation methods

放電越限現(xiàn)象可以通過電池的SOC曲線很好地表征。本文規(guī)定蓄電池和超級電容的SOC上限為80%,下限為20%,因此當(dāng)蓄電池放電使得SOC低于20%時即認(rèn)為發(fā)生了放電越限現(xiàn)象。

由于本文算例所用的HESS總功率指令較為平緩,混合儲能系統(tǒng)在整個調(diào)度日內(nèi)沒有出現(xiàn)“故障工況”。因此本文在研究故障工況時采用另外一組HESS總功率指令作為研究數(shù)據(jù),如圖14所示。

以蓄電池為例,其在故障工況時HESS功率指令分配場景及機(jī)理如下:

1)同正常工況,通過MMAF對混合儲能總功率指令進(jìn)行濾波,Pearson相關(guān)系數(shù)作為判斷濾波參數(shù)的指標(biāo),獲得蓄電池與超級電容各自的功率指令。

圖14 混合儲能總功率指令

2)蓄電池在自身分得的功率指令下進(jìn)行充放電,部分時刻出現(xiàn)深度放電致使其SOC低于下限(20%),此時蓄電池進(jìn)入“故障工況”。

3)當(dāng)蓄電池SOC首次降至下限時,即刻起令蓄電池實際輸出功率bat=0,原本由蓄電池響應(yīng)的放電功率指令全部由超級電容承擔(dān),直到后續(xù)時刻風(fēng)電出力過多時需要給HESS充電以滿足平抑需求。

4)優(yōu)先向蓄電池充電以幫助其快速恢復(fù)SOC。在恢復(fù)過程中,蓄電池的SOC會第二次升至下限,即刻起不再對蓄電池實際輸出功率做限制,恢復(fù)蓄電池與超級電容各自原本分得的功率指令。

5)蓄電池SOC一次在下降的過程中達(dá)到下限值,一次在恢復(fù)的過程中達(dá)到下限值。記蓄電池SOC兩次達(dá)到下限值的時間區(qū)間為一個完整的故障工況運行時間,作為一次故障工況修復(fù)的標(biāo)志。恢復(fù)到正常工況后,進(jìn)入步驟1)。

需要說明的是:本文僅以蓄電池放電越限為例,當(dāng)超級電容放電越限時由蓄電池對其功率指令進(jìn)行響應(yīng),機(jī)理與上述完全一致,由于篇幅有限,本文不再贅述。對于蓄電池與超級電容同時發(fā)生放電越限的情況屬于小概率事件,不屬于本文研究重點,故不再對其討論。

圖15給出了調(diào)度日內(nèi)局部時段的故障工況運行過程中蓄電池與超級電容各自SOC的變化趨勢,以此更加直觀地展示上述機(jī)理。

圖15 故障工況下蓄電池和超級電容SOC曲線

從圖15可以看出在選取的時域段內(nèi)共發(fā)生了兩次蓄電池放電越限的現(xiàn)象:A1時段和B1時段。由局部放大圖可知對于A1時段,在45 min時蓄電池的SOC降至0.2,此時若不考慮蓄電池放電越限,其SOC曲線對應(yīng)圖中綠色點畫線。可以看出,蓄電池將繼續(xù)向外輸出功率,這極大地?fù)p害了電池壽命。若此時考慮蓄電池放電越限,其SOC曲線對應(yīng)圖中藍(lán)色實線。可以看出在45 min,由于放電越限迫使此時的蓄電池實際輸出功率bat=0,因此后一小段時間內(nèi)蓄電池的SOC也一直維持在0.2。在故障工況運行時間內(nèi)蓄電池的功率指令由超級電容進(jìn)行響應(yīng),對應(yīng)的超級電容SOC曲線是A2時段。由A2時段的放大圖可知,由于此時的超級電容需要承擔(dān)部分蓄電池的放電功率指令,因此超級電容在45 min時的SOC曲線較未考慮蓄電池放電越限情況下的SOC曲線進(jìn)一步下降。48 min時需要向HESS充電,通過A1時段可以看出,蓄電池的SOC曲線逐步上升。此時從A2時段超級電容的SOC變化趨勢可以看出:在不考慮蓄電池放電越限情況下的超級電容SOC在48 min已經(jīng)開始上升,說明48 min時超級電容已經(jīng)開始充電;而考慮了蓄電池放電越限情況下的超級電容SOC在蓄電池充電過程中保持不變,說明此過程超級電容并未得到充電,其在50 min才首次充電出現(xiàn)SOC上升的趨勢,體現(xiàn)了當(dāng)蓄電池SOC過低時優(yōu)先向蓄電池充電的原則。同理對于B1和B2時段,在故障工況運行時間內(nèi)蓄電池SOC始終保持在0.2,沒有再進(jìn)行放電行為,而當(dāng)需要向HESS充電時,不考慮蓄電池放電越限情況下的超級電容SOC在124 min已經(jīng)開始上升,說明124 min時超級電容已經(jīng)開始充電;而考慮了蓄電池放電越限情況下的超級電容則是在130 min才首次充電出現(xiàn)SOC上升的趨勢,同樣體現(xiàn)了蓄電池優(yōu)先充電的原則。

故障工況時經(jīng)過這種功率分配機(jī)理后的HESS功率分配結(jié)果如圖16所示。

圖16 故障工況時功率指令分配結(jié)果

圖16中的功率分配結(jié)果符合兩種儲能電池的功率特性。需要說明的是由于超級電容在調(diào)度日內(nèi)會承擔(dān)一些蓄電池的功率指令,因此圖16中故障工況時超級電容實際響應(yīng)的功率幅值更大一些。

圖17展示了故障工況時MMAF分解后各模態(tài)的頻譜關(guān)系。

圖17 故障工況時MMAF分解后各模態(tài)的頻譜

從圖17的各模態(tài)頻譜關(guān)系中可以明顯看出,故障工況時由于超級電容承擔(dān)了相當(dāng)一部分因蓄電池放電越限而本應(yīng)該由蓄電池承擔(dān)的功率指令,因此不可避免地出現(xiàn)了模態(tài)混疊現(xiàn)象(圖中紅色圈內(nèi)),這也符合故障工況時設(shè)定的HESS功率指令分配的機(jī)理。在故障工況下雖然出現(xiàn)了模態(tài)混疊,但通過這種“互為備用”的機(jī)理避免了功率指令因電池放電越限而無法響應(yīng)的情況,從而避免了因放電越限而沒有電池響應(yīng)功率指令導(dǎo)致風(fēng)電無法平抑現(xiàn)象的發(fā)生,遏制了由于風(fēng)功率波動較大導(dǎo)致風(fēng)機(jī)脫網(wǎng)的可能。因此與風(fēng)機(jī)脫網(wǎng)帶來的經(jīng)濟(jì)損失與電力事故相比,故障工況下模態(tài)混疊導(dǎo)致的HESS成本增加就顯得微不足道了。

6 結(jié)論

本文針對平抑風(fēng)電波動的混合儲能系統(tǒng)功率分配不精確、成本過高的問題,采用引入Pearson相關(guān)系數(shù)的MMAF方法得出蓄電池與超級電容各自分配的功率指令,并基于分配結(jié)果對混合儲能進(jìn)行了定容,得到如下結(jié)論:

1)引入Pearson相關(guān)系數(shù)的MMAF為濾波參數(shù)的選取提供了理論支撐,有效地減弱了模態(tài)混疊現(xiàn)象,實現(xiàn)了高精度的功率分配。

2)采用MMAF處理混合儲能功率分配問題及定容問題時混合儲能系統(tǒng)的LCC相較于傳統(tǒng)EMD、VMD、WT方法分別降低87.6%、3.3%、56.8%,在保證蓄電池使用壽命的同時,為解決風(fēng)電場側(cè)混合儲能電站的功率分配與容量配置問題提供了新思路。

[1] 劉其輝, 逄思敏, 吳林林, 等. 大規(guī)模風(fēng)電匯集系統(tǒng)電壓不平衡機(jī)理、因素及影響規(guī)律[J]. 電工技術(shù)學(xué)報, 2022, 37(21): 5435-5450. Liu Qihui, Pang Simin, Wu Linlin, et al. The mechanism, factors and influence rules of voltage imbalance in wind power integration areas[J]. Transactions of China Electrotechnical Society, 2022, 37(21): 5435-5450.

[2] 顧雪平, 白巖松, 李少巖, 等. 考慮風(fēng)電不確定性的電力系統(tǒng)恢復(fù)全過程兩階段魯棒優(yōu)化方法[J]. 電工技術(shù)學(xué)報, 2022, 37(21): 5462-5477. Gu Xueping, Bai Yansong, Li Shaoyan, et al. Two stage robust optimization method for the whole-process power system restoration considering wind power uncertainty[J]. Transactions of China Electrotechnical Society, 2022, 37(21): 5462-5477.

[3] de Siqueira L M S, Peng Wei. Control strategy to smooth wind power output using battery energy storage system: a review[J]. Journal of Energy Storage, 2021, 35: 102252.

[4] 齊先軍, 鄭夕煒, 王曉蓉, 等. 基于時頻分析的改進(jìn)小波包風(fēng)電功率波動平抑方法[J]. 太陽能學(xué)報, 2022, 43(7): 302-309. Qi Xianjun, Zheng Xiwei, Wang Xiaorong, et al. Improved wavelet packet method of smoothing wind power fluctuations based on time-frequency analysis[J]. Acta Energiae Solaris Sinica, 2022, 43(7): 302-309.

[5] 李政, 劉宏偉, 康健, 等. 混合儲能參與自動發(fā)電控制容量優(yōu)化配置[J]. 電氣技術(shù), 2021, 22(9): 34-40. Li Zheng, Liu Hongwei, Kang Jian, et al. Participation of hybrid energy storage in capacity optimization configuration of automatic generation control system[J]. Electrical Engineering, 2021, 22(9): 34-40.

[6] Wang Huaqing, Xie Zhuoshi, Pu Lei, et al. Energy management strategy of hybrid energy storage based on Pareto optimality[J]. Applied Energy, 2022, 327: 120095.

[7] Zhou Shiyao, Chen Ziqiang, Huang Deyang, et al. Model prediction and rule based energy management strategy for a plug-in hybrid electric vehicle with hybrid energy storage system[J]. IEEE Transactions on Power Electronics, 2021, 36(5): 5926-5940.

[8] Xu Yanhui, Xu Yijia, Huang Yan. Generation of typical operation curves for hydrogen storage applied to the wind power fluctuation smoothing mode[J]. Global Energy Interconnection, 2022, 5(4): 353-361.

[9] Jin He, Liu Pei, Li Zheng. Dynamic modeling and design of a hybrid compressed air energy storage and wind turbine system for wind power fluctuation reduction[J]. Computers & Chemical Engineering, 2019, 122: 59-65.

[10] 孫玉樹, 唐西勝, 孫曉哲, 等. 基于MPC-HHT的多類型儲能協(xié)調(diào)控制策略研究[J]. 中國電機(jī)工程學(xué)報, 2018, 38(9): 2580-2588, 2826. Sun Yushu, Tang Xisheng, Sun Xiaozhe, et al. Research on multi-type energy storage coordination control strategy based on MPC-HHT[J]. Proceedings of the CSEE, 2018, 38(9): 2580-2588, 2826.

[11] 袁鐵江, 張昱, 栗磊, 等. 計及功率密度約束含氫儲能的預(yù)裝式多元儲能電站容量優(yōu)化配置研究[J]. 電工技術(shù)學(xué)報, 2021, 36(3): 496-506. Yuan Tiejiang, Zhang Yu, Li Lei, et al. Capacity optimization configuration of pre-installed multi-energy storage power station considering power density constrained hydrogen storage[J]. Transactions of China Electrotechnical Society, 2021, 36(3): 496-506.

[12] 李奇, 趙淑丹, 蒲雨辰, 等. 考慮電氫耦合的混合儲能微電網(wǎng)容量配置優(yōu)化[J]. 電工技術(shù)學(xué)報, 2021, 36(3): 486-495. Li Qi, Zhao Shudan, Pu Yuchen, et al. Capacity optimization of hybrid energy storage microgrid considering electricity-hydrogen coupling[J]. Transactions of China Electrotechnical Society, 2021, 36(3): 486-495.

[13] 錢韋廷, 趙長飛, 萬燦, 等. 基于概率預(yù)測的混合儲能平抑風(fēng)電波動隨機(jī)優(yōu)化調(diào)控方法[J]. 電力系統(tǒng)自動化, 2021, 45(18): 18-27. Qian Weiting, Zhao Changfei, Wan Can, et al. Probabilistic forecasting based stochastic optimal dispatch and control method of hybrid energy storage for smoothing wind power fluctuations[J]. Automation of Electric Power Systems, 2021, 45(18): 18-27.

[14] 李建林, 牛萌, 周喜超, 等. 能源互聯(lián)網(wǎng)中微能源系統(tǒng)儲能容量規(guī)劃及投資效益分析[J]. 電工技術(shù)學(xué)報, 2020, 35(4): 874-884. Li Jianlin, Niu Meng, Zhou Xichao, et al. Energy storage capacity planning and investment benefit analysis of micro-energy system in energy interconnection[J]. Transactions of China Electrote-chnical Society, 2020, 35(4): 874-884.

[15] Yang Yuqing, Bremner S, Menictas C, et al. Modelling and optimal energy management for battery energy storage systems in renewable energy systems: a review[J]. Renewable and Sustainable Energy Reviews, 2022, 167: 112671.

[16] 李建林, 李雅欣, 劉海濤, 等. 計及儲能電站安全性的功率分配策略研究[J]. 電工技術(shù)學(xué)報, 2022, 37(23): 5976-5986. Li Jianlin, Li Yaxin, Liu Haitao, et al. Research on power distribution strategy considering the safety of energy storage power station[J]. Transactions of China Electrotechnical Society, 2022, 37(23): 5976-5986.

[17] 郭玲娟, 魏斌, 韓肖清, 等. 基于集合經(jīng)驗?zāi)B(tài)分解的交直流混合微電網(wǎng)混合儲能容量優(yōu)化配置[J]. 高電壓技術(shù), 2020, 46(2): 527-537. Guo Lingjuan, Wei Bin, Han Xiaoqing, et al. Capacity optimal configuration of hybrid energy storage in hybrid AC/DC micro-grid based on ensemble empirical mode decomposition[J]. High Voltage Engineering, 2020, 46(2): 527-537.

[18] 李亞楠, 王倩, 宋文峰, 等. 基于變分模態(tài)分解和Hilbert變換的平滑風(fēng)電出力混合儲能容量優(yōu)化配置[J]. 電測與儀表, 2019, 56(1): 82-88, 95. Li Yanan, Wang Qian, Song Wenfeng, et al. Capacity configuration of hybrid energy storage for smoothing wind power outputs based on variational mode decomposition and Hilbert transform[J]. Electrical Measurement & Instrumentation, 2019, 56(1): 82-88, 95.

[19] 李鑫, 王娟, 邱亞, 等. 基于VMD的混合儲能容量優(yōu)化配置[J]. 太陽能學(xué)報, 2022, 43(2): 88-96. Li Xin, Wang Juan, Qiu Ya, et al. Optimal allocation of hybrid energy storage capacity based on variational mode decomposition[J]. Acta Energiae Solaris Sinica, 2022, 43(2): 88-96.

[20] 杜佳耘, 雷勇, 李永凱, 等. 基于參數(shù)優(yōu)化變分模態(tài)分解的混合儲能功率分配策略[J]. 現(xiàn)代電力, 2021, 38(1): 51-59. Du Jiayun, Lei Yong, Li Yongkai, et al. Hybrid energy storage strategy based on parameter optimized variational mode decomposition[J]. Modern Electric Power, 2021, 38(1): 51-59.

[21] 李建林, 郭斌琪, 牛萌, 等. 風(fēng)光儲系統(tǒng)儲能容量優(yōu)化配置策略[J]. 電工技術(shù)學(xué)報, 2018, 33(6): 1189-1196. Li Jianlin, Guo Binqi, Niu Meng, et al. Optimal configuration strategy of energy storage capacity in wind/PV/storage hybrid system[J]. Transactions of China Electrotechnical Society, 2018, 33(6): 1189-1196.

[22] Zhu Honglu, Hou Ruyin, Jiang Tingting, et al. Research on energy storage capacity configuration for PV power plants using uncertainty analysis and its applications[J]. Global Energy Interconnection, 2021, 4(6): 608-618.

[23] Hou Hui, Xu Tao, Wu Xixiu, et al. Optimal capacity configuration of the wind-photovoltaic-storage hybrid power system based on gravity energy storage system[J]. Applied Energy, 2020, 271: 115052.

[24] 王再闖, 陳來軍, 李笑竹, 等. 基于合作博弈的產(chǎn)銷者社區(qū)分布式光伏與共享儲能容量優(yōu)化[J]. 電工技術(shù)學(xué)報, 2022, 37(23): 5922-5932. Wang Zaichuang, Chen Laijun, Li Xiaozhu, et al. Capacity optimization of distributed PV and shared energy storage of prosumer community based on cooperative game[J]. Transactions of China Electrotechnical Society, 2022, 37(23): 5922-5932.

[25] Ye Lin, Zhang Cihang, Tang Yong, et al. Hierarchical model predictive control strategy based on dynamic active power dispatch for wind power cluster integration[J]. IEEE Transactions on Power Systems, 2019, 34(6): 4617-4629.

[26] Guo Tingting, Liu Youbo, Zhao Junbo, et al. A dynamic wavelet-based robust wind power smoothing approach using hybrid energy storage system[J]. International Journal of Electrical Power & Energy Systems, 2020, 116: 105579.

[27] 陳亞愛, 林演康, 王賽, 等. 基于濾波分配法的混合儲能優(yōu)化控制策略[J]. 電工技術(shù)學(xué)報, 2020, 35(19): 4009-4018. Chen Yaai, Lin Yankang, Wang Sai, et al. Optimal control strategy of hybrid energy storage based on filter allocation method[J]. Transactions of China Electrotechnical Society, 2020, 35(19): 4009-4018.

[28] 張明龍, 張振宇, 高源, 等. 基于變分模態(tài)分解的暫態(tài)擾動波形去噪算法[J]. 電力系統(tǒng)保護(hù)與控制, 2022, 50(8): 43-49. Zhang Minglong, Zhang Zhenyu, Gao Yuan, et al. Transient disturbance waveform denoising algorithm based on variational mode decomposition[J]. Power System Protection and Control, 2022, 50(8): 43-49.

[29] Omar A I, Ali Z M, Al-Gabalawy M, et al. Multi-objective environmental economic dispatch of an electricity system considering integrated natural gas units and variable renewable energy sources[J]. Mathematics, 2020, 8(7): 1100.

[30] 葛玉友, 尚策. 壽命約束的儲能規(guī)劃[J]. 中國電機(jī)工程學(xué)報, 2020, 40(19): 6150-6161. Ge Yuyou, Shang Ce. Energy storage planning constrained by its life[J]. Proceedings of the CSEE, 2020, 40(19): 6150-6161.

[31] 丁明, 吳杰, 張晶晶. 面向風(fēng)電平抑的混合儲能系統(tǒng)容量配置方法[J]. 太陽能學(xué)報, 2019, 40(3): 593-599. Ding Ming, Wu Jie, Zhang Jingjing. Capacity optimization method of hybrid energy storage system for wind power smoothing[J]. Acta Energiae Solaris Sinica, 2019, 40(3): 593-599.

[32] 馬蘭, 謝麗蓉, 葉林, 等. 基于混合儲能雙層規(guī)劃模型的風(fēng)電波動平抑策略[J]. 電網(wǎng)技術(shù), 2022, 46(3): 1016-1029. Ma Lan, Xie Lirong, Ye Lin, et al. Wind power fluctuation suppression strategy based on hybrid energy storage Bi-level programming model[J]. Power System Technology, 2022, 46(3): 1016-1029.

Research on Hybrid Energy Storage Power Allocation and Capacity Determination Based on Multiple Moving Average Filtering

Tian Bowen Zhang Zhiyu Yang Mengfei

(College of Electrical Engineering Xi′an University of Technology Xi′an 710048 China)

When the hybrid energy storage system (HESS) allocates power in the scenario of suppressing wind power fluctuations, the power distribution is inaccurate due to mode aliasing, which leads to the poor economy of the hybrid energy storage system. To address these issues, this paper suggests a method of multiple moving average filtering (MMAF) to decompose the total power command of hybrid energy storage.

Firstly, the minimum flat total power command of hybrid energy storage meeting the flat requirements was obtained, and the obtained total power command of hybrid energy storage was filtered through MMAF to make the battery respond to the filtered low frequency part and the supercapacitor respond to the high frequency part. Pearson correlation coefficient was used as the index to quantify modal aliasing, and Pearson correlation coefficient was used as the basis to determine the optimal filtering sliding window size and the optimal filtering times to ensure that the maximum limit of modal aliasing can be achieved under the sliding window size and filtering times. Taked the power command of the battery and the supercapacitor as the feasible region of the energy storage capacity, the minimum rated power and rated capacity of the hybrid energy storage system can be determined by comprehensively considering the battery state of charge (SOC) constraints. Secondly, based on the equivalent operating time, a quantitative model of battery life cycle was established, which provided a basis for the economic analysis of hybrid energy storage system. The cumulative effect of battery charging and discharging on its life loss was considered in the model, and the actual number of battery cycles was converted into the equivalent number of cycles at the rated discharge depth. Finally, the incremental factor of battery maintenance cost was introduced into the hybrid energy storage life cycle cost model (LCC) to describe the nonlinear increase of maintenance cost caused by the cumulative effect of battery aging every ten years, so as to better conform to the actual operation of the energy storage battery and improve the practicability of the LCC model. The effectiveness of the proposed method to reduce the comprehensive cost of hybrid energy storage was verified by simulation experiments.

According to the actual simulation experiment, the method in this paper can well smooth the power fluctuation of wind power grid connection, made it meet the grid connection requirements, and the hybrid energy storage configuration result was better than other configuration methods. The simulation results under different working conditions showed that MMAF method can ensure the battery operation life during power distribution, and can minimize the mode aliasing phenomenon and reduce the comprehensive cost of hybrid energy storage. Under normal working conditions, the LCC of MMAF was reduced by 87.6% compared with the empirical mode decomposition method (EMD). Compared with the variational mode decomposition method (VMD), the LCC of MMAF was reduced by 3.3%. Compared with the wavelet transform (WT) method, the LCC of MMAF was reduced by 56.8%. In case of failure, the battery and supercapacitor form a standby working mode for each other. The total power command of mixed energy storage was redistributed by MMAF method, which avoided the situation that the power command cannot respond due to battery failure, curbs the possibility of wind turbine disconnection, and avoided greater economic losses.

Hybrid energy storage, power distribution, mode aliasing, MMAF algorithm, quantitative model of battery life cycle

10.19595/j.cnki.1000-6753.tces.222348

TM711

國家自然科學(xué)基金聯(lián)合基金重點資助項目(U21A20485)。

2022-12-21

2023-02-12

田博文 男,1997年生,碩士研究生,研究方向為新能源側(cè)儲能優(yōu)化配置與控制。E-mail:932725835@qq.com

張志禹 男,1966年生,博士,教授,研究方向為電網(wǎng)信息處理與深度學(xué)習(xí)。E-mail:zhangzhiyu@xaut.edu.cn(通信作者)

(編輯 赫 蕾)

猜你喜歡
指令分配模態(tài)
聽我指令:大催眠術(shù)
應(yīng)答器THR和TFFR分配及SIL等級探討
遺產(chǎn)的分配
一種分配十分不均的財富
ARINC661顯控指令快速驗證方法
績效考核分配的實踐與思考
LED照明產(chǎn)品歐盟ErP指令要求解讀
電子測試(2018年18期)2018-11-14 02:30:34
國內(nèi)多模態(tài)教學(xué)研究回顧與展望
基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識別
由單個模態(tài)構(gòu)造對稱簡支梁的抗彎剛度
計算物理(2014年2期)2014-03-11 17:01:39
主站蜘蛛池模板: 色婷婷亚洲综合五月| 狠狠色噜噜狠狠狠狠色综合久| 欧美人与动牲交a欧美精品| 天堂av高清一区二区三区| 日韩精品成人网页视频在线| 美女被狂躁www在线观看| 亚洲中文字幕无码爆乳| 国产精品私拍在线爆乳| 国产成人精彩在线视频50| 精品视频第一页| 香蕉在线视频网站| 在线观看欧美精品二区| 精品国产网| 亚洲伊人电影| 亚洲欧美一区在线| 精品国产美女福到在线不卡f| 在线a视频免费观看| 成人av专区精品无码国产| 婷婷亚洲综合五月天在线| 毛片在线看网站| 青草精品视频| 亚洲AV无码乱码在线观看代蜜桃 | 红杏AV在线无码| 欧美日韩北条麻妃一区二区| 国产精品久久久久鬼色| 尤物国产在线| a在线观看免费| 亚洲欧美一区二区三区图片| 国产系列在线| 亚洲精品国产综合99| 色妺妺在线视频喷水| 国产网友愉拍精品视频| 精品色综合| 久久精品娱乐亚洲领先| 国产浮力第一页永久地址| 无码精品国产dvd在线观看9久| 强奷白丝美女在线观看| 九九九国产| 在线视频一区二区三区不卡| 成人精品免费视频| 中文字幕免费播放| 国产一区二区三区精品欧美日韩| 亚洲人成人无码www| 精品国产成人国产在线| 成人精品视频一区二区在线| 欧美成人国产| 人人爽人人爽人人片| 天堂在线www网亚洲| 国产99视频精品免费观看9e| 丁香五月激情图片| 色成人亚洲| 欧美中日韩在线| 国产人前露出系列视频| 亚洲a免费| 欧美亚洲欧美区| 成人无码区免费视频网站蜜臀| 国产在线视频二区| 日韩欧美中文| 国内自拍久第一页| 天天操精品| 亚洲日韩第九十九页| 亚洲精品无码人妻无码| 亚洲综合二区| 亚洲第一中文字幕| 91在线播放免费不卡无毒| 天堂av综合网| 99re经典视频在线| 欧美精品高清| 亚洲天堂久久久| a级毛片毛片免费观看久潮| a国产精品| 亚洲第一极品精品无码| 成人午夜视频在线| 91小视频在线观看| 亚洲色图另类| 奇米精品一区二区三区在线观看| 国产AV无码专区亚洲精品网站| 在线中文字幕网| 一级毛片免费不卡在线| 亚洲成人精品在线| 无码'专区第一页| 播五月综合|