1.一种含电池储能系统的主动配电网短期有功优化方法，其特征在于，包括以下步骤：
步骤1，建立电池储能系统有功及能量的模型，具体包括以下步骤：
步骤1-a，采集电池储能系统的基本参数，包括：电池储能系统额定有功功率Prated，最大的视在功率SBESS,max，充放电效率ηin，放电效率ηout，额定容量Wrated；设置工作电池储能系统的参数：电池荷电状态最大和最小值SOCmax、SOCmin；
步骤1-b，建立电池储能系统有功输入输出模型：
其中，It为第t时段电池储能系统的充放电状态，It＝1表示放电，It＝-1表示充电，为第t时段电池储能系统的有功出力， 和 分别表示第t时段内电池储能
系统的充电和放电有功功率；
步骤1-c，建立电池储能系统的能量模型：
其中， 为第t时段电池储能系统存储的能量， 为电池储能系统初始阶段存
储的能量，Δt为时间步长；
步骤2，建立电池储能系统等效运行成本模型；
电池储能系统等效运行成本包含固定投资成本和维护成本两部分；
设电池储能系统的固定投资成本均摊到每一次电池储能系统的循环充放电上，定义第k号电池储能系统的第m次循环周期的放电深度为DoD,k,m，所以电池储能系统的单次循环充放电的等效投资成本可表示为：
其中，CP和CW分别为电池储能系统单位功率投资成本和单位容量投资成本，其单位分别为人民币/kW，人民币/kW·h；Nctf(DoD,k,m)为放电深度DoD,k,m对应的电池储能系统循环寿命；
电池储能系统的维护成本与其额定有功功率Prated、额定容量Wrated以及工作时间有关，表示为：
CBESS＝Y(mPPrated+mWWrated)
其中，mP和mW分别为电池储能系统单位时间内单位功率和单位容量运行维护成本，其单位分别为人民币/kW·h，人民币/kW·h2；Y为工作时间；
步骤3，建立可控负荷的有功及能量关系模型，将可控负荷看作负出力的发电机：
步骤3-a，建立可直接控制负荷的有功及能量关系模型：
其中 为第t时段内可直接控制负荷的有功功率，为负值；EDMFD为调度周期内负荷总的电量需求；T为调度周期的总时段数；
步骤3-b，建立存储类可控负荷的有功及能量关系模型：
其中， 为第t时段存储类可控负荷中存储的能量， 为存储类可控负荷的
初始阶段存储的能量； 为存储类可控负荷的容量；μ为电能与存储能量之间的转换效率； 为各时段存储类可控负荷的有功注入，为负值； 为各时段从存储类可控负荷中提取热量的功率；
步骤4，在上述电池储能系统的有功及能量的模型、等效运行成本模型、可控负荷有功及能量关系模型的基础上，建立以配电网侧运行总成本最小为目标函数的短期有功优化模型；
步骤4-a，目标函数包含购售电成本、电池储能系统等效运行成本、弃风弃光成本、微型燃气轮机运行成本以及调整可控负荷的补偿成本，表示为：
min f＝C1+C2+C3+C4+C5
其中，C1为配电侧购电成本和售电收益之和，C2为电池储能系统等效运行成本，C3为优化调度过程中的弃风和弃光成本，C4为微型燃气轮机的运行成本，C5为调整可控负荷的补偿成本；
步骤4-b，建立等式约束，包含各时段节点潮流平衡方程约束、电池储能系统和可控负荷时段间耦合等式约束；
节点潮流平衡方程约束为：
其中： 和 为第t时段节点i处风电机组有功和无功出力； 和 为第t时段节
点i处光伏有功和无功出力； 和 则表示第t时段配网根节点处的有功和无功功率；
和 则表示第t时段节点i处电池储能系统的充放电有功和无功出力； 和
为第t时段节点i处可直接控制负荷有功和无功功率； 和 为第t时段节点i
处的存储类可控负荷有功和无功功率； 和 分别为第t时段i节点处微型燃气轮机有功和无功出力； 和 为第t时段节点i处固定负荷有功和无功功率；n为节点数目，Vit和为第t时段节点i的电压幅值，Gij和Bij为节点i和j之间互导纳的实部和虚部，
为第t时段i、j两节点电压的相角差；其中，风电机组和光伏的实际有功出力还应满足以下等式：
式中， 和 为第t时段节点i处风电机组和光伏的有功出力预测值； 和 分
别为第t时段节点i处的风机弃风或者光伏弃光的有功功率；
电池储能系统和存储类可控负荷能量约束分别为：
T 0
SOC＝SOC
其中，T为调度周期的总时段数；SOCT为电池储能系统最后时刻的剩余电量； 为存储类负荷最后时刻的剩余电量；
步骤4-c，建立不等式约束，包括节点电压上下限约束、支路潮流约束以及控制变量和状态变量相关约束：
节点电压上下限约束：
Vi,min＜Vit＜Vi,max
其中，Vi,min和Vi,max分别为节点i处电压允许的最小、最大值；
支路潮流约束：
其中， 为第t时段支路l的视在功率，Sl.max为支路l的视在功率最大值；
电池储能系统有功功率、视在功率和能量约束分别为：
SOCmin＜SOCt＜SOCmax
其中，SOCmax和SOCmin为电池储能系统的电池荷电状态最大和最小值， 表示第t时段电池储能系统的无功出力；
微型燃气轮机有功、无功出力限制以及爬坡约束：
其中， 为第t时段微型燃气轮机有功出力，PMT,i,min和PMT,i,max分别为微型燃气轮机有功出力的下限和上限； 为第t时段微型燃气轮机无功出力，QMT,min和QMT,max分别为微型燃气轮机无功出力的下限和上限；RMT,up和RMT,down分别表示微型燃气轮机增、减负荷时的爬坡功率上、下限；
不固定出力风电机组和光伏的有功、无功功率约束：
其中， 为第t时段风电机组的无功出力；QW,max和QW,min分别为风电机组无功出力的最大、最小值； 为第t时段光伏的无功出力，QPV,max和QPV,min分别为光伏的无功出力的最大、最小值；
可控负荷有功功率及能量约束：
其中， 为可控负荷的出力上限， 为可控负荷的出力下限；SOCSMFD,max和
SOCSMFD,min分别为存储类可控负荷荷电状态的上下限；
步骤5，考虑到上述有功优化问题为一个混合整数非线性规划问题，采用分支定界-原对偶内点法进行求解：
步骤5-a，根据风电出力和负荷的预测结果，按前一日调度情况求解各时段有功优化，得到各时段目标值，即当前时段对应目标函数的值，并将其作为分支定界法处理离散变量过程中各时段的上界；
步骤5-b，各电池储能系统在各时段内的充放电状态为离散变量，首先对所有离散变量进行松弛，然后采用原对偶内点法求解仅含连续变量的最优潮流问题，得到各控制变量、状态变量和对应的目标函数的值，判断所有离散变量结果是否为整数，若是，保存结果；否则，将该松弛问题及其目标函数值加入待分支队列，记为RP，同时将该目标函数值作为新的目标值下界；
步骤5-c，依次对RP中的子问题单个离散变量进行分支，采用原对偶内点法求解各松弛子问题，判断是否所有离散变量均取得整数值，若是，则保存结果，转5-f；若否，则将其加入下次待分支队列中，记为RRP；
步骤5-d，根据剪支准则对RRP中所有子问题进行剪支，并更新RP；
步骤5-e，判断待分支队列RP是否为空，若是，则出错，退出计算；若否，则转步骤5-c；
步骤5-f，对所有满足约束条件的调度方式进行放电深度的计算，然后得到等效运行成本，并加到每种调度方式的目标函数值中，并从中取出目标函数值最小的解作为最优解，实现主动配电网短期有功优化。
2.根据权利要求1所述的一种含电池储能系统的主动配电网短期有功优化方法，其特征在于：所述步骤2中放电深度与其循环寿命之间数学关系的建立以及放电深度的计算步骤为：
步骤2-a，储能设备的寿命与温度、峰值电流、其运行过程中的放电深度DoD密切相关，放电深度越大，储能设备的循环寿命越短，利用多项式函数法对放电深度DoD与循环寿命Nctf之间的函数关系进行拟合，放电深度DoD与循环寿命Nctf之间的函数关系为：
其中，ai为i次项对应的系数，N为多项式阶数；
上式中的放电深度DoD是一个与电池充放电状态和有功功率之间存在非线性关系的变量，即 其无法给出显式的数学表达式，因此采用雨流计数法来计算；
步骤2-b，雨流计数法包括如下步骤：
1)将电池储能系统的电池荷电状态–时间曲线顺时针转过90度，此时将曲线看做“屋檐”，雨流依次从起点和每一个峰/谷值开始；
2)雨流在流到峰/谷值处时竖直下滴，一直流到对面有一个比开始时最大值或最小值更大的峰值或更小的谷值为止；
3)根据每一个雨滴流动的起点和终点画出其流动路径，即为每个循环半周期，同时将所有循环半周期的峰谷值记录下来；
4)每一循环中雨流流动路径的水平长度即为电池储能系统在该循环周期的放电深度。
3.根据权利要求1所述的一种含电池储能系统的主动配电网短期有功优化方法，其特征在于：所述步骤4中，目标函数具体组成如下：
步骤a，考虑主配网之间可以进行功率交换且配网侧购售电价不同的情形，则配电侧购电成本和售电收益之和可表示为：
其中， 为第t时段主网向配网输送的有功功率，Ct为配网对主网在第t个调度时段的购售电价，可表示为：
其中， 和C-分别为t时段的配网购电电价和售电电价，单位均为人民币/kW·h，且步骤b，为全面考虑配网侧经济效益，除配网侧购售电成本外，还应考虑电池储能系统的等效运行成本，电池储能系统等效运行成本可表示为：
其中，NB为电池储能系统的数目，Mk为k号电池储能系统的调度周期内总循环周期数，Ce,k,m(DoD,k,m)为第k号电池储能系统的第m次放电深度为DoD,k,m的循环周期下等效投资成本；
CBESS,k为第k号BESS的等效维护成本；
步骤c，考虑到风力和光伏发电无需消耗燃料，假设其发电成本为零，从经济性和环保角度出发，弃风弃光会使得运营商过量排放，违反事先约定的排放权限，因此需要支付违约赔偿费用，则弃风弃光成本可表示为：
其中，cw和cPV分别表示单位弃风弃光成本，单位为：人民币/kW·h；Nw和Npv分别为风电和光伏个数；
步骤d，考虑到微型燃气轮机运行时需要消耗天然气类传统能源，为表示其有功出力与燃料费用的关系，微型燃气轮机的运行成本可表示为：
其中，cFP为燃料价格，单位为：人民币/kW·h； 为第t时段第r台微型燃气轮机的有功出力，NMT为微型燃气轮机数量，ηMT为微型燃气轮机的燃料转化效率；
步骤e，考虑可控负荷为工业负荷，通过与其签订相关合同，保证供电量的情况下，可对其进行适当的调整，考虑到调节可控负荷时势必影响消费者的舒适度，需对其进行一定程度的补偿，补偿成本为：
其中，ccl,l为对第l个可控负荷进行调整时的补偿价格，单位为：人民币/kW·h；Ncl为总可控负荷数目； 为第t时段对第l个可控负荷的负荷有功调整量。
4.根据权利要求1所述的一种含电池储能系统的主动配电网短期有功优化方法，其特征在于：所述步骤5中的原对偶内点法具体实现如下：
采用原对偶内点法将静态变量和动态变量解耦求解，降低求解矩阵的阶数，其中静态变量包括各节点电压幅值和相角、根节点有功和无功功率、风电机组和光伏的有功和无功功率、微型燃气轮机的无功功率、对不等式约束松弛时候的各松弛变量、各等式约束和不等式约束对应的拉格朗日乘子；动态变量包括：电池储能系统和可控负荷的有功和无功功率、电池储能系统和存储型可控负荷的荷电状态、微型燃气轮机的有功功率；上述非线性规划问题包含目标函数、等式约束、不等式约束和时段间耦合关系，形式如下：
其中，f(x)为目标函数，ht(x)＝0为等式约束，glt≤gt(x)≤gut为不等式约束，并且等式约束和不等式约束都包括静态和动态两部分，动态部分包括电池储能系统和存储型可控负荷相关的能量约束、微型燃气轮机的爬坡约束；采用原对偶内点法求解，会得到经典最优潮流修正方程的系数矩阵：
其中，WR1～WRT为对应时段下修正方程降阶后的与静态变量相关的系数矩阵，ΔXR1～ΔXRT为对应时段下所有静态变量和等式约束的拉格朗日乘子变化量，BR1～BRT为对应时段下对拉格朗日函数求静态变量的偏导，ER1～ERT为对应时段下解耦后与静态变量有关的系数矩阵，DR为与动态变量相关的系数矩阵，ΔXRd为动态变量变化量，BRD为对拉格朗日函数求动态变量的偏导；利用其分块对角带边结构，通过线性变换可得：
解耦得到：
将上式求出来的ΔXRd代入下式分别求解各时段静态变量ΔXRt：
5.根据权利要求1所述的一种含电池储能系统的主动配电网短期有功优化方法，其特征在于：所述步骤5-d中所述剪支准则具体有：
(1)该子问题无可行解；(2)所有离散变量均已取得整数解；(3)子问题目标值大于或等于上界；(4)电池储能系统存储的能量不满足上下限约束。