【WLSM、FDM状态估计】电力系统状态估计研究附Matlab代码

WLSM与FDM电力系统状态估计研究

✅作者简介:热爱科研的Matlab仿真开发者,擅长数据处理、建模仿真、程序设计、完整代码获取、论文复现及科研仿真。

🍎 往期回顾关注个人主页:Matlab科研工作室

🍊个人信条:格物致知,完整Matlab代码及仿真咨询内容私信。

🔥 内容介绍

电力系统状态估计(State Estimation, SE)是现代电力系统运行、控制和规划的基础。它利用冗余的量测数据,估计出电力系统在某一时刻的状态变量(通常指母线电压幅值和相角),为电力系统的安全分析、稳定控制、经济运行和风险评估提供准确可靠的数据支撑。随着智能电网的快速发展,电力系统的规模日益庞大,拓扑结构日趋复杂,分布式电源接入比例不断提高,这对状态估计的准确性、速度和鲁棒性提出了更高的要求。因此,对电力系统状态估计方法进行深入研究具有重要的理论意义和实际应用价值。

本文主要探讨了基于加权最小二乘法(Weighted Least Squares, WLSM)和快速解耦法(Fast Decoupled Method, FDM)的电力系统状态估计方法。这两种方法是目前应用最为广泛、最为成熟的状态估计方法,在电力工程实践中发挥着重要作用。

一、加权最小二乘法(WLSM)

WLSM是电力系统状态估计最常用的方法之一。其核心思想是通过最小化量测残差的加权平方和来估计状态变量。假设电力系统拥有m个量测值(如功率注入、支路功率潮流、电压幅值),n个状态变量(通常是母线电压幅值和相角)。量测方程可以表示为:

z = h(x) + e

其中,z是m维量测向量,x是n维状态向量,**h(x)**是m维量测函数向量,e是m维量测误差向量,通常假设服从均值为零、方差为R的正态分布。

WLSM的目标函数是:

J(x) = [z - h(x)]<sup>T</sup>R<sup>-1</sup>[z - h(x)]

其中,R是m×m的量测误差协方差矩阵,通常是一个对角矩阵,其对角元素是各个量测误差的方差。R<sup>-1</sup>称为权重矩阵,反映了各个量测的精度。精度越高的量测,权重越大,其对状态估计结果的影响也越大。

为了最小化目标函数J(x),需要求解以下方程:

g(x) = ∇J(x) = H(x)<sup>T</sup>R<sup>-1</sup>[z - h(x)] = 0

其中,H(x)是m×n的雅可比矩阵,其元素是量测函数对状态变量的偏导数。

由于量测函数h(x)通常是非线性的,因此需要采用迭代方法求解上述非线性方程组。常用的迭代方法是牛顿-拉夫逊法:

x<sub>k+1</sub> = x<sub>k</sub> - G(x<sub>k</sub>)<sup>-1</sup>g(x<sub>k</sub>)

其中,G(x)是H(x)<sup>T</sup>R<sup>-1</sup>H(x),被称为增益矩阵。

WLSM的优点在于:

  • 理论基础扎实,能够提供统计意义上的最优估计。

  • 可以处理各种类型的量测数据,例如功率注入、支路潮流、电压幅值等。

  • 能够进行量测坏数据的辨识,提高状态估计的可靠性。

WLSM的缺点在于:

  • 计算量大,需要求解大型稀疏线性方程组,对计算机的性能要求较高。

  • 对初始值敏感,如果初始值选择不当,可能会导致迭代不收敛。

  • 收敛速度较慢,在大规模电力系统中尤其明显。

二、快速解耦法(FDM)

为了克服WLSM计算量大的缺点,学者们提出了快速解耦法(FDM)。FDM是基于电力系统潮流计算的经验观察,即P-θ和Q-V具有较强的关联性,而P-V和Q-θ的关联性较弱。因此,可以将量测方程进行解耦,从而简化状态估计的计算过程。

FDM将WLSM的增益矩阵近似分解为两个对角块矩阵:

[H(x)<sup>T</sup>R<sup>-1</sup>H(x)] ≈ [B' 0; 0 B'']

其中,B'对应于P-θ的关联,B''对应于Q-V的关联。

通过这种近似,可以将状态估计的迭代过程分解为两个独立的迭代过程:

Δθ = (B')<sup>-1</sup>H<sub>P</sub><sup>T</sup>R<sub>P</sub><sup>-1</sup>(z<sub>P</sub> - h<sub>P</sub>(x))

ΔV = (B'')<sup>-1</sup>H<sub>Q</sub><sup>T</sup>R<sub>Q</sub><sup>-1</sup>(z<sub>Q</sub> - h<sub>Q</sub>(x))

其中,Δθ和ΔV分别是相角和电压幅值的修正量,H<sub>P</sub>和H<sub>Q</sub>分别是功率注入和电压幅值对状态变量的雅可比矩阵,z<sub>P</sub>和z<sub>Q</sub>分别是功率注入和电压幅值的量测向量,R<sub>P</sub>和R<sub>Q</sub>分别是功率注入和电压幅值的量测误差协方差矩阵。

FDM的优点在于:

  • 计算速度快,大大减少了计算量,适用于大规模电力系统。

  • 对初始值不敏感,收敛性较好。

  • 算法简单,易于实现。

FDM的缺点在于:

  • 精度较低,由于进行了近似处理,导致状态估计的精度有所降低。

  • 适用范围有限,只适用于径向性较强的电力系统。

  • 对量测数据的要求较高,如果量测数据质量较差,可能会导致状态估计结果不稳定。

三、WLSM和FDM的比较与改进

WLSM和FDM各有优缺点。WLSM精度高,适用范围广,但计算量大;FDM计算速度快,但精度较低,适用范围有限。为了克服各自的缺点,学者们提出了许多改进方法。

  • WLSM的改进:

     可以采用稀疏矩阵技术来提高计算效率,例如采用LU分解、Cholesky分解等方法。此外,还可以采用并行计算技术,将计算任务分配到多个处理器上并行执行,从而提高计算速度。还可以研究自适应步长算法,动态调整迭代步长,提高收敛速度。

  • FDM的改进:

     可以采用补偿技术来提高精度,例如采用一次或多次补偿迭代,对状态估计结果进行修正。此外,还可以采用自适应解耦算法,根据电力系统的运行状态动态调整解耦矩阵,从而提高状态估计的精度。还可以结合智能算法,例如粒子群算法、遗传算法等,对解耦矩阵进行优化,提高状态估计的性能。

四、应用展望

电力系统状态估计是智能电网的关键技术之一,在未来的电力系统中将发挥越来越重要的作用。随着智能传感技术、通信技术和计算技术的快速发展,状态估计的应用领域将不断拓展。

  • 实时监控与预警:

     状态估计可以实时监控电力系统的运行状态,及时发现潜在的安全隐患,为电力系统的安全稳定运行提供保障。

  • 智能决策与控制:

     状态估计可以为智能决策和控制提供准确可靠的数据支撑,例如为电压控制、无功优化、故障诊断和恢复等提供参考。

  • 需求响应与能源管理:

     状态估计可以为需求响应和能源管理提供实时数据,例如为负荷预测、分布式电源调度和储能系统控制等提供依据。

  • 配电网状态估计:

     随着分布式电源的大量接入,配电网的拓扑结构日益复杂,对配电网状态估计提出了更高的要求。需要研究更加精确、更加快速的配电网状态估计方法。

五、结论

本文对基于加权最小二乘法(WLSM)和快速解耦法(FDM)的电力系统状态估计方法进行了研究。WLSM精度高,适用范围广,但计算量大;FDM计算速度快,但精度较低,适用范围有限。为了克服各自的缺点,学者们提出了许多改进方法。随着智能电网的快速发展,电力系统状态估计的应用领域将不断拓展,对状态估计的准确性、速度和鲁棒性提出了更高的要求。未来的研究方向包括:研究更加精确、更加快速的配电网状态估计方法;研究基于人工智能的状态估计方法;研究面向分布式电源的状态估计方法;研究面向智能电网的状态估计方法。相信随着技术的不断进步,电力系统状态估计将为电力系统的安全、稳定、经济运行提供更加可靠的保障。

⛳️ 运行结果

🔗 参考文献

[1] 牛胜锁,刘颖,梁志瑞,等.基于广域测量和抗差最小二乘法的电力系统谐波状态估计[J].电力系统保护与控制, 2012, 40(8):5.DOI:10.3969/j.issn.1674-3415.2012.08.002.

[2] 刘怡.GPS功角测量在电力系统状态估计中的应用[D].华北电力大学(北京) 华北电力大学,2002.DOI:10.7666/d.Y554714.

[3] 汪瑶.基于相量量测的电力系统谐波状态估计算法研究[D].重庆大学,2009.DOI:10.7666/d.y1665245.

📣 部分代码

🎈 部分理论引用网络文献,若有侵权联系博主删除

 👇 关注我领取海量matlab电子书和数学建模资料 

🏆团队擅长辅导定制多种科研领域MATLAB仿真,助力科研梦:

🌈 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化、背包问题、 风电场布局、时隙分配优化、 最佳分布式发电单元分配、多阶段管道维修、 工厂-中心-需求点三级选址问题、 应急生活物质配送中心选址、 基站选址、 道路灯柱布置、 枢纽节点部署、 输电线路台风监测装置、 集装箱调度、 机组优化、 投资优化组合、云服务器组合优化、 天线线性阵列分布优化、CVRP问题、VRPPD问题、多中心VRP问题、多层网络的VRP问题、多中心多车型的VRP问题、 动态VRP问题、双层车辆路径规划(2E-VRP)、充电车辆路径规划(EVRP)、油电混合车辆路径规划、混合流水车间问题、 订单拆分调度问题、 公交车的调度排班优化问题、航班摆渡车辆调度问题、选址路径规划问题、港口调度、港口岸桥调度、停机位分配、机场航班调度、泄漏源定位
🌈 机器学习和深度学习时序、回归、分类、聚类和降维

2.1 bp时序、回归预测和分类

2.2 ENS声神经网络时序、回归预测和分类

2.3 SVM/CNN-SVM/LSSVM/RVM支持向量机系列时序、回归预测和分类

2.4 CNN|TCN|GCN卷积神经网络系列时序、回归预测和分类

2.5 ELM/KELM/RELM/DELM极限学习机系列时序、回归预测和分类
2.6 GRU/Bi-GRU/CNN-GRU/CNN-BiGRU门控神经网络时序、回归预测和分类

2.7 ELMAN递归神经网络时序、回归\预测和分类

2.8 LSTM/BiLSTM/CNN-LSTM/CNN-BiLSTM/长短记忆神经网络系列时序、回归预测和分类

2.9 RBF径向基神经网络时序、回归预测和分类

2.10 DBN深度置信网络时序、回归预测和分类
2.11 FNN模糊神经网络时序、回归预测
2.12 RF随机森林时序、回归预测和分类
2.13 BLS宽度学习时序、回归预测和分类
2.14 PNN脉冲神经网络分类
2.15 模糊小波神经网络预测和分类
2.16 时序、回归预测和分类
2.17 时序、回归预测预测和分类
2.18 XGBOOST集成学习时序、回归预测预测和分类
2.19 Transform各类组合时序、回归预测预测和分类
方向涵盖风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、用电量预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断
🌈图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知
🌈 路径规划方面
旅行商问题(TSP)、车辆路径问题(VRP、MVRP、CVRP、VRPTW等)、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、 充电车辆路径规划(EVRP)、 双层车辆路径规划(2E-VRP)、 油电混合车辆路径规划、 船舶航迹规划、 全路径规划规划、 仓储巡逻
🌈 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配、无人机安全通信轨迹在线优化、车辆协同无人机路径规划
🌈 通信方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化、水声通信、通信上传下载分配
🌈 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化、心电信号、DOA估计、编码译码、变分模态分解、管道泄漏、滤波器、数字信号处理+传输+分析+去噪、数字信号调制、误码率、信号估计、DTMF、信号检测
🌈电力系统方面
微电网优化、无功优化、配电网重构、储能配置、有序充电、MPPT优化、家庭用电
🌈 元胞自动机方面
交通流 人群疏散 病毒扩散 晶体生长 金属腐蚀
🌈 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合、SOC估计、阵列优化、NLOS识别
🌈 车间调度
零等待流水车间调度问题NWFSP 、 置换流水车间调度问题PFSP、 混合流水车间调度问题HFSP 、零空闲流水车间调度问题NIFSP、分布式置换流水车间调度问题 DPFSP、阻塞流水车间调度问题BFSP

👇

DATA: FILE IS D:\360MoveData\Users\asus\Desktop\纵向(id).dat; VARIABLE: NAMES ARE id school grade name number sex birth only cadre father mother A11 A12 A13 A18 A111 A112 A114 A14 A15 A16 A17 A19 A110 A113 B11 B12 B13 B14 B15 B16 B17 B18 B19 B110 B111 B112 B113 B114 B115 B116 B117 B118 C11 C13 C14 C16 C18 C19 C12 C15 C17 C110 C111 D11 D12 D13 D14 D15 D16 A21 A22 A23 A28 A211 A212 A214 A24 A25 A26 A27 A29 A210 A213 B21 B22 B23 B24 B25 B26 B27 B28 B29 B210 B211 B212 B213 B214 B215 B216 B217 B218 C21 C23 C24 C26 C28 C29 C22 C25 C27 C210 C211 D21 D22 D23 D24 D25 D26 A31 A32 A33 A38 A311 A312 A314 A34 A35 A36 A37 A39 A310 A313 B31 B32 B33 B34 B35 B36 B37 B38 B39 B310 B311 B312 B313 B314 B315 B316 B317 B318 C31 C33 C34 C36 C38 C39 C32 C35 C37 C310 C311 D31 D32 D33 D34 D35 D36 V159 V160 V161 V162 V163 V164 V165 V166 V167 V168 V169 V170 x1 x2 x3 y1 y2 y3 m1 m2 m3 w1 w2 w3; USEVARIABLES = x1 x2 x3 y1 y2 y3 m1; ANALYSIS: ESTIMATOR = MLR; MODEL: ! Step 1: 建立 x 的线性增长模型 ix sx | x1@-1 x2@0 x3@1; ! Step 2: 建立 y 的线性增长模型 iy sy | y1@-1 y2@0 y3@1; ! Step 3: 设定中介路径:x 的发展趋势 → m1 → y 的发展趋势 m1 ON ix sx; ! x 的初始水平和变化速度都可能影响 m1 sy ON m1 sx; ! y 的发展趋势受 m1 和 x 斜率共同影响 iy ON m1 ix; ! 允许潜变量协方差 ix WITH sx; iy WITH sy; sx WITH sy; MODEL INDIRECT: sy IND m1; ! 检验 m1 在 sx→sy 中的间接效应 iy IND m1; OUTPUT: STAND SAMPSTAT TECH1 TECH8 CINTERVAL(BCBOOTSTRAP); PLOT: TYPE = PLOT3; SERIES = x1-x3 (*) y1-y3 (*); (在mplus跑以上代码后变这些怎么办*** WARNING in OUTPUT command BOOTSTRAP and BCBOOTSTRAP confidence intervals require the specification of BOOTSTRAP in the ANALYSIS command. Request for CINTERVAL is ignored. *** WARNING in OUTPUT command TECH8 option is available only with analysis types MIXTURE, RANDOM, or TWOLEVEL with estimators ML, MLF, or MLR or ALGORITHM=INTEGRATION. Request for TECH8 is ignored. 2 WARNING(S) FOUND IN THE INPUT INSTRUCTIONS)
最新发布
11-11
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值