CN114580123B - 一种电力天然气互联系统分布式协同校正控制方法 - Google Patents

Abstract

本发明提出了一种电力天然气互联系统分布式协同校正控制方法。首先，基于电力系统交流潮流模型以及通过正交配置法离散的天然气系统动态管流模型构建电力‑天然气互联系统协同校正控制优化模型。其次，考虑到区域能源系统的运营独立性，将电力‑天然气互联系统协同校正优化控制模型解耦，并通过辅助问题原则法构建区域系统辅助优化子模型。最后，基于原始对偶拟牛顿算法分布式求解区域系统辅助优化子模型得到协同校正控制策略。本发明优点在于：通过正交配置法离散的天然气系统动态管流模型兼顾计算精度和效率；分布式求解协同校正控制策略，有助于维护运营独立性、缓解通信阻塞；具有有限时间收敛性，可以快速获得高精度的协同校正控制策略。

Description

一种电力天然气互联系统分布式协同校正控制方法

技术领域

本发明属于综合能源系统技术领域，具体为一种电力天然气互联系统分布式协同校正控制方法。

背景技术

燃气机组灵活调峰支撑下的可再生能源发电以及发电企业的煤改气转型，是短期内减少碳排放的有效方式。加速构建电力-天然气互联系统有助于“碳达峰、碳中和”目标的实现，但随着电力-天然气互联系统异质能流间耦合加深，不同系统间故障传播的风险随之提升。高效精准的在线协同运行与控制方法对于保障电力-天然气互联系统安全运行具有重要意义。其中，协同校正控制是确保电力-天然气互联系统安全运行的重要控制措施之一，其主要目标是消除单一预想事故场景的运行约束违反现象。当电力-天然气互联系统中发生N-1或N-k事故处于紧急状态时，部分安全限值约束条件被破坏，出现电力系统支路潮流越限、电力节点电压越限、输气管道压力越限等不安全运行状态。此时，通过电力-天然气互联系统协同校正控制进行电力系统常规机组及燃气机组发电机功率的再调整、电负荷切除、天然气系统气源再调整、气负荷切除、压缩机调整等校正控制措施，可以使得电力-天然气互联系统在允许的校正控制时间内重新过渡至安全运行状态；另一方面，电力-天然气互联系统的协同校正控制可作为预防控制的补充，对于预想故障集中重点关注的N-1或N-k事故场景提供更精准的控制策略。

本发明为一种电力-天然气互联系统分布式协同校正控制方法。本发明基于电力系统交流潮流模型以及天然气系统动态管流模型提出，尤其在对偏微分方程组描述的天然气系统动态管流模型的处理上，本发明采用正交配置法将其离散为配置点处的有限维线性代数方程模型，通过较少的离散配置点得到高精度的拟合结果，使本发明更适用于在线计算。针对以电力-天然气互联系统跨越地域广、投资运营主体多、利益关系复杂的特征，本发明能够在维持多主体运营机制，维护运营决策独立性和私密性的基础上，通过有限信息的交互，实现兼顾计算精度和计算效率的协同校正控制策略求解。本发明对电力-天然气互联系统的安全运行具有重要意义。

发明内容

本发明提出了一种电力-天然气互联系统分布式协同校正控制方法。本发明结合工程实际，适用于不同利益主体运营的电力-天然气互联系统(例如跨国能源互联网)，具体为一种针对电力-天然气互联系统中N-1或N-k故障造成的电力系统支路潮流越限、电力节点电压越限、天然气系统输气管道压力越限、天然气节点压力或流量越限事故的协同校正控制策略分布式生成方法。具体而言，本发明通过构建电力-天然气互联系统协同校正控制优化模型、按照电力-天然气互联系统的地理区划将电力-天然气互联系统协同校正控制优化模型解耦为多个区域系统辅助优化子模型、分布式地求解区域系统辅助优化子模型，最终获得可消除电力-天然气互联系统中越限事故的协同校正控制策略。

本发明中电力-天然气互联系统协同校正控制优化模型基于电力系统交流潮流模型以及天然气系统动态管流模型构建，尤其在对偏微分方程组描述的天然气系统动态管流模型的处理上，本发明通过正交配置法将其离散为配置点处的有限维线性代数方程模型。正交配置法通过较少的离散配置点可以得到高精度的计算结果，减小了本发明的求解规模和难度，增强了本发明的在线应用潜力。

本发明中区域系统辅助优化子模型的求解过程基于原始对偶拟牛顿算法设计。本发明尤其将原始对偶拟牛顿算法的应用范围由电力系统的分布式计算进一步拓展到电力-天然气互联系统的分布式计算。通过原始对偶拟牛顿算法求解区域系统辅助优化子模型的优点在于不需要中心协调机构、仅需通过有联络线连接区域能源系统间(以下简称邻居)近邻通信的协调。此外，针对区域系统辅助优化子模型条件数大、通过梯度算法求解效率低的问题，原始对偶拟牛顿算法在原始和对偶更新中通过将迭代方向由近端梯度方向改进为近端牛顿方向，从而使本发明中区域系统辅助优化子模型求解过程的收敛速度大幅提升，增强了本发明的在线应用潜力。

本发明的上述技术问题主要是通过下述技术方案得以解决的：

一种电力天然气互联系统分布式协同校正控制方法，其特征在于，包括以下步骤：

步骤1：通过以电力-天然气互联系统总体校正控制代价最小化的目标函数构建电力-天然气互联系统协同校正控制优化模型；以电力系统的常规发电机功率调整量、燃气发电机功率调整量、电负荷切除量以及天然气系统气源调整量和气负荷切除量为电力-天然气互联系统协同校正控制优化模型的连续性优化变量；通过电力-天然气互联系统协同校正控制优化模型的电力系统运行约束条件、电力-天然气互联系统协同校正控制优化模型的电力系统安全限值约束条件、电力-天然气互联系统协同校正控制优化模型的通过正交配置法离散的天然气系统运行约束条件、电力-天然气互联系统协同校正控制优化模型的天然气系统安全限值约束条件、电力-天然气互联系统协同校正控制优化模型的天然气系统边界约束条件、电力-天然气互联系统协同校正控制优化模型的电力-天然气耦合约束条件构建电力-天然气互联系统协同校正控制优化模型的约束条件，进一步构建电力-天然气互联系统协同校正控制优化模型；

步骤2：将电力-天然气互联系统协同校正控制优化模型按照电力-天然气互联系统的地理区划通过节点复制法解耦得到多个区域系统校正控制子模型；在得到的区域系统校正控制子模型中将待协调的耦合变量通过增广处理统一向量维度得到增广耦合变量；将每个区域系统校正控制子模型和每个区域系统对应的增广耦合变量，通过辅助问题原则法构建区域系统辅助优化子模型。

步骤3：以故障后电力-天然气互联系统运行状态作为步骤2的区域系统辅助优化子模型中决策变量的初值，基于原始对偶拟牛顿算法分布式求解区域系统辅助优化子模型，得到电力-天然气互联系统的协同校正控制策略；

作为优选，步骤1所述以电力-天然气互联系统总体校正控制代价最小化的目标函数，具体为：

设电力-天然气互联系统总体校正控制代价为C，可表示为N个区域电力系统控制代价C_e,i(i＝1,2,……,N)和M个区域天然气系统控制代价C_g,i(i＝1,2,……,M)之和：

式中，P_L0,n为校正前区域电力系统n的电力负荷有功功率向量，P_L,n为区域电力系统n的切负荷决策变量，M_E,n为区域电力系统n切除电负荷的成本系数；P_T0,n为校正前区域电力系统n的燃气发电机有功出力向量，P_T,n为区域电力系统n的燃气发电机调整决策变量，E_GT,n为区域电力系统n燃气发电机调整的成本系数；P_G0,n为校正前区域电力系统n的常规发电机有功出力向量，P_G,n为区域电力系统n的常规发电机调整决策变量，E_G,n为区域电力系统i常规发电机调整的成本系数；N表示区域电力系统的数量；f_S0,m为校正前区域天然气系统m气源的质量流率向量，为区域天然气系统m气源的质量流率决策变量，E_S,m区域天然气系统m气源调整的成本系数；f_L0,m为校正前区域天然气系统m气负荷的质量流率向量，为区域天然气系统m常规气负荷的质量流率决策变量，M_G,m区域天然气系统m切除常规气负荷的成本系数；f_T0,m为校正前区域天然气系统m燃气发电机气负荷的质量流率向量，为区域天然气系统m燃气发电机气负荷的质量流率决策变量，E_T,m区域天然气系统m燃气发电机气负荷调整的成本系数；M表示区域天然气系统的数量。

步骤1所述电力-天然气互联系统协同校正控制优化模型的电力系统运行约束条件，具体为：

设Λ为电力-天然气互联系统中电力系统节点的集合，且电力系统满足二阶锥交流潮流模型：

式中，P_Gi为节点i注入的常规发电机有功功率；Q_Gi为节点i注入的常规发电机无功功率；P_Ti为节点i注入的燃气发电机有功功率；Q_Ti为节点i注入的燃气发电机无功功率；P_Li为节点i流出的负荷有功功率；Q_Li为节点i流出的负荷无功功率；δ(i)为和节点i直接连接的所有节点的集合；P_ij为从节点i流向节点j的有功功率；Q_ij为从节点i流向节点j的无功功率；P_ji为从节点j流向节点i的有功功率；Q_ji为从节点j流向节点i的无功功率；g_ij为线路i-j的电导；b_ij为线路i-j的电纳；b_i,ij为线路i-j的i侧对地电纳；b_j,ij为线路i-j的j侧的对地电纳；、V_i ^(c)为节点i的电压典型值；为节点i的电压典型值；辅助变量W_i＝V_i ²；辅助变量W_j＝V_j ²；辅助变量W_cij＝V_iV_j cosθ_ij；辅助变量W_sij＝V_iV_j sinθ_ij；

步骤1所述电力-天然气互联系统协同校正控制优化模型的电力系统安全限值约束条件，具体为：

式中，Λ_G为常规发电机节点的集合，Λ_T为燃气发电机节点的集合，Λ_L为电负荷节点的集合，Λ为电力系统所有节点集合；为常规发电机i有功出力的上限；P _Gi为常规发电机i有功出力的下限；为常规发电机i无功出力的上限；Q _Gi为常规发电机i无功出力的下限；为常规发电机i爬坡率的上限；R _G,i为常规发电机i爬坡率的下限；为燃气发电机i有功出力的上限；P _Ti为燃气发电机i有功出力的下限；为燃气发电机i无功出力的上限；Q _Ti为燃气发电机i无功出力的下限；为燃气发电机i爬坡率的上限；R _T,i为燃气发电机i爬坡率的下限；为节点i有功负荷切除调整的上限，也即节点i故障前的有功负荷大小；P _Li为节点i有功负荷切除调整的下限；为节点i无功负荷切除调整的上限，也即节点i故障前的无功负荷大小；Q _Li为节点i无功负荷切除调整的下限；为系统有功网损的上界；为节点i电压幅值的上限；V _i为节点i电压幅值的下限。

步骤1所述电力-天然气互联系统协同校正控制优化模型的通过正交配置法离散的天然气系统运行约束条件，具体为：

设Υ为天然气系统中节点的集合，天然气系统满足通过正交配置法离散的动态管流模型：

式中，Δt为天然气系统待优化的时间长度；Δx_ab为输气管道ab的长度；A_ab为输气管道ab的横截面积；D_ab为输气管道ab的管壁内径；N_t为时间正交配置点的个数；N_x为空间正交配置点的个数；为空间一阶微分矩阵的第(i，k)个元素，为时间一阶微分矩阵的第(j，k)个元素；f_ab,ij为输气管道ab的质量流率在第(i，j)个时空配置点处的取值；p_ab,ij为输气管道ab的气体压力在第(i，j)个时空配置点处的取值；λ为管壁摩擦系数；为输气管道ab的平均气体流速；v_s为天然气的气体音速；f_Sa,j为天然气系统节点a气源的质量流率在第j个时间配置点处的取值；f_La,j天然气系统节点a常规气负荷的质量流率在第j个时间配置点处的取值；f_Ta,j天然气系统节点a燃气发电机气负荷的质量流率在第j个时间配置点处的取值；为压缩机ba出口处质量流率在第j个时间配置点处的取值；f_Cab,0j压缩机ab入口处质量流率在第j个时间配置点处的取值。

步骤1所述电力-天然气互联系统协同校正控制优化模型的天然气系统安全限值约束条件，具体为：

式中，i＝1,…,N_x，为空间域正交配置插值点的个数；j＝1,…,N_t，N_t为时域正交配置插值点的个数；为管道ab气体压力的上限；p_ab为管道ab气体压力的下限；为管道ab可输送气体质量流率的上限；f _ab为管道ab可输送气体质量流率的下限；为节点a气源流量的上限；f _Sa为节点a气源流量的下限；为校正控制前节点a的气负荷质量流率；f _La为节点a气负荷质量流率的下限；为节点a燃气机组气负荷质量流率可调节的上限；f _Ta分别为节点a燃气机组气负荷质量流率可调节的下限。

步骤1所述电力-天然气互联系统协同校正控制优化模型的天然气系统边界约束条件，具体为：

其中，ρ_ab,ij表示管道ab中气体密度在第(i，j)个时空配置点处的取值；和分别表示管道ab在初始时刻第i个空间配置点处的取值；f_Sa,j为天然气系统节点a气源的质量流率在第j个时间配置点处的取值；f_La,j天然气系统节点a常规气负荷的质量流率在第j个时间配置点处的取值；f_Ta,j天然气系统节点a燃气发电机气负荷的质量流率在第j个时间配置点处的取值；为初始时刻气节点a气源的质量流率；为初始时刻节点a气负荷的质量流率；为初始时刻节点a燃气发电机气负荷的质量流率；校正控制策略中节点a气源质量流率的给定值；为校正控制策略中节点a气负荷质量流率的给定值；为校正控制策略中节点a燃气机组进气量质量流率的给定值。

步骤1所述电力-天然气互联系统协同校正控制优化模型的电力-天然气耦合约束条件，具体为：

设电力节点i∈Λ的燃气发电机是天然气节点m∈Υ的气负荷，则电力节点i与天然气节点m的耦合约束条件为：

式中，为校正控制策略中天然气节点a燃气机组进气量质量流率的给定值；P_Ti为电力节点i注入的燃气发电机有功功率；η_Ti,a为燃气发电机的工作效率，H_G为天然气热值。

设压缩机ab，a,b∈Υ是电力节点i∈Λ的电负荷，则其耦合约束条件为：

式中，P_Ci为电网节点i上电压缩机消耗的电功率；η_Ci,ab压缩机ab消耗电功率与其原动机驱动流量间的比例系数；为压缩机ab出口处压力在时间配置点l处的值；p_ab,0l为为压缩机ab入口处压力在时间配置点l处的值；K_Cab为压缩机ab的增压比；f_Cab,kl为压缩机ab在第(k，l)个时空配置点处质量流率的值；f _Cab为压缩机ab质量流率的下限；为压缩机ab质量流率的上限；f_Cab,k0为压缩机ab的空间配置点k处质量流率在初始时间配置点处的值；为压缩机ab的空间配置点k处的质量流率在初始时刻值。

步骤1所述电力-天然气互联系统协同校正控制优化模型，决策变量包括所有与约束条件和目标函数有关的待确定的量。电力系统中的决策变量具体为：电力节点i注入的常规发电有功功率P_Gi；电力节点i注入的常规发电机无功功率Q_Gi；电力节点i注入的燃气发电有功功率P_Ti；电力节点i注入的燃气发电机无功功率Q_Ti；电力节点i流出的负荷有功功率P_Li；电力节点i流出的负荷无功功率Q_Li；从电力节点i流向电力节点j的有功功率P_ij；从电力节点i流向电力节点j的无功功率Q_ij；辅助变量W_i；辅助变量W_cij；辅助变量W_sij。天然气系统中的决策变量具体为：输气管道ab的质量流率在第(i，j)个时空配置点处的取值f_ab,ij；输气管道ab的气体压力在第(i，j)个时空配置点处的取值p_ab,ij；天然气系统节点a气源的质量流率在第j个时间配置点处的取值f_Sa,j；天然气系统节点a常规气负荷的质量流率在第j个时间配置点处的取值f_La,j；天然气系统节点a燃气发电机气负荷的质量流率在第j个时间配置点处的取值f_Ta,j；压缩机ba出口处质量流率在第j个时间配置点处的取值压缩机ab入口处质量流率在第j个时间配置点处的取值f_Cab,0j。

作为优选，在上述的一种电力-天然气互联系统分布式协同校正控制方法的步骤2中，所述具体方法如下：

步骤2所述将电力-天然气互联系统协同校正控制优化模型按照电力-天然气互联系统的地理区划通过节点复制法解耦得到多个区域系统校正控制子模型，具体为：

定义区域电力系统为区域内部电力节点和输电线路的集合，区域天然气系统为区域内部天然气节点和输气管道的集合。不同的区域电力系统之间通过跨区电力联络线连接，不同的区域天然气系统之间通过跨区输气管道连接，区域电力系统与区域天然气系统之间通过燃气发电机和电力驱动的加压站耦合。区域电力系统和区域天然气系统统称为区域系统。

以两个通过电力联络线连接的区域电力系统的解耦为例，对节点复制法和所得到的区域系统校正控制子模型进行说明。假设区域电力系统m与区域电力系统n通过联络线k-l连接，其中节点k位于区域电力系统m，节点l位于区域电力系统n；

基于步骤1中所定义的集中优化模型，对区域电力系统m而言联络线k-l涉及的决策变量包括y_k＝[W_k,θ_k,P_lk,Q_lk,W_ckl,W_skl]，对区域电力系统n而言联络线k-l涉及决策变量包括y_l＝[W_l,θ_l,P_lk,Q_lk,W_clk,W_slk]；由于节点l的决策变量无法直接参与区域电力系统m的优化，节点k上的决策变量无法直接参与区域电力系统n的优化，为了将步骤1所述电力-天然气互联系统协同校正控制优化模型解耦，在区域电力系统m中引入虚拟节点l’，区域电力系统n中引入节点k’。

其中区域电力系统m中虚拟节点l’上的决策变量为y_l′＝[W_l′,θ_l′,P_l′,Q_l′,W_ckl′,W_skl′]，区域电力系统n中虚拟节点k’上的决策变量为y_k′＝[W_k′,θ_k′,P_k′,Q_k′,W_ck′l,W_sk′l]；

所述y_l与所述y_l′、所述y_k与所述y_k′互为副本变量；

则区域电力系统m与区域电力系统n之间耦合变量分别定义为v_m＝[y_k,y_l′]和v_n＝[y_k′,y_l]；

在引入虚拟节点后，原联络线k-l被替换为(虚拟)输电线路k-l’和(虚拟)输电线路k’-l。其中，(虚拟)输电线路k-l’参与区域电力系统m的优化，(虚拟)输电线路k’-l参与区域电力系统n的优化。

为了使区域电力系统m和区域电力系统n分别优化得到的(虚拟)输电线路k’-l和k-l’的决策变量结果对于实际联络线k-l可行，必须使(虚拟)输电线路k-l’和(虚拟)输电线路k’-l上对应的决策变量相等，也即需要满足一致性约束条件。具体而言，在解耦之后区域电力系统m和区域电力系统n的电力-天然气互联系统协同校正控制优化模型分别为：

对于区域电力系统m，在对区域电力系统m内部的节点和输电线路相关决策变量构建如步骤1所述电力-天然气互联系统协同校正控制优化模型的基础上，还需满足一致性约束条件：

其中，为区域电力系统m通过与区域电力系统n通信获得的y_l实际值。

对于区域电力系统n，在对区域电力系统n内部的节点和输电线路相关决策变量构建如步骤1所述电力-天然气互联系统协同校正控制优化模型的基础上，还需满足一致性约束条件：

其中，为区域电力系统n通过与区域电力系统m通信获得的实际值。

电力-天然气互联系统中所有电力联络线、输气管道和电力-天然气耦合设备均可基于上述节点复制法引入虚拟节点、副本变量、耦合变量和一致性约束条件，实现电力-天然气互联系统协同校正控制优化模型的解耦，得到多个区域系统校正控制子模型。

步骤2所述在得到的多个区域系统校正控制子模型中将待协调的耦合变量通过增广处理统一向量维度得到增广耦合变量，具体为：

对各个区域系统耦合变量v_n,n＝1,…,N插入0元素构建增广耦合变量使得各个区域系统中待协调的耦合变量向量维度相同；

假设共有M个耦合变量需要协调求解，则得到的增广耦合变量

步骤2所述将每个区域系统校正控制子模型和每个区域系统对应的增广耦合变量，通过辅助问题原则法构建区域系统辅助优化子模型，具体为：

所述区域系统辅助优化子模型包含区域系统原始辅助子优化模型和区域系统对偶辅助子优化模型；

所述区域系统n的原始辅助子优化模型如下：

式中，所有变量和系数的上角标为迭代轮次标识；为根据步骤2.2得到的增广耦合变量；x_n为区域系统n的内部决策变量；构成区域系统n子优化模型的全部决策变量；Ψ_n为区域系统n决策变量的可行域，根据步骤1中定义的电力系统运行约束条件和安全限值约束条件确定；c为罚系数，一般取为较大的正数；为与区域系统n收敛性相关的缩放系数；为区域系统n增广耦合变量对应的海塞矩阵的估计矩阵；为区域系统n增广耦合变量对应的副本变量，在迭代过程中基于邻居电力系统m∈Δ_n的通信进行更新。其中ξ_n,m为区域系统n与区域系统m的关联矩阵，其定义如下：

区域系统n的对偶辅助子优化问题如下：

式中，所有变量和系数的上角标为迭代轮次标识；为区域系统n拉格朗日乘子；为区域系统n拉格朗日乘子对应的k轮迭代海塞矩阵估计值；为区域系统n增广耦合变量对应的副本变量。

作为优选，所述步骤3具体如下：

步骤3.1，所述初始化区域系统辅助优化子模型的决策变量，初始化的原始对偶拟牛顿算法；

各个区域系统n∈{1,…,N}根据故障后电力-天然气互联系统运行状态设置区域系统n内部决策变量初值区域系统n增广耦合变量初值以及区域系统n副本变量初值原始对偶拟牛顿算法的参数初始化具体为：初始化区域系统n拉格朗日乘子的初值初始化原始对偶拟牛顿算法中区域系统n原始海塞矩阵初始化区域系统n增广耦合变量变化量初始化区域系统n增广耦合变量梯度值初始化区域系统n增广耦合变量梯度变化量初始化区域系统n拉格朗日乘子的海塞矩阵初始化区域系统n拉格朗日乘子的变化量初始化区域系统n拉格朗日乘子的梯度值初始化区域系统n拉格朗日乘子的梯度变化量初始化罚系数c；设置迭代次数k＝0；设置原始对偶拟牛顿算法收敛精度门槛值ε。

步骤3.2，利用初始化的原始对偶拟牛顿算法求解初始化的原始辅助子优化模型，得到各区域系统内部决策变量和增广耦合变量的迭代值；

设置迭代次数k＝k+1，各区域系统n∈{1,...,N}并行计算如下初始化的原始辅助子优化模型，得到各区域系统内部决策变量和增广耦合变量的第k次迭代值：

其中，为区域系统n内部决策变量的第k次迭代值；为区域系统n增广耦合变量的第k次迭代值；x_n为区域系统n的内部决策变量；为区域系统n的增广耦合变量；Ψ_n为区域系统n决策变量的可行域；为区域系统n增广耦合变量初值；为原始对偶拟牛顿算法中区域系统n原始海塞矩阵初值；为区域系统n拉格朗日乘子的初值；f_n为区域系统n的校正控制代价。

步骤3.3，所有邻居区域系统之间交换增广耦合变量的迭代值，各区域系统根据交换得到的邻居区域系统耦合变量的迭代值计算各区域系统副本变量的迭代值；

所有邻居区域系统之间交换增广耦合变量各区域系统n∈{1,...,N}根据下式计算区域系统n副本变量的第k次迭代值：

式中，为区域系统n副本变量的第k次迭代值；Δ_n为区域系统n邻居的集合；ξ_n,m为区域系统n与区域系统m的关联矩阵；为邻居区域系统m的增广耦合变量的第k次迭代值。

步骤3.4，各区域系统副本变量的迭代值带入初始化的对偶辅助子优化模型，计算各区域系统的拉格朗日乘子的迭代值；

各区域系统n∈{1,...,N}并行计算如下对偶辅助子优化模型，得到区域系统n拉格朗日乘子的第k次迭代值：

式中，为区域系统n拉格朗日乘子的第k次迭代值；为区域系统n拉格朗日乘子的初值；为区域系统n拉格朗日乘子的海塞矩阵初值；为区域系统n增广耦合变量的第k次迭代值；为区域系统n副本变量第k次迭代值。

步骤3.5，各区域系统利用步骤3.4得到各区域系统的拉格朗日乘子的迭代值，计算增广耦合变量梯度值的迭代值；

各区域系统n∈{1,...,N}根据下式计算区域系统n增广耦合变量梯度值的第k次迭代值

式中，为增广耦合变量梯度值的第k次迭代值；c为罚系数；为区域系统n副本变量第k次迭代值；为区域系统n增广耦合变量的第k次迭代值。

步骤3.6，各区域系统利用各区域系统的增广耦合变量的迭代值和各区域系统的增广耦合变量梯度值的迭代值，计算各区域系统的增广耦合变量变化量的迭代值和增广耦合变量梯度变化量的迭代值；

各区域系统n∈{1,...,N}通过下式计算区域系统n增广耦合变量变化量的第k次迭代值和区域系统n增广耦合变量梯度变化量的第k次迭代值：

式中，为区域系统n增广耦合变量变化量的第k次迭代值；D_n为区域系统n的正规化矩阵；为区域系统n增广耦合变量的第k次迭代值；为区域系统n增广耦合变量的第k-1次迭代值；为区域系统n耦合变量梯度变化量的第k次迭代值；为区域系统n的增广耦合变量梯度值的第k次迭代值；为区域系统n的增广耦合变量梯度值的第k-1次迭代值；σ为正则化参数，一般设置为较小的正数。

步骤3.7，各区域系统利用步骤3.6得到的各区域系统的增广耦合变量变化量的迭代值和增广耦合变量梯度变化量的迭代值计算各区域系统的增广耦合变量原始海塞矩阵的迭代值；

各区域系统基于如下分散式(Decentralized Broyden-Fletcher-Goldfarb-Shanno，D-BFGS)迭代计算区域系统n原始海塞矩阵的第k次迭代值。若则：

式中，为区域系统n原始海塞矩阵的第k次迭代值；为区域系统n原始海塞矩阵的第k-1次迭代值；为区域系统n耦合变量梯度变化量的第k次迭代值；为区域系统n增广耦合变量变化量的第k次迭代值；I为单位矩阵；σ为正则化参数，一般设置为较小的正数。

若则：

式中，为区域系统n原始海塞矩阵的第k次迭代值；为区域系统n原始海塞矩阵的第k-1次迭代值。

步骤3.8，各区域系统基于步骤3.7得到的各区域系统的增广耦合变量海塞矩阵的迭代值，计算各区域系统的耦合变量更新方向本地值的迭代值；

各区域系统n∈{1,...,N}根据下式并行计算区域系统n耦合变量更新方向本地值的第k次迭代值：

式中，为区域系统n耦合变量更新方向本地值的第k次迭代值；为区域系统n原始海塞矩阵的第k次迭代值；Ψ为正则化参数，一般设置为较小的正数；D_n为区域系统n的正规化矩阵；为区域系统n拉格朗日乘子的第k次迭代值；c为罚系数；为区域系统n增广耦合变量的第k次迭代值；为区域系统n副本变量第k次迭代值。

步骤3.9，所有邻居区域系统之间交换由步骤3.8得到的增广耦合变量的更新方向本地值的迭代值，各区域系统根据交换得到的邻居区域系统的增广耦合变量的更新方向本地值的迭代值计算各区域系统的增广耦合变量的近端值的迭代值；

所有邻居区域系统之间交换各区域系统n∈{1,...,N}由下式并行计算区域系统n增广耦合变量近端值的第k次迭代值

式中，为区域系统n增广耦合变量近端值的第k次迭代值；为区域系统n增广耦合变量的第k次迭代值；为区域系统n迭代步长；Δ_n为区域系统n邻居的集合；为区域系统n缩放系数；为区域系统n的邻居区域系统i耦合变量更新方向本地值的第k次迭代值。其中，迭代步长具有如下形式的闭式解：

式中，为区域系统n拉格朗日乘子的第k次迭代值；Δ_n为区域系统n邻居的集合；为区域系统n增广耦合变量的第k次迭代值；为区域系统n副本变量第k次迭代值；

步骤3.10，各区域系统根据步骤3.9得到的增广耦合变量近端值，计算各区域系统的增广耦合变量的迭代值；

各区域系统n∈{1,...,N}通过求解下式并行计算内部决策变量和增广耦合变量的第k+1次迭代值：

其中，为区域系统n内部决策变量的第k+1次迭代值，为区域系统n增广耦合变量的第k+1次迭代值；x_n为区域系统n的内部决策变量；为区域系统n的增广耦合变量；Ψ_n为区域系统n决策变量的可行域；为区域系统n缩放系数；f_n为区域系统n的校正控制代价；为区域系统n增广耦合变量的第k次迭代值；为原始对偶拟牛顿算法中区域系统n原始海塞矩阵的第k次迭代值；为区域系统n增广耦合变量近端值的第k次迭代值。

步骤3.11，所有邻居区域系统交换由步骤3.10得到的增广耦合变量的迭代值，各区域系统根据交换得到的邻居区域系统的增广耦合变量的迭代值计算各区域系统的副本变量的迭代值；

各区域系统n∈{1,...,N}在邻居之间交换区域系统n增广耦合变量的迭代值各区域系统n∈{1,...,N}根据交换得到的邻居区域系统的增广耦合变量的次迭代值，利用下式计算区域系统n副本变量的第k+1次迭代值：

式中，为为区域系统n副本变量第k+1次迭代值；Δ_n为区域系统n邻居的集合；ξ_n,m为区域系统n与区域系统m的关联矩阵；为区域系统n的邻居区域系统m的增广耦合变量的第k+1次迭代值。

步骤3.12，根据步骤3.10得到的增广耦合变量和步骤3.11得到的副本变量，各区域系统计算增广耦合变量梯度值的迭代值；

各区域系统n∈{1,...,N}通过下式计算区域系统n增广耦合变量梯度值的第k次迭代值

式中，为区域系统n增广耦合变量梯度值的第k次迭代值；为区域系统n的增广耦合变量的第k+1次迭代值；为区域系统n副本变量第k+1次迭代值。

步骤3.13，各区域系统根据步骤3.4得到的拉格朗日乘子的迭代值和步骤3.12得到的拉格朗日乘子梯度值的迭代值，计算各区域系统的拉格朗日乘子变化量的迭代值和各区域系统的拉格朗日乘子梯度变化量的迭代值；

各区域系统n∈{1,...,N}根据下式计算区域系统n的拉格朗日乘子变化量的第k次迭代值和区域系统n的拉格朗日乘子梯度变化量的第k次迭代值：

式中，区域系统n的拉格朗日乘子变化量的第k次迭代值；D_n为区域系统n的正规化矩阵；为区域系统n拉格朗日乘子的第k次迭代值；为区域系统n拉格朗日乘子的第k-1次迭代值；为区域系统n的拉格朗日乘子梯度变化量的第k次迭代值；为区域系统n拉格朗日乘子梯度值的第k次迭代值；为区域系统n拉格朗日乘子梯度值的第k-1次迭代值；γ为正则化参数，一般设置为较小的正数。

步骤3.14，各区域系统根据步骤3.13得到的各区域系统的拉格朗日乘子的变化量的迭代值和各区域系统的拉格朗日乘子的梯度变化量的迭代值，计算各区域系统拉格朗日乘子的海塞矩阵的迭代值；

各区域系统n∈{1,...,N}基于如下D-BFGS迭代计算区域系统n拉格朗日乘子的海塞矩阵的第k次迭代值若则：

若则：

式中，为区域系统n拉格朗日乘子的海塞矩阵的第k次迭代值；为区域系统n拉格朗日乘子的海塞矩阵的第k-1次迭代值；为区域系统n的拉格朗日乘子梯度变化量的第k次迭代值；区域系统n的拉格朗日乘子变化量的第k次迭代值；I为单位矩阵；γ为正则化参数，一般设置为较小的正数。

步骤3.15，各区域系统根据步骤3.14得到的各区域系统的拉格朗日乘子的海塞矩阵的迭代值，计算各区域系统拉格朗日乘子的更新方向本地值的迭代值；

各区域系统n∈{1,...,N}根据下式并行计算区域系统n拉格朗日乘子更新方向的本地值的第k次迭代值：

式中，为区域系统n拉格朗日乘子更新方向的本地值的第k次迭代值；为区域系统n拉格朗日乘子的海塞矩阵的第k次迭代值；D_n为区域系统n的正规化矩阵；Υ为正则化参数，一般设置为较小的正数；为区域系统n的增广耦合变量的第k+1次迭代值；为区域系统n副本变量第k+1次迭代值。

步骤3.16，所有邻居区域系统之间交换由步骤3.13得到的拉格朗日乘子的更新方向本地值的迭代值，各区域系统根据交换得到的邻居区域系统的拉格朗日乘子的迭代值的更新方向计算各区域系统的拉格朗日乘子的迭代值；

所有邻居区域系统之间交换各区域系统n∈{1,…,N}根据各自邻居的响应信息计算区域系统n拉格朗日乘子的第k+1次迭代值：

式中，为区域系统n拉格朗日乘子的第k+1次迭代值；为区域系统n拉格朗日乘子的第k次迭代值；Δ_n为区域系统n邻居的集合；κ_n为区域系统n的拉格朗日乘子的迭代步长；为区域系统n的邻居区域系统i的拉格朗日乘子更新方向的本地值的第k次迭代值；。

步骤3.17，各区域系统根据收敛条件判断算法收敛状态，若算法收敛则停止迭代输出各区域系统的校正控制策略，若算法未收敛则进入下一轮迭代；

根据步骤3.10得到的增广耦合变量和步骤3.11得到的副本变量，各区域系统n∈{1,...,N}通过收敛条件判断算法的收敛状态。其中，收敛条件可设置为但不仅限于下式形式：

式中，N为区域系统的总数；为区域系统n的增广耦合变量的第k+1次迭代值；为区域系统n副本变量第k+1次迭代值；ε为原始对偶拟牛顿算法收敛精度门槛值。

判断是否满足收敛条件，若满足收敛条件，则停止算法，返回各区域系统n∈{1,...,N}的最优值f_n(P_n)及最优解并输出最有校正控制策略。否则，设置迭代次数k＝k+1，返回步骤3.5，进入下一轮迭代。

本发明优点在于，对天然气系统动态管流的建模方法兼顾计算精度和效率；分布式求解协同校正控制策略，有助于维护运营独立性、缓解通信阻塞；具有有限时间收敛性，可以快速获得高精度的协同校正控制策略，具有在线应用潜力。

附图说明

图1：本发明流程图

图2：改造的电力-天然气互联系统算例接线图

具体实施方式

为了便于本领域的普通技术人员理解和实施本发明，下面结合附图及具体实施例对本发明做进一步说明，应当理解，此处所描述的实施示例仅用于解释和说明本发明，并不用于限定本发明。

如图1所示，一种电力-天然气互联系统分布式协同校正控制方法，其特征在于，包括以下步骤：

下面以基于IEEE118节点电力系统和比利时20节点天然气系统改造的电力-天然气互联系统为例，对本发明的技术方案做进一步的具体说明。改造的电力-天然气互联系统包括送端电力系统、受端电力系统、送端天然气系统、受端天然气系统，共四个区域系统，具体接线图如图2所示。

以送端电力系统与受端电力系统间的电力联络线65-68和47-69发生故障跳闸为例进行说明。设算例中电力系统基准值为100MVA，上述故障造成了送端电力系统和受端电力系统电力系统支路有功功率潮流越限共计0.9563p.u.。

步骤1：通过以电力-天然气互联系统总体校正控制代价最小化的目标函数构建电力-天然气互联系统协同校正控制优化模型；以电力系统的常规发电机功率调整量、燃气发电机功率调整量、电负荷切除量以及天然气系统气源调整量和气负荷切除量为电力-天然气互联系统协同校正控制优化模型的连续性优化变量；通过电力-天然气互联系统协同校正控制优化模型的电力系统运行约束条件、电力-天然气互联系统协同校正控制优化模型的电力系统安全限值约束条件、电力-天然气互联系统协同校正控制优化模型的通过正交配置法离散的天然气系统运行约束条件、电力-天然气互联系统协同校正控制优化模型的天然气系统安全限值约束条件、电力-天然气互联系统协同校正控制优化模型的天然气系统边界约束条件、电力-天然气互联系统协同校正控制优化模型的电力-天然气耦合约束条件构建电力-天然气互联系统协同校正控制优化模型的约束条件，进一步构建电力-天然气互联系统协同校正控制优化模型；

步骤1所述以电力-天然气互联系统总体校正控制代价最小化的目标函数，具体为：

设电力-天然气互联系统总体校正控制代价为C，可表示为2个区域电力系统控制代价C_e,i(i＝1,2)和2个区域天然气系统控制代价C_g,i(i＝1,2)之和：

式中，P_L0,n为校正前区域电力系统n的电力负荷有功功率向量，P_L,n为区域电力系统n的切负荷决策变量，M_E,n＝1000为区域电力系统n切除电负荷的成本系数；P_T0,n为校正前区域电力系统n的燃气发电机有功出力向量，P_T,n为区域电力系统n的燃气发电机调整决策变量，E_GT,n＝5为区域电力系统n燃气发电机调整的成本系数；P_G0,n为校正前区域电力系统n的常规发电机有功出力向量，P_G,n为区域电力系统n的常规发电机调整决策变量，E_G,n＝10为区域电力系统i常规发电机调整的成本系数；N＝2表示区域电力系统的数量；f_S0,m为校正前区域天然气系统m气源的质量流率向量，为区域天然气系统m气源的质量流率决策变量，E_S,m＝0.4区域天然气系统m气源调整的成本系数；f_L0,m为校正前区域天然气系统m气负荷的质量流率向量，为区域天然气系统m常规气负荷的质量流率决策变量，M_G,m＝100区域天然气系统m切除常规气负荷的成本系数；f_T0,m为校正前区域天然气系统m燃气发电机气负荷的质量流率向量，为区域天然气系统m燃气发电机气负荷的质量流率决策变量，E_T,m＝1区域天然气系统m燃气发电机气负荷调整的成本系数；M＝2表示区域天然气系统的数量。

步骤1所述电力-天然气互联系统协同校正控制优化模型的电力系统运行约束条件，具体为：

设Λ为电力-天然气互联系统中电力系统节点的集合，且电力系统满足二阶锥交流潮流模型：

式中，P_Gi为节点i注入的常规发电机有功功率；Q_Gi为节点i注入的常规发电机无功功率；P_Ti为节点i注入的燃气发电机有功功率；Q_Ti为节点i注入的燃气发电机无功功率；P_Li为节点i流出的负荷有功功率；Q_Li为节点i流出的负荷无功功率；δ(i)为和节点i直接连接的所有节点的集合；P_ij为从节点i流向节点j的有功功率；Q_ij为从节点i流向节点j的无功功率；P_ji为从节点j流向节点i的有功功率；Q_ji为从节点j流向节点i的无功功率；g_ij为线路i-j的电导；b_ij为线路i-j的电纳；b_i,ij为线路i-j的i侧对地电纳；b_j,ij为线路i-j的j侧的对地电纳；、V_i ^(c)为节点i的电压典型值；为节点i的电压典型值；辅助变量W_i＝V_i ²；辅助变量W_j＝V_j ²；辅助变量W_cij＝V_iV_j cosθ_ij；辅助变量W_sij＝V_iV_j sinθ_ij；

步骤1所述电力-天然气互联系统协同校正控制优化模型的电力系统安全限值约束条件，具体为：

式中，Λ_G为常规发电机节点的集合，Λ_T为燃气发电机节点的集合，Λ_L为电负荷节点的集合，Λ为电力系统所有节点集合；为常规发电机i有功出力的上限；P _Gi为常规发电机i有功出力的下限；为常规发电机i无功出力的上限；Q _Gi为常规发电机i无功出力的下限；为常规发电机i爬坡率的上限；R _G,i＝-15％为常规发电机i爬坡率的下限；为燃气发电机i有功出力的上限；P _Ti为燃气发电机i有功出力的下限；为燃气发电机i无功出力的上限；Q _Ti为燃气发电机i无功出力的下限；为燃气发电机i爬坡率的上限；R _T,i＝-30％为燃气发电机i爬坡率的下限；为节点i有功负荷切除调整的上限，也即节点i故障前的有功负荷大小；P _Li为节点i有功负荷切除调整的下限；为节点i无功负荷切除调整的上限，也即节点i故障前的无功负荷大小；Q _Li为节点i无功负荷切除调整的下限；为系统有功网损的上界；为节点i电压幅值的上限；V _i为节点i电压幅值的下限。

步骤1所述电力-天然气互联系统协同校正控制优化模型的通过正交配置法离散的天然气系统运行约束条件，具体为：

设Υ为天然气系统中节点的集合，天然气系统满足通过正交配置法离散的动态管流模型：

式中，Δt＝3600s为天然气系统待优化的时间长度；Δx_ab为输气管道ab的长度；A_ab为输气管道ab的横截面积；D_ab为输气管道ab的管壁内径；N_t＝6为时间正交配置点的个数；N_x＝14为空间正交配置点的个数；为空间一阶微分矩阵的第(i，k)个元素，为时间一阶微分矩阵的第(j，k)个元素；f_ab,ij为输气管道ab的质量流率在第(i，j)个时空配置点处的取值；p_ab,ij为输气管道ab的气体压力在第(i，j)个时空配置点处的取值；λ＝0.01为管壁摩擦系数；为输气管道ab的平均气体流速；v_s为天然气的气体音速；f_Sa,j为天然气系统节点a气源的质量流率在第j个时间配置点处的取值；f_La,j天然气系统节点a常规气负荷的质量流率在第j个时间配置点处的取值；f_Ta,j天然气系统节点a燃气发电机气负荷的质量流率在第j个时间配置点处的取值；为压缩机ba出口处质量流率在第j个时间配置点处的取值；f_Cab,0j压缩机ab入口处质量流率在第j个时间配置点处的取值。

步骤1所述电力-天然气互联系统协同校正控制优化模型的天然气系统安全限值约束条件，具体为：

式中，i＝1,…,N_x，N_x＝14为空间域正交配置插值点的个数；j＝1,…,N_t，N_t＝6为时域正交配置插值点的个数；为管道ab气体压力的上限；p _ab为管道ab气体压力的下限；为管道ab可输送气体质量流率的上限；f _ab为管道ab可输送气体质量流率的下限；为节点a气源流量的上限；f _Sa为节点a气源流量的下限；为校正控制前节点a的气负荷质量流率；f _La为节点a气负荷质量流率的下限；为节点a燃气机组气负荷质量流率可调节的上限；f _Ta分别为节点a燃气机组气负荷质量流率可调节的下限。

步骤1所述电力-天然气互联系统协同校正控制优化模型的天然气系统边界约束条件，具体为：

其中，ρ_ab,ij表示管道ab中气体密度在第(i，j)个时空配置点处的取值；和分别表示管道ab在初始时刻第i个空间配置点处的取值；f_Sa,j为天然气系统节点a气源的质量流率在第j个时间配置点处的取值；f_La,j天然气系统节点a常规气负荷的质量流率在第j个时间配置点处的取值；f_Ta,j天然气系统节点a燃气发电机气负荷的质量流率在第j个时间配置点处的取值；为初始时刻气节点a气源的质量流率；为初始时刻节点a气负荷的质量流率；为初始时刻节点a燃气发电机气负荷的质量流率；校正控制策略中节点a气源质量流率的给定值；为校正控制策略中节点a气负荷质量流率的给定值；为校正控制策略中节点a燃气机组进气量质量流率的给定值，N_t＝6为时域正交配置插值点的个数；N_x＝14为空间域正交配置插值点的个数。

步骤1所述电力-天然气互联系统协同校正控制优化模型的电力-天然气耦合约束条件，具体为：

设电力节点i∈Λ的燃气发电机是天然气节点m∈Υ的气负荷，则电力节点i与天然气节点m的耦合约束条件为：

式中，为校正控制策略中天然气节点a燃气机组进气量质量流率的给定值；P_Ti为电力节点i注入的燃气发电机有功功率；η_Ti,a＝45％为燃气发电机的工作效率，H_G＝37.62为天然气热值。

设压缩机ab，a,b∈Υ是电力节点i∈Λ的电负荷，则其耦合约束条件为：

式中，P_Ci为电网节点i上电压缩机消耗的电功率；η_Ci,ab压缩机ab消耗电功率与其原动机驱动流量间的比例系数；为压缩机ab出口处压力在时间配置点l处的值；p_ab,0l为为压缩机ab入口处压力在时间配置点l处的值；K_Cab为压缩机ab的增压比；f _Cab,kl为压缩机ab在第(k，l)个时空配置点处质量流率的值；f _Cab为压缩机ab质量流率的下限；为压缩机ab质量流率的上限；f_Cab,k0为压缩机ab的空间配置点k处质量流率在初始时间配置点处的值；为压缩机ab的空间配置点k处的质量流率在初始时刻值。

步骤1所述电力-天然气互联系统协同校正控制优化模型，决策变量包括所有与约束条件和目标函数有关的待确定的量。电力系统中的决策变量具体为：电力节点i注入的常规发电有功功率P_Gi；电力节点i注入的常规发电机无功功率Q_Gi；电力节点i注入的燃气发电有功功率P_Ti；电力节点i注入的燃气发电机无功功率Q_Ti；电力节点i流出的负荷有功功率P_Li；电力节点i流出的负荷无功功率Q_Li；从电力节点i流向电力节点j的有功功率P_ij；从电力节点i流向电力节点j的无功功率Q_ij；辅助变量W_i；辅助变量W_cij；辅助变量W_sij。天然气系统中的决策变量具体为：输气管道ab的质量流率在第(i，j)个时空配置点处的取值f_ab,ij；输气管道ab的气体压力在第(i，j)个时空配置点处的取值p_ab,ij；天然气系统节点a气源的质量流率在第j个时间配置点处的取值f_Sa,j；天然气系统节点a常规气负荷的质量流率在第j个时间配置点处的取值f_La,j；天然气系统节点a燃气发电机气负荷的质量流率在第j个时间配置点处的取值f_Ta,j；压缩机ba出口处质量流率在第j个时间配置点处的取值压缩机ab入口处质量流率在第j个时间配置点处的取值f_Cab,0j。

步骤2：将电力-天然气互联系统协同校正控制优化模型按照电力-天然气互联系统的地理区划通过节点复制法解耦得到多个区域系统校正控制子模型；在得到的区域系统校正控制子模型中将待协调的耦合变量通过增广处理统一向量维度得到增广耦合变量；将每个区域系统校正控制子模型和每个区域系统对应的增广耦合变量，通过辅助问题原则法构建区域系统辅助优化子模型。

步骤2所述将电力-天然气互联系统协同校正控制优化模型按照电力-天然气互联系统的地理区划通过节点复制法解耦得到多个区域系统校正控制子模型，具体为：

定义区域电力系统为区域内部电力节点和输电线路的集合，区域天然气系统为区域内部天然气节点和输气管道的集合。不同的区域电力系统之间通过跨区电力联络线连接，不同的区域天然气系统之间通过跨区输气管道连接，区域电力系统与区域天然气系统之间通过燃气发电机和电力驱动的加压站耦合。区域电力系统和区域天然气系统统称为区域系统。

以两个通过电力联络线连接的区域电力系统的解耦为例，对节点复制法和所得到的区域系统校正控制子模型进行说明。

以联络线15-33为例，其中节点15位于送端电力系统，节点33位于受端电力系统。基于步骤1中所定义的电力-天然气互联系统协同校正控制优化模型，对送端电力系统而言联络线15-33涉及的决策变量包括y₁₅＝[W₁₅,θ₁₅,P_15-33,Q_15-33,W_c15-33,W_s15-33]，对受端电力系统而言联络线15-33涉及决策变量包括y_l＝[W₃₃,θ₃₃,P_33-15,Q_33-15,W_c33-15,W_s33-15]；由于节点33的决策变量无法直接参与送端电力系统的优化，节点15上的决策变量无法直接参与受端电力系统的优化，为了将步骤1所述电力-天然气互联系统协同校正控制优化模型解耦，在送端电力系统中引入虚拟节点33’，受端电力系统n中引入节点15’。

其中送端电力系统中虚拟节点33’上的决策变量为y_33′＝[W_33′,θ_33′,P_33′,Q_33′,W_c15-33′,W_s15-33′]，受端电力系统中虚拟节点15’上的决策变量为y_15′＝[W_15′,θ_15′,P_15′,Q_15′,W_c15′-33,W_s15′-33]；

所述y₁₅与所述y_15′、所述y₃₃与所述y_33′互为副本变量；

则对于电力联络线15-33而言，送端电力系统和受端电力系统之间的耦合变量分别定义为v_送端＝[y₁₅,y_33′]和v_受端＝[y_15′,y₃₃]；

在引入虚拟节点后，原联络线15-33被替换为虚拟输电线路15-33’和虚拟输电线路15’-33。其中，(虚拟)输电线路15-33’参与送端电力系统的优化，虚拟输电线路15’-33参与受端电力系统的优化。

为了使送端电力系统和受端电力系统分别优化得到的虚拟输电线路15’-33和15-33’的决策变量结果对于实际联络线15-33可行，必须使虚拟输电线路15’-33和15-33’上对应的决策变量相等，也即需要满足一致性约束条件。具体而言，在联络线15-33解耦之后送端电力系统和受端电力系统的区域系统校正控制子模型分别为：

对于送端电力系统，在对送端电力系统内部的节点和输电线路相关决策变量构建如步骤1所述电力-天然气互联系统协同校正控制优化模型的基础上，还需满足一致性约束条件：

其中，为送端电力系统通过与受端电力系统通信获得的y₃₃实际值。

对于受端电力系统，在对受端电力系统内部的节点和输电线路相关决策变量构建如步骤1所述电力-天然气互联系统协同校正控制优化模型的基础上，还需满足一致性约束条件：

其中，为受端电力系统通过与送端电力系统通信获得的y₁₅实际值。

类似地，电力-天然气互联系统中其他的电力联络线、跨区输气管道和电力-天然气耦合设备均可基于上述节点复制法引入虚拟节点、副本变量、耦合变量和一致性约束条件，实现电力-天然气互联系统协同校正控制优化模型的解耦，得到多个区域系统校正控制子模型。

步骤2所述在得到的多个区域系统校正控制子模型中将待协调的耦合变量通过增广处理统一向量维度得到增广耦合变量，具体为：

对各个区域系统耦合变量v_n,n＝1,...,4插入0元素构建增广耦合变量使得各个区域系统中待协调的耦合变量向量维度相同；

假设有M个耦合变量需要协调求解，则得到的增广耦合变量

步骤2所述将每个区域系统校正控制子模型和每个区域系统对应的增广耦合变量，通过辅助问题原则法构建区域系统辅助优化子模型，具体为：

所述区域系统辅助优化子模型包含区域系统原始辅助子优化模型和区域系统对偶辅助子优化模型；

所述区域系统n,n＝1,...,4的原始辅助子优化模型如下：

式中，所有变量和系数的上角标为迭代轮次标识；为根据步骤2.2得到的增广耦合变量；x_n为区域系统n的内部决策变量；构成区域系统n子优化模型的全部决策变量；Ψ_n为区域系统n决策变量的可行域，根据步骤1中定义的电力系统运行约束条件和安全限值约束条件确定；c为罚系数，一般取为较大的正数；为与区域系统n收敛性相关的缩放系数；为区域系统n增广耦合变量对应的海塞矩阵的估计矩阵；为区域系统n增广耦合变量对应的副本变量，在迭代过程中基于邻居电力系统m∈Δ_n的通信进行更新。其中ξ_n,m为区域系统n与区域系统m的关联矩阵，其定义如下：

区域系统n的对偶辅助子优化问题如下：

式中，所有变量和系数的上角标为迭代轮次标识；为区域系统n拉格朗日乘子；为区域系统n拉格朗日乘子对应的k轮迭代海塞矩阵估计值；为区域系统n增广耦合变量对应的副本变量。

步骤3：以故障后电力-天然气互联系统运行状态作为步骤2的区域系统辅助优化子模型中决策变量的初值，基于原始对偶拟牛顿算法分布式求解区域系统辅助优化子模型，得到电力-天然气互联系统的协同校正控制策略；

所述步骤3具体如下：

步骤3.1，所述初始化区域系统辅助优化子模型的决策变量，初始化的原始对偶拟牛顿算法；

各个区域系统n∈{1,2,3,4}根据故障后电力-天然气互联系统运行状态设置区域系统n内部决策变量初值区域系统n增广耦合变量初值以及区域系统n副本变量初值原始对偶拟牛顿算法的参数初始化具体为：初始化区域系统n拉格朗日乘子的初值初始化原始对偶拟牛顿算法中区域系统n原始海塞矩阵初始化区域系统n增广耦合变量变化量初始化区域系统n增广耦合变量梯度值初始化区域系统n增广耦合变量梯度变化量初始化区域系统n拉格朗日乘子的海塞矩阵初始化区域系统n拉格朗日乘子的变化量初始化区域系统n拉格朗日乘子的梯度值初始化区域系统n拉格朗日乘子的梯度变化量初始化罚系数c＝1500；设置迭代次数k＝0；设置原始对偶拟牛顿算法收敛精度门槛值ε＝5×10^-6。

步骤3.2，利用初始化的原始对偶拟牛顿算法求解初始化的原始辅助子优化模型，得到各区域系统内部决策变量和增广耦合变量的迭代值；

设置迭代次数k＝k+1，各区域系统n∈{1,2,3,4}并行计算如下初始化的原始辅助子优化模型，得到各区域系统内部决策变量和增广耦合变量的第k次迭代值：

其中，为区域系统n内部决策变量的第k次迭代值；为区域系统n增广耦合变量的第k次迭代值；x_n为区域系统n的内部决策变量；为区域系统n的增广耦合变量；Ψ_n为区域系统n决策变量的可行域；为区域系统n增广耦合变量初值；为原始对偶拟牛顿算法中区域系统n原始海塞矩阵初值；为区域系统n拉格朗日乘子的初值；f_n为区域系统n的校正控制代价。

步骤3.3，所有邻居区域系统之间交换增广耦合变量的迭代值，各区域系统根据交换得到的邻居区域系统耦合变量的迭代值计算各区域系统副本变量的迭代值；

所有邻居区域系统之间交换增广耦合变量各区域系统n∈{1,2,3,4}根据下式计算区域系统n副本变量的第k次迭代值：

式中，为区域系统n副本变量的第k次迭代值；Δ_n为区域系统n邻居的集合；ξ_n,m为区域系统n与区域系统m的关联矩阵；为邻居区域系统m的增广耦合变量的第k次迭代值。

步骤3.4，各区域系统副本变量的迭代值带入初始化的对偶辅助子优化模型，计算各区域系统的拉格朗日乘子的迭代值；

各区域系统n∈{1,2,3,4}并行计算如下对偶辅助子优化模型，得到区域系统n拉格朗日乘子的第k次迭代值：

式中，为区域系统n拉格朗日乘子的第k次迭代值；为区域系统n拉格朗日乘子的初值；为区域系统n拉格朗日乘子的海塞矩阵初值；为区域系统n增广耦合变量的第k次迭代值；为区域系统n副本变量第k次迭代值。

步骤3.5，各区域系统利用步骤3.4得到各区域系统的拉格朗日乘子的迭代值，计算增广耦合变量梯度值的迭代值；

各区域系统n∈{1,2,3,4}根据下式计算区域系统n增广耦合变量梯度值的第k次迭代值

式中，为增广耦合变量梯度值的第k次迭代值；c＝1500为罚系数；为区域系统n副本变量第k次迭代值；为区域系统n增广耦合变量的第k次迭代值。

步骤3.6，各区域系统利用各区域系统的增广耦合变量的迭代值和各区域系统的增广耦合变量梯度值的迭代值，计算各区域系统的增广耦合变量变化量的迭代值和增广耦合变量梯度变化量的迭代值；

各区域系统n∈{1,2,3,4}通过下式计算区域系统n增广耦合变量变化量的第k次迭代值和区域系统n增广耦合变量梯度变化量的第k次迭代值：

式中，为区域系统n增广耦合变量变化量的第k次迭代值；D_n为区域系统n的正规化矩阵；为区域系统n增广耦合变量的第k次迭代值；为区域系统n增广耦合变量的第k-1次迭代值；为区域系统n耦合变量梯度变化量的第k次迭代值；为区域系统n的增广耦合变量梯度值的第k次迭代值；为区域系统n的增广耦合变量梯度值的第k-1次迭代值；σ＝10^-2为正则化参数，一般设置为较小的正数。

步骤3.7，各区域系统利用步骤3.6得到的各区域系统的增广耦合变量变化量的迭代值和增广耦合变量梯度变化量的迭代值计算各区域系统的增广耦合变量原始海塞矩阵的迭代值；

各区域系统n∈{1,2,3,4}基于如下D-BFGS迭代计算区域系统n原始海塞矩阵的第k次迭代值。若则：

式中，为区域系统n原始海塞矩阵的第k次迭代值；为区域系统n原始海塞矩阵的第k-1次迭代值；为区域系统n耦合变量梯度变化量的第k次迭代值；为区域系统n增广耦合变量变化量的第k次迭代值；I为单位矩阵；σ＝10^-2为正则化参数。

若则：

式中，为区域系统n原始海塞矩阵的第k次迭代值；为区域系统n原始海塞矩阵的第k-1次迭代值。

步骤3.8，各区域系统基于步骤3.7得到的各区域系统的增广耦合变量海塞矩阵的迭代值，计算各区域系统的耦合变量更新方向本地值的迭代值；

各区域系统n∈{1,2,3,4}根据下式并行计算区域系统n耦合变量更新方向本地值的第k次迭代值：

式中，为区域系统n耦合变量更新方向本地值的第k次迭代值；为区域系统n原始海塞矩阵的第k次迭代值；Ψ＝10^-5为正则化参数；D_n为区域系统n的正规化矩阵；为区域系统n拉格朗日乘子的第k次迭代值；c＝1500为罚系数；为区域系统n增广耦合变量的第k次迭代值；为区域系统n副本变量第k次迭代值。

步骤3.9，所有邻居区域系统之间交换由步骤3.8得到的增广耦合变量的更新方向本地值的迭代值，各区域系统根据交换得到的邻居区域系统的增广耦合变量的更新方向本地值的迭代值计算各区域系统的增广耦合变量的近端值的迭代值；

所有邻居区域系统之间交换各区域系统n∈{1,2,3,4}由下式并行计算区域系统n增广耦合变量近端值的第k次迭代值

式中，为区域系统n增广耦合变量近端值的第k次迭代值；为区域系统n增广耦合变量的第k次迭代值；为区域系统n迭代步长；Δ_n为区域系统n邻居的集合；为区域系统n缩放系数；为区域系统n的邻居区域系统i耦合变量更新方向本地值的第k次迭代值。其中，迭代步长具有如下形式的闭式解：

式中，为区域系统n拉格朗日乘子的第k次迭代值；Δ_n为区域系统n邻居的集合；为区域系统n增广耦合变量的第k次迭代值；为区域系统n副本变量第k次迭代值；

步骤3.10，各区域系统根据步骤3.9得到的增广耦合变量近端值，计算各区域系统的增广耦合变量的迭代值；

各区域系统n∈{1,2,3,4}通过求解下式并行计算内部决策变量和增广耦合变量的第k+1次迭代值：

其中，为区域系统n内部决策变量的第k+1次迭代值，为区域系统n增广耦合变量的第k+1次迭代值；x_n为区域系统n的内部决策变量；为区域系统n的增广耦合变量；Ψ_n为区域系统n决策变量的可行域；为区域系统n缩放系数；f_n为区域系统n的校正控制代价；为区域系统n增广耦合变量的第k次迭代值；为原始对偶拟牛顿算法中区域系统n原始海塞矩阵的第k次迭代值；为区域系统n增广耦合变量近端值的第k次迭代值。

步骤3.11，所有邻居区域系统交换由步骤3.10得到的增广耦合变量的迭代值，各区域系统根据交换得到的邻居区域系统的增广耦合变量的迭代值计算各区域系统的副本变量的迭代值；

各区域系统n∈{1,2,3,4}在邻居之间交换区域系统n增广耦合变量的迭代值各区域系统n∈{1,2,3,4}根据交换得到的邻居区域系统的增广耦合变量的次迭代值，利用下式计算区域系统n副本变量的第k+1次迭代值：

式中，为为区域系统n副本变量第k+1次迭代值；Δ_n为区域系统n邻居的集合；ξ_n,m为区域系统n与区域系统m的关联矩阵；为区域系统n的邻居区域系统m的增广耦合变量的第k+1次迭代值。

步骤3.12，根据步骤3.10得到的增广耦合变量和步骤3.11得到的副本变量，各区域系统计算增广耦合变量梯度值的迭代值；

各区域系统n∈{1,2,3,4}通过下式计算区域系统n增广耦合变量梯度值的第k次迭代值

式中，为区域系统n增广耦合变量梯度值的第k次迭代值；为区域系统n的增广耦合变量的第k+1次迭代值；为区域系统n副本变量第k+1次迭代值。

步骤3.13，各区域系统根据步骤3.4得到的拉格朗日乘子的迭代值和步骤3.12得到的拉格朗日乘子梯度值的迭代值，计算各区域系统的拉格朗日乘子变化量的迭代值和各区域系统的拉格朗日乘子梯度变化量的迭代值；

各区域系统n∈{1,2,3,4}根据下式计算区域系统n的拉格朗日乘子变化量的第k次迭代值和区域系统n的拉格朗日乘子梯度变化量的第k次迭代值：

式中，区域系统n的拉格朗日乘子变化量的第k次迭代值；D_n为区域系统n的正规化矩阵；为区域系统n拉格朗日乘子的第k次迭代值；为区域系统n拉格朗日乘子的第k-1次迭代值；为区域系统n的拉格朗日乘子梯度变化量的第k次迭代值；为区域系统n拉格朗日乘子梯度值的第k次迭代值；为区域系统n拉格朗日乘子梯度值的第k-1次迭代值；γ＝10^-2为正则化参数，一般设置为较小的正数。

步骤3.14，各区域系统根据步骤3.13得到的各区域系统的拉格朗日乘子的变化量的迭代值和各区域系统的拉格朗日乘子的梯度变化量的迭代值，计算各区域系统拉格朗日乘子的海塞矩阵的迭代值；

各区域系统n∈{1,2,3,4}基于如下D-BFGS迭代计算区域系统n拉格朗日乘子的海塞矩阵的第k次迭代值若则：

若则：

式中，为区域系统n拉格朗日乘子的海塞矩阵的第k次迭代值；为区域系统n拉格朗日乘子的海塞矩阵的第k-1次迭代值；为区域系统n的拉格朗日乘子梯度变化量的第k次迭代值；区域系统n的拉格朗日乘子变化量的第k次迭代值；I为单位矩阵；γ＝10^-2为正则化参数，一般设置为较小的正数。

步骤3.15，各区域系统根据步骤3.14得到的各区域系统的拉格朗日乘子的海塞矩阵的迭代值，计算各区域系统拉格朗日乘子的更新方向本地值的迭代值；

各区域系统n∈{1,2,3,4}根据下式并行计算区域系统n拉格朗日乘子更新方向的本地值的第k次迭代值：

式中，为区域系统n拉格朗日乘子更新方向的本地值的第k次迭代值；为区域系统n拉格朗日乘子的海塞矩阵的第k次迭代值；D_n为区域系统n的正规化矩阵；Υ＝10^-3为正则化参数；为区域系统n的增广耦合变量的第k+1次迭代值；为区域系统n副本变量第k+1次迭代值。

步骤3.16，所有邻居区域系统之间交换由步骤3.13得到的拉格朗日乘子的更新方向本地值的迭代值，各区域系统根据交换得到的邻居区域系统的拉格朗日乘子的迭代值的更新方向计算各区域系统的拉格朗日乘子的迭代值；

所有邻居区域系统之间交换各区域系统n∈{1,2,3,4}根据各自邻居的响应信息计算区域系统n拉格朗日乘子的第k+1次迭代值：

式中，为区域系统n拉格朗日乘子的第k+1次迭代值；为区域系统n拉格朗日乘子的第k次迭代值；Δ_n为区域系统n邻居的集合；κ_n为区域系统n的拉格朗日乘子的迭代步长；为区域系统n的邻居区域系统i的拉格朗日乘子更新方向的本地值的第k次迭代值；。

步骤3.17，各区域系统根据收敛条件判断算法收敛状态，若算法收敛则停止迭代输出各区域系统的校正控制策略，若算法未收敛则进入下一轮迭代；

根据步骤3.10得到的增广耦合变量和步骤3.11得到的副本变量，各区域系统n∈{1,2,3,4}通过收敛条件判断算法的收敛状态。其中，收敛条件可设置为但不仅限于下式形式：

式中，N＝4为区域系统的总数；为区域系统n的增广耦合变量的第k+1次迭代值；为区域系统n副本变量第k+1次迭代值；ε＝5×10^-6为原始对偶拟牛顿算法收敛精度门槛值。

判断是否满足收敛条件，若满足收敛条件，则停止算法，返回各区域系统n∈{1,2,3,4}的最优值f_n(P_n)及最优解并输出最有校正控制策略。否则，设置迭代次数k＝k+1，返回步骤3.5，进入下一轮迭代。

在步骤3中，原始对偶拟牛顿算法的各项参数设置如表2所示。

表2 步骤3中算法参数设置

步骤3经过12次迭代计算收敛，并输出协同校正控制策略如表3所示。

表3 协同校正控制策略结果

由表3可见，本发明所述方法可以通过少量的迭代计算快速收敛，并输出可行的协同校正控制策略。基于本发明得到的协同校正控制策略可以在合理的总校正控制代价下充分调动电力-天然气互联系统中的可控资源，消除电力联络线N-2故障所引发的电力系统潮流越限事故，同时保障送受端天然气系统安全运行。

在基于本发明所得的协同校正控制策略中，燃气发电机校正控制量及校正控制后进气管道状态如表4所示。

表4 燃气发电机校正控制量及校正控制后进气管道状态

由表4可见，由于本发明构建的电力-天然气互联系统协同校正控制优化模型充分考虑电力系统和天然气系统双侧的安全运行约束，因而基于本发明所得到的协同校正控制策略既充分利用燃气发电机响应速度快、爬坡率高的灵活性优势，同时在校正控制过程中燃气机组进气口管道的气体压力和质量流率均位于安全限值范围内。

应当理解的是，本说明书未详细阐述的部分均属于现有技术。

应当理解的是，上述针对实施例的描述较为详细，并不能因此而认为是对本发明专利保护范围的限制，本领域的普通技术人员在本发明的启示下，在不脱离本发明权利要求所保护的范围情况下，还可以做出替换或变形，均落入本发明的保护范围之内，本发明的请求保护范围应以所附权利要求为准。

Claims (3)

1.一种电力天然气互联系统分布式协同校正控制方法，其特征在于，包括以下步骤：

步骤1：通过以电力-天然气互联系统总体校正控制代价最小化的目标函数构建电力-天然气互联系统协同校正控制优化模型；以电力系统的常规发电机功率调整量、燃气发电机功率调整量、电负荷切除量以及天然气系统气源调整量和气负荷切除量为电力-天然气互联系统协同校正控制优化模型的连续性优化变量；通过电力-天然气互联系统协同校正控制优化模型的电力系统运行约束条件、电力-天然气互联系统协同校正控制优化模型的电力系统安全限值约束条件、电力-天然气互联系统协同校正控制优化模型的通过正交配置法离散的天然气系统运行约束条件、电力-天然气互联系统协同校正控制优化模型的天然气系统安全限值约束条件、电力-天然气互联系统协同校正控制优化模型的天然气系统边界约束条件、电力-天然气互联系统协同校正控制优化模型的电力-天然气耦合约束条件构建电力-天然气互联系统协同校正控制优化模型的约束条件，进一步构建电力-天然气互联系统协同校正控制优化模型；

步骤2：将电力-天然气互联系统协同校正控制优化模型按照电力-天然气互联系统的地理区划通过节点复制法解耦得到多个区域系统校正控制子模型；在得到的区域系统校正控制子模型中将待协调的耦合变量通过增广处理统一向量维度得到增广耦合变量；将每个区域系统校正控制子模型和每个区域系统对应的增广耦合变量，通过辅助问题原则法构建区域系统辅助优化子模型；

步骤3：以故障后电力-天然气互联系统运行状态作为步骤2的区域系统辅助优化子模型中决策变量的初值，基于原始对偶拟牛顿算法分布式求解区域系统辅助优化子模型，得到电力-天然气互联系统的协同校正控制策略；

步骤1所述以电力天然气互联系统总体校正控制代价最小化的目标函数，具体为：

设电力-天然气互联系统总体校正控制代价为C，可表示为N个区域电力系统控制代价C_e,i(i＝1,2,……,N)和M个区域天然气系统控制代价C_g,i(i＝1,2,……,M)之和：

式中，P_L0,n为校正前区域电力系统n的电力负荷有功功率向量，P_L,n为区域电力系统n的切负荷决策变量，M_E,n为区域电力系统n切除电负荷的成本系数；P_T0,n为校正前区域电力系统n的燃气发电机有功出力向量，P_T,n为区域电力系统n的燃气发电机调整决策变量，E_GT,n为区域电力系统n燃气发电机调整的成本系数；P_G0,n为校正前区域电力系统n的常规发电机有功出力向量，P_G,n为区域电力系统n的常规发电机调整决策变量，E_G,n为区域电力系统i常规发电机调整的成本系数；N表示区域电力系统的数量；f_S0,m为校正前区域天然气系统m气源的质量流率向量，为区域天然气系统m气源的质量流率决策变量，E_S,m区域天然气系统m气源调整的成本系数；f_L0,m为校正前区域天然气系统m气负荷的质量流率向量，为区域天然气系统m常规气负荷的质量流率决策变量，M_G,m区域天然气系统m切除常规气负荷的成本系数；f_T0,m为校正前区域天然气系统m燃气发电机气负荷的质量流率向量，为区域天然气系统m燃气发电机气负荷的质量流率决策变量，E_T,m区域天然气系统m燃气发电机气负荷调整的成本系数；M表示区域天然气系统的数量；

步骤1所述电力-天然气互联系统协同校正控制优化模型的电力系统运行约束条件，具体为：

设Λ为电力-天然气互联系统中电力系统节点的集合，且电力系统满足二阶锥交流潮流模型：

式中，P_Gi为节点i注入的常规发电机有功功率；Q_Gi为节点i注入的常规发电机无功功率；P_Ti为节点i注入的燃气发电机有功功率；Q_Ti为节点i注入的燃气发电机无功功率；P_Li为节点i流出的负荷有功功率；Q_Li为节点i流出的负荷无功功率；δ(i)为和节点i直接连接的所有节点的集合；P_ij为从节点i流向节点j的有功功率；Q_ij为从节点i流向节点j的无功功率；P_ji为从节点j流向节点i的有功功率；Q_ji为从节点j流向节点i的无功功率；g_ij为线路i-j的电导；b_ij为线路i-j的电纳；b_i,ij为线路i-j的i侧对地电纳；b_j,ij为线路i-j的j侧的对地电纳；、V_i ^(c)为节点i的电压典型值；V_j ^(c)为节点i的电压典型值；辅助变量W_i＝V_i ²；辅助变量W_j＝V_j ²；辅助变量W_cij＝V_iV_jcosθ_ij；辅助变量W_sij＝V_iV_jsinθ_ij；

步骤1所述电力-天然气互联系统协同校正控制优化模型的电力系统安全限值约束条件，具体为：

式中，Λ_G为常规发电机节点的集合，Λ_T为燃气发电机节点的集合，Λ_L为电负荷节点的集合，Λ为电力系统所有节点集合；为常规发电机i有功出力的上限；P _Gi为常规发电机i有功出力的下限；为常规发电机i无功出力的上限；Q _Gi为常规发电机i无功出力的下限；为常规发电机i爬坡率的上限；R _G,i为常规发电机i爬坡率的下限；为燃气发电机i有功出力的上限；P _Ti为燃气发电机i有功出力的下限；为燃气发电机i无功出力的上限；Q _Ti为燃气发电机i无功出力的下限；为燃气发电机i爬坡率的上限；R _T,i为燃气发电机i爬坡率的下限；为节点i有功负荷切除调整的上限，也即节点i故障前的有功负荷大小；P _Li为节点i有功负荷切除调整的下限；为节点i无功负荷切除调整的上限，也即节点i故障前的无功负荷大小；Q _Li为节点i无功负荷切除调整的下限；为系统有功网损的上界；为节点i电压幅值的上限；V _i为节点i电压幅值的下限；

步骤1所述电力-天然气互联系统协同校正控制优化模型的通过正交配置法离散的天然气系统运行约束条件，具体为：

设Υ为天然气系统中节点的集合，天然气系统满足通过正交配置法离散的动态管流模型：

式中，Δt为天然气系统待优化的时间长度；Δx_ab为输气管道ab的长度；A_ab为输气管道ab的横截面积；D_ab为输气管道ab的管壁内径；N_t为时间正交配置点的个数；N_x为空间正交配置点的个数；为空间一阶微分矩阵的第(i，k)个元素，为时间一阶微分矩阵的第(j，k)个元素；f_ab,ij为输气管道ab的质量流率在第(i，j)个时空配置点处的取值；p_ab,ij为输气管道ab的气体压力在第(i，j)个时空配置点处的取值；λ为管壁摩擦系数；为输气管道ab的平均气体流速；v_s为天然气的气体音速；f_Sa,j为天然气系统节点a气源的质量流率在第j个时间配置点处的取值；f_La,j天然气系统节点a常规气负荷的质量流率在第j个时间配置点处的取值；f_Ta,j天然气系统节点a燃气发电机气负荷的质量流率在第j个时间配置点处的取值；为压缩机ba出口处质量流率在第j个时间配置点处的取值；f_Cab,0j压缩机ab入口处质量流率在第j个时间配置点处的取值；

步骤1所述电力-天然气互联系统协同校正控制优化模型的天然气系统安全限值约束条件，具体为：

式中，i＝1,…,N_x，为空间域正交配置插值点的个数；j＝1,…,N_t，N_t为时域正交配置插值点的个数；为管道ab气体压力的上限；p _ab为管道ab气体压力的下限；为管道ab可输送气体质量流率的上限；f _ab为管道ab可输送气体质量流率的下限；为节点a气源流量的上限；f _Sa为节点a气源流量的下限；为校正控制前节点a的气负荷质量流率；f _La为节点a气负荷质量流率的下限；为节点a燃气机组气负荷质量流率可调节的上限；f _Ta分别为节点a燃气机组气负荷质量流率可调节的下限；

步骤1所述电力-天然气互联系统协同校正控制优化模型的天然气系统边界约束条件，具体为：

其中，ρ_ab,ij表示管道ab中气体密度在第(i，j)个时空配置点处的取值；和分别表示管道ab在初始时刻第i个空间配置点处的取值；f_Sa,j为天然气系统节点a气源的质量流率在第j个时间配置点处的取值；f_La,j天然气系统节点a常规气负荷的质量流率在第j个时间配置点处的取值；f_Ta,j天然气系统节点a燃气发电机气负荷的质量流率在第j个时间配置点处的取值；为初始时刻气节点a气源的质量流率；为初始时刻节点a气负荷的质量流率；为初始时刻节点a燃气发电机气负荷的质量流率；校正控制策略中节点a气源质量流率的给定值；为校正控制策略中节点a气负荷质量流率的给定值；为校正控制策略中节点a燃气机组进气量质量流率的给定值；

步骤1所述电力-天然气互联系统协同校正控制优化模型的电力-天然气耦合约束条件，具体为：

设电力节点i∈Λ的燃气发电机是天然气节点m∈Υ的气负荷，则电力节点i与天然气节点m的耦合约束条件为：

式中，为校正控制策略中天然气节点a燃气机组进气量质量流率的给定值；P_Ti为电力节点i注入的燃气发电机有功功率；η_Ti,a为燃气发电机的工作效率，H_G为天然气热值；

设压缩机ab，a,b∈Υ是电力节点i∈Λ的电负荷，则其耦合约束条件为：

式中，P_Ci为电网节点i上电压缩机消耗的电功率；η_Ci,ab压缩机ab消耗电功率与其原动机驱动流量间的比例系数；为压缩机ab出口处压力在时间配置点l处的值；p_ab,0l为为压缩机ab入口处压力在时间配置点l处的值；K_Cab为压缩机ab的增压比；f_Cab,kl为压缩机ab在第(k，l)个时空配置点处质量流率的值；f _Cab为压缩机ab质量流率的下限；为压缩机ab质量流率的上限；f_Cab,k0为压缩机ab的空间配置点k处质量流率在初始时间配置点处的值；为压缩机ab的空间配置点k处的质量流率在初始时刻值；

步骤1所述电力-天然气互联系统协同校正控制优化模型，决策变量包括所有与约束条件和目标函数有关的待确定的量；电力系统中的决策变量具体为：电力节点i注入的常规发电有功功率P_Gi；电力节点i注入的常规发电机无功功率Q_Gi；电力节点i注入的燃气发电有功功率P_Ti；电力节点i注入的燃气发电机无功功率Q_Ti；电力节点i流出的负荷有功功率P_Li；电力节点i流出的负荷无功功率Q_Li；从电力节点i流向电力节点j的有功功率P_ij；从电力节点i流向电力节点j的无功功率Q_ij；辅助变量W_i；辅助变量W_cij；辅助变量W_sij；天然气系统中的决策变量具体为：输气管道ab的质量流率在第(i，j)个时空配置点处的取值f_ab,ij；输气管道ab的气体压力在第(i，j)个时空配置点处的取值p_ab,ij；天然气系统节点a气源的质量流率在第j个时间配置点处的取值f_Sa,j；天然气系统节点a常规气负荷的质量流率在第j个时间配置点处的取值f_La,j；天然气系统节点a燃气发电机气负荷的质量流率在第j个时间配置点处的取值f_Ta,j；压缩机ba出口处质量流率在第j个时间配置点处的取值压缩机ab入口处质量流率在第j个时间配置点处的取值f_Cab,0j。

2.根据权利要求1所述的电力天然气互联系统分布式协同校正控制方法，其特征在于，步骤2中，所述具体方法如下：

步骤2所述将电力-天然气互联系统协同校正控制优化模型按照电力-天然气互联系统的地理区划通过节点复制法解耦得到多个区域系统校正控制子模型，具体为：

定义区域电力系统为区域内部电力节点和输电线路的集合，区域天然气系统为区域内部天然气节点和输气管道的集合；不同的区域电力系统之间通过跨区电力联络线连接，不同的区域天然气系统之间通过跨区输气管道连接，区域电力系统与区域天然气系统之间通过燃气发电机和电力驱动的加压站耦合；区域电力系统和区域天然气系统统称为区域系统；

以两个通过电力联络线连接的区域电力系统的解耦为例，对节点复制法和所得到的区域系统校正控制子模型进行说明；假设区域电力系统m与区域电力系统n通过联络线k-l连接，其中节点k位于区域电力系统m，节点l位于区域电力系统n；

基于步骤1中所定义的集中优化模型，对区域电力系统m而言联络线k-l涉及的决策变量包括y_k＝[W_k,θ_k,P_lk,Q_lk,W_ckl,W_skl]，对区域电力系统n而言联络线k-l涉及决策变量包括y_l＝[W_l,θ_l,P_lk,Q_lk,W_clk,W_slk]；由于节点l的决策变量无法直接参与区域电力系统m的优化，节点k上的决策变量无法直接参与区域电力系统n的优化，为了将步骤1所述电力-天然气互联系统协同校正控制优化模型解耦，在区域电力系统m中引入虚拟节点l’，区域电力系统n中引入节点k’；

其中区域电力系统m中虚拟节点l’上的决策变量为y_l′＝[W_l′,θ_l′,P_l′,Q_l′,W_ckl′,W_skl′]，区域电力系统n中虚拟节点k’上的决策变量为y_k′＝[W_k′,θ_k′,P_k′,Q_k′,W_ck′l,W_sk′l]；

所述y_l与所述y_l′、所述y_k与所述y_k′互为副本变量；

则区域电力系统m与区域电力系统n之间耦合变量分别定义为v_m＝[y_k,y_l′]和v_n＝[y_k′,y_l]；

在引入虚拟节点后，原联络线k-l被替换为(虚拟)输电线路k-l’和(虚拟)输电线路k’-l；其中，(虚拟)输电线路k-l’参与区域电力系统m的优化，(虚拟)输电线路k’-l参与区域电力系统n的优化；

为了使区域电力系统m和区域电力系统n分别优化得到的(虚拟)输电线路k’-l和k-l’的决策变量结果对于实际联络线k-l可行，必须使(虚拟)输电线路k-l’和(虚拟)输电线路k’-l上对应的决策变量相等，也即需要满足一致性约束条件；具体而言，在解耦之后区域电力系统m和区域电力系统n的电力-天然气互联系统协同校正控制优化模型分别为：

对于区域电力系统m，在对区域电力系统m内部的节点和输电线路相关决策变量构建如步骤1所述电力-天然气互联系统协同校正控制优化模型的基础上，还需满足一致性约束条件：

其中，为区域电力系统m通过与区域电力系统n通信获得的y_l实际值；

对于区域电力系统n，在对区域电力系统n内部的节点和输电线路相关决策变量构建如步骤1所述电力-天然气互联系统协同校正控制优化模型的基础上，还需满足一致性约束条件：

其中，为区域电力系统n通过与区域电力系统m通信获得的实际值；

电力-天然气互联系统中所有电力联络线、输气管道和电力-天然气耦合设备均可基于上述节点复制法引入虚拟节点、副本变量、耦合变量和一致性约束条件，实现电力-天然气互联系统协同校正控制优化模型的解耦，得到多个区域系统校正控制子模型；

步骤2所述在得到的多个区域系统校正控制子模型中将待协调的耦合变量通过增广处理统一向量维度得到增广耦合变量，具体为：

对各个区域系统耦合变量v_n,n＝1,…,N插入0元素构建增广耦合变量使得各个区域系统中待协调的耦合变量向量维度相同；

假设共有M个耦合变量需要协调求解，则得到的增广耦合变量

步骤2所述将每个区域系统校正控制子模型和每个区域系统对应的增广耦合变量，通过辅助问题原则法构建区域系统辅助优化子模型，具体为：

所述区域系统辅助优化子模型包含区域系统原始辅助子优化模型和区域系统对偶辅助子优化模型；

所述区域系统n的原始辅助子优化模型如下：

式中，所有变量和系数的上角标为迭代轮次标识；为根据步骤2.2得到的增广耦合变量；x_n为区域系统n的内部决策变量；构成区域系统n子优化模型的全部决策变量；Ψ_n为区域系统n决策变量的可行域，根据步骤1中定义的电力系统运行约束条件和安全限值约束条件确定；c为罚系数，一般取为较大的正数；为与区域系统n收敛性相关的缩放系数；为区域系统n增广耦合变量对应的海塞矩阵的估计矩阵；为区域系统n增广耦合变量对应的副本变量，在迭代过程中基于邻居电力系统m∈Δ_n的通信进行更新；其中ξ_n,m为区域系统n与区域系统m的关联矩阵，其定义如下：

区域系统n的对偶辅助子优化问题如下：

式中，所有变量和系数的上角标为迭代轮次标识；为区域系统n拉格朗日乘子；为区域系统n拉格朗日乘子对应的k轮迭代海塞矩阵估计值；为区域系统n增广耦合变量对应的副本变量。

3.根据权利要求1所述的电力天然气互联系统分布式协同校正控制方法，其特征在于，所述步骤3具体如下：

步骤3.1，所述初始化区域系统辅助优化子模型的决策变量，初始化的原始对偶拟牛顿算法；

各个区域系统n∈{1,…,N}根据故障后电力-天然气互联系统运行状态设置区域系统n内部决策变量初值区域系统n增广耦合变量初值以及区域系统n副本变量初值原始对偶拟牛顿算法的参数初始化具体为：初始化区域系统n拉格朗日乘子的初值初始化原始对偶拟牛顿算法中区域系统n原始海塞矩阵初始化区域系统n增广耦合变量变化量初始化区域系统n增广耦合变量梯度值初始化区域系统n增广耦合变量梯度变化量初始化区域系统n拉格朗日乘子的海塞矩阵初始化区域系统n拉格朗日乘子的变化量初始化区域系统n拉格朗日乘子的梯度值初始化区域系统n拉格朗日乘子的梯度变化量初始化罚系数c；设置迭代次数k＝0；设置原始对偶拟牛顿算法收敛精度门槛值ε；

步骤3.2，利用初始化的原始对偶拟牛顿算法求解初始化的原始辅助子优化模型，得到各区域系统内部决策变量和增广耦合变量的迭代值；

设置迭代次数k＝k+1，各区域系统n∈{1,…,N}并行计算如下初始化的原始辅助子优化模型，得到各区域系统内部决策变量和增广耦合变量的第k次迭代值：

其中，为区域系统n内部决策变量的第k次迭代值；为区域系统n增广耦合变量的第k次迭代值；x_n为区域系统n的内部决策变量；为区域系统n的增广耦合变量；Ψ_n为区域系统n决策变量的可行域；为区域系统n增广耦合变量初值；为原始对偶拟牛顿算法中区域系统n原始海塞矩阵初值；为区域系统n拉格朗日乘子的初值；f_n为区域系统n的校正控制代价；

步骤3.3，所有邻居区域系统之间交换增广耦合变量的迭代值，各区域系统根据交换得到的邻居区域系统耦合变量的迭代值计算各区域系统副本变量的迭代值；

所有邻居区域系统之间交换增广耦合变量各区域系统n∈{1,…,N}根据下式计算区域系统n副本变量的第k次迭代值：

式中，为区域系统n副本变量的第k次迭代值；Δ_n为区域系统n邻居的集合；ξ_n,m为区域系统n与区域系统m的关联矩阵；为邻居区域系统m的增广耦合变量的第k次迭代值；

步骤3.4，各区域系统副本变量的迭代值带入初始化的对偶辅助子优化模型，计算各区域系统的拉格朗日乘子的迭代值；

各区域系统n∈{1,…,N}并行计算如下对偶辅助子优化模型，得到区域系统n拉格朗日乘子的第k次迭代值：

式中，为区域系统n拉格朗日乘子的第k次迭代值；为区域系统n拉格朗日乘子的初值；为区域系统n拉格朗日乘子的海塞矩阵初值；为区域系统n增广耦合变量的第k次迭代值；为区域系统n副本变量第k次迭代值；

步骤3.5，各区域系统利用步骤3.4得到各区域系统的拉格朗日乘子的迭代值，计算增广耦合变量梯度值的迭代值；

各区域系统n∈{1,…,N}根据下式计算区域系统n增广耦合变量梯度值的第k次迭代值

式中，为增广耦合变量梯度值的第k次迭代值；c为罚系数；为区域系统n副本变量第k次迭代值；为区域系统n增广耦合变量的第k次迭代值；

步骤3.6，各区域系统利用各区域系统的增广耦合变量的迭代值和各区域系统的增广耦合变量梯度值的迭代值，计算各区域系统的增广耦合变量变化量的迭代值和增广耦合变量梯度变化量的迭代值；

各区域系统n∈{1,...,N}通过下式计算区域系统n增广耦合变量变化量的第k次迭代值和区域系统n增广耦合变量梯度变化量的第k次迭代值：

式中，为区域系统n增广耦合变量变化量的第k次迭代值；D_n为区域系统n的正规化矩阵；为区域系统n增广耦合变量的第k次迭代值；为区域系统n增广耦合变量的第k-1次迭代值；为区域系统n耦合变量梯度变化量的第k次迭代值；为区域系统n的增广耦合变量梯度值的第k次迭代值；为区域系统n的增广耦合变量梯度值的第k-1次迭代值；σ为正则化参数，一般设置为较小的正数；

步骤3.7，各区域系统利用步骤3.6得到的各区域系统的增广耦合变量变化量的迭代值和增广耦合变量梯度变化量的迭代值计算各区域系统的增广耦合变量原始海塞矩阵的迭代值；

各区域系统n∈{1,…,N}基于如下分散式(Decentralized Broyden-Fletcher-Goldfarb-Shanno，D-BFGS)迭代计算区域系统n原始海塞矩阵的第k次迭代值；若则：

式中，为区域系统n原始海塞矩阵的第k次迭代值；为区域系统n原始海塞矩阵的第k-1次迭代值；为区域系统n耦合变量梯度变化量的第k次迭代值；为区域系统n增广耦合变量变化量的第k次迭代值；I为单位矩阵；σ为正则化参数，一般设置为较小的正数；

若则：

式中，为区域系统n原始海塞矩阵的第k次迭代值；为区域系统n原始海塞矩阵的第k-1次迭代值；

步骤3.8，各区域系统基于步骤3.7得到的各区域系统的增广耦合变量海塞矩阵的迭代值，计算各区域系统的耦合变量更新方向本地值的迭代值；

各区域系统n∈{1,…,N}根据下式并行计算区域系统n耦合变量更新方向本地值的第k次迭代值：

式中，为区域系统n耦合变量更新方向本地值的第k次迭代值；为区域系统n原始海塞矩阵的第k次迭代值；Ψ为正则化参数，一般设置为较小的正数；D_n为区域系统n的正规化矩阵；为区域系统n拉格朗日乘子的第k次迭代值；c为罚系数；为区域系统n增广耦合变量的第k次迭代值；为区域系统n副本变量第k次迭代值；

步骤3.9，所有邻居区域系统之间交换由步骤3.8得到的增广耦合变量的更新方向本地值的迭代值，各区域系统根据交换得到的邻居区域系统的增广耦合变量的更新方向本地值的迭代值计算各区域系统的增广耦合变量的近端值的迭代值；

所有邻居区域系统之间交换各区域系统n∈{1,…,N}由下式并行计算区域系统n增广耦合变量近端值的第k次迭代值

式中，为区域系统n增广耦合变量近端值的第k次迭代值；为区域系统n增广耦合变量的第k次迭代值；为区域系统n迭代步长；Δ_n为区域系统n邻居的集合；为区域系统n缩放系数；为区域系统n的邻居区域系统i耦合变量更新方向本地值的第k次迭代值；其中，迭代步长具有如下形式的闭式解：

式中，为区域系统n拉格朗日乘子的第k次迭代值；Δ_n为区域系统n邻居的集合；为区域系统n增广耦合变量的第k次迭代值；为区域系统n副本变量第k次迭代值；

步骤3.10，各区域系统根据步骤3.9得到的增广耦合变量近端值，计算各区域系统的增广耦合变量的迭代值；

各区域系统n∈{1,…,N}通过求解下式并行计算内部决策变量和增广耦合变量的第k+1次迭代值：

其中，为区域系统n内部决策变量的第k+1次迭代值，为区域系统n增广耦合变量的第k+1次迭代值；x_n为区域系统n的内部决策变量；为区域系统n的增广耦合变量；Ψ_n为区域系统n决策变量的可行域；为区域系统n缩放系数；f_n为区域系统n的校正控制代价；为区域系统n增广耦合变量的第k次迭代值；为原始对偶拟牛顿算法中区域系统n原始海塞矩阵的第k次迭代值；为区域系统n增广耦合变量近端值的第k次迭代值；

步骤3.11，所有邻居区域系统交换由步骤3.10得到的增广耦合变量的迭代值，各区域系统根据交换得到的邻居区域系统的增广耦合变量的迭代值计算各区域系统的副本变量的迭代值；

各区域系统n∈{1,…,N}在邻居之间交换区域系统n增广耦合变量的迭代值各区域系统n∈{1,…,N}根据交换得到的邻居区域系统的增广耦合变量的次迭代值，利用下式计算区域系统n副本变量的第k+1次迭代值：

式中，为为区域系统n副本变量第k+1次迭代值；Δ_n为区域系统n邻居的集合；ξ_n,m为区域系统n与区域系统m的关联矩阵；为区域系统n的邻居区域系统m的增广耦合变量的第k+1次迭代值；

步骤3.12，根据步骤3.10得到的增广耦合变量和步骤3.11得到的副本变量，各区域系统计算增广耦合变量梯度值的迭代值；

各区域系统n∈{1,…,N}通过下式计算区域系统n增广耦合变量梯度值的第k次迭代值

式中，为区域系统n增广耦合变量梯度值的第k次迭代值；为区域系统n的增广耦合变量的第k+1次迭代值；为区域系统n副本变量第k+1次迭代值；

步骤3.13，各区域系统根据步骤3.4得到的拉格朗日乘子的迭代值和步骤3.12得到的拉格朗日乘子梯度值的迭代值，计算各区域系统的拉格朗日乘子变化量的迭代值和各区域系统的拉格朗日乘子梯度变化量的迭代值；

各区域系统n∈{1,...,N}根据下式计算区域系统n的拉格朗日乘子变化量的第k次迭代值和区域系统n的拉格朗日乘子梯度变化量的第k次迭代值：

式中，区域系统n的拉格朗日乘子变化量的第k次迭代值；D_n为区域系统n的正规化矩阵；为区域系统n拉格朗日乘子的第k次迭代值；为区域系统n拉格朗日乘子的第k-1次迭代值；为区域系统n的拉格朗日乘子梯度变化量的第k次迭代值；为区域系统n拉格朗日乘子梯度值的第k次迭代值；为区域系统n拉格朗日乘子梯度值的第k-1次迭代值；γ为正则化参数，一般设置为较小的正数；

步骤3.14，各区域系统根据步骤3.13得到的各区域系统的拉格朗日乘子的变化量的迭代值和各区域系统的拉格朗日乘子的梯度变化量的迭代值，计算各区域系统拉格朗日乘子的海塞矩阵的迭代值；

各区域系统n∈{1,...,N}基于如下D-BFGS迭代计算区域系统n拉格朗日乘子的海塞矩阵的第k次迭代值若则：

若则：

式中，为区域系统n拉格朗日乘子的海塞矩阵的第k次迭代值；为区域系统n拉格朗日乘子的海塞矩阵的第k-1次迭代值；为区域系统n的拉格朗日乘子梯度变化量的第k次迭代值；区域系统n的拉格朗日乘子变化量的第k次迭代值；I为单位矩阵；γ为正则化参数，一般设置为较小的正数；

步骤3.15，各区域系统根据步骤3.14得到的各区域系统的拉格朗日乘子的海塞矩阵的迭代值，计算各区域系统拉格朗日乘子的更新方向本地值的迭代值；

各区域系统n∈{1,...,N}根据下式并行计算区域系统n拉格朗日乘子更新方向的本地值的第k次迭代值：

式中，为区域系统n拉格朗日乘子更新方向的本地值的第k次迭代值；为区域系统n拉格朗日乘子的海塞矩阵的第k次迭代值；D_n为区域系统n的正规化矩阵；Υ为正则化参数，一般设置为较小的正数；为区域系统n的增广耦合变量的第k+1次迭代值；为区域系统n副本变量第k+1次迭代值；

步骤3.16，所有邻居区域系统之间交换由步骤3.13得到的拉格朗日乘子的更新方向本地值的迭代值，各区域系统根据交换得到的邻居区域系统的拉格朗日乘子的迭代值的更新方向计算各区域系统的拉格朗日乘子的迭代值；

所有邻居区域系统之间交换各区域系统n∈{1,...,N}根据各自邻居的响应信息计算区域系统n拉格朗日乘子的第k+1次迭代值：

式中，为区域系统n拉格朗日乘子的第k+1次迭代值；为区域系统n拉格朗日乘子的第k次迭代值；Δ_n为区域系统n邻居的集合；κ_n为区域系统n的拉格朗日乘子的迭代步长；为区域系统n的邻居区域系统i的拉格朗日乘子更新方向的本地值的第k次迭代值；

步骤3.17，各区域系统根据收敛条件判断算法收敛状态，若算法收敛则停止迭代输出各区域系统的校正控制策略，若算法未收敛则进入下一轮迭代；

根据步骤3.10得到的增广耦合变量和步骤3.11得到的副本变量，各区域系统n∈{1,...,N}通过收敛条件判断算法的收敛状态；其中，收敛条件可设置为但不仅限于下式形式：

式中，N为区域系统的总数；为区域系统n的增广耦合变量的第k+1次迭代值；为区域系统n副本变量第k+1次迭代值；ε为原始对偶拟牛顿算法收敛精度门槛值；

判断是否满足收敛条件，若满足收敛条件，则停止算法，返回各区域系统n∈{1,...,N}的最优值f_n(P_n)及最优解并输出最有校正控制策略；否则，设置迭代次数k＝k+1，返回步骤3.5，进入下一轮迭代。