雷电流

雷电电流数学模型的分析与研究

（合肥工业大学电气工程与自动化学院，安徽 合肥 230009）

摘要：本文选取了常见的几个雷电电流的数学模型进行对比分析，运用Matlab 数学软件，分别做出它们的雷电电流波形图，在此基础上运用傅里叶变换分别对几种模型雷电流波的频谱进行分析，并做出了雷电流的频率分布图，为进一步深入研究雷电电磁场的计算和雷电电磁脉冲的防护提供理论基础。

关键词：雷电流 傅里叶变换 频率分布 频谱分析

Research and Analysis of the mathematic models about

the lightning current

Abstract: In this paper, it takes several familiar mathematic models about the lightning current.First of all, it can make their waveform diagrams about the lightning current using the Matlab mathematic software. On this foundation, it uses the Fourier transform to analyzed the mathematic models in frequency domain and makes their frequency distribute diagrams. Thus, it can analyze the value of the lightning current frequency diagram and offer theoretical references to further research the induction of LEMP and .the count of electromagnetism field.

Key Words : Lightning current, Fourier transform, Frequency distribute, Frequency chart analysis 0 引言

雷电电流波的数学模型是研究雷电的主要内容之一，因为一旦知道雷电电流波形，就可得到有关雷电流的参数，如雷电流的峰值，最大电流的上升率，峰值时间等，此外，通过推出的雷电流的数学表达式，将为雷电过电压保护、雷电电磁脉冲(LEMP)防护和雷电电磁场

计算等提供根本的理论基础[1]

。

由于雷电的产生具有很大的随机性，且与地质结构，土壤的电阻率等都有很大的关系，

所以没有一幅电流波形是相同的。但是人们通过长期的观测，统计出雷电流的几个特征[2]

： 1) 峰值电流：典型值为A 4102⨯左右，变化范围为3102⨯～A 5102⨯。 2) 电流上升率：典型值为1410-•s A μ左右，变化范围为310～14108-•⨯s A μ。 3) 峰值时间：典型值为2s μ左右，变化范围为1～30s μ。

4) 半峰值时间(电流随时间衰减到峰值50 %的时间)：典型值为40s μ左右，变化范围为

10 ～250s μ。 1 雷电流的解析表达式

从人们开始研究雷电流至今，已经提出了很多雷电流的数学模型，有幂级数模型，双指数函数模型，戈尔德模型，霍德勒(Heidler)函数模型，脉冲函数模型等等，本文选取较常见的经典的双指数函数模型、Heidler 函数模型和脉冲函数模型进行分析和比较。 1.1双指数函数模型

1941年，Bruce 和Golde 提出了雷电流波形的双指数函数表达式[3]

：

()()t t e e I

t i βαη

---=0 (t ≥0) （1）

上式中α为雷电流波头衰减系数；β为雷电流波尾衰减系数；p

p

t t e

e

βαη---=为峰值修正

因子； )/()/ln(αβαβ-=p t 为峰值时间。 由(1)式可知：dt t di /)(在t=0时为无穷大。

设峰值时间为m T ，峰值为m I ，半峰值时间为h T ，对(1)式两边求导并令为零：

0)(0=-=--t t e e I dt dI αβαβ 得：α

βαβln 1-=m T 将m T 代入（1）式得： )(ln

ln

0α

βα

ββα

β

α

βα--

--

-=e

e I I m

半峰值为时间为：

)(2

1

0h h T T m e e I I βα---= 从上式可以看出如果α，β和0I 已知，就可以做出雷电流波形图，通过国军标GJB2639-96

[4]

中查得不同α，β和0I ，用Matlab [5]

仿真出不同α，β和0I 时的雷电流的全波(图1)和波头(图2)的波形图：

图1 全波波形图 图2 波头波形图 其中红线表示的参数是：A I 2188100=,111354-=s α,1647265-=s β 绿线表示的参数是：A I 1094050=,122708-=s α,11294530-=s β.

1.2 Heidler 函数模型

国际电工委员会IEC 于1995年发表的文件IEC1312-1[6]中，规定了供分析用的雷电流解析表达式，即为Heidler 函数模型：

()()()

2

110

1τ

ττη

t n

n

e t t I t i -+⋅=

(t ≥0) （2）

式中：0I —峰值电流，η—峰值电流修正因子，1τ—波头时间常数，2τ—波尾时间常数； n 是电流陡度因子，一般情况下取n =10。

(2)式是基于霍德勒模型和传输线模型提出的，适应于首次雷击(10/350s μ)和后续雷击(0.25/100s μ)。通过IEC1312-1中规定的参数，用Matlab 仿真出首次雷击和后续雷击时的雷电流的全波(图3)和波头的波形图(图4)：

其中红线表示首次雷击(10/350s μ)防护级别1时的参数：KA I 2000=，930.0=η，

s μτ101=，s μτ3502=；蓝线表示首次雷击防护级别2时的参数：KA I 1500=，930.0=η，