项目:基于MATLAB的D-STATCOM电压暂降与暂升校正仿真系统
项目介绍
本项目是一个基于MATLAB脚本编写的配电静止同步补偿器(D-STATCOM)时域仿真系统。该系统不依赖Simulink图形化模块库,而是通过纯代码方式(Script)实现了电力电子系统的离散时间步进仿真。
项目旨在模拟配电网中常见的电能质量问题,特别是电压暂降(Sage)和电压暂升(Swell),并验证D-STATCOM装置通过无功功率注入或吸收来稳定公共连接点(PCC)电压的能力。仿真内核包含了电网电源、RL负载、电压源换流器(VSC)的物理模型以及基于dq坐标系的完整双闭环控制策略。
功能特性
- 纯脚本离散仿真:脱离Simulink环境,使用MATLAB脚本直接实现离散时间前向欧拉法(Forward Euler)的数值积分求解,便于理解底层物理方程。
- 多工况故障模拟:内置故障发生器,能够按时间序列自动触发电压暂降(跌落至0.5 p.u.)和电压暂升(升至1.3 p.u.)故障。
- 双闭环控制策略:
*
外环:包含直流母线电压控制(维持电容电压稳定)和交流PCC电压幅值控制(校正电网电压波动)。
*
内环:基于dq坐标系的电流解耦控制,包含前馈解耦项,确保电流快速跟踪。
- 物理系统建模:详细构建了包含线路阻抗、L型滤波器、PCC并联电容、直流侧电容在内的完整电路微分方程模型。
- 非线性限制环节:实现了控制器输出的电流限幅保护以及SPWM调制的过调制限制。
使用方法
- 确保计算机安装有MATLAB软件(基础版本即可,无需特定工具箱)。
- 打开MATLAB环境,定位到项目所在目录。
- 清除工作区变量并运行主程序。
- 程序将执行0.3秒的离散仿真循环,计算并记录电网电压、PCC电压、注入电流及直流母线电压等状态变量。
系统要求
- 软件:MATLAB R2016b 或更高版本。
- 硬件:标准PC配置,无需高性能计算资源,因为仿真时长较短且模型经过简化优化。
---
核心仿真逻辑与实现细节
以下内容基于主程序代码的实际逻辑编写,详细解析了数学模型与控制算法的实现方式。
1. 系统参数与初始化
仿真系统设定采样频率为20kHz(步长50微秒),电网标准为380V(相间有效值)、50Hz。
- 负载模型:模拟20kW/10kVar的感性RL负载,根据功率需求自动计算等效电阻和电感。
- D-STATCOM参数:设置直流母线参考电压为700V(高于线电压峰值以确保可控性),配置了滤波电感、电阻及直流侧储能电容。
- 状态变量:初始化了所有电感电流、电容电压以及积分控制器的历史状态,使用零矩阵预分配内存以记录仿真数据。
2. 故障模拟机制
在主循环中,通过时间条件判断动态改变电网电压源的幅值比例因子,实现故障模拟:
- 0.00s - 0.08s:正常运行(1.0 p.u.)。
- 0.08s - 0.14s:深度暂降模式,电压幅值降至额定值的50%。
- 0.18s - 0.24s:暂升模式,电压幅值升至额定值的130%。
其余时间保持正常运行。该机制直接作用于三相正弦波发生器。
3. 测量与信号处理
- 锁相环(PLL):代码中采用简化的开环方式,利用仿真时间直接生成电网相位角(theta),用于dq坐标变换,假设电网频率恒定且无频率波动。
- 电压幅值检测:使用Alpha-Beta变换(Clarke变换)计算PCC电压的瞬时模值。为了滤除高频开关噪声并模拟RMS检测延迟,实现了一个一阶低通滤波器(LPF)对测量幅值进行平滑处理。
4. 控制系统架构
控制逻辑在dq旋转坐标系下执行,分为三个层级:
- 直流电压外环(d轴):通过PI控制器调节直流母线电压(Vdc),输出d轴参考电流(id_ref),负责控制有功功率交换以维持电容器能量平衡。
- 交流电压校正环(q轴):通过PI控制器比较PCC电压幅值与额定值,输出q轴参考电流(iq_ref)。该环路决定了补偿电压所需的无功电流大小。
- 电流限幅:对合成的参考电流矢量幅值进行硬限制(最大500A),防止过流,优先保证系统安全。
- 电流内环(dq轴):包含两个PI控制器,分别跟踪id和iq。关键特征是引入了交叉解耦项(w*L*iq 和 -w*L*id)以及PCC电压前馈,消除了d、q轴之间的耦合影响,输出最终的逆变器参考电压指令。
5. 调制与物理求解器
- SPWM线性化处理:将内环计算出的参考电压指令从dq域变换回abc静止域。在此过程中加入过调制保护,确保发波指令不超过直流母线电压限制的线性调制区(Vdc/sqrt(3))。
- 微分方程求解(物理引擎):
该部分是仿真的核心,利用前向欧拉法(Forward Euler)手动求解电路的微分方程,更新即时状态:
1.
电网侧方程:根据电源电压、线路阻抗和PCC电压计算电网电流变化率。
2.
注入侧方程:根据逆变器输出电压、滤波阻抗和PCC电压计算VSC注入电流变化率。
3.
负载侧方程:根据PCC电压和RL负载参数计算负载电流变化率。
4.
PCC节点电压方程:基于基尔霍夫电流定律(KCL),利用进出PCC点的电流差(电网+注入-负载)对PCC并联电容充电,计算节点电压变化率。
5.
直流侧电压方程:基于功率平衡原理,假设交流侧瞬时功率等于直流侧功率(忽略损耗),计算直流电容的电压变化。
通过上述步骤,系统在每一个仿真步长内完成了从故障生成、信号采集、控制运算到物理状态更新的完整闭环过程。