基于函数波形的雷电流建模方法对比研究.docxVIP

基于函数波形的雷电流建模方法对比研究 1 电化学流数学模型的分析 随着能源系统的发展，现代能源系统的闪燃和开路事故变得越来越重要。雷电电流波的数学模型是研究雷电的主要内容之一,人们对雷电流数学模型缺乏整体性的认识。因此,在现有研究成果的基础上,对仍应用广泛的双指数函数、Heidler模型和脉冲函数模型进行分析比较愈加重要。 本文从函数波形特点、函数的一阶可导性、频谱分析等方面进行系统的比较研究,比较其优缺点,以得知其各自适用的研究领域。 2 三个雷流数学模型 2.1 热电流波形的确定 1941年,Bruce和Golde提出了双指数函数的数学表达式: i(0,t)={0,t0Ι0η[exp(-αt)-exp(-βt)],t≥0(1)i(0,t)={0,t0I0η[exp(?αt)?exp(?βt)],t≥0(1) 式中,I0—峰值电流;α—决定电流衰减的时间常数;β—决定电流上升的时间常数;η—峰值电流修正因子。 根据建筑物防雷设计规范国家标准GB50057-94选取10/350μs的雷电流波形,τ1和τ2是波头和波尾时间常数,取α=1/τ2,β=1/τ1,I0=20kA。exp(-αt)决定了函数图形的衰减部分,exp(-βt)决定其上升沿。 2.2 电流—Heidler函数模型 Heidler提出了一种较新的雷电电流的数学模型: {i0(t)=Ι0x(t)y(t)x(t)=kns/(1+kns)y(t)=exp(-t/τ2)(2)?????i0(t)=I0x(t)y(t)x(t)=kns/(1+kns)y(t)=exp(?t/τ2)(2) 这个数学模型是由两个函数组成的,x(t)—电流上升时间函数;y(t)—电流下降时间函数;I0—峰值电流;ks=t/τ1,τ1—决定电流上升(波头)的时间常数;n—电流陡度因子;τ2—决定电流下降(波尾)的时间常数。 基于Heilder模型和传输线模型,国际电工委员会(IEC)在其1995年的文件IEC1312-1中,规定了供分析用的雷电流解析表达式(n=10时): i0(t)=Ι0η[kns/(1+kns)]exp(-t/τ2)(3)i0(t)=I0η[kns/(1+kns)]exp(?t/τ2)(3) 式中,η—峰值电流修正因子。 2.3 各自适应的首效应值计算 文献用式(4)形式的脉冲函数模型来表示雷电电流。 i(0,t)={0,t0Ι0η[1-exp(-t/τ1)]nexp(-t-τ2),t≥0(4)i(0,t)={0,t0I0η[1?exp(?t/τ1)]nexp(?t?τ2),t≥0(4) 峰值修正因子η与Heidler函数电流模型的修正因子η的求解方法相同,即由di(t′)/dt=0以及i(t=t′)=I0联立求得。 η=(1-tα)ntτ1/τ2ατ1/τ2α,tα=τ1/(τ1+nτ2)。 将式(4)中的[1-exp(-t/τ1)]n展开,可得: [1-exp(-t/τ1)]n=n∑k=0(-1)kn!k!(n-k)!exp(-kt/τ1)(5) 为了计算和仿真的方便,取n=2,I0=20kA,根据建筑物防雷设计规范GB50057-94选取10/350μs的雷电流。 3 三个雷达流数学模型的比较分析 3.1 放电流波形的上升沿 为方便对函数波形对比分析,利用Matlab仿真软件将三种仿真波形的全波、波头和半峰值仿真图放在同一个坐标系里,如图1所示。 由图1可以看出,对于雷电流波形的上升沿,Heidler函数和脉冲函数相对与双指数函数有明显的改善。双指数函数和脉冲函数雷电流上升速度快,很快就达到雷电流的最大值;而Heidler函数相对与其他两个函数,雷电电流波形的上升前沿比较陡,到达最大值后,经过较长的一段时间才降到最小值,这种情况与Gerger等人的实际测量结果是一致的。因此,选用Heidler函数式(3)作为雷电流随时间变化关系比选用其它的更符合雷电电流的实际规律。 3.2 模型一阶导数仿真 为方便对比分析,利用Matlab函数仿真功能,将三幅图仿真于同一坐标系里,分别给出了这三种雷电电流数学模型在t=0时,函数的一阶导数仿真波形,如图2所示。 由图中可以看出这三条曲线中,只有双指数函数的一阶导数曲线在t=0时,di/dt达到最大值,而Heidler函数和脉冲函数在t=0时,di/dt都为0。因此,只有双指数函数在t=0时,没有连续的一阶导数。 3.3 操作在使用一定的频率 基于雷电防护的需要,对于一些精密设备进行雷电防护的主要手段就是利用滤波的方法。而要对雷电流波形信号进行滤波处理,首先应该知道雷电流波的频谱,这样才能进行有效的防护。所以,有必要对雷电流波进行频谱分析如图3～图5所示。 由图可知,雷电流所包含的频率从低频到高频都有,但主要集中在低频部分。对于双指数形式