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

摘要

本发明公开一种考虑天然气系统约束的等值最优潮流模型构建方法。在电力系统与天然气系统数据各自独立调度运行，数据不能共享的情况下，提出天然气系统的约束等值模型，将天然气系统等值约束信息提供给电力系统而不涉密；在此基础上建立考虑天然气系统约束的等值最优潮流方法，避免了电力系统与天然气系统的交替迭代，同时还有效计及了天然气系统对电力系统运行调度的影响。最后，通过IEEE14‑NGS10电‑气互联系统的算例分析验证了所提方法的有效性。

主权项

一种考虑天然气系统约束的等值最优潮流模型构建方法，其特征在于，包括以下步骤：1)输入基础数据输入电‑气互联系统基础数据：电力子系统设备参数、天然气子系统设备参数和拓扑结构2)建立天然气系统的约束等值模型2‑1)建立约束等值目标函数根据天然气系统中燃气轮机个数(N)建立多个目标函数，首先建立任意一个燃气轮机消耗天然气流量最大的目标函数，一共N个；然后建立任意两个燃气轮机消耗天然气流总流量最大的目标函数，一共个；以此类推，最后建立所有燃气轮机消耗天然气总流量最大的目标函数，一共个；即由Ci，Cij，…，Call组成。Ci=MaxFGAS,i(1)i=1,2,...,NCij=Max(FGAS,l(2)+FGAS,j(2))i,j=1,2,...,N;i>jCall=MaxΣi=1NFGAS,i(N)i=1,2,...,N---(1)]]>式中，在建立任意一个燃气轮机消耗天然气流量最大的目标函数时，Ci表示第i个燃气轮机消耗天然气流量的最大值，表示第i个燃气轮机消耗天然气流量。在建立任意两个燃气轮机消耗天然气总流量最大的目标函数时，Cij为第i个和第j个燃气轮机消耗天然气总流量的最大值，和分别表示第i个和第j个燃气轮机消耗的天然气流量。在建立所有燃气轮机消耗天然气总流量最大的目标函数时，Call为所有燃气轮机消耗天然气总流量的最大值，表示第i个燃气轮机消耗天然气流量。在求解燃气轮机消耗天然气流量最大值时，其他不相关燃气轮机消耗天然气流量为0。2‑2)建立约束等值等式约束方程天然气系统流量平衡方程为式(2)。FG，m‑FGAS，m‑FL，m‑Fm＝0，m＝1，2，...，Nm  (2)式中，FG，m、FGAS，m、FL，m、Fm分别为天然气节点m的气源流量、燃气轮机消耗的流量、气负荷和注入流量，Nm为天然气系统节点个数。其中节点注入流量Fm可以通过式(3)求得。式中，天然气节点n表示与天然气节点m相连的节点；Fmn、Fcom，mn、τcom，mn分别为管道流量、流过压缩机的流量和压缩机损耗，sgnc(m，n)和sgnτ(m，n)分别为压缩机流过流量和消耗流量的方向函数。这些量可以通过式(4)和(5)表示。Fmn=kmnsmnsmn(Πm2-Πn2)smn=+1,Πm2-Πn2≥0-1,Πm2-Πn2≤0---(4)]]>τcom,mn=ατ+βτHcom,mn+γτHcom,mn2Hcom,mn=BmnFcom,mn[(ΠnΠm)zmn-1]---(5)]]>式中，Πm是天燃气节点m的节点气压；kmn是与管道内径、长度、效率、压缩因子相关的常数；smn反应管道流量方向；Hcom，mn为压缩机消耗的功率；Bmn是压缩机系数；Zmn是压缩机压缩因子相关的常数。ατ、βτ、γτ为功率转化效率常数。2‑3)建立约束等值不等式约束方程式(6)‑(8)分别为天然气系统节点气压、气源容量、压缩机气压变比的上下限。Πmmin≤Πm≤Πmmax,m=1,2,...,Nm---(6)]]>FG,mmin≤FG,m≤FG,mmax,m=1,2,...,NS---(7)]]>Rrmin≤Rr≤Rrmax,r=1,2,...,NP---(8)]]>式中，Πm是天燃气节点m的节点气压，FG，m内气源容量，Rr为压缩机气压变比。和分别为天然气节点气压的上限和下限；和分别为气源容量的上限和下限；和分别为压缩机压变比的上限和下限。Ns和Np分别为气源和压缩机的个数。求解由公式(1)‑(8)构成的优化模型可计算得到目标函数中的通过燃气轮机耗量方程(9)可计算出对应的再通过公式(10)计算得到任意一个燃气轮机的出力上限，一共N个；得到任意两个燃气轮机总出力上限，一共个；以此类推，得到所有燃气轮机总出力上限，一共个；即由组成。FGAS,i(1)=(αi+βiPGAS,i(1)+γi(PGAS,i(1))2)/GHVi=1,2,...,NFGAS,i(2)=(αi+βiPGAS,i(2)+γi(PGAS,i(2))2)/GHVi=1,2,...,NFGAS,i(N)=(αi+βiPGAS,i(N)+γi(PGAS,i(N))2)/GHVi=1,2,...,N---(9)]]>Pimax=PGAS,i(1)i=1,2,...,NPijmax=PGAS,i(2)+PGAS,j(2)i,j=1,2,...,N;i>jPallmax=Σi=1NPGAS,i(N)i=1,2,...,N---(10)]]>式(10)中，在求解任意一个燃气轮机的出力上限时，为第i个燃气轮机输出功率的最大值，为第i个燃气轮机的输出功率。在求解任意两个燃气轮机总出力上限时，为第i个和第j个燃气轮机总输出功率的最大值，和分别表示第i个和第j个燃气轮机输出功率。在求解所有燃气轮机总出力上限时，为所有燃气轮机总输出功率的最大值，为第i个燃气轮机的输出功率。GHV为天然气高热值。3)求解天然气系统的约束等值模型式(1)‑(10)构成求解天然气约束等值参数的非线性最优化问题。通过内点法可直接求得系统中各个燃气轮机组的有功出力约束以替代天然气系统的运行约束。4)建立考虑天然气系统约束的等值最优潮流模型4‑1)等值最优潮流目标函数：Minf=Σi=1Nfi(PGAS,i)+Σj=1NGfj(PGj),i=1,...,N;j=1,...,NG---(11)]]>式中，fi(PGAS，i)为燃气轮机组的成本函数；fj(PGj)为非燃气发电机的成本函数；NG为非燃气发电机个数。4‑2)建立等值最优潮流等式约束方程：考虑电力系统功率平衡方程(12)‑(13)，构建如下等式约束：Pi-UiΣj=1NAUj(Gijcosδij+Bijsinδij)=0,i=1,2,...,NA---(12)]]>Qi-UiΣj=1NAUj(Gijsinδij-Bijcosδij)=0,i=1,2,...,NA---(13)]]>式中，Pi和Qi分别为节点i的注入有功功率和无功功率；Ui是节点i的电压幅值；θij是节点i与j的电压相角差；Gij和Bij分别为节点导纳矩阵的第i行第j列元素的实部和虚部；NA为节点个数。4‑3)建立等值最优潮流不等式约束方程：考虑天然气系统的等值约束(14)、燃气轮机自身的容量约束(15)‑(16)、非燃气常规机组容量约束(17)‑(18)、节点电压约束(19)和线路功率约束(20)构建如下不等式约束。0≤PGAS,i≤Pimaxi=1,2,...,N0≤PGAS,i+PGAS,j≤Pijmaxi,j=1,2,...,N;i>j0≤Σi=1NPGAS,i≤Pallmaxi=1,2,...,N---14)]]>PGAS,imin≤PGAS,i≤PGAS,imax,i=1,2,...,N---(15)]]>QGAS,imin≤QGAS,i≤QGAS,imax,i=1,2,...,N---(16)]]>PGimin≤PGi≤PGimax,i=1,2,...,NG---(17)]]>QGimin≤QGi≤QGimax,i=1,2,...,NG---(18)]]>Uimin≤Ui≤Uimax,i=1,2,...,NA---(19)]]>Tlmin≤Tl≤Tlmax l＝1，2，...，NL  (20)其中，PGAS，i为燃气轮机i的有功功率功率，QGAS，i为燃气轮机i的无功功率，PGi为非燃气轮机i的有功功率功率，QGi为非燃气轮机i的无功功率，Ui是节点i的电压幅值，T1是输电线路1流过的功率；NL为输电线路的总数。和分别为燃气轮机有功出力的上限和下限；和分别为燃气轮机无功出力的上限和下限；和分别为节点电压幅值的上限和下限；和为输电线路传输功率的上限和下限；和分别为非燃气常规机组有功出力的上限和下限；和分别为非燃气常规机组无功出力的上限和下限。5)求解考虑天然气系统约束的等值最优潮流模型式(11)‑(20)构成求解考虑天然气系统运行约束的最优潮流的非线性最优化问题。通过内点法可直接求得电力系统各发电机的最优出力。