最近在做某输变电工程的环评支撑时正好碰到一段三相高压电缆要从居住区地下穿过。评审专家和住户最关心的就一件事电缆埋在地下地面的人站久了到底有没有影响这个问题落到技术层面其实是两个量电缆运行时三相电流在地表附近产生的工频电场和磁场是多少如果一个人站在电缆上方或旁边人体内部会感应出多大的电场和电流密度。前者是环境评价的标准输出后者才是评估生物效应真正要用的量。这类计算用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 I0B相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里用到的核心参数清单数值为示意参数含义数值f0工频频率50 HzI_rms单相线芯电流有效值600 Ar_cond导体半径18 mmt_ins绝缘层厚度24 mmt_sheath外护套厚度6 mmcable_spacing相邻电缆中心距300 mmdepth电缆中心埋深1.2 msigma_soil土壤电导率0.005 S/m扫描基准sigma_body人体均匀电导率0.2 S/mh_body人体身高1.75 mr_air空气域半径20 m三相电流幅值用 I0 I_rms * sqrt(2) 计算相位差按120°递增。这个表里每个参数都值得留一个位置因为后期做网格收敛性验证时你会希望固定所有工况参数只改网格尺寸而参数表能帮你明确复现别人或自己之前的结果。4.2 几何搭建与网格控制二维模型里几何由三个同心圆电缆和两个同心矩形土壤和空气组成矩形间交界面一条水平线代表地面。三维模型则在此基础上将三个同心圆沿轴向拉伸一段长度我取10m两端用对称边界条件模拟无限长电缆再把人体放进去。网格是工频电磁场仿真里最容易被低估的工作量。关键有三处电缆导体和绝缘层的交界面附近电场变化梯度大需要足够细的网格分辨率至少让每个圆环截面有6到8层单元。人体表面尤其是头部、肩部、脚部这样的曲面区域需要加边界层网格否则体表电场畸变无法正确捕获。地面与空气交界面要保证两侧网格能平滑过渡我用的是地面以下加密、地面以上逐渐加粗的渐变网格。我建议第一次运行前先跑一个粗网格-中网格-细网格三档的收敛性算例监控目标量是人体头部中心点的电场模值。把网格尺寸减半后这个量变化小于5%就可以认为网格收敛。这一步看起来麻烦实际做下来是省时间的因为发现某些区域的电场峰值随网格变化剧烈时说明网格还没剖够这时候的结果不能拿去写结论。4.3 物理场与边界条件的设置顺序物理场接口选择磁场和电场频域研究然后按以下顺序设置把土壤、空气、人体、电缆各部分材料参数填入各自域。给三根电缆导体域施加外部电流密度沿电缆轴向设为z方向幅值和相位按4.1的参数表。设置求解域外边界为电绝缘和磁绝缘组合。设置地面作为电位参考零点或者用悬浮电位默认条件视求解器要求。研究步骤选频域频率填入50 Hz必要时同时填入150 Hz等谐波频率可以一次算出谐波叠加下的总场。如果使用了线圈特征还要多一步在线圈设置窗口里勾选线圈类型和激励方式并把电流方向指定为轴向。外部电流密度方式没有这一步我个人更推荐后者用于外部场分析。4.4 求解器选择直接求解器更稳大型模型要配合迭代工频频域问题求解的是复数线性方程组矩阵规模取决于网格数。二维模型几十万自由度直接求解器MUMPS或PARDISO都毫无压力。三维模型加人体后自由度可能到数百万内存吃紧时要考虑迭代求解器。我的经验是优先PARDISO它在相同内存下通常比MUMPS快而且多核扩展性好。当模型大到PARDISO内存不足时换GMRES迭代求解器配合几何多重网格或代数多重网格预处理。迭代法的调试成本较高一般通过增加迭代次数上限和调整容差来控制精度第一次跑大型模型时务必盯一下求解报告里的收敛残差不收敛就换预处理方式。另外提醒一个容易忽略的点频域求解的默认容差比静力学要求更严格因为复数解要同时满足幅值和相位。如果只是粗看云图默认容差够用但如果要提取某个点的精确数值用于和限值比较建议把求解容差调紧一档否则相位误差会被放大到幅值里。5. 结果判读与限值对照云图上最容易看错的地方5.1 输出量的选择电场模、磁感应强度模、感应电流密度结果的判读先要把单位和物理量搞对。频域解出来的场量通常是复数COMSOL后处理里默认可以显示模、实部、虚部或RMS等效。对于环评里的限值对照我习惯全部用RMS值显示这和标准限值的量纲一致。主要输出量有三个地表上方1.5m高度的电场强度ERMS环境评价标准里工程上常取该高度的代表值。磁感应强度BRMS在地表和人体周围的空间分布。人体内部的感应电流密度JRMS和体内电场强度ERMS用于公众曝露评估。磁感应强度还有另一个表达形式是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比吸收率。所以判断是否超标前先要确认拿到的仿真结果是环境场还是人体内场。我在报告里常用一张对照表把仿真值和限值摆在一起评估对象仿真提取量RMS典型仿真值示意限值对比地面电场公众活动区地面1.5m高度E数百V/m量级低于公众限值地面磁场公众活动区地面B个位数到十几µT低于公众限值人体体内感应电流密度人体最大J低于基本限值显著低于限值表里数值看的是数量级差距。关键是要把安全的结论建立在仿真工况覆盖了最不利假设的基础上不能只挑一个看起来最安全的工况下结论。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给你的是精确求解方程的能力但方程本身是不是对应你要回答的问题这件事永远得自己把关。下次做类似的电缆电磁场评估时建议先从二维模型开始快速建立直觉再谨慎地移入三维精算你会发现整个过程顺畅得多。
