基于COMSOL的输变电工程工频电磁场仿真与环评应用实践

最近在做某输变电工程的环评支撑时，正好碰到一段三相高压电缆要从居住区地下穿过。评审专家和住户最关心的就一件事：电缆埋在地下，地面的人站久了到底有没有影响？这个问题落到技术层面，其实是两个量：电缆运行时三相电流在地表附近产生的工频电场和磁场是多少；如果一个人站在电缆上方或旁边，人体内部会感应出多大的电场和电流密度。前者是环境评价的标准输出，后者才是评估生物效应真正要用的量。

这类计算用COMSOL Multiphysics做是目前比较稳妥的路线。工频电磁场本身不算复杂的物理，难在几何、材料、边界条件和激励方式的选择上，任何一个环节偷懒，结果可能差出一个数量级。这篇就从我实际搭建和调试这套分析模型的过程出发，把物理近似、建模取舍、参数设置和踩坑记录完整过一遍。

1. 从环评需求说起：这个仿真到底要回答什么问题

1.1 输变电工程环评里的电磁环境项

高压电缆和变电站的电磁环境评估，基本是所有输变电工程环评里被追问频率最高的一项。这不难理解，电缆埋在地下看不见摸不着，公众能感知的只有"听说高压线有辐射"这类模糊概念。做环评的工程师必须在报告里明确给出：线路运行时，地面及以上空间的工频电场强度、磁感应强度能达到多少数值，是否在相应导则的限值范围内，而且这个结论不能只靠定性解释，需要可复现的计算或实测数据支撑。

我之前做的模拟项目X里，有一段220kV三芯高压电缆，要沿着一条规划道路敷设，两侧几十米就是住宅区。环评大纲里明确要求给出电缆正上方、离电缆不同水平距离处的场强剖面，还要评估正常运行工况下（而不是故障工况）公众活动区域是否满足限值要求。这就意味着，常规的"算个最大值"不够，得做参数化扫描，得到电缆上方到数十米范围内的场分布曲线。

1.2 需要的输出不止是一张云图

表面上看，仿真的输出是电场强度E和磁感应强度B的分布云图，但稍往深一层想，环评关注的对象是人。工频电场和磁场对人体的作用机制完全不同：

• 工频电场作用于人体时，人体表面会感应出束缚电荷，体内反而被屏蔽得比较干净，真正有意义的量是体表外部电场畸变程度和体内通过接触电流或感应电荷形成的电流密度。
• 工频磁场能穿透人体（因为人体相对磁导率近似为1），时变磁场在体内感生出涡旋电场，进而产生感应电流密度，这是低频磁场生物效应研究的关键参量。

所以我在模型里不只看地面处的环境场强，还会把一个人体模型放进去，额外提取体内电场强度和感应电流密度。这两类数据分别对应环境评估和曝露评估，侧重点不同，但一次模型里就能同时出来。

1.3 为什么不能靠手算，也不能纯靠实测

手算解析解确实能处理无限长直导线、均匀介质、规则几何这类理想问题，工程上做快速估算够用。但实际场景里有三个麻烦：电缆三根线芯各自激励、土壤和空气是两种电导率差异巨大的介质、人体是不规则有限几何体。这三件事叠在一起，解析解基本失效，强行近似会低估人体表面电场畸变，误差可以到数倍以上。

纯靠实测也有局限。现场测量只能布点地面以上某些位置，测不到"人站在那儿时体内电场"；而且居民区环境里其他设施（水管、护栏、建筑物钢筋）都会扰动场分布，测量重现性差。仿真的价值就在于：模型里可以自由设置工况，任意提取场量，用计算结果反推现场测量点位，再结合验证实验确认模型可靠性。

2. 物理模型先厘清：准静态假设、三相激励和求解域取舍

2.1 50Hz在电缆尺度上就是准静态问题

开始建模前首先要想清楚一个问题：工频50Hz电磁场该按什么物理框架处理？电磁波在自由空间的波长约为6000km，而我们关心的电缆埋深一两米、地面观察域几十米，相比波长小了好几个数量级。这意味着在求解域内，场的相位变化几乎可以忽略，辐射效应可以忽略，位移电流相比传导电流来说也很小。

这个判断直接决定求解策略：用频域稳态分析，而不是瞬态时域分析。频域求解一次就能得到所有场量在给定频率下的幅值和相位，时域求解则需要算数十个周期再取稳态，效率和内存开销完全不是一个量级。COMSOL的AC/DC模块里，对应的接口是"磁场和电场"频域研究，内部求解的方程是包含位移电流项的频域麦克斯韦方程组，但对50Hz而言，位移电流项自动落在可忽略的数值量级上。

2.2 三相电流的复数表达与施加方式

三相高压电缆的激励是三根线芯分别流过相位相差120°的正弦电流。在频域求解器里，最直观的表达是把三相电流写成复数形式：

• A相：I_A = I0
• B相：I_B = I0 * e^(j * 2π/3)
• C相：I_C = I0 * e^(j * 4π/3)

这里的I0是电流幅值。注意一个问题：COMSOL频域研究里的电流源输入默认按幅值处理，而后处理导出时，用户经常想要的是RMS（有效值），两者相差根号2的换算系数。很多第一次做的人在这里栽跟头，云图数值和实测对比时总是对不上，回头一查就是幅值和RMS混用了。我的习惯是模型参数表里用有效值定义，施加激励时再乘根号2，这样后处理直接对照环评限值就不容易出错。

施加方式上，有两种常见选择：

• 使用"线圈"特征，把每根线芯设置为一个线圈域，指定线圈电流方向和电流大小。好处是线圈特征会自动处理净电流条件，适合后续要算损耗、电感的场景。
• 直接在导体域上施加"外部电流密度"体载荷，电流密度方向沿电缆轴向，幅值按三相复数设置。这种方式更直白，也不容易出现线圈特征和电路模块耦合带来的额外自由度问题。

我这次用的是后一种。核心原因是主题是外部场分布，不需要求解线圈端子电压和阻抗，外部电流密度方式更轻量化。

2.3 求解域截断和边界条件的思考

电缆是无限长结构（沿轴向），如果只关心一个横断面附近的场分布，可以先按二维模型算，轴向变化不进入方程，模型尺寸大幅减小。但人站在电缆正上方时，人体自身是有限长度的，会扰动原本沿轴向均匀分布的地面电场。所以完整的模型至少有两种组织方式：

• 二维横截面：适用于快速扫描电缆间距、埋深、土壤参数对地面场的影响，模型小、循环快。
• 三维模型：包含一段足够长的电缆和人体，适合精确评估人体周围的场分布和体内感应量。

无论哪种，空气域不能无限大，必须截断。截断边界设什么条件，取决于准静态问题没有向外辐射的电磁波，因此不需要用高频模拟里习惯的高吸收边界。我通常用磁绝缘和电绝缘的组合来截断，把空气域半径取到电缆埋深的5到10倍，或人体尺寸的3倍以上。边界离得太近，截断误差会通过影响远场等效电荷分布反馈到近场；离得太远，网格数量爆炸，计算时间不可接受。

3. 几何与材料的工程折中：电缆、土壤、人体都要"够用就好"

3.1 电缆几何：详细建模的边界在哪里

一根220kV单芯电缆从内到外大致有：铜导体、半导电导体屏蔽层、XLPE绝缘层、半导电绝缘屏蔽层、金属护套、外护套。如果题目的目标是算电缆内部绝缘层电场强度和局部放电风险，这些层必须全部建模，每个层界面的电场集中效应都不能忽略。但如果目标是算电缆外部（地面附近、土壤和空气中）的电磁场，电缆内部的精细层结构意义不大——从电磁场角度来看，外部场主要由线芯电流总量和各相之间的几何位置决定，导体和绝缘层内部的场分布对远场的影响可以忽略。

我这次把电缆简化为三层：导体芯（铜）、绝缘层（XLPE）、外护套，金属护套单独用一根接地导线在模型中显式建模，通过是否接地来模拟不同运行工况。对比过简化模型和完整多层模型，在地面上方1.5m处的电场磁感应强度数值差异小于1%，因此这个简化是安全的。

真正影响外部场的是三相导体之间的相对位置。高压电缆的敷设方式可能是三根单芯电缆平行排列，也可能是三芯一体成缆。我按最常见的平行排列处理，三根电缆中心距取0.3m，埋深以电缆中心到地面距离1.2m为基准。

3.2 土壤电导率：一个必须做参数扫描的材料

土壤在50Hz下的电导率范围很宽。干燥砂土可能只有0.0005 S/m，湿润黏土或回填土可达0.05 S/m，相差两个数量级。电导率影响什么？影响工频电场在地表和空气中的分布，本质上因为土壤是导体，电缆泄漏电流会通过土壤形成回流路径，土壤导电性越好，地面电场越弱；导电性越差，电缆对地相当于对大电阻充电，地面电场强。

所以我强烈建议不要把土壤电导率当固定值拍脑袋填进去，而是按高低两个数量级做参数化扫描。COMSOL里这个参数扫描可以用"扫描"研究步骤一次性跑完，每个电导率值对应一组结果，横向对比效率很高。空气和土壤的相对介电常数也建议分别设置，空气取1，土壤根据含水量取5到15之间。

3.3 人体模型：匀质体还是解剖模型

人体建到什么细节程度，取决于你要回答什么问题。如果是评估磁场在体内的感应电流密度，人体组织的磁导率基本等于空气，磁场会原样穿透人体，此时组织电导率随位置的变化是主要影响因素。50Hz下不同组织的电导率差异不小：肌肉约0.2到0.6 S/m，脂肪和骨骼约0.02到0.05 S/m。如果手边没有解剖模型，用匀质体取人体平均电导率0.2 S/m左右也能得到工程可接受的数量级估计，但体表附近的局部电场会有偏差。

如果要精确评估局部组织（比如手、脚、头部）的感应电流密度，那就得引入解剖学人体模型或者至少做皮肤、脂肪、肌肉三层的分层模型。分层模型在我的经验里精度收益和建模成本平衡得比较好。我专门做过对比：匀质人体模型在最坏工况下给出的体内最大电流密度约是分层模型的1.3倍，原因在于匀质模型的表层电导率高于真实皮肤干燥时的值，导致感应电流更集中在表面。结论要保守评估时，匀质模型反而有一定安全余量，但结论要精确时不可接受。

另外注意，人体模型在COMSOL里不需要建得太过精细。网格剖分和几何渲染精度是两回事，外轮廓用椭圆柱或稍微有些肩部过渡的造型就足够了，关键是体积和表面积要接近真实值，因为人体表面对电场的扰动主要取决于总表面积和曲率分布。

3.4 二维和三维模型的配合使用

二维横截面模型是二维的，无法放置一个真实的三维人体。我的做法是先跑二维模型做参数扫描，确定电缆间距、埋深、土壤电导率等宏观参数的最优设置，然后把固定参数带入三维模型，三维模型才引入人体。

这样安排的好处是计算资源用在刀刃上。二维扫描可以同时算几十个工况，几分钟内完成；三维模型网格量动辄数百万单元，每个工况跑几十分钟甚至几小时，只能精算少数几个代表工况。这个"二维扫参、三维精算"的策略，在时间和资源上都比上来就全三维要合理得多。

4. 从零搭模型：参数表、几何、网格和求解器的具体配置

4.1 全局参数表先立规矩

我习惯在模型最开始就把所有可调量集中到一个参数表里，这样后面做参数扫描时只需改参数名，不用动几何和物理场设置。下面是模拟项目X里用到的核心参数清单（数值为示意）：

三相电流幅值用 I0 = I_rms * sqrt(2) 计算，相位差按120°递增。这个表里每个参数都值得留一个位置，因为后期做网格收敛性验证时，你会希望固定所有工况参数只改网格尺寸，而参数表能帮你明确复现别人或自己之前的结果。

4.2 几何搭建与网格控制

二维模型里，几何由三个同心圆（电缆）和两个同心矩形（土壤和空气）组成，矩形间交界面一条水平线代表地面。三维模型则在此基础上将三个同心圆沿轴向拉伸一段长度（我取10m，两端用对称边界条件模拟无限长电缆），再把人体放进去。

网格是工频电磁场仿真里最容易被低估的工作量。关键有三处：

• 电缆导体和绝缘层的交界面附近，电场变化梯度大，需要足够细的网格分辨率，至少让每个圆环截面有6到8层单元。
• 人体表面，尤其是头部、肩部、脚部这样的曲面区域，需要加边界层网格，否则体表电场畸变无法正确捕获。
• 地面与空气交界面，要保证两侧网格能平滑过渡，我用的是地面以下加密、地面以上逐渐加粗的渐变网格。

我建议第一次运行前先跑一个"粗网格-中网格-细网格"三档的收敛性算例，监控目标量是人体头部中心点的电场模值。把网格尺寸减半后这个量变化小于5%，就可以认为网格收敛。这一步看起来麻烦，实际做下来是省时间的，因为发现某些区域的电场峰值随网格变化剧烈时，说明网格还没剖够，这时候的结果不能拿去写结论。

4.3 物理场与边界条件的设置顺序

物理场接口选择"磁场和电场"频域研究，然后按以下顺序设置：

1. 把土壤、空气、人体、电缆各部分材料参数填入各自域。
2. 给三根电缆导体域施加外部电流密度，沿电缆轴向（设为z方向），幅值和相位按4.1的参数表。
3. 设置求解域外边界为电绝缘和磁绝缘组合。
4. 设置地面作为电位参考零点（或者用悬浮电位默认条件，视求解器要求）。
5. 研究步骤选频域，频率填入50 Hz，必要时同时填入150 Hz等谐波频率，可以一次算出谐波叠加下的总场。

如果使用了"线圈"特征，还要多一步在线圈设置窗口里勾选"线圈类型"和"激励方式"，并把电流方向指定为轴向。外部电流密度方式没有这一步，我个人更推荐后者用于外部场分析。

4.4 求解器选择：直接求解器更稳，大型模型要配合迭代

工频频域问题求解的是复数线性方程组，矩阵规模取决于网格数。二维模型几十万自由度，直接求解器（MUMPS或PARDISO）都毫无压力。三维模型加人体后自由度可能到数百万，内存吃紧时要考虑迭代求解器。

我的经验是优先PARDISO，它在相同内存下通常比MUMPS快，而且多核扩展性好。当模型大到PARDISO内存不足时，换GMRES迭代求解器，配合几何多重网格或代数多重网格预处理。迭代法的调试成本较高，一般通过增加迭代次数上限和调整容差来控制精度，第一次跑大型模型时务必盯一下求解报告里的收敛残差，不收敛就换预处理方式。

另外提醒一个容易忽略的点：频域求解的默认容差比静力学要求更严格，因为复数解要同时满足幅值和相位。如果只是粗看云图，默认容差够用；但如果要提取某个点的精确数值用于和限值比较，建议把求解容差调紧一档，否则相位误差会被放大到幅值里。

5. 结果判读与限值对照：云图上最容易看错的地方

5.1 输出量的选择：电场模、磁感应强度模、感应电流密度

结果的判读先要把单位和物理量搞对。频域解出来的场量通常是复数，COMSOL后处理里默认可以显示模、实部、虚部或RMS等效。对于环评里的限值对照，我习惯全部用RMS值显示，这和标准限值的量纲一致。

主要输出量有三个：

• 地表上方1.5m高度的电场强度E（RMS），环境评价标准里工程上常取该高度的代表值。
• 磁感应强度B（RMS）在地表和人体周围的空间分布。
• 人体内部的感应电流密度J（RMS）和体内电场强度E（RMS），用于公众曝露评估。

磁感应强度还有另一个表达形式是dB/dt（磁通密度变化率），因为感应电流密度正比于dB/dt而非B本身。不过环评导则通常直接给出B限值，数据对比时用B就够了，关注机理研究时可以再导出dB/dt。

5.2 人体内外场强的数量级差异，云图要分区显示

第一次跑完三维模型，如果直接用全模型的默认色标看电场云图，几乎等于没看。原因在于人体外部电场可以达到数千V/m，而人体内部电场只有几十mV/m，相差近5个数量级，默认色标会把所有内部信息压成一个颜色。

处理办法是分两个视图：一个视图只显示空气和土壤区域，色标范围设到0到5000 V/m，看环境场；另一个视图只显示人体内部，色标范围设到0到0.5 V/m，看感应电场。如果软件版本支持，还可以用"镜像"或"表格"把不同区域的值分开统计。这个细节不处理好，写报告时非常容易给出误导性的彩色图。

磁感应强度则没有这个困扰，人体内外B几乎一致（相对磁导率差异可忽略），云图可以一张图通吃。

5.3 参数化扫描的数据提取方式

做参数扫描时，我不建议每跑完一个工况就盯着云图看，而是在模型里提前定义好探针，提取固定位置的场量。比如地面以上0.5m、1.5m、3m三个高度、水平距离从电缆正上方到20m等间距取点，把这组数据导出成表格。COMSOL的"派生值"里可以用"线图"或"剖线计算"实现。

土壤电导率的参数扫描结果很有意思：电导率从0.0005扫到0.05 S/m时，地面电场可以差出3倍以上，磁场几乎不变。这说明电场对土壤参数高度敏感，而磁场主要由电缆电流决定。做工程评估时，如果拿不准土壤的确切电导率，最稳妥的做法是把最不利（电导率最低）的结果作为报告值，并注明土壤电导率假设。

5.4 限值对照的常见误区

对照限值需要分清楚两个层面：

第一层是环境场强限值。这个值是针对公众活动空间的电场强度或磁感应强度制定的，不能直接把人体内部感应电流密度拿去和它比，因为两者的物理量纲和关注对象都不同。

第二层是人体曝露限值。这一层才看人体内部的感应电流密度或体内电场。不同频率对应的限值形式也不同，工频段通常是感应电流密度或体内电场，高频段则更多看SAR比吸收率。所以判断"是否超标"前，先要确认拿到的仿真结果是环境场还是人体内场。

我在报告里常用一张对照表把仿真值和限值摆在一起：

表里数值看的是数量级差距。关键是要把"安全"的结论建立在仿真工况覆盖了最不利假设的基础上，不能只挑一个"看起来最安全"的工况下结论。

6. 实操中踩过的坑：三相回路、网格收敛与两阶段法的取舍

6.1 三相回路不闭合，磁场会整体漂移

第一次跑三维模型时，我给三根电缆都单独施加了外部电流密度，三个电流加起来恰好等于零，物理上没有问题。但COMSOL里如果只有电流密度源，默认情况下磁场回路的边界条件由求解域外边界承担，这会引出在截断边界处磁场沿边界切向自动为零的约束，实际等效为人为设置了回流路径。

结果就是，离电缆较远处（约10m外）的磁感应强度本应衰减到非常低，但云图显示出了一个接近常数的背景B，明显不符合物理直觉。排查后发现问题是这样的：三相电流虽然净和为零，但三条线芯在空间上分开布置，局部电流密度并非处处为零，在截断边界处强制切向磁场的条件，相当于人为把回流"拉"到了远边界，导致整个求解域叠加上了一个虚假的匀强磁场。

解决办法是在求解域内增加一根"回流导线"，把这根导线设定为零电流或接地，让它承担三相不平衡电流和位移电流的返回路径。更严谨的做法是使用"线圈"特征并打开"多匝"或"回路"选项，让软件通过净电流约束自己找闭环。这个坑不踩一次很难想到，但一旦磁场云图上出现不该存在的背景值，就要首先怀疑回路的物理完整性。

6.2 网格收敛性：哪些量敏感，哪些量不敏感

网格敏感性在不同物理量上的差异非常大。我在测试中发现的规律是：

• 人体表面电场峰值对网格最敏感，边界层单元层数从3层提到6层，峰值可以变化15%左右，必须专门收敛。
• 人体内部最大感应电流密度对网格的敏感度次之，变化约5%到8%。
• 地面以上1.5m处的电场和磁场对网格相对不敏感，粗网格和细网格结果差不到3%。

这给工程实践带来一个启示：不要只收敛一个"代表性"物理量，要看最终写进报告的每个量是否都收敛。我见过有同事只看了地面电场收敛就出报告，结果人体内场其实还没收敛。评判网格是否足够的唯一标准，是被评估量本身随网格变化的幅度，而不是整体云图是否好看。

6.3 全模型内存不够时的两阶段法

三维模型加精细人体网格后，内存很容易冲破本机限制。两个减少内存的常用手段：

一是几何对称性利用。如果人体站在电缆正上方，且电缆沿轴向均匀，可以只建四分之一模型，在对称面上加对称边界条件。工频三相电流不是对称的，所以这里的"对称"必须基于时均意义下的RMS场，而不是瞬时场。

二是背景场-子模型解耦法。第一阶段先去掉人体，只算包含电缆、土壤和空气的背景场，把结果保存在数据集里；第二阶段建一个人体加周围局部空气的小模型，在子模型外边界上施加第一阶段的场值作为边界条件，只求解人体附近的扰动场。这样模型规模从全局缩小到局部，内存需求骤降。

两阶段法的物理依据是：人体和电缆之间距离至少一米以上，人体对电缆电流回路的反作用可以忽略，因此背景场不受人体存在影响。但使用时有条件：子模型边界必须离人体足够远，通常至少2倍人体高度，否则截断误差会污染体内场。我在模拟项目X里用这个方法把内存占用降低了近70%，结果与全模型对比在1%以内的偏差，工程上完全可用。

6.4 与实测对比时的单位换算细节

仿真报告最后总要面对现场测试的验证环节。现场用场强仪测出来的读数按定义是RMS有效值，而COMSOL频域云图默认显示的有时是幅值（峰值）。如果直接拿默认云图去对比实测，会高出根号2倍，这个错误真的发生过很多次。

我在后处理设置里把电场和磁场显示都改成RMS，同时在报告里注明"所有仿真数据均为RMS有效值，与现场场强仪读数口径一致"。另外要提醒一点：磁场测量探头如果装在人体附近，人体会对低频电场产生严重扰动，但磁场几乎不受人体影响，所以现场测试时电场和磁场的可信度是不同级别的，对比分析时要区别对待。

结语

如果只挑出一条最想分享的经验，那就是：这类电磁场分析项目，真正决定结果可信度的往往不是求解器本身够不够高级，而是建模者对物理问题边界的理解是否准确。三相回路要完整，土壤参数要做扫描，人体网格要单独收敛，限值对照要分清环境场和人体内场——每一条都是我在实际项目里真金白银换来的教训。COMSOL给你的是精确求解方程的能力，但方程本身是不是对应你要回答的问题，这件事永远得自己把关。下次做类似的电缆电磁场评估时，建议先从二维模型开始快速建立直觉，再谨慎地移入三维精算，你会发现整个过程顺畅得多。