CN108470220A - 考虑功率变化率限制的混合储能系统能量管理优化方法 - Google Patents
考虑功率变化率限制的混合储能系统能量管理优化方法 Download PDFInfo
- Publication number
- CN108470220A CN108470220A CN201810098996.2A CN201810098996A CN108470220A CN 108470220 A CN108470220 A CN 108470220A CN 201810098996 A CN201810098996 A CN 201810098996A CN 108470220 A CN108470220 A CN 108470220A
- Authority
- CN
- China
- Prior art keywords
- soc
- power
- bat
- dem
- electric vehicle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004146 energy storage Methods 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000005457 optimization Methods 0.000 title claims abstract description 27
- 239000003990 capacitor Substances 0.000 claims abstract description 54
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 claims abstract description 44
- 229910052744 lithium Inorganic materials 0.000 claims abstract description 44
- 230000007704 transition Effects 0.000 claims abstract description 40
- 239000011159 matrix material Substances 0.000 claims abstract description 24
- 230000008569 process Effects 0.000 claims abstract description 15
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 14
- 230000001133 acceleration Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 8
- 238000012546 transfer Methods 0.000 claims description 8
- 230000009471 action Effects 0.000 claims description 5
- 238000005096 rolling process Methods 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 6
- 230000005611 electricity Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 239000005431 greenhouse gas Substances 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- BULVZWIRKLYCBC-UHFFFAOYSA-N phorate Chemical compound CCOP(=S)(OCC)SCSCC BULVZWIRKLYCBC-UHFFFAOYSA-N 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/29—Graphical models, e.g. Bayesian networks
- G06F18/295—Markov models or related models, e.g. semi-Markov models; Markov random fields; Networks embedding Markov models
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J7/00—Circuit arrangements for charging or depolarising batteries or for supplying loads from batteries
- H02J7/007—Regulation of charging or discharging current or voltage
- H02J7/0071—Regulation of charging or discharging current or voltage with a programmable schedule
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J7/00—Circuit arrangements for charging or depolarising batteries or for supplying loads from batteries
- H02J7/34—Parallel operation in networks using both storage and other dc sources, e.g. providing buffering
- H02J7/345—Parallel operation in networks using both storage and other dc sources, e.g. providing buffering using capacitors as storage or buffering devices
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M3/00—Conversion of dc power input into dc power output
- H02M3/22—Conversion of dc power input into dc power output with intermediate conversion into ac
- H02M3/24—Conversion of dc power input into dc power output with intermediate conversion into ac by static converters
- H02M3/28—Conversion of dc power input into dc power output with intermediate conversion into ac by static converters using discharge tubes with control electrode or semiconductor devices with control electrode to produce the intermediate ac
- H02M3/325—Conversion of dc power input into dc power output with intermediate conversion into ac by static converters using discharge tubes with control electrode or semiconductor devices with control electrode to produce the intermediate ac using devices of a triode or a transistor type requiring continuous application of a control signal
- H02M3/335—Conversion of dc power input into dc power output with intermediate conversion into ac by static converters using discharge tubes with control electrode or semiconductor devices with control electrode to produce the intermediate ac using devices of a triode or a transistor type requiring continuous application of a control signal using semiconductor devices only
- H02M3/3353—Conversion of dc power input into dc power output with intermediate conversion into ac by static converters using discharge tubes with control electrode or semiconductor devices with control electrode to produce the intermediate ac using devices of a triode or a transistor type requiring continuous application of a control signal using semiconductor devices only having at least two simultaneously operating switches on the input side, e.g. "double forward" or "double (switched) flyback" converter
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Power Engineering (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Development Economics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Game Theory and Decision Science (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
- Secondary Cells (AREA)
Abstract
本发明涉及一种考虑功率变化率限制的混合储能系统能量管理优化方法，包括：根据电动汽车的行驶工况数据以及电动汽车的基本参数，确定电动汽车总需求功率；根据电动汽车的需求功率具有马尔可夫性，得到需求功率Pdem的状态转移概率矩阵，定义需求功率的状态转移概率Pdem，imj；根据锂电池和超级电容的工作特性以及实际的工况确定各参数的范围；建立系统的状态转移方程；求得状态转移概率矩阵，由需求功率Pdem和锂电池功率Pbat得到超级电容功率Puc，进而根据当前时刻的SOCuc得到下一时刻的SOCuc；以降低系统的能量损耗和维持超级电容SOC为目标，建立多目标的成本函数R；通过马尔可夫决策过程算法得到最优的策略；进一步限制锂电池功率变化率。
Description
技术领域
本发明属于能量管理系统领域，涉及到混合储能系统以及能量管理优化算法。
背景技术
由于对温室气体排放和化石燃料消费不断增长的担忧，电动汽车的发展受到了越来越多的关注，可以预见电动汽车将在未来的交通运输中发挥重要作用。锂电池作为电动汽车的主要储能器件，其必须满足电动汽车的功率和能量需求，这导致锂电池的尺寸过大。为了提高锂电池寿命，基于锂电池和超级电容的混合储能系统被用作电动汽车的储能器件。混合储能系统具有高能量以及高功率密度互补的特性，具有较大的储能和较好的功率性能。
由于存在两个相互独立的储能器件，电动汽车的需求功率在两个储能器件之间的合理分配成为一个主要的挑战，因此，能量管理优化策略成为混合储能系统能量控制的关键。能量管理优化策略的主要目的，一是通过超级电容吸收瞬时峰值功率，减少锂电池输入输出的电量，延长锂电池的寿命；二是减少混合储能系统的能量损失，提高系统的效率；三是保证超级电容的SOC始终处于合适状态，使超级电容可以及时的处理电动汽车在加速或者减速过程中产生的瞬时峰值功率。然而现有的能量管理优化策略很少将锂电池的功率变化率的限制考虑在内，使得能量控制策略并未达到最优的效果，因此能量管理优化算法的效果还有很大的提升空间。
发明内容
本发明的目的是提供一种考虑功率变化率限制的混合储能系统能量管理优化方法，在马尔可夫决策过程优化算法的基础上，通过合理的限制锂电池功率的变化率，提高混合储能系统的工作效率并延长锂电池使用寿命。方案如下：
一种考虑功率变化率限制的混合储能系统能量管理优化方法，包括如下的步骤：
(1)根据电动汽车的行驶工况数据以及电动汽车的基本参数，确定电动汽车总需求功率Pdem：
Pdem＝Pf+Pa+Pi
式中Pf是克服滚动阻力产生的功率、Pa是克服加速阻力产生的功率、Pi是克服空气阻力产生的功率、Pdem是电动汽车的需求功率。
(2)根据电动汽车的需求功率具有马尔可夫性，得到需求功率Pdem的状态转移概率矩阵，定义需求功率的状态转移概率Pdem,imj为：
式中，v为当前时刻电动汽车的速度，需求功率Pdem和电动汽车速度v被离散为有限值，即i,j＝1,2,…,Np；m＝1,2,…,Nv，当电动汽车的速度v＝vm时，汽车在t时刻的需求功率为的情况下，在t+1时刻的需求功率为的条件概率。
(3)根据锂电池和超级电容的工作特性以及实际的工况确定各参数的范围，各参数值限定范围为：
Pbat,min≤Pbat≤Pbat,max
Puc,min≤Puc≤Puc,max
SOCuc,min≤SOCuc≤SOCuc,max
式中Pbat为锂电池的功率、Puc为超级电容的功率、SOCuc为超级电容的荷电状态。
(4)已知当前时刻的状态St＝(v(t),Pdem(t),SOCuc(t))以及动作集a＝Pbat(t)，可以根据状态转移方程求出下一时刻的状态St+1＝(v(t+1),Pdem(t+1),SOCuc(t+1))，系统的状态转移方程为：
Puc(t)＝Pdem(t)-Pbat(t)
式中Iuc为超级电容的电流、Vuc为超级电容的端电压、Quc为超级电容的容量、Δt为采样时间间隔。
(5)根据需求功率的状态转移概率矩阵和系统的状态转移方程，可以求出状态S的转移概率矩阵Psoc,imj：
Psoc,imj＝P[St+1＝SOCuc,j|St＝SOCuc,i,At＝a,v＝vm]
公式的含义是当电动汽车的速度v＝vm、储能系统采取动作At＝a时，超级电容的SOC在t时刻的状态为SOCuc,i的情况下，在t+1时刻的状态为SOCuc,j的条件概率。
由需求功率Pdem和锂电池功率Pbat得到超级电容功率Puc，进而根据当前时刻的SOCuc得到下一时刻的SOCuc；
(6)以降低系统的能量损耗和维持超级电容SOC为目标，建立多目标的成本函数R，因为两个目标的优化相互矛盾，为实现均衡的优化目标，引入权重系数来实现多目标优化；
(7)已知状态S的转移概率矩阵、目标函数R和折扣因子γ，通过马尔可夫决策过程算法得到最优的策略π*(a|S)；
(8)在通过最优策略π*(a|S)确定下一时刻动作a＝Pbat(t+1)时，进一步限制锂电池功率变化率ΔPbat：当通过最优策略确定下一时刻锂电池功率时，需要与前一时刻的锂电池功率值进行比较，若功率变化率|ΔPbat|小于限定值a，则功率变化率ΔPbat无需受到限制，下一时刻锂电池功率Pbat(t+1)为有效的；若功率变化率|ΔPbat|大于限定值a，而此时SOCuc值在限定范围之内，即SOCuc,ref-b≤SOCuc≤SOCuc,ref+c，则功率变化率ΔPbat需要受到限制，下一时刻锂电池功率Pbat(t+1)＝Pbat(t)±a；若功率变化率|ΔPbat|大于限定值a，且此时SOCuc值超出限定范围，即SOCuc≤SOCuc,ref-b或SOCuc≥SOCuc,ref+c，则功率变化率ΔPbat受到的限制变小，下一时刻锂电池功率Pbat(t+1)＝Pbat(t)±a+0.2(ΔPbat±a)。其中a为功率值，b、c为SOCuc值，均为正数，ΔPbat＝Pbat(t+1)-Pbat(t)。
本发明的技术效果如下：
1.马尔科夫决策过程算法兼顾了动态规划求解全局最优和马尔可夫决策处理随机过程的优点，可以有效地实现混合储能系统的能量优化管理。
2.将锂电池功率的变化率的限制纳入优化策略中，改进了马尔科夫决策过程算法，进一步提高了混合储能系统的工作效率并延长了锂电池使用寿命。
附图说明
图1：UDDS工况速度数据。
图2：UDDS工况电动汽车需求功率。
图3：车速35Km/h时需求功率转移概率矩阵。
图4：混合储能系统结构图。
图5：双向半桥DC/DC变换器。
图6：车速35Km/h时优化策略结果。
图7：马尔可夫决策过程示意图。
图8：功率变化率的限制示意图。
图9：功率变化率限制前后电池功率。
图10：通过优化算法得到的功率分配策略。
图11：超级电容SOC变化曲线。
图12：功率变化率限制前后电池衰退对比图。
具体实施方式
本发明提供了一种考虑功率变化率限制的混合储能系统能量管理优化方法，下面先概述本发明的考虑功率变化率限制的混合储能系统能量管理优化方法，包括如下的步骤：
(1)根据电动汽车的行驶工况数据以及电动汽车的基本参数，可以确定电动汽车总需求功率。汽车在平坦路面行驶时受到滚动阻力、空气阻力和加速阻力，这三个力产生的功率的和可以得到电动汽车的需求功率：
Pdem＝Pf+Pa+Pi
式中Pf是克服滚动阻力产生的功率、Pa是克服加速阻力产生的功率、Pi是克服空气阻力产生的功率、Pdem是电动汽车的需求功率。
因为需求功率具有马尔可夫性，即下一个状态仅与当前状态有关，而与以前的状态无关，所以根据马尔可夫模型可以得到需求功率Pdem的状态转移概率矩阵，定义状态转移概率Pdem,imj为：
式中条件概率Pdem,imj是需求功率的转移概率、v为当前时刻电动汽车的速度，其中需求功率Pdem和电动汽车速度v被离散为有限值。即i,j＝1,2,…,Np；m＝1,2,…,Nv。公式的含义是当电动汽车的速度v＝vm时，汽车在t时刻的需求功率为的情况下，在t+1时刻的需求功率为的条件概率。
(3)根据锂电池和超级电容的工作特性以及实际的工况确定各参数的范围，考虑到各器件特性的限制，使系统处于最佳工作状态，各参数值必须严格限定在设定参数范围内。则各参数范围为：
Pbat,min≤Pbat≤Pbat,max
Puc,min≤Puc≤Puc,max
SOCuc,min≤SOCuc≤SOCuc,max
式中Pbat为锂电池的功率、Puc为超级电容的功率、SOCuc为超级电容的荷电状态。
(4)已知当前时刻的状态St＝(v(t),Pdem(t),SOCuc(t))以及动作集a＝Pbat(t)，可以根据状态转移方程求出下一时刻的状态St+1＝(v(t+1),Pdem(t+1),SOCuc(t+1))。系统的状态转移方程为：
Puc(t)＝Pdem(t)-Pbat(t)
式中Iuc为超级电容的电流、Vuc为超级电容的端电压、Quc为超级电容的容量、Δt为采样时间间隔。
(5)根据需求功率的状态转移概率矩阵和系统的状态转移方程，可以求出状态S的转移概率矩阵。状态S的转移概率Psoc,imj定义为：
Psoc,imj＝P[St+1＝SOCuc,j|St＝SOCuc,i,At＝a,v＝vm]
式中条件概率Psoc,imj为状态S的转移概率。公式的含义是当电动汽车的速度v＝vm、储能系统采取动作At＝a时，超级电容的SOC在t时刻的状态为SOCuc,i的情况下，在t+1时刻的状态为SOCuc,j的条件概率。
由需求功率Pdem和锂电池功率Pbat可以得到超级电容功率Puc，进而根据当前时刻的SOCuc得到下一时刻的SOCuc，可以表示为：
SOCuc(t+1)＝f(SOCuc(t),Pdem(t),Pbat(t))
(6)以降低系统的能量损耗和维持超级电容的SOC为目标，建立多目标的成本函数，因为两个目标的优化相互矛盾，为了实现均衡的优化目标，引入权重系数来实现多目标优化。目标函数的具体表达如下：
R＝ω1Rloss+ω2Rsoc
式中ω1、ω2分别代表Rloss、Rsoc的权重系数。
(7)已知状态S的转移概率矩阵、目标函数R和折扣因子γ，可以通过马尔可夫决策过程算法得到最优的策略π*(a|S)，所谓策略π(a|S)是指状态到动作的映射，即策略π(a|S)是指在每个状态S处采取一个确定的动作a。最优策略π*(a|S)定义如下：
式中Vπ(s)是状态值函数，argmax是指当∑s′∈SPsoc,imjV*(s′)取最大值时，策略π(a|S)的取值，因为状态值函数是与策略π(a|S)相对应的，所以最优策略π*(a|S)可以通过最大化的状态值函数V*(s′)来决定。累计成本具有随机性，当采用策略π(a|S)时，定义累计成本在状态S处的期望值为状态值函数：
状态值函数Vπ(s)是与策略π(a|S)相对应的，这是因为策略π(a|S)决定了累计成本的状态分布。
(8)通过马尔可夫决策过程算法得到最优策略π*(a|S)后，由于动作集a是离散的，同时成本函数中没有对锂电池功率Pbat的变化率进行限制，这导致对锂电池功率Pbat的控制结果变化过于频繁，产生一定的抖动，这对延长锂电池寿命和降低系统能量损耗十分不利。在通过最优策略π*(a|S)确定下一时刻动作a＝Pbat(t+1)时，进一步限制锂电池功率变化率ΔPbat。
式中a为功率值，b、c为SOCuc值，均为正数。ΔPbat＝Pbat(t+1)-Pbat(t)。
以下结合附图和实例对本发明进行详细的描述。本发明具体实施步骤如下：
1.根据工况确定需求功率数据曲线。以UDDS工况为例，采样时间为1S，图1给出了UDDS工况的速度数据，最大行驶速度56.7km/h，行驶距离7510m，最大加速度0.924m/s2,最小加速度-0.924m/s2。根据工况中的速度数据以及电动汽车的基本参数可以求出电动汽车的需求功率，详细公式如下：
Pf＝mgfv
式中f为滚动阻力系数、δ为汽车旋转质量系数、CD为空气阻力系数、A为汽车迎风面积。
图2给出了通过公式计算得出的电动汽车需求功率曲线。
2.确定需求功率转移概率矩阵。通过临近法，根据需求功率数据可以得到需求功率的转移概率矩阵，转移概率值由最大似然估计法确定，公式如下：
式中mimj表示在车速v＝m时，需求功率转移到发生的次数；mim表示在车速为v＝m时，需求功率发生转移的总的次数，且
采用该方法可获得不同车速下需求功率转移矩阵。此外，转移概率矩阵满足以车速35Km/h为例，图3给出了需求功率转移概率矩阵分布。
3.根据各部件即工作特性确定系统工作范围。为了防止电池和超级电容过充过放或功率过大损害混合储能系统器件，必须对SOCuc以及电池和超级电容的功率进行严格限制。
在一个完整的工况过程中，超级电容释放的能量与吸收的能量几乎相等，类似于电路中的电感和电容，因此电动汽车所需的能量必须能够全部由电池提供，电池的输出功率必须能满足在平地上电动汽车以最大速度行驶时所需的功率，因此电池最大输出功率为：
根据电池特性，设定电池最大输入功率为：
超级电容主要处理电动汽车在加速或者减速过程中的瞬时峰值功率，瞬时峰值功率主要与电动汽车的加速度有关，超级电容吸收或输出的功率应大于完整的工况中克服加速度阻力产生的功率的最大值。超级电容吸收或输出的最大功率为：
Puc,max＝-Puc,min＝max|Pa|
根据超级电容特性可知，输出电压与其SOC线性相关，为了避免超级电容输出电压过低，同时也为了防止超级电容过充过放，需要对SOCuc进行限制。本发明设定SOCuc,min＝0.3、SOCuc,max＝0.9。
4.列出系统状态转移方程。已知当前时刻的状态St＝(v(t),Pdem(t),SOCuc(t))以及动作集a＝Pbat(t)，可以根据状态转移方程求出下一时刻的状态St+1＝(v(t+1),Pdem(t+1),SOCuc(t+1))。系统的状态转移方程为：
Puc(t)＝Pdem(t)-Pbat(t)
系统中各状态值皆为离散量，例如在UDDS工况中，车速范围为0～56.7km/h，将其离散为0、5、15、25、35、45、55的网格值，即Nv＝7；需求功率的范围为-9272.25～13850.91，以500的间距离散为60个网格值，即Np＝60；SOCuc的范围为0.3～0.9，以0.1的间距离散为60个网格值。
5.确定状态转移概率矩阵。步骤2中只是确定了需求功率状态转移概率矩阵，而要通过马尔可夫决策过程算法算出最优控制策略，就必须以需求功率状态转移概率为基础，确定在某个决策a下的状态转移概率。
根据需求功率状态转移概率矩阵和系统状态转移方程，可以求出状态S的转移概率矩阵。需求功率的随机性可以看作是干扰，因此带入了SOCuc状态变量的不确定性。例如在t时刻Pdem＝1kw，t+1时刻Pdem＝1kw的概率为0.2，Pdem＝2kw的概率为0.8；可以通过状态转移方程SOCuc(t+1)＝SOCuc(t)-f(Pdem(t))求出在t时刻SOCuc为0.5时，t+1时刻SOCuc为某个值的概率。例如，当Pdem的值从1kw转移到2kw时，如果SOCuc的值从0.5转移到了0.49，那么概率Psoc,0.5m0.49＝P(SOCuc(t+1)＝0.49|SOCuc(t)＝0.5)＝P(Pdem(t+1)＝2kw|Pdem(t)＝1kw)＝0.8。
6.确定混合储能系统结构。在本发明中，选用图4所示的混合储能系统的结构，因为它提供了操作的最大灵活性，同时最小化了从电池和超级电容到负载的转换阶段的数量，电池和超级电容的电压可以独立地保持低于直流母线的电压，此外可以充分利用超级电容内储存的能量。
选用双向半桥DC-DC变换器，电路图如图5所示。双向半桥DC-DC变换器将传统的buck电路和boost电路相结合，依靠开关管和二极管之间的切换实现能量的双向流动。该DC-DC变换器具有结构简单、易于控制、可靠性高等优点，也是目前应用最多的双向DC-DC变换器。
7.建立混合储能系统的成本函数。在马尔可夫决策过程算法中，成本函数对优化目标的实现效果起着非常重要的作用。Rloss代表混合储能系统中能量损耗的成本函数，其与混合储能系统的结构相关，其表达公式为：
式中Rbat、Ruc分别为电池和超级电容的内阻，Il为DC-DC电路中的电感电流，Dboost和Dbuck为电路的占空比，RSW为功率开关的导通电阻，RD为二极管的内阻，VD为二极管的导通压降，fsw＝50HZ为开关频率，tr、tf分别为功率开关在打开和关闭的过程中电流上升和下降的时间。
Rsoc代表SOCuc与参考值的差值，为了约束超级电容过度充放、保证超级电容能够随时吸收或者释放能量，SOCuc必须保持在某一合适的值附近，其表达公式为：
Rsoc＝ω2(SOCuc-SOCuc,ref)2
式中SOCuc,ref为SOCuc的参考值，Rsoc的功能是使SOCuc保持在SOCuc,ref附近波动。
降低成本函数Rloss，虽然会使混合储能系统的损失降低，但会导致超级电容过度充放，SOCuc因此会到达甚至超过SOC限制的极限。此外，由于超级电容会到达充满电的状态，制动能量也很可能无法被储存到超级电容。另一方面，降低成本函数SOCref，虽然可以使混合储能系统满足所有的功率需求，但会导致过高的系统能量损耗。因此，需要折中处理在两个优化目标，以达到最理想的效果。
8.通过优化算法求出最优控制策略。在得到状态S的转移概率矩阵Psoc,imj、目标函数R和折扣因子γ后，可以建立状态值函数Vπ(s)：
定义最优状态值函数V*(s)为在所有策略中值最大的值函数即：最优策略π*(a|S)是与最优状态值函数V*(s)对应的策略。利用状态值函数通过使用策略迭代法求解最优策略π*：
车速35Km/h时最优策略的计算结果如图6所示。决策过程示图如图7所示。
9.对锂电池功率Pbat的变化率进行限制。通过马尔可夫决策过程算法得到最优策略π*(a|S)后，由于动作集a是离散的，同时成本函数中没有对电池功率Pbat变化率进行限制，这导致对电池功率Pbat的控制结果变化过于频繁，产生一定的抖动，这对延长电池寿命和降低系统能量损耗十分不利。在通过最优策略π*(a|S)确定下一时刻动作a＝Pbat(t+1)时，加入对电池功率变化率ΔPbat的考虑。
式中a为功率值，b、c为SOCuc值，均为正数。ΔPbat＝Pbat(t+1)-Pbat(t)。
电池功率变化率的限制示意图如图8所示。
10.下面结合Matlab仿真实验对本发明的应用效果作详细的描述。
锂电池和超级电容的功率分配仿真结果如图9、10、11、12所示,可以发现,本发明提出的一种考虑功率变化率限制的混合储能系统能量管理优化方法，用电池来满足电动汽车的功率需求,而超级电容组用于处理瞬时峰值功率并实现制动能量回收,这大大减少了整个系统的能量损耗。另一方面,超级电容处理瞬时峰值功率,使锂电池的功率在理想的范围内平稳的变化，在一定程度上延长了锂电池的使用寿命。
从仿真结果可以得出:本发明提出的考虑功率变化率限制的混合储能系统能量管理优化方法可以从全局角度进行功率分配优化,从而获得良好的综合性能,特别是在经济性能方面,从能源和锂电池两方面降低了电动汽车的运行费用。
Claims (1)
1.一种考虑功率变化率限制的混合储能系统能量管理优化方法，包括如下的步骤：
(1)根据电动汽车的行驶工况数据以及电动汽车的基本参数，确定电动汽车总需求功率Pdem：
Pdem＝Pf+Pa+Pi
式中Pf是克服滚动阻力产生的功率、Pa是克服加速阻力产生的功率、Pi是克服空气阻力产生的功率、Pdem是电动汽车的需求功率；
(2)根据电动汽车的需求功率具有马尔可夫性，得到需求功率Pdem的状态转移概率矩阵，定义需求功率的状态转移概率Pdem,imj为：
式中，v为当前时刻电动汽车的速度，需求功率Pdem和电动汽车速度v被离散为有限值，即i,j＝1,2,…,Np；m＝1,2,…,Nv，当电动汽车的速度v＝vm时，汽车在t时刻的需求功率为的情况下，在t+1时刻的需求功率为的条件概率；
(3)根据锂电池和超级电容的工作特性以及实际的工况确定各参数的范围，各参数值限定范围为：
Pbat,min≤Pbat≤Pbat,max
Puc,min≤Puc≤Puc,max
SOCuc,min≤SOCuc≤SOCuc,max
式中Pbat为锂电池的功率、Puc为超级电容的功率、SOCuc为超级电容的荷电状态；
(4)已知当前时刻的状态St＝(v(t),Pdem(t),SOCuc(t))以及动作集a＝Pbat(t)，可以根据状态转移方程求出下一时刻的状态St+1＝(v(t+1),Pdem(t+1),SOCuc(t+1))，系统的状态转移方程为：
Puc(t)＝Pdem(t)-Pbat(t)
式中Iuc为超级电容的电流、Vuc为超级电容的端电压、Quc为超级电容的容量、Δt为采样时间间隔；
(5)根据需求功率的状态转移概率矩阵和系统的状态转移方程，可以求出状态S的转移概率矩阵Psoc,imj：
Psoc,imj＝P[St+1＝SOCuc,j|St＝SOCuc,i,At＝a,v＝vm]
公式的含义是当电动汽车的速度v＝vm、储能系统采取动作At＝a时，超级电容的SOC在t时刻的状态为SOCuc,i的情况下，在t+1时刻的状态为SOCuc,j的条件概率；
由需求功率Pdem和锂电池功率Pbat得到超级电容功率Puc，进而根据当前时刻的SOCuc得到下一时刻的SOCuc；
(6)以降低系统的能量损耗和维持超级电容SOC为目标，建立多目标的成本函数R，因为两个目标的优化相互矛盾，为实现均衡的优化目标，引入权重系数来实现多目标优化；
(7)已知状态S的转移概率矩阵、目标函数R和折扣因子γ，通过马尔可夫决策过程算法得到最优的策略π*(a|S)；
(8)在通过最优策略π*(a|S)确定下一时刻动作a＝Pbat(t+1)时，进一步限制锂电池功率变化率ΔPbat：当通过最优策略确定下一时刻锂电池功率时，需要与前一时刻的锂电池功率值进行比较，若功率变化率|ΔPbat|小于限定值a，则功率变化率ΔPbat无需受到限制，下一时刻锂电池功率Pbat(t+1)为有效的；若功率变化率|ΔPbat|大于限定值a，而此时SOCuc值在限定范围之内，即SOCuc,ref-b≤SOCuc≤SOCuc,ref+c，则功率变化率ΔPbat需要受到限制，下一时刻锂电池功率Pbat(t+1)＝Pbat(t)±a；若功率变化率|ΔPbat|大于限定值a，且此时SOCuc值超出限定范围，即SOCuc≤SOCuc,ref-b或SOCuc≥SOCuc,ref+c，则功率变化率ΔPbat受到的限制变小，下一时刻锂电池功率Pbat(t+1)＝Pbat(t)±a+0.2(ΔPbat±a)；其中a为功率值，b、c为SOCuc值，均为正数，ΔPbat＝Pbat(t+1)-Pbat(t)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810098996.2A CN108470220B (zh) | 2018-01-31 | 2018-01-31 | 考虑功率变化率限制的混合储能系统能量管理优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810098996.2A CN108470220B (zh) | 2018-01-31 | 2018-01-31 | 考虑功率变化率限制的混合储能系统能量管理优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108470220A true CN108470220A (zh) | 2018-08-31 |
| CN108470220B CN108470220B (zh) | 2021-11-30 |
Family
ID=63266313
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810098996.2A Expired - Fee Related CN108470220B (zh) | 2018-01-31 | 2018-01-31 | 考虑功率变化率限制的混合储能系统能量管理优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108470220B (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109774492A (zh) * | 2018-12-29 | 2019-05-21 | 江苏大学 | 一种基于未来驱动功率需求的纯电动汽车整车功率分配方法 |
| CN110466388A (zh) * | 2019-08-23 | 2019-11-19 | 哈尔滨理工大学 | 一种复合储能电动汽车能量管理方法 |
| CN110758121A (zh) * | 2019-11-13 | 2020-02-07 | 北京理工大学 | 一种基于递阶控制的能量管理系统 |
| CN110896246A (zh) * | 2019-12-05 | 2020-03-20 | 西南交通大学 | 一种混合储能式有轨电车系统的配置优化方法 |
| CN110962684A (zh) * | 2019-11-15 | 2020-04-07 | 东华大学 | 电动汽车能源管理与分配方法 |
| CN111091249A (zh) * | 2019-12-30 | 2020-05-01 | 吉林大学 | 一种基于全局寻域算法实现车辆全局能量最优分配方法 |
| CN111459025A (zh) * | 2020-03-26 | 2020-07-28 | 河南科技大学 | 一种运用粒子算法优化的复合电源电动汽车功率分配策略 |
| CN114997467A (zh) * | 2022-05-05 | 2022-09-02 | 国网江苏省电力有限公司镇江供电分公司 | 一种考虑状态转移模型的储能优化配置方法 |
| CN115954992A (zh) * | 2023-03-14 | 2023-04-11 | 银河航天(西安)科技有限公司 | 基于马尔科夫链的蓄电池过放保护方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140244247A1 (en) * | 2013-02-28 | 2014-08-28 | Google Inc. | Keyboard typing detection and suppression |
| US20150019149A1 (en) * | 2011-12-31 | 2015-01-15 | China Electric Power Research Institute | Real-time power distribution method and system for lithium battery and redox flow battery energy storage systems hybrid energy storage power station |
| CN105140942A (zh) * | 2015-10-09 | 2015-12-09 | 国家电网公司 | 考虑核电状态偏差的混合储能优化功率分配方法 |
| CN105416077A (zh) * | 2015-12-01 | 2016-03-23 | 北京理工大学 | 一种电动汽车的能量管理系统及管理方法 |
| CN105811385A (zh) * | 2016-04-27 | 2016-07-27 | 天津大学 | 半波长输电线路电流差动保护方法 |
| CN107294122A (zh) * | 2017-04-17 | 2017-10-24 | 国网浙江省电力公司电力科学研究院 | 一种混合储能系统分层动态控制方法 |
-
2018
- 2018-01-31 CN CN201810098996.2A patent/CN108470220B/zh not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150019149A1 (en) * | 2011-12-31 | 2015-01-15 | China Electric Power Research Institute | Real-time power distribution method and system for lithium battery and redox flow battery energy storage systems hybrid energy storage power station |
| US20140244247A1 (en) * | 2013-02-28 | 2014-08-28 | Google Inc. | Keyboard typing detection and suppression |
| CN105140942A (zh) * | 2015-10-09 | 2015-12-09 | 国家电网公司 | 考虑核电状态偏差的混合储能优化功率分配方法 |
| CN105416077A (zh) * | 2015-12-01 | 2016-03-23 | 北京理工大学 | 一种电动汽车的能量管理系统及管理方法 |
| CN105811385A (zh) * | 2016-04-27 | 2016-07-27 | 天津大学 | 半波长输电线路电流差动保护方法 |
| CN107294122A (zh) * | 2017-04-17 | 2017-10-24 | 国网浙江省电力公司电力科学研究院 | 一种混合储能系统分层动态控制方法 |
Non-Patent Citations (3)
| Title |
|---|
| AMIR HOSSEIN SHAHIRINIA等: "Robust Control of Hybrid Wind/Energy Storage Power Generation System Considering Statistical Wind Characteristics", 《2012 IEEE INTERNATIONAL CONFERENCE ON POWER AND ENERGY(PECON)》 * |
| 孙纯军等: "基于SOC状态反馈的混合储能功率优化策略", 《电测与仪表》 * |
| 朱江等: "基于隐马尔可夫预测的非对称信息功率博弈机制", 《计算机应用》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109774492A (zh) * | 2018-12-29 | 2019-05-21 | 江苏大学 | 一种基于未来驱动功率需求的纯电动汽车整车功率分配方法 |
| CN110466388B (zh) * | 2019-08-23 | 2020-11-03 | 哈尔滨理工大学 | 一种复合储能电动汽车能量管理方法 |
| CN110466388A (zh) * | 2019-08-23 | 2019-11-19 | 哈尔滨理工大学 | 一种复合储能电动汽车能量管理方法 |
| CN110758121A (zh) * | 2019-11-13 | 2020-02-07 | 北京理工大学 | 一种基于递阶控制的能量管理系统 |
| CN110758121B (zh) * | 2019-11-13 | 2021-06-04 | 北京理工大学 | 一种基于递阶控制的能量管理系统 |
| CN110962684A (zh) * | 2019-11-15 | 2020-04-07 | 东华大学 | 电动汽车能源管理与分配方法 |
| CN110896246B (zh) * | 2019-12-05 | 2022-04-26 | 西南交通大学 | 一种混合储能式有轨电车系统的配置优化方法 |
| CN110896246A (zh) * | 2019-12-05 | 2020-03-20 | 西南交通大学 | 一种混合储能式有轨电车系统的配置优化方法 |
| CN111091249A (zh) * | 2019-12-30 | 2020-05-01 | 吉林大学 | 一种基于全局寻域算法实现车辆全局能量最优分配方法 |
| CN111459025A (zh) * | 2020-03-26 | 2020-07-28 | 河南科技大学 | 一种运用粒子算法优化的复合电源电动汽车功率分配策略 |
| CN114997467A (zh) * | 2022-05-05 | 2022-09-02 | 国网江苏省电力有限公司镇江供电分公司 | 一种考虑状态转移模型的储能优化配置方法 |
| CN114997467B (zh) * | 2022-05-05 | 2023-07-25 | 国网江苏省电力有限公司镇江供电分公司 | 一种考虑状态转移模型的储能优化配置方法 |
| CN115954992A (zh) * | 2023-03-14 | 2023-04-11 | 银河航天(西安)科技有限公司 | 基于马尔科夫链的蓄电池过放保护方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108470220B (zh) | 2021-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108470220A (zh) | 考虑功率变化率限制的混合储能系统能量管理优化方法 | |
| Zhao et al. | Energy transfer and utilization efficiency of regenerative braking with hybrid energy storage system | |
| CN109747625B (zh) | 一种混合动力车辆复合式能量管理方法 | |
| CN111211554B (zh) | 城轨混合储能系统的功率动态分配控制方法 | |
| CN103490494B (zh) | 一种应用于混合动力汽车车载复合电源 | |
| CN109552110B (zh) | 基于规则与非线性预测控制的电动汽车复合能量管理方法 | |
| CN107229987A (zh) | 基于光储系统运行优化的混合储能配比计算方法 | |
| CN109849693B (zh) | 一种基于自适应小波转换的电动汽车混合能源管理系统及其控制方法 | |
| CN106494328B (zh) | 一种基于电功率在线计算的燃油汽车电能控制系统及方法 | |
| CN112421745A (zh) | 一种电动汽车复合电源系统能量管理方法 | |
| CN107415730B (zh) | 一种车用燃料电池电源系统的功率控制方法 | |
| CN108879730A (zh) | 混合储能系统及基于其的风电功率波动平抑方法 | |
| CN107846043B (zh) | 一种考虑电动汽车充电影响的微网能量管理方法 | |
| CN107069782B (zh) | 应用于轨道交通车载混合储能系统的容量配置方法 | |
| CN108767872A (zh) | 一种应用于风光混合储能微网系统的模糊控制方法 | |
| CN112060982A (zh) | 一种动态平衡的燃料电池无人机能量管理方法 | |
| Zhang et al. | Genetic Algorithm based optimal component sizing for an electric vehicle | |
| CN114030392A (zh) | 燃料电池-锂电池混源供电系统效率优化能量管理策略 | |
| Wang et al. | Hybrid energy management strategy based on dynamic setting and coordinated control for urban rail train with PMSM | |
| CN203747469U (zh) | 自适应滤波器功率分流控制的混合动力车复合电源 | |
| CN102570905B (zh) | 利用太阳能的汽车供电系统和汽车供电控制方法 | |
| CN203301211U (zh) | 基于改进型滤波器功率分流控制混合动力汽车用复合电源 | |
| Li et al. | Research on multi-energy management system of fuel cell vehicle based on fuzzy control | |
| CN106849053B (zh) | 一种车载复合电源功率分配滑膜变结构控制方法 | |
| Yanzi et al. | Energy management strategy based on fuzzy logic for a new hybrid battery-ultracapacitor energy storage system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20211130 |