2016年 第23期Science and Technology & Innovation┃科技与创新
文章编号：2095－6835（2016）23－0013－03
不同模式下雷电与建筑物击距的对比分析
陈正龙，王 军
（江苏雷安检测技术有限公司，江苏 南京 210000）
摘 要：主要利用电气几何模型（EGM）、流注起始临界体积模型（CVM）和上行先导传输模型（SLIM）这三种方法，研究了雷击建筑物的引雷半径，并对三种模型的结果进行了对比分析。研究表明，在对引雷半径计算过程中，由于EGM模型忽略了上行先导的影响而具有一定的误差，且EGM的闪电距离随着高度的增加而增加，然而，对高度低于30 m的建筑物而言，使用EGM模型计算出的引雷半径误差低于20%；随着建筑物高度的增加，上行先导对引雷半径的影响越来越显著，CVM所预测的引雷半径比SLIM和EGM所预测的引雷半径大得多。 关键词：雷电；闪电击距；引雷半径；建筑物高度
中图分类号：TU5 文献标识码：A DOI：10.15913/j.cnki.kjycx.2016.23.013
28.3 m，该闪电击距结果要接近于Golde的计算结果，但是要远远地低于Whitehead等人的结果。
140 120由于雷电所带来的高峰值的雷电流和强大的电磁脉冲（LEMP）对各种物体，尤其是广泛存在的电子设备具有不可小觑的危害。因此，有必要研究闪电形成的机理。
目前国内外有很多学者对上行先导的始发条件——闪电击距进行了研究。Golde在工程实践中，首先提出了雷电流与闪电击距之间存在一定的关系。McEachron在90年代初期，对纽约帝国大厦闪电资料分析得出，当大厦周围的地面大气电场强度达到一定强度时，大厦会触发上行先导。Lee等采用EGM模型对闪电击距与回击电流大小之间的关系进行了研究。但是这个简单的EGM模型，没有考虑上行先导始发给闪击距离所带来的物理影响，同时也忽略了建筑物的高度对闪击距离的影响。
.com.cn. All Rights Reserved.魏本刚等根据高压输变电线路的结构和运行高压等情况，建立
了输变电线路表面连续上行先导起始的物理模型。清华大学科技工作者通过实际的雷电实验工程中，采用了分形理论研究上行先导的发展过程。
近年来，一些防雷公司将CVM模型运用在了建筑物防雷保护的研究工作中，但是在建筑物防雷保护的研究工作中，利用CVM方法得出结果的正确性有待考证。因此，本文主要利用电气几何模型、流注起始临界体积模型和上行先导传输模型这三种模型研究建筑物的引雷半径，对三种模型的结果进行对比分析。
1 闪电击距模型
1.1 电气几何模型（EGM）
由于闪电击距的表达式不同，所以在实际工程应用中也会出现一定的差异性。本文借助Matlab软件，研究了前人闪电击距表达式之间的差异性。由于Eriksson等人研究的闪电击距表达式考虑了建筑物的高度对其的影响，因此本文假定建筑物的高度为20 m，计算出上述表达式的闪电击距结果。图1为闪电击距随回击电流变化的趋势。从图中可以看出，在不同的回击电流作用下，不同的学者所提出的表达式计算出来的闪电击距是不一样。当回击电流大小为0～16 kA时，IEEE std计算出来的闪电击距最大，其次是Whitehead计算出来的闪电击距；当回击电流大于16 kA时，Whitehead计算出来的闪电击距最大，IEEE std计算出来的闪电击距其次。当回击电流大小在0～40 kA的研究范围内时，Eriksson计算出来的闪电击距始终最小，Golde计算出来的闪电击距结果与Eriksson较为接近。当回击电流大小为20 kA时，Whitehead计算出来的闪电击距为73.8 m；Golde计算出来的闪电击距为34.1 m，该闪电击距结果要低于Whitehead计算结果的2倍多；IEEE std计算出来的闪电击距为70.1 m，该闪电击距结果要接近于Whitehead计算结果；IEEE工作组计算出来的闪电击距为56.1 m，该闪电击距结果要低于Whitehead的计算结果17.7 m；Eriksson计算出来的闪电击距为
100闪电击距/mWhiteheadGoldeIEEE stdIEEE工作组Eriksson806040200 0510152025303540回击电流/kA
图1 闪电击距随回击电流变化趋势
由于Eriksson等人研究的闪电击距表达式考虑了建筑物的高度对其的影响，因此本文其次研究了不同建筑物高度下，闪电击距随回击电流大小的变化趋势。图2为不同高度下闪电击距随回击电流变化的趋势。从图中可以看出，在不同的高度下，利用Eriksson等人所提出的公式计算出的闪电击距具有一定差异性，说明了建筑物的高度对闪电击距具有很大的影响。总体上呈现出建筑物高度与闪电击距之间存在正相关性关系，即建筑物高度越高，闪电击距越大。当回击大小为20 kA时，建筑物高度为80 m时，计算出的闪电击距为65.1 m；当建筑物高度为50 m时，计算出的闪电击距为49.0 m，闪电击距较建筑物高为80 m时的闪击距离小16.1 m；当建筑物高度为30 m时，计算出的闪电击距为36.1 m，闪电击距较建筑物高为80 m时的闪击距离小29 m；当建筑物高度为20 m时，计算出的闪电击距为28.3 m，闪电击距较建筑物高为80 m时的闪击距离小36.8 m。
图2 不同高度下闪电击距随回击电流变化趋势
本文主要利用IEEE std 1234—1997中提出的电气几何模型（EGM）与下文的两种模型来研究建筑物的引雷半径以及连接过程。
1.2 流注起始临界体积模型（CVM）
目前在实际的工程应用中，对电极表面流注起始而需要的
·13·
科技与创新┃Science and Technology & Innovation
电场强度的判据表达式为：
2016年 第23期 1E0+=22.8[1+3](kV/cm). （1）
Regm=s for s≤h. （7）
式（6）（7）中：h是塔的高度；s是公式2～7所给出的闪电击距。
R式（1）中：R为电极的半径，取值为0.5～25 cm。
由于电极在正极性冲击电压下，所以电极头部流注的起始存在一定的时延。这主要是因为电极流注的产生需要同时满足两个条件：①电极在冲击电压的作用下，电极头部所感应出的电场强度要达到头部产生流注所需要的最小电场强度；②保证电极头部临界体积内有一定数量的自由电荷。根据此，Les提出了流注起始临界体积计算模型（CVM）。 1.3 上行先导传输模型（SLIM）
地面上建筑物产生稳定的上行先导需要3个过程：①地面建筑物表面产生流光放电；②流注向先导转化；③上行先导开始连续发展。
主要计算步骤如下：
第一步：提取梯级先导的时间或高度。在该分析中，我们假设雪崩前的阳离子数超过108时，那么电子雪崩将会转变为流光。
第二步：一旦满足了流光起始的条件，接地导体尖端物会产生一束流光。
第三步：研究这一爆发的流光是否能够产生上行先导。流光转化成上行先导，在流光所产生的电晕区内至少有1 μC的电荷量。如果流光区的电荷低于这一值，那么在之后的一小段时间内流光将不能转化为上行先导。需要注意的是，随着时间的增加，梯级先导产生的大气电场强度增加了，而流光爆发中的电荷也因此而增加了。
第四步：假设在时间t时，满足了流光向先导转化的必要条件，那么我们接下来的任务是估计所产生的上行先导部分的初始长度和半径。
第五步：在时间间隔Δt期间，背景电势中存在电荷，同时
还计算出了初始先导的长度L1，新触发的上行先导头部所产生
的流光区内电荷数量按照上述方法进行计算。
第六步：现在我们要考虑的是第n次的计算步骤，这其中有n个时间部分，且它们有各自对应的电势梯度和半径，第i个先导部分的半径和电势梯度来自于以下等式：π⋅a2Li
(t−Δt)=π⋅a2
(t)+γ
−1Liγ⋅PELi(t)⋅Δt. （2） E+Δt)=a2Li(t)
0Li(ta2t+Δt)
ELi(t). （3）
Li(式（2）（3）中：ELi（t）为电场强度；ILi（t）为时间t时先导Li的电流。通过这些，我们可以计算出电场强度随时间的变化趋势，同时先导通道的势能差值（给定时间）为：
k
ΔUL=∑ELi(t)⋅Li. （4）
i=1
计算出当先导发展到第i步时，上行先导头部的电位U（i）
tip
为：
U（i）=（tip
li）
E+xEEstrE−E∞−l(i)L∞
0∞ln[L
E−strEe/x0]. （5）
∞∞
式（5）中：（li）
L为仿真到第i步长时上行先导的长度；E∞为最终稳定上行先导的电位梯度；x0是一个常量，x0=v×θ，v为先导发展速度，θ为先导发展时间常数。 2 不同闪电击距模型计算比较
在上述分析中，电气几何模型（EGM）、上行先导自持传输模型（SLIM）和流注起始临界体积模型（CVM）所提出的闪电击距表达式仅适用于塔垂直高度（变化范围10～100 m）和回击电流（变化范围5～100 kA）。我们假设该塔的形状为圆柱形，具有半球形的尖端，该塔的半径假设为0.1 m。在EGM模型中，塔的引雷半径Regm可以利用以下表达式进行分析计算：
Regm=s2−(s−h)2for
s>h. （6）
·14·
在流注起始临界体积模型（CVM）中，利用等式2～7对闪电击距进行计算。本文模拟了塔高度分别为10 m、20m、50m、100m这4种情况下，利用上述三种模型分别对引雷半径进行计算。图3为三种模型下计算出的不同塔高度下引雷半径结果，图中的实线表示SLIM的预测值，“+”为是EGM的预测值，虚线是CVM的预测值。从图中可以看出，当塔的高度低于10 m时，SLIM的预测值与EGM预测值相一致，两条曲线重合；当塔高度大于10 m时，SLIM的预测值要稍大于EGM计算结果；在塔高10～100 m的范围内，CVM模型计算结果要远远大于SLIM、EGM计算的结果。
此外，从图中还可以看出，当塔高度的变化范围在10～50 m之间时，三种模型计算出的引雷半径的变化仅为30%.但是，从图中看出随着塔高度的增加，上行先导对引雷半径的影响越来越显著，因此SLIM模型计算出的引雷半径要比EGM预测的引雷半径大得多，主要原因是由于当塔高度增加时，上行先导变化的作用越来越大。我们发现在现如今各行各业中，利用CVM模型计算出的引雷半径值要比文中模拟出的结果要大。如果假设的临界半径低于38 cm，那么CVM和SLIM预测值之间的差异较大。在CVM模型中，如果上、下先导尖端满足一定条件时，上、下先导会发生连接，且此时CVM模型预测的引雷半径变得更大。如果忽略了上行先导的存在或其影响，那么CVM模型会预测出所有攻击地面的雷电，直到梯级先导直接达到接地导体为止。另一方面，在SLIM模型中，如果我们忽略了上行先导的存在，那么该模型将等效于EGM模型。
上述结果表明，上行先导对下行先导的发展起到了重要作用。在研究中还发现，当下行先导以一定的角度达到接地导体附近时，会发生侧击现象。
600 900 EGMEGMSLIM800SLIM500CVM700CVM600m400m//径径500半300半雷雷400引200引3002001001000 0020406080100 回击电流/kA 020406080100回击电流/kA（a）建筑物高10 m （b）建筑物高20 m
1200 1600 EGMEGM1400SLIM1000SLIMCVM1200CVMm800m1000//径径半600半800雷雷引引6004004002002000 0020406080100 0回击电流/kA 20406080100回击电流/kA
（c）建筑物高50 m （d）建筑物高100 m
图3 三种模型下计算出的不同塔高度下引雷半径结果
3 结论
本文主要利用电气几何模型（EGM）、流注起始临界体积模型（CVM）和上行先导传输模型（SLIM）这三种方法研究雷击建筑物的引雷半径，对三种模型的结果进行了对比分析。主要得出以下结论：①当塔的高度低于10 m时，SLIM的预测值与EGM没有明显差异，而CVM预测值远远大于其他两种模型计算的结果；当塔高度的变化范围在10～50 m之间时，三种模型计算出的引雷半径的变化仅为30%.②随着塔高度的增加，上行先导对引雷半径的影响越来越显著，且 SLIM模型计算出的引雷半径要比EGM预测的引雷半径大的多。③如果上、下先导尖端满足一定条件时，上、下先导会发生连接，且此时CVM模型预测的引雷半径变得更大。
.com.cn. All Rights Reserved.
2016年 第23期Science and Technology & Innovation┃科技与创新
文章编号：2095－6835（2016）23－0015－03
基于海量平台的水电厂生产实时系统建模办法
刘 浩，苏 刚，常玉红，郝国文，姜 涛，朱海峰
（1.北京许继电气有限公司，北京 100085；2.国网新源控股有限公司，北京 100053；
3.安徽响水涧抽水蓄能有限公司，安徽 芜湖 241080；4.华东桐柏抽水蓄能发电有限责任公司，浙江 台州 317200） 摘 要：基于海量平台的水电厂生产实时系统所需数据一般来源于下属电厂的监控系统、设备在线监测系统和部分离散实时系统。各厂监控厂家不一，对各自系统测点命名也没有统一的标准，数据集中到生产实时系统后，难免出现混乱的情况。同时，在这种无序状态下，接入测点数量是否与现场一致，数据是满足生产实时系统计算及统计需求，也很难管理。为此，需要一套行之有效的测点建模办法，规范建模工作，统一测点模型，使接入数据既能反馈现场实际情况，又能满足生产实时系统使用，同时方便测点管理。
关键词：海量平台；实时系统建模；描述规则；命名规范
中图分类号：TV736 文献标识码：A DOI：10.15913/j.cnki.kjycx.2016.23.015
表1 描述规则
序号1 2 3 4
层级第一层第二层第三层第四层
描述规则 举例 <源系统名><区域名> 监控系统.新源
<厂站名> 宜兴 <设备> 1号机组 <量测类型> 有功功率
说明 1
1
2
2
3
4
基于海量平台的水电厂生产实时系统需要接入多个系统不同量测类型、不同命名规则的测点信息，而随着接入规模的不断扩大，涉及多数据源问题、量纲统一问题、设备难以识辨问题等，这些测点模型的整理工作将变得非常繁重。另外，从应用角度考虑，对外需要提供统一标准的测点信息。因此，对测点信息的规范化整理工作变得异常重要。
通过制订规范化的测点模型，可以方便应用程序以一致的方式处理量测类型，最终实现测点的自动维护，以及测点和设备量测模型的自动映射。 1 概述
.com.cn. All Rights Reserved.本文所述的命名规范包含设备分层规范和测点命名规范2
部分内容，即测点标签名描述的规范。
在测点命名中，最好的方法是测点作为设备的动态属性，通过测点标签名自动识别，便于测点映射表的自动维护。但基于准实时库本身的特点，不具备分层结构的设备概念。而GAIA准实时库通过多层次的测点描述，间接建立层次性的设备概念，因此在建立GAIA准实时库的测点表时，可以按照SG-CIM的层次模型，对测点进行设备的分层描述，便于以后映射表的自动维护，同时通过分层的标签描述，按照物理及逻辑设备概念可方便进行测点的查询统计。 2 方法
通过测点的描述属性，添加测点的设备信息，实现“源系统.区域.厂站.设备.量测类型.值”，完成量测值和设备的挂接。 3 设备分层规范 3.1 综述
按照GAIA测点的设备属性及描述，将设备按照四层描述（可扩展）进行规范。具体规则如下。
参考文献
［1］R.H.Golde，The frequency of occurrence and the distribution of
lightning flashes to transmission lines［J］. Electrical Engineering，1945（）：902 910.
［2］McEachron.Lightning to the Empire State Building
［J］.Science，1939，227（4204）：149-217.
［3］Ralph H.Lee.Protection Zone for Buildings Against Lighning
Strokes Using Transmission Line Protection Practice［J］.Industry Applications IEEE Transactions on，1978，IA-14（6）：465-469. ［4］魏本刚，傅正财.雷电背景下特高压导线表面先导起始
的数值仿真［J］.中国电机工程学报，2009，29（28）： 6-12. ［5］He Jinliang，Tu Youping，Zeng Rong，et al.Numeral analysis
model for shielding failure of transmission line under
经过分层的设备描述后，可得到测点的标签描述全程。举例如下：
监控系统.新源.宜兴.1号机组.有功功率 机组在线监测系统.宜兴.1号机组.有功功率 3.2 具体说明
下面以各类设备具体说明。 3.2.1 机组类
表2为机组类描述规则。
表2 机组类描述规则
序号1 2 3 4
层级第一层第二层第三层第四层
描述规则
<源系统名><区域名>
<厂站名> <机组名称> <量测类型>
举例 监控系统.新源
宜兴 1号机组 有功功率
说明
标签描述全称：
监控系统.新源.宜兴.1号机组.有功功率 3.2.2 变压器类
表3为变压器类描述规则。
lightning stroke［J］.IEEE Trans on Power Delivery，2005，20（2）：815-822.
［6］Eriksson A J.The Incidence of Lightning Strikes to Power
Lines ［J］.IEEE Transactions on Power Delivery，1987，2（3）：859-870.
［7］Eriksson A J.An Improved Electrogeometric Model for
Transmission Line Shielding Analysis［J］.IEEE Transactions on Power Delivery，1987，PER-7（7）：67-68.
［8］A.M.Les，K.D.Srivastava.The lightning performance of unshielded
steel-structure transmission lines ［J］.IEEE Transactions on Power Apparatus and Systems，19（4）：437 445.
〔编辑：王霞〕
·15·
