基于双层代理模型的概率-区间潮流计算及灵敏度分析

王晨旭 唐 飞 刘涤尘 高 鑫 周依希

基于双层代理模型的概率-区间潮流计算及灵敏度分析

王晨旭1唐 飞1刘涤尘1高 鑫1周依希2

（1. 武汉大学电气与自动化学院 武汉 430072 2. 国网浙江省电力有限公司杭州供电公司 杭州 310000）

随着电力系统中不确定量日益复杂，同时存在的随机与区间变量使得采用概率或区间潮流计算难以准确获取系统的运行状态。为此，提出一种基于双层代理模型的概率-区间潮流计算方法。该方法仅需较少次数的确定性潮流计算便可实现上、下层代理模型的构建，进而通过代理模型求解概率-区间潮流常规求解方法中所需要的大量确定性潮流计算，可实现输出变量的快速获取。此外，该文还提出了用于描述输出变量特征的灵敏度指标，并结合所提出的双层代理模型开展灵敏度分析，以量化输入区间变量对输出变量的影响程度。在IEEE 118节点系统中进行算例分析，通过与已有方法对比验证了所提方法的精确性和快速性，借助灵敏度分析可识别对输出变量具有显著影响的关键区间变量，有助于揭示系统运行状态与区间变量之间的关系。

不确定潮流 随机变量 区间变量 代理模型 灵敏度分析

0 引言

电力系统运行中存在诸多不确定性因素[1]，例如负荷需求、线路故障等。随着以风电和光伏为代表的可再生能源大规模并网，系统运行的不确定性进一步增加[2]。因此，有必要采用适当的模型来表示各类不确定性因素，并利用不确定潮流计算对系统运行特性进行分析，从而更加准确地掌握系统的运行状态。

根据不确定量建模方式的不同，常用不确定潮流计算方法包括两种：概率潮流[3]和区间潮流[4]。概率潮流计算将输入不确定量描述为随机变量，进而采用模拟法[5]、点估计法[6]、解析法[7-8]和多项式混沌展开法[9]等获取输出变量的概率分布。区间潮流计算将输入不确定量描述为区间变量，通过仿射法[10-11]和区间泰勒展开法[12]等方法获得输出变量的上、下边界。然而，随着系统的复杂性增加，随机与区间变量往往同时存在。若采用概率潮流分析，则需假设区间变量的概率分布，导致所得结果不够准确；若采用区间潮流分析，则需丢弃随机变量的统计信息，造成对系统状态的保守估计。因此，需要研究能够同时考虑随机和区间变量的不确定潮流计算方法，即概率-区间潮流。

当随机与区间变量同时存在时，系统输出变量的波动范围为其最大、最小概率分布所围成的区域，即概率盒[13]。目前仅有少量关于概率-区间潮流计算方法的研究。文献[14]提出基于证据理论概率-区间潮流模型，并将其转换为多个区间潮流进行求解。文献[15]提出一种双层抽样法，该方法利用外层和内层抽样分别处理随机和区间变量，从而通过大量确定性潮流计算实现概率-区间潮流的求解。文献[16]将配电网中分布式电源出力和负荷需求分别视作随机和区间变量，并提出了一种基于仿射线性三相潮流的近似方法求取输出变量。文献[17]提出一种基于聚类的解析算法，该方法通过在统一最优场景下进行概率潮流计算，可快速且精确地求取单个输出变量的概率盒。针对电-气联合系统中随机与区间变量共存，文献[18]提出一种基于多项式混沌展开的概率-区间能量流计算方法。

虽然当前概率-区间潮流研究取得了一定的成果，但仍存在以下不足：一方面，算法的计算效率仍有待提升。例如双层抽样法需要进行超过106次确定性潮流计算才能得到准确的结果；基于聚类的解析算法虽然在求取单个输出变量时效率较高，但运行人员往往需要了解多个节点电压或多条支路潮流的状态。另一方面，现有研究侧重于如何获取输出变量的概率盒。然而，当系统中包含多个区间变量时会引起输出变量的波动范围较大，导致运行人员难以对系统运行状态作出准确判断。因此，有必要开展灵敏度分析来量化系统中各区间变量的影响程度，从而识别出影响系统状态的关键因素[19-20]。

针对上述不足，本文提出一种基于双层代理模型的概率-区间潮流计算方法。在所提方法中，上层代理采用径向基函数（Radial Basis Function, RBF）模型来近似系统非线性潮流方程，从而能够快速获得输出变量的边界值；基于上层代理所得结果，下层代理采用稀疏多项式混沌展开模型（sparse Polynomial Chaos Expansion, sPCE）构建输入随机变量与输出变量边界值的关系。本文所提出的双层代理模型仅需要较少次数的确定性潮流计算便可构建，进而快速求取输出变量的概率盒。此外，本文还提出了用于描述概率盒特征的灵敏度指标，并结合双层代理模型开展灵敏度分析来量化输入区间变量对输出变量概率盒的影响程度。通过在IEEE 118节点系统中进行算例分析，验证了所提方法的精确性和快速性。

4.2 产品研发投入不足 产品研发缺乏持续的资金支持，对不同土壤条件犁体适应性、关键部件耐用性、进一步提高整机作业性能等方面没有持续深入的研究；农机企业已形成“行业兴起—模仿抄袭—高速增长—产能过剩”恶性发展模式，导致企业更不愿意进行研发投入，产品更新换代慢，跟不上市场发展需求。

1 不确定变量模型

在含大规模风电的电力系统中，风电出力与负荷的不确定性对系统的优化运行和稳定性产生显著影响[2]。对于负荷，其波动范围较小，预测误差通常可以采用正态分布来表示。对于风电出力，由于受到风速、风向等天气因素的影响，其波动范围较负荷更大，且预测误差的概率分布难以获得。基于上述原因，本文将负荷需求和风电出力分别表示为随机变量和区间变量。

通常有功负荷需求的波动性可采用正态分布描述，即

图1 区间变量的相关性

对于风电出力Y和Y，若忽略二者相关性，则两个风电场出力的联合采样区域为黑色实线所围成矩形AECF。当二者相关时，所构成的区域为点画线所围成的平行四边形ABCD，且相关性系数可表示为

平行四边形ABCD所围成的联合采样区域[22]可表示为

式中，直线的斜率1、2和常数1～4由区间变量Y和Y的宽度及其相关性系数决定。

2 概率-区间潮流模型

假设系统中输入变量包含维随机变量和维区间变量。在潮流计算中，输出变量与输入变量之间的非线性关系可以简写为

式中，R=[1XX]T为维随机变量；I=[1YY]T为维区间变量；为节点电压和支路功率，也称为输出响应；(·)为由潮流方程决定的输出、输入变量之间的函数关系。

由于式（4）中同时存在随机变量与区间变量，难以直接通过概率潮流或区间潮流方法进行求解。当前，双层抽样法[15]是求解概率-区间潮流最直接的方法。该方法分别在外层和内层对随机变量和区间变量进行抽样，并将样本代入到非线性潮流方程中进行求解。双层抽样法的计算步骤总结如下：

不是吗，无限无知的宇宙，似乎天然就内在具有一种毫不犹豫的“生命指向”，在一切可能的极度艰辛中一旦有缝隙，就会“石上开花”、生命问世。没有生命的宇宙无法证明其自身的“在”与“不在”，因此，从植物到微生物到昆虫到动物等等，生命以它层层递进的宏大与渺小，让这不被思索的无限广宇在知与不知的替换中，得到思索追溯。

式中，输出响应采用标量表示，代表系统中任一节点电压或支路功率。

（2）在内层对维区间变量I进行抽样，得到样本矩阵I，其中为内层抽样次数；将随机变量的样本R与区间变量样本I结合并代入（5）中可得

在实际应用中，对于式（14）中的PCE模型，通常需要设置其最高展开阶数，将其表示为有限项截断逼近模型，即

皮肤的微细结构就像砖墙，皮脂和结缔组织就是灰泥，皮脂膜呢，就是石灰和乳胶漆。皮脂膜不足，就好像砖墙外露了。要保证这个“砖墙”外观好看且功能良好，就要避免皮脂膜的减少。

（6）对于系统中每个输入区间变量，重复步骤（3）～步骤（5），并结合式（22）得到各区间变量对输入变量概率盒的贡献度。

结晶岩孔壁安全度高，但钻遇断裂带、裂隙带、破碎带或低强度带等时，如有漏失，则漏失压力低(密度上限一般低于1 g/cm3)，调节泥浆密度难以满足压力，一般需采用空气钻井或堵漏措施；同时其坍塌压力高(密度下限一般高于2.5 g/cm3)，调节泥浆密度难以满足应力平衡要求，一般会通过坍塌扩径释放局部压力至应力重新平衡，为了维护孔壁安全，宜采用固壁或造壁技术护壁。

图2 输出变量Z的概率盒

3 基于双层代理模型的概率-区间潮流

近年来，代理模型被广泛应用于多类电力系统不确定性量化分析中。常用的代理模型包括RBF模型和sPCE模型。本文以上层代理模型近似非线性潮流方程，由于需要同时考虑输入随机与区间变量，可采用RBF模型来构建；以下层代理模型近似输入随机变量与输出变量边界值的关系，此时仅需要考虑输入随机变量，可采用sPCE模型构建。

3.1 上层代理模型

RBF代理模型也称为RBF神经网络，是一种基于多变量插值原理的前馈型神经网络，具有以任意精度逼近任意连续函数的能力，可应用于负荷与新能源预测[23-24]和潮流计算[25]。考虑输入变量为维随机变量和维区间变量，输出变量为任一节点电压或支路功率，则RBF代理模型的结构如图3所示。输入层包括+个神经元，对应于输入变量的维数；隐含层包含个神经元，采用RBF作为激励函数；输出层对应于系统节点电压或支路功率。

图3 RBF代理模型

隐含层中RBF一般选择高斯函数形式，即

6、特别注意抓好幼苗期虫害防治（播种的当天，5厘地的苗床用50g敌百虫加少量水溶解后，拌玉米粉2-2.5公斤撒于苗床上及四周防蟋蟀，施药后遇雨要重新施药，或用旧塑膜在苗床四周建1米高的围栏防虫，效果较用药更好）。

基于图3中RBF代理模型结构，可将输出响应表示为

式中，w为输出层权值，表示第个隐含层神经元对输出响应的贡献。

采用RBF代理模型来近似非线性潮流方程时，其精度取决于隐含层神经元的中心和扩展常数。为了提高构建代理模型的速度，本文以精确拟合法[25]来构建RBF模型。该方法将隐含层神经元的个数设置为输入向量的个数，从而使RBF模型能够以零误差拟合输入与输出之间的关系。

式中，()为由第个输入向量通过非线性潮流计算得到的输出响应。进一步将式（9）写为矩阵形式，即

(5)中国农村地区收入不平等显现了一定的空间聚集性，收入不平等更多地聚集在城镇化发展较快的城市集群区域，但城市群内部的要素分布并不均衡。从1986—2014年，城市群与非城市群的内部收入不平等解释了总的地区收入不平等的97%以上，且非城市群与城市群农村地区收入不平等从36.7∶62.0上升到2014年24.5∶72.9，地区内部收入不平等进一步集中在城市群地区，侧面说明当前形成的城市群并没有有效缓解区域内部的收入不平等。因此，制定政策时要避免劳动力、资本等要素过度集中到城市群核心区域，应引导资本向城市群边缘村庄流动，这是缩小农村居民收入不平等的关键。

根据式（7）～式（10）可知，当RBF模型中神经元的个数确定后，该模型输出层权值受到扩展常数的影响，因此，其拟合精度也随着选取的扩展常数不同而有所差别。扩展常数越大，RBF模型逼近曲线越光滑，但过大的扩展常数会导致模型欠缺选择性；扩展常数过小时，模型的逼近效果较差。因此，本文采用交叉验证误差（Cross-Validation Error, CVE）[26]来确定合适的扩展常数。该方法每次使用输入向量中-1个向量来构建RBF模型，并将剩余的1个输入向量作为验证计算所得模型的误差，通过次计算得到误差的期望值。CVE的定义为

在确定隐含层神经元的中心、扩展常数与输出层权值后便可得到如式（8）所示的RBF模型，从而快速地估计输出变量的边界值。采用双层抽样生成输入向量，外层对随机变量R进行抽样得到输入向量中前个元素，内层对区间变量I进行抽样得到剩余的个元素，则输入向量可表示为

式中，RBF和RBF分别为外层、内层的抽样次数。输出变量的边界值可以表示为

RS、GPS和GIS等现代技术的应用有助于提升土地利用现状的调查质量与精度，借助这些高新技术与平台所建立的土地利用变更分析与管理系统也可大大满足我国土地资源日常管理的需要。本文借助3S技术，设计并实现了基于SuperMap Object的土地利用变更分析与管理系统，该系统的应用实例表明，其能够有效满足土地利用变更调查工作需求，值得广泛应用推广。

采用RBF模型进行RBF×RBF次计算得到输出边界值URBF=[U1…U…URBF]T和L RBF=[L1…L…LRBF]T。相较于确定性潮流计算，RBF模型仅需将输入变量样本代入其中即可得到输出响应，可极大地提高计算效率。需要指出的是，虽然将RBF代理模型与双层抽样结合能够得到输出变量边界值的样本，倘若所需计算的样本数量过多将会使总计算时间增加。因此，需要构建下层代理模型以进一步提升计算效率。

3.2 下层代理模型

多项式混沌展开（Polynomial Chaos Expansion, PCE）是将输入随机变量和输出响应之间关系以一组正交多项式展开和的形式来近似表达。该模型被广泛应用于概率潮流计算[9,27]和全局灵敏度分析[19-20]。在构建PCE模型时，要求输入随机变量之间相互独立，因此首先需要通过Nataf变换[6]将具有相关性的输入随机变量R=[1…X…X]T转换为独立的标准正态分布随机变量，即R=[1…ξ…ξ]T。同样地，要将输入随机变量样本矩阵RRBF转换为样本矩阵RRBF。在此情况下，可采用各输入变量对应的正交多项式混沌展开来逼近输出响应，即

式中，OC和OD分别为线段OC和OD的长度。

多变量多项式基函数Ψ()可以通过单变量多项式的张量积求取，即

式中，φ(ξ)为变量ξ的第i阶正交多项式。通常，各输入变量对应的正交多项式取决于该变量的分布类型，本文中服从正态分布的输入变量对应的最优基函数为Hermite多项式。其他常见的连续型概率分布及其对应正交多项式体系可见文献[19]。

基于边界值U和L构建输出变量的概率盒如图2所示。输出变量最大概率分布max()和最小概率分布min()分别代表其下边界和上边界的概率分布。由图可知，输出变量的波动范围是由最大、最小概率分布所围成的区域。在概率-区间潮流中，输出变量越限事件＜limit发生的概率为区间[Prmin, Prmax]。因此，相较于概率潮流计算，概率-区间潮流的优势在于当系统中某些节点注入功率的概率分布难以准确获取时，可将其视作区间变量，从而能够估计出输出变量超过其限值的最大、最小概率。

截断后逼近模型的展开项数为

一般地，最高展开阶数取值越大，多项式展开逼近精度越高，但展开项数也会相应增加。在工程问题中，通常取2～3即可获得较高的计算精度，在本文中取=2。将式（16）采用向量形式表示为

为求解PCE模型中待求展开系数，需要利用上层RBF模型所得结果，即随机变量样本RRBF及其所对应的输出变量最大值样本U RBF，进而采用最小二乘法求解待定系数为

式中，为独立随机变量R的第个样本。为保证求解精度，通常取样本数量为2～3。由样本系数的计算表达式可知，当输入随机变量个数增大时，待定系数个数也随之快速增大，导致上层代理模型的计算量较大。由于PCE模型具有稀疏性，为了提高构建模型的效率，文献[27]提出在具有少量输入样本的情况下采用压缩感知技术实现展开系数的快速计算，从而构建sPCE模型。该方法基本思想是忽略对输出响应影响较小的展开项，仅考虑对于输出响应具有较大影响的展开项。基于展开系数向量的稀疏性，可以通过如下的优化问题进行展开系数的求解。

式中，||||0为待求系数向量中非零元素的数量。通过求解式（20），使得在给定较少的样本下能够得到一组含非零项最少的展开系数向量，从而构建sPCE模型。式（20）中优化问题可采用正交匹配追踪算法来求解，文献[20, 27]给出了该算法的详细步骤，本文不再赘述。

问卷调查结果显示，159例用药患者中吸烟、饮酒者分别为60例(37.7%)和52例(32.7%)；女性为97例(61.0%)，男性为62例(39.0%)。

3.3 灵敏度分析

在概率-区间潮流中输出变量以概率盒的形式表示，其波动范围处于最大、最小值概率分布之间。然而，当最大、最小值概率分布所涵盖的波动范围较大，则会导致其难以为运行人员提供关于系统运行状态有价值的信息。因此，需要通过开展灵敏度分析辨识出对特定输出变量影响较大的输入区间变量，以便采取措施减小其波动范围。定义描述输出变量概率盒的波动范围的平均距离指标为

这也意味着，教材尤其是高校教材绝不是抄抄写写的简单工作，而是一件复杂而有创造性的工作。甚至，编教材有时候比科研还难，因为写论文只需要明确表达自己的观点就够了，编写教材需要综合考虑多方面的因素，如读者的情趣与能力，教学的要求与师生的要求以及与相关学科的关系等。

式中，U和L为输出变量的第组上、下边界值。

本节采用IEEE 118节点系统验证所提方法的有效性。将系统中各节点有功负荷视作随机变量，且服从正态分布，期望值为原系统负荷值[28]，标准差为期望值的5%。将系统分为两个区域，区域1包含节点1～60，区域2包含节点61～118。同一区域内负荷之间的相关性系数为0.6，不同区域内负荷之间相关性系数为0.4。系统中接入八个风电场，将其分为两组，各组风电场接入节点及参数设置见表1。对于风电场的相关性，考虑以下两个场景：①场景1中不同风电出力之间相互独立；②场景2中同组风电出力之间的相关性系数为0.5，不同组风电出力之间相互独立。

对于输入区间变量相互独立的情形，区间变量Y取其中间值时，其余区间变量的波动范围未受影响，此时，ADIC仅包含变量Y对输出变量概率盒的贡献度。对于输入区间变量具有相关性的情形，根据图1所示平行四边形模型，区间变量Y取其中间值时，与其相关的区间变量波动范围需要满足式（3），此时，ADIC包含变量Y及其相关性对输出变量概率盒的贡献度。

需要指出的是，本研究虽侧重于分析输入区间变量对输出变量概率盒的影响，但所提出的双层代理模型也可用于输入随机变量的灵敏度分析。这是因为在下层代理中利用sPCE模型建立了输入随机变量与输出变量边界值之间的关系。sPCE模型可结合用于输入随机变量灵敏度计算的Sobol分解[19]，准确且快速地得到各输入随机变量对输出响应的贡献度。利用sPCE模型开展输入随机变量灵敏度分析的具体步骤可参考文献[19-20]。

4 算法流程

本文所提出的基于双层代理模型的概率-区间潮流计算及其灵敏度分析的计算流程如下：

（1）数据输入。包括电力系统潮流计算所需参数，输入随机变量R的分布参数及相关性系数，输入区间变量I的边界值参数及相关性系数。

以2016年胶东半岛耕地总面积115.95万hm2为准，依据面积比例进行平差，最终得出各耕地地力等级面积，见表3。其中，一、二级地为高产田，占耕地总面积的38.21%；三、四级地为中产田，占耕地总面积的34.61%；五、六级地为低产田，占耕地总面积的27.18%。

（2）构建RBF模型。根据输入变量R和I生成个样本向量=[1……]T，并将各向量代入非线性潮流方程(·)中得到相应的输出响应()=[(1)…()…()]T，进而通过式（9）～式（11）求解RBF代理模型s(·)中输出层的权值。

（3）基于RBF模型计算输出变量的边界值。采用拉丁超立方抽样生成对应随机变量R的维独立标准正态分布变量R的样本矩阵RRBF，并将其经逆Nataf变换为R的样本矩阵。将RRBF与区间变量样本矩阵IRBF通过式（12）进行组合，并代入式（13）中获取输出变量的边界值向量URBF和LRBF。

通过对不同囤养阶段下三种规格的雌性中华绒螯蟹体肉脂含量和脂肪酸的分析，发现中华绒螯蟹在囤养阶段体肉脂含量和脂肪酸之间均存在显著差异。这可能是由于囤养前选择的蟹本身存在生长阶段的差异，以及饲料品质、环境因素和蟹活动情况等因素密切关联。三种规格的雌性中华绒螯蟹体肉在囤养第 4阶段检测出的脂肪酸种类最低，此时已不适合继续囤养。以ω3/ω6和EPA+DHA值为评价指标，囤养第二阶段对人体健康更有利，且一级蟹囤养效果更佳；以油酸含量为评价指标，囤养第2阶段体肉营养品质较好。

（4）构建sPCE模型。基于维独立标准正态分布变量R的样本矩阵RRBF以及输出变量的边界值向量U RBF和L RBF，采用正交匹配追踪算法计算展开系数，得到sPCE模型U s(·)和L s(·)。

5 算例分析

5.1 算例设置

式中，ADI为区间变量Y的取中间值时输出变量的平均距离。

表1 风电场参数设置

Tab.1 Parameters of wind farms

为了验证本文所提方法的有效性，对如下三种算法进行对比测试：①双层抽样法（Double Layer Monte Carlo Simulation, DLMCS）[15]。以该方法所得结果作为对比依据，将其外层、内层抽样次数均设置为2 000次，因此需要4×106次非线性潮流计算。②基于聚类的解析法（Clustering-based Analytical Method, CAM）[17]。将CAM中聚类个数设置为4个，该方法能够快速且精确地获取单个输出变量的概率盒。③本文方法，即双层代理模型法（Double Layer Surrogate Method, DLSM）。算例分析利用Matlab R2020a平台进行编程实现，程序运行的硬件环境为：R7-4750U，基准频率1.7GHz，内存为16GB。

“大队胡人走了，留下了五个，他们在路边堆了个雪人，雪人有一人多高。我暗暗纳罕，难道胡人也是小孩心性，要堆雪人打雪仗？

在本文所提方法中通过上层RBF模型s(·)来获取输出变量的边界值样本，进而将其用于下层代理模型的构建。这意味着由RBF模型所得边界值样本将直接决定所提方法的精确性。定义边界值平均误差η对边界值样本的精度进行定量分析。

式中，上标“”可为上、下边界值；s,和DLMCS分别为通过代理模型和DLMCS所得输出变量的边界值样本；为所得边界值样本的数量。

跨境电商是通过网络完成产品的销售，产品从售前、售中、售后都是在网上通过客户服务人员完成的，面对全球的客户，这就要求对客户的服务意识很强。由于我国在外贸方面的人才短缺，外贸企业在电商中选拔出优秀的客户服务人员很难，只能选择懂英语、了解电商运营的人员完成客户服务工作，但是这些人员知识单纯的对产品咨询进行相应的服务，在对客户的态度或是服务上还存在一定的欠缺，缺少销售人员应当具备的素质，服务意识缺乏。

对于所得输出变量的概率盒，定义相对误差e指标来定量分析其准确性，e为

式中，为输出变量的统计矩，一般为期望值和标准差；M和DLMCS分别为通过待评估方法和DLMCS计算所得结果。由于同一类型的输出往往包含多个变量，采用同一类变量相对误差的平均值来表示结果的精确性。因此，分别采用U,mean和L,mean表示输出变量最大、最小值的相对误差。

5.2 上层代理模型精度验证

本节首先验证采用RBF模型求出输出变量边界值的精确性。所用测试系统包含99个随机变量（有功负荷），8个区间变量（风电出力）。因此，在构建RBF模型中输入向量的维度为107。将RBF代理模型中神经元的数量设置为500，即产生500个随机变量与区间变量样本作为径向基函数的中心，进而采用式（11）中CVE指标来选取合适的扩展常数。

以支路有功功率1-2为例，图4给出了不同扩展常数下该变量的CVE指标。可以看出，当扩展常数小于0.5×104时，RBF模型拟合精度较差；而当扩展常数超过6×104后，RBF模型的拟合精度具有波动性，且呈现误差逐渐增大的趋势。当扩展常数为2.5×104时，所对应的CVE指标较低，因此本文选择该值作为用作RBF模型构建的扩展常数。

图4 不同扩展常数下的CVE

利用所构建的RBF模型可获取输出变量边界值样本，为验证其在精确性上的优势，将其与3.2节所介绍的sPCE模型进行对比。选取sPCE模型作为对比方法的原因是其已被广泛地应用于多类电力系统不确定性量化问题[20, 27]，能够以较高精度近似表示非线性潮流方程中输入和输出变量之间的关系。在构建sPCE模型时需将系统中所包含的区间变量视作均匀分布的随机变量，并以潮流输出变量作为输出响应。将sPCE模型中最高展开阶数设置为2，输入随机变量的数量为107个，则待求解的系数为5 886项。与构建RBF代理模型类似，采用500个样本来求解sPCE模型的展开项系数。

利用所得RBF模型和sPCE模型来获取输出变量的边界值，并以DLMCS所得结果为标准计算两种代理模型的平均误差。图5给出了场景1中通过两种代理模型所得支路有功功率上、下边界值的平均误差（测试系统中包含186条支路）。对比两种代理模型所得结果的误差，可以发现RBF模型在计算精度上要高于sPCE模型，尤其是对于编号90～130支路，两者之间的精度具有明显差异。两种模型最大误差均发生在编号107支路，由RBF模型所得上、下边界值最大误差分别为0.126MW和0.176MW，而sPCE模型所得最大误差分别为1.074MW和1.305MW。此外，RBF模型所对应的边界值误差还呈现两个规律：一方面上、下边界值的误差相近，说明该模型对于求取输出变量上、下边界具有相似的精度；另一方面，不同输出变量的拟合精度有差异，例如编号30～55支路的误差较小，而编号90～130支路的误差较大，这是由于后者受输入区间变量的影响而波动范围较大，使得RBF模型的精度略有下降。

图5 不同代理模型所得结果的边界值误差

对于sPCE模型，由于其待求展开项系数较多，通过少量样本所构建的代理模型能够较为精确地拟合输出变量的整体分布，但难以保证所得输出变量边界值的准确性。对于RBF模型，其待求的输出层权重与样本数量一致，通过扩展常数的优化选取可使其较为精确地求取实输出变量边界值。因此，在本文中将其作为上层代理模型是合理的。

5.3 输出变量概率盒精度验证

在5.2节中已验证了RBF模型获取输出变量边界值的精确性。基于边界值样本可构建下层代理模型，进而得到双层代理模型，即DLSM。对于下层所采用的sPCE模型，设置其最高展开阶数为2，由于系统中存在99个随机变量，则待求解的系数包含5 050项。基于多项式混沌展开的稀疏特性，采用RBF模型产生500个边界值样本用于sPCE模型展开系数的求解。

为验证DLSM所得输出变量概率盒的精确性，采用CAM作为对比方法，以DLMCS所得结果作为标准，得到DLSM和CAM所得结果的相对误差见表2。由结果可知，在场景1和2中DLSM均能够取得与CAM相接近的计算精度。通过DLSM所得各类输出变量的期望值和标准差的相对误差均小于1%和2%，说明了所提方法能够较为精确地获取输出变量的概率盒。此外，DLSM的相对误差结果还呈现出三个规律：①同一场景中输出变量的上、下边界值的相对误差相近，表明了所提方法在求取输出变量上、下边界时具有相似的精度；②同一类型输出变量边界值的相对误差在场景1和场景2中非常接近，说明所提方法的计算精度不受是否考虑区间变量相关性的影响；③支路无功功率期望值的相对误差大于其余三种类型输出变量，这是因为系统中PV节点无功出力使得代理模拟在近似输入变量与支路无功功率之间的关系时精度有所降低。

图6～图8给出了在场景1和2中通过DLMCS和DLSM所得变量33、47-49和68-116的概率盒，选取上述变量进行分析的依据是：一方面这些变量与系统中输入区间变量所在节点邻近，受到区间变量的影响较大，可验证所提方法求取波动较大变量的计算精度；另一方面，变量33和68-116不仅受到输入变量的影响，还受到邻近PV节点无功出力的影响，使得其与输入变量之间的关系更为复杂，可验证所提方法在求取该类变量时的稳健性。

表2 输出变量边界值的相对误差

Tab.2 Relative errors of output variables’ boundaries

图6 节点33电压幅值的概率盒

图7 支路47-49有功功率的概率盒

图8 支路68-116无功功率的概率盒

同一变量的概率盒包含最大、最小两条概率曲线。以DLMCS所得曲线（图中实线）作为准确结果，可以发现在场景1和场景2中由DLSM所得曲线（图中虚线）均与准确曲线非常接近。该结果表明了所提方法能够以较高的精度获取输出变量概率盒。由局部放大图可知，所得变量33和68-116分布曲线的精度要低于变量47-49，其原因在于节点电压与支路无功功率不仅与输入变量有关，还与PV节点的无功出力有关，使得代理模型在曲线局部的拟合精度略有下降。此外，同一输出变量在场景2中的概率盒被包含在场景1所得结果中，说明了输出变量最大、最小概率分布曲线之间的距离与区间变量有关，且当区间变量具有相关性时会导致输出变量概率盒变窄。

概率-区间潮流所得结果可用于分析系统中节点电压或支路功率的越限概率。表3给出了不同输出变量越限概率。由结果可知，在场景1和场景2中，通过DLSM所得变量越限概率与DLMCS和CAM的结果相近。当考虑区间变量相关性时，由于输出变量的概率盒变窄，其越限概率区间也相应地变窄。例如，变量47-49小于-20MW的概率区间宽度从场景1中0.798减少至场景2中0.602。上述结果说明了在实际应用中考虑区间变量相关性的必要性，即避免对系统输出变量的保守估计，从而有利于运行人员做出更具针对性的决策。

表3 变量越限概率区间

Tab.3 Probability intervals of constraint violations

5.4 灵敏度分析结果

采用本文所提DLSM，并结合式（21）和式（22）所示指标对区间变量进行灵敏度分析。图9给出了场景1中，即区间变量相互独立时，系统中8个区间变量（以表示，下标代表该变量所在节点）对不同输出变量概率盒的影响程度。由结果可知，对于不同的输出变量，影响其概率盒的最关键因素是不同的。例如，对于变量23-24而言，对其具有最大影响的变量为27；而对于变量82-96具有最大影响的变量为100。引起该现象的主要原因是不同节点之间的电气距离存在差异。在特定节点引入的区间变量往往对与其电气联系较为紧密的节点和支路影响较大。根据灵敏度分析结果可以将区间变量按照其贡献度大小进行排序。以输出变量75-118为例，可得到排序为：100＞90＞80＞116＞27＞45＞15＞35，其中排序前4位的区间变量的贡献度占据主导。

图9 场景1中输入区间变量的平均距离指标贡献度

为了验证本文方法所得区间变量贡献度排序的正确性，将其与另外两种方法所得结果进行对比。对比方法一采用文献[29]提出的Bhattacharyya距离贡献度（以下简称BDIC）。该方法利用输出变量的最大、最小概率密度曲线的重叠面积来表征各区间变量的贡献度。某一输入区间变量所对应的BDIC越大，表明其对输出变量的影响程度越大。BDIC的定义及其具体计算步骤可见文献[29]。对比方法二采用输出变量越限概率区间宽度作为灵敏度指标。在该方法中，当某一区间变量被设置为其中间值时，输出变量越限概率区间的宽度会减少，因此，可利用其反映该区间变量的影响大小。去除某输入区间变量后输出变量的越限概率区间宽度越大，表明该区间变量对输出变量的影响程度越小。

不同方法所得区间变量贡献度结果与排序见表4。其中，对于越限概率宽度指标来说，当各输入区间变量都存在时变量75-118＜35MW的概率为[0.229, 0.968]，该区间的宽度为0.739。由结果可以看出，基于本文所提出的平均距离指标、BDIC和越限概率区间宽度所得到的区间变量排序是一致的，从而验证了本文所提出的灵敏度指标的合理性。需要指出的是，本文所提出的平均距离指标可直接通过输出变量样本计算得到，而BDIC的计算则需要利用积分，因此，前者的计算更加方便快速。

表4 不同方法的灵敏度分析结果

Tab.4 Sensitivity analysis results of different methods

图10给出了场景2下，即考虑区间变量相关性时，系统中区间变量对不同输出变量概率盒的影响程度。相较于图8中的结果，考虑相关性会使得各区间变量对输出变量概率盒的贡献度增加。这是因为当计及区间变量相关性时，区间变量Y取其中间值时，其余区间变量的波动范围需要满足图1所示平行四边形模型。此时，ADIC表示的变量Y及与其具有相关性的区间变量对输出变量概率盒的贡献度之和。

图10 场景2中输入区间变量的平均距离指标贡献度

相较于场景1，考虑区间变量相关性后各变量贡献度的排序略有变化，占据主要贡献度的变量是相同的。对于输出变量75-118，场景2中区间变量的排序为：100＞90＞80＞116＞27＞35＞45＞15，排序后3位的变量发生变化，而排序前4位的变量与场景1相同。对于场景1和场景2中占据主要贡献度的变量是相同的这一现象可结合图6～图8中曲线进行解释。虽然输出变量的概率盒在计及区间变量相关性后变窄，但其减少程度占据整体波动范围的比例仍较小。因此，相较于相关性，区间变量本身的波动范围是影响输出变量的主要因素。此外，考虑相关性后使得原本具有较大影响区间变量的贡献度进一步增加，更加突出其对输出变量的影响程度。

区间变量贡献度排序结果可为电力系统规划和运行提供新的视角。例如，为了减少变量75-118的波动范围可以提高节点100、90、80和116处风电场功率预测精度或者安装储能装置抑制其出力的不确定性。通过上述措施还可减少变量77-78和82-96的波动范围，但对于变量23-24和49-69影响较小。

5.5 计算效率对比

本小节对DLMCS、CAM和DLSM的计算效率进行对比。对于所提DLSM，其总耗时total由构建上、下层代理模型RBF和sPCE构成。RBF包含了次非线性潮流计算、扩展系数求取及进行RBF×RBF次评估的时间。sPCE包含展开系数求解与进行次评估的时间。sPCE与输出变量的个数有关，可进一步表示为sPCE=sPCE，其中sPCE表示对单个输出变量构建sPCE模型及评估所需时间；表示输出变量的个数。

以系统中所有PQ节点的电压幅值（共64个PQ节点）作为输出变量，三种方法所需计算时间见表5。由结果可知，DLMCS由于需要进行4×106次确定性潮流计算，所需计算时间最多。CAM在获取单个变量概率盒时仅需27.1s，但随着输出变量个数的增加，其计算效率也有所下降。本文所提DLSM相较于DLMCS和CAM计算效率分别提升约126倍和14倍。对于DLSM来说，仅在构建上层代理模型时需要进行少量的非线性潮流计算，并且上层代理模型得到输入随机变量和输出变量边界值样本后，下层代理模型仅需要sPCE=1.48s（94.7/64）便可获得单个输出变量的概率盒。因此，所提方法在输出变量较多时效率方面具有更为明显的优势。

表5 不同方法的计算效率

Tab.5 The computational burden of different methods

6 结论

本文针对电力系统中随机与区间变量共存时的不确定潮流计算问题，提出一种基于双层代理模型的概率-区间潮流计算方法，并基于所提方法开展了对输入区间变量的灵敏度分析。通过算例分析验证了所提方法的精确性和快速性，得到主要结论如下：

1）验证了RBF模型能够以较高的精度近似非线性潮流方程中输入变量与输出响应的关系，较为精确地获取输出响应边界值。该代理模型是本文所提方法计算精度与效率的关键所在。

2）在计算精度方面，本文所提出的DLSM能够准确地获取输出变量的概率盒，且计算误差与CAM相近；在计算效率方面，DLSM相较于DLMCS和CAM均有较大幅度的提升，且在输出变量较多时具有更为明显的优势。

3）通过灵敏度分析可识别出对输出变量影响较大的区间变量，并得到其贡献度排序。对于特定输出变量，考虑区间变量相关性与否会引起贡献度排序略有差异，而占据主要贡献度的变量仍是相同的。

[1] Jordehi A R. How to deal with uncertainties in electric power systems? a review[J]. Renewable and Sustainable Energy Reviews, 2018, 96: 145-155.

[2] 康重庆, 姚良忠. 高比例可再生能源电力系统的关键科学问题与理论研究框架[J]. 电力系统自动化, 2017, 41(9): 2-11.

Kang Chongqing, Yao Liangzhong. Key scientific issues and theoretical research framework for power systems with high proportion of renewable energy[J]. Automation of Electric Power Systems, 2017, 41(9): 2-11.

[3] Prusty B R, Jena D. A critical review on probabilistic load flow studies in uncertainty constrained power systems with photovoltaic generation and a new approach[J]. Renewable and Sustainable Energy Reviews, 2017, 69: 1286-1302.

[4] 廖小兵, 刘开培, 乐健, 等. 电力系统区间潮流计算方法综述[J]. 中国电机工程学报, 2019, 39(2): 447-458.

Liao Xiaobing, Liu Kaipei, Le Jian, et al. Review on interval power flow calculation methods in power system[J]. Proceedings of the CSEE, 2019, 39(2): 447-458.

[5] 谢桦, 任超宇, 郭志星, 等. 基于聚类抽样的随机潮流计算[J]. 电工技术学报, 2020, 35(23): 4940-4948.

Xie Hua, Ren Chaoyu, Guo Zhixing, et al. Stochastic load flow calculation method based on clustering and sampling[J]. Transactions of China Electrotechnical Society, 2020, 35(23): 4940-4948.

[6] 韩佶, 苗世洪, 李超, 等. 计及相关性的电-气-热综合能源系统概率最优能量流[J]. 电工技术学报, 2019, 34(5): 1055-1067.

Han Ji, Miao Shihong, Li Chao, et al. Probabilistic optimal energy flow of electricity-gas-heat integrated energy system considering correlation[J]. Transactions of China Electrotechnical Society, 2019, 34(5): 1055-1067.

[7] Wang Chenxu, Liu Chengxi, Tang Fei, et al. A scenario-based analytical method for probabilistic load flow analysis[J]. Electric Power Systems Research, 2020, 181: 106193.

[8] 李聪聪, 王彤, 相禹维, 等. 基于改进高斯混合模型的概率潮流解析方法[J]. 电力系统保护与控制, 2020, 48(10): 146-155.

Li Congcong, Wang Tong, Xiang Yuwei, et al. Analytical method based on improved Gaussian mixture model for probabilistic load flow[J]. Power System Protection and Control, 2020, 48(10): 146-155.

[9] 何琨, 徐潇源, 严正, 等. 基于稀疏混沌多项式展开的孤岛微电网概率潮流计算[J]. 电力系统自动化, 2019, 43(2): 67-75.

He Kun, Xu Xiaoyuan, Yan Zheng, et al. Probabilistic

power flow calculation of islanded microgrid based on sparse polynomial chaos expansion[J]. Automation of Electric Power Systems, 2019, 43(2): 67-75.

[10] 胡健, 付立军, 马凡, 等. 基于仿射算术优化的不确定系统区间潮流快速分解法[J]. 电工技术学报, 2016, 31(23): 125-131.

Hu Jian, Fu Lijun, Ma Fan, et al. Fast decoupled power flow calculation of uncertainty system based on interval affine arithmetic optimization[J]. Transactions of China Electrotechnical Society, 2016, 31(23): 125-131.

[11] 杜萍静, 杨明, 曹良晶, 等. 含电压源换流器交直流系统的仿射潮流算法[J]. 电工技术学报, 2020, 35(5): 1106-1117.

Du Pingjing, Yang Ming, Cao Liangjing, et al. Affine power flow algorithm for AC/DC systems with voltage source converter[J]. Transactions of China Electrotechnical Society, 2020, 35(5): 1106-1117.

[12] 廖小兵, 刘开培, 张亚超, 等. 基于区间泰勒展开的不确定性潮流分析[J]. 电工技术学报, 2018, 33(4): 750-758.

Liao Xiaobing, Liu Kaipei, Zhang Yachao, et al. Uncertain power flow analysis based on interval Taylor expansion[J]. Transactions of China Electrotechnical Society, 2018, 33(4): 750-758.

[13] Luo Jinqing, Shi Libao, Ni Yixin. Uncertain power flow analysis based on evidence theory and affine arithmetic[J]. IEEE Transactions on Power Systems, 2018, 33(1): 1113-1115.

[14] 鲍海波, 韦化, 郭小璇, 等. 考虑风电不确定性的概率区间潮流模型与算法[J]. 中国电机工程学报, 2017, 37(19): 5633-5642.

Bao Haibo, Wei Hua, Guo Xiaoxuan, et al. Model and algorithm of probabilistic interval power flow considering wind power uncertainty[J]. Proceedings of the CSEE, 2017, 37(19): 5633-5642.

[15] Guo Xiaoxuan, Gong Renxi, Bao Haibo, et al. Hybrid stochastic and interval power flow considering uncertain wind power and photovoltaic power[J]. IEEE Access, 2019, 7: 85090-85097.

[16] Wang Chun, Ao Xin, Fu Wenbin. Three-phase power flow calculation considering probability and interval uncertainties for power distribution systems[J]. IET Generation, Transmission & Distribution, 2019, 13(15): 3334-3345.

[17] Wang Chenxu, Liu Dichen, Tang Fei, et al. A clustering-based analytical method for hybrid probabilistic and interval power flow[J]. International Journal of Electrical Power & Energy Systems, 2021, 126: 106605.

[18] Hu Xiaoyun, Zhao Xia, Feng Xinxin. Probabilistic-interval energy flow analysis of regional integrated electricity and gas system considering multiple uncertainties and correlations[J]. IEEE Access, 2019, 7: 178209-178223.

[19] 孙鑫, 王博, 陈金富, 等. 基于稀疏多项式混沌展开的可用输电能力不确定性量化分析[J]. 中国电机工程学报, 2019, 39(10): 1-10.

Sun Xin, Wang Bo, Chen Jinfu, et al. Sparse polynomial chaos expansion based uncertainty quantification for avaiable transfer capability[J]. Proceedings of the CSEE, 2019, 39(10): 1-10.

[20] 胡潇云, 赵霞, 冯欣欣. 基于稀疏多项式混沌展开的区域电-气联合系统全局灵敏度分析[J]. 电工技术学报, 2020, 35(13): 2805-2816.

Hu Xiaoyun, Zhao Xia, Feng Xinxin. Global sensitivity analysis for regional integrated electricity and gas system based on sparse polynomial chaos expansion[J]. Transactions of China Electrotechnical Society, 2020, 35(13): 2805-2816.

[21] 鲍海波, 郭小璇. 求解含风电相关性区间潮流的仿射变换最优场景法[J]. 电力系统保护与控制, 2020, 48(18): 114-122.

Bao Haibo, Guo Xiaoxuan. Optimal scenario algorithm based on affine transformation applied interval power flow considering correlated wind power[J]. Power System Protection and Control, 2020, 48(18): 114-122.

[22] Ran Xiaohong, Leng Shipeng, Liu Kaipei. A novel affine arithmetic method with missed the triangular domain with uncertainties[J]. IEEE Transactions on Smart Grid, 2020, 11(2): 1430-1439.

[23] 陈丽娜, 张智晟, 于道林. 基于广义需求侧资源聚合的电力系统短期负荷预测模型[J]. 电力系统保护与控制, 2018, 46(15): 45-51.

Chen Lina, Zhang Zhisheng, Yu Daolin. Short-term load forecasting model of power system based on generalized demand side resources aggregation[J]. Power System Protection and Control, 2018, 46(15): 45-51.

[24] 郭茜, 匡洪海, 王建辉, 等. 单机风电功率人工智能预测模型综述[J]. 电气技术, 2020, 21(2): 1-6.

Guo Qian, Kuang Honghai, Wang Jianhui, et al. Summary of artificial intelligence prediction model for single wind power[J]. Electrical Engineering, 2020, 21(2): 1-6.

[25] Baghaee H R, Mirsalim M, Gharehpetian G B, et al. Fuzzy unscented transform for uncertainty quantification of correlated wind/PV microgrids: possibilistic-probabilistic power flow based on RBFNNs[J]. IET Renewable Power Generation, 2017, 11(6): 867-877.

[26] Li Xu, Gong Chunlin, Gu Liangxian, et al. A sequential surrogate method for reliability analysis based on radial basis function[J]. Structural Safety, 2018, 73: 42-53.

[27] Sun Xin, Tu Qingrui, Chen Jinfu, et al. Probabilistic load flow calculation based on sparse polynomial chaos expansion[J]. IET Generation, Transmission & Distribution, 2018, 12 (11): 2735-2744.

[28] Zimmerman R D, Murillo-Sanchez C E, Thomas R J. Matpower: steady-state operations, planning and analysis tools for power systems research and education[J]. IEEE Transactions on Power Systems, 2011, 26(1): 12-19.

[29] Bi Sifeng, Broggi M, Wei Pengfei, et al. The Bhattacharyya distance: enriching the P-box in stochastic sensitivity analysis[J]. Mechanical Systems and Signal Processing, 2019, 129: 265-281.

Probabilistic-Interval Power Flow and Sensitivity Analysis Using Double Layer Surrogate Method

Wang Chenxu1Tang Fei1Liu Dichen1Gao Xin1Zhou Yixi2

（1. School of Electrical Engineering and Automation Wuhan University Wuhan 430072 China 2. Hangzhou Power Supply Company State Grid Zhejiang Electric Power Co. Ltd Hangzhou 310000 China）

With the increasing complexity of uncertainties in power systems, the co-existence of random and interval variables makes it challenging to accurately obtain systems’ operating states by using probabilistic or interval power flow calculations. To cope with this issue, this paper proposes a double layer surrogate method for probabilistic-interval power flow analysis. The proposed method can obtain the upper and lower surrogate models based on a few times deterministic power flow calculations. Then, the surrogate models are used to solve a large number of deterministic power flow calculations required in the conventional probabilistic-interval power flow method, and the probability boxes of output variables can be obtained with high efficiency. In addition, a sensitivity index is proposed to describe the characteristics of output variables, and the sensitivity analysis is performed using the proposed method to identify the importance of input interval variables that affect output variables. The accuracy and efficiency of the proposed method are validated on the IEEE 118-bus test system by comparing the existing methods. The sensitivity analysis can identify the key interval variables for the output variables and reveal the relationship between the operating states and the input interval variables.

Uncertain power flow, random variable, interval variable, surrogate model, sensitivity analysis

10.19595/j.cnki.1000-6753.tces.210171

TM744

国家自然科学基金资助项目（51977157）。

2021-02-01

2021-03-21

王晨旭 男，1996年生，博士研究生，研究方向为电力系统不确定性量化分析。E-mail：Chenxu_Wang2021@163.com

唐 飞 男，1982年生，副教授，博士生导师，研究方向为电力系统运行与控制。E-mail：tangfei@whu.edu.cn.com（通信作者）

（编辑 赫蕾）