CN111523247A - 基于分立求解的综合能源系统多能流联合计算方法及装置 - Google Patents
基于分立求解的综合能源系统多能流联合计算方法及装置 Download PDFInfo
- Publication number
- CN111523247A CN111523247A CN202010395994.7A CN202010395994A CN111523247A CN 111523247 A CN111523247 A CN 111523247A CN 202010395994 A CN202010395994 A CN 202010395994A CN 111523247 A CN111523247 A CN 111523247A
- Authority
- CN
- China
- Prior art keywords
- natural gas
- node
- power
- energy
- flow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 56
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 356
- 239000003345 natural gas Substances 0.000 claims abstract description 182
- 230000008878 coupling Effects 0.000 claims abstract description 50
- 238000010168 coupling process Methods 0.000 claims abstract description 50
- 238000005859 coupling reaction Methods 0.000 claims abstract description 50
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000004088 simulation Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 34
- 238000002347 injection Methods 0.000 claims description 30
- 239000007924 injection Substances 0.000 claims description 30
- 239000007789 gas Substances 0.000 claims description 29
- 239000000243 solution Substances 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 4
- 238000012795 verification Methods 0.000 abstract 1
- 230000014509 gene expression Effects 0.000 description 6
- 238000010248 power generation Methods 0.000 description 4
- 238000002940 Newton-Raphson method Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供基于分立求解的电气耦合综合能源系统多能流联合计算方法及装置，包括：构建电气耦合综合能源系统多能流模型；使用耦合元件实现天然气系统和电力系统的解耦；解耦后分别运行电力系统潮流计算以及天然气系统潮流计算，得到各类状态变量作为多能流仿真的最终结果。本发明可通过能源耦合部分的能量等值将多能流待求信息分配到不同能源各自的求解模块，进行分立求解，物理意义明确。验证结果表明，该方法可有效计算IEGS的运行状态，为系统安全性分析提供支撑工具。
Description
技术领域
本发明属于综合能源系统技术领域，特别是涉及一种基于分立求解的电气耦合综合能源系统多能流联合计算方法及装置。
背景技术
随着经济社会的飞速发展以及环境的日益恶化，可持续发展理念深入人心。天然气作为一种优质、高效且清洁的低碳能源，在发电侧可与核能及可再生能源等其他低排放能源形成良性互补，是实现能源供应清洁化的现实选择之一。为此，天然气发电机组(Natural Gas-fired Generating Unit,NGU)由于具有污染排放少，响应速度快和发电效率高等优点，被广泛部署于电力系统。在世界范围内，天然气发电在总装机容量和发电份额中呈现稳中有进的良好“态势”，推动了电力系统同天然气系统的深度耦合。然而，电气耦合综合能源系统(Integrated Electricity-Gas System,IEGS)在提高能源利用效率，实现可再生能源充分消纳的同时，也带来了系统安全运行的整体性风险。随着电力系统同天然气系统之间的耦合日益紧密，二者间的交互影响所引起的安全问题需引起重视。将IEGS作为一个整体实现一体化安全稳定分析，对实现电力燃气系统的协同增效具有重要意义，而电气耦合的多能流计算是该过程的最基本环节。
针对IEGS的多能流计算，已有方法往往忽略天然气管网的运行特性和运行约束，以及电-气耦合环节的一些细节，仅从宏观角度分析电气与天然气系统之间的互动联系，计算结果往往存在一定误差。由于将电气系统联立增加了大量高维非线性方程与决策变量，使得统一求解法变得复杂。
发明内容
针对上述问题，考虑到天然气系统以及电力系统各自潮流计算都较成熟，本发明提出基于分立迭代的IEGS多能流计算方法，通过能源耦合部分的能量等值将多能流待求信息分配到不同能源各自的求解模块，进行分立求解，有效计算IEGS的运行状态，为系统安全性分析提供支撑工具。
为达到上述目的，本发明的技术方案是这样实现的：
基于分立求解的综合能源系统多能流联合计算方法，包括：构建电气耦合综合能源系统多能流模型；使用耦合元件实现天然气系统和电力系统的解耦；解耦后分别运行电力系统潮流计算以及天然气系统潮流计算，得到各类状态变量作为多能流仿真的最终结果。
进一步的，所述电气耦合综合能源系统多能流模型的构建过程为：
S101、建立天然气系统模型；
S102、建立电力系统模型；
S103、建立耦合元件的模型；
S104、建立电气耦合综合能源系统多能流模型：定义xeg和yeg为电气耦合综合能源系统IEGS的状态变量和注入变量；
xeg＝[V,θ,p]
yeg＝[Psp,Qsp,Lsp]
Feg(xeg,yeg)＝0
式中xeg由电力系统节点电压V、相角θ、天然气系统节点气压p组成；yeg由电力系统节点有功功率注入PSP、无功功率注入QSP、天然气系统节点注入LSP组成。
进一步的，所述耦合元件的解耦过程具体为：
将所述耦合元件在电力系统中等效为一个有功出力已知的发电机，在天然气系统中等效为一个燃气需求已知的天然气负荷，实现天然气系统和电力系统的解耦。
进一步的，所述电力系统潮流计算的方法包括：
S201、读取电网拓扑，网络参数，形成节点导纳矩阵；
S202、初始化电力系统状态变量V(0)，θ(0)；
S203、分别计算电力系统节点误差函数Fe(V，θ)；
S204、计算电力系统节点雅可比矩阵Je；
S205、计算迭代修正量ΔV，Δθ；
S206、更新电力系统状态变量，判断电力系统节点误差函数最大值是否小于阈值；若小于，迭代到此结束，进一步计算各线路潮流和平衡节点功率，反之进入步骤S203。
进一步的，所述天然气系统潮流计算的具体方法为：
S301、根据天然气网络拓扑形成全节点-支路关联矩阵以及降阶节点-支路关联矩阵，读取管道阻力系数，计算天然气系统各个节点的净负荷，初始化节点压力；
S302、以节点压力平方为变量，求解天然气节点误差函数；
S303、计算新的节点压力平方，判断天然气节点误差函数最大值是否小于阈值；若小于，得到节点压力平方真值，进入步骤S304，否则进入步骤S302；
S304、计算得到天然气流量分布情况。
本发明在另一方面，提出了基于分立求解的综合能源系统多能流联合计算装置，包括：
构建模块，用于构建电气耦合综合能源系统多能流模型；
解耦模块，用于使用耦合元件实现天然气系统和电力系统的解耦；
电力系统潮流模块，用于解耦后运行电力系统潮流计算；
天然气系统潮流模块，用于解耦后天然气系统潮流计算；
输出模块，用于得到各类状态变量作为多能流仿真的最终结果后进行输出。
进一步的，所述构建模块包括：
天然气子模块，用于建立天然气系统模型；
电力子模块，用于建立电力系统模型；
耦合元件子模块，用于建立耦合元件的模型；
综合子模块，用于建立电气耦合综合能源系统多能流模型；定义xeg和yeg为电气耦合综合能源系统IEGS的状态变量和注入变量；
xeg＝[V,θ,p]
yeg＝[Psp,Qsp,Lsp]
Feg(xeg,yeg)＝0
式中xeg由电力系统节点电压V、相角θ、天然气系统节点气压p组成；yeg由电力系统节点有功功率注入PSP、无功功率注入QSP、天然气系统节点注入LSP组成。
进一步的，所述解耦模块包括：
电力等效子模块，用于将所述耦合元件在电力系统中等效为一个有功出力已知的发电机；
天然气等效子模块，用于将所述耦合元件在天然气系统中等效为一个燃气需求已知的天然气负荷；
解耦子模块，用于实现天然气系统和电力系统的解耦。
进一步的，所述电力系统潮流模块包括：
矩阵子模块，用于读取电网拓扑，网络参数，形成节点导纳矩阵；
初始化子模块，用于初始化电力系统状态变量V(0)，θ(0)；
误差计算子模块，用于分别计算电力系统节点误差函数Fe(V，θ)；
雅可比矩阵计算子模块，用于计算电力系统节点雅可比矩阵Je；
迭代子模块，用于计算迭代修正量ΔV，Δθ；
判断子模块，用于更新电力系统状态变量，判断电力系统节点误差函数最大值是否小于阈值；若小于，迭代到此结束，进一步计算各线路潮流和平衡节点功率，反之进入误差计算子模块。
进一步的，所述天然气系统潮流模块包括：
关联矩阵子模块、用于根据天然气网络拓扑形成全节点-支路关联矩阵以及降阶节点 -支路关联矩阵，读取管道阻力系数，计算天然气系统各个节点的净负荷，初始化节点压力；
节点误差子模块，用于以节点压力平方为变量，求解天然气节点误差函数；
计算并判断子模块，用于计算新的节点压力平方，判断天然气节点误差函数最大值是否小于阈值；若小于，得到节点压力平方真值，进入分布子模块，否则进入步节点误差子模块；
分布子模块，用于计算得到天然气流量分布情况。
与现有技术相比，本发明具有如下的有益效果：
本发明提出基于分立迭代的IEGS多能流计算方法，通过能源耦合部分的能量等值将多能流待求信息分配到不同能源各自的求解模块，进行分立求解，本发明所采用的方法物理意义清晰明了，计算结果准确，该方法可有效计算并很好的模拟IEGS的运行状态，是进行IEGS安全稳定性分析的基础工作，为系统安全稳定性分析提供支撑工具。
附图说明
图1为本发明实施例所列举的一个简单的天然气系统；
图2为本发明实施例中电气耦合综合能源系统分立求解算法的流程图；
图3为本发明实施例中电气耦合综合能源系统IEGS的算例。
具体实施方式
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
为使本发明专利的目的、特征更明显易懂，下面结合附图对本发明专利的具体实施方式作进一步的说明。需说明的是，附图均采用非常简化的形式且均使用非精准的比率，仅用以方便、明晰地辅助说明本发明专利实施例的目的。
本发明提供的一种基于分立求解的电气耦合综合能源系统多能流联合计算方法包括按顺序进行的下列步骤：
步骤1)建立燃气系统与电力系统模型；
步骤2)提出燃气系统与电力系统潮流计算模型；
步骤3)一种基于分立求解的电气耦合综合能源系统多能流联合计算方法。
在步骤1)中，所建立燃气系统与电力系统模型具体步骤如下：
1.1)建立天然气系统模型：
天然气系统的变量为节点天然气注入量和节点压力。天然气系统中的节点分为两种类型。一个是天然气注入量已知的节点，另一个是压力已知的节点。假定天然气气源为平衡节点，平衡节点的压力已知，天然气注入量未知。天然气负荷节点的节点天然气注入量已知，节点气压未知。
天然气管道分为带压缩机的管道和不带压缩机的管道。本发明为了研究问题的简化，暂时不考虑带有压缩机的管道。对于不带压缩机的管道，天然气节点m到天然气节点n的管道流量fmn如式(1)、(2)所示。
式中cmn为管道阻力系数，与管道的粗糙度、直径和长度有关；pm为天然气节点m的气压；pn为天然气节点n的气压；smn为符号指示函数。
对天然气系统的任意节点，流入和流出的流量必须保持平衡，如式(3)、(4)所示。
式中
为天然气节点m的注入流量；m∈n表示天然气节点m和天然气节点n之间有管道连接；Ls,m表示天然气节点m的天然气供应量；Ll,m表示天然气节点m的天然气需求量。
1.2)建立电力系统模型：
电力系统的模型描述如式(5)-(6)所示。
式中
，
为电力系统节点i的注入有功功率和无功功率；Pg,i，Qg,i为电力系统节点i上发电机的有功功率，无功功率；Pl,i，Ql,i为电力系统节点i上负荷的有功功率，无功功率；Vi，Vj为电力系统节点i，j的电压；Gij，Bij为导纳Yij的实部，虚部；θij为节点i、j之间的相角差。
1.3)建立天然气发电机的模型：
设IEGS中的所有天然气发电机可以用集合ΩNGU表示为式(9)。
ΩNGU＝{GU1,GU2,···,GUN} (9)
式中GUi是第i个天然气发电机；N是天然气发电机的总数目。
系统中天然气发电机的有功注入向量W表示为式(10)。
W＝[PU,1,PU,2,···,PU,N] (10)
天然气发电机通过消耗天然气进行发电，其天然气消耗与发电量的关系表达为式(11)。
式中ai,bi,ci是第i个天然气发电机的消耗系数；L U,i为第i个天然气发电机的天然气消耗量。
1.4)建立电气耦合综合能源系统多能流模型：
定义xeg和yeg为IEGS的状态变量和注入变量。IEGS的多能流模型由式(1)～(11)组成，可以归纳为式(12)～(14)。
xeg＝[V,θ,p] (12)
yeg＝[Psp,Qsp,Lsp] (13)
Feg(xeg,yeg)＝0 (14)
式中xeg由电力系统节点电压、相角、天然气系统节点气压组成；yeg由电力系统节点有功功率注入、无功功率注入、天然气系统节点注入组成。
在步骤2)中，所述建立天然气和电力系统的潮流模型具体步骤如下：
2.1)建立天然气系统潮流计算模型
根据基尔霍夫第一定律，天然气系统中任何一个节点的注入流量代数和为零。在天然气系统中，任何一个节点的净负荷(负荷需求减去天然气供应)等于流入该节点的全部流量减去流出该节点的全部流量，即式(3)-(4)用矩阵的形式可以表示为式(15)。
A1f-L＝0 (15)
式中L为天然气管网中的燃气负荷矢量；A1为降阶节点-支路关联矩阵；f为支路流量矢量。
对如图1所示的一个简单的天然气系统，根据基尔霍夫第一定律，式(15)可以表达为式(16)。
A1矩阵中各元素符号规定：若支路气流流进该节点，则为正。若支路气流流出该节点，则为负。
支路两端的压力平方差用矩阵形式表示为式(17)、(18)。
ΔΠ＝-ATΠ (17)
式中△Π为天然气管网中的压力平方差；Π为节点压力平方向量；A为全节点-支路关联矩阵；Πm为天然气节点m的气压平方。
对图1所示的天然气系统，式(17)可以表达为式(19)。
根据式(17)和(18)，式(1)和(2)可以改写为式(20)和(21)。
式中i为支路编号。
在计算天然气系统潮流的时候，以气压的平方作为未知数，式(20)可以表达为式(22)
将式(22)代入式(15)得式(23)。
通过牛顿节点法求解式(23)，即可求得天然气网络的节点气压分布，再通过式(22) 计算得到天然气流量分布。首先要给出各节点压力的初始估计值，然后通过连续修正估计值直到获得最终结果。式(23)的左边在迭代过程中不等于零，其中节点压力值只是它们真值的估计值。设Fg()为天然气节点误差函数，它是所有节点压力平方的函数。除了参考节点外，整组误差可表示成
当压力平方值接近其真值时，天然气节点误差函数将会趋向于零。牛顿节点法就是求解公式组成的方程组，直到误差小于给定阈值。用于修正节点压力平方估计值的迭代过程如式(25)。
Π(k+1)＝Π(k)+(δΠ)(k) (25)
式中Π(k)为第k次迭代中的节点压力平方，(δΠ)(k)为第k次迭代中的节点压力增量。
δΠ的计算过程如式(26)。
Jg是节点雅可比矩阵，表达式如式(27)-(28)。
电力系统的状态变量初始化遵循的是“平直电压法”，即将PQ节点和平衡节点的电压设置为1，所有非平衡节点的节点相角均设为0。然而，相比于电力系统的电压相角设置，天然气系统变量的初始化需要注意初始天然气气压的设置，否则会因为潮流发散而使整个潮流计算失败。本发明的天然气节点压力设置按照以下方法进行：以气源节点为参考节点，按顺序将以后每个节点的压力减少2％进行初始化。
综上，天然气系统潮流求解过程可以描述为：
步骤一：根据天然气网络拓扑形成全节点-支路关联矩阵以及降阶节点-支路关联矩阵，读取管道阻力系数，计算天然气系统各个节点的净负荷，初始化节点压力。
步骤二：以节点压力平方为变量，代入式(24)求解天然气节点误差函数。
步骤三：根据式(25)计算新的节点压力平方，判断天然气节点误差函数最大值是否小于阈值。若小于，得到节点压力平方真值，进入步骤四，否则进入步骤二。
步骤四：通过式(22)计算得到天然气流量分布情况。
2.2)建立电力系统潮流计算模型
通过Newton-Raphson法计算电力系统潮流时，有以下的步骤：
步骤一：读取电网拓扑，网络参数，形成节点导纳矩阵。
步骤二：初始化电力系统状态变量V(0)，θ(0)。
步骤三：根据式(29)分别计算电力系统节点误差函数Fe(V，θ)。
步骤四：根据式(30)计算电力系统节点雅可比矩阵Je。
步骤五：根据式(31)计算迭代修正量ΔV，Δθ。
步骤六：更新电力系统状态变量，判断电力系统节点误差函数最大值是否小于阈值。若小于，迭代到此结束，进一步计算各线路潮流和平衡节点功率，反之进入步骤三。
max|Fe(V,θ)|＜ε (32)
在步骤3)中，所述的基于分立求解的电气耦合综合能源系统多能流联合计算方法, 并生成多能流计算结果的具体步骤如下：
本发明电气两系统的耦合环节为天然气发电机。天然气发电机在电力系统中往往作为PV节点处理，其有功功率是已知的，同时，天然气发电机在天然气系统中作为负荷，可以根据其输出的有功功率按式(11)计算其天然气消耗。此时天然气发电机在电力系统中等效为一个有功出力已知的发电机，在天然气系统中等效为一个燃气需求已知的天然气负荷，由此实现天然气系统和电力系统的解耦。解耦后，分别运行电力系统潮流计算以及天然气系统潮流计算，所得各类状态变量就是多能流仿真的最终结果，方法流程图如图2所示。
解耦模块：本模块用于将天然气系统潮流计算和电力系统系统潮流计算模块解耦。首先读取天然气发电机的有功出力信息，接着按式(11)计算其天然气消耗，最后更新天然气系统负荷，即对与天然气发电机相连接的天然气负荷节点，其天然气需求等于本地天然气需求加上天然气发电机的燃气需求。
电力系统潮流模块：本模块通过Newton-Raphson法计算电力系统潮流分布。首先输入节点和支路的相关信息，形成节点导纳矩阵；接着给定电压相角的初始值；然后依次计算电力系统雅克比矩阵，并更新节点电压相角值，不断迭代直到电力系统节点误差函数小于给定阈值；最后得到电力系统的状态变量。
天然气系统潮流模块：本模块通过牛顿节点法计算天然气系统潮流分布。首先根据节点和管道数据，形成支路-节点关联矩阵；接着给定气压的初始值；然后以节点压力平方为变量，依次计算天然气系统雅克比矩阵，并更新节点气压平方，不断迭代直到天然气系统节点误差函数小于给定阈值；最后得到天然气系统的状态变量。
输出模块：本模块用于输出分立求解算法中电-气耦合综合能源系统的状态变量，主要包括电力系统中所有母线节点的电压幅值和相角，线路的功率分布、功率损耗；天然气系统中各个节点的气压幅值，管道的流量分布。
应用例
本发明的应用算例由WSCC9节点电力系统和修改的6节点燃气网络通过2个天然气发电机耦合而成，EBi和GBi分别表示电网节点和燃气节点，其拓扑如图3所示。对于 WSCC系统而言，EB1为平衡节点，EB2和EB3为天然气发电机G2和G3所在的节点，为 PV节点。G2和G3连接天然气系统的GB6和GB3，两天然气发电机参数相同。天然气发电机的参数如表1所示，天然气网络的参数如表2所示。
表1天然气发电机的参数
表2天然气网络参数
首先，根据天然气发电机的有功出力计算其天然气消耗，如表3所示。
表3天然气发电机的有功出力以及天然气消耗
考虑天然气发电机的燃气需求，更新天然气网络的负荷如表4所示。至此，天然气系统的网络以及负荷需求分布就确定了。
表4天然气网络的负荷(计及天然气发电机燃气需求)
将多能流待求信息分配到不同能源的各自的求解模块。通过分立潮流计算，分别运行电力系统潮流计算以及天然气系统潮流计算程序，获得天然气网络的气压分布以及电网的节点电压相角等关键信息。
表5～表8给出了本算例的计算结果：天然气网络的气压分布如表5所示，天然气网络的流量分布如表6所示。电力系统的节点电压相角如表7所示，支路计算结果如表8 所示。
表5天然气网络的气压分布
表6天然气网络的流量分布
表7电力网络的节点电压相角
(注：有功负荷注入以向节点注入功率为正，因此发电机有功注入功率为正数，负荷有功注入功率为负数。无功注入同理，发电机作为无功电源向网络提供功率，负荷吸收无功功率，因此，发电机无功注入功率为正数，负荷无功注入功率为负数。联络节点无功率注入。)
表8电力网络的支路计算结果
可以看出，本发明所采用的方法物理意义清晰明了，计算结果准确，可以很好的模拟系统的运行特性，是进行IEGS安全稳定性分析的基础工作。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.基于分立求解的综合能源系统多能流联合计算方法，其特征在于，包括：构建电气耦合综合能源系统多能流模型；使用耦合元件实现天然气系统和电力系统的解耦；解耦后分别运行电力系统潮流计算以及天然气系统潮流计算，得到各类状态变量作为多能流仿真的最终结果。
2.根据权利要求1所述的基于分立求解的综合能源系统多能流联合计算方法，其特征在于，所述电气耦合综合能源系统多能流模型的构建过程为：
S101、建立天然气系统模型；
S102、建立电力系统模型；
S103、建立耦合元件的模型；
S104、建立电气耦合综合能源系统多能流模型，定义xeg和yeg为电气耦合综合能源系统IEGS的状态变量和注入变量；
xeg＝[V,θ,p]
yeg＝[Psp,Qsp,Lsp]
Feg(xeg,yeg)＝0
式中xeg由电力系统节点电压V、相角θ、天然气系统节点气压p组成；yeg由电力系统节点有功功率注入PSP、无功功率注入QSP、天然气系统节点注入LSP组成。
3.根据权利要求1所述的基于分立求解的综合能源系统多能流联合计算方法，其特征在于，所述耦合元件的解耦过程具体为：
将所述耦合元件在电力系统中等效为一个有功出力已知的发电机，在天然气系统中等效为一个燃气需求已知的天然气负荷，实现天然气系统和电力系统的解耦。
4.根据权利要求1所述的基于分立求解的综合能源系统多能流联合计算方法，其特征在于，所述电力系统潮流计算的方法包括：
S201、读取电网拓扑，网络参数，形成节点导纳矩阵；
S202、初始化电力系统状态变量V(0)，θ(0)；
S203、分别计算电力系统节点误差函数Fe(V，θ)；
S204、计算电力系统节点雅可比矩阵Je；
S205、计算迭代修正量ΔV，Δθ；
S206、更新电力系统状态变量，判断电力系统节点误差函数最大值是否小于阈值；若小于，迭代到此结束，进一步计算各线路潮流和平衡节点功率，反之进入步骤S203。
5.根据权利要求1所述的基于分立求解的综合能源系统多能流联合计算方法，其特征在于，所述天然气系统潮流计算的具体方法为：
S301、根据天然气网络拓扑形成全节点-支路关联矩阵以及降阶节点-支路关联矩阵，读取管道阻力系数，计算天然气系统各个节点的净负荷，初始化节点压力；
S302、以节点压力平方为变量，求解天然气节点误差函数；
S303、计算新的节点压力平方，判断天然气节点误差函数最大值是否小于阈值；若小于，得到节点压力平方真值，进入步骤S304，否则进入步骤S302；
S304、计算得到天然气流量分布情况。
6.基于分立求解的综合能源系统多能流联合计算装置，其特征在于，包括：
构建模块，用于构建电气耦合综合能源系统多能流模型；
解耦模块，用于使用耦合元件实现天然气系统和电力系统的解耦；
电力系统潮流模块，用于解耦后运行电力系统潮流计算；
天然气系统潮流模块，用于解耦后天然气系统潮流计算；
输出模块，用于得到各类状态变量作为多能流仿真的最终结果后进行输出。
7.根据权利要求6所述的基于分立求解的综合能源系统多能流联合计算装置，其特征在于，所述构建模块包括：
天然气子模块，用于建立天然气系统模型；
电力子模块，用于建立电力系统模型；
耦合元件子模块，用于建立耦合元件的模型；
综合子模块，用于建立电气耦合综合能源系统多能流模型；定义xeg和yeg为电气耦合综合能源系统IEGS的状态变量和注入变量；
xeg＝[V,θ,p]
yeg＝[Psp,Qsp,Lsp]
Feg(xeg,yeg)＝0
式中xeg由电力系统节点电压V、相角θ、天然气系统节点气压p组成；yeg由电力系统节点有功功率注入PSP、无功功率注入QSP、天然气系统节点注入LSP组成。
8.根据权利要求6所述的基于分立求解的综合能源系统多能流联合计算装置，其特征在于，所述解耦模块包括：
电力等效子模块，用于将所述耦合元件在电力系统中等效为一个有功出力已知的发电机；
天然气等效子模块，用于将所述耦合元件在天然气系统中等效为一个燃气需求已知的天然气负荷；
解耦子模块，用于实现天然气系统和电力系统的解耦。
9.根据权利要求6所述的基于分立求解的综合能源系统多能流联合计算装置，其特征在于，所述电力系统潮流模块包括：
矩阵子模块，用于读取电网拓扑，网络参数，形成节点导纳矩阵；
初始化子模块，用于初始化电力系统状态变量V(0)，θ(0)；
误差计算子模块，用于分别计算电力系统节点误差函数Fe(V，θ)；
雅可比矩阵计算子模块，用于计算电力系统节点雅可比矩阵Je；
迭代子模块，用于计算迭代修正量ΔV，Δθ；
判断子模块，用于更新电力系统状态变量，判断电力系统节点误差函数最大值是否小于阈值；若小于，迭代到此结束，进一步计算各线路潮流和平衡节点功率，反之进入误差计算子模块。
10.根据权利要求6所述的基于分立求解的综合能源系统多能流联合计算装置，其特征在于，所述天然气系统潮流模块包括：
关联矩阵子模块、用于根据天然气网络拓扑形成全节点-支路关联矩阵以及降阶节点-支路关联矩阵，读取管道阻力系数，计算天然气系统各个节点的净负荷，初始化节点压力；
节点误差子模块，用于以节点压力平方为变量，求解天然气节点误差函数；
计算并判断子模块，用于计算新的节点压力平方，判断天然气节点误差函数最大值是否小于阈值；若小于，得到节点压力平方真值，进入分布子模块，否则进入步节点误差子模块；
分布子模块，用于计算得到天然气流量分布情况。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010395994.7A CN111523247A (zh) | 2020-05-12 | 2020-05-12 | 基于分立求解的综合能源系统多能流联合计算方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010395994.7A CN111523247A (zh) | 2020-05-12 | 2020-05-12 | 基于分立求解的综合能源系统多能流联合计算方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111523247A true CN111523247A (zh) | 2020-08-11 |
Family
ID=71907056
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010395994.7A Pending CN111523247A (zh) | 2020-05-12 | 2020-05-12 | 基于分立求解的综合能源系统多能流联合计算方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111523247A (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111767660A (zh) * | 2020-09-02 | 2020-10-13 | 国网江西省电力有限公司电力科学研究院 | 一种电气耦合能源系统的储气装置优化布置方法及装置 |
| CN112818492A (zh) * | 2021-02-02 | 2021-05-18 | 山东大学 | 计及能源集站状态变量的电-气耦合网络能量流求解方法 |
| CN112994020A (zh) * | 2021-03-31 | 2021-06-18 | 南京信息工程大学 | 一种多能源系统潮流分解计算方法 |
| CN113887051A (zh) * | 2021-10-09 | 2022-01-04 | 清鸾科技(成都)有限公司 | 基于层次模型的综合能源系统仿真计算方法、装置、综合能源系统和计算机可读存储介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110277785A (zh) * | 2019-06-26 | 2019-09-24 | 国网浙江省电力有限公司电力科学研究院 | 基于连续多能流的电气耦合系统负荷裕度计算方法及系统 |
| CN110543695A (zh) * | 2019-08-14 | 2019-12-06 | 天津大学 | 一种电-气耦合综合能源系统可行域计算方法 |
| CN110647040A (zh) * | 2019-10-10 | 2020-01-03 | 东北电力大学 | 一种综合能源系统的安全控制方法及装置 |
-
2020
- 2020-05-12 CN CN202010395994.7A patent/CN111523247A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110277785A (zh) * | 2019-06-26 | 2019-09-24 | 国网浙江省电力有限公司电力科学研究院 | 基于连续多能流的电气耦合系统负荷裕度计算方法及系统 |
| CN110543695A (zh) * | 2019-08-14 | 2019-12-06 | 天津大学 | 一种电-气耦合综合能源系统可行域计算方法 |
| CN110647040A (zh) * | 2019-10-10 | 2020-01-03 | 东北电力大学 | 一种综合能源系统的安全控制方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| 王锡凡: "《电力工程基础》", 31 January 1998, pages: 47 - 52 * |
| 黎静华: "综合能源系统多能流潮流计算模型与方法综述", pages 1 - 9 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111767660A (zh) * | 2020-09-02 | 2020-10-13 | 国网江西省电力有限公司电力科学研究院 | 一种电气耦合能源系统的储气装置优化布置方法及装置 |
| CN112818492A (zh) * | 2021-02-02 | 2021-05-18 | 山东大学 | 计及能源集站状态变量的电-气耦合网络能量流求解方法 |
| CN112994020A (zh) * | 2021-03-31 | 2021-06-18 | 南京信息工程大学 | 一种多能源系统潮流分解计算方法 |
| CN113887051A (zh) * | 2021-10-09 | 2022-01-04 | 清鸾科技(成都)有限公司 | 基于层次模型的综合能源系统仿真计算方法、装置、综合能源系统和计算机可读存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111523247A (zh) | 基于分立求解的综合能源系统多能流联合计算方法及装置 | |
| CN104734155B (zh) | 一种获取电‑气互联能源系统可用输电能力的方法 | |
| CN111061986B (zh) | 一种多运行模式的热电综合能源系统潮流计算方法 | |
| Eltamaly et al. | Load flow analysis by gauss-seidel method; a survey | |
| CN105005940A (zh) | 计及相关性的电-气互联系统概率最优潮流计算方法 | |
| CN111428351B (zh) | 基于前推回代法的电-热综合能源系统潮流计算方法 | |
| CN105243516A (zh) | 基于主动配电网的分布式光伏发电最大消纳能力计算系统 | |
| CN102545252B (zh) | 基于三阶收敛牛顿法vsc-hvdc潮流计算方法 | |
| CN108875269A (zh) | 一种考虑电力系统多平衡机与天然气系统慢动态特性的电-气互联系统多时段能流计算方法 | |
| Zhao et al. | Distributed continuation power flow method for integrated transmission and active distribution network | |
| CN112531716A (zh) | 电-水互联系统混合能流的统一标幺计算方法 | |
| CN110532642A (zh) | 一种综合能源系统概率能流的计算方法 | |
| Zhao et al. | Holomorphic embedding power flow for AC/DC hybrid power systems using Bauer's eta algorithm | |
| CN114221346A (zh) | 一种综合能源系统潮流计算方法 | |
| CN113131513B (zh) | 考虑碳排放的电、热、气转换优化运行配置方法和存储介质 | |
| Chen et al. | Multi-stage dynamic optimal power flow in wind power integrated system | |
| Wang et al. | Optimal gas-power flow calculation with lightweight gas flow dynamics | |
| Tian et al. | Modeling and simulation for multi energy flow coupled network computing | |
| CN113627021A (zh) | 一种基于序列凸规划的电气互联系统最优能流计算方法 | |
| CN106253288A (zh) | 一种基于自动微分的含统一潮流控制器的最优潮流算法 | |
| CN103199524B (zh) | 一种适应多种分布式电源接入的潮流计算方法 | |
| Aly et al. | Voltage stability assessment for radial distribution power system with wind power penetration | |
| CN111310310A (zh) | 一种用于量调节的热力系统静态潮流快速解耦计算方法 | |
| CN117236590A (zh) | 综合能源系统运行成本和碳排放的协同优化方法及系统 | |
| CN112528445B (zh) | 一种电-气互联综合能源系统协同规划方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |