[发明专利]一种基于故障耦合传播过程中的多能流负荷裕度求解方法

权利要求书

1.一种基于故障耦合传播过程中的多能流负荷裕度求解方法，其特征在于：该方法包括如下步骤：

步骤1：准备综合能源系统的原始数据，包括电力系统、天然气系统和热力系统的各节点、支路原始数据，通过能流状态估计设备或在线能量潮流计算分析，获取综合能源系统的当前运行状态；

步骤2：计算分析综合能源系统节点、支路、耦合传播路径及能量枢纽的静态能流稳定指标，这些静态能流稳定指标包括：天然气气压、流量的相关风险指标；热力系统温度、流量的相关风险指标；能量枢纽的气压、温度、电压等的相关风险指标；

步骤3：通过步骤2计算分析确定综合能源系统的关键支路、关键节点、关键耦合传播路径、及关键能量枢纽；

步骤4：形成扩展潮流方程，计算负荷裕度；

步骤5：通过多能流负荷裕度分析，得到多能流负荷安全域；

步骤6：基于多能流负荷安全域及其边界特征，获取当前运行点到相应安全域边界的距离，得到稳态情况下负荷裕度，得到能量枢纽的负荷裕度。

2.根据权利要求1所述的一种基于故障耦合传播过程中的多能流负荷裕度求解方法，其特征在于，步骤3中综合能源系统的关键环节的确定，是通过步骤2的静态能流稳定指标计算结果分析确定的。

3.根据权利1步骤4要求所述的形成扩展潮流方程，计算负荷裕度，其特征在于基于负荷参数的多阶导数进行初值预估，得到临界点初值a₀、零特征值对应的右特征向量初值b₀，以及在步骤3中多能系统负荷裕度初值λ₀；以上述参数为基础构建非线性潮流方程组以及相对应的修正方程组，通过降阶变换将修正方程组拆分成几组同系数低维矩阵的线性方程组；基于Riccati变换的线性树形多体系统特征值求解方法，建立多能流输入输出Riccati传递矩阵递推关系；从树形系统各输入端开始沿传递路径依次求得各能量枢纽的Riccati传递矩阵，进而实现综合能源系统中能量枢纽的负荷裕度的快速求解。