CN112865132B - 双区域互联电力系统的负荷频率控制参数的处理方法 - Google Patents
双区域互联电力系统的负荷频率控制参数的处理方法 Download PDFInfo
- Publication number
- CN112865132B CN112865132B CN202011640578.5A CN202011640578A CN112865132B CN 112865132 B CN112865132 B CN 112865132B CN 202011640578 A CN202011640578 A CN 202011640578A CN 112865132 B CN112865132 B CN 112865132B
- Authority
- CN
- China
- Prior art keywords
- fractional
- pid controller
- fractional order
- controller
- power system
- 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
Links
- 238000003672 processing method Methods 0.000 title abstract description 5
- 238000000034 method Methods 0.000 claims abstract description 23
- 238000012545 processing Methods 0.000 claims description 20
- 238000004422 calculation algorithm Methods 0.000 claims description 17
- 238000004088 simulation Methods 0.000 claims description 15
- 238000012546 transfer Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 230000006870 function Effects 0.000 description 28
- 238000010586 diagram Methods 0.000 description 23
- 239000002245 particle Substances 0.000 description 15
- 230000009977 dual effect Effects 0.000 description 8
- 230000015654 memory Effects 0.000 description 8
- 238000003860 storage Methods 0.000 description 8
- UUUHXMGGBIUAPW-UHFFFAOYSA-N 1-[1-[2-[[5-amino-2-[[1-[5-(diaminomethylideneamino)-2-[[1-[3-(1h-indol-3-yl)-2-[(5-oxopyrrolidine-2-carbonyl)amino]propanoyl]pyrrolidine-2-carbonyl]amino]pentanoyl]pyrrolidine-2-carbonyl]amino]-5-oxopentanoyl]amino]-3-methylpentanoyl]pyrrolidine-2-carbon Chemical compound C1CCC(C(=O)N2C(CCC2)C(O)=O)N1C(=O)C(C(C)CC)NC(=O)C(CCC(N)=O)NC(=O)C1CCCN1C(=O)C(CCCN=C(N)N)NC(=O)C1CCCN1C(=O)C(CC=1C2=CC=CC=C2NC=1)NC(=O)C1CCC(=O)N1 UUUHXMGGBIUAPW-UHFFFAOYSA-N 0.000 description 7
- 102000004270 Peptidyl-Dipeptidase A Human genes 0.000 description 7
- 108090000882 Peptidyl-Dipeptidase A Proteins 0.000 description 7
- 238000004590 computer program Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- UDIODXADSSQKTM-GUHNCMMLSA-N (5ar,8ar,9r)-5-[[(2r,4ar,6r,7r,8r,8as)-7,8-dihydroxy-2-methyl-4,4a,6,7,8,8a-hexahydropyrano[3,2-d][1,3]dioxin-6-yl]oxy]-9-(4-hydroxy-3,5-dimethoxyphenyl)-5a,6,8a,9-tetrahydro-5h-[2]benzofuro[6,5-f][1,3]benzodioxol-8-one;(7s,9s)-7-[(2r,4s,5s,6s)-4-amino-5- Chemical compound ClCCN(CCCl)P1(=O)NCCCO1.O([C@H]1C[C@@](O)(CC=2C(O)=C3C(=O)C=4C=CC=C(C=4C(=O)C3=C(O)C=21)OC)C(=O)CO)[C@H]1C[C@H](N)[C@H](O)[C@H](C)O1.COC1=C(O)C(OC)=CC([C@@H]2C3=CC=4OCOC=4C=C3C(O[C@H]3[C@@H]([C@@H](O)[C@@H]4O[C@H](C)OC[C@H]4O3)O)[C@@H]3[C@@H]2C(OC3)=O)=C1 UDIODXADSSQKTM-GUHNCMMLSA-N 0.000 description 3
- 230000004075 alteration Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002431 foraging effect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009131 signaling function Effects 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000007723 transport mechanism Effects 0.000 description 1
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/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/241—The oscillation concerning frequency
-
- 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)
- Feedback Control In General (AREA)
Abstract
本公开是关于一种双区域互联电力系统的负荷频率控制参数的处理方法，所述方法包括：步骤1，建立分数阶PID控制器，分数阶PID控制器是基于分数阶比例积分控制器与分数阶比例微分控制器依次级联建立的控制器；步骤2，基于分数阶PID控制器，建立双区域互联电力系统模型；步骤3，基于双区域互联电力系统模型，对分数阶PID控制器参数进行优化，获得优化后的分数阶PID控制器参数；步骤4，将优化后的分数阶PID控制器参数作为双区域互联电力系统的负荷频率控制参数。本发明提供的分数阶PID控制器参数作为双区域互联电力系统的负荷频率控制参数，相比于采用其它控制器的电力系统能更好的抑止频率波动，所需要的时间更短，具有较好的动态性能。
Description
技术领域
本发明涉及电力系统调节控制技术领域，尤其涉及一种双区域互联电力系统的负荷频率控制参数的处理方法。
背景技术
电力系统的频率稳定是衡量电能质量的标准之一。随着微电网的普及，各种具有间歇性与随机性分布式能源的接入，以及随机的负荷扰动，势必会对电网的频率稳定造成一定的冲击影响，使得维持电力系统稳定的难度大大提高。在双区域互联电力系统中，任意一侧的负荷扰动都会导致两个区域的频率以及电力系统联络线上功率偏离设定的额定值。
为了解决上述问题，采用的手段是负荷频率控制(load frequency control，LFC)技术，这是一种准确高效快速的控制机制。当系统在遭受外界因素的冲击时，运用该控制机制可以快速的减少区域瞬时频率偏差，保持区域电力系统频率和联络线交换功率在设定值。
近年来，有众多专家学者在电力系统LFC问题上进行理论与技术的研究。LFC控制器涉及到两方面的内容，即控制策略的选择、控制器参数的选取与调整。许多先进的控制理论应用于LFC系统控制器中，如鲁棒控制、自适应控制、强化学习控制、模型预测控制以及结合人工智能理论的神经网络等。虽然这些方法在一定程度上提高了系统的稳态性能，但是由于其实现方法比较困难，控制器结构复杂，所以大多数的工业企业仍采用更为传统的控制器，如PID控制器。其有适用范围广，易于在数字平台上实现等优点。随着数学的发展，在科学与工程领域发现了分数阶PID控制器潜在的应用。如何快速的使系统达到稳定状态并将稳态误差保持为零，成为亟待解决的技术问题。
发明内容
为克服相关技术中存在的问题，本公开提供一种双区域互联电力系统的负荷频率控制参数的处理方法。
根据本公开实施例的一方面，提供一种双区域互联电力系统的负荷频率控制参数的处理方法，所述方法包括：
步骤1，建立分数阶PID控制器，所述分数阶PID控制器是基于分数阶比例积分控制器与分数阶比例微分控制器依次级联建立的控制器；
步骤2，基于所述分数阶PID控制器，建立双区域互联电力系统模型；
步骤3，基于所述双区域互联电力系统模型，对所述分数阶PID控制器参数进行优化，获得优化后的所述分数阶PID控制器参数；
步骤4，将优化后的所述分数阶PID控制器参数作为双区域互联电力系统的负荷频率控制参数。
其中，所述步骤1中，所述建立分数阶PID控制器包括：
步骤1-1，采用Oustaloup滤波算法构造连续滤波器；
设拟合频率段为(ωb,ωh)，则连续滤波器的传递函数为：
公式(1)中，s为拉氏变换域中的变量，K为增益，N为连续滤波器的阶次，k为连续滤波器的阶次的初始值；
其中，连续滤波器的零点为ω′k，极点为ωk，可由公式(2)计算得出：
其中，γ为分数阶的阶次，γ>0时，为微分运算，γ<0时，为积分运算；N为滤波器的阶次，ωb和ωh分别为拟合频率的上下限；
步骤1-2，利用连续滤波器对区域控制偏差进行分数阶微积分处理；
步骤1-3，将分数阶微积分处理后的区域控制偏差作为分数阶比例积分控制器的输入，将分数阶比例积分控制器的输出作为分数阶比例微分控制器的输入。
其中，所述步骤2中，基于所述分数阶PID控制器，建立双区域互联电力系统模型包括：
步骤2-1，设置所述双区域互联电力系统模型中两个区域的区域控制偏差分别为：
其中，ACE1、ACE2分别是两个区域的区域控制偏差，Δf1、Δf2分别为两个区域的频率偏差，B1、B2分别为两个区域的频率偏差因子，ΔP12为联络线功率偏差；
步骤2-2，设置所述双区域互联电力系统模型中的所述分数阶PID控制器的输出信号为：
其中，u1(s)和u2(s)分别是两个区域的控制器的输出信号，GFOPI+FOPD(s)为所述分数阶PID控制器的传递函数，其中，
其中，u(s)为所述分数阶PID控制器输出的传递函数，e(s)为所述分数阶PID控制器输入的传递函数，Kp、Ki、Kd分别为所述分数阶PID控制器的比例增益、积分增益和微分增益，s为拉氏变换域中的变量，λ与μ为所述分数阶PID控制器的积分阶次与微分阶次；
步骤2-3，控制所述分数阶PID控制器的输出信号u1(s)和u2(s)，以调整区域控制偏差ACE1和ACE2至预设范围。
其中，在所述步骤3中，所述基于所述双区域互联电力系统模型，对所述分数阶PID控制器参数进行优化，获得优化后的所述分数阶PID控制器参数包括：
步骤3-1，设置模型的约束条件
所述约束条件包括分数阶比例积分控制器的比例增益的上下阈值、分数阶比例微分控制器的比例增益的上下阈值、分数阶PID控制器的积分增益的上下阈值、分数阶PID控制器的微分增益的上下阈值、分数阶PID控制器的积分阶次的上下阈值和分数阶PID控制器的微分增益的上下阈值；
其约束条件为：
其中，Kp1为分数阶比例积分控制器的比例增益，Kp2为分数阶比例微分控制器的比例增益，Ki、Kd分别为分数阶PID控制器的积分增益和微分增益，λ与μ为分数阶PID控制器的积分阶次与微分阶次，为Kp1的上下阈值，/>为Kp2的上下阈值，是Ki的上下阈值，/> 是Kd的上下阈值，λmin,λmax是λ的上下阈值，μmin,μmax是μ的上下阈值；
步骤3-2，确定双区域互联电力系统模型的目标函数
以双区域互联电力系统的性能参数指标ITAE的目标函数，其数学表达式具体如下：
其中，Δf1、Δf2分别为两个区域的频率偏差，ΔP12为联络线功率偏差，tsim是仿真持续时间，t是时间变量；
步骤3-3，在满足约束条件的情况下，输出双区域互联电力系统模型的目标函数的最小值，确定优化后的所述分数阶PID控制器参数。
本公开的实施例提供的技术方案，基于PI与PD级联控制器的基础上提出分数阶PID控制器，通过对分数阶PID控制器参数进行优化确定最优的分数阶PID控制器参数，将优化后的分数阶PID控制器参数作为双区域互联电力系统的负荷频率控制参数，相比于采用其它控制器的电力系统能更好的抑止频率波动，所需要的时间更短，具有较好的动态性能。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本公开。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本发明的实施例，并与说明书一起用于解释本发明的原理。
图1是相关技术中示出的一种常规的分数阶PID控制的框图；
图2是根据一示例性实施例示出的一种双区域互联电力系统的负荷频率控制参数的处理方法的流程图；
图3是根据一示例性实施例示出的分数阶PID控制器的结构示意图；
图4是根据一示例性实施例示出的双区域互联电力系统模型的框图；
图5是根据一示例性实施例示出的确定优化分数阶PID控制器参数的流程图；
图6是根据一示例性实施例示出的双区域互联电力系统的仿真图；
图7是根据一示例性实施例示出的分数阶PID控制器的仿真图；
图8是根据一示例性实施例示出的Csse1频率偏差Δf1仿真图；
图9是根据一示例性实施例示出的Case2频率偏差Δf1仿真图；
图10是本根据一示例性实施例示出的Case3频率偏差Δf1仿真图；
图11是根据一示例性实施例示出的一种计算机设备的框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。
图1是根据一示例性实施例示出的一种双区域互联电力系统的负荷频率控制参数的处理方法的流程图。
如图1所示，本发明提供一种双区域互联电力系统的负荷频率控制参数的处理方法，该方法包括：
步骤101，建立分数阶PID控制器，分数阶PID控制器是基于分数阶比例积分控制器与分数阶比例微分控制器依次级联建立的控制器。
该步骤中，该分数阶PID控制器是基于PI与PD级联控制器的基础上提出的，该分数阶PID控制器既具有PID控制器的优点，又有单独使用PI控制器与PD控制器时的特征，通过调节分数阶PID控制器参数能够获得更好的系统性能，本发明提出一种新的分数阶PID控制器来稳定电力系统的频率，以抑止外界扰动及负荷扰动导致的系统频率的剧烈波动。
步骤102，基于分数阶PID控制器，建立双区域互联电力系统模型；
步骤103，基于双区域互联电力系统模型，对分数阶PID控制器参数进行优化，获得优化后的分数阶PID控制器参数；
步骤104，将优化后的分数阶PID控制器参数作为双区域互联电力系统的负荷频率控制参数。
本实施例中，分数阶PID控制器为FOPI+FOPD控制器，既具有PID控制器的优点，又有单独使用PI控制器与PD控制器时的特征，通过对分数阶PID控制器参数进行优化确定最优的分数阶PID控制器参数，将优化后的分数阶PID控制器参数作为双区域互联电力系统的负荷频率控制参数，保持双区域发电功率间的动态平衡，维持整个互联电力系统的频率恒定，相比于采用其它控制器的电力系统能更好的抑止频率波动，所需要的时间更短，具有较好的动态性能。
在示例性实施例中，如图3所示，分数阶PID控制器是基于分数阶比例积分(Fractional Order Proportional Integral，FOPI)控制器与分数阶比例微分(Fractional Order Proportional Differential，FOPD)控制器依次级联建立的FOPI+FOPD控制器。
在步骤101中，建立分数阶PID控制器包括：步骤1-1，采用Oustaloup滤波算法构造连续滤波器；
设拟合频率段为(ωb,ωh)，则连续滤波器的传递函数为：
公式(1)中，s为拉氏变换域中的变量，K为增益，N为连续滤波器的阶次，k为连续滤波器的阶次的初始值；
其中，连续滤波器的零点为ω′k，极点为ωk，可由公式(2)计算得出：
其中，γ为分数阶的阶次，γ>0时，为微分运算，γ<0时，为积分运算；N为滤波器的阶次，ωb和ωh分别为拟合频率的上下限；
其中，连续滤波器的阶次N值的大小一般不宜过大，否则会加大近似时间，从而影响整个系统的优化时间，ωb和ωh参数选取结果会影响整个系统的寻优过程，过大的ωb和ωh会导致整个系统的寻优过程过慢。
步骤1-2，利用连续滤波器对区域控制偏差进行分数阶微积分处理；
该步骤中，采用Oustaloup滤波算法构造的连续滤波器可以实现对函数进行分数阶的微分或者积分的计算，本实施例中，根据上述公式(1)和公式(2)在可以MATLAB中编写相应的函数构造出连续滤波器，然后将区域控制偏差信号函数通过构造的连续滤波器进行分数阶微积分处理。
步骤1-3，将连续滤波器处理后的区域控制偏差作为分数阶比例积分控制器的输入，将分数阶比例积分控制器的输出作为分数阶比例微分控制器的输入。
参见图3所示，该步骤中，区域控制偏差经过连续滤波器进行分数阶微积分处理后，将处理后的区域控制偏差输入到分数阶比例积分控制器(FOPI)进行处理，将分数阶比例积分控制器的输出输入到分数阶比例微分控制器(FOPD)中进行处理，实现了在比例环节对区域控制偏差进行放大，通过调节FOPD与FOP连个控制器的参数，使FOPD控制器输出的值降到最小，此时，所得到的FOPI+FOPD级联控制器为建立的分数阶PID控制器(FOPI+FOPD控制器)。
在示例性实施例中，如图4所示，双区域互联电力系统模型包括两个通过联络线互联的系统区域，两个系统区域分别存在于独立的电力系统区域，且两电力系统区域的电力输出端互联，两个独立的电力系统区域分别设为区域一及区域二中，其中，区域一中还包括调速器模型gov1、汽轮机模型tur1和发电机模型gen1。区域二中还包括调速器模型gov2、汽轮机模型tur2和发电机模型gen2。
在步骤102中，基于分数阶PID控制器，建立双区域互联电力系统模型包括：
步骤2-1，设置双区域互联电力系统模型中两个区域的区域控制偏差分别为：
公式(3)中，ACE1、ACE2分别是两个区域的区域控制偏差，Δf1、Δf2分别为两个区域的频率偏差，B1、B2分别为两个区域的频率偏差因子，ΔP12为联络线功率偏差。
步骤2-2，设置双区域互联电力系统模型中的分数阶PID控制器的输出信号为：
公式(4)中，u1(s)和u2(s)分别是两个区域的控制器的输出信号，GFOPI+FOPD(s)为分数阶PID控制器的传递函数，其中，
公式(5)中，其中u(s)为分数阶PID控制器输出的传递函数，e(s)为分数阶PID控制器输入的传递函数，Kp、Ki、Kd分别为分数阶PID控制器的比例增益、积分增益和微分增益，s为拉氏变换域中的变量，λ与μ为分数阶PID控制器的积分阶次与微分阶次。
GFOPI+FOPD(s)为分数阶PID控制器(FOPI+FOPD控制器)的传递函数，当出现负荷扰动时，两个区域的频率就会出现偏差，并且联络线上也会出现功率偏差，此时整个电力系统会对这些偏差进行抑止，即根据区域控制偏差的值所处的控制区间，对输出功率发生相应调整，使系统频率偏差保持在正常范围内。
步骤2-3，控制分数阶PID控制器的输出信号u1(s)和u2(s)，以调整区域控制偏差ACE1和ACE2至预设范围，例如调整区域控制偏差为零。
在步骤103中，基于双区域互联电力系统模型，对分数阶PID控制器参数进行优化，获得优化后的分数阶PID控制器参数包括：
步骤3-1，设置模型的约束条件
约束条件包括分数阶比例积分控制器的比例增益的上下阈值、分数阶比例微分控制器的比例增益的上下阈值、分数阶PID控制器的积分增益的上下阈值、分数阶PID控制器的微分增益的上下阈值、分数阶PID控制器的积分阶次的上下阈值和分数阶PID控制器的微分增益的上下阈值；
其约束条件为：
其中，Kp1为分数阶比例积分控制器的比例增益，Kp2为分数阶比例微分控制器的比例增益，Ki、Kd分别为分数阶PID控制器的积分增益和微分增益，λ与μ为分数阶PID控制器的积分阶次与微分阶次，为Kp1的上下阈值，/>为Kp2的上下阈值，是Ki的上下阈值，/> 是Kd的上下阈值，λmin,λmax是λ的上下阈值，μmin,μmax是μ的上下阈值。
步骤3-2，确定双区域互联电力系统模型的目标函数；
以双区域互联电力系统的性能参数指标ITAE的目标函数，其数学表达式具体如下：
其中，Δf1、Δf2分别为两个区域的频率偏差，ΔP12为联络线功率偏差，tsim是仿真持续时间，t是时间变量。
步骤3-3，在满足约束条件的情况下，输出双区域互联电力系统模型的目标函数的最小值，确定优化后的分数阶PID控制器参数。
其中，分数阶PID控制器参数包括：分数阶比例积分控制器的比例增益Kp1、分数阶比例微分控制器的比例增益Kp2、积分增益Ki、微分增益Kd、积分阶次λ与微分阶次μ。
本实施例中，选取时间乘误差绝对值积分(Integral of Time multipliedAbsolute Error，ITAE)性能指标作为目标函数，，ITAE是一种衡量系统性能的指标，目标函数为按照这个目标函数值最小为准则设计参数的系统能够达到很好的效果，该性能指标能够提供更好的时域动态响应，该目标函数包括两个区域的频率偏差以及功率偏差，选取该指标作为目标函数，可通过粒子群算法对分数阶PID控制器参数进行优化，使目标函数达到最小，则区域的频率偏差以及功率偏差亦是最小，即通过对分数阶PID控制器参数进行优化，确定最优的Kp1,Kp2,Ki,Kd,λ,μ参数，将优化后的分数阶PID控制器参数作为双区域互联电力系统的负荷频率控制参数，达到了对双区域互联电力系统频率进行控制的目的。
参照图3中的分数阶PID控制器参数Kp1,Kp2,Ki,Kd,λ,μ是可根据上述过程确定的最优分数阶PID控制器参数。
在示例性实施例中，在MATLAB中的SIMULINK环境下对双区域互联电力系统进行仿真建模，参见图3、图5、图6和图7所示，通过粒子群算法基于双区域互联电力系统模型，对分数阶PID控制器参数进行优化，获得优化后的分数阶PID控制器参数过程如下：
步骤501，对分数阶PID控制器参数进行初始化，设置分数阶PID控制器参数的上下阈值。
本实施例中，采用粒子群算法来求解分数阶PID控制器参数的最优控制器参数。初始化每个粒子的位置与速度，即对KP1,KP2,Ki,Kd,λ,μ进行初始化，初始化每个粒子的最优位置，初始化全局最优位置，计算每个粒子的适应度，并分别设置KP1,KP2,Ki,Kd,λ,μ的上下阈值
粒子群算法，也称粒子群优化算法或鸟群觅食算法(Particle SwarmOptimization)，缩写为PSO，PSO算法属于进化算法的一种，和模拟退火算法相似，它也是从随机解出发，通过迭代寻找最优解，它也是通过适应度来评价解的品质，但它比遗传算法规则更为简单，它没有遗传算法的“交叉(Crossover)”和“变异(Mutation)”操作，它通过追随当前搜索到的最优值来寻找全局最优。这种算法以其实现容易、精度高、收敛快等优点引起了学术界的重视，并且在解决实际问题中展示了其优越性。
步骤502，将初始化的分数阶PID控制器参数赋值到分数阶PID控制器中，运行双区域互联电力系统的仿真图，确定性能指标参数ITAE。
在示例性实施例中，如图7所示，为SIMULINK环境下的分数阶PID控制器的仿真图，并返回到MATLAB工作空间中运行双区域互联电力系统的仿真图，如图6所示，为MATLAB工作空间中的双区域互联电力系统的仿真图。
步骤503，判断性能指标参数ITAE是否满足预设条件，若不满足，则执行步骤504；若满足，执行步骤505。
其中，预设条件为以下条件的任意一个；
条件一，输出的性能参数指标ITAE低于所要求的第一预设值；
条件二，迭代次数已经达到第二预设值。
性能指标参数ITAE满足预设条件上述条件的任意一个，也就是目标函数达到了到最小值。第一预设值和第二预设值可以根据情况具体设置，例如第一预设值为0.1，第二预设值为100，在此不做限定。
性能指标参数ITAE不满足预设条件，也即是输出的性能参数指标ITAE未低于所要求的第一预设值并且迭代次数未达到第二预设值。
步骤504，更新分数阶PID控制器参数，执行步骤502。
该步骤中，更新分数阶PID控制器参数，即对粒子群进行更新操作，重新计算粒子群中粒子的值，包括更新每个粒子的速度与位置，评估每个粒子的适应度，更新每个粒子的历史最优位置，更新群体的全局最优位置。
步骤505，输出优化后的分数阶PID控制器参数。
该步骤中，通过PSO算法确定性能指标参数ITAE达到最小值时KP1,KP2,Ki,Kd,λ,μ值，即为获得的优化后的分数阶PID控制器参数。
结合图4、图8、图9和图10，本发明提供了双区域互联电力系统采用FOPI+FOPD控制器相对于双区域互联电力系统采用其他控制器对频率波动的抑止情况的仿真实验。
为了对比所提出的双区域互联电力系统采用FOPI+FOPD控制器对频率波动的抑止有效性，对比双区域互联电力系统分别采用PID、PI+PD、TID控制器与其比较，设计Case1、Case2和Case3三种实验仿真如下：
Case1在区域一中内加入+20％负荷扰动，系统的响应频率偏差Δf1曲线如图8所示，在Case1中，两个区域的控制器都采用FOPI+FOPD控制器，并且在区域一中加入+20％负荷扰动，观察其频率以及联络线功率偏差变化的情况。
为了与双区域互联电力系统采用FOPI+FOPD控制器对频率波动的抑止情况对比，随后又将双区域互联电力系统的控制器依次替换为PID、PI+PD、TID控制器，并且也同样在替换后的模型中的区域一中加入+20％负荷扰动，观察频率以及联络线功率偏差变化情况。
Case2区域二中加入+20％负荷扰动的，系统的响应频率偏差Δf1曲线如图9所示，在Case2中，两个区域的控制器都采用FOPI+FOPD控制器，并且在区域二中加入+20％负荷扰动，观察其频率以及联络线功率偏差变化的情况。
为了与双区域互联电力系统采用FOPI+FOPD控制器对频率波动的抑止情况，随后又将双区域互联电力系统的控制器依次替换为PID、PI+PD、TID控制器，并且也同样在替换后的模型中的区域二中加入+20％负荷扰动，观察频率以及联络线功率偏差变化情况。
Case3区域一中加入+20％负荷扰动且区域二中加入+10％负荷扰动，系统的响应频率偏差Δf1曲线如图10所示，在Case3中，两个区域的控制器都采用FOPI+FOPD控制器，在区域一中加入+20％负荷扰动同时在区域二中加入+10％负荷扰动，观察其频率以及联络线功率偏差变化的情况。
为了与双区域互联电力系统采用FOPI+FOPD控制器对频率波动的抑止情况，随后又将双区域互联电力系统的控制器依次替换为PID、PI+PD、TID控制器，并且也在区域一中加入+20％负荷扰动同时在区域二中加入+10％负荷扰动，观察频率以及联络线功率偏差变化情况。
仿真实验结果表明其中在Case1中采用FOPI+FOPD控制器的系统的Δf1调整时间为1.863s，相比于传统PID控制器的4.896s快了3s左右。由于负荷瞬时增加，导致系统频率下降，所以Δf1、Δf2的超调量为负值。故采用FOPI+FOPD控制器相比于其它控制器能更快地使系统各区域的频率波动减小到零。
在Case3中可以看出在两个区域都发生扰动时的Δf1的调整时间与超调量大多数都比只在一个区域发生扰动时要大，但是采用的FOPI+FOPD控制器的系统还是比采用其它控制器能更好的抑止频率波动，所需要的时间更短，具有较好的动态性能。
本发明实施例中提供一种计算机设备，包括：处理器；用于存储处理器可执行指令的存储器；其中，处理器被配置为执行存储器中的可执行指令以实现上述双区域互联电力系统的负荷频率控制参数的处理方法。
图11是根据一示例性实施例示出的一种计算机设备1100的框图。例如，计算机设备1100可以被提供为一服务器。参照图11，计算机设备1100包括处理器1101，处理器的个数可以根据需要设置为一个或者多个。计算机设备900还存储器1102，用于存储可由处理器901的执行的指令，例如应用程序。存储器的个数可以根据需要设置一个或者多个。其存储的应用程序可以为一个或者多个。处理器1101被配置为执行指令，以执行上述双区域互联电力系统的负荷频率控制参数的处理方法。
本发明还提供了一种包括指令的非临时性计算机可读存储介质，其上存储有可执行指令，该可执行指令被处理器执行时实现上述消除镜头撞击噪声的方法。例如，非临时性计算机可读存储介质可以是ROM、随机存取存储器(RAM)、CD-ROM、磁带、软盘和光数据存储设备等。
本领域技术人员应明白，本发明的实施例可提供为方法、装置(设备)、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质上实施的计算机程序产品的形式。计算机存储介质包括在用于存储信息(诸如计算机可读指令、数据结构、程序模块或其他数据)的任何方法或技术中实施的易失性和非易失性、可移除和不可移除介质,包括但不限于RAM、ROM、EEPROM、闪存或其他存储器技术、CD-ROM、数字多功能盘(DVD)或其他光盘存储、磁盒、磁带、磁盘存储或其他磁存储装置、或者可以用于存储期望的信息并且可以被计算机访问的任何其他的介质等。此外，本领域普通技术人员公知的是，通信介质通常包含计算机可读指令、数据结构、程序模块或者诸如载波或其他传输机制之类的调制数据信号中的其他数据，并且可包括任何信息递送介质。
本发明是参照根据本发明实施例的方法、装置(设备)和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括……”限定的要素，并不排除在包括所述要素的物品或者设备中还存在另外的相同要素。
尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明的意图也包含这些改动和变型在内。
Claims (2)
1.一种双区域互联电力系统的负荷频率控制参数的处理方法，其特征在于，所述方法包括：
步骤1，建立分数阶PID控制器，所述分数阶PID控制器是基于分数阶比例积分控制器与分数阶比例微分控制器依次级联建立的控制器；
步骤2，基于所述分数阶PID控制器，建立双区域互联电力系统模型；
步骤3，基于所述双区域互联电力系统模型，对所述分数阶PID控制器参数进行优化，获得优化后的所述分数阶PID控制器参数；
步骤4，将优化后的所述分数阶PID控制器参数作为双区域互联电力系统的负荷频率控制参数；
其中，所述步骤1中，所述建立分数阶PID控制器包括：
步骤1-1，采用Oustaloup滤波算法构造连续滤波器；
设拟合频率段为，则连续滤波器的传递函数为：
公式（1）中，s为拉氏变换域中的变量，K为增益，为连续滤波器的阶次，k为连续滤波器的阶次的初始值；
其中，连续滤波器的零点为，极点为/>，可由公式（2）计算得出：
（2）
其中，，/>为分数阶的阶次，/>时，为微分运算，/>时，为积分运算；/>为滤波器的阶次，/>和/>分别为拟合频率的上下限；
步骤1-2，利用连续滤波器对区域控制偏差进行分数阶微积分处理；
步骤1-3，将分数阶微积分处理后的区域控制偏差作为分数阶比例积分控制器的输入，将分数阶比例积分控制器的输出作为分数阶比例微分控制器的输入；
所述步骤2中，基于所述分数阶PID控制器，建立双区域互联电力系统模型包括：
步骤2-1，设置所述双区域互联电力系统模型中两个区域的区域控制偏差分别为：
（3）
其中，分别是两个区域的区域控制偏差，/>、/>分别为两个区域的频率偏差，/>分别为两个区域的频率偏差因子，/>为联络线功率偏差；
步骤2-2，设置所述双区域互联电力系统模型中的所述分数阶PID控制器的输出信号为：
（4）
其中，和/>分别是两个区域的控制器的输出信号，/>为所述分数阶PID控制器的传递函数，其中，
(/>)（5）
其中，为所述分数阶PID控制器输出的传递函数，/>为所述分数阶PID控制器输入的传递函数，/>、/>分别为所述分数阶PID控制器的比例增益、积分增益和微分增益，s为拉氏变换域中的变量，/>与/>为所述分数阶PID控制器的积分阶次与微分阶次；
步骤2-3，控制所述分数阶PID控制器的输出信号和/>，以调整区域控制偏差和/>至预设范围。
2.根据权利要求1所述的方法，其特征在于，在所述步骤3中，基于所述双区域互联电力系统模型，对所述分数阶PID控制器参数进行优化，获得优化后的所述分数阶PID控制器参数包括：
步骤3-1，设置模型的约束条件
所述约束条件包括分数阶比例积分控制器的比例增益的上下阈值、分数阶比例微分控制器的比例增益的上下阈值、分数阶PID控制器的积分增益的上下阈值、分数阶PID控制器的微分增益的上下阈值、分数阶PID控制器的积分阶次的上下阈值和分数阶PID控制器的微分增益的上下阈值；
其约束条件为：
（6）
其中，为分数阶比例积分控制器的比例增益，/>为分数阶比例微分控制器的比例增益，/>、/>分别为分数阶PID控制器的积分增益和微分增益，/>与/>为分数阶PID控制器的积分阶次与微分阶次，/>、/>为/>的上下阈值，/>、/>为/>的上下阈值，是/>的上下阈值，/>是/>的上下阈值，/>是/>的上下阈值，是/>的上下阈值；
步骤3-2，确定双区域互联电力系统模型的目标函数
以双区域互联电力系统的性能参数指标作为目标函数，其数学表达式具体如下：
（7）
其中，、/>分别为两个区域的频率偏差，/>为联络线功率偏差，/>是仿真持续时间，/>是时间变量；
步骤3-3，在满足约束条件的情况下，输出双区域互联电力系统模型的目标函数的最小值，确定优化后的所述分数阶PID控制器参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011640578.5A CN112865132B (zh) | 2020-12-31 | 2020-12-31 | 双区域互联电力系统的负荷频率控制参数的处理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011640578.5A CN112865132B (zh) | 2020-12-31 | 2020-12-31 | 双区域互联电力系统的负荷频率控制参数的处理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112865132A CN112865132A (zh) | 2021-05-28 |
| CN112865132B true CN112865132B (zh) | 2023-11-14 |
Family
ID=76000719
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011640578.5A Active CN112865132B (zh) | 2020-12-31 | 2020-12-31 | 双区域互联电力系统的负荷频率控制参数的处理方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112865132B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114123249B (zh) * | 2021-12-13 | 2023-07-14 | 华能(浙江)能源开发有限公司清洁能源分公司 | 一种基于电池储能主动响应的含风电互联电力系统负荷频率控制方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108646546A (zh) * | 2018-06-11 | 2018-10-12 | 燕山大学 | 确定分数阶pid控制器参数的方法、装置、电力系统稳定器及电力系统 |
| CN109193750A (zh) * | 2018-07-27 | 2019-01-11 | 湖北工业大学 | 基于pidd2控制器的风电集群参与agc方法 |
| CN110601222A (zh) * | 2019-08-19 | 2019-12-20 | 国网山东省电力公司电力科学研究院 | 一种基于Lévy飞蛾算法的分数阶负荷频率控制方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108803311A (zh) * | 2018-08-03 | 2018-11-13 | 佛山科学技术学院 | 一种pid控制器设计方法 |
-
2020
- 2020-12-31 CN CN202011640578.5A patent/CN112865132B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108646546A (zh) * | 2018-06-11 | 2018-10-12 | 燕山大学 | 确定分数阶pid控制器参数的方法、装置、电力系统稳定器及电力系统 |
| CN109193750A (zh) * | 2018-07-27 | 2019-01-11 | 湖北工业大学 | 基于pidd2控制器的风电集群参与agc方法 |
| CN110601222A (zh) * | 2019-08-19 | 2019-12-20 | 国网山东省电力公司电力科学研究院 | 一种基于Lévy飞蛾算法的分数阶负荷频率控制方法 |
Non-Patent Citations (2)
| Title |
|---|
| 新型分数阶PID控制器及其仿真研究;罗佑新;;哈尔滨工业大学学报(第05期);全文 * |
| 计及时滞的互联电网负荷频率控制最优分数阶PID控制器设计;左剑;王子琪;李银红;段献忠;;控制理论与应用(第09期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112865132A (zh) | 2021-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gheisarnejad et al. | Secondary load frequency control for multi-microgrids: HiL real-time simulation | |
| Li et al. | Optimal fractional order PID controller design for automatic voltage regulator system based on reference model using particle swarm optimization | |
| Arya | Automatic generation control of two-area electrical power systems via optimal fuzzy classical controller | |
| Pan et al. | Fractional-order load-frequency control of interconnected power systems using chaotic multi-objective optimization | |
| Prasad et al. | Optimal control of nonlinear inverted pendulum system using PID controller and LQR: performance analysis without and with disturbance input | |
| Panda et al. | Design and performance analysis of PID controller for an automatic voltage regulator system using simplified particle swarm optimization | |
| Li et al. | Offset-free fuzzy model predictive control of a boiler–turbine system based on genetic algorithm | |
| Kim et al. | Disturbance-observer-based model predictive control for output voltage regulation of three-phase inverter for uninterruptible-power-supply applications | |
| Singh et al. | Performance comparison of optimized controller tuning techniques for voltage stability | |
| Guha et al. | Observer‐aided resilient hybrid fractional‐order controller for frequency regulation of hybrid power system | |
| CN112865132B (zh) | 双区域互联电力系统的负荷频率控制参数的处理方法 | |
| Ławryńczuk | Nonlinear model predictive control for processes with complex dynamics: A parameterisation approach using Laguerre functions | |
| Kanchanaharuthai et al. | Fixed-time command-filtered backstepping control design for hydraulic turbine regulating systems | |
| Issa et al. | PID controller tuning parameters using meta-heuristics algorithms: comparative analysis | |
| Michalowsky et al. | Extremum control of linear systems based on output feedback | |
| Sambariya et al. | Effective PID-PSS design using Bat algorithm for SMIB power system | |
| Hsu et al. | Adaptive asymmetric fuzzy neural network controller design via network structuring adaptation | |
| Wiktorowicz | Output feedback direct adaptive fuzzy controller based on frequency-domain methods | |
| Chen et al. | Indirect model reference adaptive control for a class of linear fractional order systems | |
| Spoljarić et al. | Optimization of PID controller in AVR system by using ant lion optimizer algorithm | |
| Zhao et al. | Robust contractive economic MPC for nonlinear systems with additive disturbance | |
| Hadian et al. | A new event-based PI controller using evolutionary algorithms | |
| Wu et al. | Enhanced receding horizon optimal performance for online tuning of PID controller parameters | |
| Tajjudin et al. | Model reference input for an optimal PID tuning using PSO | |
| Rajeshwaran et al. | Hybrid Optimization Based PID Controller Design for Unstable System. |
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 |