陳博棟 逯 邁 陳小強 胡延文（蘭州交通大學(xué)自動化與電氣工程學(xué)院 蘭州 730070）

2（蘭州交通大學(xué)光電技術(shù)與智能控制教育部重點實驗室 蘭州 730070）

特高壓交流輸電線路對人體電場效應(yīng)的對比分析

陳博棟1逯 邁2陳小強1胡延文11（蘭州交通大學(xué)自動化與電氣工程學(xué)院 蘭州 730070）

2（蘭州交通大學(xué)光電技術(shù)與智能控制教育部重點實驗室 蘭州 730070）

本文通過有限元分析方法計算了1 000 kV輸電線路下方簡化人體模型的電場效應(yīng)，比較分析了不同塔型時人體中的感應(yīng)電場及其分布情況。結(jié)果表明：人體周圍電場并不是均勻分布的，而是產(chǎn)生了畸變；人體較細的部位（如頸部與腿部）電流密度較大；與普通型輸電線路相比，緊湊型線路下人體內(nèi)感應(yīng)電流密度較小；中樞神經(jīng)系統(tǒng)組織的最大感應(yīng)電流為遠小于2 mA·m-2安全限值；與緊湊型線路相比，普通型線路中樞神經(jīng)系統(tǒng)組織的電場強度較大，但其最大感應(yīng)電場仍小于20 mV·m-1的光幻視閾值。研究表明，特高壓交流輸電線路下人體感應(yīng)電流及電場均低于根據(jù)ICNIRP導(dǎo)則推薦的安全限值。

特高壓交流輸電，工頻電場，有限元法，感應(yīng)電流密度，感應(yīng)電場

低頻電磁場暴露生物效應(yīng)已成為生物電磁學(xué)研究熱點之一，如核磁共振梯度線圈電磁輻射[1-3]、手機及手機基站的電磁輻射等[4-6]。相關(guān)研究表明，長期暴露在強電磁場中，會發(fā)生中樞神經(jīng)系統(tǒng)機能障礙，出現(xiàn)頭暈頭痛、失眠多夢、記憶力衰退等癥狀，尤其是兒童白血病幾率增加[7-9]。

工頻電磁場與人作用可以在人體內(nèi)感應(yīng)出微弱的電場和電流。通過細胞和動物研究發(fā)現(xiàn)，人體內(nèi)感應(yīng)電流密度超過一定的閾值，感應(yīng)電流會對神經(jīng)和肌肉組織產(chǎn)生影響，這些影響隨著感應(yīng)電流密度增加加劇，為防止對神經(jīng)系統(tǒng)造成危害，暴露在電磁場中的人體感應(yīng)電流密度須限制在 100 mA·m-2以內(nèi)，取安全系數(shù)為10，人體職業(yè)暴露安全限值取10 mA·m-2，公眾暴露安全限值取2 mA·m-2[10]。根據(jù)ICNIRP 2010標(biāo)準(zhǔn)，即國際非電離輻射防護委員會(The International Commission for Non-Ionizing Radiation Protection, ICNIRP)發(fā)布的標(biāo)準(zhǔn)[11]，由于感應(yīng)電場是確定生物效應(yīng)的物理量，為避免感應(yīng)視網(wǎng)膜光幻視，以中央神經(jīng)系統(tǒng)組織中的感應(yīng)電場強度20 mV·m-1光幻視閾值作為公眾暴露限值基礎(chǔ)。

由于電能的需求的增加，輸電線的工作電壓也在不斷提高，特高壓輸電線路電壓都在1 000 kV及以上，這有利于緩解我國用電需求和線路走廊問題，更好地實現(xiàn)“西電東送”和“全國聯(lián)網(wǎng)”。輸電線路對環(huán)境和人體健康影響近年來倍受關(guān)注，輸電線路難免會經(jīng)過居民區(qū)，人也不可避免的會經(jīng)過輸電線路。因此，有必要對靠近人類生活輸電線的場以及與人體的相互作用進行研究。關(guān)于特高壓不同桿塔類型輸電線路在地面附近產(chǎn)生的電場和磁場進行了大量的分析[12-15]。但針對特高壓與人體相互作用的文獻較少，而且存在人體模型較為簡單，采用二維分析及邊界條件與實際不符等問題[16-18]。

本文通過研究不同塔型特高壓輸電線路下方人體與地絕緣時感應(yīng)電場對人體的影響，利用有限元法對人體感應(yīng)電流，尤其是中樞神經(jīng)的感應(yīng)電場進行計算，并與ICNIRP導(dǎo)則中的電磁暴露標(biāo)準(zhǔn)比較，對我國特高壓輸電線路的電磁暴露進行科學(xué)評估，可為生物醫(yī)學(xué)研究提供相關(guān)的參考數(shù)據(jù)，也可為國家制定相關(guān)標(biāo)準(zhǔn)提供依據(jù)。

1 分析原理與模型

1.1 分析原理

為進一步了解工頻磁場與人體健康的關(guān)系，首先必須精確計算人體各重要器官中產(chǎn)生的感應(yīng)電流及感應(yīng)電場。人體感應(yīng)電流的強度與分布不僅與外界場強有關(guān)，而且與人體的形狀、各種器官的導(dǎo)電率、介電常數(shù)等有關(guān)，因此精確計算存在一定的困難。采用有限元法解得電勢的場分布值，然后再經(jīng)過處理可得到電磁場中的其他物理量，如感應(yīng)電流密度、感應(yīng)電場等。

感應(yīng)效果的相對強弱與人體的電導(dǎo)率σ和介電常數(shù)ε有關(guān)。根據(jù)麥克斯韋方程組：

式中，H為磁場強度矢量，J為電流密度矢量，D為電通密度，E為電場強度矢量，B為磁感應(yīng)強度矢量，ρ為電荷體密度。

由電磁理論的麥克斯韋方程組構(gòu)成了3個本構(gòu)方程：

式中，μ為磁導(dǎo)率。為了便于數(shù)值求解，定義磁矢位 A和標(biāo)量電勢eΦ ，自動滿足法拉第電磁感應(yīng)電流和高斯磁通定律，定義如下：

電場可以寫成：

在靜態(tài)場中，如果所有的材料表現(xiàn)出介電性能的服從 σ＞＞ωε（歐姆電流占主導(dǎo)地位, ω為電流頻率），可以推出：

矢量勢A從電場中解耦等效于靜矢量勢A0，如果磁導(dǎo)率μ是恒定的，在整個計算區(qū)域?，由畢奧-薩伐定律可以計算出A0：

我們可以對方程(8)和(9)進行數(shù)值求解，采用有限元法解得磁勢和電勢的場分布值，然后再經(jīng)過處理可得到電磁場中的其他物理量。

1.2 分析模型

在計算特高壓輸電線路下方人體感應(yīng)電流時，首先在ANSYS中建立人體模型，各部分參數(shù)為：腳部位長0.17 m，寬0.06 m，高為0.03 m的長方體和兩個半徑為0.03 m的半圓柱組成；小腿由兩邊半徑分別為0.03 m和0.05 m，長0.42 m的椎體構(gòu)成；大腿由兩邊半徑分別為0.05 m和0.08 m，長0.43 m的椎體構(gòu)成；上身由中間長0.16 m，寬0.16 m，高0.6 m的長方體和兩邊半徑為0.08 m的半圓柱組成，內(nèi)部含一個長0.1 m，寬0.2 m，高0.25 m的內(nèi)臟，手臂由半徑0.04 m長0.64 m的圓柱構(gòu)成；頸部為半徑為0.06 m，高0.06 m的圓柱；頭部分別為半徑0.1 m、0.092 m、0.085 m的球體構(gòu)成頭皮、顱骨、和大腦的三層人頭模型，總體身高1.74 m。人體位于輸電線路正對中線位置，面向?qū)Ь€傳輸方向，如圖1所示。根據(jù)50 Hz下人體主要組織的介電常數(shù)和電導(dǎo)率[19]，假定人體各個部分由均勻介質(zhì)構(gòu)成，電導(dǎo)率及相對介電常數(shù)見表1。

圖1 人體及網(wǎng)格剖分模型Fig.1 Body and mesh generation model

表1 人體組織介電常數(shù)及電導(dǎo)率Table 1 Human tissues permittivity and conductivity

本文對當(dāng)前已運行某段1 000 kV特高壓線路模型進行分析，特高壓架空線路的普通塔型和新型緊湊型如圖2所示[17]，絕緣子串長均為11.5 m，線路運行額定線電壓為1 000 kV。

圖2 1000 kV特高壓普通型（左）和緊湊型（右）桿塔Fig.2 1000 kV UHV ordinary tower (left) and compact (right) tower

對于高壓輸電線下方人體感應(yīng)電場的計算尺寸較大， 而ANSYS可以很好地解決大型電磁仿真計算問題。在采用ANSYS三維諧性電磁場分析時，由于計算量較大，對計算機性能要求比較高，故忽略桿塔、端部效應(yīng)和導(dǎo)線弧垂的影響，線路長度取20 m，線路高度采用平均高度；選取10節(jié)點四面體實體單元SOLID 232求解三維模型的電勢分布；分別對各個實體材料進行賦值；在對模型進行網(wǎng)格剖分時采用Smartsize自動網(wǎng)格劃分工具，有利于生成形狀合理的網(wǎng)格單元，剖分后網(wǎng)格單元數(shù)量在一百多萬；對于模型的邊界條件，在導(dǎo)線上分別施加相電壓，相位差為120°，頻率為50 Hz，取地面電勢為零，其余邊界認(rèn)為是無窮遠點。由于網(wǎng)格數(shù)量較多，計算量較大，在進行求解時，計算機存儲容量需在8 GB以上，計算時間需1 h以上。

2 結(jié)果分析

2.1 線路下方人體周圍電場分布

如圖3和圖4所示，輸電線路下方人體周圍的電場，由于人體及空氣介電常數(shù)及電導(dǎo)率的變化，并不是均勻分布的，而是產(chǎn)生了畸變，致使人體周圍的場強比均勻場強增大了幾倍甚至好幾十倍。

圖3 普通塔型線路人體周圍電場分布Fig.3 Ordinary tower electric field distribution around the body

圖4 緊湊型線路人體周圍電場分布Fig.4 Compact tower electric field distribution around the body

由圖3和圖4可看出，在人體腳部和頭部周圍的電場強度較大，這是由于腳部和頭部集中了大量的異性電荷，從而產(chǎn)生了很大的場強；人體周圍最大場強出現(xiàn)在腳部位置，普通型塔桿線路與緊湊型塔桿線路下人體周圍最大場強相比緊湊型線路相對較大；但是人體頭部周圍的場強普通型桿塔線路相對較大，這可能與導(dǎo)線布置方式有關(guān)。

2.2 人體內(nèi)感應(yīng)電流密度

對于線路下方人體感應(yīng)電流的分布，如圖5所示。從總體上來看，自上而下人體內(nèi)電流密度呈逐漸增大趨勢；腳部位置的電流密度較大，而手臂位置感應(yīng)電流最??；人體較細的部位與較粗的部位相比電流密度較大，如頸部與腿部；普通型輸電線路與緊湊型線路下人體內(nèi)感應(yīng)電流密度相比緊湊型線路較小，分別為18 mA·m-2和8.6 mA·m-2，均位于腳部位置。

圖5 人體內(nèi)部感應(yīng)電流分布Fig.5 Induced current distribution inside the body

2.3 人體頭部感應(yīng)場分布

人體頭部的感應(yīng)電流和電場與周圍電場的大小以及各組織的材料有關(guān)，由圖6和圖7可知：由表及里，頭部的感應(yīng)電流呈減小趨勢，頭皮處的電流密度最大；普通型與緊湊型塔桿線路人體頭部感應(yīng)電流密度相比普通型線路較大，但均小于人體電流密度暴露限值2 mA·m-2。

由圖8和圖9可知，人體頭部電場，并非與電流密度分布一致，具體見表 2。電場強度最大值是在顱骨處；普通型線路下人體頭部最大場強為1.47 mV·m-1，低于ICNIRP導(dǎo)則推薦的20 mV·m-1的暴露限制，而緊湊型線路下人體頭部最大場強為0.32 mV·m-1遠低于這一限值，這與頭部周圍的場強以及線路布置方式有關(guān)，因此在單從公眾電磁暴露安全的角度建議采用緊湊型塔桿線路，從而對我國特高壓交流輸電電磁暴露安全評估提供理論依據(jù)。

圖6 普通塔型線路人體頭部電流密度Fig.6 Ordinary human head tower line current density

圖7 緊湊塔型線路人體頭部電流密度Fig.7 Compact tower line current density human head

圖8 普通塔型線路人體頭部電場分布Fig.8 Ordinary human head line tower electric field distribution

圖9 緊湊塔型線路人體頭部電場分布Fig.9 Compact tower line electric field distribution of human head

表2 人體頭部組織的最大感應(yīng)場強Table 2 Human head tissue maximum induction field distribution

3 結(jié)果與討論

本文分析了1 000 kV級特高壓交流輸電線路下不同塔型線路情況下的人體模型內(nèi)部感應(yīng)場的分布，得到以下結(jié)論：人體周圍的電場不是均勻分布的，而是產(chǎn)生了畸變，并且增大了幾倍甚至好幾十倍；從頭至腳，人體內(nèi)的電流密度呈增大趨勢，身體較細的部位，如頸部、小腿處的電流密度相對較大，手臂感應(yīng)電流最小。人體頭部（中樞神經(jīng)）位置的感應(yīng)電流與感應(yīng)電場分布并不一致，最大感應(yīng)電流位于頭皮處，而最大感應(yīng)電場位于顱骨處；普通塔型線路下中央神經(jīng)系統(tǒng)組織的感應(yīng)電流及感應(yīng)電場強度均略高于緊湊型線路，從公眾安全角度考慮，建議采用緊湊型塔桿線路，但人體中央神經(jīng)系統(tǒng)組織的最大感應(yīng)電流遠小于 2 mA·m-2的安全限值，感應(yīng)電場均小于20 mA·m-2的光幻視閾值。

人體簡化模型能夠很好地反映出人體內(nèi)部感應(yīng)場的變化，該結(jié)果將會推動生物醫(yī)學(xué)研究理論和應(yīng)用研究的進一步深化，為我國特高壓電磁暴露評估提供參考，但還存在不足之處，需進一步對真實人體模型的電場效應(yīng)進行研究。

1 胡格麗, 倪志鵬, 王秋良. 結(jié)合振動控制的柱面縱向梯度線圈目標(biāo)場設(shè)計方法 [J]. 物理學(xué)報, 2014, 63(1): 018301. HU Geli, NI Zhipeng, WANG Qiuliang. A target feld method for designing cylindrical z-gradient coil combined with vibration control [J]. Acta Physica Sinica, 2014, 63(1): 018301.

2 MAI Lu, Shoogo Ueno. Dosimetry of exposure of patients to pulsed gradient magnetic fields in MRI [J]. IEEE Transactions on Magnetics, 2011, 47(10): 3841-3844.

3 劉虎, 尤克增. 磁共振掃描中受檢者電磁輻射現(xiàn)狀與防護對策 [J]. 環(huán)境與職業(yè)醫(yī)學(xué), 2012, 29(9): 576-579. LIU Hu, YOU Kezeng. Electromagnetic radiation and protection in magnetic resonance scanning [J]. Journal of Environmental & Occupational Medicine, 2012, 29(9): 576-579.

4 高翔, 劉群芳, 姚勇. GSM 移動通信基站電磁輻射水平防護距離理論預(yù)測與現(xiàn)場實測對比分析 [J]. 輻射防護, 2011, 31(1): 46-49. GAO Xiang, LIU Qunfang, YAO Yong. Analysis of theoretical prediction and field measurement about GSM mobile base station electromagnetic radiation protection level [J]. Radiation Protection, 2011, 31(1): 46-49.

5 李旸, 逯貴禎. 手機輻射對DNA分子作用的機理研究[J]. 武漢大學(xué)學(xué)報(理學(xué)版), 2009, 55(4): 460-464. LI Yang, LU Guizhen. Research of mechanism of interaction between the microwaves radiated from the mobile phone and DNA molecules [J]. Journal of Wuhan University: Natural Science Edition, 2009, 55(4): 460-464.

6 MAI Lu, Shoogo Ueno. Comparison of specific absorption rate induced in brain tissues of a child and an adult using mobile phone [J]. Journal of Applied Physics, 2012, 111(7): 07B311.

7 Attwell D. Interaction of low frequency electric fields with the nervous system: the retina as a model system [J]. Radiation Protection Dosimetry, 2003, 106(4): 341-348.

8 Feychting M, Floderus B, Ahlbom A. Parental occupational exposure to magnetic fields and childhood cancer (Sweden) [J]. Cancer Causes Control, 2000, 11(2): 151-156.

9 Davis S, Mirick D K, Stevens R G. Residential magnetic fields and the risk of breast cancer [J]. American Journal of Epidemiology, 2002, 155(5): 446-454.

10 International Commission on Non-Ionizing Radiation Protection. Guidelines for limiting exposure to time-varying electric, magnetic and electromagnetic fields (up to 300GHz) [J]. Health Physics, 1998, 74(4): 494-522.

11 ICNIRP Guidelines for limiting exposure to time-varying electric and magnetic fields (1Hz to 100 kHz) [J]. Health Physics, 2010, 99(6): 818-836.

12 彭迎, 阮江軍. 模擬電荷法計算特高壓架空線路3維工頻電場 [J]. 高電壓技術(shù). 2006, 32(12): 69-74. PENG Ying, RUAN Jiangjun. Calculation of three-dimensional harmonic electric field around ultra high voltage overhead line based on the charge simulation method [J]. High Voltage Engineering, 2006, 32(12): 69-74.

13 肖冬萍, 何為, 張占龍, 等. 特高壓輸電線工頻磁場三維優(yōu)化模型 [J]. 中國電機工程學(xué)報, 2009, 29(12): 116-120. XIAO Dongping, HE Wei, ZHANG Zhanlong, et al. Three-dimension optimization model for power frequency magnetic field of UHV transmission lines [J]. Proceedings of the CSEE, 2009, 29(12): 116-120.

14 蘭生, 張振興, 原永濱. 考慮弧垂的交流特高壓輸電線三維電磁場 [J]. 電機與控制學(xué)報, 2012, 16(12): 42-46. LAN Sheng, ZHANG Zhenxing, YUAN Yongbin. The 3D electromagnetic field of AC EHV transmission lines with sag [J]. Electric Machines and Control, 2012, 16(12): 42-46.

15 黃韜, 呂建紅, 文遠芳, 等. 基于矩量法的輸電線路工頻電磁場實測與計算分析 [J]. 高壓電器. 2013, 49(2): 30-36. HUANG Tao, LYU Jianhong, WEN Yuanfang, et al. Calculation of power frequency electromagnetic field around transmission line based on moment method [J]. High Voltage Apparatus, 2013, 49(2): 30-36.

16 王建華, 文武, 阮江軍. 特高壓交流輸電線路工頻磁場在人體內(nèi)的感應(yīng)電流密度計算分析 [J]. 電網(wǎng)技術(shù). 2007, 31(13): 7-10. WANG Jianhua, WEN Wu, RUAN Jiangjun. Calculation and analysis on power frequency magnetic field induced current density within human body under uhvac transmission line [J]. Power System Technology, 2007, 31(13): 7-10.

17 胡宇. 超高壓輸電線環(huán)境中人體電磁場分析 [D]. 沈陽.沈陽工業(yè)大學(xué). 2003. HU Yu. The electromagnetic field analysis in human body in the environment of the extra-high voltage transmission line [D]. Shenyang. Shenyang University of Technology. 2003.

18 田子山. 交流架空輸電線路附近工頻電場及其人體內(nèi)感應(yīng)電流計算研究 [D]. 重慶. 重慶大學(xué). 2013. TIAN Zishan. Calculation of power frequency electric field from ac overhead transmission lines and induced current in human body [D]. Chongqing. Chongqing University. 2013.

19 Gabriel C, Gabriel S, Corthout E. The dielectric properties of biological tissues: I. literature survey [J]. Physics in Medicine and Biology, 1996, 41(11): 2231.

Comparative analysis of electric field effects on the human body induced by ultra high voltage AC transmission line

CHEN Bodong1LU Mai2CHEN Xiaoqiang1HU Yanwen11(College of Automation and Electrical Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China)
2(Key Laboratory of Opto-Electronic Technology and Intelligent Control, Lanzhou Jiaotong University, Lanzhou 730070, China)

To study the influence of ultra high voltage (UHV) AC transmission line frequency electric field on the human body, finite element analysis method was used to calculate electric field effect of the simplified human body model under 1000 kV transmission line. Then, comparative analysis of different towers was carried out to show the differences in the induced electric field and its distribution in the human body. The results indicated that electric field around the body is not evenly distributed, and distortion was produced; current density in the small parts of human body (e.g. the neck and legs) is high; compared with compact transmission lines, normal transmission lines induce smaller current density in the human body; the maximum induced current in the central nervous system is much smaller than safety limits of 2 mA·m-2; electric filed strength in the central nervous system under normal transmission lines is high, but the maximum induced electric field is still less than 20 mV·m-1, which is the phosphene threshold value. Studies demonstrate that induced current and electric field in the human body under UHV AC transmission line are both less than the reference levels from the International Commission on Non-Ionizing Radiation Protection (ICNIRP) guidelines. The results can not only provide the basis for calculation of electromagnetic exposure of the real human body, but also provide reference for electromagnetic exposure assessment standard of UHV power transmission engineering in China.

Ultra high Voltage (UHV) AC transmission, Frequency electric field, Finite element method, Induced current density

TM81，TL72DOI: 10.11889/j.1000-3436.2014.rrj.32.060501

國家自然科學(xué)基金(51267010)和甘肅省杰出青年基金(1308RJDA013)資助

陳博棟，男，1988年2月出生，2010年畢業(yè)于蘭州交通大學(xué)，現(xiàn)為該校在讀碩士研究生，研究方向為電力系統(tǒng)及其自動化

逯邁，博士，教授，博士生導(dǎo)師，E-mail: mai.lu@126.com

初稿2014-07-23；修回2014-08-28

CLCTM81, TL72