3. 配电网潮流计算的数学基础
配电网潮流计算的核心任务，是在给定网络拓扑、支路参数和节点注入功率的条件下，求出各节点的电压幅值与相角，以及各支路的功率分布。要理解潮流计算为什么能算、怎么算、为什么配电网算起来更困难，必须回到四个数学基础：节点导纳矩阵、功率方程、雅可比矩阵和牛顿-拉夫逊法。本节按“模型—方程—求解—适用性”的逻辑展开。
一、节点导纳矩阵：网络的线性代数描述
配电网由节点（母线）和支路（线路、变压器、开关）构成。把每条支路的阻抗折算为导纳后，可以按基尔霍夫电流定律写出每个节点的注入电流与各节点电压的关系。对含 n 个节点的网络，写成矩阵形式为：
I = Y · U
其中 I 是节点注入电流列向量，U 是节点电压列向量，Y 即节点导纳矩阵。Y 的元素有明确物理含义：对角元 Yii 是连接到节点 i 的所有支路导纳之和，称为自导纳；非对角元 Yij 是节点 i 与 j 之间支路导纳的负值，称为互导纳。若 i、j 之间无直接支路，则 Yij=0。
配电网的 Y 矩阵有两个显著特点。第一，高度稀疏：每条馈线只与相邻节点相连，一个节点平均只与两三个节点有连接，因此每行非零元很少。第二，结构上接近辐射状：正常开环运行时，网络是一棵以变电站母线为根的树，Y 矩阵经适当排序后可呈现近似三角结构。这两点决定了后续求解算法可以利用稀疏性大幅降低计算量。
二、功率方程：从线性到非线性
实际配电网中，负荷和分布式电源给出的通常是功率而非电流。节点注入功率与电压、电流的关系为 Si=UiIi*，代入 I=YU 并展开，得到节点 i 的复功率方程：
S_i = U_i · Σ_j (Y_ij · U_j)*
把 Si=Pi+jQi、Ui=Vi∠θi、Yij=Gij+jBij 代入，分离实部与虚部，得到每个节点的有功、无功功率方程：
P_i = V_i · Σ_j V_j (G_ij·cosθ_ij + B_ij·sinθ_ij)
Q_i = V_i · Σ_j V_j (G_ij·sinθ_ij − B_ij·cosθ_ij)
其中 θij=θi−θj。这两个方程是潮流计算的基本方程，其关键特征是电压以乘积形式出现、相角以三角函数形式出现，因此对未知量 V 和 θ 而言是非线性的。一个 n 节点网络有 2n 个方程，但每个节点的 P、Q、V、θ 四个量中已知两个、未知两个，故实际待求量为 2n 个，方程数与未知数相等。
按已知量不同，节点分为三类：PQ 节点（已知 P、Q，求 V、θ），是负荷节点和大多数分布式电源并网点的常见类型；PV 节点（已知 P、V，求 Q、θ），对应具备无功电压调节能力的机组；平衡节点（已知 V、θ，求 P、Q），一般取变电站高压侧母线，承担全网网损。配电网中平衡节点通常只有一个，PQ 节点占绝大多数，PV 节点较少甚至没有。
三、雅可比矩阵与牛顿-拉夫逊法
非线性方程组无法直接求解，牛顿-拉夫逊法的思路是：在当前估计点把非线性方程线性化，用切线近似曲线，解线性方程得到修正量，反复迭代逼近真解。设待求向量 x=[θ; V]，不平衡量向量 Δf=[ΔP; ΔQ]，则线性化方程为：
J · Δx = Δf
J 即雅可比矩阵，其元素是功率方程对电压幅值和相角的偏导数，按四块组织：
J = [ ∂P/∂θ ∂P/∂V ]
[ ∂Q/∂θ ∂Q/∂V ]
对角块 ∂P/∂θ、∂Q/∂V 反映同一节点自身功率对自身相角、幅值的敏感度；非对角块 ∂P/∂V、∂Q/∂θ 反映节点间的耦合。雅可比矩阵与节点导纳矩阵同维、同稀疏结构，但它是随迭代变化的——每轮迭代电压更新后，J 的元素都要重新计算。这正是牛顿法计算量大的根源。
完整迭代流程为：给定电压初值（通常取平启动 V=1.0、θ=0）→ 计算 ΔP、ΔQ → 判断最大不平衡量是否小于收敛精度 ε（如 1e-5）→ 若未收敛则形成 J、解线性方程得 Δx → 更新 θ、V → 回到第二步。牛顿法具有平方收敛特性，初值合理时通常 3～5 次迭代即可收敛。
四、牛顿-拉夫逊法在配电网中的适用性分析
牛顿-拉夫逊法在输电网中是标准算法，但直接搬到配电网会遇到几个现实问题，需要具体分析。
第一，收敛性受 R/X 比影响。输电网线路电抗远大于电阻，雅可比矩阵对角优势明显，牛顿法收敛稳定。配电网线路截面小、电压等级低，电阻与电抗接近甚至电阻更大，雅可比矩阵的对角优势减弱，初值偏离真解时可能出现迭代振荡或不收敛。工程上常采用“平启动 + 阻尼因子”或先用前推回代法给出较好初值，再切换到牛顿法。
第二，三相不平衡问题。配电网存在大量单相负荷和单相分布式电源，三相电压、电流不对称。若按单相模型计算，误差可能达到百分之几，无法满足低压配电网分析要求。此时需要建立三相节点导纳矩阵，把每相作为一个独立节点，节点数变为原来的三倍，雅可比矩阵维数随之扩大，计算量显著增加。
第三，网络规模与稀疏性利用。一个中等城市配电网可能有数千个节点，若按满矩阵处理，雅可比矩阵存储和求逆的开销不可接受。实际实现中必须采用稀疏存储（如按行压缩）和稀疏线性方程求解（如稀疏 LU 分解），并利用辐射状网络的近似三角结构进行节点编号优化，减少注入元。
第四，与替代算法的对比。配电网潮流还有前推回代法、Zbus 高斯法、回路阻抗法等。前推回代法利用辐射状结构，无需形成雅可比矩阵，每次迭代只做一次前推和一次回代，单次计算量小、内存占用低，对辐射状配电网收敛可靠，是配电网潮流最常用的算法之一。但它的收敛速度是线性的，且对含环网或 PV 节点的网络处理不便。牛顿-拉夫逊法虽然单次迭代开销大，但收敛快、对环网和 PV 节点适应性强，适合作为配电网中需要高精度、多类型节点混合场景的求解器。