雷击风险评估参数优化探讨.pdfVIP

雷击风险评估参数优化探讨 1 1 1 1 2 1 1 李 霞 樊 荣 肖稳安 王振会 曾山泊 植耀玲 徐鸣一 （1.江苏省南京信息工程大学 南京 210044；2.江苏省苏州市气象局 苏州 215021 ） 摘 要 雷击风险评估参数的取值决定了评估结果的精确性。本文依据GB/T 21714.2对雷击风险评估参数进行 了优化和探讨：（1）利用闪电定位仪监测资料，提出雷击大地密度（Ng）的计算方法；（2）根据国标中的相关数据，研 究提出雷击建筑物附近引起内部系统故障概率（Pms）的既快又准的计算方法；（3）研究提出多线路系统损失分量（Lo） 的细化方法。通过这些参数的优化，达到优化风险值的目的。 关键词 ： 雷击风险 评估参数 计算方法 引 言 雷击灾害是一种频发性的自然灾害，其危害随着社会信息化和电子化的发展变得更加易发、明显和 广泛。开展雷击风险评估，是防止和减少雷击灾害损失的重要手段，也是进行科学合理雷电防护的必要 前提。目前，我国开展雷击风险评估工作依据的标准为《GB/T 21714.2-2008/62305-2：2006，雷电防 护第2部分：风险管理》[1]，其评估结果的准确性取决于各计算参数的精确取值。而在规范中，很多参 数的取值不够具体，特别是针对某个区域、某个建筑时，无法给出精确的参数值，这就严重影响了计算 结果的精准性。因此，本文依据《GB/T 21714.2-2008/62305-2：2006，雷电防护第2部分：风险管理》， 对其参数进行研究和分析，探究雷击风险参数的计算方法，以便更为准确、真实地反映出被评估建筑物 的真实情况。 1 雷击大地密度（Ng）计算方法的优化 1.1 利用人工观测资料计算Ng的局限性 雷击大地密度（Ng）是进行雷击风险评估的重要参数之一。雷击风险评估业务的传统方法中，通 常采用当地气象部门提供的所在城市的年平均雷暴日(Td)来计算评估点的Ng值。雷暴日Td值一般是由 该城市全天候的观测站给出的，或者是该城市多个观测站（不同观测站之间的距离一般有几十公里）观 测数据的平均。一个城市往往以当地气象部门发布的Td为标准[2]。得到了Td值，然后用GB/T 21714.2 提供的Ng的计算公式：Ng=0.1Td，计算Ng。使用这种方法计算Ng非常明显的缺陷是同一个行政辖区 内的Ng值完全相同，这就在很大程度上忽略了一个城市不同地区因地理、气象、建筑结构布局对雷电 发生条件的不同而引起的雷击次数不同的影响，与实际的Ng分布情况通常有较大的差异。 1.2 利用闪电定位仪观测资料计算Ng 根据雷击大地年平均密度的含义（即一年里某一确定地区的地闪次数），通过雷电定位系统提供的 一年中某地区的地闪次数，再除以该地区的面积，可以得到更为精确的雷击大地年平均密度[3]。计算公 式为： Ng=D/S D——某地区一年中的地闪次数（次/a） 2 S——该地区的面积（km） 根据目前的技术水平和条件，D和S都可以得到较为精确的数值，这样就可以计算出精确的Ng值。 1.3 软件算法 为了在计算机上实现Ng的快捷计算，需要将上述方法转变成软件算法。如图1所示，黑色格子的 中心位置为评估点所在经纬度。每个格子边长为3km，每边共有9个格子。算法具体内容如下： 气象与减灾 2011专刊 (1) 在程序中声明四个二维数组：① ArrayNg；②ArrayAvgLigCurt；③ ArrayBoxsLng；④ArrayBoxsLat。分别用来存 放图1中每个格子的①Ng值；②雷电流平均 强度；③格子中心点经度；④格子中心点纬度。 (2) 已知用户输入的评估点经纬度。而且 根据地球上任意两点的经纬度又可以得到两 点间距离。计算公式如下： d(x1,y1,x2,y2)=r*arccos(sin(x1)*sin (x2)+cos(x1)*cos(x2)*cos(y1-y2)) 其中：x1,y1是纬度、经度的弧度单位， r为地球半径。 图1评估点附近网格分布图 据此，可以经过计算得出每个格子的经纬度范围，ArrayBox