基于多源数据的电网一次调频能力平行计算研究

电网频率是电力系统运行中最重要的关键参数之一, 系统频率的异常波动不仅影响电网的电能质量, 制约电力设备的正常使用, 严重的频率失衡甚至将导致并网设备脱网及系统崩溃, 对电力系统的稳定运行带来严重危害[1-2]. 因此, 系统频率的稳定控制是电力系统稳定不可或缺的重要组成. 电力系统频率的异常波动是由并网机组发电功率和用电负荷功率不平衡引起的, 因而, 电力系统能够根据系统频率的变化主动快速进行发用电控制以维持频率稳定就显得尤为重要[3]. 电力系统的一次调频就是当系统频率发生异常波动时, 自动迅速地作出响应, 改变并网机组的实时出力, 将系统频率稳定在一个合理安全的水平线上. 这对提高电能质量及电网频率控制水平, 迅速平息电网频率波动起到了重要作用[4-5].

电力系统是一个非线性时变强耦合的复杂系统, 其一次调频特性与每台并网运行的发电机调速器特性息息相关, 是每台并网运行机组调速器特性在当前电网工况下的综合特性[6]. 对于单台发电机组而言, 其一次调频特性与调速器自身的特性显著相关[7]. 而当机组并入电力系统时, 电力系统的一次调频特性不仅与机组自身的调速器特性有关, 也和所有并列运行的其他发电机组特性相关, 甚至还与机组所处的实时工况和机组间的相互作用有关[8]. 另一方面, 电网中各机组的发电类型、容量、设备结构(例如, 凝汽式机组的负荷响应速度明显高于背压式机组), 以及机组运行工况(例如, 滑压或定压、不同季节乃至一天内不同时间段的负荷率高低)等各类因素共同作用影响着电网的一次调频性能[9]. 所以, 在不同的电网运行工况下, 非常有必要对电网的一次调频性能进行在线估计, 使电网面对不断变化的负荷调节需求均能够具备足够的一次调频性能.

目前, 国内外对一次调频能力的研究偏向单台机组侧较多, 对电网一次调频性能的研究和计算比较分散, 整体性考虑电网一次调性能各个因素的综合性研究较少, 也未见实际的工程化应用. 此外, 由于很难、甚至无法对电网一次调频特性进行精确建模, 以至于无法有效估计电网真实的一次调频性能.

本文基于平行控制理论框架下由人工社会 (Artificial societies, A)、计算实验(Computational experiments, C)和平行执行 (Parallel execution, P) 构成的ACP理论体系, 构建了多源数据的电网一次调频性能平行计算平台. 本文提出的计算方法能够有效兼顾机组类型的静态特性和运行工况的动态特性, 并以平行执行方式完成人工估算系统与实际电力系统的滚动优化, 实现了电网一次调频性能的在线全面估计, 从而有效弥补其无法进行精确建模的不足, 有助于提高电网主动应对负荷突变的能力, 也是对电网一次调频性能管控的一种全新的可行性方式.

1.   平台整体框架设计

所谓平行系统, 是指由某一个自然的现实系统和对应的一个或多个虚拟或理想的人工系统所组成的共同系统[10]. 平行系统主要包括实际系统和人工系统, 其基本框架如图1所示[10].

图 1  平行系统运行的基本框架

Fig. 1  Basic framework of parallel systems

设计平行系统的最主要目的是通过人工系统与实际系统的相互连接, 对二者之间的行为进行比对和分析, 完成对各自未来状况的“借鉴”和“预估”, 相应地调节各自的管理与控制方式, 达到实施有效解决方案以及学习和培训的目的[10-11].

人工系统试图尽可能地模拟实际系统, 对其行为进行预估, 从而为寻找对实际系统有效的解决方案或对当前方案进行改进提供依据[12]. 进一步, 通过观察实际系统与人工系统评估的状态之间的不同, 生成误差反馈信号, 对人工系统的评估方式或参数进行修正, 并重复进行新一轮的优化和评估[12].

目前, 平行系统理论已经在自动化码头[13]、智能交通[14]等领域中开展了应用研究, 但尚未用于现代电力系统. 现代电力系统是一个巨大的参数海量、广域分布、模型复杂的非线性复杂系统. 现阶段, 绝大多数电力仿真系统是基于物理模型进行的, 并没有深入地考虑人的行为、自然因素以及社会环境的影响, 这些因素限制了对复杂现代电力系统运行影响的评估, 尤其是无法定性、定量地对电力系统应对突发事件的能力问题进行分析评估[15].

为了克服采用传统仿真方法导致的上述局限, 本文将平行系统理论应用到电力系统的实际工程中, 改变传统的被动、离线的电力系统仿真方法, 建立人工电力平行计算平台, 进而实现管理与控制、实验与评估、学习与优化及故障诊断等功能. 以现代电网的海量数据为驱动, 电力平行计算平台为载体, 该平台可对正在运行的电力系统进行滚动式在线改进与不间断优化, 实现电网一次调频性能估算达到最优解, 为电力系统和复杂电网的能力优化与管理开辟一条新的途径[16].

基于上述平行系统理念, 本文设计了电网调频能力计算平台, 其整体框架如图2所示, 主要包括实际系统和人工平行系统两部分. 实际电力系统的频率控制本质上是一个以电力调控中心为中枢, 实现电源与负荷的电力供需平衡的闭环控制系统. 电力调控中心不仅需要满足区域内机组发电及负荷消纳的控制要求, 也肩负着电力系统实际运行产生的海量数据信息采集的重担. 人工平行系统中, 基于D5000电网调度系统的一次调频在线测试系统不仅能够在正常运行工况下测试控制区内的机组调频性能, 保证机组一次调频性能的测试结果更具真实性, 同时将历次测试数据存储、分析及评价, 实现机组一次调频的静态参数、动态性能的数据积累和量化管理. 将电力调控中心的海量数据和一次调频测试系统的特性数据应用于电网一次调频能力在线估计模型中, 通过实际电力系统和人工平行系统的彼此作用、相互修正, 既可以实现对实际电力系统的管理与控制,也能够完成对电网能力、行为与决策的评估与实验.

图 2  基于平行系统的计算平台整体框架

Fig. 2  Overall framework of computing platform based on parallel systems

2.   发电机组一次调频数据采集体系

当前, 电力系统已经发展成一个具有多源信息融合的综合性电力控制系统. 根据现有的SCADA(Supervisory control and data acquisition)、WAMS(Wide area measurement system)、OMS (Operations management system)等系统, 建立机组一次调频数据采集体系, 如图3所示. SCADA的一次调频在线测试系统能够在一次调频在线测试期间自动记录测试机组的有功PMU (Phasor measurement unit)数据, 整编录入历史库存储[17]. WAMS的机组性能考核系统能够不间断采集机组实际运行的一次调频动作、响应等广域量测数据, 实现机组转速不等率、积分电量、响应时间等一次调频性能指标的在线计算[18]. 基于OMS的网源协调信息管理系统, 通过设置统一的一次调频试验台账, 维护和存储机组的历次原始试验数据和参数, 统一存储于数据库中.

图 3  发电机组一次调频数据采集体系

Fig. 3  Primary frequency control data acquisition system of generator sets

2.1   机组一次调频在线测试系统

如图4所示, 在D5000系统中部署试验触发、结果展示、参数维护及统计分析功能. 试验触发功能可实现选定机组进行试验, 试验开始时将测试机组投入一次调频扰动测试模式, 调控中心向电厂侧下发调频测试开始信号和频率偏差信号, 机组根据接收的信号进行一次调频动作响应, 具体控制信息流程见图5所示.

图 4  机组一次调频在线测试与评价功能结构图

Fig. 4  Functional structure diagram of unit primary frequency control online test and evaluation

图 5  机组一次调频在线测试系统图

Fig. 5  System diagram of unit primary frequency control online test

在一次调频在线测试期间, 系统自动记录测试机组的有功PMU数据, 整编录入D5000历史库, 永久保存, 供后期调频指标分析与数据展示使用[17]. 在结果展示功能中, 将频率扰动信号及电厂侧PMU装置上传的有功出力信号进行展示, 并根据有功出力曲线自动计算一次调频动作性能指标. 在浙江某厂#2机组进行一次调频在线测试系统联调试验时, 调控中心主站侧的操作及结果展示如图6所示, 发电厂机组一次调频动作响应曲线如图7所示.

图 6  机组一次调频在线测试主站侧界面图

Fig. 6  Interface diagram of main station side of unit primary frequency control online test

图 7  浙江某厂#2机组一次调频在线测试动作曲线(75%负荷点 + 11 r/min转速偏差)

Fig. 7  Action curves of primary frequency control online test of unit 2 of a power plant in Zhejiang (+11 r/min speed deviation at 75% load point)

通过PMU实时采集高频度的功率、相位、相角等机组数据, 并利用PMU的GPS对时优势, 计算机组在一次调频在线测试中不同时刻、不同工况下的性能指标, 反向推算出机组实际的调频特性参数, 以此形成海量的机组一次调频特征数据, 实现机组一次调频静态参数、动态性能等特征数据的分析及管理[19].

2.2   网源协调信息管理系统

基于OMS的网源协调信息管理系统(如图8所示)能够实现设备台账、电气参数、涉网试验报告等技术资料的及时更新及有效管理, 打通厂网之间的涉网信息共享交互渠道, 更加科学、系统、高效地开展网源协调管理工作. 以各机组台账信息为基础, 实现对网内各并网机组数据信息的实时统计汇总; 以并网机组各类信息关键字为基础, 实现对特定类型、特定机组相关信息的快速查询和准确统计, 该系统所部署的功能模块如图9所示. 同时, 建立网源协调数据中心, 这是网源协调信息管理平台的核心内容, 可实现参数模型库、文档资料库、督办信息库及信息变更库的数据获取与整合、数据质量校核、合格率督办和数据查询等服务.

图 8  基于OMS的网源协调信息管理系统主界面

Fig. 8  Main interface of grid power coordination information management system based on OMS

图 9  网源协调信息管理系统的功能模块

Fig. 9  Functional modules of grid power coordination information management system

在网源协调信息管理系统中, 通过设置统一的一次调频试验台账模板, 采集机组的历次试验数据, 以维护和存储机组的原始试验数据. 图10为浙江某百万千瓦机组A修后一次调频试验的数据台账.

图 10  浙江某厂#1机组A修后一次调频试验台账

Fig. 10  Ledger of primary frequency control test after class a maintenance of unit 1 of a power plant in Zhejiang

根据设置统一的数据台账, 可以追溯机组每一次试验工况, 分析不同工况下机组一次调频试验时的真实性能.

2.3   发电机组一次调频性能考核系统

发电机组一次调频性能考核系统是基于WAMS部署于D5000上的一次调频数据量化及评价系统[18]. 该系统通过PMU将现场采集的机组出力数据送至电网调控中心EMS (Energy management system)系统中, 并由EMS采集电网中的频率值及时间, 打包形成一个数据包存储于该系统, 以便进行量化处理.

机组一次调频贡献电量是机组一次调频多项指标的综合, 它反映机组对电网频率的实际贡献大小. 这一参数的引入使量化机组一次调频功能的动作力度和快速性成为可能.

在系统设计中, 机组对象被抽象封装成一个类, 每个机组都是该类的实例化. 类中数据为机组的实时信息, 包括机组基本参数(代号、装机容量、积分间隔、频率超前时间等)和积分电量以及其他计算辅助变量; 类中对数据的操作包含积分电量的计算和对数据库的读写操作等. 其中最为核心的函数为贡献电量积分函数, 用来计算和累积贡献电量, 以反映机组在实际运行时的一次调频动作性能.

一次调频性能考核系统通过采集各机组一次调频投退信号、有功功率和频率信息等广域量测信息, 实现机组转速不等率、积分电量、响应时间等一次调频性能指标的在线计算, 使系统能够考核机组是否贡献了满足要求的一次调频服务.

该系统依托于调度D5000系统, 能够不间断实时获取所有并网机组的一次调频动作性能. 据系统统计, 每个月一台机组的一次调频动作次数可达到近5000次(如表1所示), 大量原始数据的累积奠定了电网一次调频性能计算的数据基础.

3.   电网一次调频能力数据建模

目前, 国内外各类指标都是针对机组调频能力, 比如控制性能标准(Control performance standard, CPS)、投运率、贡献电量、调频效果等[20]. 现有的评价指标本质为事后评价, 不能预测电网总的调频裕度. 本文构建的数据平行驱动体系(如图11所示), 依靠电网运行的大数据, 通过不断演算、迭代优化, 形成调频能力的概率区间来衡量电网调频性能, 解决了电网一次调频能力的在线估计问题, 从而使电网的运行状态具有更高的可控性及透明度, 提高大电网应对负荷突变的能力.

图 11  电网一次调频能力估算技术路线图

Fig. 11  Technical roadmap of grid primary frequency control capacity estimation

3.1   同步电网一次调频能力建模

定义 1. PFCA (Primary frequency control ability)表示仅在一次调频作用情况下(即没有二次调频作用时), 某时间段内负荷变化与电网频率变化的比值[8], 可表示为

同步发电机并列运行的电网调频能力模型如图12所示. 在图12中, Ri$R_i$为二次调频的给定值; δi${\delta }_i$为电网中第i台机组的不等率; αi${\alpha }_i$为电网中第i台机组装机容量/电力系统装机容量; Gi(s)$G_i(s)$为第i台汽轮机的传递函数; PNL(s)$P_{NL}(s)$为负载的标幺值; φs(s)${\phi }_s(s)$为电网转速变化的标幺值. TαΣ=∑Mi=1αiTαi$T_{\alpha \mathrm{\Sigma }}={\sum }_{i=1}^M{\alpha }_i{T}_{\alpha i}$为电网的惯性时间常数; Tαi$T_{\alpha i}$为第i台汽轮机的转子时间常数; βΣ${\beta }_{\mathrm{\Sigma }}$为电网负载频率特性系数.

由图12可推出电网的周波变化为[8]

图 12  M台同步发电机并列运行的电网调频能力模型图

Fig. 12  Model diagram of power grid frequency regulation capacity of M synchronous generators operating in parallel

其中,

各变量前加Δ$\mathrm{\Delta }$作为静态标记, 静态下电网周波变化可表示为

由于汽轮机传递函数可以表示为

故

由此得到, 在静态情况下电网的一次调频能力可表示为

RG=∑Mi=1(αi/δi)$R_G={\sum }_{i=1}^M({\alpha }_i/{\delta }_i)$定义为电网的等值调差率, 等效为所有挂网同步发电机组一次调频出力总和的最大值. RL=β∑$R_L={\beta }_{\sum }$为电网负载频率特性系数, 等效为电网负荷对频率变化的贡献能力.

对于某台额定功率Pei$P_{ei}$, 不等率δi${\delta }_i$的机组而言, 其在Δf$\mathrm{\Delta }f$频率扰动下的一次调频能力的实际贡献为

显然, 在Δf$\mathrm{\Delta }f$频率扰动下的同步电网中所有挂网运行机组的一次调频实际出力总和为

3.2   机组一次调频能力数据模型

从式(11)中可以发现电网一次调频能力主要取决于各挂网运行机组的转速不等率.

机组在不同功频下的转速不等率不尽相同, 将每台机组的转速不等率设定为常数势必难以准确描述机组的一次调频能力. 调度D5000系统具有各类丰富的一次调频数据(如图13所示), 其中, 机组性能考核系统实时采集的机组一次调频动作数据量最大, 每月达几千次. 而网源协调信息管理系统中存储的一次调频参数及试验数据, 则表征机组一次调频性能的理论出力和实际能力上下限值. 一次调频在线测试系统可以探测机组在不同工况下的大频差真实能力, 完善整个机组的一次调频数据集. 本文将采用极大似然估计及数值拟合等算法对调度D5000系统中的多源一次调频数据进行有效整合, 以获取机组真实的一次调频性能功频特性图谱.

图 13  发电机组一次调频动作历史数据库

Fig. 13  Historical database of primary frequency control actions of generator sets

可以从D5000系统中获取不同时刻机组一次调频能力的数据, 从而构造出机组转速不等率的样本空间. 我们用随机变量δi${\delta }_i$表示单台机组在第i个时间段内的转速不等率. 利用D5000系统中不同时间的转速不等率数据, 可以构造用于描述转速不等率的随机变量δ1,⋯,δn${\delta }_1,\cdots ,{\delta }_n$. 由于机组的短期特性通常不会发生变化, 可以认为δ1,⋯,δn${\delta }_1,\cdots ,{\delta }_n$独立同分布. 按照中心极限定理, 随机变量δ1,⋯,δn${\delta }_1,\cdots ,{\delta }_n$独立同分布, 当n很大时, (∑ni=1δi−nμ)/(n−−√σ)$({\sum }_{i=1}^n{\delta }_i-n\mu )/(\sqrt{n}\sigma )$近似服从标准正态分布N(0,1)$\mathrm{N}(0,1)$. 设X=∑ni=1δi$X={\sum }_{i=1}^n{\delta }_i$, 则X服从正态分布N(nμ,nσ2)$\mathrm{N}(n\mu ,n{\sigma }^2)$.

本文通过极大似然估计来估计μ$\mu$和σ2${\sigma }^2$, 似然函数可以表示为

那么, 似然函数的对数为

则似然方程组为

求解(14a), 得

代入(14b), 得

那么, μ$\mu$和σ2${\sigma }^2$的最大似然估计为

本文用大写字母表示所有涉及的样本, 因为最大似然估计μ^$\hat{\mu }$和σ^2${\hat{\sigma }}^2$均为统计量, 脱离具体的试验或观测, 它们都是随机的. 依据3套系统得到3个维度的真实观测数据, 采用上述的估计算法, 即可得到机组功频坐标下的一次调频性能估计值, 记为K$K$, 显然K=μ^$K=\hat{\mu }$, 从而得到功频坐标下的机组一次调频性能网格如表2所示, 对浙江某660万千瓦级超临界机组进行原始数据采集和数值迭代模拟, 可得到该机组的一次调频性能三维分布, 如图14所示, 其中, 横轴f$f$和P$P$分别为电网频率和机组实际出力值, 纵轴K$K$表示为机组一次调频性能值.

图 14  发电机组一次调频性能功频图谱

Fig. 14  Power frequency spectrum of primary frequency control performance of generator sets

4.   电网一次调频性能数据分析及估计

4.1   浙江电网一次调频性能分析

2018年10月22日10 : 36时, 浙江省全社会口径用电负荷39753.02 MW, 受电17849.91 MW (占全省用电负荷的52.9%), 发电21903.11 MW (占全省用电负荷的47.1%). 其中, 浙江省调统调机组共并网45台(不含新能源), 其中纯凝汽燃煤机组39台, 背压供热燃煤机组2台, 燃气机组4台, 并网机组平均负荷率约为58.25%, 总发电负荷16221.59 MW, 占当时全省用电负荷的40.8%.

据此, 按照式(11)进行计算, 得到浙江2018年10月22日10 : 36时省调统调并网机组的实际调频能力为12.2327, 根据当时的用电负荷进行折算后约为6 762.23 MW/Hz, 略高于浙江电网ACE系统中自然频率特性系数5 979 MW/Hz[21]. 推此及彼, 采用本文的方法可以得到浙江统调机组每天任意时刻的实际调频能力值, 如图15所示.

图 15  2018年10月22日浙江统调机组实际调频能力时刻图

Fig. 15  Time map of actual primary frequency control performance of Zhejiang dispatching unit on October 22, 2018

显然, 由图15可以明显得知, 电网的实际调频性能与并网机组台数、并网机组类型、机组平均负荷率等因素息息相关, 是一项电网拓扑结构的综合表征.

4.2   大频差时浙江电网一次调频容量分析

电网频率是描述电力系统发电与用电平衡关系的一个重要变量. 电网中用电负荷的随机特性使电网频率也成为随机变量. 电力负荷的随机变化部分是相互独立、不相关的[22]. 由中心极限定理可知, 连续的独立随机变量叠加后是正态分布的, 因此, 电力系统中负荷的随机性分量也是正态分布的. 在通常情况下, 电网频率的概率分布也就可近似认为是正态分布[23], 因此, 假设发生直流闭锁时的浙江电网系统最低频率近似正态分布, 并结合近三年的历次直流闭锁数据分析, 可得表3数据.

根据表3, 可知电网频率的最大变化差值为

最大概率功率缺额为

按浙江历次直流闭锁数据, 可知在不少于63台同步发电机组, 总装机不低于34 658.58 MW的条件下, 浙江统调机组一次调频能力为9.5856 (95%以上), 折算后为6 644.46 MW/Hz, 在最大概率频差0.159 Hz下, 浙江电网可提供一次调频最大能力至少约为1 056.47 MW (3.1% MCR (Maximum continuous rating)), 能够满足功率缺额950.66 MW的要求, 也就是说, 浙江统调机组一次调频能力可以保证电网在直流闭锁等大频差工况下系统频率不低于49.841 Hz.

5.   结束语

机组的转速不等率与电网频率往往呈非线性关系, 采用固定的转速不等率来计算机组的一次调频能力误差较大. 本文基于数据平行驱动体系, 依靠电网运行的大数据, 建立了更为准确的电网一次调频模型, 对提升电网频率控制水平具有非常重要的意义.