最近一直在整理高压直流输电系统的小信号稳定性问题，趁着手头这个DR-MMC串联型拓扑的项目刚结题，把阻抗建模和稳定性分析的完整思路梳理一遍。DR-MMC（Diode Rectifier–Modular Multilevel Converter）这个概念很多做柔直的人都不陌生——送端用不控二极管整流器替代可控换流器，受端用MMC，整套系统在成本和损耗上比双MMC方案有质的优势。但便宜是有代价的：二极管整流器没有主动控制能力，系统的稳定边界完全靠受端MMC和送端电源特性撑着，这时候阻抗模型就不再是“锦上添花”，而是判断整个系统会不会在某个频率点悄悄失稳的关键工具。
这篇内容我按实际操作顺序来写，从DR-MMC为什么要这么搭、阻抗建模在哪一步容易踩坑，到怎么用阻抗比判据判断稳定性、控制器参数怎么调，最后整理几个我排查过程中遇到的典型问题。适合正在做柔直控制保护或者新能源经柔直送出仿真的工程师，也适合刚接触阻抗分析法、想找个实例入门的同学。过程中会给出具体的建模流程、参数选取思路和仿真验证方法，尽量做到能直接抄作业。
1. DR-MMC系统的结构与建模难点
1.1 为什么选择DR-MMC拓扑
传统的两端柔性直流输电通常是双MMC结构，两个换流站都具备独立控制能力，运行灵活但代价高——每个MMC都需要大量子模块电容、阀控系统和辅助供电设备，一个站的造价动辄数亿元。在远距离大容量输电场景下，送端如果只是把交流电整流成直流送出去，不需要频繁调节功率方向，那用不控二极管整流器（DR）替代送端MMC就能省掉一多半设备成本。
DR-MMC串联系统的基本结构是：送端交流源经过变压器接到二极管整流器，整流器输出的直流电经过直流线路送到受端MMC，由MMC完成逆变并接入受端交流电网。这里说“串联”是指直流侧拓扑结构，整流器和MMC在直流侧串联成一条完整的直流回路，而不是像并联型多端系统那样存在多个直流节点。
这个拓扑的工程意义很直接。预计一下全生命周期成本，送端DR没有子模块电容和IGBT阀组，维护量比MMC小得多；额定工况下整流器的损耗占传输功率的比例能做到0.5%以内，而MMC的损耗通常在1%以上。对海上风电送出这类对损耗和运维成本敏感的场景，DR-MMC几乎是为它们量身定做的。
1.2 站在阻抗法角度看系统交互
但DR-MMC也有个绕不开的问题：系统的可控自由度变少了。双MMC系统里两端都能调节电压、功率，稳定性问题可以由两端协调控制解决；DR-MMC里送端完全是“被动”的，整流器输出的直流电压随交流侧电压和电流变化而波动，送端电源（比如风电场的等值机组）和整流器之间、整流器和直流线路之间、受端MMC和受端交流电网之间的耦合关系会更直接地暴露出来。
阻抗分析法的思路是把每个子系统都当作一个“端口网络”——从某个公共连接点看进去，用一个频域阻抗（或导纳）来描述它的外部特性。如果一个系统可以自然划分成“源”和“负载”两部分，那么只要这两部分的阻抗比值在满足Nyquist判据，系统就是小信号稳定的。
对DR-MMC系统来说，阻抗模型的难点集中在三个地方：
- 二极管整流器本身是强非线性器件，导通换相过程带来大量谐频分量，平均化处理需要特别小心；
- 整流器直流侧的电压是六脉动或十二脉动的，纹波会沿直流线路传到受端MMC，受端为了抑制这些纹波会调整控制带宽，控制参数和线路参数之间存在频率耦合；
- MMC的控制系统包含电流内环、功率外环、锁相环和内环本身的三相调制延迟，这些环节的截止频率相互接近时，阻抗幅频特性会出现明显“峰谷跳动”，稳定性判定容易受模型精确度影响。
2. 阻抗建模的完整流程
2.1 从时域状态空间到频域阻抗
我习惯的建模路径不是直接在频域里猜传递函数，而是先建立系统在时域的状态空间方程，再做线性化、拉普拉斯变换得到频域阻抗表达式。这条路径的好处是物理关系清楚，什么地方做了近似自己心里有数，后续要加控制环节也方便。
以受端MMC为例，建立阻抗模型分四步：
- 写出MMC在abc三相静止坐标系下的电气动态方程，包括桥臂电感、子模块电容电压和桥臂环流方程；
- 将abc坐标变换到dq旋转坐标系，得到直角坐标下的状态方程，此时基波正序分量变成直流分量，方便做稳态工作点求解；
- 在稳态工作点做小信号线性化，得到状态变量和输入变量之间的偏差传递关系；
- 对偏差方程做拉普拉斯变换，整理成端口电压与端口电流之比，即阻抗表达式。
相比直接把开关平均模型写进仿真里扫频，这种解析模型能让你拿到一个可以写成多项式的传递函数，后面做Nyquist判据和控制器参数灵敏度分析都会方便很多。
但这里有个容易掉进去的坑：dq坐标系下的阻抗是二阶矩阵，含耦合项。你要么接受一个2×2矩阵形式的阻抗，分析时处理交叉耦合；要么利用对称分量法把正序、负序阻抗分别提取出来，变成两个独立的标量阻抗。后面我会讲两种方式的取舍。
2.2 小信号线性化的关键细节
线性化这一步看着是纯数学操作，实操时几个细节直接决定模型可用不可用：
第一，工作点变量必须自洽。不能一边用额定电压算增益、一边用半载电流算稳态值，那样线性化点都不一致，出来的阻抗曲线一定是歪的。正确的做法是先让模型稳定运行到目标工况，记录各状态变量的稳态值，再在稳态解上叠加小扰动。
第二，MMC的桥臂环流和内模控制器的状态量必须包含在状态向量里。很多人做MMC阻抗建模时只保留交流侧电流状态，把环流抑制环节直接简化成比例项，结果高频段的阻抗特性完全对不上仿真。环流控制器的带宽通常设计在二倍频附近，它和直流侧阻抗特性直接相关，漏不得。
第三，锁相环（PLL）的增广模型要保留。对受端MMC来说，PLL决定了控制器的参考角度，PLL的动态会通过坐标变换耦合到阻抗表达式中。常规做法是把PLL的输入、输出作为两个增广状态变量，分别描述相角偏差对电流环给定值的扰动传递。
2.3 序阻抗与dq阻抗的选择
做稳定性分析时需要选一种阻抗表达形式。dq阻抗适合两个子系统都采用同一旋转坐标系、且坐标系可以全局共享的分析场景；序阻抗则更贴近电网故障分析的习惯，适合在公共连接点处观察正序、负序的耦合特性。
从工程效率上讲，我强烈推荐在验证阶段两者都算一遍：先用dq阻抗做控制器参数设计，因为它表达式的系数物理意义明确；再用正序阻抗做系统级稳定性校验，因为交流电网的等值阻抗在序分量下更容易获得。
具体到DR侧，二极管整流器没有控制系统，它的阻抗计算不需要增广PLL和控制器，只需要把变压器漏抗、整流器换相电抗和直流侧滤波器折算到交流侧，然后根据换相重叠角计算动态换相阻抗。注意，二极管整流器在dq坐标系下会有明显的耦合项，因为换相过程本质上是三相之间的相间短路过程，等效阻抗在不同旋转坐标之间并不完全对称。
2.4 参数扫描与验证
解析模型算完一定要用扫频仿真做对照。我常用的方法是在电磁暂态仿真软件（如PSCAD/EMTDC或MATLAB/Simulink Simscape Electrical）里搭建详细开关模型，在稳态工作点上注入一个小幅正弦扰动：
- 给电流参考值叠加一个频率从1Hz到5kHz的扰动信号，幅值取额定电流的1%—2%；
- 在每个频率点等待系统进入稳态（通常3-5个周波），采样电压和电流响应；
- 利用FFT提取该频率分量的幅值和相位，得到频率点的阻抗实测值；
- 将扫频结果与解析模型计算出的阻抗曲线叠放在同一张图上，检查幅频和相频的偏差。
如果两者偏差超过3dB或者相位偏差超过10度，先回到线性化步骤检查状态变量是否完备，而不是急着调判据。我见过不少项目在建模阶段误差很大，后面稳定性分析结论完全失真，根子都在这里。
3. 稳定性分析与判据应用
3.1 阻抗比与Nyquist判据
阻抗法判断稳定性有一个基本前提：系统能划分为源侧子系统和负载侧子系统，且它们在公共连接点处的交互频率范围内阻抗特性是线性的。对于DR-MMC直流输电系统，最常见的划分方式是：
- 从送端交流系统及DR整流器组成的“源”；
- 从直流线路及受端MMC组成的“负载”。
在直流端口处，源阻抗和负载阻抗满足关系。如果两侧阻抗在某个频点幅值相近且相位差达到180°，就可能出现负阻尼振荡。将阻抗比定义为源阻抗与负载阻抗之比，用Nyquist判据判断稳定性：只要在复平面上不包围(-1,0)点，系统稳定；如果包围，则该频点附近存在失稳风险。
这个判据直接易懂，但实际应用时有三个注意点：
- 判据使用的是“开环传递函数”形式的比值，必须先确认两侧阻抗各自本身是稳定的，再用阻抗比判断交互稳定性；
- 如果两侧都存在宽频带的多个谐振峰，需要扫描不止一个频率区间，否则可能漏掉高频段的高次谐波谐振；
- 判据给出的是稳定与否的定性结论，没有直接给出稳定裕度。工程上我还是习惯把阻抗比的幅值和相位裕度曲线也画出来，取幅值裕度大于6dB、相位裕度大于30度作为可接受的边界。
3.2 计及DR侧的约束
二极管整流器没有控制自由度，它的阻抗特性主要由主电路参数决定：交流侧等值阻抗、换相电抗、平波电抗器以及直流滤波电容。分析时要注意DR侧的阻抗并不是固定在某一频率下不变的，它随传输功率水平变化——功率越大，换相重叠角越大，动态换相阻抗中关于电流导数的项就越明显。
具体操作时，我会把DR侧的阻抗当成“不可控的负载特性”看待：它的低频段（几十Hz以下）表现为电感和电阻的串联特性，中频段（100Hz到2倍频）会出现由平波电抗器与直流滤波电容谐振带来的负阻尼区间，高频段则主要被线路分布电容和整流器阀侧杂散电容主导。
这个特性带来的直接后果是：DR-MMC系统的稳定性工作点往往是“窄带”的，不能只校验额定工况。我建议至少做三组工况的稳定性校验——10%轻载、额定工况、1.1倍过载工况——每一组都要重新计算DR侧的阻抗曲线和受端MMC的阻抗曲线，再重新做阻抗比判据。只算额定工况的模型可能在现场轻载启动时直接振荡。
3.3 控制器参数灵敏度
既然受端MMC是可调的，那它就是整个系统稳定性的“调节阀”。常见的控制器结构包括交流电压外环、电流内环、PLL和环流抑制控制器。参数调整对阻抗特性影响很大，但规律不复杂：
- 电流内环比例增益加大，中频段阻抗幅值抬升，相频曲线穿越-180度的频率点会向高频移动；
- PLL带宽增大，低频段的阻抗相位会变得更“软”，容易在几十Hz附近制造负阻尼；
- 环流抑制控制器带宽设置过高时，会在100Hz附近引入一个阻抗凹坑，如果这个凹坑恰好在DR侧等效阻抗的谐振峰附近，两者就会形成明显的交互风险。
工程中我的习惯是先固定PLL带宽（通常取10-20Hz），再整定电流内环（截止频率取150-300Hz），最后调环流抑制（截止频率取50-100Hz），每调完一个环节就重新看一次阻抗比曲线。控制器参数设计的核心目标不是让单一换流器的阻抗很好看，而是让阻抗比在全频段都离(-1,0)点足够远。
4. 仿真与实测对照
4.1 仿真模型搭建
仿真模型的精度等级决定了稳定性分析的结论可信度。我不建议直接把一个等值开关模型用于全部稳定性分析，太高自由度的模型反而难以定位问题。
我的做法是分两级搭建：
第一级是详细电磁暂态模型，用于验证解析阻抗模型以及观察失稳时的真实波形。MMC每个桥臂用工程实际子模块数量（例如200个半桥子模块），阀控采用最近电平逼近调制，控制器按实际参数写。这个模型计算量大，但能揭示高频谐振等细节。
第二级是平均值模型，用于系统级稳定性统计和参数扫描。MMC用受控电压源代替，DR用理想整流桥加换相阻抗等效，直流线路用集中参数π型等效。平均模型在几百Hz以下的精度完全够用，计算仿真速度可以提升几十倍，适合做大批量的工况扫描。
搭建过程中要注意整流变压器接线组别和换相电抗的处理，DR的换相过程严重依赖变压器漏感参数，这里如果直接用理想变压器，扫频结果会在特征谐波处失真。
4.2 关键波形与结果
以一个典型的DR-MMC串联系统为例，送端交流电压等级220kV，受端MMC直流额定电压±320kV，传输功率1000MW，直流线路长度200km。在额定工况下，解析阻抗模型和详细模型的扫频结果曲线能很好地吻合，特别是在0.1Hz到1kHz区间，幅值偏差基本在1dB以内。
稳定性校验场景下，我特意把受端MMC的PLL带宽从15Hz逐步提高到35Hz，观察阻抗比曲线变化。PLL带宽为15Hz时，阻抗比在30Hz附近虽然有峰值隆起，但距离Nyquist临界点还有10dB以上的裕度，系统稳定；当PLL带宽提高到35Hz，阻抗比的峰值显著抬高，同时相位裕度跌到20度以下。再结合时域仿真，此时直流电流出现30Hz附近的振荡分量，振幅经过几十个周波缓慢增长，虽然不会立刻导致换流器闭锁，但已经逼近保护定值——这就是典型的稳定裕度不足表现。
这个对照实验说明，阻抗比曲线上的“隆起峰值”并非只是一个理论量，它直接对应时域里可以观察到的振荡模态。反过来，如果你在实测数据里看到直流电流存在某种频率的振荡，也可以通过阻抗比曲线快速定位是哪一侧、哪个控制环节在“制造”负阻尼。
4.3 控制器参数修正
针对上述失稳风险，我在项目中采用了两步修正：
第一步，把PLL带宽从35Hz降回20Hz，同时调整PLL阻尼比以减小30Hz附近的相位滞后。这一步直接改善了阻抗比的低频段特性。
第二步，修改电流内环的交叉解耦设置，把dq轴耦合项前馈系数从理论值稍微放宽到实际值的95%，适度解耦后，中频段的阻抗凹坑变浅，和DR侧谐振峰的交互被削弱。
修正后重新扫频，阻抗比峰值下降约6dB，相位裕度恢复到35度左右；时域中30Hz振荡分量不再增长，直流电流平稳，系统通过小信号稳定校验。
很多文章到这里就结束了，但实际项目中还有第三步：把修正后的参数在详细电磁暂态模型里做一次三相短路故障和功率阶跃测试，确保控制器参数调整没有牺牲暂态响应能力。实测下来，PLL带宽降下来之后暂态功角恢复略慢，但换流器没有出现直流过压或电流越限，满足工程验收标准。
5. 常见问题排查与工程建议
5.1 阻抗拟合偏差超限
解析阻抗模型与扫频仿真不匹配时，先检查线性化点的稳态值是否由同一个运行工况计算得来——我踩过一次坑，功率从1000MW改到500MW后只更新了电流工作点，没有更新电容电压平均值和桥臂环流稳态值，结果在50Hz附近阻抗相位差了15度。
另一个高频打包因素来自子模块电容的等效串联电阻损耗。解析模型里如果忽略这部分损耗，高频段的阻抗实部会偏低，而实部恰恰决定了阻尼特性。建议在频域表达式里给每个子模块电容串联一个等效电阻，阻值取厂家数据或经验值。
5.2 谐振点偏移
扫频仿真得到的阻抗前段谐振峰的位置和解析模型对不上，先检查直流线路模型。400km以上的长线路阶段需要按分布参数建模或至少分成两个π型节，集中参数π模型会把线路谐振频率计算偏低，这对稳定性判据的结论影响是方向性的。
如果谐振点整体偏移不多但幅值偏高，排查平波电抗器和直流滤波器之间的杂散电容。实际工程现场的杂散参数无法忽略，我通常用扫频仪的实测结果反向校核仿真里的杂散参数设置。
5.3 工程建议
根据这几个项目经验，给后面做类似系统的人几条实用的建议：
- 先把阻抗模型做成脚本化，参数变化时能够一键重新生成阻抗曲线和稳定性报告，不然后续做工况批量校验时会被重复劳动拖死；
- 稳定性校验报告里一定要包含阻抗比幅频、相频曲线和Nyquist图三者，只给一个判据结论无法应对评审专家的追问；
- 控制器参数优化时用自动化搜索比手工调试更快，我用过粒子群算法对PLL带宽、电流内环截止频率做多目标优化，同时约束阻抗比裕度和暂态响应指标，收敛速度令人满意；
- DR-MMC系统的稳定性问题往往在高频段（几百Hz级）出现，这部分问题在平均值模型里可能被掩盖，务必在详细模型里做最终验证。
回头再看这个项目，阻抗建模和稳定性分析并不是一个“算完交差”的任务，它的价值在于把一连串看似孤立的工程决策——拓扑选择、控制器参数、直流线路设计——串在了一条明确的分析链条上：建立模型、验证模型、分析交互、优化参数、回归验证。这套流程走通之后，做别的拓扑（比如全半桥混合结构或者双极配置）的稳定性分析也只是换模型、不换思路的事。后面我打算把DR侧稳态工作点的迭代计算部分做一个更稳定的求解器，现在遇到轻载工况下工作点迭代偶尔不收敛的问题，等有新进展再来更新。
到O(2ⁿ)的工程直觉)
