本发明一种基于供需互动和dg运行特性的主动配电网网架规划方法，属于基于供需互动和dg运行特性的主动配电网网架规划
技术领域：
。
背景技术：
：配电网是电力系统的重要组成部分，是输电网和用户之间的重要中间环节；在传统配电系统中，电力潮流一般由上端变电站单一流向负荷节点，其运行方式和规划准则相对简单，配电网与用户侧保持着简单而稳定的供需关系，但随着分布式能源、电动汽车、储能不断增大的普及率以及可控负荷接入配电网，用户需求侧和供应侧稳定的供需关系发生改变，二者之间的频繁互动对系统潮流分布、电压水平、短路容量等原有的电气特性造成显著的影响；而传统配电网在设计阶段并未考虑上述因素，因此难以满足低碳经济背景下高渗透率可再生能源发电接入与高效利用的要求。在此背景下，传统依据潮流单向流动运行模式的配电网网架规划方法不再适用，供需互动给配电网网架规划带来新的问题和挑战，主动配电网应运而生，相比传统配电网它能够主动对配电网运行进行控制管理；因此，需要一种新的一种基于供需互动和dg运行特性的主动配电网网架规划方法，以提高配电网各类能源资源的使用效率，降低投资及运营的成本，同时引导用户有序用电，全面提升配电网经济效益、运行效益以及风险控制能力。技术实现要素：本发明为克服现有技术中存在的不足，所要解决的技术问题为：提供一种基于供需互动和dg运行特性的主动配电网网架规划方法；为解决该技术问题，本发明采用的技术方案为：一种基于供需互动和dg运行特性的主动配电网网架规划方法，包括以下步骤：步骤一：运用基于wasserstein距离的最优场景生成方法处理分布式电源dg如风电、光伏出力的不确定性，生成风光分布式发电联合运行的若干确定性场景；步骤二：构建以年网络综合费用最优为主动配电网网架规划的总目标函数；步骤三：将步骤二的总目标函数分解为上层的网架规划问题和下层的配电网运行优化问题分别进行规划，所述网架规划问题以年网络损耗费用、网架投资费用等年值和年dg固定投资费用之和最小为目标函数，所述配电网运行优化问题以年dg运行维护费用、年主动管理费用和年需求侧管理费用之和最小为目标函数；步骤四：利用交替迭代方法对步骤三的上层的网架规划问题和下层的配电网运行优化问题进行求解：即先优化出上层的最优的一组网架结果，然后将最优的网架结果逐个代入到下层对dg出力和可中断负荷量进行优化，然后将下层优化结果传递给上层计算总目标函数，得到最终优化网架结果。所述步骤二中基于wasserstein距离的最优场景生成方法，由zs(s＝1,2,…,s)表示最优分点位置，表达式如下：所述zs点对应的概率ps数学表达式如下：式中，s表示将f(x)离散化后得到的离散场景内离散点个数，f(x)是变量x的连续概率密度函数；所述步骤一中风电、光伏出力的不确定性分别采用两参数beta分布和weibull分布进行描述：式中，i为光照强度；imax表示i的最大值；α和β是beta分布的两个参数，都取为0.95；γ(·)是伽马函数；c和k分别为尺度和形状两个参数，分别取8.92、2.30。所述步骤一中运用基于wasserstein距离的最优场景生成方法处理分布式电源dg如风电、光伏出力的不确定性包括以下步骤：(1)、根据风速、光照强度的概率密度分布函数以及风速-风电功率、光强-光伏功率函数关系，得到风力发电功率、光伏发电功率的概率密度函数；(2)、采用权利要求2中基于wasserstein距离的最优场景生成方法对风电出力和光伏出力分别进行场景生成，分别将两者生成场景数、对应场景概率值进行交叉相乘得到风光联合运行场景及对应场景的概率。所述步骤三中所述网架规划问题的目标函数表达式如下：minc1＝closs+cline+cf；式中，closs为年损耗费用；cline为网架投资费用等年值；cf折算到每年的dg固定投资费；所述步骤三中所述配电网运行优化问题的目标函数如下所示：minc2＝com+cam+cdsm；式中，com为dg年运行维护费；cam为dg年主动管理费；cdsm为需求侧管理成本。本发明相对于现有技术具备的有益效果为：本发明的提供的一种基于供需互动和dg运行特性的主动配电网网架规划方法不仅简化了规划，并充分考虑了风、光分布式电源出力的不确定性，同时考虑了能够优化运行的需求侧管理措施，是更符合实际的一种配电网网架规划方法；针对网架和运行同时优化可能有维数灾以及收敛性较差问题，本发明提出分解协调联合规划的思想，即对网架规划和运行优化不使用联合编码的形式，而是将其进行转化和分解，变为两个问题进行求解，但是二者之间通过相互迭代求得最优解。附图说明下面结合附图对本发明做进一步说明：图1为双层规划模型上下层规划传递关系图；图2规划模型求解流程图；图3改进29节点配电网区域图；图4方案一的最优网架方案图；图5方案二的最优网架方案图；图6方案三的最优网架方案图；图7三种规划方案各项费用比较图。具体实施方式下面结合附图和具体实施例，进一步阐明本发明，应理解这些实施例仅是本发明的优选实施方式，在阅读了本发明之后，本领域技术人员对本发明的各种不脱离本发明原理的修改均落于本申请所附权利要求所限定的范围。本发明的原理为：针对目前的主动配电网网架规划研究较少，相关规划方法存在求解方法较复杂、未考虑需求侧管理措施、目标函数单一等不足，本发明根据分解协调的思想将主动配电网网架规划模型简化成双层模型，并采用基于wasserstein距离的最优场景生成方法处理了风、光分布式电源出力的不确定性，同时考虑了运行过程中适当中断可中断负荷的需求侧管理措施；双层规划模型分为上层的网架规划问题和下层的配电网运行优化问题分别进行规划，所述网架规划问题以年网络损耗费用、网架投资费用等年值和年dg固定投资费用之和最小为目标函数进行网架的优化配置，所述配电网运行优化问题以年dg运行维护费用、年主动管理费用和年需求侧管理费用之和最小为目标函数进行运行优化；最后采用改进的29节点配电系统作为具体实施例验证了所提模型的合理性和有效性。本发明的一种基于供需互动和dg运行特性的主动配电网网架规划方法，包括以下步骤：步骤一：运用基于wasserstein距离的最优场景生成方法处理分布式电源(dg)如风电、光伏出力的不确定性，生成风光分布式发电联合运行的若干确定性场景，具体为：对于风电、光伏分布式电源的不确定性处理，首先风速大小、光照强度的不确定性分别采用两参数beta分布和weibull分布进行描述：式中，i为光照强度；imax表示i的最大值；α和β是beta分布的两个参数，都取为0.95；γ(·)是伽马函数；c和k分别为尺度和形状两个参数，分别取8.92、2.30。光伏发电功率ppv与光照强度i可以处理为正比关系，由此，可以推导出光伏发电功率的概率密度函数：式中，ppvn是imax对应的光伏发电功率。风力发电功率和风速关系函数可以表示如下：式中，vci、vco、vr是切入风速、切出风速和额定风速；pwtn是风力发电的额定功率。结合上式和风速大小概率密度函数推导出风力发电功率pwt的概率密度函数：式中，δ(·)为冲击函数。基于wasserstein距离的最优场景生成方法是用离散的场景及对应离散场景的概率值来代替连续场景下概率密度函数：式中，zs(s＝1,2,…,s)表示最优分点位置；s表示将f(x)离散化后得到的离散场景内离散点个数；f(x)是变量x的连续概率密度函数。zs点对应的概率ps数学表达式如下：最后，结合得到的风力发电功率、光伏发电功率的概率密度函数，采用基于wasserstein距离的最优场景生成方法对风电出力和光伏出力分别进行场景生成，分别将两者生成场景数、对应场景概率值进行交叉相乘得到风光联合运行场景及对应场景的概率。步骤二：构建以年网络综合费用最优为主动配电网网架规划的总目标函数，将总目标函数分解为上下两层目标函数，上层以年网络损耗费用、网架投资费用等年值和年dg固定投资费用之和最小为目标函数，下层以年dg运行维护费用、年主动管理费用和年需求侧管理费用之和最小为目标函数，具体为：上层网架规划问题的目标函数如下所示：minc1＝closs+cline+cf；式中，closs为年损耗费用；cline为网架投资费用等年值；cf折算到每年的dg固定投资费。上式中各费用的计算表达式如下所示：closs＝αeloss；式中，α为电价；eloss为年损耗量；fline是起始网架建设费用；r为折现率，取8％；n是线路寿命，架空线路取30，电缆线路取40；rwtg和rpvg分别表示wtg和pvg的现值转等年值系数；nnode为配电网节点数；和为安装在节点i的wtg和pvg的单位容量固定投资成本；pwtg,i和ppvg,i表示安装在节点i的wtg和pvg的额定容量。所述上层网架规划问题的约束条件如下所示：1)、网络辐射状约束：n＝m+1；式中，m是网络的支路数；n是网络的节点数。2)、网络连通性约束：下层配电网运行优化问题的目标函数如下所示：minc2＝com+cam+cdsm；式中，com为dg年运行维护费；cam为dg年主动管理费；cdsm为需求侧管理成本。上式中各费用的计算表达式如下所示：式中，ns为场景数；和为wtg和pvg的单位发电量运行维护成本；pwtg,i,s和ppvg,i,s分别表示安装在节点i的wtg和pvg在场景s的有功出力；和分别为安装在节点i的wtg和pvg的单位发电量的主动管理成本；ndsm表示可以采取需求侧管理的负荷节点数；pdsm,m,s是m节点负荷在场景s的负荷切除量；ξj是m节点负荷单位中断成本。所述下层配电网运行优化问题的约束条件如下所示：1)、wtg在每个场景的运行约束：式中：表示场景s下节点i处wtg的有功出力上限；ωwtg,i为节点i处wtg的有功切除比例；为最大切除比例；为场景s下wtg功角；和分别是节点i处wtg功角的最小和最大值。2)、pvg在每个场景的运行约束：式中：是场景s节点i处pvg的有功出力上限；ωpvg,i是节点i处pvg的有功切除比例；为最大切除比例；为场景s下pvg的功角；和分别是节点i处pvg功角的最小和最大值。3)、可中断负荷中断量在每个场景的约束：式中：表示其最大允许中断量。4)、潮流方程约束：式中：pi,s和qi,s是场景s下注入节点i的有功功率和无功功率；ui,s和uj,s是场景s下节点i和节点j的电压幅值；θij,s是场景s下节点i和j之间的相角差；gij和bij分别是节点i和j之间的互电导和互电纳。5)、每个场景的节点电压约束：uimin≤ui,s≤uimax；式中：uimin和uimax分别为节点i电压的上下限。6)、每个场景的支路功率约束：0≤sk,s≤skmax；式中：skmax为第k条支路上允许传输的最大视在功率。7)、每个场景的dg切除量约束：利用交替迭代方法对上层的网架规划问题和下层的配电网运行优化问题进行求解：即先优化出上层的最优的一组网架结果，然后将最优的网架结果逐个代入到下层对dg出力和可中断负荷量进行优化，然后将下层优化结果传递给上层计算总目标函数，得到最终优化网架结果。所述上下层规划模型之间的传递关系如图1所示。步骤三：对于上层规划模型和下层规划模型分别采用改进的最小生成树prim算法和原对偶内点法进行求解。由于配电网是辐射状的连通性网络，可用图论中的树表示，因此，上层模型采用改进的最小生成树算法求解，算法的核心思路可以阐述如下：1)、对于树状图m＝(g，f)，g是所有顶点全集，f是所有边的全集；2)、令gnew＝{z}，z为集合g中某一个起始节点，令e起始为空集；3)、重复以下的流程，直到gnew＝g：a)、在集合f中选取权值最小的边<m,n>，其中m为集合gnew中的元素，而n不属于gnew，但n属于g，如若存在多条边权值相同，那么取其中任意一条边即可；b)、往集合gnew中加入新元素n点，同时往集合fnew中添加新元素<m,n>边；4)、mnew＝(vnew，enew)即为最终得到的最小生成树。下层模型既是运行优化问题，也可以说是最优潮流问题，本文采用已有文献提出的原对偶内点法进行求解，本发明不进行详细介绍。本文分别采用两种智能算法对上下层规划模型进行求解，规划模型求解流程框图如图2所示，其中，下层规划模型内嵌在上层规划模型之中，最终以上下层目标函数值加和最小为依据确定最优网架方案。下面以改进的一个配电系统作为实施例，验证本发明所提规划方法的有效性和合理性。选择改进的29节点配电网算例，规划区域如图3所示；该规划区域电压等级10kv，规划年限10a，各负荷节点的坐标、有功功率、无功功率期望值以及重要程度见表1；为充分考虑分布式电源的运行特性以及计算方便，在4、5、6节点加装了400kw额定容量的wtg，在9、10节点装了400kw额定容量的pvg；14、15、16节点负荷为可中断负荷；额定风速、切入风速和切出风速分别为12、3.5和20m/s；分布参数和形状参数k、c、α、β分别取2.30、8.92、0.85、0.85；wtg、pvg运行维护费分别设为0.3元/kw·h和0.2元/kw·h，向上级电网购电费用取0.39元/kw·h，售卖电价0.5元/kw·h，主动管理成本0.08元/kw·h，可中断负荷补偿成本0.2元/kw·h，dg最大切除比例20％，可中断负荷最大中断量100％。表1配电网节点数据节点号横坐标/km纵坐标/km有功负荷/kw无功负荷/kw负荷重要程度11.9761.090---21.0561.0261881170.530.4801.3041801120.541.9281.798136840.850.1961.0761841140.663.6400.474160990.570.5240.9141721070.382.8761.8081641020.890.1841.6022441510.3101.0081.5862521560.9110.6641.8221801120.4123.3600.904204990.6130.5480.4302481540.8140.9160.182160990.8153.4241.1921961220.5162.8560.182144890.7172.4880.2721721070.6183.2721.7381881170.2192.8761.5601721070.9203.1121.3942921810.7212.3480.1121921190.9222.1280.3342081290.5233.3000.474136840.5243.4401.4901681040.5252.3041.556118730.2261.1720.35468420.9272.3880.506134830.3282.9441.19677480.5293.6160.71868420.91)、算例结果风力和光伏发电采用前面得到的出力概率密度函数，通过基于wasserstein距离的场景生成方法得到如下25个运行场景见表2。表2风光联合运行场景功率-概率分布情况为了进行对比分析，本文考虑了三种方案进行网架的规划，见表3。方案一，考虑主动管理和需求侧管理措施；方案二，只考虑主动管理不考虑需求侧管理措施(可中断负荷)；方案三，不考虑任何管理措施。表3三种规划方案情况采用基于最小生成树的prim算法和原对偶内点法对方案一、方案二和方案三进行求解，得到最优网架规划方案如图4、图5和图6所示。三种规划方案的各项成本费用对比如表4所示。为了更直观的对比三种规划方案各项费用，将表4中的数据处理成折线图如图7所示。表4三种规划方案各项费用2)结果分析在年网络损耗方面，方案1比方案2少10.49万元，方案2比方案3少13.06万元，这说明采取主动管理和需求侧管理措施之后，可以更好发挥dg在延缓电网投资、减少系统网损以及改善系统运行的积极作用。在向上级购电成本方面，方案1比方案2少44.41万元，方案2比方案3少52.14万元，同时在dg运行维护费上，方案1比方案2多29.19万元，方案2比方案3多14.03万元，这两者说明采取主动管理和需求侧管理之后，配电网络向上级电网购电量大大减少，而dg发电出力更多，证明了主动管理和需求侧管理措施能够促进dg的消纳。在年综合费用方面，方案1＜方案2＜方案3，说明同时采取主动管理和需求侧管理措施进行规划得到的网架方案综合最优。可见，本发明提出的一种基于供需互动和dg运行特性的主动配电网网架规划方法，能够增加dg的消纳，更好发挥dg在减少系统网损、改善系统运行方面的作用。同时，在规划阶段考虑主动管理措施，将运行纳入规划之中，能够科学地得到网架整体最优的规划方案，更加符合未来配电网中dg大规模接入的发展趋势，为未来配电网规划的研究提供一定的参考意义。当前第1页12