各位做电力电子控制的朋友应该都有体会有源电力滤波器APF样机刚上电最先要确认的不是“补偿算法够不够高级”而是电网电压的相位有没有锁准。相位锁不准后面的谐波提取、dq 坐标变换、直流侧电压控制都会跟着偏移严重时甚至会造成过流。很多初学者会在 MATLAB 里搭一个理想三相电压源再配合理想 PLL 做仿真结果一到三相不平衡、电压畸变工况就暴露问题。本文针对 APF 控制系统里常用的电网同步方案之一——DSOGI-PLLDual Second Order Generalized Integrator PLL双二阶广义积分器软件锁相环做一次原理建模分析。因为内容较多本文先完成“原理部分 模型搭建思路”也就是标题里说的“上半”讲清楚它由哪些模块组成、每个模块解决什么问题、为什么比单纯的 SRF-PLL 更适合 APF 工况以及数字化实现时的参数设计与离散化注意事项。有基础的开发者可以直接跳到第 5、第 7 节看建模思路和排错内容新手建议从头往后看把坐标变换、SOGI、正负序分离这条链路理清楚。1. DSOGI-PLL 是什么为什么 APF 需要它1.1 APF 对相位同步的要求有源电力滤波器的基本任务是检测出负载谐波电流然后通过逆变器输出一个与谐波电流幅值相等、方向相反的补偿电流使电网侧的电流恢复为标准正弦波。但问题在于“检测谐波电流”并不是单纯地从采样波形里做减法和傅里叶变换。实际控制器通常会先把三相电流变换到与电网电压同步的旋转 dq 坐标系里在 dq 坐标系下分别处理基波有功、无功和谐波。此时电网电压的相位信息就变成了整个控制器的基准坐标系电网相位偏移谐波检测结果就会错误电网相位抖动指令电流会出现扰动电网相位在暂态中发生跳变很可能引发逆变器过流或直流侧电压波动。因此APF 对同步环节的要求可以概括为在电网电压不对称、含谐波、频率波动等工况下仍然能快速准确地输出基波正序电压的相位、频率并且具备一定的抗扰动能力。1.2 传统同步方案为什么不够用早期最简单的方法是硬件过零检测。通过比较器检测电网电压的过零点再用定时器计算周期从而得到相位和频率。硬件过零检测在理想电网中比较简单但在工程中却有三个明显问题第一当电网电压含有谐波时波形可能在一个周期内多次穿越零点产生错误的相位跳变。第二过零检测只能得到“过零”这个瞬间的信息半个周期内相位只能靠开环计时一旦频率波动中段相位误差会逐渐积累。第三在三相不平衡系统中某一相电压幅值很低时过零点可能不明显甚至消失导致检测失败。后来控制系统普遍改用软件锁相环Software PLL。锁相环本质上是一个闭环反馈系统它不再关注某个电压“是否过零”而是持续估计电网电压的相位和频率并通过反馈控制使估计相位不断逼近真实相位。这样既没有硬件过零检测的脆弱性又能对频率变化进行跟踪。不过根据前端滤波器结构的不同软件锁相环的性能差异很大DSOGI-PLL 就是其中一种针对“不平衡 畸变”工况设计的方案。1.3 DSOGI-PLL 的整体位置从系统角度理解DSOGI-PLL 并不是一个独立算法而是多个模块的组合三相 abc 电压经过 Clarke 变换到 αβ 静止坐标系两个 SOGI 模块分别对 vα 和 vβ 做滤波并生成两路正交信号通过正负序分离公式得到基波正序分量再经过 Park 变换进入 SRF-PLL 完成相位锁定锁相环输出的角频率反馈给两个 SOGI让 SOGI 的中心频率始终跟随电网频率。这种“DSOGI 前端 SRF-PLL 后端”的结构既保留了 SRF-PLL 锁相精度高的优点又把负序分量提前滤除避免了普通 SRF-PLL 在三相不平衡时出现二倍频锁相误差的问题。所以 DSOGI-PLL 在 APF、光伏逆变器、风电变流器等需要适应复杂电网的应用中非常常见。2. SRF-PLL从坐标变换到闭环锁相2.1 Clarke 坐标变换先来复习最基础的坐标变换。假设三相电网电压为 va、vb、vc在不考虑零序的前提下可以变换到两相静止 αβ 坐标系vα (2/3) * (va - 0.5*vb - 0.5*vc) vβ (2/3) * ((√3/2)*vb - (√3/2)*vc)这里的系数取 2/3 是为了保证变换前后幅值不变也就是“等幅值 Clarke 变换”。如果取 √(2/3)则变成等功率变换。实际工程中APF 的电流控制常使用等幅值变换因为 dq 轴指令值可以直接对应到相电压幅值概念上更直观。Clarke 变换把三相互成 120° 的交流量变成了 αβ 坐标系中两个互相垂直的交流量。在理想三相平衡电压下vα 和 vβ 是一对幅值相同、相位相差 90° 的正弦波。这个特性非常关键后面 DSOGI 之所以能工作也依赖于 αβ 坐标系中的正交关系。2.2 Park 变换与鉴相过程有了 αβ 轴分量后再经过 Park 变换进入同步旋转 dq 坐标系vd vα * cosθ vβ * sinθ vq -vα * sinθ vβ * cosθ这里的 θ 是锁相环估计出的电网电压相位。如果 θ 正好等于电网电压的真实相位且电网电压平衡对称那么 d 轴分量会等于电压幅值q 轴分量会等于 0。反过来想我们也可以把 vq 看成“相位误差指示器”当 θ 比真实相位小时vq 0当 θ 比真实相位大时vq 0当 θ 等于真实相位时vq 0。于是SRF-PLL 的闭环思路就出来了用 Park 变换得到 vq再用 PI 调节器把 vq 调节到 0。PI 调节器输出的角频率经过积分得到锁相环输出的相位 θ而 θ 又返回给 Park 变换形成闭环。这个结构在理想电网下性能很好响应速度也可以设计得比较快。2.3 SRF-PLL 面对不平衡与畸变时的问题虽然 SRF-PLL 结构简单但它有一个前提假设输入的三相电压是平衡对称的基波电压。一旦电网电压不平衡电压就会同时存在正序、负序分量。负序分量在同步旋转坐标系下会表现为 2 倍频的交流分量。为了理解这一点可以这样记忆正序分量在 dq 坐标系中是直流量负序分量相当于以相反方向旋转的空间矢量所以折算到正序 dq 坐标系中会变成 2 倍频脉动。此时 vq 不再是一个稳定的直流量而是一个叠加了 2 倍频分量的波动信号。PI 调节器很难在抑制 2 倍频波动的同时保持较快的动态响应。如果电网电压本身还含有 5 次、7 次等谐波SRF-PLL 的 vq 里还会出现更高次的脉动导致输出相位发生周期性抖动。这种相位抖动对 APF 的影响很大因为 APF 要输出的补偿电流本身就是谐波如果同步相位不干净谐波提取环节和各次指令电流的坐标变换就会互相干扰。因此我们必须在 SRF-PLL 的前面增加预处理环节把不平衡带来的负序分量和部分谐波滤除DSOGI 就是承担这个预处理的方案。3. SOGI 与正交信号发生器3.1 二阶广义积分器的传递函数SOGI 的完整名字是 Second Order Generalized Integrator中文常译为二阶广义积分器。它的基本结构是输入一个交流量 v输出两个信号 v 和 qv其中 v 是与输入同相位的滤波信号qv 是滞后 v 90° 的正交信号。理想情况下SOGI 相当于一个中心频率为 ω 的带通滤波器同时对中心频率附近的基波信号进行“提取”和“正交化”。它的典型传递函数如下D(s) v(s) / v(s) k * ω * s / (s^2 k*ω*s ω^2) Q(s) qv(s) / v(s) k * ω^2 / (s^2 k*ω*s ω^2)其中ω 是中心角频率在锁相环应用中通常来自 PLL 反馈频率估计值k 是阻尼系数也叫增益系数决定滤波器在中心频率附近的带宽D(s) 是“同相输出”的传递函数呈现带通滤波特性Q(s) 是“正交输出”的传递函数相当于在中心频率处让输出滞后输入 90°。如果输入信号频率正好等于中心频率 ω那么 D(jω) 的幅值增益约为 1相位偏移约为 0Q(jω) 的幅值增益约为 1相位偏移约为 -90°。也就是说SOGI 能从输入波形中把中心频率附近的基波“取出来”同时生成一对高质量的正交信号。3.2 同相输出与正交输出的关系从传递函数上看Q(s) 和 D(s) 之间存在一个非常经典的关系Q(s) (ω / s) * D(s)也就是说正交输出 qv 可以看作把同相输出 v 再经过一个积分环节得到的信号。在稳态正弦条件下“积分”等价于滞后 90°所以 qv 是 v 的正交信号。但要注意这里的 90° 滞后只在输入频率等于中心频率 ω 时严格成立。如果输入频率偏离中心频率qv 和 v 之间的相位关系就会偏离理想 90°同时 v 的幅值也会下降。因此在 PLL 应用中必须把 ω 做成自适应量让 SOGI 中心频率跟随电网频率变化。这也是为什么 DSOGI-PLL 一定要把 PLL 的频率输出反馈回 SOGI而不是把中心频率固定为 314 rad/s 的原因。3.3 参数 k 对动态与滤波性能的影响SOGI 里只有一个可调参数 k但它的作用却很关键。当 k 较大时SOGI 的带宽较宽动态响应更快但滤除附近谐波的能力会下降当 k 较小时SOGI 的选择性更好对谐波抑制能力强但动态响应会变慢而且中心频率附近的相位和幅值误差更容易出现。工程上比较常见的做法是把 k 取在 0.5 到 √2 之间经典取值是 √2对应比较好的阻尼特性。实际调试时不建议一上来就改 k 值而是先固定一个初始值结合整个锁相环的开环 Bode 图来观察稳定裕度。如果只单独调大 SOGI 的 k会发现锁相环对谐波的抑制变差vq 中的高频分量明显增大如果 k 调得太小电网电压发生相位跳变时PLL 跟踪速度又会变慢。所以要记住一个原则DSOGI-PLL 的整个滤波延迟由 SOGI 和 PI 调节器共同决定单独看任何一个参数都不完整。3.4 SOGI 的频率自适应思想固定中心频率的 SOGI 只适合电网频率非常稳定的场合。一旦电网频率从 50Hz 偏移到 49.5Hz固定中心频率为 50Hz 的 SOGI 就会在输出中引入幅值误差和相位误差进而影响正负序分离结果。自动适应频率的做法是把 PLL 闭环估计出的角频率 ω 直接接到 SOGI 的中心频率端口上。这样当电网频率发生变化时PLL 先检测到相位误差变化调整输出角频率SOGI 中心频率也随之改变。虽然这条反馈回路上存在延迟但由于电网频率变化通常比较缓慢工程上完全够用。更重要的是频率自适应让 SOGI 在动态过程中也能保持正交性不至于因为频率偏差产生额外的负序分量残留。4. DSOGI 与正负序分离4.1 双 SOGI 结构前面的 SOGI 只能处理一路单相交流信号。三相电网电压经过 Clarke 变换后有两路交流信号 vα 和 vβ所以要对每一路都配一个 SOGI这就是“Dual SOGI”的含义。DSOGI 的结构如下vα 进入 SOGI1输出 vα 和 qvαvβ 进入 SOGI2输出 vβ 和 qvβ两个 SOGI 共用同一个中心角频率 ωω 由后级 PLL 的反愤输出提供。DSOGI 的输出总共有四路信号vα、qvα、vβ、qvβ。其中 vα 是与 vα 同相位的基波滤波结果qvα 是滞后 vα 90° 的正交信号vβ 与 qvβ 的关系同理。这些信号会进一步送入正负序分离环节而不会直接送到 Park 变换因为只有把这些信号组合起来才能真正分离出正序分量和负序分量。4.2 用正交信号提取基波正序分量三相不平衡电压可以分解为正序分量、负序分量和零序分量。在 APF 等三相三线系统中通常不考虑零序分量因此主要任务是分离正序和负序。在 αβ 静止坐标系下设 α 轴和 β 轴信号分别经过 90° 滞后处理正序分量的计算公式可以写为vα 0.5 * (vα - qvβ) vβ 0.5 * (qvα vβ)这里的 qvβ 可以理解成 β 轴信号滞后 90° 的结果。这个公式的物理含义是把正序和负序分量在 αβ 平面上做一个“对称分量滤波”最后把负序分量抵消掉只留下正序。如果三相电压本来就是理想平衡的那么负序分量为 0DSOGI 输出 vα 约等于 vαvβ 约等于 vβ整个 DSOGI-PLL 的工作状态相当于普通的 SRF-PLL。但在三相不平衡时通过这个正序提取步骤负序分量会被有效抑制进入后级 Park 变换的只有基波正序因此 vq 中的二倍频波动会大幅下降锁相结果也就更加稳定。4.3 与直接在 dq 域加滤波器的方案对比如果不想用 DSOGI也可以保留普通 SRF-PLL然后在 vq 反馈通道加低通滤波器或陷波器来滤除二倍频分量。这种方法实现起来更简单但有一个明显问题低通滤波器本身会引入相位滞后等效地降低了锁相环的相位裕度动态响应变差陷波器又很难同时适应频率变化电网频率一旦偏离设计值陷波效果就会下降。DSOGI 方案把滤波功能放在 αβ 静止坐标系中完成而且中心频率随 PLL 输出自适应所以它既能抑制负序又能适应频率偏移。与此同时DSOGI 还天然提供了两路正交信号可以直接用于正序提取不需要额外再构造 90° 延迟。这种“滤波 正交 对称分量提取”集成在一起的特点是 DSOGI-PLL 在复杂电网工况下受欢迎的根本原因。5. DSOGI-PLL 整体建模分析5.1 整体结构把前面各个模块串联起来就得到了完整的 DSOGI-PLL 控制结构。从信号流上看整个系统可以分为以下几个部分abc 到 αβ 的 Clarke 变换模块αβ 轴双 SOGI 正交信号发生器模块正序分量计算模块αβ 到 dq 的 Park 变换模块由 q 轴误差、PI 调节器和积分器组成的锁相环模块频率反馈到两个 SOGI 的自适应回路。把这个结构画成控制框图后可以看到一条清晰的主信号路径电压采样 - Clarke - DSOGI - 正序计算 - Park - PI - 相位积分。频率 ω 则从“相位积分”前引出同时反馈给 SOGI 和用于计算。在 Simulink 里建模时最推荐的方式是把每个模块做成独立的 Subsystem这样既方便单独调试滤波参数也方便观察中间信号。比如把 SOGI 模块封装好后可以直接复制一份给 β 轴使用减少重复搭建的工作量。5.2 正序 dq 模型与相位误差经过正序分离后基波正序电压 vα 和 vβ 可以看成一对幅值恒定的正交交流量。再经过 Park 变换后vd V vq ≈ V * (θ - θ_hat)这里的 V 是正序电压幅值θ 是真实相位θ_hat 是 PLL 估计相位。当相位误差较小时sin(θ - θ_hat) 可以近似为 θ - θ_hat所以 vq 近似正比于相位误差。这正是 PI 调节器可以工作的基础锁相环把电压量变成了一个可以控制的“角度误差量”。如果希望 PI 参数不随电网电压幅值变化可以在 vq 乘以 1/V相当于对误差做了归一化这样当电压跌落或抬高时锁相环的动态特性