本发明涉及一种电-气耦合系统概率连续潮流建模及其求解方法，属于多能流系统安全稳定运行领域。
背景技术：
在能源互联的背景下，电力系统、天然气系统彼此之间的耦合不断加强，电力、天然气等能源需求量日益增大。风电作为可再生清洁能源，被广泛地应用于地区发电，但由于风机的有功出力受风速的影响而具有随机性和不确定性等特征，为了保证多能流系统的安全高效稳定运行，在含风电电力系统的电-气耦合系统的静态安全稳定分析的问题中，对风电不确定性的分析以及电-气耦合系统概率连续潮流建模亟待研究。基于此，提出一种电-气耦合系统概率连续潮流建模及其求解方法对多能流系统的安全稳定运行具有现实意义。
文献《计及电-气互联能源系统安全约束的可用输电能力计算》考虑了电气耦合系统安全约束，提出一种基于线性预测的连续潮流法进行求解，用线性预测法辨识制约电力系统可输电能力的关键约束，对电力系统可输电能力进行了计算；文献《一种基于连续潮流的在线静态稳定综合评估方法》通过对电压可控节点的优化，来达到提高系统负荷裕度的要求；文献《newparameterizationschemesforthecontinuationloadflowmethod》考虑了电力系统保护和自动装置动作对于负荷裕度计算的影响，进而应用了一种在线计算的连续潮流算法，来计算系统的稳定裕度。目前，电力系统的静态安全分析方面已有广泛研究，但随着电-气系统的高度耦合，考虑电-气耦合系统的静态安全稳定分析的鲜有研究，所以，需要研究一种电-气耦合系统概率连续潮流建模及其求解方法。
为了保证多能流系统安全稳定运行，研究电-气耦合系统概率连续潮流建模及其求解方法十分必要，在能源互联网背景下有着重要意义。
技术实现要素：
针对现有技术的不足，本发明“一种电-气耦合系统概率连续潮流建模及其求解方法”，建立电-气耦合系统概率连续潮流模型，应用预测-校正的方法来逐渐逼近系统临界运行点来计算电-气耦合系统中电力网络及天然气网络的最大负荷裕度。
本发明采用如下技术方案：一种电-气耦合系统概率连续潮流建模及其求解方法，该方法包括如下步骤：
步骤1：建立电-气耦合系统概率连续潮流模型；
步骤2：应用预测-校正的方法来逐渐逼近系统临界运行点来计算电-气耦合系统中电力网络及天然气网络的最大负荷裕度。
附图说明
图1是本发明的算法流程图。
具体实施方式
本发明包括以下步骤：
步骤1：建立电-气耦合系统概率连续潮流模型；
本发明研究风电场出力的不确定性，采用二参数的威布尔分布模型来描述风速v的随机变化规律，风力发电机的输出功率pw。关于风电机组的建模已被广泛研究，本发明将直接应用，细节不再赘述。
下式表示的是下式表示的是含风电电-气耦合系统含参数潮流方程组：
式中，pgi、pli、qgi、qli分别代表初始点发电机发出的有功功率、有功负荷、无功功率以及无功负荷；vi、θi分别表示节点电压和节点相角；a为去掉压缩机管道的支路节点关联矩阵；u为压缩机节点关联矩阵；ω为天然气气体输入量；tτ为压缩机气体需求量；ηtrans、ηchp.e、ηchp.th、ηfur、ηex，分别为能源集线器(energyhub，eh)中的变压器、chp发电部分、chp供热部分、热锅炉和热交换机的转化效率；pe、pg分别为eh输入的电功率和天然气。
电-气耦合系统概率连续潮流模型简记为：
步骤2：应用预测-校正的方法来逐渐逼近系统临界运行点来计算电气耦合系统中电力网络及天然气网络的最大负荷裕度。
本发明在连续潮流方程中计及负荷参数λe,λg。使得λe∈[0,λec]，λg∈[0,λgc],其中λec和λgc分别为电力网络和天然气网络在临界点时的负荷增长系数。在λ＝0时对应于原始运行点的负荷水平，在λ＝λc时对应于在给定负荷增长方式下静态稳定临界点的负荷水平。
(1)预测计算
1)对含参数潮流方程组求微分得到：
由于潮流方程中引入了参数λe,λg，增加了两个未知状态变量，而需增设的两个一维向量对上式构成的雅可比矩阵加入两个对应的行向量
式中x表示电-气耦合系统的状态变量v，θ，fg，p，pe，pg，λe,λg，连续潮流预测步的计算公式为：
切向量求出后，下一步的预测值由下式给出。
式中x％、λe％、λg％为预估值，x0、λe0、λg0为初始值，σ为预测步长。
(2)校正计算
在校正步计算中，采用垂直校正的方法，以[x％，λe％，λg％]t为初值，应用牛顿-拉夫逊法求解下述校正方程。
式中δx、δλe、δλg为修正量，在校正计算中，通过牛顿-拉夫逊计算得到校正值作为下一次预测部分计算的初值，一直迭代至系统临界运行点，由此来计算电-气耦合系统中电力网络及天然气网络的最大负荷裕度。
据此就得到一种电-气耦合系统概率连续潮流建模及其求解方法。
以上实施方案仅用于说明本发明，而并非对本发明的限制，有关技术领域的普通技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型，因此所有等同的技术方案也属于本发明的保护范畴。