CN109245090B - 一种大功率缺失下频率最低点预测的解析模型的建模方法 - Google Patents
一种大功率缺失下频率最低点预测的解析模型的建模方法 Download PDFInfo
- Publication number
- CN109245090B CN109245090B CN201811036058.6A CN201811036058A CN109245090B CN 109245090 B CN109245090 B CN 109245090B CN 201811036058 A CN201811036058 A CN 201811036058A CN 109245090 B CN109245090 B CN 109245090B
- Authority
- CN
- China
- Prior art keywords
- frequency
- lowest point
- fitting
- unit
- response
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000007812 deficiency Effects 0.000 title abstract description 6
- 230000004044 response Effects 0.000 claims abstract description 80
- 230000009471 action Effects 0.000 claims abstract description 4
- 230000006735 deficit Effects 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 230000014509 gene expression Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000007142 ring opening reaction Methods 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000007123 defense Effects 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012360 testing method 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/382—
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种大功率缺失下频率最低点预测的解析模型的建模方法,步骤为:S1:动态频率响应模型解环及模拟频差输入;S2:在频域内计算模拟频差输入作用下各机组的频率响应出力;S3:求解调速器的频率响应功率;S4:对各机组调速器阶跃响应曲线进行离散积分运算得到其斜坡响应曲线;S5:利用多项式对斜坡响应曲线进行拟合;S6:利用最小二乘法确定多项式拟合系数;S7:计算机组在系统频率到达最低点时刻的频率响应出力;S8:计算系统在系统频率到达最低点时刻的频率响应出力;S9:建立大功率缺失下频率最低点预测的解析模型。本发明基于各调速器阶跃响应曲线,可建立大功率缺失下频率最低点预测的解析模型,并准确快速计算大功率缺失下的最大频率偏移和频率最低点时间。
Description
技术领域
本发明属于电力系统分析领域,涉及一种大功率缺失下频率最低点预测的解析模型的建模方法。
背景技术
随着高比例可再生能源的接入,风力、光伏等新能源机组将逐步成为电力系统的主力电源。相比传统发电机组,新能源机组的惯性响应和频率响应能力先天不足,其大规模并网将使得系统频率稳定问题日益凸显。另一方面,随着特高压交直流输电工程的发展,由于直流闭锁故障造成的大功率缺失,易引发直流馈入受端电网的低频事故。上述两方面因素的作用,将导致频率控制“攻防”两端出现此消彼长的态势,从而使电力系统的频率稳定面临严峻的挑战。
频率响应作为电力系统频率控制的第一道防线,特别是大功率缺失下阻止频率大幅下降的关键措施,一直以来备受重视。在频率响应过程中,频率最低点位置既是频率监视的关键参量,亦是频率控制的决策依据,如快速频率响应资源调用和低频减载控制等。因此,准确估计最大频率偏移和频率最低点时间尤为重要。
发明内容
根据现有技术存在的问题,本发明公开了一种大功率缺失下频率最低点预测的解析模型的建模方法。
本发明的技术方案为:
一种大功率缺失下频率最低点预测的解析模型的建模方法,包括以下步骤:
S1:采用电力系统仿真软件,如PSS/E,BPA或PowerWorld对大规模电力系统进行动态频率仿真需建立高阶网络结构模型和元件模型,涉及诸多模型参数,其仿真过程一般耗时较长,通常难以满足频率最低点在线预测的需要。为此,建立了大功率缺失下频率最低点预测的解析模型,其核心思想是通过频率反馈解环和模拟二次频差输入实现系统频率和频率响应功率的解耦计算。
所述的模拟二次频差输入如公式(1)所示:
式中:Pdeficit为功率缺额;Hsys为系统等效惯性时间常数;tnadir为系统频率到达最低点的时间;t为时间变量。
S2:在大功率缺失下的频率响应过程中,各机组依据频率偏差进行功率调整。因此,在模拟频差输入ΔfFNL(t)的作用下,各机组的频率响应出力为:
式中:Gunit,i(s)为调速器i的传递函数;Pdeficit为功率缺额;Hsys为系统等效惯性时间常数;tnadir为系统频率到达最低点的时间;s为Laplace算子;i为机组编号。
S5:在平均系统频率模型中,由于保留了机组调速器动态,使得该模型阶数随机组数量增长而提高。针对该问题,采用统一数学形式描述不同类型调速器的响应特性以进一步降低系统模型阶数,即利用n次多项式对斜坡响应曲线在t∈[0,tfit]内进行拟合。
所述的n次多项式如公式(4)和公式(5)所示:
式中:m为时域拟合区间内的样本数量;k为样本数据编号;i为机组编号;tk为采样时间;tfit为时域拟合时长。
S7:在系统发生大功率缺失故障后,各机组频率响应出力不断增加,在系统频率到达最低点时,其频率响应出力为:
式中:Pdeficit为功率缺额;Hsys为系统等效惯性时间常数;tnadir为系统频率到达最低点的时间;l为多项式拟合阶数编号;i为机组编号。
S8:系统统频率响应是所辖区域内各机组频率响应的综合。在发生大功率缺失故障后,所有机组共同参与频率响应过程,其频率响应出力是前述各机组频率响应出力的叠加。因此,包含NG台机组的系统在系统频率到达最低点时刻的频率响应出力为:
其中:
式中:Ki为标幺值转换系数,即从机组容量标幺值转化为系统容量标幺值。NG为系统内频率响应机组的数量;i为机组编号;tnadir为系统频率到达最低点的时间;n为拟合阶数;l为拟合阶数编号;为机组调速器i斜坡响应的多项式拟合系数向量;为系统等效调速器斜坡响应的多项式拟合系数向量。
S9:在频率到达最低点时刻,系统频率响应出力与功率缺额大小相等。基于此,建立大功率缺失下频率最低点预测的解析模型,并利用数值方法求解最大频率偏移和频率最低点时间。
所述大功率缺失下频率最低点预测的解析模型如公式(12)所述:
与现有技术相比,本发明的有益效果为:
(1)数据需求少。大功率缺失下频率最低点预测的解析模型无须已知各调速器结构及参数,仅对调门扰动试验曲线进行多项式拟合便可建立系统的频率响应模型。
(2)模型阶数低。大功率缺失下频率最低点预测的解析模型不随机组数量增多而进一步升高,有效提高了模型的求解计算速度。
(3)计算精度高。多项式拟合阶数和拟合时长对模型求解结果影响不大,采用低阶模型便可获取较高的计算精度。
附图说明
图1为大功率缺失下频率最低点预测解析建模的流程图;
图2为本发明动态频率响应模型解环及模拟频差输入示意图。
具体实施方式
为使本发明的技术方案和优点更加清楚,下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚完整的描述,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
如图1所示的一种大功率缺失下频率最低点预测的解析建模流程,包括以下步骤:
S1:采用电力系统仿真软件,如PSS/E,BPA或PowerWorld对大规模电力系统进行动态频率仿真需建立高阶网络结构模型和元件模型,涉及诸多模型参数,其仿真过程一般耗时较长,通常难以满足频率最低点在线预测的需要。为此,建立了大功率缺失下频率最低点预测的解析模型,其核心思想是通过频率反馈解环和模拟二次频差输入实现系统频率和频率响应功率的解耦计算。
频率反馈解环及模拟二次频差输入过程如图2所示。
所述的模拟二次频差输入如公式(1)所示:
式中:Pdeficit为功率缺额;Hsys为系统等效惯性时间常数;tnadir为系统频率到达最低点的时间;t为时间变量。
S2:在大功率缺失下的频率响应过程中,各机组依据频率偏差进行功率调整。因此,在模拟频差输入ΔfFNL(t)的作用下,各机组的频率响应出力为:
式中:Gunit,i(s)为调速器i的传递函数;Pdeficit为功率缺额;Hsys为系统等效惯性时间常数;tnadir为系统频率到达最低点的时间;s为Laplace算子;i为机组编号。
S5:在平均系统频率模型中,由于保留了机组调速器动态,使得该模型阶数随机组数量增长而提高。针对该问题,采用统一数学形式描述不同类型调速器的响应特性以进一步降低系统模型阶数,即利用2次多项式对斜坡响应曲线在0~5s内进行拟合。
所述的2次多项式如公式(4)和公式(5)所示:
式中:m为时域拟合区间内的样本数量;k为样本数据编号;i为机组编号;tk为采样时间。
利用最小二乘法计算所得G1~G10调速器5s的二阶多项式拟合结果如表1所示。
表1各机组二阶多项式拟合结果
S7:在系统发生大功率缺失故障后,各机组频率响应出力不断增加,在系统频率到达最低点时,其频率响应出力为:
式中:Pdeficit为功率缺额;Hsys为系统等效惯性时间常数;tnadir为系统频率到达最低点的时间;l为多项式拟合阶数编号;i为机组编号。
S8:系统统频率响应是所辖区域内各机组频率响应的综合。在G6发生机组跳闸故障后,其余机组共同参与频率响应过程,其频率响应出力是前述各机组频率响应出力的叠加。因此,包含10台机组的系统在系统频率到达最低点时刻的频率响应出力为:
其中:
其中,标幺值转换系数K=[0.5263 0.0526 0.0526 0.0526 0.0526 0.05260.0526 0.0526 0.0526 0.0526];P deficit=3.42%;Hsys=3.9363。
S9:在频率到达最低点时刻,系统频率响应出力与功率缺额大小相等。基于此,建立大功率缺失下频率最低点预测的解析模型,并利用数值方法求解最大频率偏移和频率最低点时间。
所述大功率缺失下频率最低点预测的解析模型如公式(12)所述:
以上所述实施例仅表达了本发明的实施方式,但并不能因此而理解为对本发明专利的范围的限制,应当指出,对于本领域的技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些均属于本发明的保护范围。
Claims (1)
1.一种大功率缺失下频率最低点预测的解析模型的建模方法,其特征在于以下步骤:
S1:首先,通过频率反馈解环和模拟二次频差输入实现系统频率和频率响应功率的解耦计算;
所述的模拟二次频差输入如公式(1)所示:
式中:Pdeficit为功率缺额;Hsys为系统等效惯性时间常数;tnadir为系统频率到达最低点的时间;t为时间变量;
S2:在频域内计算模拟二次频差输入ΔfFNL(t)作用下,各机组的频率响应出力为:
式中:Gunit,i(s)为调速器i的传递函数;s为Laplace算子;i为机组编号;Pdeficit为功率缺额;Hsys为系统等效惯性时间常数;tnadir为系统频率到达最低点的时间;
S3:获取各机组频率响应出力的时域解析表达
S5:在平均系统频率模型中保留机组调速器动态,使得该模型阶数随机组数量增长而提高;采用统一数学形式描述不同类型调速器的响应特性以进一步降低系统模型阶数,即利用n次多项式对斜坡响应曲线在t∈[0,tfit]内进行拟合;
所述的n次多项式如公式(4)和公式(5)所示:
S6:求解各机组调速器斜坡响应曲线的多项式拟合参数
式中:m为时域拟合区间内的样本数量;k为样本数据编号;tk为采样时间;tfit为时域拟合时长;
S7:在系统发生大功率缺失故障后,各机组频率响应出力不断增加,在系统频率到达最低点时,其频率响应出力为:
式中:Pdeficit为功率缺额;Hsys为系统等效惯性时间常数;tnadir为系统频率到达最低点的时间;l为多项式拟合阶数编号;
S8:在发生大功率缺失故障后,所有机组共同参与频率响应过程,其频率响应出力是前述各机组频率响应出力的叠加;因此,包含NG台机组的系统在系统频率到达最低点时刻的频率响应出力为:
其中:
式中:Ki为标幺值转换系数,即从机组容量标幺值转化为系统容量标幺值;NG为系统内频率响应机组的数量;tnadir为系统频率到达最低点的时间;n为拟合阶数;l为拟合阶数编号;为机组调速器i斜坡响应的多项式拟合系数向量;为系统等效调速器斜坡响应的多项式拟合系数向量;
S9:建立大功率缺失下频率最低点预测的解析模型,并利用数值方法求解最大频率偏移和频率最低点时间;
所述大功率缺失下频率最低点预测的解析模型如公式(12)所述:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811036058.6A CN109245090B (zh) | 2018-09-06 | 2018-09-06 | 一种大功率缺失下频率最低点预测的解析模型的建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811036058.6A CN109245090B (zh) | 2018-09-06 | 2018-09-06 | 一种大功率缺失下频率最低点预测的解析模型的建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109245090A CN109245090A (zh) | 2019-01-18 |
CN109245090B true CN109245090B (zh) | 2021-04-20 |
Family
ID=65060997
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811036058.6A Expired - Fee Related CN109245090B (zh) | 2018-09-06 | 2018-09-06 | 一种大功率缺失下频率最低点预测的解析模型的建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109245090B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112865071B (zh) * | 2020-12-29 | 2023-04-14 | 四川大学 | 直流受端电网分布式光伏高渗透下频率偏移预测方法 |
CN114583746B (zh) * | 2022-02-15 | 2022-11-25 | 四川大学 | 一种预测新能源接入电网频率最低点的通用建模方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105138739A (zh) * | 2015-07-30 | 2015-12-09 | 南方电网科学研究院有限责任公司 | 一种计及死区效应的电力系统最低频率快速计算方法 |
CN106130017A (zh) * | 2016-08-01 | 2016-11-16 | 国电南瑞科技股份有限公司 | 一种电网整体一次调频响应时间的检测方法 |
CN106169115A (zh) * | 2016-06-23 | 2016-11-30 | 西南交通大学 | 一种基于解析‑数值结合的两区域电网最低频率估计算法 |
WO2017062914A1 (en) * | 2015-10-08 | 2017-04-13 | Johnson Controls Technology Company | Photovoltaic energy system with value function optimization |
CN106849092A (zh) * | 2017-03-27 | 2017-06-13 | 南方电网科学研究院有限责任公司 | 一种交直流电网扰动最大频率偏差的计算方法及装置 |
CN107026460A (zh) * | 2017-04-27 | 2017-08-08 | 大连理工大学 | 基于新能源大规模并网下频率偏差系数的确定方法 |
CN107546754A (zh) * | 2017-09-25 | 2018-01-05 | 国网河南省电力公司电力科学研究院 | 特高压大功率缺失下互联电网一次频率响应能力评价方法 |
CN107834610A (zh) * | 2017-11-29 | 2018-03-23 | 西南交通大学 | 一种考虑水轮机水锤效应的电网频率动态分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8694172B2 (en) * | 2011-07-12 | 2014-04-08 | General Electric Company | Systems and devices for controlling power generation |
-
2018
- 2018-09-06 CN CN201811036058.6A patent/CN109245090B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105138739A (zh) * | 2015-07-30 | 2015-12-09 | 南方电网科学研究院有限责任公司 | 一种计及死区效应的电力系统最低频率快速计算方法 |
WO2017062914A1 (en) * | 2015-10-08 | 2017-04-13 | Johnson Controls Technology Company | Photovoltaic energy system with value function optimization |
CN106169115A (zh) * | 2016-06-23 | 2016-11-30 | 西南交通大学 | 一种基于解析‑数值结合的两区域电网最低频率估计算法 |
CN106130017A (zh) * | 2016-08-01 | 2016-11-16 | 国电南瑞科技股份有限公司 | 一种电网整体一次调频响应时间的检测方法 |
CN106849092A (zh) * | 2017-03-27 | 2017-06-13 | 南方电网科学研究院有限责任公司 | 一种交直流电网扰动最大频率偏差的计算方法及装置 |
CN107026460A (zh) * | 2017-04-27 | 2017-08-08 | 大连理工大学 | 基于新能源大规模并网下频率偏差系数的确定方法 |
CN107546754A (zh) * | 2017-09-25 | 2018-01-05 | 国网河南省电力公司电力科学研究院 | 特高压大功率缺失下互联电网一次频率响应能力评价方法 |
CN107834610A (zh) * | 2017-11-29 | 2018-03-23 | 西南交通大学 | 一种考虑水轮机水锤效应的电网频率动态分析方法 |
Non-Patent Citations (4)
Title |
---|
"Frequency Deviations and Generation Scheduling in the Nordic System";Z. W. Li,et al;《2011 IEEE Trondheim PowerTech》;20111231;第1-6页 * |
"Maximum Frequency Deviation Calculation in Small Isolated Power Systems";Ignacio Egido, et al.;《IEEE TRANSACTIONS ON POWER SYSTEMS》;20091130;第24卷(第4期);第1731-1738页 * |
"Estimation of maximum wind power penetration level to maintain an adequatefrequency response in a power system";Nahid-Al-Masood, et al.;《8th International Conference on Electrical and Computer Engineering》;20141222;第587-590页 * |
"基于 v-SVR 的电力系统扰动后最低频率预测";薄其滨,等;《电力自动化设备》;20150731;第35卷(第7期);第83-88页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109245090A (zh) | 2019-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110120686B (zh) | 一种基于电力系统在线惯量估计的新能源承载力预警方法 | |
CN103715718A (zh) | 网源联合仿真及其多级调度闭环控制系统 | |
CN109638879A (zh) | 基于性能指标评估的一次调频动态补偿调整系统及方法 | |
CN104201700B (zh) | 计及风电不确定性波动的区域电网火电调频机组配置方法 | |
CN110429648B (zh) | 考虑风速随机波动的小干扰稳定裕度概率评估方法 | |
CN110492524A (zh) | 无通信海上风电场柔直并网的有功-频率控制方法及系统 | |
CN109245090B (zh) | 一种大功率缺失下频率最低点预测的解析模型的建模方法 | |
CN109149566B (zh) | 一种大功率缺失下频率最低点预测的仿真模型的建模方法 | |
CN114123344B (zh) | 基于自适应递推最小二乘的电力系统惯量评估方法及装置 | |
CN110061521A (zh) | 一种考虑频率累积作用的最大风电渗透率快速评估方法 | |
CN104167758B (zh) | 基于超调补偿的一次调频优化控制方法及系统 | |
CN111311021A (zh) | 一种风电场理论功率预测方法、装置、设备和存储介质 | |
CN110649596B (zh) | 一种考虑系统初始状态的频率全响应解析模型 | |
CN108879721A (zh) | 一种基于风电消纳的频率控制方法 | |
CN106021754B (zh) | 考虑vsc无功越限调整策略的混联电网概率潮流算法 | |
CN104037772A (zh) | 一种互联电网agc调频统一控制器及控制方法 | |
Zhang et al. | Frequency-constrained unit commitment for power systems with high renewable energy penetration | |
CN107979112B (zh) | 一种风机控制方法、系统、终端及可读存储介质 | |
CN116488267B (zh) | 一种基于建模的风电场无功容量极限仿真计算方法及装置 | |
US11262713B2 (en) | Method for calculating control parameters of heating supply power of heating network | |
CN110048428B (zh) | 基于概率守恒原理的电力系统概率潮流计算方法 | |
CN107346900B (zh) | 基于考核数据的一次调频参数自调整控制方法及系统 | |
CN111682552B (zh) | 数据驱动的无功电压控制方法、装置、设备及存储介质 | |
CN112290584B (zh) | 一种农村风光柴储微网频率控制方法、装置及终端设备 | |
CN114520521A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210420 Termination date: 20210906 |