[发明专利]计及电-气耦合系统约束的燃气电厂储气容量计算方法

权利要求书

1.一种计及电-气耦合系统约束的燃气电厂储气容量计算方法，其特征在于该方法包括以下步骤：

(1)建立一个电-气耦合系统中的燃气电厂发电成本的优化目标函数如下：

C_coal(G_coal)＝P_coal×G_coal

C_gas(L_gas)＝P_gas×L_gas

其中，C为电力系统运行的总成本，C_coal(G_coal)为燃煤成本，G_coal为耗煤量，P_coal为煤炭的价格，C_gas(L_gas)为燃气成本，L_gas为耗气量，P_gas为天然气的购买价格，C_storage(V)为储气成本，V为储气设备的容积，T_life为储气设备的使用年限；

(2)建立电-气耦合系统稳态安全运行的等式约束方程，包括：

(2-1)一个电-气耦合系统中电力系统的潮流方程如下：

其中，Pⁱ为电力系统中第i个节点的注入有功功率，Qⁱ为电力系统中第i个节点的注入无功功率，G_ij为与电力系统的节点导纳矩阵Y中第i行、第j列相对应的电导，B_ij为与电力系统的节点导纳矩阵Y中第i行、第j列相对应的电纳，电力系统节点导纳矩阵Y从电力系统调度中心获取；

(2-2)一个电-气耦合系统中天然气网中管道的水力方程如下：

其中，f_km为天然气网中第k个节点和第m个节点之间的管道中天然气体积流量，p_k，p_m分别为第k个节点和第m个节点的压强，D^km和L^km分别为第k个节点和第m个节点之间管道km的管道直径和管道长度，F为管道内壁的摩擦系数，F由公式计算得到，E_f为管道的效率系数，E_f取值为0.92，Re为雷诺数，由公式计算得到，ρ为天然气密度，μ_m为天然气的动力粘度系数，与天然气密度有关，查表获得，γ_G为天然气比重，0＜γ_G＜1，T_a为天然气平均温度，T_n和p_n分别为标准状态下天然气的温度和压力，Z_g为天然气平均可压缩系数，0.9＜Z_g＜1.5，上述天然气网中管道的水力方程中，当时，上式中的sgn_p(p_k,p_m)＝1，当时，sgn_p(p_k,p_m)＝-1；

(2-3)一个电-气耦合系统的天然气网中压缩机的能量消耗方程如下：

其中，p_k，p_m分别为天然气网中第k个节点和第m个节点的压强，BHP^km为第k个节点和第m个节点之间的压缩机的能量消耗，为该压缩机的入口体积流量，η_c为该压缩机的总效率，c_k为该压缩机的多变系数，η_c和c_k从压缩机的出厂说明书获取；

(2-4)一个电-气耦合系统中通过燃气轮机耦合的电力系统与天然气网之间的耦合方程如下：

μP_Tur＝H_gas×[f_Tur+(V_t-1-V_t)]，

其中，f_Tur为燃气轮机由管道输入的燃气的体积流量，P_Tur为燃气轮机的有功功率输出，V_t为储气设备在t时刻的天然气储量，H_gas为天然气的燃烧热值，取值为37.59MJ/m³，μ为燃气轮机的效率系数，由燃气轮机的出厂说明书获取；

(2-5)一个电-气耦合系统中天然气网的节点气流平衡方程如下：

A_Gf＝L，

其中，A_G为天然气网的节点-支路矩阵，f为天然气网支路的体积流量，L为天然气网节点的气体负荷，L根据天然气网历史运行数据得到；

(3)设定电-气耦合系统稳态安全运行的不等式约束条件，包括：

(3-1)电力系统中发电机组的输出功率P_i^gen大于或等于0，小于或等于该发电机组出厂铭牌上给出的最大功率即：

(3-2)电力系统第i个节点的电压幅值U_i在设定的电力系统安全运行电压的上限值和下限值U_i之间运行，U_i为第i个节点额定电压的0.95倍，为第i个节点额定电压的1.05倍，即：

(3-3)电力系统中第l条线路的传输容量S_l小于或等于设定的电力系统安全运行传输容量的最大值即：

(3-4)天然气网中第k个节点的压强p_k在设定的管道安全运行气压的上限值、下限值p_k、内，即：

(3-5)天然气网中第b个管道的流量f_b在设定的管道安全运行流量的上限值、下限值f_b、内，即：

(3-6)天然气网中气源供气量f_s小于或等于该气源能提供天然气流的最大值f_s,max，即：

f_s≤f_s,max；

(3-7)天然气网中的压缩机安全运行的约束条件：

1S＜S_max，p^out≤p_c,max，

其中：S是该压缩机的升压比，S_max是该压缩机的最大升压比，S_max由压缩机的出厂铭牌获取，为该压缩机的入口的体积流量，为该压缩机的入口最大允许体积流量，由压缩机的出厂铭牌获取，p^out为压缩机的出口压力，p_c,max为压缩机的出口最大允许压力，p_c,max由压缩机的出厂铭牌获取；

(3-8)燃气电厂中储气设备在t时刻的储量应大于等于0，小于等于储气设备的最大容积：

0≤V_t≤V_max；

(4)利用优化算法，求解由上述步骤(1)为目标函数、上述步骤(2)和步骤(3)为约束条件的优化模型，得到储气设备在预测年负荷下的月最优储气容量曲线，最优储气容量曲线中的横坐标为12个月，最优储气容量曲线中的纵坐标为月储气容量；

(5)根据上述步骤(4)的最优储气容量曲线，求取容量的期望值，将最接近期望值的容量作为燃气电厂中储气设备的最优容量。