一种N-1下新能源渗透率的确定方法

一种n-1下新能源渗透率的确定方法
技术领域
1.本发明涉及一种n-1下新能源渗透率的确定方法，属于新能源渗透率评估技术领域。


背景技术：

2.新能源并网是当前关注的热点，在西北、西南等地区拥有丰富的风力、光伏等资源，近年来加大了对新能源的扶持力度，风力发电和光伏发电也得到了快速发展。然而，新能源资源分布不均，需要通过电网进行大容量、长距离的输送，同时也带来了一系列稳定安全性问题。有效评估大规模新能源消纳能力是提高新能源消纳水平的有效途径。新能源渗透率是限制新型电力系统调节能力和抗扰动性能的重要界限，为了提高风光资源利用效率，降低大规模新能源接入对系统运行的影响，准确评估系统可承载的新能源极限规模对于判断新能源发展的规模边界和合理增速具有重要意义。
3.针对大规模新能源并网所呈现的明显地域性，以及新能源临界渗透率受多类稳定问题约束，本发明将多区域新能源进行聚合，以负载率作为指标快速搜寻最严重的n-1状态，再以新能源消纳率最大值作为优化目标，刻画n-1静态安全域，将该静态安全域映射至新能源出力的低维空间，可获取最大新能源消纳能力曲线，能够直接、快速、准确获取在n-1故障或检修状态下的最大性能源渗透率的计算结果。


技术实现要素：

4.本发明要解决的技术问题是提供一种n-1下新能源渗透率的确定方法，用以大规模新能源并网并考虑n-1故障或检修状态下的最大新能源渗透率进行评估。
5.本发明的技术方案是：一种n-1下新能源渗透率的确定方法，具体步骤为：
6.step1：对新能源多区域进行聚合，将新能源分为两区域或三区域。
7.所述step1具体步骤为：
8.step1.1：基于k-mean聚类将多区域新能源场站进行聚类。
9.step1.2：根据聚类结果，按照容量加权法，将多区域新能源场站进行聚合，等值为两区域或三区域。
10.所述step1.1具体步骤为：
11.step1.1.1：在待分类样本中随机选择k个样本数据作为初始凝聚点，每个凝聚点代表1类。
12.step1.1.2：计算其余样本与凝聚点之间的欧式距离，将每个样本归入与它欧式距离最小的凝聚点所在的类。
13.step1.1.3：所有数据归类完成后，将每一类的凝聚点更新为该类中当前数据的均值。
14.step1.1.4：重复step1.1.2，直至前后2次凝聚点不再变化，即实现了数据的分类。
15.step2：对新能源接入的电力系统发生n-1故障或检修的所有情况进行潮流计算。
16.所述step2中潮流计算的约束条件具体为：
[0017][0018][0019][0020][0021][0022][0023][0024]
式中，v
gk
、p
gk
、q
gk
分别表示第k台发电机的机端电压、有功功率、无功功率，ui、ii、pi、qi分别表示第i个节点的电压幅值、电流幅值、有功功率及无功功率，下标min、max表示对应参量的上限、下限。
[0025]
所述step2中的约束条件，除基本约束条件外，还包括：
[0026]
新能源出力约束，具体如下所示：
[0027][0028]
其中，为新能源的最大出力值。
[0029]
n-1断线输电线路传输功率约束，具体表示如下：
[0030][0031]
其中，表示系统发生n-1后仍运行线路的传输功率，表示对应线路的最大允许传输值，正负号表示传输方向。
[0032]
step3：根据输电线路的负载率，得到系统的不平衡度，以此作为最严重n-1的选择标准并对其进行选择。
[0033]
所述step3具体步骤为：
[0034]
step3.1：定义输电线路的负载率为：
[0035][0036]
其中，r
ij
表示节点i、j所相连输电线路的负载率，p
lij
为对应线路的实际功率传输值，为基态时对应线路所能允许功率传输的最大值。
[0037]
step3.2：采用系统的不平衡度作为表征n-1下严重程度的指标，记为r，根据step3.1求取出不同n-1情况下的各条输电线路负载率，进行求和并计算出输电线路平均负载率。继续计算出各条输电线路与平均负载率的绝对偏差，最终得到不同n-1情况下系统的不平衡度，具体如下所示：
[0038]
[0039]
其中，l表示两节点所相连的输电线路集合，n
l
表示输电线路数量。
[0040]
step3.3：根据不同n-1情况下的系统不平衡度，最大系统不平衡度即为最严重的n-1情况，映射关系如下：
[0041][0042]
其中，c表示与最大系统不平衡度相对应的最严重故障场景。
[0043]
step4：根据最严重的n-1状态，设置新能源渗透率最大值作为优化目标，实现n-1静态安全域的刻画。
[0044]
所述step4具体步骤为：
[0045]
step4.1：根据step3在不同情况下的n-1状态下计算系统不平衡度，根据最大系统不平衡度可以确定系统最严重的n-1情况，以最严重n-1状态作为系统安全运行约束，在此基础上，设置新能源渗透率最大值作为优化目标，具体如下所示：
[0046][0047]
其中，为第i个n-1静态安全域边界点所对应的新能源出力增长系数，为第i个n-1静态安全域边界点所对应的新能源功率初始值，为系统总出力，由新能源出力与传统机组出力共同构成。
[0048]
step4.2：确定新能源聚合后的初始功率与注入方向，如新能源功率沿着某方向增长，根据step4.1中的目标函数，以及step2中的约束条件，基于最优潮流的方法，求解获得，从而可得到新能源渗透率最大临界值；
[0049]
step4.3：更新新能源的注入方向，根据step4.2中的方法继续求解不同新能源区域组合下的临界点，直至达到两区域或三区域的坐标边界，连接所有临界点，形成能够表征新能源最大渗透率的n-1静态安全域。
[0050]
step5：将静态安全域映射至新能源出力的低维空间，得到系统的最大消纳能力曲线，实现系统n-1下的满足系统运行安全约束的最大新能源渗透率的计算：
[0051]
设两区域i、j的新能源初始功率分别为、，将其映射至二维坐标轴的空间内，分别对应着横轴与竖轴的坐标点、，根据坐标点对应的全系统新能源功率注入方向，与坐标点（，）的关系为：
[0052][0053]
式中，为功率增长角，是与横轴的夹角，、分别为区域i、j的新能源出力值。
[0054]
从而将静态安全域映射至新能源出力的二维或三维空间，得到系统的最大消纳能力曲线，新能源消纳能力曲线维度与聚合区域有关，两区域为二维关系，三区域为三维关
系，最大消纳能力曲线坐标轴信息均表示不同区域下新能源的出力情况，计算系统的新能源总出力情况，所得新能源总出力情况与系统总功率的比值即为该系统所能允许的新能源最大渗透率，表达为：
[0055][0056]
式中，为系统所能允许的新能源最大渗透率，为新能源总出力，为系统总出力。
[0057]
现有技术对于大规模新能源并网的承载能力评估，主要关注于规划层面的中长期时段，而从短期或超短期的调度角度来看，需要引入逐段降低负荷、灵活运行调度等方法，以提高电网承载新能源的能力。同时，考虑到新能源产量的不确定性以及外部因素对电网的影响，对于大规模新能源并网所呈现的明显地域性，以及新能源临界渗透率受多类稳定问题约束，本发明将大规模新能源集中接入的多区域新能源进行聚合，以负载率作为指标快速搜寻最严重的n-1状态，在各类电网安全约束下以新能源消纳率最大值作为优化目标，刻画n-1静态安全域，将该静态安全域映射至新能源出力的低维空间，可获取最大新能源消纳能力曲线，最大消纳能力曲线坐标轴信息均表示不同区域下新能源的出力情况，通过多区域不同场景下新能源的出力情况能够直接、快速、准确获取在n-1故障或检修状态下最大新能源渗透率的计算结果，具有较高的实际应用价值。
[0058]
本发明的有益效果是：
[0059]
1、本发明不需要对所有的n-1情况都进行静态安全域的刻画，即可求出满足所有安全要求的最大渗透率。
[0060]
2、本发明适用于任何新能源接入的电力系统，不受电网结构限制。
[0061]
3、本发明可实现n-1的快速搜寻，在保证计算精度的同时提高计算速度。
[0062]
4、本发明计及n-1故障或检修状态等电网安全约束，更加适用于对安全性要求更高的调度系统。
[0063]
5、本发明可实现大规模新能源集中接入的多区域不同场景下的全天不同时段新能源渗透率计算，具有较高的实际应用价值。
附图说明
[0064]
图1是本发明提供的刻画新能源最大消纳能力的n-1静态安全域方法流程图；
[0065]
图2是本发明提供的确定二维最大新能源出力边界方法示意图。
具体实施方式
[0066]
下面结合附图和具体实施方式，对本发明作进一步说明。
[0067]
实施例1：如图1所示，一种n-1下新能源渗透率的确定方法，具体步骤为：
[0068]
step1：对新能源多区域进行聚合，将新能源分为两区域或三区域。
[0069]
具体的，基于k-mean聚类将多区域新能源场站进行聚类，根据聚类结果，按照容量加权法，将多区域新能源场站进行聚合，等值为两区域或三区域。
[0070]
step2：对新能源接入的电力系统发生n-1故障或检修的所有情况进行潮流计算。
[0071]
具体的，潮流计算的约束条件具体如下所示：
[0072][0073][0074][0075][0076][0077][0078][0079]
上式中，v
gk
、p
gk
、q
gk
分别表示第k台发电机的机端电压、有功功率、无功功率，ui、ii、pi、qi分别表示第i个节点的电压幅值、电流幅值、有功功率及无功功率。下标min、max表示对应参量的上限、下限。
[0080]
step3：根据输电线路的负载率，得到系统的不平衡度，以此作为最严重n-1的选择标准并对其进行选择。
[0081]
具体的，定义输电线路的负载率为：
[0082][0083]
其中，r
ij
表示节点i、j节点所相连输电线路的负载率，p
lij
为对应线路的实际功率传输值，为基态时对应线路所能允许功率传输的最大值。
[0084]
采用系统的不平衡度作为表征n-1下严重程度的指标，记为r。根据步骤step3求取出不同n-1情况下的各条输电线路负载率，进行求和并计算出输电线路平均负载率；继续计算出各条输电线路与平均负载率的绝对偏差，最终得到不同n-1情况下系统的不平衡度，具体如下所示：
[0085][0086]
其中，l表示两节点所相连的输电线路集合，n
l
表示输电线路数量。
[0087]
根据不同n-1情况下的系统不平衡度，最大系统不平衡度即为最严重的n-1情况，映射关系如下：
[0088][0089]
其中，c表示与最大系统不平衡度相对应的最严重故障场景。
[0090]
step4：根据最严重的n-1状态，设置新能源渗透率最大值作为优化目标，实现n-1静态安全域的刻画。
[0091]
根据step3在不同情况下的n-1状态下计算系统不平衡度，根据最大系统不平衡度可以确定系统最严重的n-1情况，以最严重n-1状态作为系统安全运行约束，在此基础上，设置新能源渗透率最大值作为优化目标，具体如下所示：
[0092][0093]
其中，为第i个n-1静态安全域边界点所对应的新能源出力增长系数；为第i个n-1静态安全域边界点所对应的新能源功率初始值；为系统总出力，由新能源出力与传统机组出力共同构成。
[0094]
确定新能源聚合后的初始功率与注入方向，如新能源功率沿着某方向增长，根据step4中的目标函数，以及step2中的约束条件，基于最优潮流的方法，求解获得，从而可得到新能源渗透率最大临界值；
[0095]
进一步，更新新能源的注入方向，继续求解不同新能源区域组合下的临界点，直至达到两区域或三区域的坐标边界，连接所有临界点，形成能够表征新能源最大渗透率的n-1静态安全域。
[0096]
step5：将静态安全域映射至新能源出力的低维空间，得到系统的最大消纳能力曲线，实现系统n-1下的满足系统运行安全约束的最大新能源渗透率的计算。
[0097]
具体的，如图2所示，设两区域i、j的新能源初始功率分别为、，将其映射至二维坐标轴的空间内，分别对应着横轴与竖轴的坐标点、，根据坐标点对应的全系统新能源功率注入方向，与坐标点（，）的关系为：
[0098][0099]
式中，为功率增长角，是与横轴的夹角，、分别为区域i、j的新能源出力值。
[0100]
从而将静态安全域映射至新能源出力的二维或三维空间，得到系统的最大消纳能力曲线，新能源消纳能力曲线维度与聚合区域有关，两区域为二维关系，三区域为三维关系，最大消纳能力曲线坐标轴信息均表示不同区域下新能源的出力情况，计算系统的新能源总出力情况，所得新能源总出力情况与系统总功率的比值即为该系统所能允许的新能源最大渗透率，表达为：
[0101][0102]
式中，为系统所能允许的新能源最大渗透率，为新能源总出力，为系统总出力。
[0103]
如图2所示，当系统运行在新能源最大消纳能力的n-1静态安全域边界上点时，通过映射得到对应的二维新能源出力的最大边界值点，进而该点对应的新能源总出力与系统总功率的比值即可得到该系统所能允许的新能源最大渗透率，同理，对应的在安全边界内，则对应的在安全边界外，此时系统新能源消纳能力超过边界值。
[0104]
以上结合附图对本发明的具体实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。

技术特征：
1.一种n-1下新能源渗透率的确定方法，其特征在于：step1：对新能源多区域进行聚合，将新能源分为两区域或三区域；step2：对新能源接入的电力系统发生n-1故障或检修的所有情况进行潮流计算；step3：根据输电线路的负载率，得到系统的不平衡度，以此作为最严重n-1的选择标准并对其进行选择；step4：根据最严重的n-1状态，设置新能源渗透率最大值作为优化目标，实现n-1静态安全域的刻画；step5：将静态安全域映射至新能源出力的低维空间，得到系统的最大消纳能力曲线，实现系统n-1下的满足系统运行安全约束的最大新能源渗透率的计算。2.根据权利要求1所述的n-1下新能源渗透率的确定方法，其特征在于，所述step1具体步骤为：step1.1：基于k-mean聚类将多区域新能源场站进行聚类；step1.2：根据聚类结果，按照容量加权法，将多区域新能源场站进行聚合，等值为两区域或三区域。3.根据权利要求2所述的n-1下新能源渗透率的确定方法，其特征在于，所述step1.1具体步骤为：step1.1.1：在待分类样本中随机选择k个样本数据作为初始凝聚点，每个凝聚点代表1类；step1.1.2：计算其余样本与凝聚点之间的欧式距离，将每个样本归入与它欧式距离最小的凝聚点所在的类；step1.1.3：所有数据归类完成后，将每一类的凝聚点更新为该类中当前数据的均值；step1.1.4：重复step1.1.2，直至前后2次凝聚点不再变化，即实现了数据的分类。4.根据权利要求1所述的n-1下新能源渗透率的确定方法，其特征在于，所述step2中潮流计算的约束条件具体为：；；；；；；；式中，v
gk
、p
gk
、q
gk
分别表示第k台发电机的机端电压、有功功率、无功功率，u
i
、i
i
、p
i
、q
i
分别表示第i个节点的电压幅值、电流幅值、有功功率及无功功率，下标min、max表示对应参量的上限、下限；新能源出力约束，具体如下所示：；
其中，为新能源的最大出力值；n-1断线输电线路传输功率约束，具体表示如下：；其中，表示系统发生n-1后仍运行线路的传输功率，表示对应线路的最大允许传输值，正负号表示传输方向。5.根据权利要求1所述的n-1下新能源渗透率的确定方法，其特征在于，所述step3具体步骤为：step3.1：定义输电线路的负载率为：；其中，r
ij
表示节点i、j所相连输电线路的负载率，p
lij
为对应线路的实际功率传输值，为基态时对应线路所能允许功率传输的最大值；step3.2：采用系统的不平衡度作为表征n-1下严重程度的指标，记为r，根据step3.1求取出不同n-1情况下的各条输电线路负载率，进行求和并计算出输电线路平均负载率；继续计算出各条输电线路与平均负载率的绝对偏差，最终得到不同n-1情况下系统的不平衡度，具体如下所示：；其中，l表示两节点所相连的输电线路集合，n
l
表示输电线路数量；step3.3：根据不同n-1情况下的系统不平衡度，最大系统不平衡度即为最严重的n-1情况，映射关系如下：；其中，c表示与最大系统不平衡度相对应的最严重故障场景。6.根据权利要求1所述的n-1下新能源渗透率的确定方法，其特征在于，所述step4具体步骤为：step4.1：根据step3在不同情况下的n-1状态下计算系统不平衡度，根据最大系统不平衡度可以确定系统最严重的n-1情况，以最严重n-1状态作为系统安全运行约束，在此基础上，设置新能源渗透率最大值作为优化目标，具体如下所示：；其中，为第i个n-1静态安全域边界点所对应的新能源出力增长系数，为第i个n-1静态安全域边界点所对应的新能源功率初始值，为系统总出力，由新能源出力与传统机组出力共同构成；step4.2：确定新能源聚合后的初始功率与注入方向，如新能源功率沿着某方向增
长，根据step4.1中的目标函数，以及step2中的约束条件，基于最优潮流的方法，求解获得，从而可得到新能源渗透率最大临界值；step4.3：更新新能源的注入方向，根据step4.2中的方法继续求解不同新能源区域组合下的临界点，直至达到两区域或三区域的坐标边界，连接所有临界点，形成能够表征新能源最大渗透率的n-1静态安全域。7.根据权利要求1所述的n-1下新能源渗透率的确定方法，其特征在于，所述step5具体为：假设已知两区域i、j的新能源初始功率为和，将其映射至二维坐标轴的空间内，分别对应着横轴与竖轴的坐标点、，根据坐标点对应的全系统新能源功率注入方向，与坐标点（，）的关系为：；式中，为功率增长角，是与横轴的夹角，、分别为区域i、j的新能源出力值；从而将静态安全域映射至新能源出力的二维或三维空间，得到系统的最大消纳能力曲线，新能源消纳能力曲线维度与聚合区域有关，两区域为二维关系，三区域为三维关系，最大消纳能力曲线坐标轴信息均表示不同区域下新能源的出力情况，计算系统的新能源总出力情况，所得新能源总出力情况与系统总功率的比值即为该系统所能允许的新能源最大渗透率，表达为：；式中，为系统所能允许的新能源最大渗透率，为新能源总出力。

技术总结
本发明涉及一种N-1下新能源渗透率的确定方法，属于新能源渗透率评估技术领域。首先对新能源区域进行聚合，聚合后完成新能源接入的电力系统发生N-1故障或检修的所有情况搜寻，根据系统不平衡度最大的标准实现最严重情况的N-1选择，确定最严重N-1情况后，以新能源渗透率最大值作为优化目标，实现N-1静态安全域的刻画。将静态安全域映射至新能源出力的低维空间，得到系统的最大消纳能力曲线，实现新能源最大渗透率的计算。同时，该方法不需要对所有的N-1情况都进行静态安全域的刻画即可求出满足所有安全要求的最大渗透率，且不受系统拓扑的限制。扑的限制。扑的限制。


技术研发人员：束洪春 何业福 王广雪 廖孟黎 董俊 韩一鸣 张焦婕 朱柳青 陈靖 时波涛 李建男
受保护的技术使用者：昆明理工大学
技术研发日：2023.07.19
技术公布日：2023/8/16