CN114094568A - 一种发电-储能系统的波动压力发电控制方法及系统 - Google Patents

一种发电-储能系统的波动压力发电控制方法及系统 Download PDF

Info

Publication number
CN114094568A
CN114094568A CN202111264203.8A CN202111264203A CN114094568A CN 114094568 A CN114094568 A CN 114094568A CN 202111264203 A CN202111264203 A CN 202111264203A CN 114094568 A CN114094568 A CN 114094568A
Authority
CN
China
Prior art keywords
gas
air
pressure
liquid mixing
liquid
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.)
Granted
Application number
CN202111264203.8A
Other languages
English (en)
Other versions
CN114094568B (zh
Inventor
寇攀高
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
State Grid Hunan Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
State Grid Hunan Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd, State Grid Hunan Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202111264203.8A priority Critical patent/CN114094568B/zh
Publication of CN114094568A publication Critical patent/CN114094568A/zh
Application granted granted Critical
Publication of CN114094568B publication Critical patent/CN114094568B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/002—Flicker reduction, e.g. compensation of flicker introduced by non-linear load
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D16/00—Control of fluid pressure
    • G05D16/02—Modifications to reduce the effects of instability, e.g. due to vibrations, friction, abnormal temperature, overloading or unbalance
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D16/00—Control of fluid pressure
    • G05D16/20—Control of fluid pressure characterised by the use of electric means
    • G05D16/2093—Control of fluid pressure characterised by the use of electric means with combination of electric and non-electric auxiliary power
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J15/00—Systems for storing electric energy
    • H02J15/006—Systems for storing electric energy in the form of pneumatic energy, e.g. compressed air energy storage [CAES]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28—Arrangements for balancing of the load in a network by storage of energy
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/16—Mechanical energy storage, e.g. flywheels or pressurised fluids
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E70/00—Other energy conversion or management systems reducing GHG emissions
    • Y02E70/30—Systems combining energy storage with energy generation of non-fossil origin

Abstract

本发明公开了一种发电‑储能系统的波动压力发电控制方法及系统，本发明包括分别获取发电‑储能系统的测量气压、实测功率值以及高压气系统、气液混合系统之间气压调节阀的实测流量值，并以此生成有功功率误差信号、气压误差信号以及生成流量误差信号；将气压误差信号、有功功率误差信号以及流量误差信号求和生成多变量控制信号并控制高压气系统、气液混合系统之间气压调节阀的流量大小以维持气液混合系统中气体气压恒定。本发明通过多变量控制信号控制高压气系统、气液混合系统之间气压调节阀的流量大小以维持气液混合系统中气体气压恒定，能够提升发电机组功率的稳定性、机组的安全稳定运行，防止电力系统的低频振荡。

Description

一种发电-储能系统的波动压力发电控制方法及系统
技术领域
本发明涉及发电-储能系统的发电控制技术，具体涉及一种发电-储能系统的波动压力发电控制方法及系统。
背景技术
随着新能源发电的超常规发展与电网建设相对滞后的矛盾日益明显，大规模具有随机性、问歇性、反调节性及出力波动大等特点的风电/光伏能源接入电网对系统的电压稳定、暂态稳定和频率稳定都有较大的影响，风电/光伏能源并网难、并网后消纳难等问题严重制约着能源结构的变革。常规水电厂、抽液蓄能电厂在大规模新能源存储、能量转化方面作用有限，不能吸收丰沛的风电、太阳能等大规模可再生新能源电力，且对地势、地质有一定的要求。
申请号为202020451871.6的中国专利文献公开了一种发电-储能系统，该发电-储能系统包括高压气系统和气液混合系统，高压气系统包括储气容器和用于给储气容器供气的气体压缩装置，气液混合系统包括分别连接有补液循环系统和液力发电机组的气液混合容器，储气容器、气液混合容器之间通过调节阀相连，该发电-储能系统可通过储气容器、气液混合容器来实现储能、发电，也可以选择多种运行模式。该发电-储能系统作为一种新型发电系统，具有清洁无污染、布置灵活、运行寿命长等优点，对于大规模新能源消纳及提升电力系统的灵活性具有重要价值。然而，发电过程中气液混合系统内液体体积减小、压力趋向于下降，高压气系统通过压力调节阀对气液混合系统内补气，试图维持气压稳定，然而气压调节阀输出端气压难以维持恒定值，实践表明，气压调节阀输出气压波动幅度与输出气压平均值相比不小于5％。气液混合系统内气压的波动不仅严重影响了发电机组功率的稳定、机组的安全稳定运行，严重情况下会导致电力系统的低频振荡，因此提出维持气液混合系统气压稳定的控制方法对理论研究及工程实践具有重要的科学意义及实践价值。
发明内容
本发明要解决的技术问题：针对现有技术的上述问题，提供一种发电-储能系统的波动压力发电控制方法及系统，本发明通过多变量控制信号控制高压气系统、气液混合系统之间气压调节阀的流量大小以维持气液混合系统中气体气压恒定，能够提升发电机组功率的稳定性、机组的安全稳定运行，防止电力系统的低频振荡。
为了解决上述技术问题，本发明采用的技术方案为：
一种发电-储能系统的波动压力发电控制方法，包括：
1)分别获取发电-储能系统的测量气压Pmeas、实测功率值Nmeas以及高压气系统、气液混合系统之间气压调节阀的实测流量值
Figure BDA0003326454760000021
2)根据测量气压Pmeas、设定气压值Pref之间的差值生成气压误差信号，根据实测功率值Nmeas、目标功率值Nref之间的差值生成有功功率误差信号，根据目标流量值
Figure BDA0003326454760000022
实测流量值
Figure BDA0003326454760000023
之间的差值生成流量误差信号；
3)将气压误差信号、有功功率误差信号以及流量误差信号求和生成多变量控制信号；
4)通过多变量控制信号控制高压气系统、气液混合系统之间气压调节阀的流量大小以维持气液混合系统中气体气压恒定。
可选地，所述设定气压值Pref、目标功率值Nref均为常数值，所述目标流量值
Figure BDA0003326454760000024
为根据气液混合系统的动态特性计算得到。
可选地，步骤2)中还包括根据气液混合系统的动态特性计算目标流量值
Figure BDA0003326454760000025
的步骤：
S1)将高压气系统、气液混合系统之间的气体视为两个气体控制体，气液混合系统中的液体视为液体控制体，气液混合系统中的器壁视为器壁控制体，分别建立发电阶段气体控制体和液体控制体的温度、体积、压力动态变化特性方程、液体控制体的方程、与气体接触器壁控制体的方程以及与液体接触器壁控制体的方程；
S2)针对发电阶段气体控制体和液体控制体的温度、体积、压力动态变化特性方程、液体控制体的方程、与气体接触器壁控制体的方程以及与液体接触器壁控制体的方程采用“乒乓”策略求解得到最优的目标流量值
Figure BDA0003326454760000026
可选地，步骤S1)中建立的发电阶段气体控制体和液体控制体的温度、体积、压力动态变化特性方程为：
Figure BDA0003326454760000031
式(1)中，TB,gas,mB,gas分别为气液混合系统内空气控制体B的温度、质量，TA,gas,mA,gas分别为高压气系统内空气控制体A的温度、质量，t为时间，k为空气比热容的比值，
Figure BDA0003326454760000032
为气液混合系统内空气控制体B的质量流入速率，hl表示气液混合系统的液体控制体与器壁换热系数，Al表示气液混合系统的液体控制体与器壁接触的面积，Tl表示气液混合系统的液体控制体的温度，
Figure BDA0003326454760000033
表示气液混合系统内空气控制体B的空气定容比热容，U为与环境接触器壁的换热系数，AG为气液混合系统内空气控制体B与气体接触的器壁面积，TB,gas,T为气液混合系统内空气控制体B接触的器壁控制体温度，pB,gas为气液混合系统内空气控制体B的压力，
Figure BDA0003326454760000034
为高压气系统内空气控制体A的质量流出速率，ρw为汽液混合容器内液体控制体的密度，TA,gas,T为气液混合系统内空气控制体A接触的器壁控制体温度，
Figure BDA0003326454760000035
表示气液混合系统内空气控制体A的空气定容比热容，
Figure BDA0003326454760000036
为汽液混合容器内液体控制器的质量流出流率，
Figure BDA0003326454760000037
为高压气系统内空气控制体A、气液混合系统内空气控制体B之间的交换速率，且满足
Figure BDA0003326454760000038
Rg为空气气体常数；
液体控制体的方程为：
Figure BDA0003326454760000039
式(2)中，ρl为液体密度，Al为液体控制体水平面的面积，cl为液体控制体的比热容，Zl为液体控制体的液位，pgas为气液混合容器内空气压力,pgas＝pB,gas，
Figure BDA0003326454760000041
为液体控制体出口液体温度，
Figure BDA0003326454760000042
为液体控制体出口的质量流率，hgas,l为液体控制体与气体控制体之间的换热系数，hl,T为液体与汽液混合容器器壁之间的换热系数，Tgas为气液混合容器的空气温度，Al,T为液体控制体与汽液混合容器器壁接触的面积，Tl,T为与液体控制体接触的容器器壁温度，
Figure BDA0003326454760000043
表示抽液过程中的质量流量；
与气体接触器壁控制体的方程为：
Figure BDA0003326454760000044
式(3)中，mgas,T为与空气接触的容器器壁控制体质量，cT为气液混合容器器壁的比热容，Tgas,T为与汽液混合容器器壁接触的气体控制体的温度，Agas,T为气体控制体与汽液混合容器器壁接触的面积，Tgas为气液混合容器的空气温度，ho为气液混合容器器壁与外部环境之间换热系数，Ao,gas为与气液混合容器内部空气接触的器壁暴露在外部环境的面积，Tamb为外部环境温度；
发电过程中与液体接触器壁控制体的方程为：
Figure BDA0003326454760000045
式(4)中，ml,T为与液体接触的容器器壁控制体质量，Ao,l为与液体接触的容器器壁和外部环境的接触面积。
可选地，步骤S2)包括：
S2.1)初始化迭代变量n为1，初始化目标流量
Figure BDA0003326454760000046
以及状态变量；根据nΔt时刻的目标流量
Figure BDA0003326454760000047
以及状态变量，结合式(1)中气液混合系统内空气控制体B的压力pB,gas的计算函数表达式计算得到nΔt时刻的气液混合系统内空气控制体B的压力
Figure BDA0003326454760000048
S2.2)将nΔt时刻的状态变量代入式(1)～(4)，并采用数值差分的方法计算出(n+1)Δt时刻的状态变量；
S2.3)根据(n+1)Δt时刻的状态变量结合式(1)中气液混合系统内空气控制体B的压力pB,gas的计算函数表达式计算得到(n+1)Δt时刻的气液混合系统内空气控制体B的压力
Figure BDA0003326454760000049
S2.4)判断nΔt时刻的气液混合系统内空气控制体B的压力
Figure BDA0003326454760000051
(n+1)Δt时刻的气液混合系统内空气控制体B的压力
Figure BDA0003326454760000052
两者之间的差值的绝对值小于预设阈值是否成立，若成立，则将当前的目标流量
Figure BDA0003326454760000053
作为求解得到的最优的目标流量值
Figure BDA0003326454760000054
输出，结束并退出；否则，根据下式更新当前的目标流量
Figure BDA0003326454760000055
Figure BDA0003326454760000056
式(5)中，
Figure BDA0003326454760000057
为更新后的目标流量，Δ为质量流量的搜索步长，a为质量流量
Figure BDA0003326454760000058
搜索范围的宽度，σ为质量流量
Figure BDA0003326454760000059
搜索范围宽度的缩放因子，m为第m次动态调整质量流量
Figure BDA00033264547600000510
的搜索范围。
可选地，步骤2)中根据测量气压Pmeas、设定气压值Pref之间的差值生成气压误差信号是指将测量气压Pmeas、设定气压值Pref之间的差值输入气压闭环控制器得到气压误差信号。
可选地，步骤2)中根据实测功率值Nmeas、目标功率值Nref之间的差值生成有功功率误差信号是指将实测功率值Nmeas、目标功率值Nref之间的差值输入有功功率闭环控制器得到有功功率误差信号。
可选地，步骤2)中根据目标流量值
Figure BDA00033264547600000511
实测流量值
Figure BDA00033264547600000512
之间的差值生成流量误差信号是指将目标流量值
Figure BDA00033264547600000513
实测流量值
Figure BDA00033264547600000514
之间的差值输入流量闭环控制器得到流量误差信号。
此外，本发明还提供一种发电-储能系统的波动压力发电控制系统，包括相互连接的微处理器和存储器，该微处理器被编程或配置以执行所述发电-储能系统的波动压力发电控制方法的步骤。
此外，本发明还提供一种计算机可读存储介质，该计算机可读存储介质中存储有被编程或配置以执行所述发电-储能系统的波动压力发电控制方法的计算机程序。
和现有技术相比，本发明具有下述优点：本发明包括分别获取发电-储能系统的测量气压、实测功率值以及高压气系统、气液混合系统之间气压调节阀的实测流量值，并以此生成有功功率误差信号、气压误差信号以及生成流量误差信号；将气压误差信号、有功功率误差信号以及流量误差信号求和生成多变量控制信号并控制高压气系统、气液混合系统之间气压调节阀的流量大小以维持气液混合系统中气体气压恒定。本发明通过多变量控制信号控制高压气系统、气液混合系统之间气压调节阀的流量大小以维持气液混合系统中气体气压恒定，能够提升发电机组功率的稳定性、机组的安全稳定运行，防止电力系统的低频振荡。
附图说明
图1为本发明实施例方法的控制原理图。
图2为本发明实施例中求解目标流量值的流程示意图。
具体实施方式
为了解决发电-储能系统中与高压气系统、气液混合系统相联结的压力调节动态调节导致气液混合系统内气压波动问题，如图1所示，本实施例发电-储能系统的波动压力发电控制方法包括：
1)分别获取发电-储能系统的测量气压Pmeas、实测功率值Nmeas以及高压气系统、气液混合系统之间气压调节阀的实测流量值
Figure BDA0003326454760000061
2)根据测量气压Pmeas、设定气压值Pref之间的差值生成有功功率误差信号，根据实测功率值Nmeas、目标功率值Nref之间的差值生成气压误差信号，根据目标流量值
Figure BDA0003326454760000062
实测流量值
Figure BDA0003326454760000063
之间的差值生成流量误差信号；
3)将气压误差信号、有功功率误差信号以及流量误差信号求和生成多变量控制信号；
4)通过多变量控制信号控制高压气系统、气液混合系统之间气压调节阀的流量大小以维持气液混合系统中气体气压恒定。
本实施例中，设定气压值Pref、目标功率值Nref均为常数值，所述目标流量值
Figure BDA0003326454760000064
为根据气液混合系统的动态特性计算得到。
本实施例中，步骤2)中还包括根据气液混合系统的动态特性计算目标流量值
Figure BDA0003326454760000065
的步骤：
S1)将高压气系统、气液混合系统之间的气体视为两个气体控制体，气液混合系统中的液体视为液体控制体，气液混合系统中的器壁视为器壁控制体，分别建立发电阶段气体控制体和液体控制体的温度、体积、压力动态变化特性方程、液体控制体的方程、与气体接触器壁控制体的方程以及与液体接触器壁控制体的方程；
S2)针对发电阶段气体控制体和液体控制体的温度、体积、压力动态变化特性方程、液体控制体的方程、与气体接触器壁控制体的方程以及与液体接触器壁控制体的方程采用“乒乓”策略求解得到最优的目标流量值
Figure BDA0003326454760000066
本实施例中，步骤S1)中建立的发电阶段气体控制体和液体控制体的温度、体积、压力动态变化特性方程为：
Figure BDA0003326454760000071
式(1)中，TB,gas,mB,gas分别为气液混合系统内空气控制体B的温度、质量，TA,gas,mA,gas分别为高压气系统内空气控制体A的温度、质量，t为时间，k为空气比热容的比值，
Figure BDA0003326454760000072
为气液混合系统内空气控制体B的质量流入速率，hl表示气液混合系统的液体控制体与器壁换热系数，Al表示气液混合系统的液体控制体与器壁接触的面积，Tl表示气液混合系统的液体控制体的温度，
Figure BDA0003326454760000073
表示气液混合系统内空气控制体B的空气定容比热容，U为与环境接触器壁的换热系数，AG为气液混合系统内空气控制体B与气体接触的器壁面积，TB,gas,T为气液混合系统内空气控制体B接触的器壁控制体温度，pB,gas为气液混合系统内空气控制体B的压力，
Figure BDA0003326454760000074
为高压气系统内空气控制体A的质量流出速率，ρw为汽液混合容器内液体控制体的密度，TA,gas,T为气液混合系统内空气控制体A接触的器壁控制体温度，
Figure BDA0003326454760000075
表示气液混合系统内空气控制体A的空气定容比热容，
Figure BDA0003326454760000076
为汽液混合容器内液体控制器的质量流出流率，
Figure BDA0003326454760000077
为高压气系统内空气控制体A、气液混合系统内空气控制体B之间的交换速率，且满足
Figure BDA0003326454760000078
Rg为空气气体常数；
液体控制体的方程为：
Figure BDA0003326454760000079
式(2)中，ρl为液体密度，Al为液体控制体水平面的面积，cl为液体控制体的比热容，Zl为液体控制体的液位，pgas为气液混合容器内空气压力,pgas＝pB,gas，
Figure BDA0003326454760000081
为液体控制体出口液体温度，
Figure BDA0003326454760000082
为液体控制体出口的质量流率，hgas,l为液体控制体与气体控制体之间的换热系数，hl,T为液体与汽液混合容器器壁之间的换热系数，Tgas为气液混合容器的空气温度，Al,T为液体控制体与汽液混合容器器壁接触的面积，Tl,T为与液体控制体接触的容器器壁温度，
Figure BDA0003326454760000083
表示抽液过程中的质量流量；
与气体接触器壁控制体的方程为：
Figure BDA0003326454760000084
式(3)中，mgas,T为与空气接触的容器器壁控制体质量，cT为气液混合容器器壁的比热容，Tgas,T为与汽液混合容器器壁接触的气体控制体的温度，Agas,T为气体控制体与汽液混合容器器壁接触的面积，Tgas为气液混合容器的空气温度，ho为气液混合容器器壁与外部环境之间换热系数，Ao,gas为与气液混合容器内部空气接触的器壁暴露在外部环境的面积，Tamb为外部环境温度；
发电过程中与液体接触器壁控制体的方程为：
Figure BDA0003326454760000085
式(4)中，ml,T为与液体接触的容器器壁控制体质量，Ao,l为与液体接触的容器器壁和外部环境的接触面积。
本实施例控制方法的控制目标为气液混合系统中气体气压恒定；本实施例控制方法为多变量控制方法，所述控制方法考虑了发电过程中高压气系统、气液混合系统中气体、液体、容器器壁温度的动态变化特性，将气液混合系统内气压信号、有功功率信号、开度信号、气压调节阀流量信号引入到控制器中形成了流量闭环控制、气压闭环控制、有功功率闭环控制器，实现了发电过程功率的稳定控制。控制方法中的流量闭环控制设定流量通过乒乓求解方法求出。压力误差信号、功率误差信号、流量误差信号分别经过各自控制器形成控制信号输出到气压调节阀，实现对气压调节阀过流量的控制。采用“乒乓”策略求解得到最优的目标流量值
Figure BDA0003326454760000086
的基本原理为：首先假设出一个
Figure BDA0003326454760000087
然后根据(1)-(4)求解出mA,gas和mB,gas和气压值，判断气压值是否满足误差需要，如果气压误差满足要求则假设质量流量值满足要求；否则修正
Figure BDA0003326454760000088
数值。修正时采用“乒乓”策略：搜索到了第一个质量流量
Figure BDA0003326454760000091
则以该质量流量为中心，在有限区域内遍历搜索下一个质量流量，当在该区域内寻找不到解答时，动态扩展该求解区域，直到搜索到解答。如图2所示，本实施例步骤S2)包括：
S2.1)初始化迭代变量n为1，初始化目标流量
Figure BDA0003326454760000092
以及状态变量；根据nΔt时刻的目标流量
Figure BDA0003326454760000093
以及状态变量，结合式(1)中气液混合系统内空气控制体B的压力pB,gas的计算函数表达式计算得到nΔt时刻的气液混合系统内空气控制体B的压力
Figure BDA0003326454760000094
S2.2)将nΔt时刻的状态变量代入式(1)～(4)，并采用数值差分的方法计算出(n+1)Δt时刻的状态变量；
S2.3)根据(n+1)Δt时刻的状态变量结合式(1)中气液混合系统内空气控制体B的压力pB,gas的计算函数表达式计算得到(n+1)Δt时刻的气液混合系统内空气控制体B的压力
Figure BDA0003326454760000095
S2.4)判断nΔt时刻的气液混合系统内空气控制体B的压力
Figure BDA0003326454760000096
(n+1)Δt时刻的气液混合系统内空气控制体B的压力
Figure BDA0003326454760000097
两者之间的差值的绝对值小于预设阈值是否成立，若成立，则将当前的目标流量
Figure BDA0003326454760000098
作为求解得到的最优的目标流量值
Figure BDA0003326454760000099
输出，结束并退出；否则，根据下式更新当前的目标流量
Figure BDA00033264547600000910
Figure BDA00033264547600000911
式(5)中，
Figure BDA00033264547600000912
为更新后的目标流量，Δ为质量流量的搜索步长，a为质量流量
Figure BDA00033264547600000913
搜索范围的宽度，σ为质量流量
Figure BDA00033264547600000914
搜索范围宽度的缩放因子，m为第m次动态调整质量流量
Figure BDA00033264547600000915
的搜索范围。
本实施例中，步骤2)中根据测量气压Pmeas、设定气压值Pref之间的差值生成气压误差信号是指将测量气压Pmeas、设定气压值Pref之间的差值输入气压闭环控制器得到气压误差信号。作为一种可选的实施方式，本实施例中有功功率闭环控制器采用PID控制器，此外也可以根据需要采用模糊控制器等其他控制器。
本实施例中，步骤2)中根据实测功率值Nmeas、目标功率值Nref之间的差值生成有功功率误差信号是指将实测功率值Nmeas、目标功率值Nref之间的差值输入有功功率闭环控制器得到有功功率误差信号。作为一种可选的实施方式，本实施例中压闭闭环控制器采用PID控制器，此外也可以根据需要采用模糊控制器等其他控制器。
本实施例中，步骤2)中根据目标流量值
Figure BDA00033264547600000916
实测流量值
Figure BDA00033264547600000917
之间的差值生成流量误差信号是指将目标流量值
Figure BDA0003326454760000101
实测流量值
Figure BDA0003326454760000102
之间的差值输入流量闭环控制器得到流量误差信号。作为一种可选的实施方式，本实施例中流量闭环控制器采用PID控制器，此外也可以根据需要采用模糊控制器等其他控制器。
综上所述，本实施例方法针对高压气系统和气液混合系统共同发电阶段气液混合系统内气压波动导致发电机组功率难以稳定的问题，考虑了气体、液体、容器器壁温度的动态变化，提出了高压气系统和气液混合系统动态模型的乒乓求解方法，将气液混合系统内气压控制问题转化为气压调节阀的流量控制，并将气液混合系统内气压信号、有功功率信号、开度信号、气压调节阀流量信号引入到控制器中，通过闭环控制和控制策略实现了发电过程功率的稳定控制，为该发电-储能系统实现压力波动下发电机组功率恒定输出提供了一种控制方法，为提升电力系统灵活性提供了重要的理论支撑与实践依据。
此外，本实施例还提供一种发电-储能系统的波动压力发电控制系统，包括相互连接的微处理器和存储器，该微处理器被编程或配置以执行前述发电-储能系统的波动压力发电控制方法的步骤。
此外，本实施例还提供一种计算机可读存储介质，该计算机可读存储介质中存储有被编程或配置以执行前述发电-储能系统的波动压力发电控制方法的计算机程序。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可读存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (10)

1.一种发电-储能系统的波动压力发电控制方法，其特征在于，包括：
1)分别获取发电-储能系统的测量气压Pmeas、实测功率值Nmeas以及高压气系统、气液混合系统之间气压调节阀的实测流量值
Figure FDA0003326454750000011
2)根据测量气压Pmeas、设定气压值Pref之间的差值生成气压误差信号，根据实测功率值Nmeas、目标功率值Nref之间的差值生成有功功率误差信号，根据目标流量值
Figure FDA0003326454750000012
实测流量值
Figure FDA0003326454750000013
之间的差值生成流量误差信号；
3)将气压误差信号、有功功率误差信号以及流量误差信号求和生成多变量控制信号；
4)通过多变量控制信号控制高压气系统、气液混合系统之间气压调节阀的流量大小以维持气液混合系统中气体气压恒定。
2.根据权利要求1所述的发电-储能系统的波动压力发电控制方法，其特征在于，所述设定气压值Pref、目标功率值Nref均为常数值，所述目标流量值
Figure FDA0003326454750000014
为根据气液混合系统的动态特性计算得到。
3.根据权利要求2所述的发电-储能系统的波动压力发电控制方法，其特征在于，步骤2)中还包括根据气液混合系统的动态特性计算目标流量值
Figure FDA0003326454750000015
的步骤：
S1)将高压气系统、气液混合系统之间的气体视为两个气体控制体，气液混合系统中的液体视为液体控制体，气液混合系统中的器壁视为器壁控制体，分别建立发电阶段气体控制体和液体控制体的温度、体积、压力动态变化特性方程、液体控制体的方程、与气体接触器壁控制体的方程以及与液体接触器壁控制体的方程；
S2)针对发电阶段气体控制体和液体控制体的温度、体积、压力动态变化特性方程、液体控制体的方程、与气体接触器壁控制体的方程以及与液体接触器壁控制体的方程采用“乒乓”策略求解得到最优的目标流量值
Figure FDA0003326454750000016
4.根据权利要求3所述的发电-储能系统的波动压力发电控制方法，其特征在于，步骤S1)中建立的发电阶段气体控制体和液体控制体的温度、体积、压力动态变化特性方程为：
Figure FDA0003326454750000021
式(1)中，TB,gas,mB,gas分别为气液混合系统内空气控制体B的温度、质量，TA,gas,mA,gas分别为高压气系统内空气控制体A的温度、质量，t为时间，k为空气比热容的比值，
Figure FDA0003326454750000022
为气液混合系统内空气控制体B的质量流入速率，hl表示气液混合系统的液体控制体与器壁换热系数，Al表示气液混合系统的液体控制体与器壁接触的面积，Tl表示气液混合系统的液体控制体的温度，
Figure FDA0003326454750000023
表示气液混合系统内空气控制体B的空气定容比热容，U为与环境接触器壁的换热系数，AG为气液混合系统内空气控制体B与气体接触的器壁面积，TB,gas,T为气液混合系统内空气控制体B接触的器壁控制体温度，pB,gas为气液混合系统内空气控制体B的压力，
Figure FDA0003326454750000024
为高压气系统内空气控制体A的质量流出速率，ρw为汽液混合容器内液体控制体的密度，TA,gas,T为气液混合系统内空气控制体A接触的器壁控制体温度，
Figure FDA0003326454750000025
表示气液混合系统内空气控制体A的空气定容比热容，
Figure FDA0003326454750000026
为汽液混合容器内液体控制器的质量流出流率，
Figure FDA0003326454750000027
为高压气系统内空气控制体A、气液混合系统内空气控制体B之间的交换速率，且满足
Figure FDA0003326454750000028
Rg为空气气体常数；
液体控制体的方程为：
Figure FDA0003326454750000029
式(2)中，ρl为液体密度，Al为液体控制体水平面的面积，cl为液体控制体的比热容，Zl为液体控制体的液位，pgas为气液混合容器内空气压力,pgas＝pB,gas，
Figure FDA0003326454750000031
为液体控制体出口液体温度，
Figure FDA0003326454750000032
为液体控制体出口的质量流率，hgas,l为液体控制体与气体控制体之间的换热系数，hl,T为液体与汽液混合容器器壁之间的换热系数，Tgas为气液混合容器的空气温度，Al,T为液体控制体与汽液混合容器器壁接触的面积，Tl,T为与液体控制体接触的容器器壁温度；
与气体接触器壁控制体的方程为：
Figure FDA0003326454750000033
式(3)中，mgas,T为与空气接触的容器器壁控制体质量，cT为气液混合容器器壁的比热容，Tgas,T为与汽液混合容器器壁接触的气体控制体的温度，Agas,T为气体控制体与汽液混合容器器壁接触的面积，Tgas为气液混合容器内的空气温度，ho为气液混合容器器壁与外部环境之间换热系数，Ao,gas为与气液混合容器内部空气接触的器壁暴露在外部环境的面积，Tamb为外部环境温度；
发电过程中与液体接触器壁控制体的方程为：
Figure FDA0003326454750000034
式(4)中，ml,T为与液体接触的容器器壁控制体质量，Ao,l为与液体接触的容器器壁和外部环境的接触面积。
5.根据权利要求4所述的发电-储能系统的波动压力发电控制方法，其特征在于，步骤S2)包括：
S2.1)初始化迭代变量n为1，初始化目标流量
Figure FDA0003326454750000035
以及状态变量；根据nΔt时刻的目标流量
Figure FDA0003326454750000036
以及状态变量，结合式(1)中气液混合系统内空气控制体B的压力pB,gas的计算函数表达式计算得到nΔt时刻的气液混合系统内空气控制体B的压力
Figure FDA0003326454750000037
S2.2)将nΔt时刻的状态变量代入式(1)～(4)，并采用数值差分的方法计算出(n+1)Δt时刻的状态变量；
S2.3)根据(n+1)Δt时刻的状态变量结合式(1)中气液混合系统内空气控制体B的压力pB,gas的计算函数表达式计算得到(n+1)Δt时刻的气液混合系统内空气控制体B的压力
Figure FDA0003326454750000038
S2.4)判断nΔt时刻的气液混合系统内空气控制体B的压力
Figure FDA0003326454750000041
(n+1)Δt时刻的气液混合系统内空气控制体B的压力
Figure FDA0003326454750000042
两者之间的差值的绝对值小于预设阈值是否成立，若成立，则将当前的目标流量
Figure FDA0003326454750000043
作为求解得到的最优的目标流量值
Figure FDA0003326454750000044
输出，结束并退出；否则，根据下式更新当前的目标流量
Figure FDA0003326454750000045
Figure FDA0003326454750000046
式(5)中，
Figure FDA0003326454750000047
为更新后的目标流量，Δ为质量流量的搜索步长，a为质量流量
Figure FDA0003326454750000048
搜索范围的宽度，σ为质量流量
Figure FDA0003326454750000049
搜索范围宽度的缩放因子，m为第m次动态调整质量流量
Figure FDA00033264547500000410
的搜索范围。
6.根据权利要求1所述的发电-储能系统的波动压力发电控制方法，其特征在于，步骤2)中根据测量气压Pmeas、设定气压值Pref之间的差值生成气压误差信号是指将测量气压Pmeas、设定气压值Pref之间的差值输入气压闭环控制器得到气压误差信号。
7.根据权利要求6所述的发电-储能系统的波动压力发电控制方法，其特征在于，步骤2)中根据实测功率值Nmeas、目标功率值Nref之间的差值生成有功功率误差信号是指将实测功率值Nmeas、目标功率值Nref之间的差值输入有功功率闭环控制器得到有功功率误差信号。
8.根据权利要求7所述的发电-储能系统的波动压力发电控制方法，其特征在于，步骤2)中根据目标流量值
Figure FDA00033264547500000411
实测流量值
Figure FDA00033264547500000412
之间的差值生成流量误差信号是指将目标流量值
Figure FDA00033264547500000413
实测流量值
Figure FDA00033264547500000414
之间的差值输入流量闭环控制器得到流量误差信号。
9.一种发电-储能系统的波动压力发电控制系统，包括相互连接的微处理器和存储器，其特征在于，该微处理器被编程或配置以执行权利要求1～8中任意一项所述发电-储能系统的波动压力发电控制方法的步骤。
10.一种计算机可读存储介质，其特征在于，该计算机可读存储介质中存储有被编程或配置以执行权利要求1～8中任意一项所述发电-储能系统的波动压力发电控制方法的计算机程序。
CN202111264203.8A 2021-10-28 2021-10-28 一种发电-储能系统的波动压力发电控制方法及系统 Active CN114094568B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111264203.8A CN114094568B (zh) 2021-10-28 2021-10-28 一种发电-储能系统的波动压力发电控制方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111264203.8A CN114094568B (zh) 2021-10-28 2021-10-28 一种发电-储能系统的波动压力发电控制方法及系统

Publications (2)

Publication Number Publication Date
CN114094568A true CN114094568A (zh) 2022-02-25
CN114094568B CN114094568B (zh) 2023-06-09

Family

ID=80298219

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111264203.8A Active CN114094568B (zh) 2021-10-28 2021-10-28 一种发电-储能系统的波动压力发电控制方法及系统

Country Status (1)

Country Link
CN (1) CN114094568B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023071809A1 (zh) * 2021-10-28 2023-05-04 国网湖南省电力有限公司 一种发电-储能系统的波动压力发电控制方法及装置

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB700098A (en) * 1949-07-29 1953-11-25 British Thomson Houston Co Ltd Improvements in and relating to control systems for air compressor plants
US20050049757A1 (en) * 2003-08-29 2005-03-03 Masami Funakura Cogeneration system, operation controller for cogeneration facility, and operation program for cogeneration facility
US20100089063A1 (en) * 2008-04-09 2010-04-15 Sustainx, Inc. Systems and Methods for Energy Storage and Recovery Using Rapid Isothermal Gas Expansion and Compression
US20100229544A1 (en) * 2009-03-12 2010-09-16 Sustainx, Inc. Systems and Methods for Improving Drivetrain Efficiency for Compressed Gas Energy Storage
CN103791690A (zh) * 2014-02-20 2014-05-14 北京市燃气集团有限责任公司 利用管道压力能发电-制冷的液化天然气生产装置及方法
CN106677966A (zh) * 2016-12-28 2017-05-17 华北电力大学 气液两相联合储能发电系统及其储能发电方法
US20200332709A1 (en) * 2018-01-12 2020-10-22 Kabushiki Kaisha Kobe Seiko Sho (Kobe Steel, Ltd.) Cold-heat power generation device
CN211975319U (zh) * 2020-03-31 2020-11-20 国网湖南省电力有限公司 一种发电系统
CN112596394A (zh) * 2020-12-28 2021-04-02 华电郑州机械设计研究院有限公司 热电联产机组电、热负荷调节的协调控制方法及系统

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB700098A (en) * 1949-07-29 1953-11-25 British Thomson Houston Co Ltd Improvements in and relating to control systems for air compressor plants
US20050049757A1 (en) * 2003-08-29 2005-03-03 Masami Funakura Cogeneration system, operation controller for cogeneration facility, and operation program for cogeneration facility
US20100089063A1 (en) * 2008-04-09 2010-04-15 Sustainx, Inc. Systems and Methods for Energy Storage and Recovery Using Rapid Isothermal Gas Expansion and Compression
US20100229544A1 (en) * 2009-03-12 2010-09-16 Sustainx, Inc. Systems and Methods for Improving Drivetrain Efficiency for Compressed Gas Energy Storage
CN103791690A (zh) * 2014-02-20 2014-05-14 北京市燃气集团有限责任公司 利用管道压力能发电-制冷的液化天然气生产装置及方法
CN106677966A (zh) * 2016-12-28 2017-05-17 华北电力大学 气液两相联合储能发电系统及其储能发电方法
US20200332709A1 (en) * 2018-01-12 2020-10-22 Kabushiki Kaisha Kobe Seiko Sho (Kobe Steel, Ltd.) Cold-heat power generation device
CN211975319U (zh) * 2020-03-31 2020-11-20 国网湖南省电力有限公司 一种发电系统
CN112596394A (zh) * 2020-12-28 2021-04-02 华电郑州机械设计研究院有限公司 热电联产机组电、热负荷调节的协调控制方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张涛;吕卫国;: "三峡地下电站700MW水轮发电机蒸发冷却系统环流、温升与稳定性提高", 水电站机电技术, no. 04 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023071809A1 (zh) * 2021-10-28 2023-05-04 国网湖南省电力有限公司 一种发电-储能系统的波动压力发电控制方法及装置

Also Published As

Publication number Publication date
CN114094568B (zh) 2023-06-09

Similar Documents

Publication Publication Date Title
CN107168055B (zh) 一种基于多变量预测控制的氨法脱硫优化控制方法
CN102608911B (zh) 一种基于多参数预测的火电机组协调控制方法
EP2939726B1 (en) Carbon dioxide capturing system and method of operating same
CN105114242A (zh) 基于模糊自适应的dfpso算法的水轮机调速器参数优化方法
CN101587357A (zh) 温度控制方法及其温度控制系统
CN114094568A (zh) 一种发电-储能系统的波动压力发电控制方法及系统
CN105258099B (zh) 一种加热炉汽包液位控制系统及控制方法
CN112147891A (zh) 一种火电机组协调系统全局非线性优化控制方法
CN114336673B (zh) 一种基于模型预测控制的风储联合电站一次调频控制策略
CN102543232A (zh) 压水堆核电站稳压器水位和压力复合控制方法
CN113904588B (zh) 一种发电-储能系统的波动压力发电控制方法及装置
CN102736640A (zh) 用于除氧器水位控制的系统和方法
JP4841848B2 (ja) ポンプの最適運転方法、情報処理システム、ポンプの最適運転プログラム、複数ポンプの最適流量配分方法
CN114059084A (zh) 电解制氢系统及其温度控制方法和装置
CN112713288A (zh) 一种燃料电池鼓泡加湿器控制系统及控制方法
CN113739447B (zh) 变频热泵温度自适应控制方法、装置、设备及存储介质
Dibiasio et al. Phase‐plane analysis of feedback control of unstable steady states in a biological reactor
Yue et al. Active disturbance rejection control strategy for PEMFC oxygen excess ratio based on adaptive internal state estimation using unscented Kalman filter
CN102723468A (zh) 酸循环极板化成方法及其设备
Gupta Control of integrating processes using dynamic matrix control
CN112563541A (zh) 一种改进粒子群pid的燃料电池阴极压力控制方法
CN107763890A (zh) 一种基于高压储液罐控制的高温热泵系统及控制方法
CN112632757A (zh) 一种基于混合模糊控制的质子交换膜燃料电池的过氧量实时切换控制装置设计方法
Ouyang et al. Optimization and realization of the continuous reactor with improved automatic disturbance rejection control
CN112653157B (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