一种电‑气耦合多能流系统的优化调度方法

摘要

本发明涉及一种电‑气耦合多能流系统优化调度方法，属于含多种能源形式的电力系统运行和控制技术领域。本方法考虑电‑气系统的相互影响，实现了电‑气耦合多能流系统的优化调度。相比独立地对供电、供气系统进行优化调度分析，不仅能得到更优的调度方案(总运行成本或网络损耗更小等)，还提高了调度的灵活性。该方法可以应用于电‑气耦合多能流系统的调度计划制定，有利于提高电‑气耦合多能流系统的用能效率，减少运行成本。

主权项

一种电‐气耦合多能流系统的优化调度方法，其特征在于，该方法包括以下步骤：(1)建立一个电‑气耦合多能流系统优化调度的目标函数：$min\underset{x\in GU}{\Sigma }{W}_x^{gas}+\underset{y\notin GU}{\Sigma }{F}_y^C\left(P_y\right)$其中，GU为电‐气耦合多能流系统中的燃气机组编号的集合，为电‐气耦合多能流系统中第x台燃气机组的燃料成本，为电‐气耦合多能流系统中第y台非燃气机组的运行成本，P_y为电‐气耦合多能流系统中第y台非燃气机组的有功功率，非燃气机组的运行成本是机组有功出力的非线性函数，从非燃气机组的出厂说明书获得；(2)设定电‑气耦合多能流系统中电网与天然气网稳态安全运行的等式约束条件，包括：(2‑1)电‑气耦合多能流系统中电网潮流方程如下：$\begin{array}{c}{P}^i=U_i\underset{j\in i}{\Sigma }{U}_j\left(G_{ij}\cos \right({\theta }_i-{\theta }_j)+B_{ij}\sin ({\theta }_i-{\theta }_j)),i,j=1,2,...n\\ Q^i=U_i\underset{j\in i}{\Sigma }{U}_j\left(G_{ij}\sin \right({\theta }_i-{\theta }_j)-B_{ij}\cos ({\theta }_i-{\theta }_j)),i,j=1,2,...n\end{array},$其中，Pⁱ为电网中节点i的注入有功功率，Qⁱ为电网中节点i的注入无功功率，θ_i、θ_j分别为节点i、节点j的电压相角，U_i、U_j分别为节点i、节点j的电压幅值，G_ij为电网节点导纳矩阵Y第i行、第j列元素的实部，B_ij为电网节点导纳矩阵Y第i行、第j列元素的虚部，电网节点导纳矩阵Y从电‑气耦合多能流系统的能量管理系统中获取；(2‑2)电‑气耦合系统的天然气网中管道的水力方程如下：$f_{km}={\mathrm{sgn}}_p(p_k,p_m)\times \left(\frac{77.54T_n}{p_n}\right){\left(D^{km}\right)}^{2.5}{\left(\frac{1}{L^{km}{\gamma }_G{T}_a^{km}{\mathrm{FZ}}_g}\right)}^{0.5}\sqrt{(p_k^2-p_m^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为雷诺数，Re＝ρvd/μ，ρ为天然气密度，v为天然气流速，μ为天然气黏性系数，取值为0＜u＜100帕秒，d为天然气管道直径，γ_G为天然气比重，0＜γ_G＜1，T_a为天然气平均温度，T_n、p_n分别为标准状态下天然气的温度和压力，T_n、p_n的取值分别为288K、0.1Mpa，Z_g为天然气平均可压缩系数，取值为0＜Z_g＜2，天然气网管道水力方程中，当时，上式中的sgn_p(p_k,p_m)＝1，当时，sgn_p(p_k,p_m)＝‑1；(2‑3)电‑气耦合系统中的天然气网中压缩机的能量消耗方程如下：${\mathrm{BHP}}^{km}=\frac{f_C^{in}{c}_k}{{\eta }_C(c_k-1)}[{\left(\frac{p_m}{p_k}\right)}^{\frac{c_k-1}{c_k}}-1],$其中，p_k，p_m分别为天然气网中第k个节点和第m个节点的压强，BHP^km为第k个节点和第m个节点之间的压缩机的能量消耗，为该压缩机的入口体积流量，η_c为该压缩机的总效率，c_k为该压缩机的多变系数，η_c和c_k从压缩机的出厂说明书获取；(2‑4)通过燃气轮机耦合的电网与天然气网之间的耦合方程如下：$P_{Tur}=C_1{f}_{Tur}^2+C_2{f}_{Tur}+C_3,$其中，f_Tur为燃气轮机的燃气的体积流量，P_Tur为燃气轮机的有功功率输出，C₁、C₂和C₃分别为燃气轮机的燃气系数，由燃气轮机的出厂说明书获取；(2‑5)电‑气耦合系统中天然气网的节点气流平衡方程如下：A_Gf＝L，其中：A_G为天然气网的节点‑支路矩阵，f为天然气网的所有支路的体积流量，L为天然气网的所有节点的天然气注入体积流量，正数表示气源注入天然气网的体积流量，负数表示负荷从天然气网获得的体积流量，L根据天然气网历史运行数据得到；(3)设定电‑气耦合系统中电网与天然气网稳态安全运行的不等式约束条件，包括：(3‑1)电‑气耦合系统的电网中第i个节点的电压幅值U_i在设定的电力系统安全运行电压的上、下限值U_i、之间运行，U_i为第i个节点额定电压的0.95倍，为第i个节点额定电压的1.05倍：${\underset{‾}{U}}_i\le U_i\le {\bar{U}}_i;$(3‑2)电‑气耦合系统的电网中第l条线路的传输容量小于或等于设定的电网安全运行传输容量的最大值$S_l\le {\bar{S}}_l;$(3‑3)电‑气耦合系统的电网中燃气机组有功功率的爬坡约束：$-{\mathrm{RAMP}}_x^{down}·\Delta t\le P_{x,t}-P_{x,t-1}\le {\mathrm{PAMP}}_x^{up}·\Delta t;$其中，和分别为第x台燃气机组有功功率的向上和向下爬坡速率，和从燃气机组的出厂说明书中获得，Δt为相邻两个调度时段的时间间隔，P_x,t和P_x,t‑1分别为第x台燃气机组在第t个调度时段和第t‑1个调度时段的有功功率；(3‑4)电‑气耦合系统的电网中非燃气机组有功功率的爬坡约束：$-{\mathrm{ramp}}_y^{down}·\Delta t\le P_{y,t}-P_{y,t-1}\le {\mathrm{ramp}}_y^{up}·\Delta t;$其中，和分别为第y台非燃气机组有功功率的向上和向下爬坡速率，和从非燃气机组的出厂说明书中获得，Δt为相邻两个调度时段的时间间隔，P_y,t和P_y,t‑1分别为第y台非燃气机组在第t个调度时段和第t‑1个调度时段的有功功率；(3‑5)电‑气耦合系统的电网中第x台燃气机组的有功功率P_x在设定的电网安全运行第x台燃气机组有功功率的上、下限值P_x之间：${\underset{‾}{P}}_x\le P_x\le {\bar{P}}_x;$(3‑6)电‑气耦合系统的电网中第y台非燃气机组的有功功率P_y在设定的电网安全运行第y台非燃气机组有功功率的上、下限值P_y之间：${\underset{‾}{P}}_y\le P_y\le {\bar{P}}_y;$(3‑7)电‑气耦合系统的天然气网中第k个节点的压强p_k在设定的管道安全运行气压的上、下限值p_k、内：${\underset{‾}{p}}_k\le p_k\le {\bar{p}}_k;$(3‑8)电‑气耦合系统的天然气网中第b个管道的流量f_b在设定的管道安全运行流量的上、下限值f_b、内：$0\le f_b\le {\bar{f}}_b;$(3‑9)电‑气耦合系统的天然气网中气源供气量f_s小于或等于该气源能提供天然气流的最大值f_s,max：f_s≤f_s,max；(3‑10)电‑气耦合系统的天然气网中的压缩机安全运行时需要满足以下约束：$1<S<S_{max}(S=\frac{p^{out}}{p^{in}}),f_{comp}^{in}\le f_{c,max}^{in},p^{out}\le p_{c,max},$其中：S是该压缩机的升压比，S_max是该压缩机的最大升压比，S_max由压缩机的出厂铭牌获取，为该压缩机的入口的体积流量，为该压缩机的入口最大允许体积流量，由压缩机的出厂铭牌获取，p^out为压缩机的出口压力，p_c,max为压缩机的出口最大允许压力，p_c,max由压缩机的出厂铭牌获取；(4)采用内点法，将步骤(1)中的方程作为目标函数，将上述步骤(2)和步骤(3)的所有方程作为约束条件，求解得到电‑气耦合多能流系统中燃气轮机的有功功率、非燃气轮机的有功功率和气源注入天然气网的体积流量，作为电‐气耦合多能流系统的优化调度方案。