一种风电集中接入后的互联电网分布式AGC控制方法 |
|||||||
申请号 | CN201410769057.8 | 申请日 | 2014-12-12 | 公开(公告)号 | CN104467015A | 公开(公告)日 | 2015-03-25 |
申请人 | 国家电网公司; 中国电力科学研究院; 国网浙江省电力公司; | 发明人 | 单茂华; 田伟; 鄢蜜昉; 严春华; 姚建国; | ||||
摘要 | 本 发明 提供一种 风 电集中接入后的互联 电网 分布式AGC控制方法,包括以下步骤:风电-非风电控制区协调决策器确定风电 波动 平衡责任系数矩阵;风电控制区 控制器 根据风电控制区承担自身控制区的风电波动平衡责任系数,确定风电控制区的AGC总调节量,并分配至各AGC机组执行;非风电控制区控制器根据非风电控制区承担的风电控制区的风电波动平衡责任系数,确定非风电控制区的AGC总调节量,并分配至各AGC机组执行。本发明充分发扬了互联电网控制区之间的功率互援作用,解决了局部电网调节能 力 欠缺而造成的风电无法消纳的问题。 | ||||||
权利要求 | 1.一种风电集中接入后的互联电网分布式AGC控制方法,其特征在于,所述方法包括以下步骤: |
||||||
说明书全文 | 一种风电集中接入后的互联电网分布式AGC控制方法技术领域[0001] 本发明涉及一种控制方法,具体设计一种风电集中接入后的互联电网分布式AGC控制方法。 背景技术[0002] 目前,风力发电在我国已进入快速发展阶段,并网容量和装机容量均已位居世界第一。大规模风电并网后,风电出力的随机性和间歇性特性使得互联电网负荷频率(LFC)控制不仅需要平衡负荷波动功率,同时也需要平衡来自风电的“负”负荷不规则的波动,电网的源荷双侧波动给电网的负荷频率(LFC)的AGC控制带来新的挑战。 [0003] 目前的AGC控制采用分区的控制模式。在互联电网中,各控制区的AGC控制模式通常采用联络线频率偏差控制(TBC),以维持本控制区有功功率的就地平衡。但是风电受到自然条件的束缚,风电场一般分布在互联电网中风资源丰富的地区,风电的分布也就往往集中在少数的某些控制区,例如中国的风电主要集中在三北地区以及沿海的一些地区。风电接入互联电网局部区域后,局部区域的区域控制偏差ACE变化范围大,使得局部电网的AGC调节资源显得相对紧张,依靠本局部控制区内的AGC往往很难平衡风电的大幅度波动功率,从而给风电的消纳以及互联电网频率带来很大的影响。并且更特别的,风电的快速大范围波动需要快速的调节资源与之匹配,但是在一些国家如中国的风电集中区域往往缺少例如燃气等快速的AGC机组,因此更加剧了风电接入后的AGC控制难度。因此必须重新审视当前高风电渗透下的电网AGC控制方法,以提高AGC系统的控制性能。 [0004] 目前,如何确定互联电网控制区域AGC总调节功率,提升互联电网LFC的有功AGC控制性能研究备受关注。早在20世纪50年代,Kirchmayer根据经典控制理论中的传递函数原理,提出了互联系统LFC的数学模型,研究了PI控制方法;至今目前工程上仍然主要应用针对控制区域偏差ACE采用PI控制算法来确定互联电网各控制区的AGC调节功率。随着1997年,NERC提出新的CPS1和CPS2的AGC自动发电控制性能评价标准(简称CPS标准),增强了互联电网控制区之间的功率支援能力,互联电网控制区域AGC总调节功率的确定还包括了区域间功率支援增量。但是这种功率支援量是极其有限的,很难满足风电大规模集中电网下的互联电网有功功率控制的需求。因此。文献[高宗和,滕贤亮,张小白.适应大规模风电接入的互联电网有功功率控制[J].电力系统自动化,2010,34(17):37-41.]针对我国互联电网提出了集中式AGC控制的构想,通过平等地调动网内优质调频电源,以应对大规模风电并网的有功控制挑战,但该文并没有讨论具体的控制模型与算法。并且集中式的AGC控制对数据通信的可靠性提出了更高的要求,控制中心需要收集整个互联电网内的联络线功率、稳定断面潮流、AGC机组调节信息等,通讯数据量大为增加,信息交换较以往变得更为复杂。 发明内容[0005] 为了克服上述现有技术的不足,本发明提供一种风电集中接入后的互联电网分布式AGC控制方法,各控制区只需要与本控制区相关的风电信息,并不需要整个互联电网的信息。本发明充分发扬了互联电网控制区之间的功率互援作用,解决了局部电网调节能力欠缺而造成的风电无法消纳的问题。 [0006] 为了实现上述发明目的,本发明采取如下技术方案: [0007] 本发明提供一种风电集中接入后的互联电网分布式AGC控制方法,所述方法包括以下步骤: [0008] 步骤1:风电-非风电控制区协调决策器确定风电波动平衡责任系数矩阵; [0009] 步骤2:风电控制区控制器根据风电控制区承担自身控制区的风电波动平衡责任系数,确定风电控制区的AGC总调节量,并分配至各AGC机组执行; [0010] 步骤3:非风电控制区控制器根据非风电控制区承担的风电控制区的风电波动平衡责任系数,确定非风电控制区的AGC总调节量,并分配至各AGC机组执行。 [0011] 所述步骤1中,风电波动平衡责任系数矩阵用n×m矩阵用α表示,有: [0012] [0013] α由上半部分和下半部分两部分组成,上半部分矩阵为风电控制区的风电波动平衡责任系数矩阵,其为对角方阵,αmm为风电控制区m用于自身的风电波动平衡责任系数; [0014] 下半部分矩阵为非风电控制区参与平衡风电波动功率的责任系数矩阵,其为(n-m)×m矩阵,αn,m为非风电控制区n参与平衡风电控制区m的风电波动功率的责任系数。 [0015] 在每个调度周期内,α按照如下步骤获得: [0016] (1)确定上半部分矩阵中的元素,具体过程如下: [0017] 1-1)针对风电控制区i,预测下个调度周期内的风电出力期望值 和最大预测误差ΔPi,i=1,2,…,m; [0018] 1-2)寻找风电控制区i的自身控制区的风电波动平衡责任系数αii,使得|(1-αii)ΔPi|≤Li,其中αii∈[0,1],Li为风电控制区i的功率转移容许设定阈值,且: [0019] [0020] 其中,ε10为一年时段内系统实际频率与标准频率偏差的10min平均值的均方根值;Bs为互联电力系统总频率偏差系数;Bi为风电控制区i的频率偏差系数; [0021] 1-3)重新执行1-2),直至确定上半部分矩阵中的所有元素; [0022] (2)确定下半部分矩阵中的元素,具体过程如下: [0023] 3-1)逐一确定每一列元素,使每一列元素之和等于1;3-2)针对非风电控制区j,预测下个调度周期的该非风电控制区j内的负荷需求Plj,其中j=m+1,m+2,…,n;3-3)针对非风电控制区j,检验第j行元素是否满足下列约束: [0024] [0025] [0026] [0027] 其中,ΔPi为风电控制区i的风电功率预测误差极值, 和 为非风电控制区j中所有机组的最大出力和最小出力,β为旋转备用水平,Lj为非风电控制区j的功率转移容许设定阈值,且: [0028] [0029] 其中,Bj为非风电控制区j的频率偏差系数; [0030] 3-4)针对非风电控制区j,如果第j行元素出现不满足上述条件,则返回3-1),重新调整下半部分矩阵的元素,即非风电控制区j的风电波动平衡责任系数; [0031] 3-5)如果非风电控制区j的风电波动平衡责任系数整定发生超时,则增大上半部分矩阵的元素,即加大风电控制区i的风电波动平衡责任系数后,返回3-1),重新调整下半部分矩阵的元素,即非风电控制区j的风电波动平衡责任系数。 [0032] 所述步骤2具体包括以下步骤: [0033] 步骤2-1:计算风电控制区i的区域控制偏差ACEi,有: [0034] ACEi=BiΔfi+(ΣPti-ΣI0i) [0035] 其中,Bi为风电控制区i的频差系数;Δfi为风电控制区i的频率偏差,且Δfi=fi-f0,fi为当前控制周期内风电控制区i的电网频率实际测量值,f0为标称频率;ΣPti为风电控制区i中所有联络线交换功率的实际量测值之和,ΣI0i为风电控制区i与外控制区的交易计划之和; [0036] 步骤2-2:计算风电控制区i的基本调节量Pbase,有: [0037] [0038] 其中,kp为比例系数,ki为积分系数; [0039] 步骤2-3:计算风电控制区i的平衡区内风电波动的调节量Pwind; [0040] 步骤2-4:计算风电控制区i参与平衡频率波动的支援调节量Pf; [0041] 步骤2-5:计算风电控制区i的AGC总调节量Ptotal,有Ptotal=Pbase+Pwind+Pf。 [0042] 所述步骤2-3具体包括以下步骤: [0043] 步骤2-3-1:首先接收风电-非风电控制区协调决策器下发的风电控制区i的自身控制区的风电波动平衡责任系数αii; [0044] 步骤2-3-2:计算风电控制区i的风电控制偏差WCEi,有: [0045] [0046] 其中, 为在当前控制周期内风电控制区i的风电实际出力, 为在当前控制周期内风电控制区i的风电计划出力; [0047] 步骤2-3-3:计算风电控制区i的平衡区内风电波动的调节量Pwind,有: [0048] [0049] 所述步骤2-4具体包括以下步骤: [0050] 步骤2-4-1:计算风电控制区i的频率偏差平均值 有: [0051] [0052] 其中,Ni为风电控制区i的频率偏差采集样本个数; [0053] 步骤2-4-2:根据 分级确定风电控制区i参与平衡频率波动的支援调节量Pf,具体有: [0054] A)如果 则Pf=0,其中ΔFmin为频率调节门槛设定值; [0055] B)如果 则 其中kf为频率增益系数,具体有: [0056] B-1)如果Pf>Pmaxup,则Pf=Pmaxup,其中Pmaxup为风电控制区i向上调节的最大值; [0057] B-2)如果Pf<-Pmaxdown,则Pf=Pmaxdown,其中Pmaxdown为风电控制区i向下调节的最大值。 [0058] 所述步骤3具体包括以下步骤: [0059] 步骤3-1:计算非风电控制区j的区域控制偏差ACEj,有 [0060] ACEj=BjΔfj+(ΣPtj-ΣI0j) [0061] 其中,Bj为非风电控制区j的频差系数;Δfj为非风电控制区j的频率偏差,且Δfj=fj-f0,fj为当前控制周期内非风电控制区j的电网频率实际测量值;ΣPtj为非风电控制区j中所有联络线交换功率的实际量测值之和,ΣI0j为非风电控制区j与外控制区的交易计划之和; [0062] 步骤3-2:计算非风电控制区j的基本调节量P′base,有: [0063] [0064] 步骤3-3:计算非风电控制区j的平衡区内风电波动的调节量P′wind; [0065] 步骤3-4:计算非风电控制区j参与平衡频率波动的支援调节量P′f; [0066] 步骤3-5:计算非风电控制区j的AGC总调节量P′total,有P′total=P′base+P′wind+P′f。 [0067] 所述步骤3-3具体包括以下步骤: [0068] 步骤3-3-1:首先接收风电-非风电控制区协调决策器下发的α的第j行元素,组成非风电控制区j参与风电控制区的风电波动平衡责任系数αj,有αj=[αj1 αj2 ... αjm],其中,αjm为非风电控制区j参与的风电控制区m的风电波动平衡责任系数; [0069] 步骤3-3-2:接收各风电控制区发送来的风电控制偏差值,形成非风电控制区j的风电控制偏差WCEj,有: [0070] [0071] 其中, 在当前控制周期内风电控制区m的风电实际出力, 为在当前控制周期内风电控制区m的风电计划出力; [0072] 步骤3-3-3:计算非风电控制区j的平衡区内风电波动的调节量P′wind,有: [0073] [0074] 其中,WCE′j为WCEj的转置向量。 [0075] 所述步骤3-4具体包括以下步骤: [0076] 步骤3-4-1:计算频率偏差平均值 有: [0077] [0078] 其中,Nj为非风电控制区j的频率偏差采集样本个数; [0079] 步骤3-4-2:根据 分级确定非风电控制区j参与平衡频率波动的支援调节量P′f,具体有: [0080] A)如果 则P′f=0; [0081] B)如果 则 具体有: [0082] B-1)如果P′f>P′ maxup,则P′f=P′ maxup,其中P′maxup为非风电控制区j向上调节的最大值; [0083] B-2)如果P′f<-P′ maxdown,则P′f=P′ maxdown,其中P′maxdown为非风电控制区j向下调节的最大值。 [0084] 与现有技术相比,本发明的有益效果在于: [0085] 1、综合考虑了频率恢复与风电波动功率平衡的双重问题,在当风电集中的控制区域风电波动恶化了电网频率品质时,本发明方法利用非风电集中控制区的调节资源去支援风电集中的控制区域;当风电集中的控制区风电波动增强电网频率品质时,本方明方法利用风电集中的控制区域的有益波动功率去支援非风电集中的控制区域。 [0086] 2、鼓励调动全网调节资源去平抑风电的波动,提升了电网对风电的消纳能力。 [0087] 3、便于电网选择性价比好的机组参与风电波动的平抑,降低了电网的平抑风电波动的成本,提升了电网运行的经济性。 [0089] 图1是本发明实施例中分布式AGC控制器结构示意图; [0090] 图2是本发明实施例中风电波动平衡责任系数矩阵流程图。 具体实施方式[0091] 下面结合附图对本发明作进一步详细说明。 [0092] 如图1,分布式AGC控制器主要由风电-非风电控制区协调决策器、风电控制区控制器、非风电控制区控制器等三部分组成。 [0093] 如图2,本发明提供一种风电集中接入后的互联电网分布式AGC控制方法,所述方法包括以下步骤: [0094] 步骤1:风电-非风电控制区协调决策器确定风电波动平衡责任系数矩阵; [0095] 步骤2:风电控制区控制器根据风电控制区承担自身控制区的风电波动平衡责任系数,确定风电控制区的AGC总调节量,并分配至各AGC机组执行; [0096] 步骤3:非风电控制区控制器根据非风电控制区承担的风电控制区的风电波动平衡责任系数,确定非风电控制区的AGC总调节量,并分配至各AGC机组执行。 [0097] 所述步骤1中,风电波动平衡责任系数矩阵用n×m矩阵用α表示,有: [0098] [0099] α由上半部分和下半部分两部分组成,上半部分矩阵为风电控制区的风电波动平衡责任系数矩阵,其为对角方阵,αmm为风电控制区m用于自身的风电波动平衡责任系数; [0100] 下半部分矩阵为非风电控制区参与平衡风电波动功率的责任系数矩阵,其为(n-m)×m矩阵,αn,m为非风电控制区n参与平衡风电控制区m的风电波动功率的责任系数。 [0101] 在每个调度周期内,α按照如下步骤获得: [0102] (1)确定上半部分矩阵中的元素,具体过程如下: [0103] 1-1)针对风电控制区i,预测下个调度周期内的风电出力期望值 和最大预测误差ΔPi,i=1,2,…,m; [0104] 1-2)寻找风电控制区i的自身控制区的风电波动平衡责任系数αii,使得|(1-αii)ΔPi|≤Li,其中αii∈[0,1],Li为风电控制区i的功率转移容许设定阈值,且: [0105] [0106] 其中,ε10为一年时段内系统实际频率与标准频率偏差的10min平均值的均方根值;Bs为互联电力系统总频率偏差系数;Bi为风电控制区i的频率偏差系数; [0107] 1-3)重新执行1-2),直至确定上半部分矩阵中的所有元素; [0108] (2)确定下半部分矩阵中的元素,具体过程如下: [0109] 3-1)逐一确定每一列元素,使每一列元素之和等于1;3-2)针对非风电控制区j,预测下个调度周期的该非风电控制区j内的负荷需求Plj,其中j=m+1,m+2,…,n;3-3)针对非风电控制区j,检验第j行元素是否满足下列约束: [0110] [0111] [0112] [0113] 其中,ΔPi为风电控制区i的风电功率预测误差极值, 和 为非风电控制区j中所有机组的最大出力和最小出力,β为旋转备用水平,Lj为非风电控制区j的功率转移容许设定阈值,且: [0114] [0115] 其中,Bj为非风电控制区j的频率偏差系数; [0116] 3-4)针对非风电控制区j,如果第j行元素出现不满足上述条件,则返回3-1),重新调整下半部分矩阵的元素,即非风电控制区j的风电波动平衡责任系数; [0117] 3-5)如果非风电控制区j的风电波动平衡责任系数整定发生超时,则增大上半部分矩阵的元素,即加大风电控制区i的风电波动平衡责任系数后,返回3-1),重新调整下半部分矩阵的元素,即非风电控制区j的风电波动平衡责任系数。 [0118] 所述步骤2具体包括以下步骤: [0119] 步骤2-1:计算风电控制区i的区域控制偏差ACEi,有: [0120] ACEi=BiΔfi+(ΣPti-ΣI0i) [0121] 其中,Bi为风电控制区i的频差系数;Δfi为风电控制区i的频率偏差,且Δfi=fi-f0,fi为当前控制周期内风电控制区i的电网频率实际测量值,f0为标称频率;ΣPti为风电控制区i中所有联络线交换功率的实际量测值之和,ΣI0i为风电控制区i与外控制区的交易计划之和; [0122] 步骤2-2:计算风电控制区i的基本调节量Pbase,有: [0123] [0124] 其中,kp为比例系数,ki为积分系数; [0125] 步骤2-3:计算风电控制区i的平衡区内风电波动的调节量Pwind; [0126] 步骤2-4:计算风电控制区i参与平衡频率波动的支援调节量Pf; [0127] 步骤2-5:计算风电控制区i的AGC总调节量Ptotal,有Ptotal=Pbase+Pwind+Pf。 [0128] 所述步骤2-3具体包括以下步骤: [0129] 步骤2-3-1:首先接收风电-非风电控制区协调决策器下发的风电控制区i的自身控制区的风电波动平衡责任系数αii; [0130] 步骤2-3-2:计算风电控制区i的风电控制偏差WCEi,有: [0131] [0132] 其中, 为在当前控制周期内风电控制区i的风电实际出力, 为在当前控制周期内风电控制区i的风电计划出力; [0133] 步骤2-3-3:计算风电控制区i的平衡区内风电波动的调节量Pwind,有: [0134] [0135] 所述步骤2-4具体包括以下步骤: [0136] 步骤2-4-1:计算风电控制区i的频率偏差平均值 有: [0137] [0138] 其中,Ni为风电控制区i的频率偏差采集样本个数; [0139] 步骤2-4-2:根据 分级确定风电控制区i参与平衡频率波动的支援调节量Pf,具体有: [0140] A)如果 则Pf=0,其中ΔFmin为频率调节门槛设定值; [0141] B)如果 则 其中kf为频率增益系数,具体有: [0142] B-1)如果Pf>Pmaxup,则Pf=Pmaxup,其中Pmaxup为风电控制区i向上调节的最大值; [0143] B-2)如果Pf<-Pmaxdown,则Pf=Pmaxdown,其中Pmaxdown为风电控制区i向下调节的最大值。 [0144] 所述步骤3具体包括以下步骤: [0145] 步骤3-1:计算非风电控制区j的区域控制偏差ACEj,有 [0146] ACEj=BjΔfj+(ΣPtj-ΣI0j) [0147] 其中,Bj为非风电控制区j的频差系数;Δfj为非风电控制区j的频率偏差,且Δfj=fj-f0,fj为当前控制周期内非风电控制区j的电网频率实际测量值;ΣPtj为非风电控制区j中所有联络线交换功率的实际量测值之和,ΣI0j为非风电控制区j与外控制区的交易计划之和; [0148] 步骤3-2:计算非风电控制区j的基本调节量P′base,有: [0149] [0150] 步骤3-3:计算非风电控制区j的平衡区内风电波动的调节量P′wind; [0151] 步骤3-4:计算非风电控制区j参与平衡频率波动的支援调节量P′f; [0152] 步骤3-5:计算非风电控制区j的AGC总调节量P′total,有P′total=P′base+P′wind+P′f。 [0153] 所述步骤3-3具体包括以下步骤: [0154] 步骤3-3-1:首先接收风电-非风电控制区协调决策器下发的α的第j行元素,组成非风电控制区j参与风电控制区的风电波动平衡责任系数αj,有αj=[αj1 αj2 ... αjm],其中,αjm为非风电控制区j参与的风电控制区m的风电波动平衡责任系数; [0155] 步骤3-3-2:接收各风电控制区发送来的风电控制偏差值,形成非风电控制区j的风电控制偏差WCEj,有: [0156] [0157] 其中, 在当前控制周期内风电控制区m的风电实际出力, 为在当前控制周期内风电控制区m的风电计划出力; [0158] 步骤3-3-3:计算非风电控制区j的平衡区内风电波动的调节量P′wind,有: [0159] [0160] 其中,WCE′j为WCEj的转置向量。 [0161] 所述步骤3-4具体包括以下步骤: [0162] 步骤3-4-1:计算频率偏差平均值 有: [0163] [0164] 其中,Nj为非风电控制区j的频率偏差采集样本个数; [0165] 步骤3-4-2:根据 分级确定非风电控制区j参与平衡频率波动的支援调节量P′f,具体有: [0166] A)如果 则P′f=0; [0167] B)如果 则 具体有: [0168] B-1)如果P′f>P′ maxup,则P′f=P′ maxup,其中P′maxup为非风电控制区j向上调节的最大值; [0169] B-2)如果P′f<-P′ maxdown,则P′f=P′ maxdown,其中P′maxdown为非风电控制区j向下调节的最大值。 |