多区互联电力系统负荷频率的约束GPC优化控制方法
专利汇可以提供多区互联电力系统负荷频率的约束GPC优化控制方法专利检索，专利查询，专利分析的服务。并且本 发明 公开一种多区互联电 力 系统负荷 频率 的约束GPC优化控制方法，通过小 信号 机理分析建模方法和Z变换建立多区互联电力系统负荷频率的受控自回归积分滑动平均模型，采用约束广义预测控制作为控制方法，将系统预测输出与参考轨迹之间误差的平方值、 控制器 输出增量的平方值和控制输出的平方值这三者的加权值作为优化目标函数，并设计一种高效的自适应进化优化求解器实现约束GPC的在线滚动优化，实时更新系统各个变量从而实现负荷频率优化控制。采用本发明可实现多区互联电力系统负荷频率优化控制效果，相比 现有技术 具有更优良的稳态性能和动态性能以及更健壮的鲁棒性。，下面是多区互联电力系统负荷频率的约束GPC优化控制方法专利的具体信息内容。
1.一种多区互联电力系统负荷频率的约束GPC优化控制方法，其特征在于，该方法包括以下步骤：
(1)通过小信号机理分析建模方法建立多区互联电力系统负荷频率的状态空间模型：
Y＝CX (2)
其中，状态向量X＝[X1,X2,…,Xn]，系统输入U＝[U1,U2,…,Un]，负荷扰动W＝[W1,W2,…,Wn]，系统输出Y＝[Y1,Y2,…,Yn]，n表示分区的数量，第i个状态向量为Xi＝[Δfi,ΔPmi,ΔPvi,ΔPti]T，Δfi表示第i个分区的频率偏差，ΔPmi表示第i个分区的发电功率偏差，ΔPvi表示第i个分区的调节阀位置偏差，ΔPti表示第i个分区的联络线功率偏差；第i个系统输入为Ui＝ΔPri，ΔPri表示第i个分区的负荷参考设定值；第i个负荷扰动为Wi＝ΔPLi，ΔPLi表示第i个分区的负荷扰动偏差；第i个系统输出为Yi＝ACEi＝β1iΔfi+β2iΔPti，ACEi表示第i个区域控制偏差，β1i和β2i分别表示第i个分区频率偏差和联络线功率偏差对应的权重系数；
A、B、C、E分别表示多区互联电力系统的状态矩阵、输入矩阵、输出矩阵和扰动矩阵，具体描述如下：
其中，Di表示第i个分区发电机的负荷频率变化量比值系数，Mai表示第i个分区发电机组的总惯性时间系数，TGi表示第i个分区的调速器时间常数，TCHi表示第i个分区的汽轮机时间常数，Tij表示第i个分区和第j个分区之间联络线功率同步系数，Rfi表示第i个分区的发电机组的下垂系数；
(2)通过Z变换建立多区互联电力系统负荷频率的受控自回归积分滑动平均模型：
其中y(k)表示系统在k时刻的输出信号，u(k-1)表示在k-1时刻的控制输出信号，ω(k-
1)表示在k-1时刻的负荷扰动信号，ε(k)表示均值为0且方差为δ2的n维白噪声，Δd表示差分因子，Δd＝1-z-1； 表示关于z-1的n×n维对角矩阵， 表示关于z-1的n×m维矩阵，表示关于z-1的n×p维矩阵，其中 的第i个对角元素 的第i行第
j列元素 的第i行第j列元素 的表达式如下：
其中， 表示 各阶次的系数， 表示 各阶次的系数，
表示 各阶次的系数，nai和nbij分别表示 和 的最高阶次；
(3)设置参数值：预测时域长度N，控制域长度Nu，权重系数矩阵Q、R1和R2,柔化系数矩阵αj，种群规模p，变异系数b，最大迭代次数Imax；
(4)读取当前时刻k的历史信息，k＝1,2,…,tmax，包括系统输出信号y(k-1),…,y(k-na)和控制器输出信号u(k-1),…,u(k-nb)，tmax表示系统运行时窗的最大值；
(5)随机产生满足约束条件的实数编码的种群P＝{S1,S2,…,Sp}，其中第i个体Si表示待优化的控制增量序列{Δu(k),Δu(k+1),…,Δu(k+Nu-1)}，具体产生过程如下：
Si＝Δumin+ri(Δumax-Δumin),i＝1,2,...,p
其中，Δumin和Δumax分别表示控制增量序列的下限和上限，ri表示一组在0和1之间产生的随机数；
(6)按照公式(4)～(10)所示的约束GPC优化目标对种群P进行适应度函数评价，获得当前时刻k种群P的最好适应度函数Fbest＝min{F(P,k)}，将对应的个体设置为最好个体Sbest，其中F(P,k)表示当前时刻k种群P的适应度集合；
s.t.Δumin≤Δu(k+i-1)≤Δumax,i＝1,2,...,Nu (5)
umin≤u(k+i-1)≤umax,i＝1,2,...,Nu (6)
ymin≤y(k+j)≤ymax,j＝1,2,...,N (7)
Δu(k+i-1)＝u(k+i-1)-u(k+i-2) (8)
ω(k+j)＝αjω(k+j-1)+(1-αj)yr(k),ω(k)＝yp(k),j＝1,2,...,N (9)y(k+j)＝Gj(z-1)Δu(k+j-1)+Hj(z-1)Δu(k-1)+Fj(z-1)y(k),j＝1,2,...,N (10)其中，y(k+j)表示系统的预测输出信号，ω(k+j)表示参考轨迹，yr(k)表示设定输入信号，yp(k)为系统实际输出信号，Q、R1、R2表示权重系数矩阵，umin和umax分别表示控制序列的下限和上限，ymin和ymax分别表示预测输出信号的下限和上限，公式(10)所示的y(k+j)的计算方法为求取如下丢番图方程(11)～(12)获得：
其中，l0,l1,…,lj-1表示多项式Lj(z-1)各阶次的系数， 表示多项式Fj(z-1)各阶次的系数，g0,g1,…,gj-1表示多项式Gj(z-1)各阶次的系数， 表示多项式Hj(z-1)各阶次的系数，na和nb分别表示Fj(z-1)和Hj(z-1)的最高阶次；
(7)按照公式(13)～(14)所示的实数变异操作因子产生新的种群Pm＝{Sm1,Sm2,…,Smp}，并无条件地接受P＝Pm；
其中r、r1是[0,1]范围内产生的随机数，t表示当前迭代次数，b为变异操作的可调系数，Imax为用户设定的最大迭代次数，γ(t,α)表示关于t和α的函数，α的具体描述如下：当r<0.5时，α＝Δumax-Si，否则α＝Si- Δ umin ；
(8)重复步骤(6)～(7)直到优化求解器的迭代次数达到Imax，并保存Sbest＝{Δu*(k+j),j＝0,1,…,Nu-1}；
(9)计算在最优控制增量Δu*(k)作用下的当前时刻k对应的系统输出y*(k)和u*(k)，并设置为Δu(k)＝Δu*(k)，y(k)＝y*(k)，u(k)＝u*(k)；
(10)重复步骤(4)～(9)直到系统运行时窗达到tmax，输出多区互联电力系统最优系统输出曲线、最优控制增量信号曲线和最优控制信号曲线。
多区互联电力系统负荷频率的约束GPC优化控制方法
技术领域
背景技术
发明内容
第i行第j列元素 的第i行第j列元素 的表达式如下：
阶次。(3)设置参数值：预测时域长度N，控制域长度Nu，权重系数矩阵Q、R1和R2,柔化系数矩阵αj，种群规模p，变异系数b，最大迭代次数Imax；
具体实施方式
第i行第j列元素 的第i行第j列元素 的表达式如下：
阶次。(3)设置参数值：预测时域长度N＝4，控制域长度Nu＝8，权重系数矩阵Q＝2×I4×4、R1j
＝3×I8×8和R2＝1.5×I8×8，I表示单位矩阵，柔化系数矩阵α＝0.3×I4×4，种群规模p＝20，变异系数b＝20，最大迭代次数Imax＝30；
多项式Hj(z-1)各阶次的系数，na和nb分别表示Fj(z-1)和Hj(z-1)的最高阶次。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 用于运行由燃料运行的车辆加热器的方法 | 2020-06-29 | 0 |
| 一种获取角度测量误差曲线的方法 | 2020-09-28 | 2 |
| 电梯油箱油量远程监控系统及方法 | 2020-07-13 | 0 |
| 一种预应力混凝土连续刚构桥温度作用有限元分析方法 | 2021-01-21 | 0 |
| 转向系统中方向盘扭振的补偿方法和装置 | 2021-08-07 | 0 |
| 一种基于模糊控制的虚拟同步机虚拟惯性自适应调节方法及频率控制方法 | 2021-01-30 | 1 |
| Temperature compensated timing signal generator | 2021-12-03 | 1 |
| GNSS signal processing with clock difference modeling | 2023-06-02 | 1 |
| Apparatus and method for changing frequency deviation | 2021-06-10 | 1 |
| ELECTROSURGICAL ULTRASONIC VESSEL SEALING AND DISSECTING SYSTEM | 2021-06-29 | 1 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
