CN113065244B - 一种火电储能联合调频磨损的量化方法 - Google Patents

一种火电储能联合调频磨损的量化方法 Download PDF

Info

Publication number
CN113065244B
CN113065244B CN202110306811.4A CN202110306811A CN113065244B CN 113065244 B CN113065244 B CN 113065244B CN 202110306811 A CN202110306811 A CN 202110306811A CN 113065244 B CN113065244 B CN 113065244B
Authority
CN
China
Prior art keywords
frequency modulation
thermal power
abrasion
energy storage
power
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
Application number
CN202110306811.4A
Other languages
English (en)
Other versions
CN113065244A (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.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
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 Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN202110306811.4A priority Critical patent/CN113065244B/zh
Publication of CN113065244A publication Critical patent/CN113065244A/zh
Application granted granted Critical
Publication of CN113065244B publication Critical patent/CN113065244B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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/24—Arrangements for preventing or reducing oscillations of power in networks
    • 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
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/04—Power grid distribution networks
    • 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]
    • 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
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/40—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation wherein a plurality of decentralised, dispersed or local energy generation technologies are operated simultaneously

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明属于电气工程领域，具体涉及一种火电储能联合调频磨损的量化方法，包括：构建火电机组‑储能联合调频模型，基于量化时段内的实际负荷数据仿真，将模型中油动机传递函数的输出量在待量化时段内积分得到该时段对应的油动机活塞与缸壁之间的磨损行程；采用磨损行程和材料损耗体积之间的计算模型，求取该调频过程中的材料损耗体积，并结合长期磨损区域表面积得到该调频过程中的缸壁内直径增加量；求解该增加量与表征油动机机缸报废的内直径增加阈值的比值，将该比值与各调频机构元件成本求乘积，实现火电机组调频磨损的量化。本发明通过研究联合调频火电机组自身的磨损降低，有效解决了储能电站辅助火电机组调频的经济性难以衡量的问题。

Description

一种火电储能联合调频磨损的量化方法
技术领域
本发明属于电气工程领域，更具体地，涉及一种火电储能联合调频磨损的量化方法。
背景技术
新能源机组的大规模接入给电力系统带来较大的消纳压力，新能源的出力具有波动性与随机性，容易引发电网的频率稳定问题。目前，我国的调频机组仍旧为火电机组，受到其响应偏差大、响应时间长的固有特性限制，调频性能不够理想。为了应对这一挑战，储能系统被广泛用于改善电网的电能质量。但由于储能造价成本高，其经济性难以衡量。绝大多数研究从储能系统入手，通过减少储能寿命的衰减以提升储能电站的经济性，却忽略了储能电站在降低火电机组自身调频机构磨损方面的重要作用。
由于调频过程中调速系统需要不断动作来响应外部负荷的频繁变化，长期的调频动作会对火电机组的油动机以及调节汽门造成严重磨损，从而造成元件损坏引发非计划停机。因此，为了储能电站更好地推广使用，保障电力系统的安全稳定运行，迫切需要一种从火电机组自身调频机构出发，量化火电机组元件由于调频产生磨损的方法，进而实现联合调频经济效益的量化。
发明内容
针对现有技术的缺陷和改进需求，本发明提供了一种火电储能联合调频磨损的量化方法，旨在解决当储能电站辅助火电机组进行调频时火电机组自身磨损降低的经济性难以衡量的问题。
为实现上述目的，按照本发明的一个方面，提供了一种火电储能联合调频磨损的量化方法，包括：
构建火电机组-储能联合调频模型,用于在外部负荷实时功率输入下联合调整系统实时输出功率；基于量化时段内的实际负荷数据进行仿真，将模型中的油动机传递函数输出量在待量化时段内积分，计算得到该时段对应的调频过程中油动机活塞与缸壁之间的磨损行程；
采用已建立的磨损行程和材料损耗体积之间的关系，计算该调频过程中油动机缸壁造成的材料损耗体积，并结合长期磨损区域表面积得到该调频过程中摩擦导致的缸壁内直径增加量；求解该增加量与表征油动机机缸报废的内缸壁直径阈值的比值，作为经过该调频行为后油动机机缸的损耗系数，将该损耗系数与各调频机构元件成本进行乘积运算，实现各调频机构元件磨损的量化。
进一步，通过Archard磨损模型建立所述磨损行程和材料损耗体积之间的关系，则材料损耗体积
Figure BDA0002987882700000021
式中，K表示摩擦系数，物理意义为单位负荷作用下滑动单位距离所引起的体积磨损；N为摩擦物体所受到的法向载荷；L为磨损行程；H为材料硬度。
进一步，所述直径增加量表示为：
Figure BDA0002987882700000022
式中，V为材料损耗体积，S为长期磨损区域表面积。
进一步，所述火电机组-储能联合调频模型进行联合调整系统实时输出功率的具体方式为：
当储能电站调频仿真模型的功率发出时为正、吸收时为负，与火电机组调频仿真模型发出的功率进行加和，作为联合调频模型的实时输出功率；该输出功率与外部输入的负荷实时功率求差，作为系统实时功率缺额；该实时功率缺额经过发电机转子传递函数得到系统频率波动量；该波动量经过数据传输通道合成系统出力变化量的指令，通过联合调频仿真模型对火电机组和储能电站进行出力变化量的分配，改变系统的实时输出功率，实现联合系统实时输出功率的调节；其中，所述储能电站调频仿真模型通过采用编程方式模拟储能电站控制策略建模得到。
进一步，所述乘积运算包括：
将所述损耗系数与油动机机缸更换成本求解乘积，作为待量化时段内的调频行为后油动机损耗成本。
进一步，所述乘积运算还包括：
通过油动机损耗系数来间接衡量汽门阀杆的机械损耗，获得汽门阀杆的磨损成本CFG，表示为：
Figure BDA0002987882700000031
式中，CFGT表示汽门阀杆的总成本，ηFGoil表示火电厂更换汽门阀杆和油动机的平均个数比，D表示所述直径增加量，Dthreshold表示所述表征油动机机缸报废的内缸壁直径阈值。
本发明还提供一种火电储能联合调频磨损的量化设备，包括处理器和存储器，所述存储器包含计算机可读指令，所述计算机可读指令在由所述处理器执行时使所述量化设备执行如上所述的一种火电储能联合调频磨损的量化方法。
本发明还提供一种计算机可读存储介质，所述计算机可读存储介质包括存储的计算机程序，其中，在所述计算机程序被处理器运行时控制所述存储介质所在设备执行如上所述的一种火电储能联合调频磨损的量化方法。
总体而言，通过本发明所构思的以上技术方案，能够取得以下有益效果：
(1)储能电站辅助火电机组调频能够极大的改善系统电能质量，但由于储能造价成本高、使用寿命有限，联合调频的经济效益难以衡量，阻碍了火电机组运营商对储能的投资。当前研究主要集中在如何降低储能设备的寿命损耗，却忽视了火电机组调频元件在动作过程中也会产生相应的磨损，储能电站的辅助作用可以一定程度降低火电机组自身的磨损成本。从火电机组自身元件入手，将储能装置降低的火电机组磨损成本进行量化，为储能电站的成本回收提供了新思路，使联合调频的经济性得以体现，有利于引导储能电站的推广使用，保障电力系统的安全稳定运行。
(2)长期的调频动作使火电机组调频元件老化从而产生元件更换成本，当前研究缺乏将调频动作与元件老化程度的关联研究。将仿真得到的油动机磨损行程结合Archard磨损模型计算磨损体积进而求得缸壁内径增加量后除以内径增加阈值，可以实现对火电机组调频元件磨损程度的量化，并将此量化量记为损耗系数，即元件使用寿命的损失程度。当损耗系数达到100％时元件需要进行更换从而造成更换成本，因此将损耗系数乘以元件的更换成本可以作为火电机组调频过程中的磨损成本。
(3)对建立的联合调频仿真模型中油动机的输出量进行积分可以得到油动机油缸内壁和活塞的磨损行程，如何将磨损行程转化为元件磨损程度是磨损量化的关键。通过借助摩擦学中经典的Archard磨损模型，可以在已知磨损行程的情况下，结合金属材料其他已知参数求得磨损体积，进而对磨损程度进行量化。这种跨学科的创新可以实现对火电机组元件的磨损成本清晰的量化计算。
(4)汽门阀杆的动作依靠油动机活塞的上下位移带动，因此阀杆与油动机的磨损老化程度存在关联性。通过定义比例系数ηFGoil表征火电厂更换汽门阀杆和油动机的平均个数比，将此系数乘以油动机的损耗系数，即可得到汽门阀杆的损耗系数。
附图说明
图1为本发明实施例提供的一种火电储能联合调频磨损的量化方法框图；
图2为本发明实施例提供的汽轮机传递函数图；
图3是本发明实施例提供的联合调频思路示意图；
图4是本发明实施例提供的联合调频仿真模型控制框图；
图5是本发明实施例提供的火电机组调频机构结构图；
图6是本发明实施例提供的一种能油动机的结构示意图；
图7是本发明实施例提供的通过仿真获得的实时频率图；
图8是本发明实施例提供的通过仿真获得的油动机磨损行程图；
图9是本发明实施例提供的通过仿真获得的油动机磨损成本图；
图10是本发明实施例提供的通过仿真获得的汽门磨损成本图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
实施例一
一种火电储能联合调频磨损的量化方法，如图1所示，包括：
构建火电机组-储能联合调频模型,用于在外部负荷实时功率输入下联合调整系统实时输出功率；基于量化时段内的实际负荷数据进行仿真，将模型中的油动机传递函数输出量在待量化时段内积分，计算得到该时段对应的调频过程中油动机活塞与缸壁之间的磨损行程；
采用已建立的磨损行程和材料损耗体积之间的关系，计算该调频过程中油动机缸壁造成的材料损耗体积，并结合长期磨损区域表面积得到该调频过程中摩擦导致的缸壁内直径增加量；
求解该增加量与表征油动机机缸报废的内缸壁直径阈值的比值，作为经过该调频行为后油动机机缸的损耗系数，将该损耗系数与各调频机构元件成本进行乘积运算，实现各调频机构元件磨损的量化。
其中，在火电机组调频仿真模型基础上添加储能电站调频仿真模型，构成火电机组-储能联合调频模型，用于在外部负荷实时功率输入下联合调整系统实时输出功率，具体的：
(1)关于火电机组的建模，采用传统模型，如下：
力矩马达是一个电气-机械转换器，外部电压信号ug经过放大器产生差动控制电流Δi，从而输出电磁力矩Td，使衔铁产生偏转。力矩马达电路的基本电压方程如下式所示。
Figure BDA0002987882700000061
式中，Ku表示伺服放大器的增益，ug表示伺服放大器的输入电压信号，Rc表示线圈的电阻，rp表示线圈回路中放大器的内阻，Kb表示线圈反动电动势常数，Lc表示线圈的自感系数，Δi表示输入力矩马达的控制电流，θ表示衔铁的转角。
根据力矩马达中存在的力矩平衡关系对力矩马达的电压方程进行简化变形，引入由机械阻尼和电磁阻尼产生的阻尼比，以及忽略力矩马达控制线圈的动态和滑阀的动态，可以得到伺服阀以Ug(s)为输入、以阀芯位移Xv(s)为输出的传递函数。通常情况下伺服阀二阶环节的固有频率ωmf很高，阻尼很小，则其传递函数可简化为一阶惯性，设Ty是伺服阀时间常数，无量纲化后伺服阀的传递函数可以写为：
Figure BDA0002987882700000062
油动机的进油量可以用其体积变化描述，如下式所示。
Figure BDA0002987882700000063
其中，μ为油口的流量系数，XV为滑阀的阀芯位移，b为油口的宽度，ρ为压力油的密度，P0为伺服阀的进油压力，Pd为伺服阀的出油压力。A为进油侧活塞有效面积，Z为油动机活塞行程变化量。
通过对上式进行拉普拉斯变换，并引入油动机的时间常数Tc，意义为油动机活塞在最大进油量条件的行程时间，得到油动机模型的传递函数为：
Figure BDA0002987882700000071
功率比较大、蒸汽参数比较高的汽轮机，由于容积比较大，储存在蒸汽容积中的蒸汽质量比较大，容积效应的影响显著。因此汽轮机应分为蒸汽容积和转子两个环节进行建模。
首先分析蒸汽容积的动态过程。qm1为通过节流阀的蒸汽流量，qm1的大小决定于调节汽阀开度Sz；qm2是流入汽轮机的蒸汽流量，qm2取决于调节级的压力P1，而压力变化取决于容积内积存工质的多少。
(qm2-qm1)dt＝Vdρ
根据调节汽阀的特性和汽轮机的流量特性对上式进行变形，并用小偏差线性化的方法将其转化为一个近似的线性方程。将方程中的各个变量用额定值的百分比表示并进行化简，可以得到输入量为高压调阀开度，输出为容积的压力，中间在热机组喷嘴室进汽容积方程为：
Figure BDA0002987882700000072
式中，TRH为再热器容积的时间常数。
以上为蒸汽容积的动态过程分析，接下来进行转子方程的分析。转子的功率平衡方程为：
Figure BDA0002987882700000073
对转子运动方程进行小偏差线性化，省略高阶项并进行一系列化简后，可将转子的运动方程变形为如下形式：
α1χP1(s)+α2χP2(s)-χL(s)-βχn(s)＝Tasχn(s)
式中，α1为高压缸所占的功率比，α2为中、低压缸所占的功率比，χP1为喷嘴室容积中压力的相对变化量，χP2为中间再热容积中压力的相对变化量。Ta称为转子的飞升时间常数，对于发电机用汽轮机Ta＝6～15s。β称为自平衡系数，是一无量纲量，表示与附加损耗有关的系数，一般为0.03～0.05，χL(s)为负荷扰动相对变化量，χn(s)表示转速相对变化量。
通过以上分析，即可搭建出火电机组调频的仿真模型。
(2)关于储能部分的建模，优选的，根据实际研究的火电厂配置的储能装置所采取的控制策略进行的建模，是用编程的方式模拟控制策略进行建模。
通过编程对现场运行中的储能控制策略进行仿真，并将参数设置为可以手动调节。如图4所示，储能电站辅助调频时相当于一个功率源，当储能电站的功率发出时为正，吸收时为负，与两台火电机组发出的功率进行加和即为系统的输出功率。系统的输出功率与外部负荷的实时功率进行求差可得系统实时的功率缺额。实时功率缺额经过发电机和负荷的调频模型可得系统频率的波动量。系统频率的波动量经过数据传输通道合成系统出力变化量的指令，通过控制策略对火电机组和储能电站进行分配，以改变此时系统的输出功率。整体的建模思路如图3所示，按照以上思路，本实施例搭建了储能电站和火电机组联合调频的仿真模型。其中，图2为汽轮机传递函数图，汽轮机为图4中的从汽门压力到发电机和负荷调频模型(即发电机转子传递函数)部分。图2中，χSZ表示调节阀位移相对变化量；T0称为容积时间常数，表示进汽容积效应对进汽的影响；α1为高压缸所占的功率比；α2为中、低压缸所占的功率比；TRH为再热器容积的时间常数；χP1为喷嘴室容积中压力的相对变化量；χP2为中间再热容积中压力的相对变化量；Ta称为转子的飞升时间常数；β称为自平衡系数，是一无量纲量，表示与附加损耗有关的系数；χL(s)为负荷扰动相对变化量，χn(s)表示转速相对变化量。图4中，R为火电机组的下垂控制系数；T1表示DEH系统的时间常数；Toil表示油动机的时间常数；K1为高压缸所占的功率比；K2为中、低压缸所占的功率比；Th，Tm，Ts分别表示高压缸、中压缸及低压缸的时间常数；M为转子的飞升时间常数；D为自平衡系数；ΔPL、ΔPG、ΔPess分别为负荷、发电机输出功率及储能输出功率的波动量；Δf为系统频率的变化量。
在构建联合调频的仿真模型之后，基于联合调频仿真模型通过带入实际负荷实际数据进行仿真分析，并在油动机传递函数
Figure BDA0002987882700000091
后加入积分环节进行油动机活塞与滑阀相对位移的累加计算，即可获取调频过程中油动机活塞与缸壁间的磨损行程。
关于材料损耗体积的计算。油动机是将油压信号转换为实际阀门开度变化的重要部件，在调频过程中，油动机活塞需要不断的进往复运动从而控制控制汽轮机阀门的开度变化，调频机构如图5所示，图5中的油动机的具体结构图如图6所示。因此，油动机的活塞和油动机的内缸壁之间由于长期的高速相对运动会产生金属材料损耗，金属材料损耗用磨损体积加以衡量。长期的磨损使油动机缸壁内直径增大，在固定磨损区域与其他区域间形成内壁的直径差，造成油动机在调频过程中的卡涩行为，导致元件损坏。
优选的，磨损行程与油动机缸壁材料的损耗体积之间的关系可以通过经典的Archard磨损模型建立，由此磨损模型可得材料损耗体积和滑动摩擦持续的行程之间的关系如下式所示：
Figure BDA0002987882700000092
式中，K表示摩擦系数，物理意义为单位负荷作用下滑动单位距离所引起的体积磨损；N为摩擦物体所受到的法向载荷；L为磨损行程；H为材料硬度。
通过将前述仿真模型获得的油动机活塞与缸壁的磨损行程带入上式进行计算，即可求得调频过程中油动机缸壁由于磨损造成的材料损耗体积。
关于磨损成本的量化。优选的，通过前述计算得到的材料损耗体积以及长期磨损区域的表面积可以得到长期摩擦导致的缸壁内直径的增加量。直径的增加量计算如下式所示：
Figure BDA0002987882700000101
其中，V为磨损体积，S为长期磨损区域表面积。得到内缸壁直径的增加量后，与表征油动机机缸报废的内缸壁直径阈值求解比值可得经过该调频行为后油动机机缸的损耗情况，将此比值与更换油动机机缸所造成的成本求解乘积即为调频行为对火电机组造成的经济损失。
基于以上分析，油动机在调频过程中的磨损成本计算公式为：
Figure BDA0002987882700000102
式中，Coil和CoilT表示调频过程油动机损耗成本和油动机总成本，Dthreshold表示油动机发生活塞卡塞的缸壁直径增量阈值，一般取1.4mm，D表示一个计费周期内的缸壁直径增量。
由于汽门阀杆直接受油动机影响，实现对汽门开度的大小调节，通过油动机活塞行程磨损造成油动机更换来间接衡量汽门阀杆的机械损耗，获得阀杆的磨损成本：
Figure BDA0002987882700000103
式中，CFG和CFGT表示调频汽门阀杆损耗成本和汽门阀杆总成本，ηFGoil表示某个时段内(比如一年内)火电厂更换汽门阀杆和油动机的平均个数比。
通过本实施例所构思的以上技术方案，与现有技术相比，公开了一种火电储能联合调频磨损的量化方法。基于调频机构的动作机理以及金属摩擦损耗的磨损模型对火电机组元件的磨损情况进行了量化。结合元件更换成本得到了不同元件的磨损成本计算公式。
实施案例
以贵州某电厂2019年9月7日的系统实时运行工况和实际数据为例，基于电厂实际火电储能联合调频控制策略，通过MATLAB仿真，分析火电机组单独调频与火电储能联合调频时，调频机构元件磨损成本的差异。
火电单独调频和火电储能联合调频模式下，系统日频率波动曲线如图7所示。火电储能联合调频情况下，系统频率低于火电单独调频模式，储能辅助火电调频的优势得以体现，系统频率毛刺变小、频率波动幅度变小，系统频率偏差变小，频率更为平稳。
火电单独调频和火电储能联合调频油动机活塞行程如图8所示。可以发现，火电储能联合调频情况下活塞行程低于火电单独调频的情况。火电储能联合调频情况下油动机活塞的变化更小，有利于减少对油动机缸壁的磨损。假设火电厂油动机为70万元/个，平均2年进行一更换，调频过程油动机寿命损耗成本如图9所示。可以发现，火电储能联合调频的油动机寿命损耗成本低于火电单独调频的情况。以2019年9月7日为例，该电厂火电单独调频对油动机的损耗约910元/日，折合33.215万元/年；火电储能联合调频对油动机的损耗约510元/日，折合18.615万元/年；储能辅助火电调频将节省油动机设备成本14.5万元/年。
假设火电机组阀杆2.5万元/套，平均每年更换2套，火电单独调频和火电储能联合调频情况下，汽门阀杆的寿命成本如图10所示。以2019年9月7日为例，清水河电厂火电单独调频对汽门阀杆的损耗约180元/日，折合6.57万元/年；火电储能联合调频对汽门阀杆的损耗约100元/日，折合3.65万元/年；储能辅助火电调频将节省汽门阀杆设备成本2.92万元/年。
综上，本实施例公开的一种火电机组调频的磨损量化方法，基于火电机组调频机构的动作机理以及金属材料的摩擦损耗原理，对火电机组调频机构在调频工况下的磨损进行了量化分析。首先根据火电机组各部件的物理提醒进行仿真模型的搭建，并基于储能电站的控制策略建立了火电储能联合调频模型，然后根据Archard磨损模型进行了磨损体积的计算，基于易损件的物理参数对磨损程度进行了量化。通过分析储能电站辅助号调频与火电机组单独调频时的磨损成本的差异，得到了储能电站在降低火电机组调频磨损方面的重要作用。本实施例方法为回收储能电站投资成本提供了一种新思路，有利于引导火电机组运营商进行储能电站的投资，保障电力系统的电能质量和运行稳定性。
实施例二
一种火电储能联合调频磨损的量化设备，包括处理器和存储器，所述存储器包含计算机可读指令，所述计算机可读指令在由所述处理器执行时使所述量化设备执行如上实施例一所述的一种火电储能联合调频磨损的量化方法。相关技术方案同实施例一，在此不再赘述。
实施例三
一种计算机可读存储介质，所述计算机可读存储介质包括存储的计算机程序，其中，在所述计算机程序被处理器运行时控制所述存储介质所在设备执行如上实施例一所述的一种火电储能联合调频磨损的量化方法。相关技术方案同实施例一，在此不再赘述。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种火电储能联合调频磨损的量化方法，其特征在于，包括：
构建火电机组-储能联合调频模型，用于在外部负荷实时功率输入下联合调整系统实时输出功率；基于量化时段内的实际负荷数据进行仿真，将模型中的油动机传递函数输出量在待量化时段内积分，计算得到该时段对应的调频过程中油动机活塞与缸壁之间的磨损行程；
采用已建立的磨损行程和材料损耗体积之间的关系，计算该调频过程中油动机缸壁造成的材料损耗体积，并结合长期磨损区域表面积得到该调频过程中摩擦导致的缸壁内直径增加量；求解该增加量与表征油动机机缸报废的内缸壁直径阈值的比值，作为经过该调频行为后油动机机缸的损耗系数，将该损耗系数与各调频机构元件成本进行乘积运算，实现各调频机构元件磨损的量化；
其中，所述火电机组-储能联合调频模型进行联合调整系统实时输出功率的具体方式为：
当储能电站调频仿真模型的功率发出时为正、吸收时为负，与火电机组调频仿真模型发出的功率进行加和，作为联合调频模型的实时输出功率；该输出功率与外部输入的负荷实时功率求差，作为系统实时功率缺额；该实时功率缺额经过发电机转子传递函数得到系统频率波动量；该波动量经过数据传输通道合成系统出力变化量的指令，通过联合调频仿真模型对火电机组和储能电站进行出力变化量的分配，改变系统的实时输出功率，实现联合系统实时输出功率的调节；其中，所述储能电站调频仿真模型通过采用编程方式模拟储能电站控制策略建模得到；
所述该调频行为为所述待量化时段对应的调频过程。
2.根据权利要求1所述的一种火电储能联合调频磨损的量化方法，其特征在于，通过Archard磨损模型建立所述磨损行程和材料损耗体积之间的关系，则材料损耗体积
Figure FDA0003773027370000021
式中，K表示摩擦系数，物理意义为单位负荷作用下滑动单位距离所引起的体积磨损；N为摩擦物体所受到的法向载荷；L为磨损行程；H为材料硬度。
3.根据权利要求1所述的一种火电储能联合调频磨损的量化方法，其特征在于，所述直径增加量表示为：
Figure FDA0003773027370000022
式中，V为材料损耗体积，S为长期磨损区域表面积。
4.根据权利要求1所述的一种火电储能联合调频磨损的量化方法，其特征在于，所述乘积运算包括：
将所述损耗系数与油动机机缸更换成本求解乘积，作为待量化时段内的调频行为后油动机损耗成本。
5.根据权利要求4所述的一种火电储能联合调频磨损的量化方法，其特征在于，所述乘积运算还包括：
通过油动机的损耗系数来间接衡量汽门阀杆的机械损耗，获得汽门阀杆的磨损成本CFG，表示为：
Figure FDA0003773027370000023
式中，CFGT表示汽门阀杆的总成本，ηFGoil表示火电厂更换汽门阀杆和油动机的平均个数比，D表示所述直径增加量，Dthreshold表示所述表征油动机机缸报废的内缸壁直径阈值。
6.一种火电储能联合调频磨损的量化设备，其特征在于，包括处理器和存储器，所述存储器包含计算机可读指令，所述计算机可读指令在由所述处理器执行时使所述量化设备执行权利要求1至5任一项所述的一种火电储能联合调频磨损的量化方法。
7.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质包括存储的计算机程序，其中，在所述计算机程序被处理器运行时控制所述存储介质所在设备执行如权利要求1至5任一项所述的一种火电储能联合调频磨损的量化方法。
CN202110306811.4A 2021-03-23 2021-03-23 一种火电储能联合调频磨损的量化方法 Active CN113065244B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110306811.4A CN113065244B (zh) 2021-03-23 2021-03-23 一种火电储能联合调频磨损的量化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110306811.4A CN113065244B (zh) 2021-03-23 2021-03-23 一种火电储能联合调频磨损的量化方法

Publications (2)

Publication Number Publication Date
CN113065244A CN113065244A (zh) 2021-07-02
CN113065244B true CN113065244B (zh) 2022-09-20

Family

ID=76562966

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110306811.4A Active CN113065244B (zh) 2021-03-23 2021-03-23 一种火电储能联合调频磨损的量化方法

Country Status (1)

Country Link
CN (1) CN113065244B (zh)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106998079B (zh) * 2017-04-28 2020-05-05 东南大学 一种热电联合优化调度模型的建模方法
CN110021968B (zh) * 2019-04-29 2021-03-23 国网陕西省电力公司电力科学研究院 一种火电机组一次调频控制器、控制方法及系统
CN111756074B (zh) * 2020-06-29 2021-10-08 中国电力科学研究院有限公司 电池储能辅助火电机组深度调峰的优化方法、系统及设备
CN112103980B (zh) * 2020-09-27 2022-04-12 中国科学院电工研究所 一种联合火电机组agc调频的混合储能系统能量管理方法

Also Published As

Publication number Publication date
CN113065244A (zh) 2021-07-02

Similar Documents

Publication Publication Date Title
CN103670536B (zh) 一种火力发电厂汽轮机调门流量的调节方法
CN100545772C (zh) 一种燃煤锅炉系统混合控制方法
CN113489024B (zh) 热电联供机组多模式抽汽辅助调峰调频控制系统及方法
CN101709869B (zh) 燃煤锅炉过热蒸汽温度系统混合控制方法
CN109783936B (zh) 压水堆核电站核岛变工况蓄热增量的计算方法
Xu et al. Stability of speed regulating system of hydropower station with surge tank considering nonlinear turbine characteristics
CN105700380A (zh) 二次再热机组汽轮机调速系统仿真模型及其建模方法
CN104500150B (zh) 一种汽轮机调速器执行机构分段线性模型及参数获取方法
Gao et al. A new linear model of fossil fired steam unit for power system dynamic analysis
CN101709863B (zh) 燃煤锅炉炉膛压力系统混合控制方法
Haji et al. Adaptive model predictive control design for the speed and temperature control of a V94. 2 gas turbine unit in a combined cycle power plant
CN113065244B (zh) 一种火电储能联合调频磨损的量化方法
Ren et al. Application of improved fuzzy immune PID controller to bending control system
CN108131238A (zh) 一种抑制水锤压力波动的新型pid控制方法
CN111927576B (zh) 一种汽轮发电机组空载能耗测试方法
CN110489703B (zh) 先进绝热压缩空气储能热电联供运行可行域分析方法
CN112036009A (zh) 汽轮机调速系统宽工况精细化仿真模型及参数辨识方法
CN105804819B (zh) 一种提高烧结余热回收量及发电系统运行稳定性的调控方法
Albanesi et al. Optimization of the start-up procedure of a combined cycle power plant
CN113914950B (zh) 一种超超临界二次再热多抽汽机组及热电解耦控制方法
CN214660397U (zh) 双低压缸汽轮机低负荷运行时提升供热抽汽参数的装置
CN114963298A (zh) 一种蒸汽管网储能在线定量计算方法和系统
JP5091712B2 (ja) ガバナフリー制御装置およびガバナフリー制御方法
Zhang et al. Modeling and generalized predictive control strategy of a micro gas turbine
CN109696828B (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