本发明涉及一种高比例可再生能源电力系统的骨干网架优化方法，来考虑高比例可再生能源在电力系统核心骨干网架中的作用和对电力系统的影响。基本实施流程包括步骤：输入电力系统的拓扑结构和电气特性参数；计算负荷节点、风电场节点、光伏电站节点、常规机组节点的相对重要度；采用改进的熵权‑理想解法来求取不同类型节点间的相对重要度；采用分阶段骨干网架构建策略，第一阶段模型以优化负荷节点为主要任务，采用蜂群算法求解；第二阶段模型以优化电源节点为主要任务，采用蜂群算法结合贪心策略求解。本发明方法优化得到的核心骨干网架，能够保障对重要负荷的供电，并且具有较高的间歇性能源渗透率与利用率。
一种高比例可再生能源电力系统的骨干网架优化方法
技术领域
[0001] 本发明涉及电力系统的技术领域，更具体地，涉及一种高比例可再生能源电力系统的骨干网架优化方法。
背景技术
[0002] 传统化石能源的过度开采和利用，对人类社会的发展产生很多不利的影响，比如环境污染和气候变化等。可再生能源是全球能源转型和低碳发展的重要解决方案。当前，可再生能源的开发和利用在世界主要国家均处于快速发展阶段，以水能、风能和太阳能为主要代表的可再生能源发电的大规模接入，是未来电力系统发展的必然趋势。我国能源转型与发展的目标是实现2020年和2030年非化石能源分别占一次能源消费比重15％和20％的，并力争在2050年实现可再生能源占比达60％。可见，高比例可再生能源电力系统也是我国电力系统发展和规划的重要研究内容。
[0003] 电力系统的核心骨干网架是指在发生灾害或故障时，在保证电力系统安全运行并保障对重要负荷供电的前提下满足拓扑连通性的支路和节点的集合。电力系统核心骨干网架的研究，对于提高系统的抗灾能力，保障其对重要负荷的供电具有重要意义。目前核心骨干网架的研究均是基于传统的无可再生能源的电力系统而开展的，未考虑高比例可再生能源在电力系统核心骨干网架中的作用和对电力系统的影响。
[0004] 研究高比例可再生能源电力系统的核心骨干网架构建，对提高未来电力系统的稳定性、可靠性和抗灾能力具有重要的意义。然而，间歇性电源如风电、光伏电站等的大规模接入，显著地增加了电力系统的不确定性和复杂程度，导致原本适用于常规电力系统的核心骨干网架构建方法，并不适用于高比例可再生能源电力系统。
[0005] 在此背景下，亟需研究高比例可再生能源电力系统的核心骨干网架优化策略。
发明内容
[0006] 本发明为克服上述现有技术所述的至少一种缺陷，提供一种高比例可再生能源电力系统的骨干网架优化方法。
[0007] 本发明所采用的技术方案如下：
[0008] 高比例可再生能源电力系统的骨干网架优化方法，该方法的实现流程包括以下步骤：
[0009] S1：输入高比例可再生能源电力系统的网络拓扑结构和电气运行参数；
[0010] S2：根据改进的负荷节点的拓扑重要度、负荷节点的电能和潮流通量，计算负荷节点的重要度；
[0011] S3：根据风电场在电力网络中的接近度、等效利用小时数、平均风速和储能配置，计算风电场节点的重要度；
[0012] S4：根据光伏电站在电力网络中的接近度、等效利用小时数和系统效率，计算光伏电站节点的重要度；
[0013] S5：根据常规机组的接近度、机组的备用容量和机组爬坡率，计算常规机组节点的重要度；
[0014] S6：对S2～S5中得到的不同量纲的指标进行标准化，采用改进的熵权-理想解法，来求取不同类型节点间的相对重要度；
[0015] S7：采取两阶段核心骨干网架优化策略，第一阶段模型以优化负荷节点为目标，构建一个包含全部重要负荷节点以及部分支路的连通子网架；
[0016] 根据第一阶段所得子网架，第二阶段模型以优化电源节点为目标，在满足间歇性电源的穿透功率极限约束、备用约束和安全约束的前提下，投入相对重要度超过阈值的各类电源节点并优化得到最终的核心骨干网架。
[0017] 作为优选，采用蜂群算法进行第一阶段模型的优化；采用蜂群算法结合贪心策略求解核心骨干网架优化的第二阶段模型，即在满足约束的前提下按照各类型电源节点的排序逐次加入电源节点，直至目标函数最大或者达到迭代次数为止。
[0018] 作为优选，S2中，根据改进的负荷节点的拓扑重要度、负荷节点的电能和潮流通量来计算负荷节点的相对重要度，具体实现方法如下：
[0019] 设高比例可再生能源电力系统的原网架对应的拓扑图为G0＝(V,E)，其中V＝{V1,2
V2,…,Vn}为节点集，n为节点个数，E为V上的一个子集；节点i的度ki定义为与节点i直接相连的其他节点的数目；聚集系数为描述复杂网络中节点的聚集程度的特征量；计及邻居节点规模的节点拓扑重要度表示为：
[0020]
[0021] 式中：εi为节点i自身的度与其所有邻居节点的度之和，归一化函数θi表示为：
[0022]
[0023] 式中：ci为节点i的聚集系数，表示为：
[0024]
[0025] 式中：ki为节点i的邻居节点个数；si为节点i的ki个邻节点之间实际存在的边数；
[0026] 计及不同时段负荷功率的不同，负荷节点在时段T内的电能表示为：
[0027]
[0028] 式中：PL,i(t)为t时刻负荷节点i的有功功率，PB为功率基准值，SL为负荷节点的集合；
[0029] 潮流通量的定义为正常运行状态下负荷节点的最大视在功率，负荷节点i的潮流通量表示为：
[0030]
[0031] 式中：QL,i(t)为t时刻负荷节点i的无功功率。
[0032] 作为优选，S3中，根据风电场在电力网络中的接近度、等效利用小时数、平均风速和储能配置，来计算风电场节点的相对重要度，具体实现方法如下：
[0033] 对于任意线路i*-j*，节点j*到节点i*的有效长度表示为：
[0034]
[0035] 式中： 为节点j*到节点i*的有效长度， 为节点j*到节点i*之间的边的参数， 表示节点j*与节点i*直接相连；当线路两端为负荷节点i*和j*时，计及负荷节点i*与j*的相对重要度， 表示为：
[0036]
[0037] 式中， 为负荷节点i*的相对重要度， 为负荷节点i*与负荷节点j*之间的线路电抗；
[0038] 结合有效距离，改进后节点i的接近度中心性表示为：
[0039]
[0040] 式中：Dji为任意两节点j到i之间的有效距离；
[0041] 定义风电场节点的接近度中心性为：
[0042]
[0043] 式中：nL为第一阶段网架的负荷节点的数量，SW为风电场节点集合；
[0044] 一定时间T内的风电场等效利用小时数表示为：
[0045]
[0046] 式中：PW,i(t)为风电场i在时刻t的出力， 为该风电场的总装机容量；
[0047] 将一定时间T内风电机组的平均风速作为风电场节点的指标之一，该平均风速表示为：
[0048]
[0049] 式中：vi,t为风电场i在第t个时刻的风速，NT为T时间段内对风速进行测量的次数；
[0050] 将风电场的储能装置指标表示为：
[0051]
[0052] 式中： 为风电场i的储能装置容量，单位为MWh， 为风电场i的储能装置的最大放电功率。
[0053] 作为优选，S4中，根据光伏电站在电力网络中的接近度、等效利用小时数和系统效率，计算光伏电站节点的相对重要度，具体实现方法如下：
[0054] 根据风电节点的接近度中心性，定义光伏节点的接近度中心性为：
[0055]
[0056] 式中：SP为光伏电站节点集合；
[0057] 一定时间T内的光伏电站的等效利用小时数定义为：
[0058]
[0059] 式中：PP,i(t)为光伏电站i的出力关于时间的函数， 为该光伏电站的装机总容量；
[0060] 系统效率表示一段时间内光伏电站的满发小时数与理论发电小时数之比，即光伏电站的实际交流发电量与理论直流发电量之比，表示为：
[0061]
[0062] 式中：Yf为满发小时数，表示一段时间内并网光伏电站的上网电量与光伏电站的额定功率之比；Yr为理论发电小时数，表示单位面积的光伏阵列倾斜面总辐射量与光伏电池在标准测试条件下的标准辐照度之比。
[0063] 作为优选，S5中，根据常规机组的接近度、机组的备用容量和机组爬坡率，计算常规机组节点的相对重要度，具体实现方法如下：
[0064] 根据风电节点的接近度中心性，定义常规机组节点的接近度中心性为：
[0065]
[0066] 式中：SC为常规机组节点集合；
[0067] 用常规机组的最大旋转备用容量代替其容量大小作为常规机组的相对重要度的一个评价指标，将常规机组的最大旋转备用容量表示为：
[0068]
[0069] 式中： 为常规机组i的额定功率或最大技术出力， 为常规机组i的最小技术出力；
[0070] 以火电机组爬坡率的大小反映提供旋转备用的快慢程度：
[0071]
[0072] 式中：Ri为火电机组i的爬坡率，ST为火电机组集合，
[0073] 作为优选，S6中，采用改进的熵权-理想解法对节点进行相对重要度评价，具体实现方法如下：
[0074] 步骤1：对不同量纲的指标进行标准化，设节点i的第m个指标为r'im，则正向指标标准化表示为：
[0075]
[0076] 式中：rim为标准化后的节点i的第m个指标，M为指标个数，NN为节点的个数；逆向指标的标准化与正向指标标准化的形式相同，记标准化后的指标矩阵为
[0077]
[0078] 步骤2：根据标准化后的数据求取第m个指标的熵：
[0079]
[0080] 步骤3：根据熵值求得第m个指标的熵权为：
[0081]
[0082] 记指标熵权向量为w＝[w1,w2,…,wM]；
[0083] 步骤4：根据指标熵权和标准化指标矩阵构建加权指标矩阵 其中，zim＝wmrim；
[0084] 步骤5：定义“理想解”z+与“负理想解”z-，其中
[0085]
[0086] 式中： 其中m＝1,2,…,M；
[0087] 步骤6：计算节点i与“理想解”的欧式距离 以及其与“负理想解”的欧式距离其中：
[0088]
[0089] 步骤7：计算节点i的贴近度，贴近度表示节点i远离“负理想解”并靠近“理想解”的程度，贴近度越大，节点越重要；节点i的贴近度表示为：
[0090]
[0091] 步骤8：按照各个节点与理想解的贴近度大小，对负荷节点和电源节点进行重要性程度排序，从而得到节点重要度的评估结果；此外，为了使网架优化模型满足尽可能多地投入间歇性电源节点的要求，将所有间歇性电源节点的重要度设定为Δi＝Const+Δi，其中Const为一常数，i∈SW∪SP。
[0092] 作为优选，S7中，两阶段核心骨干网架优化策略的第一阶段模型以优化负荷节点为目标，构建一个包含全部重要负荷节点以及部分支路的连通子网架，具体实现方法如下：
[0093] 将第一阶段的网架优化模型表示为：
[0094]
[0095] 式中：(i,j)∈G0表示线路i-j属于网架G0；αi-j为线路i-j的投入状态，当网架G包含线路i-j时，αi-j＝1，否则αi-j＝0；βk为负荷节点k的切除状态，当网架G包含负荷节点k时，βk＝0，否则βk＝1；λ1为折衷系数；为了满足连通性约束，则第一阶段的网架优化模型表示为：
[0096]
[0097] 式中：IC(G)为指示函数，当网架G连通时，IC(G)＝0；否则IC(G)＝+∞。
[0098] 作为优选，第二阶段模型以优化电源节点为目标，投入相对重要度超过阈值的各类电源节点并优化得到最终的核心骨干网架，具体实现方法如下：
[0099] 在投入电源节点的同时，兼顾电源节点的位置，即考虑电源节点的投入所导致的线路数目的增加量；构造第二阶段核心骨干网架优化的目标函数为：
[0100]
[0101] 式中：p∈G0表示电源节点p属于原网架G0；FS,p为电源节点p的重要度；γp为电源节点的投入状态，若电源节点被选入核心骨干网架，则γp＝1，否则，γp＝0；λ2为线路权重调整系数，用于调节对电源节点到G的线路数量的重视程度；lp-G为电源节点p到子网架G的最少线路数；
[0102] 第二阶段的网架优化应满足如下约束
[0103] 1)风电的穿透功率极限约束：
[0104]
[0105] 式中：SWF和SLF分别为核心骨干网架的风电场节点集合和负荷节点集合，γW％为核心骨干网架中风电的穿透功率极限，γW0％为原网架中风电的穿透功率极限；
[0106] 2)光伏发电的穿透功率极限约束：
[0107]
[0108] 式中：SPF为核心骨干网架的光伏电站节点集合，γP％为核心骨干网架的光伏发电的穿透功率极限，γP0％为原网架的光伏发电的穿透功率极限；
[0109] 3)备用约束：
[0110]
[0111] 式中：SCF为核心骨干网架的常规机组节点集合，PR(t)为t时刻系统的负荷备用功率；
[0112] 4)等式约束g(x)和不等式约束h(x)统一表示为：
[0113]
[0114] 式中：x为状态变量向量，包括节点电压、发电机功率。
[0115] 该方法优化得到的核心骨干网架能保障对重要负荷的供电，并且其具有较高的间歇性能源渗透率与利用率。
附图说明
[0116] 图1是本发明的整体流程示意图。
[0117] 图2是本发明的IEEE-118节点系统图。
[0118] 图3是本发明的夏季典型场景的源荷曲线。
[0119] 图4是本发明的IEEE-118系统重要负荷分布图。
[0120] 图5是本发明的IEEE-118节点系统的骨干网架第一阶段网架优化结果。
[0121] 图6是本发明的IEEE-118节点系统的骨干网架第二阶段网架优化结果。
具体实施方式
[0122] 附图仅用于示例性说明，不能理解为对本专利的限制；对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。附图中描述位置关系仅用于示例性说明，不能理解为对本专利的限制。
[0123] 高比例可再生能源电力系统的核心骨干网架具有以下特征：
[0124] 1)不中断对重要负荷的供电；
[0125] 2)核心骨干网架的拓扑结构配置合理，并满足网络的拓扑连通性约束；
[0126] 3)在可再生能源的任何出力状况下，核心骨干网架均能满足高比例可再生能源电力系统的安全运行要求；
[0127] 4)在满足1)到3)的情况下，间歇性电源渗透率和利用率尽可能高；
[0128] 5)在满足1)到4)的情况下，核心骨干网架支路总数尽可能少。
[0129] 本发明的高比例可再生能源电力系统的骨干网架优化方法，其实现流程包括以下步骤：
[0130] S1：输入高比例可再生能源电力系统的网络拓扑结构和电气运行参数。
[0131] S2：根据改进的负荷节点的拓扑重要度、负荷节点的电能和潮流通量，计算负荷节点的重要度。本步骤的具体实现方法如下：
[0132] 设高比例可再生能源电力系统的原网架对应的拓扑图为G0＝(V,E)，其中V＝{V1,2
V2,…,Vn}为节点集，n为节点个数，E为V上的一个子集；节点i的度ki定义为与节点i直接相连的其他节点的数目；聚集系数为描述复杂网络中节点的聚集程度的特征量；计及邻居节点规模的节点拓扑重要度表示为：
[0133]
[0134] 式中：εi为节点i自身的度与其所有邻居节点的度之和，归一化函数θi表示为：
[0135]
[0136] 式中：ci为节点i的聚集系数，表示为：
[0137]
[0138] 式中：ki为节点i的邻居节点个数；si为节点i的ki个邻节点之间实际存在的边数；
[0139] 计及不同时段负荷功率的不同，负荷节点在时段T内的电能表示为：
[0140]
[0141] 式中：PL,i(t)为t时刻负荷节点i的有功功率，PB为功率基准值，SL为负荷节点的集合；
[0142] 潮流通量的定义为正常运行状态下负荷节点的最大视在功率，负荷节点i的潮流通量表示为：
[0143]
[0144] 式中：QL,i(t)为t时刻负荷节点i的无功功率。
[0145] S3：根据风电场在电力网络中的接近度、等效利用小时数、平均风速和储能配置，计算风电场节点的重要度。本步骤的具体实现方法如下：
[0146] 对于任意线路i*-j*，节点j*到节点i*的有效长度表示为：
[0147]
[0148] 式中： 为节点j*到节点i*的有效长度， 为节点j*到节点i*之间的边的参数， 表示节点j*与节点i*直接相连；当线路两端为负荷节点i*和j*时，计及负荷节点i*与j*的相对重要度， 表示为：
[0149]
[0150] 式中， 为负荷节点i*的相对重要度， 为负荷节点i*与负荷节点j*之间的线路电抗；
[0151] 结合有效距离，改进后节点i的接近度中心性表示为：
[0152]
[0153] 式中：Dji为任意两节点j到i之间的有效距离；
[0154] 定义风电场节点的接近度中心性为：
[0155]
[0156] 式中：nL为第一阶段网架的负荷节点的数量，SW为风电场节点集合；
[0157] 一定时间T内的风电场等效利用小时数表示为：
[0158]
[0159] 式中：PW,i(t)为风电场i在时刻t的出力， 为该风电场的总装机容量；
[0160] 将一定时间T内风电机组的平均风速作为风电场节点的指标之一，该平均风速表示为：
[0161]
[0162] 式中：vi,t为风电场i在第t个时刻的风速，NT为T时间段内对风速进行测量的次数；
[0163] 将风电场的储能装置指标表示为：
[0164]
[0165] 式中： 为风电场i的储能装置容量，单位为MWh， 为风电场i的储能装置的最大放电功率。
[0166] S4：根据光伏电站在电力网络中的接近度、等效利用小时数和系统效率，计算光伏电站节点的重要度。本步骤的具体实现方法如下：
[0167] 根据风电节点的接近度中心性，同理定义光伏节点的接近度中心性为：
[0168]
[0169] 式中：SP为光伏电站节点集合；
[0170] 一定时间T内的光伏电站的等效利用小时数定义为：
[0171]
[0172] 式中：PP,i(t)为光伏电站i的出力关于时间的函数， 为该光伏电站的装机总容量；
[0173] 系统效率表示一段时间内光伏电站的满发小时数与理论发电小时数之比，即光伏电站的实际交流发电量与理论直流发电量之比，表示为：
[0174]
[0175] 式中：Yf为满发小时数，表示一段时间内并网光伏电站的上网电量与光伏电站的额定功率之比；Yr为理论发电小时数，表示单位面积的光伏阵列倾斜面总辐射量与光伏电池在标准测试条件下的标准辐照度之比。
[0176] S5：根据常规机组的接近度、机组的备用容量和机组爬坡率，计算常规机组节点的重要度。本步骤的具体实现方法如下：
[0177] 根据风电节点的接近度中心性，定义常规机组节点的接近度中心性为：
[0178]
[0179] 式中：SC为常规机组节点集合；
[0180] 用常规机组的最大旋转备用容量代替其容量大小作为常规机组的相对重要度的一个评价指标，将常规机组的最大旋转备用容量表示为：
[0181]
[0182] 式中： 为常规机组i的额定功率或最大技术出力， 为常规机组i的最小技术出力；
[0183] 以火电机组爬坡率的大小反映提供旋转备用的快慢程度：
[0184]
[0185] 式中：Ri为火电机组i的爬坡率，ST为火电机组集合，
[0186] S6：S2～S5中计算得到的各指标即对应其重要度，对S2～S5中得到的不同量纲的指标进行标准化，采用改进的熵权-理想解法，来求取不同类型节点间的相对重要度，改进的熵权-理想解法对节点进行相对重要度评价的具体实现方法如下：
[0187] 步骤1：对不同量纲的指标进行标准化，设节点i的第m个指标为r'im，则正向指标标准化表示为：
[0188]
[0189] 式中：rim为标准化后的节点i的第m个指标，M为指标个数，NN为节点的个数；逆向指标的标准化与正向指标标准化的形式相同，记标准化后的指标矩阵为
[0190] 步骤2：根据标准化后的数据求取第m个指标的熵：
[0191]
[0192] 步骤3：根据熵值求得第m个指标的熵权为：
[0193]
[0194] 记指标熵权向量为w＝[w1,w2,…,wM]；
[0195] 步骤4：根据指标熵权和标准化指标矩阵构建加权指标矩阵 其中，zim＝wmrim；
[0196] 步骤5：定义“理想解”z+与“负理想解”z-，其中
[0197]
[0198] 式中： 其中m＝1,2,…,M；
[0199] 步骤6：计算节点i与“理想解”的欧式距离 以及其与“负理想解”的欧式距离其中：
[0200]
[0201] 步骤7：计算节点i的贴近度，贴近度表示节点i远离“负理想解”并靠近“理想解”的程度，贴近度越大，节点越重要；节点i的贴近度表示为：
[0202]
[0203] 步骤8：按照各个节点与理想解的贴近度大小，对负荷节点和电源节点进行重要性程度排序，从而得到节点重要度的评估结果；此外，为了使网架优化模型满足尽可能多地投入间歇性电源节点的要求，将所有间歇性电源节点的重要度设定为Δi＝Const+Δi，其中Const为一常数，i∈SW∪SP。
[0204] S7：采取两阶段核心骨干网架优化策略，第一阶段模型以优化负荷节点为目标，构建一个包含全部重要负荷节点以及部分支路的连通子网架，具体实现方法如下：
[0205] 将第一阶段的网架优化模型表示为：
[0206]
[0207] 式中：(i,j)∈G0表示线路i-j属于网架G0；αi-j为线路i-j的投入状态，当网架G包含线路i-j时，αi-j＝1，否则αi-j＝0；βk为负荷节点k的切除状态，当网架G包含负荷节点k时，βk＝0，否则βk＝1；λ1为折衷系数；为了满足连通性约束，则第一阶段的网架优化模型表示为：
[0208]
[0209] 式中：IC(G)为指示函数，当网架G连通时，IC(G)＝0；否则IC(G)＝+∞。
[0210] 根据第一阶段所得子网架，第二阶段模型以优化电源节点为目标，在满足间歇性电源的穿透功率极限约束、备用约束和安全约束的前提下，投入相对重要度超过阈值的各类电源节点并优化得到最终的核心骨干网架，具体实现方法如下：
[0211] 在投入电源节点的同时，兼顾电源节点的位置，即考虑电源节点的投入所导致的线路数目的增加量；构造第二阶段核心骨干网架优化的目标函数为：
[0212]
[0213] 式中：p∈G0表示电源节点p属于原网架G0；FS,p为电源节点p的重要度；γp为电源节点的投入状态，若电源节点被选入核心骨干网架，则γp＝1，否则，γp＝0；λ2为线路权重调整系数，用于调节对电源节点到G的线路数量的重视程度；lp-G为电源节点p到子网架G的最少线路数；
[0214] 第二阶段的网架优化应满足如下约束
[0215] 1)风电的穿透功率极限约束：
[0216]
[0217] 式中：SWF和SLF分别为核心骨干网架的风电场节点集合和负荷节点集合，γW％为核心骨干网架中风电的穿透功率极限，γW0％为原网架中风电的穿透功率极限；
[0218] 2)光伏发电的穿透功率极限约束：
[0219]
[0220] 式中：SPF为核心骨干网架的光伏电站节点集合，γP％为核心骨干网架的光伏发电的穿透功率极限，γP0％为原网架的光伏发电的穿透功率极限；
[0221] 3)备用约束：
[0222]
[0223] 式中：SCF为核心骨干网架的常规机组节点集合，PR(t)为t时刻系统的负荷备用功率；
[0224] 4)等式约束g(x)(包括潮流约束、功率平衡约束)和不等式约束h(x)(包括发电机出力约束、线路热稳定极限约束等)统一表示为：
[0225]
[0226] 式中：x为状态变量向量，包括节点电压、发电机功率。
[0227] 采用蜂群算法进行第一阶段模型的优化；采用蜂群算法结合贪心策略求解核心骨干网架优化的第二阶段模型，即在满足约束的前提下按照各类型电源节点的排序逐次加入电源节点，直至目标函数最大或者达到迭代次数为止。
[0228] 实施例
[0229] 如图1所示，一种高比例可再生能源电力系统的骨干网架优化方法，其基本实施流程包括如下详细步骤：
[0230] 步骤1、输入高比例可再生能源电力系统的拓扑结构参数(线路连接关系)、电气特性参数(节点参数、线路参数)；
[0231] 步骤2、计算负荷节点的各项参数值，具体参数包括：
[0232] 1)改进拓扑重要度：
[0233]
[0234]
[0235] 式中：εi为节点i自身的度与其所有邻居节点的度之和，θi为归一化函数[0236] 2)节点的电能：
[0237]
[0238] 式中：PL,i(t)为t时刻负荷节点i的有功功率，PB为功率基准值，SL为负荷节点的集合。
[0239] 3)潮流通量：
[0240]
[0241] 式中：QL,i(t)为t时刻负荷节点i的无功功率。
[0242] 步骤3、计算风电场节点的各项参数值，具体参数包括：
[0243] 1)电力网络中的接近度：
[0244]
[0245]
[0246]
[0247]
[0248] 式中：Dij为任意两节点j到i之间的有效距离，h0为构成节点j和节点i之间有向路径的中间节点。 为负荷节点i*的相对重要度， 为负荷节点i*与负荷节点j*之间的线路电抗； 为节点j*到节点i*的有效长度， 为节点j*到节点i*之间的边的参数，表示节点j*与节点i*直接相连。
[0249] 2)等效利用小时数：
[0250]
[0251] 式中：PW,i(t)为风电场i在时刻t的出力， 为该风电场的总装机容量。
[0252] 3)平均风速：
[0253]
[0254] 式中：vi,t为风电场i在第t个时刻的风速，NT为T时间段内对风速进行测量的次数。
[0255] 4)储能配置：
[0256]
[0257] 式中： 为风电场i的储能装置容量，单位为MWh， 为风电场i的储能装置的最大放电功率。
[0258] 步骤4、计算光伏电站节点的各项参数值，具体参数包括：
[0259] 1)电力网络中的接近度：
[0260]
[0261] 式中：SP为光伏电站节点集合。
[0262] 2)等效利用小时数：
[0263]
[0264] 式中：PP,i(t)为光伏电站i的出力关于时间的函数， 为该光伏电站的装机总容量。
[0265] 3)系统效率：
[0266]
[0267] 式中：Yf为满发小时数，表示一段时间内并网光伏电站的上网电量与光伏电站的额定功率之比；Yr为理论发电小时数，表示单位面积的光伏阵列倾斜面总辐射量与光伏电池在标准测试条件下的标准辐照度之比。
[0268] 步骤5、计算常规机组节点的各项参数值，具体参数包括：
[0269] 1)电力网络中的接近度：
[0270]
[0271] 式中：SC为常规机组节点集合。
[0272] 2)机组的备用容量：
[0273]
[0274] 式中： 为常规机组i的额定功率(或最大技术出力，本发明认为最大技术出力等于额定功率)， 为常规机组i的最小技术出力。
[0275] 3)爬坡率：
[0276]
[0277] 式中：Ri为火电机组i的爬坡率，ST为火电机组集合，
[0278] 步骤6、对不同量纲的指标进行标准化：设节点i的第m个指标为r'im，则正向指标标准化可表示为
[0279]
[0280] 式中：rim为标准化后的节点i的第m个指标，M为指标个数，NN为节点的个数。逆向指标的标准化与正向指标标准化的形式相同，记标准化后的指标矩阵为
[0281]
[0282] 步骤7、根据标准化后的数据求取第m个指标的熵权，记指标熵权向量为w＝[w1,w2,…,wM]：
[0283]
[0284]
[0285] 步骤8、根据指标熵权和标准化指标矩阵构建加权指标矩阵 其中，zim＝wmrim。
[0286] 步骤9、步骤6：计算节点i与“理想解”的欧式距离 以及其与“负理想解”的欧式距离
[0287]
[0288] 其中，z+为“理想解”，z-为“负理想解”：
[0289]
[0290] 式中： 其中m＝1,2,…,M。
[0291] 步骤10、计算节点i的贴近度，按照各个节点与理想解的贴近度大小，对负荷节点和电源节点进行重要性程度排序，从而得到节点重要度的评估结果。节点i的贴近度为：
[0292]
[0293] 步骤11、以优化负荷节点为主要任务，采用蜂群算法进行核心骨干网架的第一阶段模型优化，构建一个包含全部重要负荷节点以及部分支路的连通子网架。模型的目标函数为：
[0294]
[0295] 式中：(i,j)∈G0表示线路i-j属于网架G0；αi-j为线路i-j的投入状态，当网架G包含线路i-j时，αi-j＝1，否则αi-j＝0；βk为负荷节点k的切除状态，当网架G包含负荷节点k时，βk＝0，否则βk＝1；λ1为折衷系数，表示对节点重要度和线路数目的折衷程度，考虑节点相对重要度的大小，λ1一般在[2,5]之间选取。
[0296] 若采用指示函数表示网架的连通性，则第一阶段的网架优化模型可表示为[0297]
[0298] 式中：IC(G)为指示函数，当网架G连通时，IC(G)＝0；否则IC(G)＝+∞。
[0299] 步骤12、以优化电源节点为主要任务，采用蜂群算法结合贪心策略求解第二阶段模型，投入相对重要度较大的各类电源节点并优化得到最终的核心骨干网架。模型的目标函数为：
[0300]
[0301] 式中：p∈G0表示电源节点p属于原网架G0；FS,p为电源节点p的重要度；γp为电源节点的投入状态，若电源节点被选入核心骨干网架，则γp＝1，否则，γp＝0；λ2为线路权重调整系数，lp-G为电源节点p到子网架G的最少线路数。
[0302] 第二阶段的网架优化应满足如下约束
[0303] 1)风电的穿透功率极限约束：
[0304]
[0305] 式中：SWF和SLF分别为核心骨干网架的风电场节点集合和负荷节点集合，γW％为核心骨干网架中风电的穿透功率极限，γW0％为原网架中风电的穿透功率极限；
[0306] 2)光伏发电的穿透功率极限约束：
[0307]
[0308] 式中：SPF为核心骨干网架的光伏电站节点集合，γP％为核心骨干网架的光伏发电的穿透功率极限，γP0％为原网架的光伏发电的穿透功率极限；
[0309] 3)备用约束：
[0310]
[0311] 式中：SCF为核心骨干网架的常规机组节点集合，PR(t)为t时刻系统的负荷备用功率；
[0312] 4)等式约束g(x)和不等式约束h(x)统一表示为：
[0313]
[0314] 式中：x为状态变量向量，包括节点电压、发电机功率。
[0315] 为了进一步理解本发明，以下以修改后的IEEE-118节点为例，来解释本发明的实际应用。修改后的IEEE-118节点包括18个火电机组，2个水电机组，14个风电场，15个光伏电站，系统的拓扑结构如附图2所示。源荷参数如表1至表5所示。
[0316] 表1负荷参数及其节点相对重要度
[0317]
[0318]
[0319]
[0320]
[0321] 表2风电场参数及其节点相对重要度
[0322]
[0323]
[0324] 表3光伏电站参数及其节点相对重要度
[0325]
[0326] 表4常规机组参数及其节点相对重要度
[0327]
[0328]
[0329] 注：带*号的为水电机组，如4*和85*表示节点4和节点85为水电机组节点。
[0330] 表5夏季典型场景日各小时负荷
[0331]
[0332] 注：各节点负荷随总负荷变化的比例而变化
[0333] 附图3所示为夏季某典型场景下系统的负荷、风电场出力、光伏电站出力的有功功率曲线，本发明以夏季典型场景为例说明所提出的核心骨干网架优化策略。重要负荷节点和水电机组节点在核心骨干网架中必须予以保留，其在拓扑图中的位置如附图4所示。原网架的风电和光伏的穿透功率极限分别为35％和10％，核心骨干网架的风电和光伏的穿透功率极限分别为25％和7％。折衷系数λ1＝4，线路权重调整系数λ2＝0.5。负荷备用功率取为负荷功率的0.2倍。蜂群算法求解时，其种群数量为30，迭代次数为100次。
[0334] 根据IEEE-118节点的网架参数和夏季典型场景的源荷参数可求得各类型节点的相对重要度，表6所示为部分节点的相对重要度，完整的节点相对重要度结果如表1至表4所示。表6中，节点编号由两部分组成，即“字母+数字”，字母L、W、P和T分别表示负荷节点、风电场节点、光伏电站节点和火电机组节点，数字为图2中的节点编号。RL表示相对重要度。
[0335] 表6部分节点相对重要度
[0336]
[0337]
[0338] 根据表6的节点相对重要度结果，采用本发明方法优化得到如附录图5所示的第一阶段子网架，该网架含70条线路、65个节点。第一阶段网架优化结果如表7所示。由于本发明第一阶段优化模型的目标函数同时考虑了网架的线路数目和节点重要度，其意义更加明确，且满足核心骨干网架的结构配置合理的要求，有利于依靠核心骨干网架恢复灾害后的其余负荷节点和电源节点。因此，本发明提出的高比例可再生能源电力系统的第一阶段核心骨干网架优化策略更加合理，其能够保证所得到的网架的线路数最少且负荷节点的相对重要度最大。
[0339] 表7第一阶段网架优化结果
[0340]
[0341] 根据电源节点相对重要度和核心骨干网架的第二阶段优化模型可得到最终的核心骨干网架如图6所示。图6中，对角线条纹阴影的图元和菱形阴影图元代表核心骨干网架的重要负荷或电源，加粗虚线的边代表核心骨干网架中的线路，部分线路连接的节点并未标红，如节点30和节点37等，这里称之为冗余节点。冗余节点虽然保留在核心骨干网架中，但是并不保障不中断对其负荷的供电，冗余节点主要承担着连接网架的作用。当高比例可再生能源电力系统的功率充足时，可向冗余节点的负荷供电。此外，由于在核心骨干网架优化时充分考虑到了节点相对重要度信息，将冗余节点保留在核心骨干网架中亦有利于灾害过后依靠核心骨干网架恢复整个网络。
[0342] 在本发明所构建的夏季典型场景的高比例可再生能源电力系统核心骨干网架中，风电、光伏的最大功率占比分别为21.56％和6.71％，其电能占比分别为11.18％和2.71％。可见，所构建的高比例可再生能源电力系统的核心骨干网架具有较高的可再生能源渗透率和利用率，进一步验证了本发明两阶段网架优化策略的合理性与有效性。
[0343] 显然，上述实施例仅是为清楚地说明本发明而所作的举例，并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。
