1.一种模块化多电平电池储能系统电池混用的方法，其特征在于，包括以下步骤：
S1：获取模块化多电平电池储能系统各个子模块电池电压、SOC、SOH、SOF、额定容量信息；模块化多电平电池储能系统中，每个子模块包括电池单元和功率单元，电池单元由电池管理系统管理，功率单元作为功率转换系统的一部分由PCS控制器控制，PCS控制器定时从BMS获取各个功率单元对应的电池单元的SOC状态和SOH状态；
S2：计算各个子模块的可充电能量和可放电能量：根据步骤S1获取的子模块电池SOC、SOH以及电池额定容量信息，分别计算出各个子模块电池的可充电能量和可放电能量；
S3：计算各相和整个系统总的可充电电量和可放电电量；
S4：各相充放电功率的预分配：按照三相的可充电能量SOCEa，SOCEb，SOCEc的比例分配充电功率，按照三相的可放电能量SOCEa，SOCEb，SOCEc的比例分配放电功率；
S5：子模块功率的预分配：在各相功率预分配的基础上分配各个子模块的功率，各个子模块功率按照与子模块的可充/放电能量的成正比例的原则分配；
S6：子模块功率的校核；
S7：子模块功率和三相功率的再分配：根据电池电压和电池配合电流，重新分配各个子模块功率，并且重新分配三相功率；
S8：零序电压注入实现三相间均衡；
S9：子模块交流电压的分配；同相各个子模块交流电压相位相同，按照各子模块的可充/放电能量的比例分配交流电压，即控制了交流侧功率。
2.根据权利要求1所述的模块化多电平电池储能系统电池混用的方法，其特征在于，步骤S1中，PCS控制器定时从BMS获取各个功率单元对应的电池单元的SOC状态和SOH状态，时间间隔根据电池储能系统的状态刷新速率确定，取0.1s-1min。
3.根据权利要求1所述的模块化多电平电池储能系统电池混用的方法，其特征在于，步骤S2中，计算各个子模块的可充电能量和可放电能量，具体方法为：
可放电能量：
SODE(x,n)＝[(SOCx,n-SOCdown)×SOHx,n×CN]×VN
可充电能量：
SOCE(x,n)＝[(SOCup-SOCx,n)×SOHx,n×CN]×VN
式中，SOCup和SOCdown分别代表电池运行的SOC上下边界，0≤SOCdown4.根据权利要求1所述的模块化多电平电池储能系统电池混用的方法，其特征在于，步骤S3中，计算各个子模块的可充电能量和可放电能量，具体方法为：
计算每相的可最大放电能量：
式中，下标x表示a、b、c三相之一，n表示该相的第n个子模块，N为每相的子模块数；
计算三相总的最大可放电电量：
SODE＝SODEa+SODEb+SODEc
计算每相的充电能量：
式中，下标x表示a、b、c三相之一，n表示该相的第n个子模块，N为每相的子模块数；
计算三相总的可充电能量：
SOCE＝SOCEa+SOCEb+SOCEc。
5.根据权利要求4所述的模块化多电平电池储能系统电池混用的方法，其特征在于，步骤S4中，各相充放电功率的预分配具体包括：
总功率指令为Psum，放电时，abc三相分配的放电功率指令分别为：
充电时，abc三相的充电功率指令分别为：
6.根据权利要求5所述的模块化多电平电池储能系统电池混用的方法，其特征在于，步骤S5中，子模块功率的预分配，具体方法为：
充电时，abc三相各个子模块的功率分配如下：
放电时，abc三相各个子模块的功率分配如下：
式中，Pa,n，Pb,n，Pc,n分别表示a、b、c三相的第n个子模块的功率，下标a、b、c表示abc三相，n表示子模块的编号，N表示每相的子模块数量。
7.根据权利要求6所述的模块化多电平电池储能系统电池混用的方法，其特征在于，步骤S6中，子模块功率的校核具体方法为：
根据能量守恒，电池侧功率等于交流功率和直流功率之和，只有交流侧控制的情况下，直流侧功率为零；
三相中各个子模块电池中的电流平均值分别为：
根据步骤S1获取的各个子模块电池的SOF包含的可充电电流和可放电电流信息，校核上述电流平均值是否在SOF允许的可充电电流和可放电电流范围之内；
充电时，对于a相子模块的电池：
如Ibat,a,n≤Ichg,a,n,则保持Ibat,a,n不变；
如Ibat,a,n>Ichg,a,n,则重置Ibat,a,n＝Ichg,a,n；
同理，对于b相子模块的电池：
如Ibat,b,n≤Ichg,b,n,则保持Ibat,b,n不变；
如Ibat,b,n>Ichg,b,n,则重置Ibat,b,n＝Ichg,b,n；
对于c相子模块的电池：
如Ibat,c,n≤Ichg,c,n,则保持Ibat,c,n不变；
如Ibat,c,n>Ichg,c,n,则重置Ibat,c,n＝Ichg,c,n；
式中，Ibat,a,n，Ibat,b,n，Ibat,c,n分别表示a、b、c三相的第n个子模块的电池平均电流，下标a、b、c表示abc三相，Ichg,a,n，Ichg,b,n，Ichg,c,n分别表示a、b、c三相的第n个子模块的电池的允许充电电流，n表示子模块的编号,N表示每相的子模块数量；
放电时，对于a相子模块的电池：
如Ibat,a,n≤Idchg,a,n,则保持Ibat,a,n不变；
如Ibat,a,n>Idchg,a,n,则重置Ibat,a,n＝Idchg,a,n；
同理，对于b相子模块的电池：
如Ibat,b,n≤Idchg,b,n,则保持Ibat,b,n不变；
如Ibat,b,n>Idchg,b,n,则重置Ibat,b,n＝Idchg,b,n；
对于c相子模块的电池：
如Ibat,c,n≤Idchg,c,n,则保持Ibat,c,n不变；
如Ibat,c,n>Idchg,c,n,则重置Ibat,c,n＝Idchg,c,n；
式中，Ibat,a,n，Ibat,b,n，Ibat,c,n分别表示a、b、c三相的第n个子模块的电池平均电流，下标a、b、c表示abc三相，Idchg,a,n，Idchg,b,n，Idchg,c,n分别表示a、b、c三相的第n个子模块的电池的允许充电电流，n表示子模块的编号,N表示每相的子模块数量。
8.根据权利要求7所述的模块化多电平电池储能系统电池混用的方法，其特征在于，步骤S7中，子模块功率和三相功率的再分配的具体方法为：
根据电池电压和电池配合电流，重新分配各个子模块功率如下：
Pa,n＝Ubat,a,n×Ibat,a,n
Pb,n＝Ubat,b,n×Ibat,b,n
Pc,n＝Ubat,c,n×Ibat,c,n
重新分配三相功率如下：
9.根据权利要求8所述的模块化多电平电池储能系统电池混用的方法，其特征在于，步骤S8中，零序电压注入实现三相间均衡，具体方法为：
abc三相的功率与平均值的偏差为：
三相功率的偏差通过注入零序电压的方式来实现；
注入的零序电压的相对A相正序电流相位为：
零序电压的大小为：
A相正序电流的相位通过对电网电压锁相得到；
a、b、c三相电压相量分别为：
三相电流相量及大小分别为：
式中，Us为电网系统电压有效值。
10.根据权利要求9所述的模块化多电平电池储能系统电池混用的方法，其特征在于，步骤S9中，子模块交流电压的分配，具体方法为：
同相各个子模块交流电压相位相同，按照各子模块的可充/放电能量的比例分配交流电压，即控制了交流侧功率；
充电时，abc三相各个子模块交流电压分配如下：
放电时，abc三相各个子模块交流电压分配如下：
式中，Ua,n，Ub,n，Uc,n分别表示a、b、c三相的第n个子模块的交流侧电压，下标a、b、c表示abc三相，n表示子模块的编号，N表示每相的子模块数量。
技术领域
[0001] 本
发明涉及电池储能系统领域，更具体地，涉及一种模块化多电平电池储能系统电池混用的方法。
背景技术
[0002] 模块化多电平电池储能系统以其具有可同时连接交直流
电网的端口、模块化的结构、谐波小的特点，被认为是最有前景的电池储能系统拓扑结构之一。该储能系统模块化的结构使其对于接入的电池具有差异化控制的能
力，可用于实现不同电池的混用。
[0003] 目前国内外对同一个电池储能系统内不同电池混用的研究尚未见报道。相近的研究是对电池的梯次利用。张大宝等在《变电站智能并联
蓄电池系统新技术研究与应用[J].科技与创新，2018(09)：148-149.》中提出通过DC/DC或DC/AC变换之后在进行并联，从而实现新旧电池的混用。姚克宇在《关于存量电池利旧和新旧电池混用模式的研究[J].移动通信，2016，40(15)：64-67.》中采用电池合路器实现新旧电池的混用，电池合路器的本质是多路切换
开关。上述研究和应用虽名称叫混用，但只是同类电池的混用，而且也仅仅限于移动基站备用电源领域。
[0004] 对于电池储能领域的梯次利用，目前有两个技术路线。第一，通过增强
单体间的均衡适应梯次电池应用的需求。如赵光金等在《主被动均衡技术及其在电池梯次利用中的应用[J].电源技术，2018，42(07)：983-986+1075.》中采用主动和被动组合实现梯次电池的均衡，主动和被动均衡之间分时切换；朱运征《集装箱式储能系统用梯次利用锂
电池组的一致性管理研究[J].电源学报，2018，16(04)：80-86.》和李晓依《锂离子动力电池梯次利用的主动均衡系统研究与实现[D].武汉理工大学，2018.》均通过增加
电池管理系统的主动均衡能力来消减梯次电池的差异性的影响。第二，通过电力
电子模块实现对电池模块的差异化控制，适应梯次电池应用的需求。如李丹在《模块化独立控制梯次利用电池储能系统[D].北京交通大学，2018.》中利用级联H桥拓扑结构的相间和相内均衡功能实现梯次电池的利用；李圆圆等在《快速收敛的梯次电池成组技术[J].高
电压技术，2018，44(01)：169-175.》中将电池模组通过DC/DC变换后
串联实现梯次电池的应用。
[0005] 以梯次电池利用为目标的控制方法，其目标均是使均衡功率与SOC的偏差值成正比相关。但这种以消除SOC偏差为直接目标的控制策略对于储能系统中存在不同电池混用时就存在一些不足：1)同一储能系统中不同电池的容量可能存在较为显著的差异，同样SOC的电池却存储着不同大小的电量，存储着相同电量的电池处于不同的SOC状态。2)同一储能系统中不同电池的平台电压可能存在较为显著的差异，同样SOC和容量的电池却存储着不同大小的
能量，存储着相同能量的电池具有不同的容量和处于不同的SOC状态。单纯以SOC和SOC偏差进行的控制在不同电池混用的储能系统中将不再适用，不仅不能取得好的效果，甚至可能适得其反。
发明内容
[0006] 本发明针对
现有技术存在的不足，提出一种模块化多电平电池储能系统电池混用的方法。该方法以在模块化多电平电池储能系统中实现不同电池的混用为目标，利用交流侧的功率控制实现对容量和平台电压不同电池模块的差异化控制，实现不同电池模块的安全运行，起到最大化电池利用的目的。
[0007] 为解决上述技术问题，本发明的技术方案如下：
[0008] 一种模块化多电平电池储能系统电池混用的方法，包括以下步骤：
[0009] S1：获取模块化多电平电池储能系统各个子模块电池电压、SOC(state of charge，电池剩余电量百分比)、SOH(state of health，电池健康度)、SOF(state of
function，电池的功能状态)、额定容量信息；模块化多电平电池储能系统中，每个子模块包括电池单元和功率单元，电池单元由电池管理系统(Battery Management System，BMS)管理，功率单元作为功率转换系统(Power Conversion System，PCS)的一部分由PCS
控制器控制，PCS控制器定时从BMS获取各个功率单元对应的电池单元的SOC状态和SOH状态；
[0010] S2：计算各个子模块的可充
电能量和可放电能量：根据步骤S1获取的子模块电池SOC、SOH以及电池额定容量信息，分别计算出各个子模块电池的可充电能量和可放电能量；
[0011] S3：计算各相和整个系统总的可充电电量和可放电电量；
[0012] S4：各相充放电功率的预分配：按照三相的可充电能量SOCEa，SOCEb，SOCEc的比例分配充电功率，按照三相的可放电能量SOCEa，SOCEb，SOCEc的比例分配放电功率；
[0013] S5：子模块功率的预分配：在各相功率预分配的
基础上分配各个子模块的功率，各个子模块功率按照与子模块的可充/放电能量的成正比例的原则分配；
[0014] S6：子模块功率的校核；
[0015] S7：子模块功率和三相功率的再分配：根据电池电压和电池配合
电流，重新分配各个子模块功率，并且重新分配三相功率；
[0016] S8：零序电压注入实现三相间均衡；
[0017] S9：子模块交流电压的分配；同相各个子模块交流电压
相位相同，按照各子模块的可充/放电能量的比例分配交流电压，即控制了交流侧功率。
[0018] 在其中一个
实施例中，步骤S1中，PCS控制器定时从BMS获取各个功率单元对应的电池单元的SOC状态和SOH状态，时间间隔根据电池储能系统的状态刷新速率确定，取0.1s-1min。获取方式通常通讯，具体由PCS和BMS之间的
接口规范和协议确定。
[0019] 在其中一个实施例中，步骤S2中，计算各个子模块的可充电能量和可放电能量，具体方法为：
[0020] 可放电能量：
[0021] SODE(x，n)＝[(SOCx，n-SOCdown)×SOHx，n×CN]×VN
[0022] 可充电能量：
[0023] SOCE(x，n)＝[(SOCup-SOCx，n)×SOHx，n×CN]×VN
[0024] 式中，SOCup和SOCdown分别代表电池运行的SOC上下边界，0≤SOCdown＜SOCup≤1，x表示a、b、c三相之一，n表示某一相中的子模块编号，CN为电池额定容量，VN为电池标称电压。
[0025] 在其中一个实施例中，步骤S3中，计算各个子模块的可充电能量和可放电能量，具体方法为：
[0026] 计算每相的可最大放电能量：
[0027]
[0028] 式中，下标x表示a、b、c三相之一，n表示该相的第n个子模块，N为每相的子模块数；
[0029] 计算三相总的最大可放电电量：
[0030] SODE＝SODEa+SODEb+SODEc
[0031] 计算每相的充电能量：
[0032]
[0033] 式中，下标x表示a、b、c三相之一，n表示该相的第n个子模块，N为每相的子模块数；
[0034] 计算三相总的可充电能量：
[0035] SOCE＝SOCEa+SOCEb+SOCEc。
[0036] 在其中一个实施例中，步骤S4中，各相充放电功率的预分配具体包括：
[0037] 总功率指令为Psum，放电时，abc三相分配的放电功率指令分别为：
[0038]
[0039]
[0040]
[0041] 充电时，abc三相的充电功率指令分别为：
[0042]
[0043]
[0044]
[0045] 在其中一个实施例中，步骤S5中，子模块功率的预分配，具体方法为：
[0046] 充电时，abc三相各个子模块的功率分配如下：
[0047]
[0048]
[0049]
[0050] 放电时，abc三相各个子模块的功率分配如下：
[0051]
[0052]
[0053]
[0054] 式中，Pa，n，Pb，n，Pc，n分别表示a、b、c三相的第n个子模块的功率，下标a、b、c表示abc三相，n表示子模块的编号，N表示每相的子模块数量。
[0055] 在其中一个实施例中，步骤S6中，子模块功率的校核具体方法为：
[0056] 根据能量守恒，电池侧功率等于交流功率和直流功率之和，只有交流侧控制的情况下，直流侧功率为零；
[0057] 三相中各个子模块电池中的电流平均值分别为：
[0058]
[0059]
[0060]
[0061] 根据步骤S1获取的各个子模块电池的SOF包含的可充电电流和可放电电流信息，校核上述电流平均值是否在SOF允许的可充电电流和可放电电流范围之内；
[0062] 充电时，对于a相子模块的电池：
[0063] 如Ibat，a，n≤Ichg，a，n，则保持Ibat，a，n不变；
[0064] 如Ibat，a，n＞Ichg，a，n，则重置Ibat，a，n＝Ichg，a，n；
[0065] 同理，对于b相子模块的电池：
[0066] 如Ibat，b，n≤Ichg，b，n，则保持Ibat，b，n不变；
[0067] 如Ibat，b，n＞Ichg，b，n，则重置Ibat，b，n＝Ichg，b，n；
[0068] 对于c相子模块的电池：
[0069] 如Ibat，c，n≤Ichg，c，n，则保持Ibat，c，n不变；
[0070] 如Ibat，c，n＞Ichg，c，n，则重置Ibat，c，n＝Ichg，c，n；
[0071] 式中，Ibat，a，n，Ibat，b，n，Ibat，c，n分别表示a、b、c三相的第n个子模块的电池平均电流，下标a、b、c表示abc三相，Ichg，a，n，Ichg，b，n，Ichg，c，n分别表示a、b、c三相的第n个子模块的电池的允许充电电流，n表示子模块的编号，N表示每相的子模块数量；
[0072] 放电时，对于a相子模块的电池：
[0073] 如Ibat，a，n≤Idchg，a，n，则保持Ibat，a，n不变；
[0074] 如Ibat，a，n＞Idchg，a，n，则重置Ibat，a，n＝Idchg，a，n；
[0075] 同理，对于b相子模块的电池：
[0076] 如Ibat，b，n≤Idchg，b，n，则保持Ibat，b，n不变；
[0077] 如Ibat，b，n＞Idchg，b，n，则重置Ibat，b，n＝Idchg，b，n；
[0078] 对于c相子模块的电池：
[0079] 如Ibat，c，n≤Idchg，c，n，则保持Ibat，c，n不变；
[0080] 如Ibat，c，n＞Idchg，c，n，则重置Ibat，c，n＝Idchg，c，n；
[0081] 式中，Ibat，a，n，Ibat，b，n，Ibat，c，n分别表示a、b、c三相的第n个子模块的电池平均电流，下标a、b、c表示abc三相，Idchg，a，n，Idchg，b，n，Idchg，c，n分别表示a、b、c三相的第n个子模块的电池的允许充电电流，n表示子模块的编号，N表示每相的子模块数量。
[0082] 在其中一个实施例中，步骤S7中，子模块功率和三相功率的再分配的具体方法为：
[0083] 根据电池电压和电池配合电流，重新分配各个子模块功率如下：
[0084] Pa，n＝Ubat，a，n×Ibat，a，n
[0085] Pb，n＝Ubat，b，n×Ibat，b，n
[0086] Pc，n＝Ubat，c，n×Ibat，c，n
[0087] 重新分配三相功率如下：
[0088]
[0089]
[0090]
[0091] 在其中一个实施例中，步骤S8中，零序电压注入实现三相间均衡，具体方法为：
[0092] abc三相的功率与平均值的偏差为：
[0093]
[0094]
[0095]
[0096] 三相功率的偏差通过注入零序电压的方式来实现；
[0097] 注入的零序电压的相对A相正序电流相位为：
[0098]
[0099] 零序电压的大小为：
[0100]
[0101] A相正序电流的相位通过对电网电压
锁相得到；
[0103]
[0104]
[0105]
[0107]
[0108]
[0109]
[0110]
[0111] 式中，Us为电网系统电压有效值。
[0112] 在其中一个实施例中，步骤S9中，子模块交流电压的分配，具体方法为：
[0113] 同相各个子模块交流电压相位相同，按照各子模块的可充/放电能量的比例分配交流电压，即控制了交流侧功率；
[0114] 充电时，abc三相各个子模块交流电压分配如下：
[0115]
[0116]
[0117]
[0118] 放电时，abc三相各个子模块交流电压分配如下：
[0119]
[0120]
[0121]
[0122] 式中，Ua，n，Ub，n，Uc，n分别表示a、b、c三相的第n个子模块的交流侧电压，下标a、b、c表示abc三相，n表示子模块的编号，N表示每相的子模块数量。
[0123] 与现有技术相比，本发明技术方案的有益效果是：本发明提供了一种模块化多电平电池储能系统电池混用的方法，考虑了子模块中不同电池的标称电压的不同、
健康状态的不同、容量的不同、可充放电能量的不同，实现了对不同电池混用时能量的充分利用；考虑到电池的SOF，体现了不同电池充放电电流的约束，电池的寿命和安全性更有保证。最终达到在模块化多电平电池储能系统中安全、合理、充分利用不同电池的目的。
附图说明
[0124] 图1为实施例1模块化多电平电池储能系统示意图。
[0125] 图2为实施例1模块化多电平电池储能系统电池混用的方法的
流程图。
具体实施方式
[0126] 下面结合附图和实施例对本发明的技术方案做进一步的说明。
[0127] 实施例1
[0128] 如图1所示，本实施例为150kW电池储能系统，交流侧额定电压380V，
星形连接，每个桥臂有N＝10个子模块，整个系统共60个子模块。桥臂电抗0.6mH，交流并网电抗0.3mH。电池储能系统收到100kW放电指令。
[0129] 本实施例中，其中57个子模块电池为由32节3.2V/100Ah
磷酸铁锂电池单体串联组成标称电压102.4V，标称容量40Ah的储能电池模块。另有2个子模块电池为42节2.3V/70A
钛酸锂电池单体串联组成的标称电压96.6V，标称容量30A的储能电池模块，位于A相和C相最后。有1个子模块电池为42节2.3V/70A钛酸锂电池单体串联组成的标称电压96.6V，标称容量35A的储能电池模块，位于B相最后。电池SOC的运行上限均设为0.9，运行下限均设为0.1。
[0130] 如图2所示，本实施例的过程如下：
[0131] S1：获取模块化多电平储能系统各个子模块电池电压、SOC、SOH、SOF、额定容量信息
[0132] 能量转换系统通过通讯方式定时每1s从电池管理系统获取三相共120个子模块电池的信息。设置SOC运行上下限分别为0.9和0.1。获取的信息如下：
[0133] A相信息：
[0134]
荷电状态SOCa＝[0.59，0.65，0.65，0.56，0.61，0.56，0.63，0.60，0.57，0.61，0.62，0.62，0.60，0.55，0.58，0.56，0.57，0.60，0.59，0.64]
[0135] 健康状态SOHa＝[0.91，0.93，0.92，0.93，0.90，0.91，0.90，0.92，0.93，0.92，0.91，0.95，0.95，0.92，0.94，0.91，0.93，0.90，0.91，0.92]
[0136] 电池实测电压Ua＝[102.3，101.7，100.6，97.1，98.6，103.4，97.7，97.2，96.1，100.7，97.8，96.4，100.6，103.3，103.2，103.9，102.5，102.6，97.9，96.1]，单位V；
[0137] 电池标称电压UN＝[102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，96.6]，单位V；
[0138] 电池标称容量CN＝[40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，30]，单位AH；
[0139] 可放电电流Idchg＝[40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，20]，单位A；
[0140] 可充电电流Ichg＝[40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，20]，单位A。
[0141] B相信息：
[0142] 荷电状态SOCb＝[0.63，0.58，0.65，0.62，0.64，0.58，0.60，0.57，0.63，0.62，0.57，0.56，0.64，0.60，0.62，0.56，0.58，0.58，0.62，0.63]；
[0143] 健康状态SOHb＝[0.92，0.95，0.92，0.90，0.94，0.95，0.92，0.90，0.90，0.94，0.93，0.94，0.91，0.95，0.95，0.94，0.90，0.90，0.90，0.93]；
[0144] 电池实测电压Ub＝[100.9，100.4，98.4，97.4，96.7，98.8，98.0，102.2，98.8，103.6，98.6，96.2，104.0，97.7，103.7，101.3，99.9，103.2，102.1，97.5]，单位V；
[0145] 电池标称电压UN＝[102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，96.6]，单位V；
[0146] 电池标称容量CN＝[40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，35]，单位AH；
[0147] 可放电电流Idchg＝[40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，20]，单位A；
[0148] 可充电电流Ichg＝[40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，25]，单位A。
[0149] C相信息：
[0150] 荷电状态SOCc＝[0.56，0.61，0.61，0.63，0.63，0.63，0.60，0.61，0.65，0.63，0.60，0.63，0.61，0.59，0.65，0.62，0.57，0.56，0.55，0.64]；
[0151] 健康状态SOHc＝[0.94，0.94，0.90，0.94，0.92，0.93，0.95，0.91，0.92，0.95，0.93，0.90，0.92，0.95，0.94，0.94，0.93，0.94，0.93，0.90]；
[0152] 电池实测电压Uc＝[102.6，103.4，98.7，102.8，96.3，96.7，100.0，102.9，101.0，102.9，99.0，96.6，96.4，99.7，96.6，98.0，98.6，103.9，96.5，95.9]，单位V；
[0153] 电池标称电压UN＝[102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，102.4，96.6]，单位V；
[0154] 电池标称容量CN＝[40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，30]，单位AH；
[0155] 可放电电流Idchg＝[40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，20]，单位A；
[0156] 可充电电流Ichg＝[40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，40，20]，单位A。
[0157] S2：计算各个子模块的可充电能量和可放电能量
[0158] 根据步骤一获取的信息，分别计算出各个子模块的可充电能量和可放电能量。
[0159] A相：
[0160] 可放电能量SODE＝[2.22，2.47，2.45，2.14，2.23，2.10，2.33，2.25，2.18，2.28，2.31，2.43，2.35，2.08，2.22，2.10，2.16，2.20，2.19，1.70]，单位kWh；
[0161] 可充电能量SOCE＝[1.51，1.34，1.32，1.67，1.46，1.63，1.36，1.52，1.63，1.49，1.42，1.46，1.55，1.69，1.63，1.63，1.65，1.48，1.54，0.96]，单位kWh。
[0162] B相：
[0163] 可放电能量SODE＝[2.37，2.27，2.43，2.30，2.46，2.25，2.25，2.12，2.34，2.38，2.17，2.17，2.38，2.33，2.40，2.15，2.13，2.15，2.27，1.98]，单位kWh；
[0164] 可充电能量SOCE＝[1.40，1.62，1.34，1.39，1.39，1.64，1.52，1.57，1.35，1.47，1.64，1.68，1.35，1.56，1.49，1.70，1.56，1.53，1.42，1.16]，单位kWh。
[0165] C相：
[0166] 可放电能量SODE＝[2.16，2.34，2.25，2.44，2.37，2.42，2.33，2.27，2.44，2.44，2.30，2.32，2.32，2.28，2.50，2.38，2.17，2.14，2.10，1.66]，单位kWh；
[0167] 可充电能量SOCE＝[1.69，1.51，1.44，1.41，1.39，1.39，1.56，1.46，1.33，1.45，1.50，1.36，1.45，1.61，1.35，1.47，1.64，1.71，1.71，0.94]，单位kWh。
[0168] S3：计算各相和整个系统总的可充电电量和可放电电量
[0169] 根据各个子模块充放电能量，求和得到：
[0170] A相总的可放电能量SODE＝44.39kWh；总的可充电能量SOCE＝29.92kWh；
[0171] B相总的可放电能量SODE＝45.29kWh；总的可充电能量SOCE＝29.78kWh；
[0172] C相总的可放电能量SODE＝45.63kWh；总的可充电能量SOCE＝29.40kWh；
[0173] 电池储能系统总的可放电能量135.30kWh，总的可充电能量89.10kWh。
[0174] S4：各相充放电功率的预分配
[0175] 根据三相的可放电容量，预分配三相功率如下：A相功率Pa为32.81kW；B相功率Pb为33.47kW；C相功率Pc为33.72kW。
[0176] S5：子模块功率的预分配
[0177] 预分配各个子模块功率如下：
[0178] Pa＝[1.64，1.83，1.81，1.58，1.65，1.55，1.72，1.66，1.61，1.69，1.71，1.79，1.73，1.54，1.64，1.55，1.60，1.63，1.62，1.26]，单位kW；
[0179] Pb＝[1.75，1.68，1.80，1.70，1.82，1.66，1.66，1.56，1.73，1.76，1.60，1.60，1.76，1.72，1.78，1.59，1.57，1.59，1.68，1.46]，单位kW；
[0180] Pc＝[1.59，1.73，1.66，1.80，1.76，1.79，1.72，1.67，1.80，1.80，1.70，1.72，1.71，1.69，1.85，1.76，1.60，1.58，1.55，1.23]，单位kW。
[0181] 根据预分配功率和电池电压，计算得到电池平均电流：
[0182] Ia＝[16.02，17.95，17.98，16.29，16.71，15.03，17.61，17.11，16.76，16.74，17.46，18.62，17.23，14.91，15.88，14.92，15.60，15.87，16.52，13.10]，单位A；
[0183] Ib＝[17.37，16.69，18.27，17.44，18.80，16.83，16.97，15.30，17.49，16.95，16.23，16.66，16.91，17.60，17.13，15.67，15.75，15.42，16.43，15.01]，单位A；
[0184] Ic＝[15.53，16.70，16.82，17.52，18.22，18.45，17.22，16.28，17.86，17.51，17.21，17.77，17.75，16.93，19.13，17.96，16.24，15.24，16.11，12.83]，单位A。
[0185] S6：子模块功率的校核
[0186] 与第一步得到的电池可放电电流校核发现，预分配电流小于电池可放电电流。电池储能系统各个子模块可以按照预分配功率运行。
[0187] S7：子模块功率和三相功率的再分配
[0188] 子模块可以按照预分配功率运行，则三相均可按照预分配功率运行。
[0189] S8：零序电压注入实现三相间均衡
[0190] 三相功率与三相功率平均值的偏差：
[0191] ΔPa＝Pa-Psum/3＝33.53kW-33.33kW＝0.20kW；
[0192] ΔPb＝Pb-Psum/3＝33.09kW-33.33kW＝-0.24kW；
[0193] ΔPc＝Pc-Psum/3＝33.39kW-33.33kW＝0.06kW；
[0194] 三相电流Ia＝Ib＝Ic＝100kW/380V/1.732＝152A
[0195] 三相功率的偏差通过注入零序电压的方式来实现。
[0196] 注入的零序电压的相对A相正序电流相位为(以A相电压相位为参考)：
[0197]
[0198]
[0199] 零序电压的大小为：
[0200]
[0201] 故零序电压相位为38.9°，幅度为1.69V。
[0202] 本实例中，根据系统描述，计算得到交流等效电感为0.6mH，100kW时电感上工频电压为28.6V。由此，注入零序电压得到的
中性点偏移的三相电压为：
[0203]
[0204]
[0205]
[0206] S9：子模块交流电压的分配
[0207] 在三相电压确定的基础上分配各个子模块的电压，按照各子模块的可放电能量的比例分配子模块电压即实现了按此比例分配功率。
[0208] 三相各个子模块电压：
[0209] Ua＝[10.98，10.81，12.47，11.75，12.21，11.40，11.87，10.81，11.55，10.78，10.62，11.03，10.84，12.38，11.53，11.12，11.54，10.42，10.77，8.41]，单位V；
[0210] Ub＝[11.28，10.77，11.63，10.73，10.39，10.69，11.60，10.85，11.04，11.23，11.31，11.33，11.06，11.46，10.72，11.30，11.44，10.98，11.42，9.15]，单位V；
[0211] Uc＝[10.97，11.30，10.40，11.55，10.14，10.80，11.63，12.11，12.03，11.74，11.24，11.09，10.26，11.25，11.90，11.61，11.90，11.25，10.84，7.89]，单位V。
[0212] 控制各个子模块的交流侧电压比例即实现了按可放电容量比例控制各个子模块电池的放电功率。
[0213] 尽管本发明的内容已经通过上述优选实施例作了详细介绍，但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后，对于本发明的多种
修改和替代都将是显而易见的。因此，本发明的保护范围应由所附的
权利要求来限定。
