一种基于多维评估指标体系的配电网DG最大准入容量计算方法与流程

一种基于多维评估指标体系的配电网dg最大准入容量计算方法
技术领域
1.本发明涉及电力技术领域，特别涉及一种基于多维评估指标体系的配电网dg最大准入容量计算方法。


背景技术：

2.目前，以分布式光伏为主的分布式电源(distributed generation,dg)并网给配电网的运行带来了一系列的影响，在改善潮流，提升网络电压质量并降低网络损耗的同时，过高的渗透率也将影响电网的安全运行。研究配电网对dg的接纳能力对于合理规划利用dg实现高效节能配电具有十分重要的意义。借助仿真软件如matlab、pscad、opendss进行dg的接纳能力分析，其优点是可以实现配电网各运行参数指标的可视化展示，原理和操作过程相对简单，但仿真过程需重复验算，仿真验证工作量较大，传统的对于dg接入配电网最大渗透率的计算评估，是构建一个最优化的规划计算模型，利用优化算法求解一个以接入容量最大为目标函数，计及多种安全约束，控制变量为dg的最优接入方式的规划问题。对dg接纳能力评估方法主要有解析法、进化算法及随机场景模拟法三种方式，其中解析法包含内点法、线性规划方法、牛顿法等求解非线性优化模型的数学方法；进化算法包含遗传算法、粒子群算法、深度神经网络等人工智能方法；随机场景模拟法主要通过蒙特卡洛抽样依照对象的概率密度函数和分布函数生成抽样样本，并计及各种安全约束计算其接纳能力。此种方法的优势在于其数学建模准确，计算结果可靠，但由于需要计及多个约束条件，建模复杂且计算量极大，对优化算法的计算精度要求高。
3.现有技术中，研究往往着眼于限制dg有效接纳的某一因素进行深入分析建模，并用数学优化方法进行求解，然而影响dg接纳能力的因素众多，若将所有影响因素均在数学模型中进行考虑则会造成建模复杂、计算量大的问题，难以求解。
4.例如，一种在中国专利文献上公开的“一种区域电网新能源消纳容量计算方法和系统”，其公告号：cn108390415a，其申请日：2018年01月30日，该发明根据电网运行方式计算考虑调峰约束后的区域电网新能源消纳容量；将不大于区域电网新能源消纳容量的新能源接入容量分解至区域电网的各分区电网；当各分区电网的新能源接入容量均满足安全稳定要求时，根据各分区电网的新能源接入容量确定区域电网新能源消纳容量，但是存在影响维度片面，模型计算结果准确度低的问题。


技术实现要素：

5.针对现有技术影响维度片面，模型计算结果准确度低的不足，本发明提出了一种基于多维评估指标体系的配电网dg最大准入容量计算方法，综合考虑了影响配电网接纳dg能力的多维影响因素，采用刚性约束阶段和柔性约束阶段进行求解，保证了配电网接纳dg的准确度和潜力。
6.以下是本发明的技术方案，一种基于多维评估指标体系的配电网dg最大准入容量
计算方法，包括以下步骤：
7.s1：选取分析评估维度，所述评估维度包括配电网运行指标维度、电气设备约束维度、dg消纳能力维度和投资及运行经济性维度，基于所述评估维度建立配电网接纳dg能力评估指标体系；
8.s2：分析建立基于加权平均电压偏移率、平均静态电压稳定性、平均供电可用率、期望缺供电能量、加权节点电压谐波总畸变率和电负荷均负荷率的约束目标函数，基于所述约束目标函数建立配电网dg最大准入容量计算模型；
9.s3：基于所述配电网dg最大准入容量计算模型，通过刚性约束阶段和柔性约束阶段计算dg最大准入容量。
10.本方案中，从配电网限制dg接入的四个维度建立了配电网接纳dg能力评估指标体系，结合指标体系建立dg最大准入容量计算模型，综合考虑了影响配电网接纳dg能力的多维影响因素，实现了影响因素的多维融合保证了计算结果的有效性，分析建立基于加权平均电压偏移率、平均静态电压稳定性、平均供电可用率、期望缺供电能量、加权节点电压谐波总畸变率和电负荷均负荷率的约束目标函数，基于所述约束目标函数建立配电网dg最大准入容量计算模型，本模型通过刚性约束、柔性约束的两阶段计算策略改进基于分解的多目标进化算法进行计算模型的求解，旨在减少计算的复杂程度和计算用时，增强模型的实用性。
11.作为优选，步骤s1中，所述配电网运行指标维度包括节点电压偏差维度、静态电压稳定性维度、供电可靠性维度和谐波维度，在所述节点电压偏差维度，选取加权平均电压偏移率指标，其算式如下：
[0012][0013]
式中，u
cs
为加权平均电压偏移率指标，λi为第i个节点的电压偏移权重系数，u
i,t
为第i个节点在t时刻的电压幅值，un为系统标称电压，k为配电网总节点个数，t为整个计算周期；
[0014]
在所述静态电压稳定性维度，选取平均静态电压稳定性指标，所述静态电压稳定性指标为周期内静态电压最大值；
[0015]
在所述供电可靠性维度，选取平均供电可用率指标和期望缺供电能量指标，其算式如下：
[0016][0017][0018]
式中：asai为平均供电可用率指标，eens为缺供电能量指标，c
ij
为第i次停电事件中负荷点j的切负荷总量，t为模拟周期时间，s为模拟周期时间内发生停电的事件集合，n为模拟周期时间内停电次数，d为模拟周期时间内总停电时间，d
ij
为第i次停电事件中负荷点j的停电时间；
[0019]
在所述谐波维度，选取加权平均谐波偏移率指标，其算式如下：
[0020][0021][0022]
式中，thd
μ
为加权节点电压谐波总畸变率指标，thdi为第i个节点在计算周期内的电压总谐波畸变率，u
h,i
为第i节点的h次谐波电压均方根值，m为所考虑的最高谐波次数，μi为第i个节点的电压谐波权重系数，k为配电网总节点个数。
[0023]
本方案中，主要针对以上四点建立配电网运行指标维度的指标集，选取了共计五个指标：在节点电压偏差方面，选取加权平均电压偏移率指标；在静态电压稳定性方面，选取平均静态电压稳定性指标；在供电可靠性方面，选取平均供电可用率指标和期望缺供电能量；在谐波方面，选取加权平均谐波偏移率指标，提高指标维度的广度和数据的有效性。
[0024]
作为优选，步骤s1中，dg消纳能力维度满足发电、用电的动态平衡关系和调节能力的约束，算式如下：
[0025]
pa(t)＝p
l
(t)+p
t
(t)-pr(t)-p
pcc
(t)，
[0026]
p
t,min
(t)≤p
t
(t)≤p
t,max
(t)，
[0027][0028][0029]
式中，v为电负荷均负荷率指标，pa(t)为最大消纳dg电力，p
l
(t)为t时刻的负荷功率，p
t
(t)为t时刻该配网通过联络线向其他电网外送的功率，以下简称联络线功率，送出方向为正，pr(t)为t时刻配网内储能装置的总出力，输出功率为正，蓄能为负，p
pcc
(t)为上级电网的pcc传输功率，e1为计算周期t内平均负荷量，p
t,min
(t)、p
t,max
(t)为联络线在t时刻输送功率的最大值和最小值限制。
[0030]
作为优选，步骤s1中，基于投资成本费用、运维成本费用、运行线损和等价现值计算投资回收期，算式如下：
[0031][0032]
式中，t
inv
为投资回收期，c
inv
为各项治理措施的总投资成本，c
ope
(t)为电网规划治理投资项目的运行维护成本，cr(t)为储能装置运行成本，c
dr
(t)为需求侧响应投入成本，c
t
(t)为通过联络线传输外送功率收益，c
l
(t)为系统线损成本，c
ori
(t)为治理前原电网的运行成本。
[0033]
本方案中，投资及运行经济性指标用于提升配电网dg接纳能力的各项措施的效能分析，为电网规划提出最佳辅助决策。
[0034]
作为优选，步骤s2中，基于加权平均电压偏移率指标、平均静态电压稳定性指标、平均供电可用率指标、期望缺供电能量指标、加权节点电压谐波总畸变率指标和电负荷均
负荷率指标构建约束目标函数，步骤如下：
[0035]
s21：计算控制变量的指标值；
[0036]
s22：计算指标的标准值的偏移量，经过线性变换和归一化处理得到指标的计算值；
[0037]
s23：采用最小二乘法和dare-熵权法确定指标的权重；
[0038]
s24：组合指标的计算值和权重，形成约束目标函数。
[0039]
本方案中，配电网运行指标维度涵盖了dg并网对配电网各节点电压、支路静态稳定性、供电可靠性和谐波畸变四个方面的影响，其偏离额定值的绝对值越小越好。由于负荷及dg出力的不确定性，用电负荷往往达不到最小负荷的极限状态，dg出力也极少达到理论最大出力的情况，部分设备对电压升高等运行状态也有短暂的过载承受能力，配电网运行指标维度在dg最大准入容量计算模型可以量化为约束目标函数，作为柔性约束进行建模，相比将运行指标设置为约束条件的刚性约束计算方法可以更加充分地挖掘配电网dg接纳潜力，在工程实践中具有更好的适应性。
[0040]
作为优选，步骤s2中，所述配电网dg最大准入容量计算模型如下：
[0041][0042]
式中，f1为主目标函数，f2为约束目标函数，p
dg,i
为第i个节点配置的dg最大准入功率，n为所有可以接入dg的节点，t
inv
为投资回收期，ai分别为第i个指标的计算值，λi为ai对应的权重，k1，k2体现了dg最大准入容量同经济性指标的权衡关系，在计算理论最大准入容量时可以将k2取0，具体工程项目中可以根据实际情况对k1，k2进行调整。
[0043]
本方案中，选取dg配置的位置及对应节点的配置容量为控制变量，以dg最大配置容量和投资及运行经济性指标的组合函数作为主目标函数，计算模型在连续时间序列上建立约束目标函数，并将其作为柔性约束，采用基于刚性约束、柔性约束的两阶段计算方法求解，充分挖掘了配电网接纳dg的潜力。
[0044]
作为优选，所述配电网dg最大准入容量计算模型满足配电网运行潮流约束、系统内电气设备约束、各时刻消纳平衡约束、计算周期内消纳平衡约束、计算周期内储能装置充放电容量约束、储能装置充放电功率约束和联络线传输功率约束。
[0045]
本方案中，从配电网运行潮流约束、系统内电气设备约束、各时刻消纳平衡约束、计算周期内消纳平衡约束、计算周期内储能装置充放电容量约束、储能装置充放电功率约束和联络线传输功率约束的角度对模型进行约束，提高模型的实用性。
[0046]
作为优选，所述配电网dg最大准入容量计算模型采用moea/d算法进行求解，将多目标最优前沿的逼近问题转化为多个子目标模型的求解问题，对每一个子问题，在迭代过程中利用与其相邻的子问题的信息实现求解。
[0047]
作为优选，步骤s3中，通过刚性约束阶段和柔性约束阶段计算dg最大准入容量的步骤如下：
[0048]
s31：准备数据，对所述数据进行预处理；
[0049]
s32：建立基于刚性约束的单目标模型，求解得到计算初值；
[0050]
s33：建立基于刚性约束和柔性约束的双目标模型，将所述双目标模型分解为两个子目标模型进行求解，得到所述计算初值临近区域内的多目标帕累托前沿解集。
[0051]
本方案中，通过刚性约束、柔性约束的两阶段计算策略改进基于分解的多目标进化算法进行计算模型的求解，减少计算的复杂程度和计算用时，增强模型的实用性。
[0052]
本发明的有益效果是：综合考虑了影响配电网接纳dg能力的多维影响因素，采用刚性约束阶段和柔性约束阶段进行求解，保证了配电网接纳dg的准确度和潜力。
附图说明
[0053]
图1本发明一种基于多维评估指标体系的配电网dg最大准入容量计算方法的流程图。
[0054]
图2本发明一种基于多维评估指标体系的配电网dg最大准入容量计算方法的dg消纳空间限制示意图。
[0055]
图3本发明一种基于多维评估指标体系的配电网dg最大准入容量计算方法的计算模型指标的处理方法。
[0056]
图4本发明一种基于多维评估指标体系的配电网dg最大准入容量计算方法的moea/d算法计算流程图。
[0057]
图5本发明一种基于多维评估指标体系的配电网dg最大准入容量计算方法的两阶段计算策略流程示意图.
[0058]
图6本发明一种基于多维评估指标体系的配电网dg最大准入容量计算方法的改进ieee33节点拓扑图。
[0059]
图7本发明一种基于多维评估指标体系的配电网dg最大准入容量计算方法的典型日总负荷曲线图。
[0060]
图8本发明一种基于多维评估指标体系的配电网dg最大准入容量计算方法的理论pv出力曲线图。
[0061]
图9本发明一种基于多维评估指标体系的配电网dg最大准入容量计算方法的场景模拟评估结果比例示意图。
[0062]
图10本发明一种基于多维评估指标体系的配电网dg最大准入容量计算方法的节点电压越限时间比例示意图。
具体实施方式
[0063]
下面通过实施例，并结合附图，对本发明的技术方案作进一步具体的说明。
[0064]
实施例：如图1所示，一种基于多维评估指标体系的配电网dg最大准入容量计算方法，包括以下步骤：
[0065]
步骤1：选取分析评估维度，所述评估维度包括配电网运行指标维度、电气设备约束维度、dg消纳能力维度和投资及运行经济性维度，基于所述评估维度建立配电网接纳dg能力评估指标体系；
[0066]
步骤2：分析建立基于加权平均电压偏移率、平均静态电压稳定性、平均供电可用率、期望缺供电能量、加权节点电压谐波总畸变率和电负荷均负荷率的约束目标函数，基于所述约束目标函数建立配电网dg最大准入容量计算模型；
[0067]
步骤3：基于所述配电网dg最大准入容量计算模型，通过刚性约束阶段和柔性约束阶段计算dg最大准入容量。
[0068]
大量分布式光伏电源(distributed photovoltaic generation，本发明简称dg)接入配电网后，配电网拓扑将由简单的单向潮流模式转换为多点电源供电、双向潮流运行的复杂结构。配电网的各节点电压、电压波动、供电可靠性、谐波含量等受到很大的影响，在直接影响系统各运行指标的同时，间接影响了电网一次设备、二次设备的选型约束。配电网接纳dg能力进行评估的基础，是要对dg并网后对配电网各运行指标的影响及其出力特征进行充分分析和把握。综合以上评估维度结合投资及运行经济性指标维度选取特征指标，建立了配电网接纳dg能力评估指标体系，作为构建dg最大准入容量计算模型的基础。分别从配电网运行指标维度、电气设备约束维度、dg消纳能力维度、投资及运行经济性维度四个维度对配电网接纳dg能力评估指标体系进行指标选取分析。
[0069]
配电网运行指标维度。对于配电网稳态运行，dg的并网接入主要对节点电压偏差、静态电压稳定性、供电可靠性、谐波产生影响，针对以上四点影响建立配电网运行指标维度的指标集，选取了共计五个指标：在节点电压偏差方面，选取加权平均电压偏移率指标；在静态电压稳定性方面，选取平均静态电压稳定性指标；在供电可靠性方面，选取平均供电可用率指标asai和期望缺供电能量指标eens；在谐波方面，选取加权平均谐波偏移率指标。运行指标维度中，根据配电网内各负荷对电压偏移、谐波畸变的敏感程度不同，应用权重向量在模型中对各节点进行差异化分析及处理，得到更好的指导和评估效果。
[0070]
加权平均电压偏移率指标。由于dg并网，配电网的供电形式由传统配电网单向辐射转变为多点网状的结构，潮流不再是由单电源端以辐射状的单一方向流向用电负荷，各节点电压幅值和相角将随着dg接入而产生变化。系统轻载运行时，潮流方向发生改变，可能导致电压降落方向发生变化，使得负载侧电压逐步升高。高渗透率的dg并网影响了系统电压稳定性，甚至可能导致脱网事故。10kv及以下三相供电电压允许偏差为标称电压的+7％至-7％，偏离的绝对值越小越好。系统电压偏差可以表征各节点电压偏差绝对值的均值，其计算方法如下：
[0071][0072]
式中，us为电压偏差，ui为配电网内第i个节点电压幅值，un为系统标称电压，k为配电网总节点个数。
[0073]
在配电网中，各个负荷的重要程度及其对电压偏移的敏感程度不同。考虑重要负荷的电压偏移率往往能得到更好的指导和评价效果，引入加权平均电压偏移率指标。用赋权方法依照各节点负荷的重要程度和其对电压的敏感程度形成节点电压权重系数向量[λ1,λ2,...,λk]，基于负荷类别的节点电压权重系数向量确定方法，类比供电可靠性的三类负荷划分，将居民用电负荷等对供电电压不敏感类型负荷定义为c类，节点电压权重取0.2；商业用电、办公楼等公用设施用电负荷定义为b类，节点电压权重取0.5；精密工厂等对供电电压要求严格的用电负荷定义为a类，节点电压权重取1.0。用该系数向量对系统电压偏移指标进行修正，可以计算某时刻的加权平均电压偏移率。同时考虑到计算周期内整体的平均电压偏移指标，建立计算周期内的系统加权平均电压偏移率指标u
cs
，其计算方法如下：
[0074][0075]
式中，u
cs
为加权平均电压偏移率指标，λi为第i个节点的电压偏移权重系数，u
i,t
为第i个节点在t时刻的电压幅值，un为系统标称电压，k为配电网总节点个数，t为整个计算周期。
[0076]
平均静态电压稳定性指标。静态电压稳定性指标反映了系统负荷电压稳定的关系，用等效阻抗变换方法将多节点配电网转换为两节点配电网实现指标计算，但配电网拓扑复杂此方法误差较大且忽略了各节点的节点电压和负荷分布对静态电压稳定的影响。对静态电压稳定性的修正计算方法，定义t时刻支路b
ij
的静态电压稳定性指标为l
ij,t
，其计算方法如下：
[0077][0078]
l
t
＝max{l
b,t
}
ꢀꢀꢀ
(4)
[0079]
式中，l
ij,t
为t时刻支路b
ij
的静态电压稳定性，p
j,t
为节点j在t时刻的注入有功功率，x
ij
为支路b
ij
的线路电抗值，q
j,t
为为节点j在t时刻的注入无功功率，r
ij
为支路b
ij
的线路电阻值，u
i,t
为第i个节点在t时刻的电压幅值，l
t
为平均静态电压稳定性指标，l
b,t
为t时刻该配电网中全部支路的静态电压稳定指标集合。
[0080]
静态电压稳定性指标l越小表示网络的电压稳定性越好；l越大表示网络的稳定性越差，当l接近于1时系统具有较大的电压崩溃风险。在计算周期内，网络的静态电压稳定性指标l同样选取周期内l
t
最大值来表示，其计算方法如下：
[0081]
l＝max{l
t
},t＝1,2,...,t
ꢀꢀꢀ
(5)
[0082]
式中，l为静态电压稳定性指标，l
t
为选取周期t内电压指标，t为周期。
[0083]
供电可靠性指标。可再生能源发电占比较低时，传统的发电机组能够完全满足电力电量的平衡要求，可再生能源出力仅作为电力系统的备用补充。在电网传统机组发生故障导致供电能力不足形成计划孤岛时，孤岛内的dg可以起到在孤岛内传输电力的作用，将提高配电网的供电可靠性。随着分布式可再生能源并网比例的进一步提高，源荷将呈现较大的波动性和不确定性特点，电网的运行状态将更为复杂。传统机组不再能够独立满足负荷需求，以水电、光伏、风电为主的可再生能源将承担一部分负荷平衡的责任。此时，dg由于其自身的不确定性造成出力下降无法满足负荷需求将会带来部分负荷断电或系统供电电压质量下降现象，降低了系统可靠性。随着dg接入比例的逐渐增加，配电网供电可靠性将呈现先上升后下降的特点。选取平均供电可用率指标asai和期望缺供电能量指标eens来表征配电系统的整体供电可靠性水平，其计算方法如下：
[0084][0085]
式中，asai为平均供电可用率指标，t为模拟周期时间，s为模拟周期时间内发生停电的事件集合，其中n为模拟周期时间内停电次数，d为模拟周期时间内总停电时间，d
ij
为第i次停电事件中负荷点j的停电时间。asai主要体现的是平均停电时间的指标。
[0086][0087]
式中，eens为缺供电能量指标，c
ij
为第i次停电事件中负荷点j的切负荷总量，t为模拟周期时间，s为模拟周期时间内发生停电的事件集合，其中n为模拟周期时间内停电次数，d为模拟周期时间内总停电时间，d
ij
为第i次停电事件中负荷点j的停电时间。eens主要体现缺供电能量的指标。
[0088]
电压谐波总畸变率指标。dg接入配电网，并网接口的电力电子装置将向配电网侧注入大量谐波，谐波引起的过电压、谐波畸变率过大问题将影响电气设备的正常运行，各个谐波源的谐波不同且其具有很强的随机性，应用概率谐波潮流分析方法充分考虑dg谐波的波动性和不确定性。节点畸变波形偏离基波波形的程度，一般以总谐波畸变率thd来表示，其计算方法如下：
[0089][0090]
式中，thdi为第i个节点在计算周期内的电压总谐波畸变率，u
h,i
为第i节点的h次谐波电压均方根值，m为所考虑的最高谐波次数。
[0091]
与电压指标相同，各个节点负荷对谐波畸变率的敏感程度不同，可以引入节点电压谐波权重系数向量[μ1,μ2,...,μk]对系统所有节点的平均电压总谐波畸变率进行修正，节点电压谐波权重系数的取值同节点电压权重系数，形成加权节点电压谐波总畸变率指标thd
μ
，其计算方法如下：
[0092][0093]
式中，thd
μ
为加权节点电压谐波总畸变率指标，thdi为第i个节点在计算周期内的电压总谐波畸变率，u
h,i
为第i节点的h次谐波电压均方根值，m为所考虑的最高谐波次数，μi为第i个节点的电压谐波权重系数，k为配电网总节点个数。
[0094]
电气设备约束维度。dg接入后将引起电流分布的变化，由于电网一次设备的载流限制，必须要求配电网支路电流不超过该支路上工作的电气一次设备的工作电流耐受能力，其计算方法如下：
[0095]
ik≤i
xk
ꢀꢀꢀ
(10)
[0096]
式中，k为配电网中一次设备编号；ik为对应的支路电流；i
xk
为该一次设备允许的最大工作电流。
[0097]
与此同时要求最大短路电流不超过该支路上工作的电气一次设备的短路电流耐受能力，其计算方法如下：
[0098]idk
≤i
dxk
ꢀꢀꢀ
(11)
[0099]
式中，i
dxk
为该一次设备允许最大短路电路；i
dk
为最大短路电流。
[0100]
显然提高一次设备的载流限制可以有效提高配电网接纳dg的能力，提高dg最大准入容量。
[0101]
dg消纳能力维度。对于配电网系统，实现dg消纳需要满足发电、用电的动态平衡关
系并满足调节能力的约束。dg消纳空间限制示意图如图1所示。在系统t时刻，忽略网损的情况下最大消纳dg电力pa(t)，其计算方法如下：
[0102]
pa(t)＝p
l
(t)+p
t
(t)-pr(t)-p
pcc
(t)
ꢀꢀꢀ
(12)
[0103]
式中，pa(t)为最大消纳dg电力，p
l
(t)为t时刻的负荷功率，p
t
(t)为t时刻该配网通过联络线向其他电网外送的功率，以下简称联络线功率，送出方向为正，pr(t)为t时刻配网内储能装置的总出力，输出功率为正，蓄能为负，p
pcc
(t)为上级电网的pcc传输功率。
[0104]
其中联络线功率p
t
(t)必须满足通道能力的约束，其计算方法如下：
[0105]
p
t,min
(t)≤p
t
(t)≤p
t,max
(t)
ꢀꢀꢀ
(13)
[0106]
式中，p
t,min
(t)、p
t,max
(t)为联络线在t时刻输送功率的最大值和最小值限制。
[0107]
由于储能装置的充放电特性限制以及配电区域间供需电能关系的限制，系统dg消纳电量的分析应建立在长时间序列上，以一个长时间序列周期t为例(通常可以取一天24h、或一周7天)，系统dg消纳电量ea为最大消纳dg电力pa(t)的积分，其计算方法如下：
[0108][0109]
式中，ea为系统dg消纳电量，pa(t)为最大消纳dg电力，p
l
(t)为t时刻的负荷功率，p
t
(t)为t时刻该配网通过联络线向其他电网外送的功率，以下简称联络线功率，送出方向为正，pr(t)为t时刻配网内储能装置的总出力，输出功率为正，蓄能为负，p
pcc
(t)为上级电网的pcc传输功率。
[0110]
量化分析由于负荷峰谷差对系统消纳dg能力的制约和影响，引入电负荷均负荷率表示一个计算周期t内平均负荷量与最大负荷量的比值关系，电负荷均负荷率用v表示，其计算方法如下：
[0111][0112]
式中，v为电负荷均负荷率指标，p
l
(t)为t时刻的负荷功率，e1为计算周期t内平均负荷量。
[0113]
投资及运行经济性维度。如图2所示，增加配电网消纳dg的能力，通过风、光等可再生能源的出力特性实现多能互补，增加生物质能发电、燃气轮机发电等灵活调节的分布式电源进行调节；从网侧灵活性可以考虑基于dg源出力的时间特性安排外送，利用联络线转供至其他电压等级或相同电压等级的其他电网，增加本区域dg消纳空间，增强电网的互联互通能力。在可再生资源丰富的地区，可以将过剩的能源通过外送到其他区域增加消纳，并实现调节的灵活性；同时，从储侧考虑储能装置的配置和灵活控制方式可以实现负荷削峰填谷，系统增大消纳能力，但与此同时储能装置具有充放电特性、容量限制等多方面限制，需要从长时间序列进行分析考虑，同时其投资成本和运行成本较高，在实际运用时需要重点探讨其规划配置方式和控制策略，以期达到最好的效果；从荷侧考虑需求侧响应技术，例如新型负荷具有时空可转移特性，在实现负荷曲线形状改变的同时也可调节系统潮流，建立友好的辅助市场交易体制可以更好地引导用户用电行为，增强负荷用能曲线和dg出力供能曲线的匹配程度，提高配电网消纳dg的能力。
[0114]
投资及运行经济性指标则可以用于以上提升配电网dg接纳能力的各项措施的效能分析，为电网规划提出最佳辅助决策，选取投资回收期指标来进行电网规划治理投资项目经济性的评估。
[0115]
投资回收期t
inv
是涉及了投资成本费用、运维成本费用、运行线损和等价现值等项目全过程的综合经济性指标，其表征具体电网规划项目实施投产后获得的收益总额达到该项目投资成本所需要的时间，其计算方法如下：
[0116][0117]
式中，t
inv
为投资回收期，c
inv
为各项治理措施的总投资成本，c
ope
(t)为电网规划治理投资项目的运行维护成本，cr(t)为储能装置运行成本，c
dr
(t)为需求侧响应投入成本，c
t
(t)为通过联络线传输外送功率收益，c
l
(t)为系统线损成本，c
ori
(t)为治理前原电网的运行成本。
[0118]
基于评估指标体系的配电网dg最大准入容量计算模型。为了求解配电网的dg最大准入容量，可以通过建立规划模型求解最优解的形式来实现，为此应选取dg接入的最大容量及其对应的配置位置为控制变量，以配电网各节点负荷、电气设备及拓扑信息等其他运行参数为扰动变量。为了在目标函数和约束条件的构造中考虑配电网接纳dg能力的约束限制，基于配电网接纳dg能力评估指标体系建立了dg最大准入容量的计算模型。在建立模型前首先对指标体系的各项指标进行详细分析。
[0119]
如图3所示，评估体系中各指标的处理方法。配电网运行指标维度涵盖了dg并网对配电网各节点电压、支路静态稳定性、供电可靠性和谐波畸变四个方面的影响。在以往研究中大多将运行指标作为刚性约束条件进行限制，并以最小负荷场景的dg最大理论出力进行校验，保证各时刻各项运行指标均不超过限制。由于负荷及dg出力的不确定性，用电负荷往往达不到最小负荷的极限状态，dg出力也极少达到理论最大出力的情况。同时，部分设备对电压升高等运行状态也有短暂的过载承受能力。配电网运行指标维度在dg最大准入容量计算模型可以量化为约束目标函数，作为柔性约束进行建模。此方法相比将运行指标设置为约束条件的刚性约束计算方法可以更加充分地挖掘配电网dg接纳潜力，在工程实践中具有更好的适应性。dg消纳能力维度中的电负荷均负荷率指标与配电网运行指标极为相似，同样也在约束目标函数中进行考虑。
[0120]
建立了包含加权平均电压偏移率(公式(2))、平均静态电压稳定性(公式(3)、(4)、(5))、平均供电可用率asai(公式(6))、期望缺供电能量eens(公式(7))、加权节点电压谐波总畸变率(公式(9))、电负荷均负荷率(公式(15))共六项指标组成的约束目标函数。
[0121]
构建基于以上六个指标的约束目标函数，采用归一化方法和主客观赋权方法方法，其步骤如下：
[0122]
步骤21：计算当前控制变量取值情形下各指标值；
[0123]
步骤22：计算各指标与其对应标准值的偏移量并线性变换到[0,1]之间，得到经归一化后第i个指标计算值ai。以10kv配电网为研究对象，节点电压标准值取10kv，静态电压稳定性标准值取0，asai标准值取1，eens标准值取0，电压谐波畸变率标准值取0，电负荷均负荷率标准值取1。
[0124]
步骤23：采用最小二乘法综合dare-熵权法的主客观赋权方法确定指标权重λi。
[0125]
步骤24：组合权重λi及指标计算值ai形成约束目标函数f2，其计算方法如下：
[0126][0127]
式中，f2为约束目标函数，ai分别为第i个指标的计算值，λi为ai对应的权重。
[0128]
电气设备约束维度指标集表征了系统电气设备热稳定性、动稳定性对各运行指标的限制。由于过载会对降低电气设备的电气寿命，影响系统安全运行，故电气设备约束维度指标无法参照配电网运行指标维度构建约束目标函数作为柔性约束，而是需要作为刚性约束直接加入计算模型的约束条件。
[0129]
dg消纳能力维度指标集中，电负荷均负荷率已在约束目标函数中进行考虑，运行消纳平衡条件需要在计算中的每一时刻进行校验，由于储能装置的充放电容量限制同样也要在整个计算周期内进行校验，需要作为刚性约束加入计算模型的约束条件进行限制。
[0130]
投资及运行经济性指标用于配电网分布式电源接纳能力提升措施的效能分析，应加入主目标函数与dg接入容量进行比较。
[0131]
作为多目标规划模型，选取dg配置的位置及对应节点的配置容量为控制变量，以dg最大配置容量和投资及运行经济性指标的组合函数作为主目标函数，计算模型在连续时间序列上建立约束目标函数，并将其作为柔性约束，采用基于刚性约束、柔性约束的两阶段计算方法求解，充分挖掘了配电网接纳dg的潜力。构造的计算模型目标函数如下所示：
[0132][0133]
式中，f1为主目标函数，f2为约束目标函数，p
dg,i
为第i个节点配置的dg最大准入功率，n为所有可以接入dg的节点，t
inv
为投资回收期，ai分别为第i个指标的计算值，λi为ai对应的权重，k1，k2体现了dg最大准入容量同经济性指标的权衡关系，在计算理论最大准入容量时可以将k2取0，具体工程项目中可以根据实际情况对k1，k2进行调整。
[0134]
配电网dg最大准入容量计算模型满足以下约束条件。计算模型建立在配电网正常运行状态下，故其需要满足配电网运行潮流的等式约束：
[0135][0136][0137]
式中，p
dg,i
为节点i处dg向节点注入的有功功率，p
l,i
为节点i处负荷向节点注入的有功功率，ui为节点i处电压幅值，g
ij
为节点i与节点j间线路电导值，uj为节点j处电压幅值，θ
ij
为节点i与节点j间电压相位差，b
ij
为节点i与节点j间线路电纳值，q
dg,i
为节点i处dg向节点注入的有功功率，q
l,i
为节点i处负荷向节点注入的无功功率，n为配电网总节点个数。
[0138]
同时dg并网后需要满足系统内电气设备约束，校验公式如公式(12)、公式(13)所示。
[0139]
各时刻消纳平衡约束：
[0140]
pa(t)＝p
l
(t)+p
t
(t)-pr(t)-p
pcc
(t)
ꢀꢀꢀ
(21)
[0141]
式中，pa(t)为最大消纳dg电力，p
l
(t)为t时刻的负荷功率，p
t
(t)为t时刻该配网通过联络线向其他电网外送的功率，以下简称联络线功率，送出方向为正，pr(t)为t时刻配网内储能装置的总出力，输出功率为正，蓄能为负，p
pcc
(t)为上级电网的pcc传输功率。
[0142]
计算周期内消纳平衡约束：
[0143][0144]
式中，pa(t)为最大消纳dg电力，p
l
(t)为t时刻的负荷功率，p
t
(t)为t时刻该配网通过联络线向其他电网外送的功率，以下简称联络线功率，送出方向为正，pr(t)为t时刻配网内储能装置的总出力，输出功率为正，蓄能为负，p
pcc
(t)为上级电网的pcc传输功率，t为计算周期。
[0145]
计算周期内储能装置充放电容量约束：
[0146][0147]
式中，pr(t)为t时刻配网内储能装置的总出力，t为计算周期。
[0148]
储能装置充放电功率约束：
[0149]
p
r,i,min
≤p
r,i
(t)≤p
r,i,max
ꢀꢀꢀ
(24)
[0150]
式中，p
r,i
(t)为i节点所接储能装置在t时刻的充放电功率，p
r,i,min
、p
r,i,max
分别为其充放电功率限制下限及上限。
[0151]
联络线传输功率约束：
[0152]
p
t,i,min
≤p
t,i
(t)≤p
t,i,max
ꢀꢀꢀ
(25)
[0153]
式中，p
t,i
(t)为第i条联络线在t时刻的传输功率，p
t,i,min
、p
t,i,max
分别为输送功率的最小值和最大值限制。
[0154]
dg最大准入容量计算模型包含主目标函数和约束目标函数两个目标函数，属含约束的多目标规划问题。选用改进的基于分解的多目标进化算法(moea/d算法)进行计算模型的求解。在电力行业中，moea/d算法在运行优化协调控制技术原理是将多目标模型通过组合权重转化为多个单目标子模型，具有高搜索效率、高优化效率、能够提高最优解性能及计算效率等多方面优点。moea/d算法由分解算法和进化算法两部分组成。分解算法将多目标最优前沿的逼近问题转化为多个子目标模型的求解问题；进化算法是对每一个子问题，在迭代过程中利用与其相邻的子问题的信息实现求解。moea/d算法在迭代过程中的每一个解均是由各个子问题产生的最优解组成的前沿解集，利用相邻的子问题共享优化信息特性，可以避免进行重复的规划计算，从而在保障算法寻优能力的同时避免陷入局部最优的问题。如图4所示，moea/d算法的计算流程图。
[0155]
moea/d的分解算法有加权和法、切比雪夫法、边界交集法。其中切比雪夫法对最优前沿形状不敏感，对凸规划和非凸规划均有很好的适应性。moea/d的进化算法运用最多的是子问题差分进化算法。与遗传算法等传统进化算法相比，差分进化算法参数少，鲁棒性和全局优化能力强。
[0156]
moea/d算法是无约束问题的求解算法，但本模型需要考虑多种约束条件，将约束
条件利用惩罚函数法构造进目标函数进行处理，对不可行解进行惩罚保证优化过程在可行域内进行。如图5所示，计算流程示意图，设计了基于刚性约束、柔性约束的两阶段计算策略进行模型的求解。
[0157]
第一阶段为基于刚性约束模型的求解，将约束目标函数中各指标建立为刚性约束条件，基于主目标函数和约束条件构建单目标规划模型进行求解，得到dg接入容量的保守解，此保守解将保证所有控制变量遵从约束条件的可行域。
[0158]
第二阶段建立约束目标函数作为柔性约束，以moea/d算法进行多目标模型的求解，应用切比雪夫分解法将主目标函数、约束目标函数分解为两个子目标模型，并通过差分进化算法对两个子目标模型进行求解。以第一阶段刚性约束模型的计算结果作为计算初值，在其临近区域内寻找多目标帕累托前沿解集。此方法确定计算初值可以迅速得到最优解节省运算时间同时提升解集质量。
[0159]
具体步骤如下：
[0160]
步骤31：准备数据，生成计算场景集，对数据进行预处理；
[0161]
步骤32：将运行指标构建为刚性约束，建立为传统单目标模型，用差分进化算法对单目标模型进行求解，得到计算结果；
[0162]
步骤33：构建约束目标函数作为柔性约束，建立为双目标模型，用moea/d算法进行多目标模型的求解，用切比雪夫分解法分解为两个子目标模型，用差分进化算法求解两个子目标模型，以步骤32的计算结果作为计算初值，在其临近区域内寻找多目标帕累托前沿解集。
[0163]
模型具有极佳的可拓展性，在具体解决实际问题时可基于当地供电负荷特征、用电需求对模型进行修改，使其更佳贴合实际情形，更佳具有工程应用价值。本计算模型也为提升配电网接纳dg能力的方法研究提供了模型仿真方法，为电网规划提出最佳辅助决策。
[0164]
实施方案：
[0165]
本发明示范例选用改进的ieee33节点算例模拟南方某条10kv供电出线，运用本发明构建模型计算其dg最大准入容量研究其接纳dg能力。如图6所示，改进ieee33节点拓扑图，由于模型建立在连续时间序列上，为了让示范例仿真更加具有实际意义，如图7所示，典型日总负荷曲线，选取电网实测的负荷曲线对示范例负荷进行修正，建立一天96时段的连续负荷数据场景。光伏理论出力曲线采用当地一实际光伏电站出力曲线如图8所示。
[0166]
计算结果分析。首先在第一阶段计算过程中，运用刚性约束的单目标规划方法，以dg接入总量为目标函数，约束目标函数转化为硬约束条件进行计算，运用差分进化算法进行求解得到模型计算初值，计算得到总接入容量为880kw。本模型在基于电压约束的同时也包含了可靠性指标和谐波畸变率指标，运行环境更为复杂。综合考虑了影响配电网接纳dg能力的多维影响因素，实现了影响因素的多维融合保证了计算结果的有效性。
[0167]
根据原理分析根据刚性约束模型计算的初值应略小于最终计算结果，在此初值的基础上，第二阶段运用moea/d算法结合本发明所提出计算模型构建约束目标函数计算dg最大准入容量，计算得到总接入容量为1380kw。两个阶段计算的接入节点位置及各节点接入容量详见表1所示。相较刚性约束的单目标模型计算得到的总接入容量提升了56.8％。
[0168]
表1两阶段计算dg最大准入容量计算结果及各节点接入容量对比
[0169][0170]
为了验证计算模型的效果，在上述dg最大接入配置的情形下，模拟全年系统运行仿真获得8760组dg运行场景，在每一个场景下通过潮流计算方法得到系统内各节点的电压幅值，其计算结果如图8所示。
[0171]
根据落入电压不合格区域内点的比例判定方法的可行性。如图9所示，从场景模拟评估结果比例示意图可以看出，在全部潮流共289080组节点电压运行数据中，有12524组节点电压幅值落在不合格区域内，仅占总体的4.33％，同时电压越限值小于9.8％的有844组，仅占总体的0.29％。电压越限集中出现在节点12～节点17主干线末端负荷中，在这些节点处若安装无功调节装置将可以有效解决电压越限问题。
[0172]
如图10所示，仿真全年数据的节点电压越限时间的比例示意图。全年各节点仅出现434次电压越限状况，最大越限幅值为10.92kv，76.5％的越限时间在2个小时之内。结果表明，本模型可以充分挖掘配电网接纳dg的潜力，在建立柔性约束的情形下显著提升了dg总接入容量，同时在大样本的数据支持下系统内仅有极少且短暂的电压越限现象出现，充分证明了其有效性和优越性。