CN115204549A - 一种区域电气综合能源系统小扰动稳定性评估方法 - Google Patents
一种区域电气综合能源系统小扰动稳定性评估方法 Download PDFInfo
- Publication number
- CN115204549A CN115204549A CN202210235874.XA CN202210235874A CN115204549A CN 115204549 A CN115204549 A CN 115204549A CN 202210235874 A CN202210235874 A CN 202210235874A CN 115204549 A CN115204549 A CN 115204549A
- Authority
- CN
- China
- Prior art keywords
- gas
- stability
- comprehensive energy
- variable
- state
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 16
- 238000004088 simulation Methods 0.000 claims abstract description 38
- 239000007789 gas Substances 0.000 claims description 71
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 62
- 239000003345 natural gas Substances 0.000 claims description 31
- 238000002485 combustion reaction Methods 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 17
- 238000000034 method Methods 0.000 claims description 16
- 230000010355 oscillation Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000000446 fuel Substances 0.000 claims description 4
- 230000003993 interaction Effects 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 7
- 230000006641 stabilisation Effects 0.000 abstract description 2
- 238000011105 stabilization Methods 0.000 abstract description 2
- 238000012795 verification Methods 0.000 abstract description 2
- 230000008878 coupling Effects 0.000 description 6
- 238000010168 coupling process Methods 0.000 description 6
- 238000005859 coupling reaction Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000087 stabilizing effect Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004540 process dynamic Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- 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)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Public Health (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种电‑气综合能源系统小扰动稳定性评估方法，建立了考虑燃机压力控制和功率调整的多燃机区域电‑气稳定性分析模型，基于李雅普诺夫稳定理论有效评估区域电‑气综合能源系统小扰动稳定性，并能够通过时序仿真来进行对比验证。采用以上技术方案与现有技术相比，能够基于该稳定性评估方法对系统稳定性影响因素进行对比分析，具有实践意义。
Description
技术领域
本发明涉及一种电-气综合能源系统小扰动稳定性评估方法。特别是涉及一种区域电气综合能源系统小扰动稳定性评估方法。
背景技术
在能源转型和能效提升的时代背景下，燃气发电机因其快速灵活的调节特性和低碳排的优势，广泛应用于平抑可再生能源波动和提高区域供能效率等场景，成为未来一段时间内能源系统中不可或缺的重要环节之一。然而，燃气发电机的大量应用将导致电力系统和天然气系统的耦合程度逐步加深，电力系统或天然气系统的扰动会通过耦合设备影响另一个系统安全稳定运行，有可能引发跨系统的严重故障。电-气两种能源系统相互耦合影响及安全稳定运行成为亟待研究的问题。
现有电-气综合能源系统稳定性的研究，主要根据天然气网络模型的不同而分为三类，即不考虑天然气网络、天然气网络稳态模型和动态模型。其中，不考虑天然气网络，即针对燃气发电机和电力系统构成的系统展开研究，主要包括燃气发电机容量、位置、接入数量等因素对电力系统电压稳定性的影响。但在实际运行场景中，燃气发电机的耦合作用使电力系统的稳定性在一定程度上还要依赖天然气系统的运行稳定性，因此还需考虑天然气网络的影响。基于天然气网络稳态模型的稳定性分析中，多是参照电力系统安全域针对电-气综合能源系统提出其安全域的概念并加以分析，但稳态模型假设条件是天然气网络某节点状态发生变化会立刻反映到全网中,忽略了天然气网络固有的慢变过程动态特性，在诸多应用场景中并不适用。基于天然气网络动态模型的研究中，针对考虑连接单根天然气管道燃气发电机的电-气综合能源系统展开分析，但并没有考虑天然气网络对系统稳定性的影响。目前已有对电-气综合能源系统的安全稳定问题的研究中，大多是从耦合设备或天然气系统的角度分析电力系统的稳定性，且多为电力和天然气的输送系统，而很少有从天然气的角度对区域配电、配气耦合系统稳定性进行研究，但在电-气互联下，天然气系统稳定性的破坏不仅会影响居民、工业用气负荷，同时也会通过燃气发电机等耦合设备传播，进而影响电力系统安全稳定运行。因此，亟需发展区域电-气综合能源系统小扰动稳定性评估方法。
发明内容
本发明所要解决的技术问题是，提出了一种适用于计及微燃机的区域电-气综合能源系统稳定性分析方法，用于评估区域电-气综合能源系统小扰动稳定性。
本发明所采用的技术方案是：
步骤1：输入系统参数(天然气系统拓扑结构，管道长度、内径与摩擦系数，天然气负荷，天然气气源压力，微燃机额定功率，微燃机额定流量，阀门PI控制参数)。
步骤2：确定空间步长，进行管道空间维度有限差分。
步骤3：根据系统拓扑结构、负荷位置及类型和空间离散结果确定系统状态变量x、控制变量和中间变量u，并建立区域电-气综合能源系统非线性微分代数方程组。
步骤4：整理线性代数方程组，用x表示u并带入到非线性微分方程组中，消去控制变量和中间变量u，得到形如
只含有状态变量x的系统状态方程。
步骤5：根据给定数据求解非线性微分代数方程得到系统运行点(x0,u0)，并在运行点处对系统状态方程进行线性化，得到系统状态矩阵。
步骤6：计算系统状态矩阵特征值，并根据其主导特征值在复平面的分布情况对区域电- 气综合能源系统小扰动稳定性进行评估。同时可以根据某一参数的变化情况绘制其根轨迹或稳定域情况，分析参数对区域电-气综合能源系统小扰动稳定性的影响进行评估。
步骤7：基于步骤3中建立的区域电-气综合能源系统非线性微分代数方程组，即时域仿真模型，设置仿真参数进行时域仿真后可将系统仿真结果的收敛、振荡及发散情况与步骤6 中的稳定性评估结果进行对比，验证区域电-气综合能源系统小扰动稳定性评估方法的正确性和有效性。
上述技术方案的进一步特征在于，所述步骤3中根据系统拓扑结构、负荷位置及类型和空间离散结果确定系统状态变量x、控制变量和中间变量u的方法如下：
由于在区域配电系统中，微燃机主要以接受区域配电调度或平抑可再生能源波动模式运行，电力负荷波动由上级电网平抑，配电系统的拓扑结构和负荷变化对系统稳定性的影响不大。因此在区域电-气综合能源系统小扰动稳定性评估过程中主要考虑微燃机和配气系统两部分的影响。
步骤3-1：针对空间维度有限差分后的配气系统，选取状态变量xGas由管道末端压力pout和首端质量流量Min构成，选取控制变量及中间变量uGas由管道首端压力pin和末端质量流量 Mout构成；
步骤3-2：针对计及天然气系统相互作用及动态特性的微燃机，选取状态变量xMT由转速 w、燃料流量Wf等构成，选取控制变量及中间变量uMT由参考流量Fref、实际流量Fmt等构成；
步骤3-3：系统状态变量由xGas和xMT两部分构成，系统控制变量及中间变量由uGas和uMT两部分构成。
上述技术方案的进一步特征在于，所述步骤5中在运行点处对系统状态方程进行线性化得到系统状态矩阵的方法如下：
步骤5-1：将系统状态方程中状态变量微分项置零，联立计算非线性代数方程组得到系统运行点；
步骤5-2：在状态方程的基础上将状态变量写作增量形式Δx，并将状态方程Fi依次对状态变量xj进行求导，其结果为状态矩阵中对应位置Aij元素构成；
步骤5-3：将运行点(x0,u0)带入并依次计算Aij的值，得到状态矩阵A。
上述技术方案的进一步特征在于，所述步骤7中区域电-气综合能源系统时域仿真并根据仿真结果与小扰动稳定性评估结果进行对比的方法如下：
步骤7-1：设置变量初值、系统扰动及仿真参数(仿真时长，仿真步长，收敛精度等)；
步骤7-2：选择常微分方程数值解法进行迭代求解，如各阶龙格库塔法(Runge-Kutta)等；
步骤7-3：根据仿真结果及步骤6中系统特征值的计算结果进行对比，可对系统小扰动稳定性进行评估，具体分为以下三种情况：
1)若所有的特征值都处在复平面虚轴左侧，则运行点(x0,u0)处系统是稳定运行的，且系统仿真结果收敛；
2)若至少一个特征值处在复平面虚轴右侧，则运行点(x0,u0)处系统运行是不稳定的，且系统仿真结果发散；
3)若存在特征值在处在复平面虚轴上，而其他特征值都在虚轴左侧，则系统是临界稳定的，且系统仿真结果呈现等幅振荡的趋势。
本发明公开一种区域电气综合能源系统小扰动稳定性评估方法，采用以上技术方案与现有技术相比，具有以下技术效果：本发明所设计计及微燃机的区域电-气综合能源系统小扰动稳定性评估方法，建立了考虑燃机压力控制和功率调整的多燃机区域电-气稳定性分析模型，基于李雅普诺夫稳定理论有效评估区域电-气综合能源系统小扰动稳定性，并能够通过时序仿真来进行对比验证。同时，还能够基于该稳定性评估方法对系统稳定性影响因素进行对比分析，具有实践意义。
附图说明
图1是本发明计及微燃机的区域电-气综合能源系统小扰动稳定性评估方法流程图；
图2是双机三节点配气系统算例结构图；
图3是微燃机控制环节比例系数一定、积分系数逐渐增大时系统主导特征值根轨迹图；
图4是微燃机控制参数组合在不同情况下系统主导特征值实部的数值情况；
图5是微燃机控制参数按表1所示进行取值时G3节点压力时域仿真结果；
图6是微燃机控制参数按表1所示进行取值时MT3燃机流量时域仿真结果；
具体实施方式
为了使本发明的目的、技术方案及有点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
如图1所示，本发明所述一种计及微燃机的区域电-气综合能源系统小扰动稳定性评估方法，详细说明，该方法包括以下步骤：
步骤1：输入系统参数(天然气系统拓扑结构，管道长度、内径与摩擦系数，天然气负荷，天然气气源压力，微燃机额定功率，微燃机额定流量，阀门PI控制参数)，本发明采用双机三节点天系统作为分析算例，如图2所示。微燃机采用Rowen简化模型，并对其中的燃料供应环节加以改进，天然气管道采用如式(1)所示动态模型：
式中，x和t分别代表空间位置和时间；M为质量流量；A为管道横截面积；p为气体压力；d为管道直径；λ为摩擦系数；c为天然气中的声速。
步骤2：确定空间步长Δx，进行管道空间维度有限差分，差分格式如式(2)所示：
以空间步长Δx进行空间差分后有：
式中，针对离散后的某一段管道，pin为管道首端压力；pout为管道末端压力；Min为管道首端质量流量；Mout为管道末端质量流量。
步骤3：根据系统拓扑结构、负荷位置及类型和空间离散结果确定系统状态变量x、控制变量和中间变量u:
步骤3-1：针对空间维度有限差分后的配气系统，选取状态变量xGas由管道末端压力pout和首端质量流量Min构成，选取控制变量及中间变量uGas由管道首端压力pin和末端质量流量 Mout构成，即xGas＝[pout,Min]T，uGas＝[pin,Mout]T；
步骤3-2：针对计及天然气系统相互作用及动态特性的微燃机，选取状态变量xMT由转速 w、燃料流量Wf等构成，选取控制变量及中间变量uMT由参考流量Fref、实际流量Fmt等构成，即xMT＝[w,Tmpv,Wf]T，uMT＝[Fref,Fmt]T；
步骤3-3：系统状态变量由xGas和xMT两部分构成，系统控制变量及中间变量由uGas和uMT两部分构成，即x＝[xGas,xMT]T，u＝[uGas,uMT]T；
建立区域电-气综合能源系统非线性微分代数方程组如式(4)所示：
步骤4：整理线性代数方程组0＝g(x,u)，用x表示u可得到u＝G(x)的形式，将其带入到非线性微分方程组中，消去控制变量和中间变量u，得到只含有状态变量x的系统状态方程如式(5)所示。
步骤5：根据给定数据求解非线性微分代数方程得到系统运行点(x0,u0)，并在运行点处对系统状态方程进行线性化，得到系统状态矩阵。
步骤5-1：将系统状态方程中状态变量微分项置零，即令
联立计算非线性代数方程组0＝F(x)得到系统运行点(x0,u0)；
步骤5-2：在状态方程的基础上将状态变量写作增量形式Δx并在运行点处对式(5)进行线性化，即可得到式(6)所示形式；
式中，A(x)为系统在运行点处的状态矩阵，且有
步骤5-3：将运行点(x0,u0)带入并依次计算Aij(x)的值，得到状态矩阵A。
步骤6：计算系统状态矩阵A特征值，并根据其主导特征值在复平面的分布情况对区域电-气综合能源系统小扰动稳定性进行评估。同时可以根据某一参数的变化情况绘制其根轨迹或稳定域情况，分析参数对区域电-气综合能源系统小扰动稳定性的影响进行评估。
特别的，当两台微燃机PI控制参数分别按照表1中三组参数进行设置时，其主导特征值实部分别为0.000226、0、-0.009815，分别对应系统不稳定、系统临界稳定及系统稳定。
表1三组微型燃气轮机PI控制参数
在此基础上，令两台燃机比例系数均为0.003，积分系数由0.0001取至0.1，分别计算状态矩阵的特征值后，可将其主导特征值绘制为如图3所示的根轨迹，可以看出随着积分系数的逐渐增大，根轨迹反复、多次穿越虚轴，系统稳定性呈现非线性特征。
类似的，同时将PI控制参数组合范围扩大，并保持MT2与MT3控制参数一致，比例系数由0.0001取至0.02，积分系数由0.0001取至0.1，分别计算不同控制参数组合下系统状态矩阵特征值后，可将其主导特征值实部绘制为如图4所示形式，其中黑线所示为主导特征值为0时对应的控制参数组合，此时系统临界稳定，黑线左侧黄绿色区域，系统主导特征值实部均大于零，系统是不稳定的；黑线右侧蓝色区域，系统主导特征值实部均小于零，系统是稳定的。
步骤7：基于区域电-气综合能源系统时域仿真模型，设置仿真参数进行时域仿真后可将系统仿真结果的收敛、振荡及发散情况与步骤6中的稳定性评估结果进行对比，验证区域电- 气综合能源系统小扰动稳定性评估方法的正确性和有效性。
步骤7-1：设置变量初值、系统扰动及仿真参数(仿真时长，仿真步长，收敛精度等)；
步骤7-2：选择常微分方程数值解法进行迭代求解，此处采用4阶龙格库塔法进行计算；
步骤7-3：根据仿真结果及步骤6中系统特征值的计算结果进行对比，可对系统小扰动稳定性进行评估。
如当两台微燃机PI控制参数分别按照表1中三组参数进行设置时，其G3节点压力、MT3 流量的时域仿真结果如图5、6所示，积分系数取值0.0007、0.000757、0.0008显示为系统发散、系统等幅振荡及系统收敛三种情况，分别对应其特征值为0.000226的系统不稳定、特征值为0的系统临界稳定及特征值为-0.009815的系统稳定，可说明区域电-气综合能源系统小扰动稳定性评估方法的正确性和有效性。
最后应当说明的是：上述仅为本发明的较佳实施例及所运用技术原理，并不用以限制本发明。本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，对于本发明做出的改进和修改都应在本发明的保护范围之内。
Claims (4)
1.一种区域电气综合能源系统小扰动稳定性评估方法，包括如下步骤：
步骤1：输入系统参数(天然气系统拓扑结构，管道长度、内径与摩擦系数，天然气负荷，天然气气源压力，微燃机额定功率，微燃机额定流量，阀门PI控制参数)。
步骤2：确定空间步长，进行管道空间维度有限差分。
步骤3：根据系统拓扑结构、负荷位置及类型和空间离散结果确定系统状态变量x、控制变量和中间变量u，并建立区域电-气综合能源系统非线性微分代数方程组。
步骤4：整理线性代数方程组，用x表示u并带入到非线性微分方程组中，消去控制变量和中间变量u，得到形如
只含有状态变量x的系统状态方程。
步骤5：根据给定数据求解非线性微分代数方程得到系统运行点(x0,u0)，并在运行点处对系统状态方程进行线性化，得到系统状态矩阵。
步骤6：计算系统状态矩阵特征值，并根据其主导特征值在复平面的分布情况对区域电-气综合能源系统小扰动稳定性进行评估。同时可以根据某一参数的变化情况绘制其根轨迹或稳定域情况，分析参数对区域电-气综合能源系统小扰动稳定性的影响进行评估。
步骤7：基于步骤3中建立的区域电-气综合能源系统非线性微分代数方程组，即时域仿真模型，设置仿真参数进行时域仿真后可将系统仿真结果的收敛、振荡及发散情况与步骤6中的稳定性评估结果进行对比，验证区域电-气综合能源系统小扰动稳定性评估方法的正确性和有效性。
2.根据权利要求1所述的一种区域电气综合能源系统小扰动稳定性评估方法，其特征在于，所述的步骤3中的根据系统拓扑结构、负荷位置及类型和空间离散结果确定系统状态变量x、控制变量和中间变量u的具体步骤如下：
步骤3-1：针对空间维度有限差分后的配气系统，选取状态变量xGas由管道末端压力pout和首端质量流量Min构成，选取控制变量及中间变量uGas由管道首端压力pin和末端质量流量Mout构成；
步骤3-2：针对计及天然气系统相互作用及动态特性的微燃机，选取状态变量xMT由转速w、燃料流量Wf等构成，选取控制变量及中间变量uMT由参考流量Fref、实际流量Fmt等构成；
步骤3-3：系统状态变量由xGas和xMT两部分构成，系统控制变量及中间变量由uGas和uMT两部分构成。
3.根据权利要求1所述的一种区域电气综合能源系统小扰动稳定性评估方法，其特征在于，所述的步骤5中的在运行点处对系统状态方程进行线性化得到系统状态矩阵的具体步骤如下：
步骤5-1：将系统状态方程中状态变量微分项置零，联立计算非线性代数方程组得到系统运行点；
步骤5-2：在状态方程的基础上将状态变量写作增量形式Δx，并将状态方程Fi依次对状态变量xj进行求导，其结果为状态矩阵中对应位置Aij元素构成；
步骤5-3：将运行点(x0,u0)带入并依次计算Aij的值，得到状态矩阵A。
4.根据权利要求1所述的一种区域电气综合能源系统小扰动稳定性评估方法，其特征在于，所述的步骤7中的区域电-气综合能源系统时域仿真并根据仿真结果与小扰动稳定性评估结果进行对比的具体步骤如下：
步骤7-1：设置变量初值、系统扰动及仿真参数(仿真时长，仿真步长，收敛精度等)；
步骤7-2：选择常微分方程数值解法进行迭代求解，如各阶龙格库塔法(Runge-Kutta)等；
步骤7-3：根据仿真结果及步骤6中系统特征值的计算结果进行对比，可对系统小扰动稳定性进行评估，具体分为以下三种情况：
1)若所有的特征值都处在复平面虚轴左侧，则运行点(x0,u0)处系统是稳定运行的，且系统仿真结果收敛；
2)若至少一个特征值处在复平面虚轴右侧，则运行点(x0,u0)处系统运行是不稳定的，且系统仿真结果发散；
3)若存在特征值在处在复平面虚轴上，而其他特征值都在虚轴左侧，则系统是临界稳定的，且系统仿真结果呈现等幅振荡的趋势。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210235874.XA CN115204549A (zh) | 2022-03-11 | 2022-03-11 | 一种区域电气综合能源系统小扰动稳定性评估方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210235874.XA CN115204549A (zh) | 2022-03-11 | 2022-03-11 | 一种区域电气综合能源系统小扰动稳定性评估方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115204549A true CN115204549A (zh) | 2022-10-18 |
Family
ID=83574918
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210235874.XA Pending CN115204549A (zh) | 2022-03-11 | 2022-03-11 | 一种区域电气综合能源系统小扰动稳定性评估方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115204549A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116150969A (zh) * | 2022-12-22 | 2023-05-23 | 国网甘肃省电力公司电力科学研究院 | 一种光储-虚拟同步发电机稳定性分析方法 |
-
2022
- 2022-03-11 CN CN202210235874.XA patent/CN115204549A/zh active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116150969A (zh) * | 2022-12-22 | 2023-05-23 | 国网甘肃省电力公司电力科学研究院 | 一种光储-虚拟同步发电机稳定性分析方法 |
| CN116150969B (zh) * | 2022-12-22 | 2024-01-30 | 国网甘肃省电力公司电力科学研究院 | 一种光储-虚拟同步发电机稳定性分析方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107291990B (zh) | 基于电-气互联综合能源系统暂态模型的能量流仿真方法 | |
| Sun et al. | A hybrid paradigm combining model-based and data-driven methods for fuel cell stack cooling control | |
| CN110647040B (zh) | 一种综合能源系统的安全控制方法及装置 | |
| CN115204549A (zh) | 一种区域电气综合能源系统小扰动稳定性评估方法 | |
| Mirrashid et al. | Robust control design for air breathing proton exchange membrane fuel cell system via variable gain second‐order sliding mode | |
| Tafreshi et al. | Linear dynamic parameter‐varying sliding manifold for air–fuel ratio control in lean‐burn engines | |
| Giménez et al. | Design and experimental validation of the temperature control of a PEMFC stack by applying multiobjective optimization | |
| CN113656941B (zh) | 一种面向综合能源系统的静态电压稳定性指标计算方法 | |
| Petkovic et al. | Modeling the transient behavior of gas turbines | |
| Wang et al. | Modeling and control for PEMFC hydrogen management subsystem based on neural network compensation and prescribed tracking accuracy | |
| Nayak et al. | A quantitative approach for convergence analysis of a singularly perturbed inverter-based microgrid | |
| Chen et al. | Observer based fuel delivery control for PEM fuel cells with a segmented anode model | |
| CN110277785B (zh) | 基于连续多能流的电气耦合系统负荷裕度计算方法及系统 | |
| CN116357463A (zh) | 氢混燃气轮机回火控制方法、装置、电子设备及存储介质 | |
| Palmieri et al. | An innovative sliding mode load controller for gas turbine power generators: Design and experimental validation via real-time simulation | |
| Bai et al. | The nonlinear single controller of DGEN380 aero engine design | |
| CN113885335B (zh) | 一种存在部分解耦扰动的网络化系统容错控制方法 | |
| Liu et al. | Compensators design for bumpless switching in aero‐engine multi‐loop control system | |
| Tang et al. | Off-equilibrium linearisation-based nonlinear control of turbojet enginese with sum-of-squares programming | |
| Lou et al. | Stabilisation for a class of switched nonlinear systems and its application to aero‐engines | |
| CN111796814A (zh) | 反应堆模块化微分方程模型的图形界面控制系统设计方法 | |
| Xie et al. | Research on Air Mass Flow and Pressure Control Method for the Multi-Stack Fuel Cell System Based on Model Predictive Control | |
| CN109002419B (zh) | 天然气管网的动态分析方法及装置 | |
| Messina et al. | Aergia number: A new non‐dimensional group for gas turbine power estimation | |
| Xi et al. | Robust Nonlinear Adaptive Backstepping Coordinated Control for Boiler-Turbine Units |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication |