CN113221358B - 基于可靠性参数的电-气耦合系统的备用出力优化方法 - Google Patents
基于可靠性参数的电-气耦合系统的备用出力优化方法 Download PDFInfo
- Publication number
- CN113221358B CN113221358B CN202110522290.6A CN202110522290A CN113221358B CN 113221358 B CN113221358 B CN 113221358B CN 202110522290 A CN202110522290 A CN 202110522290A CN 113221358 B CN113221358 B CN 113221358B
- Authority
- CN
- China
- Prior art keywords
- natural gas
- coal
- gas
- time
- electric
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 230000008878 coupling Effects 0.000 title claims abstract description 50
- 238000010168 coupling process Methods 0.000 title claims abstract description 50
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 50
- 238000005457 optimization Methods 0.000 title claims abstract description 19
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 365
- 239000007789 gas Substances 0.000 claims abstract description 188
- 239000003345 natural gas Substances 0.000 claims abstract description 175
- 238000004519 manufacturing process Methods 0.000 claims abstract description 57
- 230000008569 process Effects 0.000 claims abstract description 9
- 239000003245 coal Substances 0.000 claims description 68
- 238000003860 storage Methods 0.000 claims description 42
- 230000005540 biological transmission Effects 0.000 claims description 14
- 230000009194 climbing Effects 0.000 claims description 12
- 238000005520 cutting process Methods 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 9
- 238000010248 power generation Methods 0.000 claims description 8
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 6
- 239000003990 capacitor Substances 0.000 claims description 6
- 230000005611 electricity Effects 0.000 claims description 5
- 230000006835 compression Effects 0.000 claims description 4
- 238000007906 compression Methods 0.000 claims description 4
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 3
- 239000001569 carbon dioxide Substances 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000002360 preparation method Methods 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
- 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
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/06—Wind turbines or wind farms
-
- 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/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- 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)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Game Theory and Decision Science (AREA)
- Mathematical Physics (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Development Economics (AREA)
- Pure & Applied Mathematics (AREA)
- Quality & Reliability (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于可靠性参数的电‑气耦合系统的备用出力优化方法。本发明包括以下步骤:1)通过各类互感器采集电‑气耦合系统中各个设备的运行状态参数;2)构建电‑气耦合系统的可靠性模型;3)计算在故障传播过程中天然气系统与电力系统的可靠性参数:4)构建基于可靠性参数的电‑气耦合系统的备用出力优化模型;5)通过内点法求解电‑气耦合系统的备用出力优化模型,获得在故障下天然气系统的备用出力和电力系统的备用出力,调度服务器调整各产能设备的出力,使得各产能设备的出力满足电力负荷与天然气负荷需求。本发明可有效保障电‑气耦合系统中备用出力的充裕度,对于电‑气耦合系统的安全可靠运行具有重要作用。
Description
技术领域
本发明涉及了属于综合能源系统可靠性管理及优化领域的一种电-气耦合系统的备用出力优化方法,尤其是涉及了一种基于可靠性参数的电-气耦合系统的备用出力优化方法。
背景技术
近年来,为建设清洁、高效、低碳和可持续的能源体系,电、热和气等耦合的综合能源系统已成为国内外研究的热点。天然气机组因具备发电效率高、灵活性强和排放低等优点,被世界各国广泛采用,其发电容量比例也逐年增加。天然气机组的广泛建设使得电力系统与天然气系统间的耦合日渐紧密,电-气耦合系统日渐成型。
电力与天然气耦合交互在促进能源行业提质增效的同时,也为故障跨系统传播提供了路径,带来诸多可靠性问题。具体来说,天然气系统中出现气源供应中断、管道泄漏等随机故障均可能导致天然气机组供气中断,由此机组不得不降低出力,使得电力系统运行面临诸多挑战。美国德克萨斯州等大停电事故表明,无法有效应对多能耦合引起的风险问题将严重威胁电-气耦合系统的安全运行。
作为系统可靠性管理的重要资源,备用容量可为调度机构提供充足的调节裕度,以应对负荷变化和可再生能源波动等因素。考虑到电-气耦合系统的交互作用,电力备用与天然气备用间存在支撑特性。然而,现有的研究多集中于研究电力与天然气备用间如何协调优化,忽略了故障传播引起的可靠性问题。因此,如何在考虑故障传播的同时,协同配置电力与天然气备用,对于保障电-气耦合系统的安全可靠运行至关重要。
发明内容
针对上述背景技术中存在的问题,本发明提出了一种基于可靠性参数的电-气耦合系统的备用出力优化方法,针对电-气耦合系统,本发明首先通过电压、电流和气压等各类互感器获取各设备的运行状态及物理参数,以此作为备用出力优化模型的输入;紧接着,本发明根据各产能设备提供备用的可靠性模型,考虑故障传播过程依次构建天然气系统与电力系统的可靠性参数;最后,本发明构建考虑可靠性参数的备用出力优化模型,确定电-气耦合系统中气源、发电机组和储气装置的备用安排。本发明提出的基于可靠性参数的电-气耦合系统的备用出力优化方法,可有效保障系统备用容量的充裕度,对于电-气耦合系统的安全可靠运行具有重要作用。本发明可以直接应用于电力调度平台,在考虑故障传播的情况下确定各发电机组提供备用计划,所得的结果可直接用于电力系统调度过程,以保证系统的安全可靠运行。
本发明的技术方案如下:
本发明包括以下步骤:
1)通过各类互感器采集电-气耦合系统中各个产能设备、传输设备和用能设备的运行状态参数;
2)根据部分产能设备的运行状态参数,构建可靠性模型;
3)基于可靠性模型和剩余产能设备的运行状态参数,计算在故障传播过程中天然气系统与电力系统的可靠性参数:
4)基于电力系统和天然气系统的可靠性参数以及各个传输设备和用能设备的运行状态参数,构建基于可靠性参数的电-气耦合系统的备用出力优化模型;
5)通过内点法求解电-气耦合系统的备用出力优化模型,获得在故障下天然气系统的备用出力和电力系统的备用出力,调度服务器根据天然气系统的备用出力和电力系统的备用出力调整各产能设备的出力,使得各产能设备的出力满足电力负荷与天然气负荷需求,保证电-气耦合系统的安全可靠运行。
所述步骤1)中产能设备包括气源和发电机组,发电机组包括天然气机组、风电机组和煤电机组,传输设备包括天然气管道和电力线路,用能设备包括天然气负荷和电力负荷。
所述步骤2)具体为:
所述产能设备包括气源和发电机组,发电机组包括天然气机组、风电机组和煤电机组;可靠性模型包括气源可靠性模型和煤电机组可靠性模型;
2.1)根据气源的运行状态参数,构建气源可靠性模型,气源可靠性模型主要由气源的最大产气量模型和气源状态概率模型组成;
气源状态概率模型具体为各时刻t气源w处于气源状态v的概率prwv(t)的集合;
气源的最大产气量模型具体为:
2.2)根据煤电机组的运行状态参数,构建煤电机组可靠性模型,煤电机组可靠性模型主要由煤电机组的最大出力模型和煤电机组状态概率模型组成;
煤电机组状态概率模型具体为各时刻t煤电机组g处于煤电机组状态h的概率prgh(t)的集合;
煤电机组的最大出力模型具体为:
其中,Pgt表示时刻t的调度服务器中煤电机组g的正常出力,Rgt表示时刻t的调度服务器中煤电机组g的备用出力,表示时刻t处于煤电机组状态h的煤电机组g的装机容量;煤电机组状态h满足h=1,2,...,Kg;
所述步骤3)具体为:
3.1)基于气源可靠性模型,计算天然气系统的可靠性参数;
所述天然气系统的可靠性参数具体为天然气负荷切除量,天然气负荷切除量通过以下公式进行设置:
其中,EGNSt表示t时刻的天然气负荷切除量,DGt表示t时刻的总天然气负荷;GSgs,t表示t时刻储气装置gs中天然气的存储量;prws(t)表示t时刻气源w处于天然气系统状态s的概率,表示t时刻气源w处于天然气系统状态s下的最大产气量;Δt为时间间隔;为0-1变量,当处于天然气系统状态s时的负荷不足量表示天然气系统状态的总个数;
3.2)基于煤电机组可靠性模型以及剩余产能设备的运行状态参数,计算电力系统的可靠性参数:
所述电力系统的可靠性参数具体为电力负荷切除量,电力负荷切除量通过以下公式进行设置:
其中,DEt表示t时刻的总电力负荷;PGLe,t表示在电力系统状态ε下天然气机组的功率损失量;Pgg,t表示t时刻的天然气机组gg的出力;表示t时刻在电力系统状态ε下煤电机组g的最大出力;表示t时刻在电力系统状态ε下节点m处风电机组wd的出力,prgε(t)和分别表示t时刻煤电机组g和风电机组wd处于电力系统状态ε的概率;电力系统状态ε满足表示电力系统状态的总个数;为0-1变量,电力系统负荷切除量大于0时,为1,反之则为0;ηg2e表示天然气转化为电力的转化系数。
所述步骤4)具体为:
4.1)建立电-气耦合系统的备用出力目标函数:
min TC=OCE+OCG
其中,TC表示电-气耦合系统的总计量参数;OCE表示电力系统运行的总计量参数;OCG表示天然气系统运行的总计量参数;
所述电力系统运行的总计量参数OCE主要包括煤电机组的发电计量参数、备用计量参数和启停计量参数:
其中,ECg(·)表示煤电机组g的发电计量参数函数;Pgt表示时刻t的调度服务器中煤电机组g的正常出力,Rgt表示时刻t的调度服务器中煤电机组g的备用出力,RCgt表示煤电机组g的备用计量参数;UCg和DCeg分别表示煤电机组g的启动计量参数和停机计量参数;αgt和βgt为0-1变量,分别表示煤电机组g的启动状态和停机状态;
所述天然气系统运行的总计量参数OCG主要包括气源的产气计量参数、备用计量参数和储气装置的运行计量参数:
其中,GCw(·)表示气源w的产气计量参数函数;RCwt表示气源w的备用计量参数;Cgs表示t时刻储气装置gs的运行计量参数;GCgs,t和GDgs,t分别表示t时刻储气装置gs的进气量与出气量;Wwt表示时刻t的调度服务器中气源w的正常产气量;Rwt表示时刻t的调度服务器中气源w的备用产气量;
4.2)基于电力系统和天然气系统的可靠性参数以及各个传输设备和用能设备的运行状态参数,构建电-气耦合系统的约束条件;
I.构建天然气系统运行约束:
天然气管道上的天然气流量与天然气管道的两端气压有关,天然气管道上的天然气流量表示为:
其中,QD,t表示t时刻天然气管道D的天然气流量;πDi,t表示t时刻天然气管道D节点i处的气压;πDj,t表示t时刻天然气管道D节点j处的气压;CD为天然气管道D的气流传输参数;
天然气管道D的天然气流量满足以下约束:
天然气系统的各节点的天然气流入量与流出量相同,表示为:
天然气系统的各节点气压满足以下约束:
压缩机模型满足以下约束:
t时刻与t+1时刻储气装置gs的进气量、出气量和存量满足以下约束:
SGgs,t+1=SGgs,t+GCgs,t-GDgs,t
其中,SGgs,t表示t时刻储气装置gs的存量;
此外,储气装置还满足以下约束:
气源的正常产气量与备用产气量之和的约束表示为:
II.构建电力系统运行约束;
电力系统的节点功率平衡约束可表示为:
其中,Pmgt表示t时刻节点m的煤电机组g的出力,Pmgg,t表示t时刻节点m处天然气机组gg的出力;Pmwd,t表示t时刻节点m处风电机组wd的出力;DEmt表示t时刻节点m处的电力负荷;flt表示时刻t电力线路l的电力流量;为节点的电力流入量,为节点的电力流出量;
电力线路的电力流量约束表示为:
煤电机组的正常出力与备用出力满足以下约束:
其中,xmgt表示t时刻节点m处煤电机组g的运行状态,xmgt为0-1变量,0表示煤电机组g处于停机状态,1表示煤电机组g处于开机状态;和分别表示节点m处煤电机组g的最大与最小出力;Rmgt表示t时刻节点m的煤电机组g的备用出力;
煤电机组的备用出力还小于煤电机组的上爬坡速率:
相邻时刻煤电机组的正常出力的差值还小于下爬坡速率,表示为:
煤电机组从开机状态转化为停机状态时满足最小开机时间约束:
煤电机组从停机状态转化为开机状态时满足最小停机时间约束:
电力系统的节点m的相角约束表示为:
θmin≤θmt≤θmax
其中,θmax和θmin分别表示节点m相角的最大值和最小值;
III.构建可靠性约束;
电力系统和天然气系统的可靠性参数满足以下约束:
本发明的有益效果是:
本发明的针对现阶段电、气备用优化方法忽略了故障传播而使得各产能设备的备用安排无法满足可靠性要求,本发明在备用安排过程中将故障传播引起的可靠性问题考虑在内。本发明相较于传统方法更为有效可靠,所计算确定的备用安排更能保证电-气耦合系统的安全运行。
此外,本发明还可直接应用到现阶段电力系统调度系统中,通过研究天然气系统与电力系统间的故障传播过程,精确高效地确定系统中发电机组、气源和储气装置的备用计划,对于电-气耦合系统的可靠性提升具有重要意义。
附图说明
图1为本发明方法的逻辑框图。
图2为本发明方法效果与传统方法效果的对比图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
按照本发明内容完整方法实施的具体实施例如下:
如图1所示,本发明包括以下步骤:
1)通过电压、电流和气压等各类互感器采集电-气耦合系统中各个产能设备、传输设备和用能设备的运行状态参数;
步骤1)中产能设备包括气源和发电机组,发电机组包括天然气机组、风电机组和煤电机组,传输设备包括天然气管道和电力线路,用能设备包括天然气负荷和电力负荷。
2)根据部分产能设备的运行状态参数,部分产能设备具体为气源和煤电机组,构建可靠性模型;
步骤2)具体为:
产能设备包括气源和发电机组,发电机组包括天然气机组、风电机组和煤电机组;可靠性模型包括气源可靠性模型和煤电机组可靠性模型;
2.1)根据气源的运行状态参数,构建气源可靠性模型,气源可靠性模型主要由气源的最大产气量模型和气源状态概率模型组成;
气源状态概率模型具体为各时刻t气源w处于气源状态v的概率prwv(t)的集合;
气源的最大产气量模型具体为:
具体实施中,v=1时表示气源w处于正常运行状态,时刻t后气源w的最大产气量为时刻t的调度服务器中气源w的正常产气量Wwt与备用产气量Rwt之和。
2.2)根据煤电机组的运行状态参数,构建煤电机组可靠性模型,煤电机组可靠性模型主要由煤电机组的最大出力模型和煤电机组状态概率模型组成;
煤电机组状态概率模型具体为各时刻t煤电机组g处于煤电机组状态h的概率prgh(t)的集合;
煤电机组的最大出力模型具体为:
其中,Pgt表示时刻t的调度服务器中煤电机组g的正常出力,Rgt表示时刻t的调度服务器中煤电机组g的备用出力,表示时刻t处于煤电机组状态h的煤电机组g的装机容量;煤电机组状态h满足h=1,2,...,Kg;
具体实施中,h=1时表示煤电机组g处于正常运行状态,时刻t后煤电机组g的最大出力为时刻t的调度服务器中煤电机组g的正常出力Pgt与备用出力Rgt之和。
3)基于可靠性模型和剩余产能设备的运行状态参数,计算在故障传播过程中天然气系统与电力系统的可靠性参数:
步骤3)具体为:
3.1)基于气源可靠性模型,计算天然气系统的可靠性参数;
天然气系统的可靠性参数具体为天然气负荷切除量,天然气负荷切除量为天然气系统的总最大产气量与总天然气负荷的差值,天然气负荷切除量通过以下公式进行设置:
其中,EGNSt表示t时刻的天然气负荷切除量,DGt表示t时刻的总天然气负荷;GSgs,t表示t时刻储气装置gs中天然气的存储量;prws(t)表示t时刻气源w处于天然气系统状态s的概率,表示t时刻气源w处于天然气系统状态s下的最大产气量;Δt为时间间隔;为0-1变量,当处于天然气系统状态s时的负荷不足量表示天然气系统状态的总个数;
3.2)基于煤电机组可靠性模型以及剩余产能设备的运行状态参数,计算电力系统的可靠性参数:剩余产能设备具体为风电机组和天然气机组;
电力系统的可靠性参数具体为电力负荷切除量,电力负荷切除量为发电机组的出力总和与总电力负荷的差值,发电机组的出力总和为风电机组的出力、煤电机组的最大出力和天然气机组的出力减去天然气机组功率损失量的差值;电力负荷切除量通过以下公式进行设置:
其中,DEt表示t时刻的总电力负荷;PGLε,t表示在电力系统状态ε下天然气机组的功率损失量;Pgg,t表示t时刻的天然气机组gg的出力;表示t时刻在电力系统状态ε下煤电机组g的最大出力;表示t时刻在电力系统状态ε下节点m处风电机组wd的出力,prgε(t)和分别表示t时刻煤电机组g和风电机组wd处于电力系统状态e的概率;电力系统状态e满足表示电力系统状态的总个数;为0-1变量,电力系统负荷切除量大于0时,为1,反之则为0;ηg2e表示天然气转化为电力的转化系数。
4)基于电力系统和天然气系统的可靠性参数以及各个传输设备和用能设备的运行状态参数,构建基于可靠性参数的电-气耦合系统的备用出力优化模型;
步骤4)具体为:
4.1)建立电-气耦合系统的备用出力目标函数:
min TC=OCE+OCG
其中,TC表示电-气耦合系统的总计量参数;OCE表示电力系统运行的总计量参数;OCG表示天然气系统运行的总计量参数;
电力系统运行的总计量参数OCE主要包括煤电机组的发电计量参数、备用计量参数和启停计量参数:
其中,ECg(·)表示煤电机组g的发电计量参数函数;Pgt表示时刻t的调度服务器中煤电机组g的正常出力,Rgt表示时刻t的调度服务器中煤电机组g的备用出力,RCgt表示煤电机组g的备用计量参数;UCg和DCeg分别表示煤电机组g的启动计量参数和停机计量参数;agt和βgt为0-1变量,分别表示煤电机组g的启动状态和停机状态;
天然气系统运行的总计量参数OCG主要包括气源的产气计量参数、备用计量参数和储气装置的运行计量参数:
其中,GCw(·)表示气源w的产气计量参数函数;RCwt表示气源w的备用计量参数;Cgs表示t时刻储气装置gs的运行计量参数;GCgs,t和GDgs,t分别表示t时刻储气装置gs的进气量与出气量;Wwt表示时刻t的调度服务器中气源w的正常产气量;Rwt表示时刻t的调度服务器中气源w的备用产气量;
4.2)基于电力系统和天然气系统的可靠性参数以及各个传输设备和用能设备的运行状态参数,构建电-气耦合系统的约束条件;
I.构建天然气系统运行约束:
在天然气系统运行中,气源、天然气管道和储气装置需要满足相应的运行约束。
根据Weymouth方程,天然气管道上的天然气流量与天然气管道的两端气压有关,天然气管道上的天然气流量表示为:
其中,QD,t表示t时刻天然气管道D的天然气流量;πDi,t表示t时刻天然气管道D节点i处的气压;πDj,t表示t时刻天然气管道D节点j处的气压;CD为天然气管道D的气流传输参数;
本实施例中,气流传输参数取值为6。
此外,天然气管道D的天然气流量满足以下约束:
本实施例中,气流传输参数取值为10×105m3/h和0。
天然气系统的各节点的天然气流入量与流出量相同,压缩机设置在天然气管道中,表示为:
天然气系统的各节点气压满足以下约束:
压缩机模型满足以下约束:
本实施例中,压缩机压缩系数的上限和下限取值为2和1。
考虑到充放气过程,t时刻与t+1时刻储气装置gs的进气量、出气量和存量满足以下约束:
SGgs,t+1=SGgs,t+GCgs,t-GDgs,t
其中,SGgs,t表示t时刻储气装置gs的存量;
此外,储气装置还满足以下约束:
本实施例中,储气装置的存量设置为存量上限值的50%,进气量与出气量的最大值设置为存量上限值的50%。
气源的正常产气量与备用产气量之和的约束表示为:
II.构建电力系统运行约束;
电力系统运行约束包括电力线路和发电机组约束。
电力系统的节点功率平衡约束可表示为:
其中,Pmgt表示t时刻节点m的煤电机组g的出力,Pmgg,t表示t时刻节点m处天然气机组gg的出力;Pmwd,t表示t时刻节点m处风电机组wd的出力;DEmt表示t时刻节点m处的电力负荷;flt表示时刻t电力线路l的电力流量;为节点的电力流入量,为节点的电力流出量;
电力线路的电力流量约束表示为:
实际运行中,煤电机组的正常出力与备用出力满足以下约束:
其中,xmgt表示t时刻节点m处煤电机组g的运行状态,xmgt为0-1变量,0表示煤电机组g处于停机状态,1表示煤电机组g处于开机状态;和分别表示节点m处煤电机组g的最大与最小出力;Rmgt表示t时刻节点m的煤电机组g的备用出力;
此外,煤电机组的备用出力还小于煤电机组的上爬坡速率:
相邻时刻煤电机组的正常出力的差值还小于下爬坡速率,表示为:
煤电机组从开机状态转化为停机状态时满足最小开机时间约束:
煤电机组从停机状态转化为开机状态时满足最小停机时间约束:
电力系统的节点m的相角约束表示为:
θmin≤θmt≤θmax
其中,θmax和θmin分别表示节点m相角的最大值和最小值;
III.构建可靠性约束;
电-气耦合系统在运行时需保证可靠性参数在给定范围内。为此,电力系统和天然气系统的可靠性参数满足以下约束:
5)通过内点法求解电-气耦合系统的备用出力优化模型,获得在故障下天然气系统的备用出力和电力系统的备用出力,调度服务器根据天然气系统的备用出力和电力系统的备用出力调整各产能设备的出力,使得各产能设备的出力满足电力负荷与天然气负荷需求,保证电-气耦合系统的安全可靠运行。
Claims (4)
1.一种基于可靠性参数的电-气耦合系统的备用出力优化方法,其特征在于:包括以下步骤:
1)通过各类互感器采集电-气耦合系统中各个产能设备、传输设备和用能设备的运行状态参数;
2)根据部分产能设备的运行状态参数,构建可靠性模型;
3)基于可靠性模型和剩余产能设备的运行状态参数,计算在故障传播过程中天然气系统与电力系统的可靠性参数:
4)基于电力系统和天然气系统的可靠性参数以及各个传输设备和用能设备的运行状态参数,构建基于可靠性参数的电-气耦合系统的备用出力优化模型;
5)通过内点法求解电-气耦合系统的备用出力优化模型,获得在故障下天然气系统的备用出力和电力系统的备用出力,调度服务器根据天然气系统的备用出力和电力系统的备用出力调整各产能设备的出力,使得各产能设备的出力满足电力负荷与天然气负荷需求,保证电-气耦合系统的安全可靠运行;
所述步骤3)具体为:
3.1)基于气源可靠性模型,计算天然气系统的可靠性参数;
所述天然气系统的可靠性参数具体为天然气负荷切除量,天然气负荷切除量通过以下公式进行设置:
其中,EGNSt表示t时刻的天然气负荷切除量,DGt表示t时刻的总天然气负荷;GSgs,t表示t时刻储气装置gs中天然气的存储量;prws(t)表示t时刻气源w处于天然气系统状态s的概率,表示t时刻气源w处于天然气系统状态s下的最大产气量;Δt为时间间隔;为0-1变量,当处于天然气系统状态s时的负荷不足量大于0时,为1,反之当处于天然气系统状态s时的负荷不足量小于等于0时,则为0;天然气系统状态s满足表示天然气系统状态的总个数;
3.2)基于煤电机组可靠性模型以及剩余产能设备的运行状态参数,计算电力系统的可靠性参数:
所述电力系统的可靠性参数具体为电力负荷切除量,电力负荷切除量通过以下公式进行设置:
2.根据权利要求1所述的一种基于可靠性参数的电-气耦合系统的备用出力优化方法,其特征在于:所述步骤1)中产能设备包括气源和发电机组,发电机组包括天然气机组、风电机组和煤电机组,传输设备包括天然气管道和电力线路,用能设备包括天然气负荷和电力负荷。
3.根据权利要求1所述的一种基于可靠性参数的电-气耦合系统的备用出力优化方法,其特征在于:所述步骤2)具体为:
所述产能设备包括气源和发电机组,发电机组包括天然气机组、风电机组和煤电机组;可靠性模型包括气源可靠性模型和煤电机组可靠性模型;
2.1)根据气源的运行状态参数,构建气源可靠性模型,气源可靠性模型由气源的最大产气量模型和气源状态概率模型组成;
气源状态概率模型具体为各时刻t气源w处于气源状态v的概率prwv(t)的集合;
气源的最大产气量模型具体为:
2.2)根据煤电机组的运行状态参数,构建煤电机组可靠性模型,煤电机组可靠性模型由煤电机组的最大出力模型和煤电机组状态概率模型组成;
煤电机组状态概率模型具体为各时刻t煤电机组g处于煤电机组状态h的概率prgh(t)的集合;
煤电机组的最大出力模型具体为:
4.根据权利要求1所述的一种基于可靠性参数的电-气耦合系统的备用出力优化方法,其特征在于:所述步骤4)具体为:
4.1)建立电-气耦合系统的备用出力目标函数:
min TC=OCE+OCG
其中,TC表示电-气耦合系统的总计量参数;OCE表示电力系统运行的总计量参数;OCG表示天然气系统运行的总计量参数;
所述电力系统运行的总计量参数OCE主要包括煤电机组的发电计量参数、备用计量参数和启停计量参数:
其中,ECg(·)表示煤电机组g的发电计量参数函数;Pgt表示时刻t的调度服务器中煤电机组g的正常出力,Rgt表示时刻t的调度服务器中煤电机组g的备用出力,RCgt表示煤电机组g的备用计量参数;UCg和DCeg分别表示煤电机组g的启动计量参数和停机计量参数;αgt和βgt为0-1变量,分别表示煤电机组g的启动状态和停机状态;
所述天然气系统运行的总计量参数OCG主要包括气源的产气计量参数、备用计量参数和储气装置的运行计量参数:
其中,GCw(·)表示气源w的产气计量参数函数;RCwt表示气源w的备用计量参数;Cgs表示t时刻储气装置gs的运行计量参数;GCgs,t和GDgs,t分别表示t时刻储气装置gs的进气量与出气量;Rwt表示时刻t的调度服务器中气源w的备用产气量;
4.2)基于电力系统和天然气系统的可靠性参数以及各个传输设备和用能设备的运行状态参数,构建电-气耦合系统的约束条件;
I.构建天然气系统运行约束:
天然气管道上的天然气流量与天然气管道的两端气压有关,天然气管道上的天然气流量表示为:
其中,QD,t表示t时刻天然气管道D的天然气流量;πDi,t表示t时刻天然气管道D节点i处的气压;πDj,t表示t时刻天然气管道D节点j处的气压;CD为天然气管道D的气流传输参数;
天然气管道D的天然气流量满足以下约束:
天然气系统的各节点的天然气流入量与流出量相同,表示为:
其中,DGit表示t时刻节点i的天然气负荷;tDci,t表示t时刻天然气管道D节点i处压缩机c的天然气注入量;为节点i的天然气流入量,节点i的天然气流出量;Wiwt表示时刻t节点i的调度服务器中气源w的正常产气量;
天然气系统的各节点气压满足以下约束:
压缩机模型满足以下约束:
t时刻与t+1时刻储气装置gs的进气量、出气量和存量满足以下约束:
SGgs,t+1=SGgs,t+GCgs,t-GDgs,t
其中,SGgs,t表示t时刻储气装置gs的存量;
此外,储气装置还满足以下约束:
气源的正常产气量与备用产气量之和的约束表示为:
II.构建电力系统运行约束;
电力系统的节点功率平衡约束可表示为:
其中,Pmgt表示t时刻节点m的煤电机组g的出力,Pmgg,t表示t时刻节点m处天然气机组gg的出力;Pmwd,t表示t时刻节点m处风电机组wd的出力;DEmt表示t时刻节点m处的电力负荷;flt表示时刻t电力线路l的电力流量;为节点的电力流入量,为节点的电力流出量;
电力线路的电力流量约束表示为:
flt=|Bl(θmt-θnt)|≤Fl max
其中,θmt表示t时刻节点m的相角;Bl表示电力线路l的电纳;Fl max表示电力线路l的最大电力流量;
煤电机组的正常出力与备用出力满足以下约束:
其中,xmgt表示t时刻节点m处煤电机组g的运行状态,xmgt为0-1变量,0表示煤电机组g处于停机状态,1表示煤电机组g处于开机状态;和分别表示节点m处煤电机组g的最大与最小出力;Rmgt表示t时刻节点m的煤电机组g的备用出力;
煤电机组的备用出力还小于煤电机组的上爬坡速率:
相邻时刻煤电机组的正常出力的差值还小于下爬坡速率,表示为:
煤电机组从开机状态转化为停机状态时满足最小开机时间约束:
煤电机组从停机状态转化为开机状态时满足最小停机时间约束:
电力系统的节点m的相角约束表示为:
θmin≤θmt≤θmax
其中,θmax和θmin分别表示节点m相角的最大值和最小值;
III.构建可靠性约束;
电力系统和天然气系统的可靠性参数满足以下约束:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110522290.6A CN113221358B (zh) | 2021-05-13 | 2021-05-13 | 基于可靠性参数的电-气耦合系统的备用出力优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110522290.6A CN113221358B (zh) | 2021-05-13 | 2021-05-13 | 基于可靠性参数的电-气耦合系统的备用出力优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113221358A CN113221358A (zh) | 2021-08-06 |
CN113221358B true CN113221358B (zh) | 2022-09-27 |
Family
ID=77095367
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110522290.6A Active CN113221358B (zh) | 2021-05-13 | 2021-05-13 | 基于可靠性参数的电-气耦合系统的备用出力优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113221358B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116362478B (zh) * | 2023-02-15 | 2024-03-22 | 浙江大学 | 考虑综合能源枢纽灵活性的电-气耦合系统风险调度方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108921404A (zh) * | 2018-06-15 | 2018-11-30 | 浙江大学 | 一种基于通用生成函数的电-气耦合系统的可靠性评估方法 |
WO2019200891A1 (zh) * | 2018-04-16 | 2019-10-24 | 清华大学 | 计及电-气系统耦合的电力系统电压稳定裕度计算方法 |
WO2019200662A1 (zh) * | 2018-04-20 | 2019-10-24 | 东北大学 | 电-热-气综合能源系统的稳定评估与静态控制方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103066620B (zh) * | 2012-12-24 | 2014-10-22 | 中国电力科学研究院 | 间歇能源并网下的自动发电控制模型的设计方法 |
CN108767852B (zh) * | 2018-06-15 | 2020-05-08 | 浙江大学 | 一种电-气耦合系统节点可靠性评估方法 |
CN109615248B (zh) * | 2018-12-14 | 2023-02-10 | 浙江大学 | 一种考虑天然气系统影响的电力系统连锁故障估计方法 |
CN112260261B (zh) * | 2020-09-08 | 2022-05-27 | 浙江大学 | 一种电-气综合能源系统脆弱性评估与改善方法 |
CN112182907B (zh) * | 2020-10-19 | 2022-12-27 | 贵州电网有限责任公司 | 一种基于可靠性约束的电-气耦合系统储能装置规划方法 |
CN112288592B (zh) * | 2020-10-20 | 2022-10-25 | 东南大学 | 一种气热电耦合系统scuc优化调度方法、装置及存储介质 |
CN112736923B (zh) * | 2020-12-08 | 2023-06-20 | 华南理工大学 | 一种考虑电网故障影响的天然气网连锁故障评估控制方法 |
-
2021
- 2021-05-13 CN CN202110522290.6A patent/CN113221358B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019200891A1 (zh) * | 2018-04-16 | 2019-10-24 | 清华大学 | 计及电-气系统耦合的电力系统电压稳定裕度计算方法 |
WO2019200662A1 (zh) * | 2018-04-20 | 2019-10-24 | 东北大学 | 电-热-气综合能源系统的稳定评估与静态控制方法 |
CN108921404A (zh) * | 2018-06-15 | 2018-11-30 | 浙江大学 | 一种基于通用生成函数的电-气耦合系统的可靠性评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113221358A (zh) | 2021-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108110804B (zh) | 一种含抽水蓄能机组的电网优化调度方法 | |
CN109167383B (zh) | 基于精确线性化电力网络模型的电力系统调峰优化方法 | |
CN111769553A (zh) | 考虑风电不确定性的气电互联系统优化运行方法 | |
CN102493915B (zh) | 风电场的风机调用方法及风机调用系统 | |
CN210111621U (zh) | 一种风电场风力发电机组储能式后备电源系统 | |
CN106856333B (zh) | 一种风光火打捆外送系统调峰容量分配确定方法 | |
CN106779377A (zh) | 风电消纳过程评估方法 | |
CN106600022B (zh) | 基于多目标优化的风-光-燃气-海水抽水蓄能孤立电力系统容量优化配置方法 | |
CN113937819B (zh) | 一种多能源短期优化调度方法 | |
CN110492524A (zh) | 无通信海上风电场柔直并网的有功-频率控制方法及系统 | |
CN105978008A (zh) | 一种具有风场黑启动功能的液流电池储能系统及其工作方法 | |
CN109992818A (zh) | 大规模风电参与一次调频的机组组合模型及求解方法 | |
CN113221358B (zh) | 基于可靠性参数的电-气耦合系统的备用出力优化方法 | |
CN114583766B (zh) | 面向风电消纳的海上油气平台能源系统低碳运行调度方法 | |
CN114142536B (zh) | 一种考虑容量备用的多类型机组协调方法 | |
CN109615125B (zh) | 一种考虑特高压调峰的多区域随机生产模拟方法及应用 | |
CN102496962A (zh) | 调峰调频约束下电力系统可消纳风电的辨识与控制方法 | |
CN112491091B (zh) | 一种大规模区域电网日前安全约束机组组合快速求解方法 | |
CN111342453B (zh) | 计及多类型备用资源的电气综合能源系统备用决策方法 | |
CN110380447B (zh) | 一种风机失效下的电-气互联能源系统降风险调度方法 | |
CN104484772A (zh) | 电量计划的可行性校验方法和系统 | |
CN116562563A (zh) | 一种基于最小惯量需求评估的电力系统优化运行方法 | |
CN106026197B (zh) | 风电与火电耦合外送优化控制方法 | |
CN108596521B (zh) | 一种区外来电对受端电网辅助服务影响的定量分析方法 | |
CN112434841A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |