雷电流数学模型仿真分析
雷电流数学模型仿真分析
摘 要:简要分析了雷电的产生机理和雷电的危害,比较选择出双指数函数、Heidler函数和脉冲函数三种能较好模拟雷电流变化的数学模型,并对其数学推导及计算过程进行了详细分析。在对数学模型详细的理论分析的基础上,结合实际中雷电的幅值,频率等参考值,对这三种雷电流数学模型进行了合理的参数选择。在MATLAB中分别建立了双指数函数、Heidler函数和脉冲函数三种雷电流模型,详细分析了三种模型产生雷电波形的波形变化趋势。
关键词:双指数函数;Heidler函数;雷电流 1.引言 自古以来,雷电就是一种令人生畏的自然现象。人们从很早就致力于对它的研究,但是目前,国内外关于雷电理论研究及计算方法尚欠完善,所以研究雷电的产生机理和对雷电的仿真分析已经成为对雷电研究的重点。为了更方便的对雷电进行仿真分析,用数学模型对其进行模拟已经得到了广泛??应用,但是人们对雷电流数学模型缺乏整体性的认识,在现有研究成果的基础上,对雷电流数学模型进行系统的比较研究,并对雷击输电线路产生的影响进行了分析。 2.雷电电流的数学模型 2.1 双指数函数 1941年,Bruce和Golde提出了双指数函数的数学表达式: ■ (1) 式中α是决定电流衰减的时间常数;β是决定电流上升的时间常数;α、β的大小可按照当时已知的知识,由闪电的三个特性推得。这三个特性就是:沿先导通道的电荷密度,回击速度,以及回击过程中先导电荷的复合率。 ■,为峰值电流的修正因子; tp=■,为峰值时间(通过求解■=0,可以得到tp值)。 把tp带入(4-1)式中可求得峰值电流为: ■(2) 假设半峰值时间为■,则有: ■(3) 从(1)、(2)和(3)式子中可以看出■不仅与■有关,而且与α、β有关。tp和■与α、β也有着复杂的数学关系(将(2)式带入式(3)中可求得■的数学表达式)。 2.2 霍德勒函数 Heidler提出了一种较新的雷电电流的数学模型: ■(4) ■(5) 其中■,■是决定电流上升(波头)的时间常数,n是电流陡度因子。 ■(6) 其中■是决定电流下降(波尾)的时间常数。 这个??学模型是由两个函数组成的,x(t)是电流上升时间函数,y(t)是电流下降时间函数,■是电流峰值。 假设 : 在电流下降时间内,x(t)≈1; 在电流上升时间内,y(t)≈1。 由式(5)可知当t >0时,x(t)<1,则电流的最大值小于■,因此,需要一个电流最大值的修正系数η。 国际电工委员会(IEC)在其1995年的文件IEC1312-1中,规定了供分析用的雷电流解析表达式(n=10时): ■(7) 式(7)是基于Heilder模型和传输线模型提出的,适用于首次雷击和后续雷击。 2.3 脉冲函数 雷电的放电回击过程可以由简单的天线模型来模拟,由这个模型得到关于雷电的电磁场表达式中的静电场项与回击电流的时间积分有关。因此,在进行雷电电磁场的计算中,都涉及到复杂的重积分运算,工作量巨大。双指数函数在t=0时有连续的一阶导数,而Heidler函数又没有明显的积分式。为了克服这种困难,有人提出了用式(8)形式的脉冲函数模型[2]来表示雷电电流 ■(8) 其中峰值修正因子■,■。 修正因子η与Heidler函数电流模型的修正因子η的求解方法相同,即由■以及■联立求得。 将式(8)中的■展开,可以得到[2]: ■(9) 有关文献认为:脉冲函数展开式(9)中的第一项(k=0时)是决定脉冲函数衰减的主要项。 3.三种雷电流数学模型仿真的对比分析 3.1函数的一阶可导性 为了方便对比分析,这里我们利用Matlab函数仿真功能,将三幅图仿真于同一坐标系里,分别给出了这三种雷电电流数学模型在t=0时,函数的一阶导数仿真波形。如图所示: 从图中可以看出这三条曲线中,只有双指数函数的一阶导数曲线在t=0时,di/dt达到最大值,而Heidler函数和脉冲函数在t=0时,di/dt都为0。因此,只有只有双指数函数在t=0时,没有连续的一阶导数 3.2三种函数仿真波形的特点 为了方便对函数波形对比分析,我们将利用Matlab仿真软件把用PSCAD/EMTDC的仿真出的三种函数仿真波形的全波、波头和半峰值仿真图放在同一个坐标系里,如图所示: 图2中的各函数其实是这样得到的:由于它们的衰减项主要决定于参数■,而脉冲函数和Heidler函数对应衰减项相同,对双指数函数令■,■。我们就可以看出它们的波形是比较接近的。特别是双指数函数和脉冲函数,在n值不大的情况下是非常接近的(图2(a)中可以得以体现)。实际上,在n=1的情况下脉冲函