THE DISTRIBUTION CHARACTERISTICS OF MULTIPLE STROKES OF CLOUD-TO-GROUND LIGHTNING AND THEIR LIGHTNING CURRENT MAGNITUDE BASED ON LLS
-
摘要: 为进一步研究多回击地闪参数分布特征, 以便为雷电防护工程设计和雷电物理研究提供参考, 根据湖北省雷电定位系统(LLS)2007年1月—2018年12月监测资料, 采用计算机编程处理和数理统计方法, 对多回击地闪次数、多重回击次数和不同类型多回击地闪雷电流幅值等参数进行了统计分析。结果表明: 多回击正地闪、负地闪和总地闪次数占其地闪总数的百分比分别为2.06%、34.76%和32.64%, 多重回击次数分别占其回击总数的0.01%、0.42%和0.40%, 多回击负地闪回击次数占多回击总地闪回击总数的99.69%。首次回击强度大于后续回击强度的多回击正地闪和负地闪分别占多回击地闪总数的82.52%和57.87%;在多回击地闪后续回击中, 正地闪约有9%的后续回击强度大于首次回击强度, 负地闪约有20%的后续回击强度大于首次回击强度。多回击正地闪和负地闪中值电流分别为59.30 kA和35.10 kA, 首次回击分别为90.90 kA和40.00 kA, 后续回击中分别为43.90 kA和33.00 kA。首次回击中, 多回击正地闪和负地闪雷电流幅值大于100 kA的累积概率分别为44.06%和4.64%, 首次回击强度大于后续回击强度的多回击正地闪和负地闪雷电流幅值大于100 kA的累积概率最大分别为52.21%和7.94%;后续回击中, 多回击正地闪和负地闪雷电流幅值小于等于40 kA的累积概率分别为41.80%和69.92%, 首次回击强度大于后续回击强度的多回击负地闪, 雷电流幅值小于等于40 kA的累积概率最大为77.71%。多回击正地闪和负地闪后续回击与首次回击中值电流的比值分别为0.48和0.83。拟合得出的不同类型的多回击正地闪和负地闪雷电流幅值累积概率公式, 拟合效果显著; 拟合公式中a值附近的雷电流幅值累积概率与b值呈显著正相关关系。Abstract: To further investigate the distribution characteristics of multiple strokes of cloud-to-ground(CG)lightning, and to provide reference for the design of lightning protection engineering and lightning physics research, the present study analyzes the parameters such as multiple strokes lightning frequency, multiple return-strokes lightning frequency and lightning current magnitude for different types of multiple strokes lightning using computer programing and mathematical statistics method. These data were collected by the Lightning Location System(LLS) of Hubei Meteorological Service from January 2007 to December 2018.The results show that the percentage for multiple strokes of positive, negative and total CG lightning was 2.06%, 34.76% and 32.64% respectively, and the number of multiple return-strokes was 0.01%, 0.42% and0.40% of the total number of return strokes respectively. The number of negative multiple strokes accounted for 99.69% of the multiple strokes of total CG lightning. The percentage of multiple strokes of positive and negative CG lightning in which the first return stroke intensity was greater than the subsequent was 82.52% and 57.87% of the total number of return strokes respectively. In the subsequent return stroke of multiple strokes, the percentage of the subsequent return stroke intensity of positive and negative CG lightning being greater than that of the first return stroke was about 9% and 20%, respectively. The median current for multiple strokes of positive and negative CG lightning was 59.30 kA and 35.10 kA, respectively. The first return stroke intensity was 90.90 kA and 40.00 kA, respectively, and the subsequent return stroke intensity was 43.90 kA and 33.00 kA, respectively. In the first return stroke, the cumulative probability for multiple strokes of positive and negative CG lightning in which current magnitude was greater than 100 kA was 44.06% and 4.64%, respectively. When the first return stroke intensity was greater than the subsequent return stroke intensity, the maximum cumulative probability of the positive and negative CG lightning in which current magnitude was greater than 100 kA was 52.21%and 7.94%, respectively. In the subsequent return stroke, the cumulative probability of multiple strokes of positive and negative CG lightning with current magnitude being less than or equaling to 40 kA was 41.80% and 69.92%, respectively. Besides, when the first return stroke intensity was greater than the subsequent return stroke intensity, the maximum cumulative probability of the negative CG lightning with current magnitude being less than or equaling to 40 kA was 77.71%. In the multiple strokes of positive and negative CG lightning, the ratio of the subsequent return stroke to the first return stroke was 0.48 and0.83, respectively. The fitting results show that the cumulative probability formula of the current magnitude for multiple strokes of the positive and negative CG lightning with different types was significant. There was a significant positive correlation between the cumulative probability of lightning current magnitude near value a and value b in the fitting formula.
-
表 1 2007—2018年多回击正地闪、负地闪和总地闪次数及其百分比随回击次数变化
回击次数 多回击正地闪次数 多回击负地闪次数 多回击总地闪次数 多回击地闪次数占其地闪总数的百分比/% 多回击地闪次数占其多回击地闪总数的百分比/% 正地闪 负地闪 总地闪 正地闪 负地闪 总地闪 2 2 610 384 619 387 229 2.01 20.55 19.35 97.28 59.12 59.28 3 66 148 196 148 262 0.05 7.92 7.41 2.46 22.78 22.70 4 7 65 006 65 013 0.01 3.47 3.25 0.26 9.99 9.95 5 0 29 949 29 949 0.00 1.60 1.50 0.00 4.60 4.58 6 0 13 172 13 172 0.00 0.70 0.66 0.00 2.02 2.02 7 0 5 737 5 737 0.00 0.31 0.29 0.00 0.88 0.88 8 0 2 402 2 402 0.00 0.13 0.12 0.00 0.37 0.37 9 0 922 922 0.00 0.05 0.05 0.00 0.14 0.14 10 0 350 350 0.00 0.02 0.02 0.00 0.05 0.05 11 0 140 140 0.00 0.01 0.01 0.00 0.02 0.02 12 0 42 42 0.00 0.00 0.00 0.00 0.01 0.01 13 0 11 11 0.00 0.00 0.00 0.00 0.00 0.00 14 0 4 4 0.00 0.00 0.00 0.00 0.00 0.00 15 0 2 2 0.00 0.00 0.00 0.00 0.00 0.00 16 0 1 1 0.00 0.00 0.00 0.00 0.00 0.00 合计 2 683 650 553 653 236 2.06 34.76 32.64 100.00 100.00 100.00 表 2 2007—2018年负地闪和总地闪多重回击次数及其百分比随回击次数变化
回击次数 后续回击次数/次 多重回击次数/次 多重回击次数占其地闪回击总次数的百分比/% 多重回击次数占其后续回击次数的百分比/% 负地闪 负地闪 总地闪 负地闪 总地闪 负地闪 总地闪 2 384 619 4 415 4 426 0.15 0.14 1.15 1.14 3 296 392 3 280 3 280 0.11 0.10 1.11 1.11 4 195 018 2 135 2 135 0.07 0.07 1.09 1.09 5 119 796 1 275 1 275 0.04 0.04 1.06 1.06 6 65 860 729 729 0.02 0.02 1.11 1.11 7 34 422 372 372 0.01 0.01 1.08 1.08 8 16 814 193 193 0.01 0.01 1.15 1.15 9 7 376 74 74 0.00 0.00 1.00 1.00 10 3 150 32 32 0.00 0.00 1.02 1.02 11 1 400 13 13 0.00 0.00 0.93 0.93 合计 1 124 847 12 518 12 529 0.42 0.40 1.11(平均) 1.11(平均) 表 3 2007—2018年不同类型的多回击地闪次数、回击次数及其百分比分布
项目 T型多回击地闪 F型多回击地闪 S型多回击地闪 后续回击强度大于首次 地闪总数/次 回击总数/次 次数/次 占其地闪总数百分比/% 次数/次 占其地闪总数百分比/% 次数/次 占其回击总数百分比/% 正地闪 2 683 5 446 2 214 82.52 469 17.48 481 8.83 负地闪 650 553 1 776 089 376 448 57.87 274 105 42.13 353 245 19.89 总地闪 653 236 1 781 535 378 662 57.97 274 574 42.03 353 726 19.86 表 4 2007—2018年不同多回击地闪、回击类型雷电流幅值累积概率(%)
地闪极性 雷电流幅值/kA 全回击 首次回击 后续回击 F型 S型 T型 F型 S型 T型 F型 S型 T型 正地闪 > 100 27.65 12.38 24.95 52.21 5.54 44.06 3.70 18.90 6.41 > 150 10.64 3.54 9.38 21.05 1.28 17.59 0.48 5.69 1.41 > 200 3.26 0.94 2.85 6.55 0.00 5.40 0.04 1.83 0.36 ≤40 25.86 36.63 27.76 5.56 49.89 13.31 45.66 23.98 41.80 负地闪 > 100 3.07 0.57 1.97 7.94 0.11 4.64 0.11 0.82 0.43 > 150 0.63 0.05 0.37 1.66 0.01 0.96 0.01 0.07 0.03 > 200 0.14 0.01 0.08 0.37 < 0.01 0.22 < 0.01 0.01 0.01 ≤40 60.07 66.01 62.68 31.01 76.42 50.15 77.71 60.37 69.92 表 5 2007—2018年不同多回击地闪、回击类型雷电流幅值累积概率拟合公式中的a、b参数值
地闪极性 拟合公式
a、b参数值全回击 首次回击 后续回击 F型 S型 T型 F型 S型 T型 F型 S型 T型 正地闪 a 62.01 50.12 59.24 103.83 40.07 90.56 41.94 60.33 43.97 b 2.52 3.13 2.58 3.31 3.39 2.88 4.15 3.22 3.61 负地闪 a 35.92 34.27 35.08 49.11 31.52 39.94 30.90 36.19 33.04 b 3.86 5.07 4.14 4.02 5.90 3.79 5.85 4.95 5.15 -
[1] 黎勋, 郄秀书, 刘昆, 等. 基于高时间分辨率快电场变化资料的北京地区地闪回击统计特征[J]. 气候与环境研究, 2017, 22(2): 231-241. [2] 张阳, 张义军, 孟青, 等. 北京地区正地闪时间分布及波形特征[J]. 应用气象学报, 2010, 21(4): 442-449. [3] 陈绿文, 吕伟涛, 张义军, 等. 不同高度建筑物上的下行地闪回击特征[J]. 应用气象学报, 2015, 26(3): 311-318. [4] 付茂金, 阮小飞, 王州龙, 等. 高速铁路通信信号综合防雷技术[M]. 北京: 科学出版社, 2014: 1-8. [5] 孙萍. 220kV新杭线雷电流幅值实测结果的统计分析[J]. 中国电力, 2000, 33(3): 72-75. [6] 张义军, 吕伟涛, 郑栋, 等. 负地闪先导-回击过程的光学观测和分析[J]. 高电压技术, 2008, 34(10): 2 022-2 029. [7] COORAY V, PEREZ H. Some features of lightning flashes observed in Sweden[J]. J Geophy Res, 1994, 99(D5): 10 683-10 688. [8] WILLIAM C V, PHILIP K E. Statistics and characteristics of cloud-to-ground lightning with multiple ground contacts[J]. J Geophy Res, 2002, 107(20), 4441, doi:10.1029/2001JD001360. [9] 王学良, 张科杰, 张义军, 等. 雷电定位系统与人工观测雷暴日数统计比较[J]. 应用气象学报, 2014, 25(6): 741-750. [10] 王学良, 刘学春, 黄小彦, 等. 湖北地区云地闪电时空分布特征分析[J]. 气象, 2010, 36(10): 91-96. [11] 王学良, 张科杰, 余田野, 等. 湖北地区云地闪电频次及雷电流幅值时间分布特征[J]. 电瓷避雷器, 2017(3): 1-9. [12] 王学良, 张科杰, 余田野, 等. 湖北省山区与平原雷电分布及其参数特征[J]. 气象科技, 2019, 47(2): 337-348. [13] 王学良, 余田野, 贺姗. 湖泊与陆地雷电参数分布特征及其对比分析[J]. 气象与环境学报, 2019, 35(6): 132-138. [14] 王学良, 余田野, 贺姗, 等. 区域海拔高度对云地闪电参数分布的影响[J]. 高电压技术, 2020, 46(4): 1 206-1 215. [15] 费蕾蕾, 毕新慧, 刘永林, 等. 香港地区地闪时空分布特征及其影响因素[J]. 热带气象学报, 2017, 33(5): 617-626. [16] 郭润霞, 王迎春, 张文龙, 等. 基于VLF/LF三维闪电监测定位系统的北京闪电特征分析[J]. 热带气象学报, 2018, 34(3): 393-400. [17] HUSSEIN A M, JANISCHEWSKY J W, CHANG J S, et al. Simultaneous measurement of lightning parameters for strokes to the Toronto Canadian National Tower[J]. J Geophy Res, 1995, 100(D5): 8 853-8 861. [18] 朱时阳, 周文俊, 邓雨荣, 等. 可记录多重雷击的雷电流波形在线监测装置研制[J]. 高电压技术, 2013, 39(11): 2 699-2 705. [19] 孔祥贞, 郄秀书, 张广庶, 等. 多接地点闪电的梯级先导与回击过程的研究[J]. 中国电机工程学报, 2005, 25(22): 142-147. [20] 曾金全, 朱彪, 曾颖婷, 等. 福建省多回击闪电特征参数的统计分析[J]. 暴雨灾害, 2017, 36(6): 573-578. [21] 许维伟, 李在光, 祝宝友, 等. 基于不间断闪电波形采集的一次皖北雷暴负地闪特征观测[J]. 高原气象, 2015, 34(3): 850-862. [22] 陈家宏, 赵淳, 谷山强, 等. 我国电网雷电监测与防护技术现状及发展趋势[J]. 高电压技术, 2016, 42(11): 3 361-3 375. [23] 田芳, 肖稳安, 冯民学, 等. 闪电定位仪观测结果的修订分析[J]. 华东电力, 2008, 36(6): 38-42. [24] 曾庆峰, 张其林, 赖鑫, 等. 深圳市闪电定位资料误差分析及其优化[J]. 气象科技, 2015, 43(3): 530-536. [25] 中国科学院空间科学与应用研究中心. ADTD雷电探测仪用户手册[R]. 北京: 中国科学院, 2004: 2-3. [26] 秦建峰, 刘梦雨, 吴昊. ADTD雷电探测系统典型故障分析[J]. 气象科技, 2012, 40(2): 180-184. [27] 谷山强, 王剑, 张义军, 等. 基于雷电定位系统(LLS)的地闪密度总则(GB/T 37047-2018, IEC 62858: 2015, MOD)[S]. 北京: 中国标准出版社, 2015: 1-4. [28] WIKIPEDIA. Haversine formula[EB/OL]. [2019-01-06]. https://en.wikipedia.org/wiki/Haversine_formula. [29] 王胜开, 徐志洁, 张健钦, 等. 逆向热力图的绘制方法[J]. 地球信息科学学报, 2018, 20(4): 515-522. [30] 郄秀书, 郭昌明, 刘欣生. 北京与兰州地区的地闪特征[J]. 高原气象, 1990, 9(4): 388-394. [31] RAKOV V A. THOTTAPPILLIL R. Review of lightning properties from electric field and TV observations[J]. J Geophy Res, 1994, 99(D5): 10 745-10 750. [32] 李京校, 郭凤霞, 沈永海, 等. 利用SAFIR资料对北京及其周边地区地闪参数的特征分析[J]. 高原气象, 2013, 32(5): 1 450-1 459. [33] 郄秀书, 余晔, 王怀斌, 等. 中国内陆高原地闪特征的统计分析[J]. 高原气象, 2001, 20(4): 395-401. [34] 林维勇, 黄有根, 焦兴学, 等. 建筑物防雷设计规范(GB 50057-2010)[S]. 北京: 中国计划出版社, 2011: 43-58. [35] 陈家宏, 童雪芳, 谷山强, 等. 雷电定位系统测量的雷电流幅值分布特征[J]. 高电压技术, 2008, 34(9): 1 893-1 897. [36] 李瑞芳, 吴广宁, 曹晓斌, 等. 雷电流幅值概率计算公式[J]. 电工技术学报, 2011, 26(4): 161-167. [37] 刘刚, 唐军, 孙雷雷, 等. 不同地形地貌的雷电流幅值概率分布对输电线路雷击跳闸的影响[J]. 高电压技术, 2013, 39(1): 17-23. [38] 豆朋, 黄松波. 佛山地区雷电流幅值概率分布研究[J]. 电瓷避雷器, 2011, 3: 43-45. [39] GRANT I S, ANDERSON J G, HILEMAN A R, et al. A simplified method for estimating lightning performance of transmission lines[J]. IEEE Transaction on Power Apparatus and Systems, 1985, 104(4): 919-932. [40] TAKAMI J, OKABE S. Observational results of lightning current on transmission towers[J]. IEEE Transactions on Power Delivery, 2007, 22(1): 547-556. [41] IEEE STANDARDS BOARD. IEEE Std 1243-1997 IEEE Guide for Improving the Lightning Performance of Transmission Lines[S]. New York IEEE Inc, 1997: 8. [42] 李家启, 王劲松, 申双和, 等. 基于ADTD系统监测的雷电流幅值累积概率特征分析[J]. 气象, 2011, 37(2): 226-231. [43] 赵淳, 雷梦飞, 王剑, 等. 雷电流幅值累积概率曲线拟合方法[J]. 高电压技术, 2018, 44(5): 1 598-1 604.