最近在做某輸變電工程的環(huán)評支撐時正好碰到一段三相高壓電纜要從居住區(qū)地下穿過。評審專家和住戶最關心的就一件事電纜埋在地下地面的人站久了到底有沒有影響這個問題落到技術層面其實是兩個量電纜運行時三相電流在地表附近產生的工頻電場和磁場是多少如果一個人站在電纜上方或旁邊人體內部會感應出多大的電場和電流密度。前者是環(huán)境評價的標準輸出后者才是評估生物效應真正要用的量。這類計算用COMSOL Multiphysics做是目前比較穩(wěn)妥的路線。工頻電磁場本身不算復雜的物理難在幾何、材料、邊界條件和激勵方式的選擇上任何一個環(huán)節(jié)偷懶結果可能差出一個數(shù)量級。這篇就從我實際搭建和調試這套分析模型的過程出發(fā)把物理近似、建模取舍、參數(shù)設置和踩坑記錄完整過一遍。1. 從環(huán)評需求說起這個仿真到底要回答什么問題1.1 輸變電工程環(huán)評里的電磁環(huán)境項高壓電纜和變電站的電磁環(huán)境評估基本是所有輸變電工程環(huán)評里被追問頻率最高的一項。這不難理解電纜埋在地下看不見摸不著公眾能感知的只有聽說高壓線有輻射這類模糊概念。做環(huán)評的工程師必須在報告里明確給出線路運行時地面及以上空間的工頻電場強度、磁感應強度能達到多少數(shù)值是否在相應導則的限值范圍內而且這個結論不能只靠定性解釋需要可復現(xiàn)的計算或實測數(shù)據(jù)支撐。我之前做的模擬項目X里有一段220kV三芯高壓電纜要沿著一條規(guī)劃道路敷設兩側幾十米就是住宅區(qū)。環(huán)評大綱里明確要求給出電纜正上方、離電纜不同水平距離處的場強剖面還要評估正常運行工況下而不是故障工況公眾活動區(qū)域是否滿足限值要求。這就意味著常規(guī)的算個最大值不夠得做參數(shù)化掃描得到電纜上方到數(shù)十米范圍內的場分布曲線。1.2 需要的輸出不止是一張云圖表面上看仿真的輸出是電場強度E和磁感應強度B的分布云圖但稍往深一層想環(huán)評關注的對象是人。工頻電場和磁場對人體的作用機制完全不同工頻電場作用于人體時人體表面會感應出束縛電荷體內反而被屏蔽得比較干凈真正有意義的量是體表外部電場畸變程度和體內通過接觸電流或感應電荷形成的電流密度。工頻磁場能穿透人體因為人體相對磁導率近似為1時變磁場在體內感生出渦旋電場進而產生感應電流密度這是低頻磁場生物效應研究的關鍵參量。所以我在模型里不只看地面處的環(huán)境場強還會把一個人體模型放進去額外提取體內電場強度和感應電流密度。這兩類數(shù)據(jù)分別對應環(huán)境評估和曝露評估側重點不同但一次模型里就能同時出來。1.3 為什么不能靠手算也不能純靠實測手算解析解確實能處理無限長直導線、均勻介質、規(guī)則幾何這類理想問題工程上做快速估算夠用。但實際場景里有三個麻煩電纜三根線芯各自激勵、土壤和空氣是兩種電導率差異巨大的介質、人體是不規(guī)則有限幾何體。這三件事疊在一起解析解基本失效強行近似會低估人體表面電場畸變誤差可以到數(shù)倍以上。純靠實測也有局限?，F(xiàn)場測量只能布點地面以上某些位置測不到人站在那兒時體內電場而且居民區(qū)環(huán)境里其他設施水管、護欄、建筑物鋼筋都會擾動場分布測量重現(xiàn)性差。仿真的價值就在于模型里可以自由設置工況任意提取場量用計算結果反推現(xiàn)場測量點位再結合驗證實驗確認模型可靠性。2. 物理模型先厘清準靜態(tài)假設、三相激勵和求解域取舍2.1 50Hz在電纜尺度上就是準靜態(tài)問題開始建模前首先要想清楚一個問題工頻50Hz電磁場該按什么物理框架處理電磁波在自由空間的波長約為6000km而我們關心的電纜埋深一兩米、地面觀察域幾十米相比波長小了好幾個數(shù)量級。這意味著在求解域內場的相位變化幾乎可以忽略輻射效應可以忽略位移電流相比傳導電流來說也很小。這個判斷直接決定求解策略用頻域穩(wěn)態(tài)分析而不是瞬態(tài)時域分析。頻域求解一次就能得到所有場量在給定頻率下的幅值和相位時域求解則需要算數(shù)十個周期再取穩(wěn)態(tài)效率和內存開銷完全不是一個量級。COMSOL的AC/DC模塊里對應的接口是磁場和電場頻域研究內部求解的方程是包含位移電流項的頻域麥克斯韋方程組但對50Hz而言位移電流項自動落在可忽略的數(shù)值量級上。2.2 三相電流的復數(shù)表達與施加方式三相高壓電纜的激勵是三根線芯分別流過相位相差120°的正弦電流。在頻域求解器里最直觀的表達是把三相電流寫成復數(shù)形式A相I_A I0B相I_B I0 * e^(j * 2π/3)C相I_C I0 * e^(j * 4π/3)這里的I0是電流幅值。注意一個問題COMSOL頻域研究里的電流源輸入默認按幅值處理而后處理導出時用戶經常想要的是RMS有效值兩者相差根號2的換算系數(shù)。很多第一次做的人在這里栽跟頭云圖數(shù)值和實測對比時總是對不上回頭一查就是幅值和RMS混用了。我的習慣是模型參數(shù)表里用有效值定義施加激勵時再乘根號2這樣后處理直接對照環(huán)評限值就不容易出錯。施加方式上有兩種常見選擇使用線圈特征把每根線芯設置為一個線圈域指定線圈電流方向和電流大小。好處是線圈特征會自動處理凈電流條件適合后續(xù)要算損耗、電感的場景。直接在導體域上施加外部電流密度體載荷電流密度方向沿電纜軸向幅值按三相復數(shù)設置。這種方式更直白也不容易出現(xiàn)線圈特征和電路模塊耦合帶來的額外自由度問題。我這次用的是后一種。核心原因是主題是外部場分布不需要求解線圈端子電壓和阻抗外部電流密度方式更輕量化。2.3 求解域截斷和邊界條件的思考電纜是無限長結構沿軸向如果只關心一個橫斷面附近的場分布可以先按二維模型算軸向變化不進入方程模型尺寸大幅減小。但人站在電纜正上方時人體自身是有限長度的會擾動原本沿軸向均勻分布的地面電場。所以完整的模型至少有兩種組織方式二維橫截面適用于快速掃描電纜間距、埋深、土壤參數(shù)對地面場的影響模型小、循環(huán)快。三維模型包含一段足夠長的電纜和人體適合精確評估人體周圍的場分布和體內感應量。無論哪種空氣域不能無限大必須截斷。截斷邊界設什么條件取決于準靜態(tài)問題沒有向外輻射的電磁波因此不需要用高頻模擬里習慣的高吸收邊界。我通常用磁絕緣和電絕緣的組合來截斷把空氣域半徑取到電纜埋深的5到10倍或人體尺寸的3倍以上。邊界離得太近截斷誤差會通過影響遠場等效電荷分布反饋到近場離得太遠網格數(shù)量爆炸計算時間不可接受。3. 幾何與材料的工程折中電纜、土壤、人體都要夠用就好3.1 電纜幾何詳細建模的邊界在哪里一根220kV單芯電纜從內到外大致有銅導體、半導電導體屏蔽層、XLPE絕緣層、半導電絕緣屏蔽層、金屬護套、外護套。如果題目的目標是算電纜內部絕緣層電場強度和局部放電風險這些層必須全部建模每個層界面的電場集中效應都不能忽略。但如果目標是算電纜外部地面附近、土壤和空氣中的電磁場電纜內部的精細層結構意義不大——從電磁場角度來看外部場主要由線芯電流總量和各相之間的幾何位置決定導體和絕緣層內部的場分布對遠場的影響可以忽略。我這次把電纜簡化為三層導體芯銅、絕緣層XLPE、外護套金屬護套單獨用一根接地導線在模型中顯式建模通過是否接地來模擬不同運行工況。對比過簡化模型和完整多層模型在地面上方1.5m處的電場磁感應強度數(shù)值差異小于1%因此這個簡化是安全的。真正影響外部場的是三相導體之間的相對位置。高壓電纜的敷設方式可能是三根單芯電纜平行排列也可能是三芯一體成纜。我按最常見的平行排列處理三根電纜中心距取0.3m埋深以電纜中心到地面距離1.2m為基準。3.2 土壤電導率一個必須做參數(shù)掃描的材料土壤在50Hz下的電導率范圍很寬。干燥砂土可能只有0.0005 S/m濕潤黏土或回填土可達0.05 S/m相差兩個數(shù)量級。電導率影響什么影響工頻電場在地表和空氣中的分布本質上因為土壤是導體電纜泄漏電流會通過土壤形成回流路徑土壤導電性越好地面電場越弱導電性越差電纜對地相當于對大電阻充電地面電場強。所以我強烈建議不要把土壤電導率當固定值拍腦袋填進去而是按高低兩個數(shù)量級做參數(shù)化掃描。COMSOL里這個參數(shù)掃描可以用掃描研究步驟一次性跑完每個電導率值對應一組結果橫向對比效率很高。空氣和土壤的相對介電常數(shù)也建議分別設置空氣取1土壤根據(jù)含水量取5到15之間。3.3 人體模型勻質體還是解剖模型人體建到什么細節(jié)程度取決于你要回答什么問題。如果是評估磁場在體內的感應電流密度人體組織的磁導率基本等于空氣磁場會原樣穿透人體此時組織電導率隨位置的變化是主要影響因素。50Hz下不同組織的電導率差異不小肌肉約0.2到0.6 S/m脂肪和骨骼約0.02到0.05 S/m。如果手邊沒有解剖模型用勻質體取人體平均電導率0.2 S/m左右也能得到工程可接受的數(shù)量級估計但體表附近的局部電場會有偏差。如果要精確評估局部組織比如手、腳、頭部的感應電流密度那就得引入解剖學人體模型或者至少做皮膚、脂肪、肌肉三層的分層模型。分層模型在我的經驗里精度收益和建模成本平衡得比較好。我專門做過對比勻質人體模型在最壞工況下給出的體內最大電流密度約是分層模型的1.3倍原因在于勻質模型的表層電導率高于真實皮膚干燥時的值導致感應電流更集中在表面。結論要保守評估時勻質模型反而有一定安全余量但結論要精確時不可接受。另外注意人體模型在COMSOL里不需要建得太過精細。網格剖分和幾何渲染精度是兩回事外輪廓用橢圓柱或稍微有些肩部過渡的造型就足夠了關鍵是體積和表面積要接近真實值因為人體表面對電場的擾動主要取決于總表面積和曲率分布。3.4 二維和三維模型的配合使用二維橫截面模型是二維的無法放置一個真實的三維人體。我的做法是先跑二維模型做參數(shù)掃描確定電纜間距、埋深、土壤電導率等宏觀參數(shù)的最優(yōu)設置然后把固定參數(shù)帶入三維模型三維模型才引入人體。這樣安排的好處是計算資源用在刀刃上。二維掃描可以同時算幾十個工況幾分鐘內完成三維模型網格量動輒數(shù)百萬單元每個工況跑幾十分鐘甚至幾小時只能精算少數(shù)幾個代表工況。這個二維掃參、三維精算的策略在時間和資源上都比上來就全三維要合理得多。4. 從零搭模型參數(shù)表、幾何、網格和求解器的具體配置4.1 全局參數(shù)表先立規(guī)矩我習慣在模型最開始就把所有可調量集中到一個參數(shù)表里這樣后面做參數(shù)掃描時只需改參數(shù)名不用動幾何和物理場設置。下面是模擬項目X里用到的核心參數(shù)清單數(shù)值為示意參數(shù)含義數(shù)值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°遞增。這個表里每個參數(shù)都值得留一個位置因為后期做網格收斂性驗證時你會希望固定所有工況參數(shù)只改網格尺寸而參數(shù)表能幫你明確復現(xiàn)別人或自己之前的結果。4.2 幾何搭建與網格控制二維模型里幾何由三個同心圓電纜和兩個同心矩形土壤和空氣組成矩形間交界面一條水平線代表地面。三維模型則在此基礎上將三個同心圓沿軸向拉伸一段長度我取10m兩端用對稱邊界條件模擬無限長電纜再把人體放進去。網格是工頻電磁場仿真里最容易被低估的工作量。關鍵有三處電纜導體和絕緣層的交界面附近電場變化梯度大需要足夠細的網格分辨率至少讓每個圓環(huán)截面有6到8層單元。人體表面尤其是頭部、肩部、腳部這樣的曲面區(qū)域需要加邊界層網格否則體表電場畸變無法正確捕獲。地面與空氣交界面要保證兩側網格能平滑過渡我用的是地面以下加密、地面以上逐漸加粗的漸變網格。我建議第一次運行前先跑一個粗網格-中網格-細網格三檔的收斂性算例監(jiān)控目標量是人體頭部中心點的電場模值。把網格尺寸減半后這個量變化小于5%就可以認為網格收斂。這一步看起來麻煩實際做下來是省時間的因為發(fā)現(xiàn)某些區(qū)域的電場峰值隨網格變化劇烈時說明網格還沒剖夠這時候的結果不能拿去寫結論。4.3 物理場與邊界條件的設置順序物理場接口選擇磁場和電場頻域研究然后按以下順序設置把土壤、空氣、人體、電纜各部分材料參數(shù)填入各自域。給三根電纜導體域施加外部電流密度沿電纜軸向設為z方向幅值和相位按4.1的參數(shù)表。設置求解域外邊界為電絕緣和磁絕緣組合。設置地面作為電位參考零點或者用懸浮電位默認條件視求解器要求。研究步驟選頻域頻率填入50 Hz必要時同時填入150 Hz等諧波頻率可以一次算出諧波疊加下的總場。如果使用了線圈特征還要多一步在線圈設置窗口里勾選線圈類型和激勵方式并把電流方向指定為軸向。外部電流密度方式沒有這一步我個人更推薦后者用于外部場分析。4.4 求解器選擇直接求解器更穩(wěn)大型模型要配合迭代工頻頻域問題求解的是復數(shù)線性方程組矩陣規(guī)模取決于網格數(shù)。二維模型幾十萬自由度直接求解器MUMPS或PARDISO都毫無壓力。三維模型加人體后自由度可能到數(shù)百萬內存吃緊時要考慮迭代求解器。我的經驗是優(yōu)先PARDISO它在相同內存下通常比MUMPS快而且多核擴展性好。當模型大到PARDISO內存不足時換GMRES迭代求解器配合幾何多重網格或代數(shù)多重網格預處理。迭代法的調試成本較高一般通過增加迭代次數(shù)上限和調整容差來控制精度第一次跑大型模型時務必盯一下求解報告里的收斂殘差不收斂就換預處理方式。另外提醒一個容易忽略的點頻域求解的默認容差比靜力學要求更嚴格因為復數(shù)解要同時滿足幅值和相位。如果只是粗看云圖默認容差夠用但如果要提取某個點的精確數(shù)值用于和限值比較建議把求解容差調緊一檔否則相位誤差會被放大到幅值里。5. 結果判讀與限值對照云圖上最容易看錯的地方5.1 輸出量的選擇電場模、磁感應強度模、感應電流密度結果的判讀先要把單位和物理量搞對。頻域解出來的場量通常是復數(shù)COMSOL后處理里默認可以顯示模、實部、虛部或RMS等效。對于環(huán)評里的限值對照我習慣全部用RMS值顯示這和標準限值的量綱一致。主要輸出量有三個地表上方1.5m高度的電場強度ERMS環(huán)境評價標準里工程上常取該高度的代表值。磁感應強度BRMS在地表和人體周圍的空間分布。人體內部的感應電流密度JRMS和體內電場強度ERMS用于公眾曝露評估。磁感應強度還有另一個表達形式是dB/dt磁通密度變化率因為感應電流密度正比于dB/dt而非B本身。不過環(huán)評導則通常直接給出B限值數(shù)據(jù)對比時用B就夠了關注機理研究時可以再導出dB/dt。5.2 人體內外場強的數(shù)量級差異云圖要分區(qū)顯示第一次跑完三維模型如果直接用全模型的默認色標看電場云圖幾乎等于沒看。原因在于人體外部電場可以達到數(shù)千V/m而人體內部電場只有幾十mV/m相差近5個數(shù)量級默認色標會把所有內部信息壓成一個顏色。處理辦法是分兩個視圖一個視圖只顯示空氣和土壤區(qū)域色標范圍設到0到5000 V/m看環(huán)境場另一個視圖只顯示人體內部色標范圍設到0到0.5 V/m看感應電場。如果軟件版本支持還可以用鏡像或表格把不同區(qū)域的值分開統(tǒng)計。這個細節(jié)不處理好寫報告時非常容易給出誤導性的彩色圖。磁感應強度則沒有這個困擾人體內外B幾乎一致相對磁導率差異可忽略云圖可以一張圖通吃。5.3 參數(shù)化掃描的數(shù)據(jù)提取方式做參數(shù)掃描時我不建議每跑完一個工況就盯著云圖看而是在模型里提前定義好探針提取固定位置的場量。比如地面以上0.5m、1.5m、3m三個高度、水平距離從電纜正上方到20m等間距取點把這組數(shù)據(jù)導出成表格。COMSOL的派生值里可以用線圖或剖線計算實現(xiàn)。土壤電導率的參數(shù)掃描結果很有意思電導率從0.0005掃到0.05 S/m時地面電場可以差出3倍以上磁場幾乎不變。這說明電場對土壤參數(shù)高度敏感而磁場主要由電纜電流決定。做工程評估時如果拿不準土壤的確切電導率最穩(wěn)妥的做法是把最不利電導率最低的結果作為報告值并注明土壤電導率假設。5.4 限值對照的常見誤區(qū)對照限值需要分清楚兩個層面第一層是環(huán)境場強限值。這個值是針對公眾活動空間的電場強度或磁感應強度制定的不能直接把人體內部感應電流密度拿去和它比因為兩者的物理量綱和關注對象都不同。第二層是人體曝露限值。這一層才看人體內部的感應電流密度或體內電場。不同頻率對應的限值形式也不同工頻段通常是感應電流密度或體內電場高頻段則更多看SAR比吸收率。所以判斷是否超標前先要確認拿到的仿真結果是環(huán)境場還是人體內場。我在報告里常用一張對照表把仿真值和限值擺在一起評估對象仿真提取量RMS典型仿真值示意限值對比地面電場公眾活動區(qū)地面1.5m高度E數(shù)百V/m量級低于公眾限值地面磁場公眾活動區(qū)地面B個位數(shù)到十幾μT低于公眾限值人體體內感應電流密度人體最大J低于基本限值顯著低于限值表里數(shù)值看的是數(shù)量級差距。關鍵是要把安全的結論建立在仿真工況覆蓋了最不利假設的基礎上不能只挑一個看起來最安全的工況下結論。6. 實操中踩過的坑三相回路、網格收斂與兩階段法的取舍6.1 三相回路不閉合磁場會整體漂移第一次跑三維模型時我給三根電纜都單獨施加了外部電流密度三個電流加起來恰好等于零物理上沒有問題。但COMSOL里如果只有電流密度源默認情況下磁場回路的邊界條件由求解域外邊界承擔這會引出在截斷邊界處磁場沿邊界切向自動為零的約束實際等效為人為設置了回流路徑。結果就是離電纜較遠處約10m外的磁感應強度本應衰減到非常低但云圖顯示出了一個接近常數(shù)的背景B明顯不符合物理直覺。排查后發(fā)現(xiàn)問題是這樣的三相電流雖然凈和為零但三條線芯在空間上分開布置局部電流密度并非處處為零在截斷邊界處強制切向磁場的條件相當于人為把回流拉到了遠邊界導致整個求解域疊加上了一個虛假的勻強磁場。解決辦法是在求解域內增加一根回流導線把這根導線設定為零電流或接地讓它承擔三相不平衡電流和位移電流的返回路徑。更嚴謹?shù)淖龇ㄊ鞘褂镁€圈特征并打開多匝或回路選項讓軟件通過凈電流約束自己找閉環(huán)。這個坑不踩一次很難想到但一旦磁場云圖上出現(xiàn)不該存在的背景值就要首先懷疑回路的物理完整性。6.2 網格收斂性哪些量敏感哪些量不敏感網格敏感性在不同物理量上的差異非常大。我在測試中發(fā)現(xiàn)的規(guī)律是人體表面電場峰值對網格最敏感邊界層單元層數(shù)從3層提到6層峰值可以變化15%左右必須專門收斂。人體內部最大感應電流密度對網格的敏感度次之變化約5%到8%。地面以上1.5m處的電場和磁場對網格相對不敏感粗網格和細網格結果差不到3%。這給工程實踐帶來一個啟示不要只收斂一個代表性物理量要看最終寫進報告的每個量是否都收斂。我見過有同事只看了地面電場收斂就出報告結果人體內場其實還沒收斂。評判網格是否足夠的唯一標準是被評估量本身隨網格變化的幅度而不是整體云圖是否好看。6.3 全模型內存不夠時的兩階段法三維模型加精細人體網格后內存很容易沖破本機限制。兩個減少內存的常用手段一是幾何對稱性利用。如果人體站在電纜正上方且電纜沿軸向均勻可以只建四分之一模型在對稱面上加對稱邊界條件。工頻三相電流不是對稱的所以這里的對稱必須基于時均意義下的RMS場而不是瞬時場。二是背景場-子模型解耦法。第一階段先去掉人體只算包含電纜、土壤和空氣的背景場把結果保存在數(shù)據(jù)集里第二階段建一個人體加周圍局部空氣的小模型在子模型外邊界上施加第一階段的場值作為邊界條件只求解人體附近的擾動場。這樣模型規(guī)模從全局縮小到局部內存需求驟降。兩階段法的物理依據(jù)是人體和電纜之間距離至少一米以上人體對電纜電流回路的反作用可以忽略因此背景場不受人體存在影響。但使用時有條件子模型邊界必須離人體足夠遠通常至少2倍人體高度否則截斷誤差會污染體內場。我在模擬項目X里用這個方法把內存占用降低了近70%結果與全模型對比在1%以內的偏差工程上完全可用。6.4 與實測對比時的單位換算細節(jié)仿真報告最后總要面對現(xiàn)場測試的驗證環(huán)節(jié)。現(xiàn)場用場強儀測出來的讀數(shù)按定義是RMS有效值而COMSOL頻域云圖默認顯示的有時是幅值峰值。如果直接拿默認云圖去對比實測會高出根號2倍這個錯誤真的發(fā)生過很多次。我在后處理設置里把電場和磁場顯示都改成RMS同時在報告里注明所有仿真數(shù)據(jù)均為RMS有效值與現(xiàn)場場強儀讀數(shù)口徑一致。另外要提醒一點磁場測量探頭如果裝在人體附近人體會對低頻電場產生嚴重擾動但磁場幾乎不受人體影響所以現(xiàn)場測試時電場和磁場的可信度是不同級別的對比分析時要區(qū)別對待。結語如果只挑出一條最想分享的經驗那就是這類電磁場分析項目真正決定結果可信度的往往不是求解器本身夠不夠高級而是建模者對物理問題邊界的理解是否準確。三相回路要完整土壤參數(shù)要做掃描人體網格要單獨收斂限值對照要分清環(huán)境場和人體內場——每一條都是我在實際項目里真金白銀換來的教訓。COMSOL給你的是精確求解方程的能力但方程本身是不是對應你要回答的問題這件事永遠得自己把關。下次做類似的電纜電磁場評估時建議先從二維模型開始快速建立直覺再謹慎地移入三維精算你會發(fā)現(xiàn)整個過程順暢得多。