CN113765157A - 适用于调频辅助服务市场的发电单元调频性能测算方法和系统 - Google Patents
适用于调频辅助服务市场的发电单元调频性能测算方法和系统 Download PDFInfo
- Publication number
- CN113765157A CN113765157A CN202010488689.2A CN202010488689A CN113765157A CN 113765157 A CN113765157 A CN 113765157A CN 202010488689 A CN202010488689 A CN 202010488689A CN 113765157 A CN113765157 A CN 113765157A
- Authority
- CN
- China
- Prior art keywords
- power generation
- generation unit
- frequency modulation
- time
- calculating
- 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
Links
- 238000010248 power generation Methods 0.000 title claims abstract description 218
- 238000000034 method Methods 0.000 title claims abstract description 130
- 230000008569 process Effects 0.000 claims abstract description 70
- 238000004364 calculation method Methods 0.000 claims abstract description 61
- 238000005259 measurement Methods 0.000 claims abstract description 44
- 238000012795 verification Methods 0.000 claims abstract description 14
- 230000004044 response Effects 0.000 claims description 52
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 claims description 8
- 230000001105 regulatory effect Effects 0.000 claims description 5
- 230000009471 action Effects 0.000 claims description 4
- 230000002123 temporal effect Effects 0.000 claims description 4
- 230000006872 improvement Effects 0.000 description 10
- 238000012360 testing method Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000010926 purge Methods 0.000 description 1
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/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/48—Controlling the sharing of the in-phase component
-
- 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/004—Generation forecast, e.g. methods or systems for forecasting future energy generation
-
- 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/007—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
- H02J3/0075—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source according to economic or energy efficiency considerations, e.g. economic dispatch
-
- 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/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明公开了一种适用于调频辅助服务市场的发电单元调频性能测算方法和系统,所述方法包括获取发电单元在出清时段内的原始数据;对发电单元在出清时段内的各个AGC指令调节过程进行校验;对通过校验的各个AGC指令调节过程进行调频性能测算;统计发电单元在出清时段内的调频性能指标,完成发电单元调频性能测算。本发明能够对发电单元在各个出清时段的调频性能进行数值计算,具有良好的准确性,可适用于单机或多机发电单元,可兼容水电、火电、新能源等不同类型发电单元,适合程序化、规范化、定值化地自动测算,计算结果稳定合理。
Description
技术领域
本发明属于电力市场和电力系统控制技术领域,具体涉及一种适用于调频辅助服务市场的发电单元调频性能测算方法和系统,能够适用于调频辅助服务市场。
背景技术
调频是电力系统中重要的辅助服务项目,对于电网安全和电能质量有着重要意义。传统的非市场环境电力系统中,调频辅助服务提供者包括水电厂、火电厂和新能源电厂等,在电厂接入、AGC优化调用和考核结算等环节都需要考虑它们的调频性能。调度机构通常使用手动或自动的方法向发电单元发送测试指令对其调频性能进行测算。
在调频辅助服务市场中,调频市场技术支持系统根据调频辅助服务提供者的报价和调频性能等数据制定出清结果,AGC按照调频市场技术支持系统提供的出清结果,根据电网调频需求向中标或者临时征用的发电单元发送调节指令,从而实现调频资源的调用,根据发电单元调频性能、调频里程等数据对辅助服务提供者进行考核结算。发电单元的调频性能数据在调频辅助服务市场的准入、出清、考核和结算等过程中都至关重要。
在调频辅助服务市场中,调频性能影响着发电单元是否能够中标、能够获得多少收益。以往使用人工或自动进行响应测试的性能测算方法,只对发电单元少数时段或个别指令过程进行测算,不能满足自动对发电单元各出清时段性能进行测算的要求,而使用扰动测试或测试指令进行调频性能测算会影响市场环境下调频辅助服务提供者之间的公平性,因此,传统的性能测算方法并不适用于调频辅助服务市场环境。同时,传统的性能测算方法对于发电单元的调节速率进行测速时,一般按照一个固定周期计算发电单元调节速率,得到的速率指标计算结果不稳定也不准确。市场环境下,需要一种准确的数值计算方法,能够兼容水电、火电、新能源等不同种类的发电单元,在发电单元按照AGC指令提供调频辅助服务时,自动地对发电单元在各出清时段的调频性能进行测算。由于调频市场环境下,发电单元是按照AGC下发的调节指令发挥调频作用的,因此,AGC指令应是发电单元调频性能测算的重要依据。对于很多电全厂发电单元,在AGC指令调节量较小时,会将指令分配给单独机组执行,在AGC调节量较大时,指令以一定策略被分配给多个机组执行,区分不同指令大小情况下发电单元的性能测算方法,才能使测算结果更加稳定合理。
发明内容
针对上述问题,本发明提出一种适用于调频辅助服务市场的发电单元调频性能测算方法和系统,能够对发电单元在各个出清时段的调频性能进行数值计算,具有良好的准确性,可适用于单机或多机发电单元,可兼容水电、火电、新能源等不同类型发电单元。
为了实现上述技术目的,达到上述技术效果,本发明通过以下技术方案实现:
第一方面,本发明提供了一种适用于调频辅助服务市场的发电单元调频性能测算方法,包括:
获取发电单元在出清时段内的原始数据;
对发电单元在出清时段内的各个AGC指令调节过程进行校验;
对通过校验的各个AGC指令调节过程进行调频性能测算;
统计发电单元在出清时段内的调频性能指标,完成发电单元调频性能测算。
作为本发明的进一步改进,所述获取发电单元在出清时段内的原始数据,包括以下步骤:
获取出清时段内发电单元的所有AGC指令数据,包括各个AGC指令的发出时刻、目标有功功率和发出指令时发电单元所处的AGC控制模式;
获取出清时段内发电单元包含的各个机组每秒的实时有功功率;
基于各个机组每秒的实时有功功率,累计得到出清时段内发电单元每秒的实时有功功率。
作为本发明的进一步改进,所述对发电单元在出清时段内的各个AGC指令调节过程进行校验,包括以下步骤:
将各个AGC指令调节过程与预设的校验条件进行比对,当且仅当AGC指令调节过程满足预设的校验条件时通过校验;
所述预设的校验条件包括:
AGC指令发出时刻,发电单元的AGC控制模式为调频模式;
AGC指令发出时的实时有功功率与目标有功功率之间偏差绝对值不小于指令测算门槛设定值。
作为本发明的进一步改进,所述对通过校验的各个AGC指令调节过程进行调频性能测算,包括以下步骤:
计算AGC指令调节过程中的限值参数;
计算AGC指令调节过程中发电单元的调节速率指标、响应速度指标和调节精度指标。
作为本发明的进一步改进,定义Max/Min(a,b)运算符表示根据设定对a和b取大值或小值;记出清时段开始时间为Tstart,结束时间额为Tend;对于按照AGC指令发出时间先后顺序排列的指令调节过程,记发电单元在t时刻的实时有功功率为P(t),第i个AGC指令调节过程中,指令发出时刻为Tstart(i),发出时刻有功功率Pstart(i)=P(Tstart(i)),指令目标有功功率为Ptarget(i),L1为测速最短时长设定值,按照以下方法计算该AGC指令调节过程中的计算限值参数:
计算Tstart(i)时刻的发电单元开机容量Son(i);
计算Tstart(i)时刻的发电单元下开机的最大单机容量Smax(i);
计算指令调节过程结束时间Tend(i)=min(Tstart(i+1),Tend);
计算指令调节量Preg(i)=Ptarget(i)-P(Tstart(i));
计算动作死区Psdeadb(i)=Max/Min(Son(i)×R1,B1),其中E1和B1为设定的与发电单元无关的常数;
计算目标死区Pddeadb(i)=Max/Min(Smax(i)×R2,B2),其中R2和B2为设定的与发电单元无关的常数;
计算指令分担判定门槛Ps(i)=Max/Min(Son(i)×R3,B3),其中R3和B3为设定的与发电单元类型相关的常数;
计算调节速率计算门槛PT(i),若|Preg(i)|<Ps(i),PT(i)=Max/Min(Smax(i)×R3,B3),否则,PT(i)=Max/Min(Son(i)×R3,B3),其中R4和B4为设定的与发电单元类型相关的常数;
计算调节速率起始计算门槛Psd(i),若|Preg(i)|<Ps(i),Psd(i)=Max/Min(Smax(i)×R5,B5),否则,Psd(i)=Max/Min(Son(i)×R5,B5),其中R5和B5为设定的与发电单元类型相关的常数。
作为本发明的进一步改进,所述第i个AGC指令调节过程中发电单元的调节速率指标为k1(i)的计算步骤包括:
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的测速起始时刻t0(i),并记P0(i)=P(t0(i)),若找不到满足所述条件的时刻,则该AGC指令调节过程的调节速率指标k1(i)=-1,结束调节速率指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Pstart(i)|>Psd(i);
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的测速结束时刻t1(i),并记P1(i)=P(t1(i)),若找不到满足所述条件的时刻,则本指令调节过程的调节速率指标k1(i)=-1,结束调节速率指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Pstart(i)|>Preg(i)·D,其中D为与发电单元无关的测速终点调节量门槛设定值;
按照所述方法计算本调节过程的调节速率值V(i),此处所述方法为若|P0(i)-P1(i)|≥PT(i)且t1(i)-t0(i)≥L1,则否则指标k1(i)=-1,结束调节速率指标计算;此处L1为与发电单元无关的测算最小时长设定值;
作为本发明的进一步改进,所述第i个指令调节过程中发电单元的响应速度指标k2(i)的计算步骤包括:
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的时刻,记录为可靠响应时刻Tact(i),若找不到满足所述条件的时刻,则Tact(i)=-1,结束响应速度指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Pstart(i)|>Psdeadb(i)且维持U3秒,此处U3为可靠动作时间门槛,是与发电单元无关的设定常数;
计算此调节过程中的响应时间Tresponse(i)=Tact(i)-Tstart(i);
作为本发明的进一步改进,所述第i个指令调节过程中发电单元的调节精度指标k3(i)的计算步骤包括:
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的时刻,记录为进入调节死区时刻Tdes(i),若找不到满足所述条件的时刻,则k3(i)=-1,结束响应速度指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Ptarget(i)|>Pddeadb(i);
若Tend(i)-Tdes(i)<L2,则k3(i)=-1,结束响应速度指标计算,此处L2为精度测算最小时长,是与发电单元无关的设定常数;
计算此调节过程中的精度测算时长Taccu(i)=min(Tend(i)-Tdes(i),L3),L3为精度测算最大时长,是与发电单元无关的设定常数;
作为本发明的进一步改进,所述统计发电单元在出清时段内的调频性能指标,包括以下步骤:
记第i个AGC指令调节过程中发电单元的调节速率指标为k1(i),响应速度指标为k2(i),调节精度指标为k3(i);
计算发电单元在出清时段内的调节速率性能指标k1为所有大于0的k1(i)的算术平均值,若该出清时段没有大于0的k1(i),则该出清时段内的调节速率性能指标k1为最近一次测得的大于0的出清时段调节速率性能指标;
计算发电单元在出清时段内的响应速度性能指标k2为所有大于0的k2(i)的算术平均值,若该出清时段没有大于0的k2(i),则该出清时段内的响应速度性能指标k2为最近一次测得的大于0的出清时段响应速度性能指标;
计算发电单元在出清时段内的调节精度性能指标k3为所有大于0的k3(i)的算术平均值,若该出清时段没有大于0的k3(i),则该出清时段内的调节精度性能指标k3为最近一次测得的大于0的出清时段调节精度性能指标;
第二方面,本发明提供了一种适用于调频辅助服务市场的发电单元调频性能测算系统,包括:
获取单元,用于获取发电单元在出清时段内的原始数据;
校验单元,用于对发电单元在出清时段内的各个AGC指令调节过程进行校验;
测算单元,用于对通过校验的各个AGC指令调节过程进行调频性能测算;
统计单元,用于统计发电单元在出清时段内的调频性能指标,完成发电单元调频性能测算。
与现有技术相比,本发明的有益效果:
本发明能够被应用于某一大型区域调频辅助服务市场中,实现了市场环境下对发电单元在各个出清时段调频性能的自动化、规范化、定值化测算,具有良好的准确性,可适用于单机或多机发电单元,可兼容水电、火电、新能源等不同类型发电单元,计算结果稳定合理。
附图说明
为了使本发明的内容更容易被清楚地理解,下面根据具体实施例并结合附图,对本发明作进一步详细的说明,其中:
图1为本发明一种实施例的发电单元性能测算整体流程图;
图2为本发明一种实施例的发电单元AGC指令调节过程示意图;
图3为本发明一种实施例的发电单元调节速率计算方法流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明的保护范围。
下面结合附图对本发明的应用原理作详细的描述。
实施例1
本发明实施例中提供了一种适用于调频辅助服务市场的发电单元调频性能测算方法,如图1所示,具体包括以下步骤:
获取发电单元在出清时段内的原始数据;
对发电单元在出清时段内的各个AGC指令调节过程进行校验;
对通过校验的各个AGC指令调节过程进行调频性能测算;
统计发电单元在出清时段内的调频性能指标,完成发电单元调频性能测算。
在本发明实施例的一种具体实施方式中,所述获取发电单元在出清时段内的原始数据,包括以下步骤:
从商用数据库获取出清时段内发电单元的所有AGC指令数据,包括各个AGC指令的发出时刻、目标有功功率和发出指令时发电单元所处的AGC控制模式,图1中CMD表示按照指令发出时刻先后顺序存储的AGC指令数据集合,CMD(i)表示集合中第i条指令数据,card(CMD)表示CMD中的指令数目;
从时序数据库获取出清时段内发电单元包含的各个机组每秒的实时有功功率;
基于各个机组每秒的实时有功功率,累计得到出清时段内发电单元每秒的实时有功功率。
在本发明实施例的一种具体实施方式中,所述对发电单元在出清时段内的各个AGC指令调节过程进行校验,包括以下步骤:
将各个AGC指令调节过程与预设的校验条件进行比对,当且仅当AGC指令调节过程满足预设的校验条件时通过校验;
所述预设的校验条件包括:
AGC指令发出时刻,发电单元的AGC控制模式为调频模式;
AGC指令发出时的实时有功功率与目标有功功率之间偏差绝对值不小于指令测算门槛设定值。
在本发明实施例的一种具体实施方式中,所述对通过校验的各个AGC指令调节过程进行调频性能测算,包括以下步骤:
计算AGC指令调节过程中的限值参数;
计算AGC指令调节过程中发电单元的调节速率指标、响应速度指标和调节精度指标。
作为本发明的进一步改进,定义Max/Min(a,b)运算符表示根据设定对a和b取大值或小值;记出清时段开始时间为Tstart,结束时间额为Tend;对于按照AGC指令发出时间先后顺序排列的指令调节过程,记发电单元在t时刻的实时有功功率为P(t),第i个AGC指令调节过程中,指令发出时刻为Tstart(i),发出时刻有功功率Pstart(i)=P(Tstart(i)),指令目标有功功率为Ptarget(i),L1为测速最短时长设定值,按照以下方法计算该AGC指令调节过程中的计算限值参数:
计算Tstart(i)时刻的发电单元开机容量Son(i);
计算Tstart(i)时刻的发电单元下开机的最大单机容量Smax(i);
计算指令调节过程结束时间Tend(i)=min(Tstart(i+1),Tend);
计算指令调节量Preg(i)=Ptarget(i)-P(Tstart(i));
计算动作死区Psdeadb(i)=Max/Min(Son(i)×R1,B1),其中R1和B1为设定的与发电单元无关的常数;
计算目标死区Pddeadb(i)=Max/Min(Smax(i)×R2,B2),其中R2和B2为设定的与发电单元无关的常数;
计算指令分担判定门槛Ps(i)=Max/Min(Son(i)×R3,B3),其中R3和B3为设定的与发电单元类型相关的常数;
计算调节速率计算门槛PT(i),若|Preg(i)|<Ps(i),PT(i)=Max/Min(Smax(i)×R3,B3),否则,PT(i)=Max/Min(Son(i)×R3,B3),其中R4和B4为设定的与发电单元类型相关的常数;
计算调节速率起始计算门槛Psd(i),若|Preg(i)|<Ps(i),Psd(i)=Max/Min(Smax(i)×R5,B5),否则,Psd(i)=Max/Min(Son(i)×R5,B5),其中R5和B5为设定的与发电单元类型相关的常数。
如图3所示,所述第i个AGC指令调节过程中发电单元的调节速率指标为k1(i)的计算步骤包括:
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的测速起始时刻t0(i),并记P0(i)=P(t0(i)),若找不到满足所述条件的时刻,则该AGC指令调节过程的调节速率指标k1(i)=-1,结束调节速率指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Pstart(i)|>Psd(i);
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的测速结束时刻t1(i),并记P1(i)=P(t1(i)),若找不到满足所述条件的时刻,则本指令调节过程的调节速率指标k1(i)=-1,结束调节速率指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Pstart(i)|>Preg(i)·D,其中D为与发电单元无关的测速终点调节量门槛设定值;
按照所述方法计算本调节过程的调节速率值V(i),此处所述方法为若|P0(i)-P1(i)|≥PT(i)且t1(i)-t0(i)≥L1,则否则指标k1(i)=-1,结束调节速率指标计算;此处L1为与发电单元无关的测算最小时长设定值;
所述第i个指令调节过程中发电单元的响应速度指标k2(i)的计算步骤包括:
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的时刻,记录为可靠响应时刻Tact(i),若找不到满足所述条件的时刻,则Tact(i)=-1,结束响应速度指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Pstart(i)|>Psdeadb(i)且维持U3秒,此处U3为可靠动作时间门槛,是与发电单元无关的设定常数;
计算此调节过程中的响应时间Tresponse(i)=Tact(i)-Tstart(i);
所述第i个指令调节过程中发电单元的调节精度指标k3(i)的计算步骤包括:
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的时刻,记录为进入调节死区时刻Tdes(i),若找不到满足所述条件的时刻,则k3(i)=-1,结束响应速度指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Ptarget(i)|>Pddeadb(i);
若Tend(i)-Tdes(i)<L2,则k3(i)=-1,结束响应速度指标计算,此处L2为精度测算最小时长,是与发电单元无关的设定常数;
计算此调节过程中的精度测算时长Taccu(i)=min(Tend(i)-Tdes(i),L3),L3为精度测算最大时长,是与发电单元无关的设定常数;
在本发明实施例的一种具体实施方式中,所述统计发电单元在出清时段内的调频性能指标,包括以下步骤:
记第i个AGC指令调节过程中发电单元的调节速率指标为k1(i),响应速度指标为k2(i),调节精度指标为k3(i);
计算发电单元在出清时段内的调节速率性能指标k1为所有大于0的k1(i)的算术平均值,若该出清时段没有大于0的k1(i),则该出清时段内的调节速率性能指标k1为最近一次测得的大于0的出清时段调节速率性能指标;
计算发电单元在出清时段内的响应速度性能指标k2为所有大于0的k2(i)的算术平均值,若该出清时段没有大于0的k2(i),则该出清时段内的响应速度性能指标k2为最近一次测得的大于0的出清时段响应速度性能指标;
计算发电单元在出清时段内的调节精度性能指标k3为所有大于0的k3(i)的算术平均值,若该出清时段没有大于0的k3(i),则该出清时段内的调节精度性能指标k3为最近一次测得的大于0的出清时段调节精度性能指标;
实施例2
基于与实施例1相同的发明构思,本发明实施例中提供了一种适用于调频辅助服务市场的发电单元调频性能测算系统,包括:
获取单元,用于获取发电单元在出清时段内的原始数据;
校验单元,用于对发电单元在出清时段内的各个AGC指令调节过程进行校验;
测算单元,用于对通过校验的各个AGC指令调节过程进行调频性能测算;
统计单元,用于统计发电单元在出清时段内的调频性能指标,完成发电单元调频性能测算。
在本发明实施例的一种具体实施方式中,所述获取发电单元在出清时段内的原始数据,包括以下步骤:
从商用数据库获取出清时段内发电单元的所有AGC指令数据,包括各个AGC指令的发出时刻、目标有功功率和发出指令时发电单元所处的AGC控制模式;
从时序数据库获取出清时段内发电单元包含的各个机组每秒的实时有功功率;
基于各个机组每秒的实时有功功率,累计得到出清时段内发电单元每秒的实时有功功率。
在本发明实施例的一种具体实施方式中,所述对发电单元在出清时段内的各个AGC指令调节过程进行校验,包括以下步骤:
将各个AGC指令调节过程与预设的校验条件进行比对,当且仅当AGC指令调节过程满足预设的校验条件时通过校验;
所述预设的校验条件包括:
AGC指令发出时刻,发电单元的AGC控制模式为调频模式;
AGC指令发出时的实时有功功率与目标有功功率之间偏差绝对值不小于指令测算门槛设定值。
在本发明实施例的一种具体实施方式中,所述对通过校验的各个AGC指令调节过程进行调频性能测算,包括以下步骤:
计算AGC指令调节过程中的限值参数;
计算AGC指令调节过程中发电单元的调节速率指标、响应速度指标和调节精度指标。
作为本发明的进一步改进,定义Max/Min(a,b)运算符表示根据设定对a和b取大值或小值;记出清时段开始时间为Tstart,结束时间额为Tend;对于按照AGC指令发出时间先后顺序排列的指令调节过程,记发电单元在t时刻的实时有功功率为P(t),第i个AGC指令调节过程中,指令发出时刻为Tstart(i),发出时刻有功功率Pstart(i)=P(Tstart(i)),指令目标有功功率为Ptarget(i),L1为测速最短时长设定值,按照以下方法计算该AGC指令调节过程中的计算限值参数:
计算Tstart(i)时刻的发电单元开机容量Son(i);
计算Tstart(i)时刻的发电单元下开机的最大单机容量Smax(i);
计算指令调节过程结束时间Tend(i)=min(Tstart(i+1),Tend);
计算指令调节量Preg(i)=Ptarget(i)-P(Tstart(i));
计算动作死区Psdeadb(i)=Max/Min(Son(i)×R1,B1),其中R1和B1为设定的与发电单元无关的常数;
计算目标死区Pddeadb(i)=Max/Min(Smax(i)×R2,B2),其中R2和B2为设定的与发电单元无关的常数;
计算指令分担判定门槛Ps(i)=Max/Min(Son(i)×R3,B3),其中R3和B3为设定的与发电单元类型相关的常数;
计算调节速率计算门槛PT(i),若|Preg(i)|<Ps(i),PT(i)=Max/Min(Smax(i)×R3,B3),否则,PT(i)=Max/Min(Son(i)×R3,B3),其中R4和B4为设定的与发电单元类型相关的常数;
计算调节速率起始计算门槛Psd(i),若|Preg(i)|<Ps(i),Psd(i)=Max/Min(Smax(i)×R5,B5),否则,Psd(i)=Max/Min(Son(i)×R5,B5),其中R5和B5为设定的与发电单元类型相关的常数。
如图3所示,所述第i个AGC指令调节过程中发电单元的调节速率指标为k1(i)的计算步骤包括:
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的测速起始时刻t0(i),并记P0(i)=P(t0(i)),若找不到满足所述条件的时刻,则该AGC指令调节过程的调节速率指标k1(i)=-1,结束调节速率指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Pstart(i)|>Psd(i);
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的测速结束时刻t1(i),并记P1(i)=P(t1(i)),若找不到满足所述条件的时刻,则本指令调节过程的调节速率指标k1(i)=-1,结束调节速率指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Pstart(i)|>Preg(i)·D,其中D为与发电单元无关的测速终点调节量门槛设定值;
按照所述方法计算本调节过程的调节速率值V(i),此处所述方法为若|P0(i)-P1(i)|≥PT(i)且t1(i)-t0(i)≥L1,则否则指标k1(i)=-1,结束调节速率指标计算;此处L1为与发电单元无关的测算最小时长设定值;
所述第i个指令调节过程中发电单元的响应速度指标k2(i)的计算步骤包括:
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的时刻,记录为可靠响应时刻Tact(i),若找不到满足所述条件的时刻,则Tact(i)=-1,结束响应速度指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Pstart(i)|>Psdeadb(i)且维持U3秒,此处U3为可靠动作时间门槛,是与发电单元无关的设定常数;
计算此调节过程中的响应时间Tresponse(i)=Tact(i)-Tstart(i);
所述第i个指令调节过程中发电单元的调节精度指标k3(i)的计算步骤包括:
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的时刻,记录为进入调节死区时刻Tdes(i),若找不到满足所述条件的时刻,则k3(i)=-1,结束响应速度指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Ptarget(i)|>Pddeadb(i);
若Tend(i)-Tdes(i)<L2,则k3(i)=-1,结束响应速度指标计算,此处L2为精度测算最小时长,是与发电单元无关的设定常数;
计算此调节过程中的精度测算时长Taccu(i)=min(Tend(i)-Tdes(i),L3),L3为精度测算最大时长,是与发电单元无关的设定常数;
在本发明实施例的一种具体实施方式中,所述统计发电单元在出清时段内的调频性能指标,包括以下步骤:
记第i个AGC指令调节过程中发电单元的调节速率指标为k1(i),响应速度指标为k2(i),调节精度指标为k3(i);
计算发电单元在出清时段内的调节速率性能指标k1为所有大于0的k1(i)的算术平均值,若该出清时段没有大于0的k1(i),则该出清时段内的调节速率性能指标k1为最近一次测得的大于0的出清时段调节速率性能指标;
计算发电单元在出清时段内的响应速度性能指标k2为所有大于0的k2(i)的算术平均值,若该出清时段没有大于0的k2(i),则该出清时段内的响应速度性能指标k2为最近一次测得的大于0的出清时段响应速度性能指标;
计算发电单元在出清时段内的调节精度性能指标k3为所有大于0的k3(i)的算术平均值,若该出清时段没有大于0的k3(i),则该出清时段内的调节精度性能指标k3为最近一次测得的大于0的出清时段调节精度性能指标;
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (10)
1.一种适用于调频辅助服务市场的发电单元调频性能测算方法,其特征在于,包括:
获取发电单元在出清时段内的原始数据;
对发电单元在出清时段内的各个AGC指令调节过程进行校验;
对通过校验的各个AGC指令调节过程进行调频性能测算;
统计发电单元在出清时段内的调频性能指标,完成发电单元调频性能测算。
2.根据权利要求1所述的一种适用于调频辅助服务市场的发电单元调频性能测算方法,其特征在于,所述获取发电单元在出清时段内的原始数据,包括以下步骤:
获取出清时段内发电单元的所有AGC指令数据,包括各个AGC指令的发出时刻、目标有功功率和发出指令时发电单元所处的AGC控制模式;
获取出清时段内发电单元包含的各个机组每秒的实时有功功率;
基于各个机组每秒的实时有功功率,累计得到出清时段内发电单元每秒的实时有功功率。
3.根据权利要求2所述的一种适用于调频辅助服务市场的发电单元调频性能测算方法,其特征在于,所述对发电单元在出清时段内的各个AGC指令调节过程进行校验,包括以下步骤:
将各个AGC指令调节过程与预设的校验条件进行比对,当且仅当AGC指令调节过程满足预设的校验条件时通过校验;
所述预设的校验条件包括:
AGC指令发出时刻,发电单元的AGC控制模式为调频模式;
AGC指令发出时的实时有功功率与目标有功功率之间偏差绝对值不小于指令测算门槛设定值。
4.根据权利要求1所述的一种适用于调频辅助服务市场的发电单元调频性能测算方法,其特征在于,所述对通过校验的各个AGC指令调节过程进行调频性能测算,包括以下步骤:
计算AGC指令调节过程中的限值参数;
计算AGC指令调节过程中发电单元的调节速率指标、响应速度指标和调节精度指标。
5.根据权利要求4所述的一种适用于调频辅助服务市场的发电单元调频性能测算方法,其特征在于:定义Max/Min(a,b)运算符表示根据设定对a和b取大值或小值;记出清时段开始时间为Tstart,结束时间额为Tend;对于按照AGC指令发出时间先后顺序排列的指令调节过程,记发电单元在t时刻的实时有功功率为P(t),第i个AGC指令调节过程中,指令发出时刻为Tstart(i),发出时刻有功功率Pstart(i)=P(Tstart(i)),指令目标有功功率为Ptarget(i),L1为测速最短时长设定值,按照以下方法计算该AGC指令调节过程中的计算限值参数:
计算Tstart(i)时刻的发电单元开机容量Son(i);
计算Tstart(i)时刻的发电单元下开机的最大单机容量Smax(i);
计算指令调节过程结束时间Tend(i)=min(Tstart(i+1),Tend);
计算指令调节量Preg(i)=Ptarget(i)-P(Tstart(i));
计算动作死区Psdeadb(i)=Max/Min(Son(i)×R1,B1),其中R1和B1为设定的与发电单元无关的常数;
计算目标死区Pddeadb(i)=Max/Min(Smax(i)×R2,B2),其中R2和B2为设定的与发电单元无关的常数;
计算指令分担判定门槛Ps(i)=Max/Min(Son(i)×R3,B3),其中R3和B3为设定的与发电单元类型相关的常数;
计算调节速率计算门槛PT(i),若|Preg(i)|<Ps(i),PT(i)=Max/Min(Smax(i)×R3,B3),
否则,PT(i)=Max/Min(Son(i)×R3,B3),其中R4和B4为设定的与发电单元类型相关的常数;
计算调节速率起始计算门槛Psd(i),若|Preg(i)|<Ps(i),Psd(i)=Max/Min(Smax(i)×R5,B5),否则,Psd(i)=Max/Min(Son(i)×R5,B5),其中R5和B5为设定的与发电单元类型相关的常数。
6.根据权利要求5所述的一种适用于调频辅助服务市场的发电单元调频性能测算方法,其特征在于:所述第i个AGC指令调节过程中发电单元的调节速率指标为k1(i)的计算步骤包括:
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的测速起始时刻t0(i),并记P0(i)=P(t0(i)),若找不到满足所述条件的时刻,则该AGC指令调节过程的调节速率指标k1(i)=-1,结束调节速率指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Pstart(i)|>Psd(i);
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的测速结束时刻t1(i),并记P1(i)=P(t1(i)),若找不到满足所述条件的时刻,则本指令调节过程的调节速率指标k1(i)=-1,结束调节速率指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Pstart(i)|>Preg(i)·D,其中D为与发电单元无关的测速终点调节量门槛设定值;
按照所述方法计算本调节过程的调节速率值V(i),此处所述方法为若|P0(i)-P1(i)|≥PT(i)且t1(i)-t0(i)≥L1,则否则指标k1(i)=-1,结束调节速率指标计算;此处L1为与发电单元无关的测算最小时长设定值;
7.根据权利要求5所述的一种适用于调频辅助服务市场的发电单元调频性能测算方法,其特征在于:所述第i个指令调节过程中发电单元的响应速度指标k2(i)的计算步骤包括:从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的时刻,记录为可靠响应时刻Tact(i),若找不到满足所述条件的时刻,则Tact(i)=-1,结束响应速度指标计算;此处所述条件为Pref(i)·(P(t)-Pstart(i))>0且|P(t)-Pstart(i)|>Psdeadb(i)且维持U3秒,此处U3为可靠动作时间门槛,是与发电单元无关的设定常数;
计算此调节过程中的响应时间Tresponse(i)=Tact(i)-Tstart(i);
8.根据权利要求5所述的一种适用于调频辅助服务市场的发电单元调频性能测算方法,其特征在于:所述第i个指令调节过程中发电单元的调节精度指标k3(i)的计算步骤包括:
从指令发出时刻Tstart(i)开始到Tend(i),按时间顺序查找发电单元实时出力P(t)首次满足所述条件的时刻,记录为进入调节死区时刻Tdes(i),若找不到满足所述条件的时刻,则k3(i)=-1,结束响应速度指标计算;此处所述条件为Preg(i)·(P(t)-Pstart(i))>0且|P(t)-Ptarget(i)|>Pddeadb(i);
若Tend(i)-Tdes(i)<L2,则k3(i)=-1,结束响应速度指标计算,此处L2为精度测算最小时长,是与发电单元无关的设定常数;
计算此调节过程中的精度测算时长Taccu(i)=min(Tend(i)-Tdes(i),L3),L3为精度测算最大时长,是与发电单元无关的设定常数;
9.根据权利要求4所述的一种适用于调频辅助服务市场的发电单元调频性能测算方法,其特征在于:所述统计发电单元在出清时段内的调频性能指标,包括以下步骤:
记第i个AGC指令调节过程中发电单元的调节速率指标为k1(i),响应速度指标为k2(i),调节精度指标为k3(i);
计算发电单元在出清时段内的调节速率性能指标k1为所有大于0的k1(i)的算术平均值,若该出清时段没有大于0的k1(i),则该出清时段内的调节速率性能指标k1为最近一次测得的大于0的出清时段调节速率性能指标;
计算发电单元在出清时段内的响应速度性能指标k2为所有大于0的k2(i)的算术平均值,若该出清时段没有大于0的k2(i),则该出清时段内的响应速度性能指标k2为最近一次测得的大于0的出清时段响应速度性能指标;
计算发电单元在出清时段内的调节精度性能指标k3为所有大于0的k3(i)的算术平均值,若该出清时段没有大于0的k3(i),则该出清时段内的调节精度性能指标k3为最近一次测得的大于0的出清时段调节精度性能指标;
10.一种适用于调频辅助服务市场的发电单元调频性能测算系统,其特征在于,包括:
获取单元,用于获取发电单元在出清时段内的原始数据;
校验单元,用于对发电单元在出清时段内的各个AGC指令调节过程进行校验;
测算单元,用于对通过校验的各个AGC指令调节过程进行调频性能测算;
统计单元,用于统计发电单元在出清时段内的调频性能指标,完成发电单元调频性能测算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010488689.2A CN113765157B (zh) | 2020-06-02 | 2020-06-02 | 适用于调频辅助服务市场的发电单元调频性能测算方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010488689.2A CN113765157B (zh) | 2020-06-02 | 2020-06-02 | 适用于调频辅助服务市场的发电单元调频性能测算方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113765157A true CN113765157A (zh) | 2021-12-07 |
CN113765157B CN113765157B (zh) | 2024-07-30 |
Family
ID=78782727
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010488689.2A Active CN113765157B (zh) | 2020-06-02 | 2020-06-02 | 适用于调频辅助服务市场的发电单元调频性能测算方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113765157B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105207242A (zh) * | 2015-09-17 | 2015-12-30 | 山东大学 | 储能装置参与机组调频的优化控制与容量规划系统及方法 |
CN105529748A (zh) * | 2016-01-11 | 2016-04-27 | 中国南方电网有限责任公司 | 一种适用于电力系统动态仿真的自动发电控制系统及方法 |
CN105760998A (zh) * | 2016-02-16 | 2016-07-13 | 国电南瑞科技股份有限公司 | 一种引入竞争机制的调频资源调用优化决策系统及方法 |
CN107579515A (zh) * | 2017-07-31 | 2018-01-12 | 东南大学 | 一种改善电厂调频性能的容量配置方法 |
CN107681695A (zh) * | 2017-10-30 | 2018-02-09 | 华泰慧能(北京)能源技术有限公司 | 一种储能辅助火电机组调频的容量配置方法 |
CN108711891A (zh) * | 2018-04-12 | 2018-10-26 | 广东电网有限责任公司 | 一种考虑机组调频性能的调频辅助服务优化系统及方法 |
JP2019115249A (ja) * | 2017-12-25 | 2019-07-11 | 赫普科技発展(北京)有限公司 | 電力網に周波数変調を行う方法、装置及びそのシステム |
CN110380440A (zh) * | 2019-08-19 | 2019-10-25 | 阳光电源股份有限公司 | 一种火储系统及其储能系统agc调频方法和装置 |
CN110994648A (zh) * | 2019-12-27 | 2020-04-10 | 广东电网有限责任公司 | 一种应用于调频辅助服务中的发电单元的调用方法 |
CN111082455A (zh) * | 2019-12-13 | 2020-04-28 | 中国电力科学研究院有限公司 | 一种风电场并网性能监测方法及系统 |
-
2020
- 2020-06-02 CN CN202010488689.2A patent/CN113765157B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105207242A (zh) * | 2015-09-17 | 2015-12-30 | 山东大学 | 储能装置参与机组调频的优化控制与容量规划系统及方法 |
CN105529748A (zh) * | 2016-01-11 | 2016-04-27 | 中国南方电网有限责任公司 | 一种适用于电力系统动态仿真的自动发电控制系统及方法 |
CN105760998A (zh) * | 2016-02-16 | 2016-07-13 | 国电南瑞科技股份有限公司 | 一种引入竞争机制的调频资源调用优化决策系统及方法 |
CN107579515A (zh) * | 2017-07-31 | 2018-01-12 | 东南大学 | 一种改善电厂调频性能的容量配置方法 |
CN107681695A (zh) * | 2017-10-30 | 2018-02-09 | 华泰慧能(北京)能源技术有限公司 | 一种储能辅助火电机组调频的容量配置方法 |
JP2019115249A (ja) * | 2017-12-25 | 2019-07-11 | 赫普科技発展(北京)有限公司 | 電力網に周波数変調を行う方法、装置及びそのシステム |
CN108711891A (zh) * | 2018-04-12 | 2018-10-26 | 广东电网有限责任公司 | 一种考虑机组调频性能的调频辅助服务优化系统及方法 |
CN110380440A (zh) * | 2019-08-19 | 2019-10-25 | 阳光电源股份有限公司 | 一种火储系统及其储能系统agc调频方法和装置 |
CN111082455A (zh) * | 2019-12-13 | 2020-04-28 | 中国电力科学研究院有限公司 | 一种风电场并网性能监测方法及系统 |
CN110994648A (zh) * | 2019-12-27 | 2020-04-10 | 广东电网有限责任公司 | 一种应用于调频辅助服务中的发电单元的调用方法 |
Non-Patent Citations (2)
Title |
---|
曹伟 等: "新电改条件下电力市场需求侧响应交易机制分析", 电网与清洁能源, vol. 35, no. 10, 31 October 2019 (2019-10-31), pages 50 - 55 * |
袁桂丽,苏伟芳: "计及电动汽车不确定性的虚拟电厂参与AGC 调频服务研究", 电网技术, 29 April 2020 (2020-04-29), pages 1 - 13 * |
Also Published As
Publication number | Publication date |
---|---|
CN113765157B (zh) | 2024-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112182720B (zh) | 一种基于建筑能源管理应用场景的建筑能耗模型评价方法 | |
CN109873455B (zh) | 一种储能辅助火电机组agc调频方法及系统 | |
CN112332404B (zh) | 采暖服务智能管理系统及方法 | |
CN112904266A (zh) | 一种电能表寿命预测方法及装置 | |
CN109508895A (zh) | 烟草制丝设备的控制性能评估装置、方法以及存储介质 | |
CN115860644A (zh) | 一种基于互联网的云仓储智慧拣选系统 | |
CN110649638B (zh) | 一种弥补风电功率预测误差的储能系统的优化方法 | |
CN117406844B (zh) | 一种基于神经网络的显卡风扇控制方法及相关装置 | |
CN118117581A (zh) | 针对特定时段的区域发电量预测方法及装置 | |
CN113765157B (zh) | 适用于调频辅助服务市场的发电单元调频性能测算方法和系统 | |
CN117674142A (zh) | 一种基于分时电量核算的电力调度方法及调度装置 | |
US20240097444A1 (en) | Hybrid system and method for distributed virtual power plants integrated intelligent net zero | |
CN114399166A (zh) | 一种基于电价引导及意愿的源荷储效益评估系统及方法 | |
CN116432510A (zh) | 一种基于集成学习-天牛须优化算法的柔性负荷预测方法 | |
CN115898787A (zh) | 一种风电机组静态偏航误差动态识别方法及装置 | |
CN111614105B (zh) | Agc控制时段不合格原因分析方法、系统和存储介质 | |
CN111817339B (zh) | 一种量化电化学储能提升agc调节性能效果的方法 | |
CN105116723B (zh) | 一种用于需量响应的能耗削减评估算法 | |
CN113964820B (zh) | 一种可控负荷资源参与电力系统调频的参与因子计算方法 | |
CN117949886B (zh) | 互感器校验仪智能调控方法、系统、电子设备及存储介质 | |
CN111900730B (zh) | 一种输电通道热稳定限额区间在线计算方法和装置 | |
CN113177692B (zh) | 一种国家指令性电量的年计划完成风险评估方法 | |
CN112712105A (zh) | 一种基于场景划分法的风电鲁棒调度区间优化方法 | |
CN117767261A (zh) | 一种多元资源调节性能及效果评价方法、系统、设备、存储介质 | |
CN117674234A (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 |