一种电力系统热稳定安全域构建的方法与流程

本发明涉及电力系统领域，尤其涉及一种基于bootstrap抽样^[1]的电力系统热稳定安全域快速构建方法。

背景技术：

当前，电力系统n-1热稳定校验完全基于逐点极限计算模式，该方法简单易用，已在电网的运行规划中发挥着巨大作用。但随着电网规模不断发展、随着高渗透率可再生能源大规模并网，该方法已暴露出越来越多的局限性，主要体现在：①在n-1热稳定极限之下，电力系统也可能存在热稳定越限；②极限计算结果受电力系统运行方式变化而难以确定，且某一个固定的极限值难以确保电网的n-1热稳定安全。随着电网规模不断扩大、可再生能源大规模并网，影响电力系统支路潮流分布的因素已经越来越多、呈现了明显的多维特性，仅仅采用过去简单的、一维的极限计算方法来监视、控制电网支路热稳定性已不能满足大电网安全运行、精细化运行要求，亟需发展新技术手段来解决该问题，而静态热稳定安全域(staticthermalsecurityregion,stsr)就是一种最理想、最实用的新的解决方案。

stsr定义为在功率注入空间内，满足电力系统n-1热稳定的所有运行点集合。通过判断电力系统当前运行点与stsr边界的相对位置，可定性评估电力系统的n-1热稳定运行状态，实现全面、直观地监视和控制电力系统的n-1热稳定性，可有效评估电力系统在波动性、随机性和间歇性等诸多不确定性因素影响下的热稳定性，但stsr边界搜索是构建stsr的关键。stsr边界呈现高维非线性拓扑特性，求取复杂度高且难以形成精确的边界解析表达式。通过大量离线计算可获得处于stsr边界上的热稳定临界运行点集，进而形成stsr边界，通过对比电力系统运行点相对stsr边界的位置，可定性评估电力系统电压稳定性，但这类方法难以形成描述stsr边界的解析表达式，不能定量分析电力系统运行点到stsr边界的距离，直观给出电力系统n-1热稳定裕度和制定合理、可行的控制方案以改善系电力统n-1的热稳定性。

大量科学研究和工程实践发现：在工程关心的范围内，电力系统的热稳定安全域边界可用一个或少数几个超平面来近似描述。为获得热稳定安全域边界的近似超平面，目前常采用重复调用潮流方程，搜索大量不同功率增长方向下处于所研究支路n-1热稳定临界的运行点，然后根据这些运行点，采用最小二乘求解对应的近似超平面的系数，进而获取超平面的近似表达式。理论上，搜索的热稳定临界点越多，所构建的近似超平面精度越高，可更高精度逼近支路热稳定安全域边界。但在实际应用过程中，临界点的搜索是热稳定安全域构建最耗时的部分。高精度的近似超平面意味着高耗时的临界点搜索，这显然与热稳定安全域的在线应用相矛盾。

为避免陷入计算精度与在线应用的矛盾中，目前工程应用时，在牺牲一定的计算精度前提下，通过搜索有限个临界点来构建热稳定安全域边界的近似超平面。为保证所构建的热稳定安全域内部不存在空洞，通常还需将所得的超平面进一步平移，以保证热稳定安全域内的运行点都是安全的。显而易见，该方法在保证热稳定安全域内部不存在空洞同时，也将部分稳定运行点排除在热稳定安全域外，带来了较高的保守性。

参考文献

[1]b.efron,r.tibshirani,“bootstrapmethods:anotherlookatthejackknife,”ann.statist.,vol.7,no.1,pp.1–26,1979.

技术实现要素：

本发明提供了一种电力系统热稳定安全域快速构建方法，本发明基于bootstrap抽样，实现准确、快速、高精度地构建电力系统热稳定安全域，详见下文描述：

一种电力系统热稳定安全域构建的方法，所述方法基于bootstrap抽样，所述方法包括以下步骤：

1)根据初始近似超平面解析表达式、结合支路k在n-1下的原始热稳定临界点集合的子集合，获取近似超平面预测值集合、残差集合及残差的经验分布函数；

2)根据预测值集合和残差集合，获取重构值集合；基于重构值集合、原始热稳定临界点集合的子集合，计算n-1下热稳定安全边界近似超平面的解析表达式的系数集合；

3)根据系数集合，计算各系数的经验分布函数的初始均值和初始标准差；

4)重复执行步骤2)-3)若干次，获取基于bootstrap抽样的经验分布函数的bootstrap均值和bootstrap标准差；

5)根据bootstrap均值和bootstrap标准差，结合概率论和统计学中的t分布、χ²分布近似得支路k所对应的经验分布函数的均值和标准差，并计算相应的置信区间，进而求取支路k所对应的近似超平面表达式。

所述根据初始近似超平面解析表达式、结合支路k在n-1下的原始热稳定临界点集合的子集合，获取近似超平面预测值集合、残差集合及残差的经验分布函数具体为：

根据初始近似超平面解析表达式、原始热稳定临界点子集合计算近似超平面预测值集合；

由列向量的实际值集合和近似超平面预测值集合的对应向量之差，获取残差集合；

由求得的残差集合，计算残差的经验分布函数，即均值μ(ξ′)和标准差σ(ξ′)。

所述近似超平面预测值集合为：

其中，

为由第j个原始热稳定临界所得的第j个预测值，β1＝-βk,2/βk,1和βn-1＝-βk,n/βk,1，为第n个描述性参数在第j个原始热稳定临界点处的值，βk,n为第n个描述性参数在第k条支路n-1下的热稳定安全边界近似超平面的系数。

所述计算残差的经验分布函数，即均值μ(ξ′)和标准差σ(ξ′)为：

其中，ξ′为残差集合中的元素。

进一步地，所述根据预测值集合和残差集合，获取重构值集合具体为：

基于残差的经验分布函数，从残差集合中任意抽取一残差ξj′，与的预测值集合中任一元素计算的重构值

重复执行上述步骤z次，获取到的重构值集合。

其中，所述bootstrap均值和bootstrap标准差具体为：

其中，均为bootstrap均值；均为bootstrap标准差；μ′i(α′)、μ′i(β′1)、μ′i(β′n-1)为经验分布函数的均值。

本发明提供的技术方案的有益效果是：

1、本发明在原始热稳定临界点基础上，通过bootstrap抽样生成大量的热稳定临界点，提高了电力系统热稳定安全域边界的近似效率；

2、本发明通过估计超平面系数的经验分布函数，确定各超平面系数的值，并引入超平面系数的置信区间和标准化，以有效评估所得超平面近似电力系统热稳定安全域边界的效果和可靠性；

3、本发明基于统计学的方法求解电力系统热稳定安全域边界的近似超平面，可有效计及了运行方式的不确定性对超平面计算结果的影响；

4、本发明可为电网运行调度人员提供更加快速、准确、高精度的电力系统支路n-1热稳定信息，有利于电网运行调度人员及时掌握电力系统热稳定信息，采用合理措施，改善电网运行的n-1热稳定性。

附图说明

图1为一种基于bootstrap抽样的电力系统热稳定安全域快速构建方法的流程图；

图2为总有功功率的初始预测值图；

图3为初始预测残差图；

图4为某一单线n-1的热稳定安全域及其边界近似超平面图的示意图；

图5为d1的预测值与实际值的对比图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面对本发明实施方式作进一步地详细描述。

为了解决背景技术中关于快速构建电力系统热稳定安全域的不足，本发明实施例提出一种基于bootstrap抽样的电力系统热稳定安全域快速构建新方法，以实现电力系统热稳定安全域边界的快速、准确近似，进而提高电力系统热稳定安全域的构建效率和精度，为电网运行调度人员提供快速、准确、可靠的n-1静态热稳定安全评估，以改善电力系统n-1静态热稳定安全的监视和控制效果。

实施例1

一种基于bootstrap抽样的电力系统热稳定安全域快速构建方法，参见图1，该方法包括以下步骤：

101：根据初始近似超平面解析表达式、结合支路k在n-1下的原始热稳定临界点集合的子集合，获取近似超平面预测值集合、残差集合及残差的经验分布函数；

102：根据预测值集合和残差集合，获取重构值集合；基于重构值集合、原始热稳定临界点集合的子集合，计算n-1下热稳定安全边界近似超平面的解析表达式的系数集合；

103：根据系数集合，计算各系数的经验分布函数的初始均值和初始标准差；

104：重复执行步骤102-103若干次，获取基于bootstrap抽样的经验分布函数的bootstrap均值和bootstrap标准差；

105：根据bootstrap均值和bootstrap标准差，近似得支路k所对应的经验分布函数的均值和标准差，并计算相应的置信区间，进而求取支路k所对应的近似超平面表达式。

其中，上述步骤101中的根据初始近似超平面解析表达式、结合支路k在n-1下的原始热稳定临界点集合的子集合，获取近似超平面预测值集合、残差集合及残差的经验分布函数具体为：

根据初始近似超平面解析表达式、原始热稳定临界点子集合计算近似超平面预测值集合；

由列向量的实际值集合和近似超平面预测值集合的对应向量之差，获取残差集合；

由求得的残差集合，计算残差的经验分布函数，即均值μ(ξ′)和标准差σ(ξ′)。

其中，上述步骤102中的根据预测值集合和残差集合，获取重构值集合具体为：

基于残差的经验分布函数，从残差集合中任意抽取一残差ξj′，与的预测值集合中任一元素计算的重构值

重复执行上述步骤z次，获取到的重构值集合

综上所述，本发明实施例在原始热稳定临界点基础上，通过bootstrap抽样生成大量的热稳定临界点，提高了电力系统热稳定安全域边界的近似效率。

实施例2

下面结合具体的计算公式、实例对实施例1中的方案进行进一步地介绍，详见下文描述：

201：获取电力系统热稳定安全域边界的初始近似超平面解析表达式，包括：

1)将式(1)描述的热稳定安全域边界的近似超平面，改写为式(2)所示的形式：

式中，βk,0、βk,1、…、βk,n为满足第k条支路n-1约束的热稳定安全域边界近似的超平面系数；为电力系统满足第k条支路n-1约束的热稳定安全域边界上的临界运行点，中元素为电力系统的描述性参数，可为支路电流、支路有功功率、支路相位差、区域总负荷、断面有功功率等(具体实现时，本发明实施例对此不做限制)；为第k条支路的热稳定限值；n为安全域的空间维数。

式中，β1＝-βk,2/βk,1和βn-1＝-βk,n/βk,1。

2)采用最优潮流或重复潮流计算得到支路k在n-1下的原始热稳定临界点集合；

其中，原始热稳定临界点集合为列向量，该列向量中元素数量对应于采用最优潮流或重复潮流计算所得的支路k在n-1下的原始热稳定临界点个数，按式(3)计算式(2)所示的支路k所对应的n-1下热稳定安全边界近似超平面的解析表达式的初始系数α⁰,如下：

式中，t为转置。

3)根据式(3)所求解的支路k对应的n-1下热稳定安全边界近似超平面的解析表达式的初始系数α⁰,获取支路k对应的n-1下热稳定安全边界初始近似超平面表达式为：

式中，为的预测值。

202：根据所求得的支路k在n-1下的原始热稳定临界点集合的子集合和n-1下热稳定安全边界初始近似超平面表达式(4)，求由式(4)计算所得的近似超平面预测值集合残差集合ξ′及残差的经验分布函数，包括：

1)根据式(4)和原始热稳定临界点子集合计算近似超平面预测值集合

式中，为由第j个原始热稳定临界所得的第j个预测值。

2)由的实际值集合和预测值集合

求式(4)所示的热稳定安全域边界近似超平面的残差集合

3)由求得的残差集合计算残差的经验分布函数的均值μ(ξ′)和标准差σ(ξ′)：

203：根据残差集合ξ′＝[ξ1′ξ2′…ξm′]和的预测值集合计算的重构值集合包括：

1)基于残差ξ′＝[ξ1′ξ2′…ξj′…ξm′]的经验分布函数，从ξ′＝[ξ1′ξ2′…ξj′…ξm′]中任意抽取一残差ξj′，与的预测值集合

中任一元素按式(8)计算的重构值

2)重复执行步骤1)z次(z≥m×n)，获取到的重构值集合

204：基于的重构值集合及支路k的原始热稳定临界点子集合按式(9)计算式(2)所示的支路k所对应的n-1下热稳定安全边界近似超平面的解析表达式的系数集合α′＝[α′1α′2…α′n]、β′1＝[β′1，1β′1,2…β′1,n]、…及β′n-1＝[β′n-1，1β′n-1,2…β′n-1,n]，包括：

1)基于和由式(9)计算支路k所对应的n-1下热稳定安全边界近似超平面的解析表达式的系数α′1,β′1，1,…,β′n-1,1：

2)重复执行步骤203、及步骤204中的步骤1)n次，获取支路k所对应的n-1下热稳定安全边界近似超平面的解析表达式的系数集合α′＝[α′1α′2…α′n]、β′1＝[β′1，1β′1,2…β′1,n]、…、β′n-1＝[β′n-1，1β′n-1,2…β′n-1,n]。

205：根据步骤204所得的超平面的解析表达式的各系数集合α′＝[α′1α′2…α′n]、β′1＝[β′1，1β′1,2…β′1,n]、…、β′n-1＝[β′n-1，1β′n-1,2…β′n-1,n]，计算近似超平面的解析表达式的各系数的经验分布函数的均值μ′和标准差σ′；

即，以式(10)计算各系数的经验分布函数的均值、和以式(11)计算各系数的经验分布函数的标准差：

206：基于bootstrap抽样，重复执行步骤203、204和205l次，获取基于bootstrap抽样的经验分布函数的均值μ^*和标准差σ^*；；

1)重复步骤203、204和205l次，计算每次的近似超平面的解析表达式的各系数的经验分布函数的均值μ′和标准差σ′，得式(12)所示的各系数α′1,β′1，1,…,β′n-1,1的经验分布函数的均值集合和式(13)所示的各系数α′1,β′1，1,…,β′n-1,1的经验分布函数的标准差集合；

2)针对式(12)和式(13)所得的经验分布函数的均值集合和标准差集合，分别根据式(14)和(15)计算bootstrap抽样的近似超平面各系数的经验分布函数的均值μ^*和标准差σ^*。

其中，均为均值；均为标准差。

207：根据步骤206中计算的均值μ^*和标准差σ^*，近似得支路k所对应的经验分布函数的均值μ和标准差σ，并按式(16)和(17)计算均值μ和标准差σ的η(η通常取值为95％、90％等)的置信区间；。

式中，t为概率论和统计学中的t-分布，χ²为概率论和统计学中的χ²-分布。

208：根据步骤207近似所得的支路k所对应的经验分布函数，得式(18)所示的支路k所对应的近似超平面表达式。

式中，

由式(18)所得在高维空间第一象限内的近似超平面及坐标轴，构建满足支路k出现n-1后仍满足热稳定运行要求的热稳定安全域。

上述具体构建的步骤为本领域技术人员所公知，本发明实施例对此不做赘述。

综上所述，本发明实施例通过上述步骤201-步骤208实现了基于bootstrap抽样的电力系统热稳定安全域边界的近似超平面快速估计，在不降低电力系统热稳定安全域边界近似精度的前提下，实现了基于bootstrap抽样的热稳定安全域边界的快速近似。可为电网运行和控制人员改善电力系统的热稳定安全性提供快速、准确、可靠的热稳定信息，有利于提升电力系统的n-1热稳定安全监视与控制能力。

实施例3

下面结合具体的实例、图2-图5、以及表1，对实施例1和2中的方案进行可行性验证，详见下文描述：

本实例是以在三维有功功率空间中构建位于南方电网广东交流输入断面的某一单线n-1热稳定安全域为例，验证本发明实施例1和2的有效性。

显然，三维空间中式(2)所描述的单线n-1热稳定安全域边界的近似超平面表达式应为通过长期运行经验和研究表明：影响单线n-1热稳定的关键性发电机功率注入节点为yx和ts。

为此，在三维有功功率空间中坐标轴分别选取为流过广东交流输入断面总有功功率，以及发电机功率注入节点yx和ts的有功注入，即式(2)中选择为单线处于n-1热稳定临界时流过广东交流输入断面的总有功功率；和分别选择为单线处于n-1热稳定临界时发电机功率注入节点yx和ts注入到电力系统中的有功功率。基于d2和d3不同的功率增长方向，通过重复潮流计算获取到60个单线n-1的热稳定临界点，进而构成步骤201中的单线n-1热稳定临界点集合其中为列向量，且列向量中元素数量均为60。根据单线n-1热稳定临界点集合获取子集合和

进而根据式(3)求得单线n-1的热稳定安全边界近似超平面的解析表达式的初始系数α⁰、β1⁰和分别为9343.7499、0.3224和0.9135，即式(4)所示的单线n-1的热稳定安全域边界初始近似超平面表达式为

由求得的单线n-1的热稳定安全域边界初始近似超平面表达式及根据式(5)计算的预测值集合其中m＝60，求得的初始预测值如图2所示。

根据重复潮流计算所得的原始实际值集合为和预测值集合

由式(6)计算该热稳定安全域边界近似超平面的残差集合ξ′＝[ξ1′ξ2′…ξm′]，结果如图3所示。

由图3中的残差集合ξ′＝[ξ1′ξ2′…ξm′]，根据式(7)计算残差ξ′的经验分布函数的均值和标准差。进而根据步骤203，由残差集合ξ′＝[ξ1′ξ2′…ξm′]和的预测值集合

按式(8)计算的重构值，重复该过程3000次，得的重构值集合

根据步骤204，由的重构值集合及按式(9)计算式(2)所示的单线n-1的热稳定安全边界近似超平面的解析表达式的系数集合α′＝[α′1α′2…α′n]、β′1＝[β′1，1β′1,2…β′1,n]和β′2＝[β′2，1β′2,2…β′2,n]，此处n取值为1000。由求得的系数集合α′、β′1和β′2，进一步由步骤205的式(10)和(11)计算系数集合α′、β′1和β′2的经验分布函数的均值和标准方差，其中α′的经验分布函数的均值和标准方差分别为9343.6801和3.1703；β′1的经验分布函数的均值和标准方差分别为0.3224和0.0047；β′2经验分布函数的均值和标准方差分别为0.9135和0.0069。

表1某一单线n-1的热稳定安全域边界近似超平面系数计算结果

根据步骤206重复执行步骤203、204和205共l次(算例中l取值为200，根据实际应用设定)，获取到基于bootstrap抽样的单线n-1的均值和标准差。进一步再根据步骤207，获取到单线n-1的经验分布函数的均值和标准差，及均值和标准差获得95％置信区间，如表1所示。显然，由表1中结果可知：采用本方法计算的单线n-1的系数α、β1和β2分别为9343.7427、0.3224和0.9135。

进而由步骤208获取单线n-1的解析表达式为d1＝9343.7426+3.1171d2+0.9135d3。由该解析表达式进而可获取到单线n-1的热稳定安全域。

图4进一步给出了通过本方法求得的近似ym单线n-1的热稳定安全边界的超平面，显然，由图4中结果可见：图4中的ym单线n-1的热稳定临界点均位于本方法所得的超平面上，表明本方法所得的超平面具有较高的安全域边界近似精度。

图5进一步对比了在ym单线n-1的热稳定临界点处，根据本方法所得的近似超平面解析表达式预测的d1值和d1的实际值，由图5中结果可知：根据本方法所预测的d1值与d1的实际值近似相等。进一步，根据图5中的d1的预测值和实际值，计算本方法所得单线n-1的热稳定安全域边界近似超平面的近似误差，综上由图5中的结果可知：

本方法所得的近似超平面最大相对近似误差为0.3211％，平均相对近似误差为0.0805％，结果表明本方法所得的近似超平面对单线n-1的热稳定安全域边界具有较高的近似精度。

本领域技术人员可以理解附图只是一个优选实施例的示意图，上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。