[发明专利]考虑天然气系统约束的等值最优潮流模型构建方法

权利要求书

1.一种考虑天然气系统约束的等值最优潮流模型构建方法，其特征在于，包括以下步骤:

1)输入基础数据

输入电-气互联系统基础数据：电力子系统设备参数、天然气子系统设备参数和拓扑结构

2)建立天然气系统的约束等值模型

2-1)建立约束等值目标函数

根据天然气系统中燃气轮机个数N建立多个目标函数，首先建立任意一个燃气轮机消耗天然气流量最大的目标函数，一共N个；然后建立任意两个燃气轮机消耗天然气流总流量最大的目标函数，一共C_2N个；以此类推，最后建立所有燃气轮机消耗天然气总流量最大的目标函数，一共C_NN个；即由C_i，C_ij，…，C_all组成；

式中，在建立任意一个燃气轮机消耗天然气流量最大的目标函数时，C_i表示第i个燃气轮机消耗天然气流量的最大值，表示第i个燃气轮机消耗天然气流量；在建立任意两个燃气轮机消耗天然气总流量最大的目标函数时，C_ij为第i个和第j个燃气轮机消耗天然气总流量的最大值，和分别表示第i个和第j个燃气轮机消耗的天然气流量；在建立所有燃气轮机消耗天然气总流量最大的目标函数时，C_all为所有燃气轮机消耗天然气总流量的最大值，表示第i个燃气轮机消耗天然气流量；在求解燃气轮机消耗天然气流量最大值时，其他不相关燃气轮机消耗天然气流量为0；

2-2)建立约束等值等式约束方程

天然气系统流量平衡方程为式(2)：

F_G,m-F_GAS,m-F_L,m-F_m＝0,m＝1,2,...,N_m (2)

式中，F_G,m、F_GAS,m、F_L,m、F_m分别为天然气节点m的气源流量、燃气轮机消耗的流量、气负荷和注入流量，N_m为天然气系统节点个数；其中节点注入流量F_m通过式(3)求得；

式中，天然气节点n表示与天然气节点m相连的节点；F_mn、F_com,mn、τ_com,mn分别为管道流量、流过压缩机的流量和压缩机损耗，sgn_c(m,n)和sgn_τ(m,n)分别为压缩机流过流量和消耗流量的方向函数；这些量通过式(4)和(5)表示：

式中，Π_m是天然气节点m的节点气压；k_mn是与管道内径、长度、效率、压缩因子相关的常数；s_mn反应管道流量方向；H_com,mn为压缩机消耗的功率；B_mn是压缩机系数；Z_mn是压缩机压缩因子相关的常数；α_τ、β_τ、γ_τ为功率转化效率常数；

2-3)建立约束等值不等式约束方程

式(6)-(8)分别为天然气系统节点气压、气源容量、压缩机气压变比的上下限；

式中，Π_m是天然气节点m的节点气压，F_G,m为气源容量，R_r为压缩机气压变比；和分别为天然气节点气压的上限和下限；

和分别为气源容量的上限和下限；和分别为压缩机压变比的上限和下限；N_s和N_p分别为气源和压缩机的个数；

求解由公式(1)-(8)构成的优化模型，进而计算得到目标函数中的通过燃气轮机耗量方程(9)计算出对应的再通过公式(10)计算得到任意一个燃气轮机的出力上限，一共N个；得到任意两个燃气轮机总出力上限，一共C_2N个；以此类推，得到所有燃气轮机总出力上限，一共C_NN个；即由组成；

式(10)中，在求解任意一个燃气轮机的出力上限时，为第i个燃气轮机输出功率的最大值，为第i个燃气轮机的输出功率；在求解任意两个燃气轮机总出力上限时，为第i个和第j个燃气轮机总输出功率的最大值，和分别表示第i个和第j个燃气轮机输出功率；在求解所有燃气轮机总出力上限时，为所有燃气轮机总输出功率的最大值，为第i个燃气轮机的输出功率；GHV为天然气高热值；

3)求解天然气系统的约束等值模型

式(1)-(10)构成求解天然气约束等值参数的非线性最优化问题；通过内点法直接求得系统中各个燃气轮机组的有功出力约束以替代天然气系统的运行约束；

4)建立考虑天然气系统约束的等值最优潮流模型

4-1)等值最优潮流目标函数：

式中，f_i(P_GAS,i)为燃气轮机组的成本函数；f_j(P_Gj)为非燃气发电机的成本函数；N_G为非燃气发电机个数；

4-2)建立等值最优潮流等式约束方程：

考虑电力系统功率平衡方程(12)-(13)，构建如下等式约束：

式中，P_i和Q_i分别为节点i的注入有功功率和无功功率；U_i是节点i的电压幅值；θ_ij是节点i与j的电压相角差；G_ij和B_ij分别为节点导纳矩阵的第i行第j列元素的实部和虚部；N_A为节点个数；

4-3)建立等值最优潮流不等式约束方程：

考虑天然气系统的等值约束燃气轮机自身的容量约束、非燃气常规机组容量约束、节点电压约束和线路功率约束构建如下不等式约束：

T_l^min≤T_l≤T_l^max l＝1,2,...,N_L (20)

其中，P_GAS,i为燃气轮机i的有功功率，Q_GAS,i为燃气轮机i的无功功率，P_Gi为非燃气轮机i的有功功率，Q_Gi为非燃气轮机i的无功功率，U_i是节点i的电压幅值，T_l是输电线路l流过的功率；N_L为输电线路的总数；和分别为燃气轮机有功出力的上限和下限；和分别为燃气轮机无功出力的上限和下限；和分别为节点电压幅值的上限和下限；和为输电线路传输功率的上限和下限；和分别为非燃气常规机组有功出力的上限和下限；和分别为非燃气常规机组无功出力的上限和下限；

5)求解考虑天然气系统约束的等值最优潮流模型

式(11)-(20)构成求解考虑天然气系统运行约束的最优潮流的非线性最优化问题；通过内点法直接求得电力系统各发电机的最优出力。