DB43_T 2133-2021 雷电致灾风险区划技术规范.docx

《DB43_T 2133-2021 雷电致灾风险区划技术规范.docx》由会员分享，可在线阅读，更多相关《DB43_T 2133-2021 雷电致灾风险区划技术规范.docx（18页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、ICS07.060CCSA4743湖南省地方标准DB43/T21332021雷电致灾风险区划技术规范Technicalspecificationforlightningdisasterriskzoning2021-08-03发布2021-10-03实施湖南省市场监督管理局发布DB43/T21332021目次前言1范围12规范性引用文件13术语和定义14区划流程25资料收集26资料处理37区划模型和方法4附录A（资料性）归一化方法7附录B（资料性）雷电落点误差计算8附录C（资料性）权重大小确定方法10附录D（资料性）自然断点法13参考文献14IDB43/T21332021前言本文件按照GB/T1

2、.12020标准化工作导则第1部分：标准化文件的结构和起草规则的规定起草。请注意本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。本文件由湖南省气象局提出。本文件由湖南省气象标准化技术委员会归口。本文件起草单位：长沙市气象局、怀化市气象局、湖南省气象灾害防御技术中心。本文件主要起草人：钟金莲、薛志明、邓战满、黄柱坚、薛明、刘越屿、周璐、薛德锋、唐瑶、陈梦醒、杨加艳、黄钰杰、谢露、黄正才、吴琏、吴瑕、卢荣琼、邱庆栋、陈力强、汤光玉。IIIDB43/T21332021雷电致灾风险区划技术规范1范围本文件给出了雷电致灾风险区划的流程、资料收集与处理、区划模型和方法。本文件适用于雷电

3、致灾风险区划。2规范性引用文件下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用文件，仅该日期对应的版本适用于本文件；不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。GB/T370472018基于雷电定位系统（LLS）的地闪密度总则QX/T4052017雷电灾害风险区划技术指南DB43/T21342021地闪等级划分3术语和定义下列术语和定义适用于本文件。3.1雷电致灾风险Riskoflightningdisaster雷电灾害发生的可能性。3.2致灾因子hazard可能造成人员伤亡、财产损失、资源与环境破坏、社会系统混乱等的异变因子。QX/T40

4、52017，定义3.23.3雷电致灾风险指数lightningdisasterriskindex根据致灾因子危险性对雷电致灾风险进行评定的量化指标。3.4雷电致灾风险区划lightningdisasterriskzoning根据雷电致灾风险指数大小，对雷电致灾风险的空间范围进行区域划分。3.5地闪cloud-to-groundlightning；CG雷暴云与大地之间的放电现象。GB/T370472018,定义3.1.13.6雷电定位系统lightninglocationsystem；LLS由多个雷电传感器组成的用于监测和定位其覆盖区域内雷电事件的探测网。1注：单位为次每平方千米年次/(kma)

5、。DB43/T21332021雷电定位系统又称闪电定位系统。GB/T370472018,定义3.1.103.7地闪密度groundflashdensity单位面积、单位时间的平均地闪次数。2GB/T370472018,定义3.1.73.8地闪日数lightningdays网格内年平均地闪日数。注：单位为天每网格年。3.9雷电强度lightningIntensity按地闪等级划分中的地闪强度等级加权平均得到的强度。注：单位为千安培（kA）。4区划流程雷电致灾风险区划流程，见图1。雷电资料灾情资料地理信息资料资料收集归一化处理归一化处理栅格化处理GIS数据库雷电致灾风险区划模型雷电致灾风险区划结果

6、图1雷电致灾风险区划流程资料处理风险区划5资料收集5.1雷电资料2等级等级名称地闪强度（kA）条件关系地闪陡度（kA/S）1级弱雷电0，15)或52级中等雷电15，25)且53级强雷电25，50)且54级超强雷电50且5DB43/T21332021宜收集5年以上的雷电定位系统资料，包括雷击的时间、地理位置、雷电流幅值、陡度值等资料。5.2地理信息资料按QX/T4052017规定，收集分辨率不低于1250000的数字高程模型（DEM）数据。5.3灾情资料宜收集5年以上雷灾起数。6资料处理6.1处理方法6.1.1归一化处理参照QX/T4052017给出的方法，对收集的资料进行归一化处理，归一化处理

7、方法参见附录A。6.1.2栅格化处理按2000国家大地坐标参数，将区域划分为5km5km的网格。6.2地闪资料6.2.1落点误差处理对收集的雷电资料进行落点误差计算，剔除雷电最大落点误差大于2倍网格边长的数据，雷电落点误差计算参见附录B。6.2.2陡度值处理雷电定位系统数据中，当陡度值大于雷电流幅值时，陡度值取雷电流幅值。6.2.3地闪强度等级划分根据雷电流幅值、陡度值大小，参照DB43/T21342021给出的方法，按表1将地闪强度划分为4个等级。表1地闪强度等级6.2.4地闪日数统计各网格内地闪强度2级（含）以上地闪日数，除以资料年限，得到各网格内的年地闪日数，并进行归一化处理，形成地闪日

8、数栅格数据。6.2.5地闪密度3雷灾起数（起）1-1011-2021雷电致灾等级（级）123iLn=(Fi)（1）DB43/T21332021统计各网格内地闪强度2级（含）以上地闪频次，除以资料年限，得到各网格内的年地闪密度，并进行归一化处理，形成地闪密度栅格数据。6.2.6雷电强度根据表1，按式（1）计算各网格内的地闪强度2级（含）以上雷电强度并形成栅格数据。4i=29式中：Ln雷电强度；i地闪强度等级；Fi地闪强度为i等级的地闪频次的归一化值。6.3地理信息资料6.3.1土壤电导率对土壤电导率资料进行归一化处理，形成归一化的土壤电导率栅格数据。6.3.2海拨高度对数字高程模型（DEM）资料

9、进行归一化处理，形成归一化的海拨高度栅格数据。6.3.3地形起伏计算以目标栅格为中心、大小为所定栅格大小的正方形范围内高程的标准差，并进行归一化处理，形成归一化的地形起伏栅格数据。6.4致灾等级参照雷灾起数，按表2统计栅格面积上雷电致灾等级，并进行归一化处理,形成归一化的雷电致灾等级栅格数据。表2雷电致灾等级表6.5建立数据库将地闪资料、灾情资料和地理信息资料处理成相同空间分辨率和空间投影坐标系统，建立地理信息系统（GIS）风险区划数据库。7区划模型和方法7.1区划模型雷电致灾风险区划模型由地闪因子和触发因子数据组成。区划模型见图2。4雷电致灾风险区划模型雷电致灾风险指数RHDB43/T213

10、32021雷电致灾风险等级划分地闪因子触发因子极高高风一般地闪日数地闪密度雷电强度致灾等级土壤电导海拨高度地形起伏风险等级险等级风险等级LrLdLnCl率EsTrSc图2雷电致灾风险区划模型7.2区划方法7.2.1雷电致灾风险指数计算雷电致灾因子危险性按式（2）计算：chRH=(Lwr+Lwd+Lwn)(Clwc+Sws+Ewe+Trwt)（2）式中：RH雷电致灾因子危险性；Lr地闪日数；wr地闪日数权重；Ld地闪密度；wd地闪密度权重；Ln雷电强度；wn雷电强度权重；Cl致灾等级；wc致灾等级权重；Sc土壤电导率；ws土壤电导率权重；Eh海拨高度；we海拨高度权重；Tr地形起伏；wt地形起伏

11、权重。各权重的计算参见附录C。5DB43/T213320217.2.2雷电致灾风险等级划分依据雷电致灾风险大小，采用自然断点法，将雷电致灾风险等级划分为极高风险等级（级）、高风险等级（级）、和一般风险等级（级）。自然断点法参见附录D。6DB43/T21332021附录A（资料性）归一化方法归一化是将有量纲的数值经过变换，化为无量纲的数值，进而消除各指标的量纲差异。计算公式为：Dij=0.5+0.5Aij-minimaxi-mini（A.1）式中：DijAijmini、maxij站（格）点第i个指标的归一化值；j站（格）点第i个指标值；分别是第i个指标值中的最小值和最大值。7DB43/T2133

12、2021附录B（资料性）雷电落点误差计算B.1前言雷电落点探测误差直接影响雷电灾害的定位确定。雷电位置是根据定位方式原理确定的，雷电定位方式包含时差定位和磁方向定位。雷电定位方式包含二站振幅、二站混合、三站混合、四站算法四类，不同的定位方式可产生不同的定位误差。通过对湖南雷电数据的落点分析，雷电落点误差主要来自于二站振幅。二站混合虽有一定误差，但比二站振幅小很多，且影响范围较小，基本可以接受。三站混合或四站算法误差极小，可以忽略。故本文件中只对二站振幅落点进行误差计算。B.2原理简介磁定向法（MDF）是多站定位法中的一种。MDF是利用两个相互正交且垂直的环（平面分别指向东西向和南北向）测量来自

13、近似垂直的闪电回击通道的辐射磁场，每个环上的感应电压与回击磁场的强度和入射方向有关，通过比较两个环上的感应电压，来确定闪电辐射源的方向。某次闪电被二个测站同时探测到，则二个测站磁方向的交点就是闪电落点。事实上，磁方向并不是一条线，而是存在一个误差角（一般小于1），两个测站M、N以误差角方向与地面的交叉区域（图B.1中HABC所围成部分）不是唯一的一个点，而是一个面，所以磁定向法普遍存在一定误差。图B.1二站振幅误差计算示意图B.3计算方法设H（H,H）点为闪电落点，M（M,M）,N（N,N）为二个探测站点，三点组成的HMN对边设为h、m、n，为误差角，MHN设为，利用经纬坐标差计算各边距离。h

14、=dx(lM-lN)2cos2FM+(FM-FN)2(1-f)2)1/2；m=dx(lH-lN)2cos2FH+(FH-FN)2(1-f)2)1/2；n=dx(lM-lH)2cos2FM+(FM-FN)2(1-f)2)1/2。8DB43/T21332021式中：dx=2pRd/360赤道上相距一个经度的距离；Rd地球赤道半径；f地球扁率。应用余弦定理计算角值：；b=arccoscosb=n2+m2-h22mnn2+m2-h22mn。应用正弦定理：在HMC中，HC=nsinq/sin(b+q)；在HAN中，HA=msinq/sin(b+q)。应用余弦定理：AC2=HC2+HA2-2HCHAcos

15、b；AC=sinqsin(b+q)(m2+n2-2mncosb)1/2=hsinqsin(b+q)（B.1）因误差角为1或0.5，将四边形HABC视为近似平行四边形，则：HAN=HCM=𝜋𝜋-,HABC，有：HB2=HC2+HA2+2HCHAcosb；HB=sinqsin(b+q)(m2+n2+2mncosb)1/2=sinqsin(b+q)(2m2+2n2-h2)1/2（B.2）选择AC、HB中的最大值作为二站振幅定位方式可能的雷电落点最大误差，并用（B.3）式表示：gz=12(max(abs(AC),abs(HB)（B.3）9rij=（C.1）DB43/T2

16、1332021附录C（资料性）权重大小确定方法C.1熵值法C.1.1原理熵是系统无序程度的度量，可用于度量已知数据所包含的有效信息量和确定权重。通过对熵的计算确定权重，即根据各项指标值的差异程度，确定各指标的权重。当评价对象的某项指标值相差较大时，熵值较小，说明指标提供的有效量较大，其权重也应较大；反之，若某项指标值相差较小，熵值较大，说明该指标提供的信息量较小，其权重也应较小。当各被评价对象的某项指标值完全相同时，熵值达到最大，意味着该指标无有用信息，可以从评价指标体系中去除。C.1.2步骤C.1.2.1原始数据矩阵归一化设m个评价指标i、n个被评价对象j，构成原始数据矩阵（aij）mn，对

17、其归一化后，对大者为优的评价指标，归一化公式为：aij-minjaijmaxjaij-minjaij式中：rijaijj站（格）点第i个指标的归一化值；j站（格）点第i个指标值。rij=（C.2）对小者为优的评价指标，归一化公式为：maxjaij-aijmaxjaij-minjaij式中：rijaijj站（格）点第i个指标的归一化值；j站（格）点第i个指标值。hi=-kfijlnfij（C.3）C.1.2.2定义熵在有m个评价指标i、n个被评价对象j的评估案例中，第i个评价指标的熵的计算公式为：mj=1式中：hi第i个评价指标的熵；k与样本m有关的常数，一般为1/lnm；fij第j项评价指标下

18、第i个评估方案占该评价指标的比重，计算公式见（C.4）。10标度含义aij=1因素与因素具有相等的重要性aij=3因素比稍显重要aij=5因素比明显重要aij=7因素比强烈重要aij=9因素比极度重要aij=2、4、6、8因素与因素相比，介于结果的中间值倒数aij=1/aijfij=rij/rij（C.4）DB43/T21332021mj=1式中：rij归一化后的值。C.1.2.3定义熵权定义了第i个评价指标的熵之后，按照公式（C.5）可得到第i个评价指标的熵权：m-hi(0wi1,wi=1)（C.5）wi=1-himi=1mi=1式中：wi第i个评价指标的熵权；hi第i个评价指标的熵；m评价

19、指标数量。C.2层次分析法C.2.1原理把一个复杂系统中的每个指标都分解为若干个有序层次，每一层次中的元素具有大致相等的地位，并且每一层与上一层的某个指标和下一层的若干指标有着一定的联系，每一个层次之间按照隶属关系组建成一个有序的递阶层次结构模型。在这个层次结构模型中，根据客观事实的判断，通过两两比较判断的方式确定同一层次中每个指标的相对重要性，以数字的方式建立判断矩阵，然后利用向量的计算方法得出同一层次中每个指标的相对重要性权重系数，最后通过组合计算所有层次的相对权重系数得到每个最底层指标相对于目标的重要性权重系数。C.2.2步骤C.2.2.1构造判断矩阵采用19标度法对各指标进行成对比较，

20、确定各指标间的相对重要性并给出相应的比值，见表C.1。表C.1两两比较赋值表11判断矩阵的阶数R.I.102030.5240.951.1261.2671.36DB43/T21332021上述过程得出的判断矩阵A为：a=21A=(aij)mna11.am1a12a22.am2.a1n.a2n.amn（C.6）式中：aij=1，aij=1/aij。C.2.2.2计算相对权重通过求解判断矩阵A的最大特征值lmax及最大特征值对应的特征向量W，得出同一层次各指标的相对权重系数。C.2.2.3一致性检验用平均随机一致性指标（R.I.）对各指标重要程度比较链上的相容性进行检验，当成对比较得出的判断矩阵的阶

21、数大于或等于3时，则需要进行一致性检验。根据判断矩阵得出一致性指标（C.I.）：C.I.=lmax-nn-1（C.7）根据判断矩阵阶数，按照表C.2找出对应的R.I.。表C.2平均随机一致性指标值C.R.=（C.8）根据C.I.和R.I.的值，计算一致性比例（C.R.）:C.I.R.I.当C.R.小于等于0.1时，则判断矩阵A的一致性是符合要求的，反之，需要对判断矩阵A的两两比较值作调整，直到计算出符合一致性要求的值。C.2.2.4计算合成权重当所有层次的相对权重计算得出后，利用各层次指标的层次单排序结果，进一步计算递阶层次结构模型中最底层指标相对于总目标的组合权重，由下而上逐层进行，进行层次

22、总排序。12DB43/T21332021附录D（资料性）自然断点法自然断点法（Jenksnaturalbreaksmethod）是一种地图分级法。该算法认为数据本身有断点，可利用数据这一特点进行分级。算法原理是一个小聚类，聚类结束条件是组间方差最大、组内方差最小。计算方法见式（D.1）：SSDij=Ak-(Ak)2jk=12jk=1j-i+1（D.1）式中：SSD方差；i、j第i、j个元素；A长度为N的数组；Ki、j中间的数，表示A组中的第k个元素。13DB43/T21332021参考文献1GB/500572010建筑物防雷设计规范2GB/T370472018基于雷电定位系统（LLS）的地闪密度总则3QX/T4052017雷电灾害风险区划技术指南4QX/T792007闪电监测定位系统第1部分：技术条件5MZ/T0272011自然灾害风险管理基本术语6DB43/T18022020雷电预警等级划分技术规范6章国材.自然灾害风险评估与区划原理和方法M.北京：气象出版社，20148刘凤姣,薛明,薛德锋，等.基于闪电数据的雷电风险优化研究模型与应用J.灾害学,2020,35（4）：48-53，8414