CN113204903A - 一种预测雷电的方法 - Google Patents
一种预测雷电的方法 Download PDFInfo
- Publication number
- CN113204903A CN113204903A CN202110474699.5A CN202110474699A CN113204903A CN 113204903 A CN113204903 A CN 113204903A CN 202110474699 A CN202110474699 A CN 202110474699A CN 113204903 A CN113204903 A CN 113204903A
- Authority
- CN
- China
- Prior art keywords
- data
- lightning
- intermediate data
- prediction
- thunder
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000008859 change Effects 0.000 claims abstract description 13
- 238000012545 processing Methods 0.000 claims abstract description 13
- 238000007781 pre-processing Methods 0.000 claims abstract description 12
- 238000013528 artificial neural network Methods 0.000 claims abstract description 11
- 238000002310 reflectometry Methods 0.000 claims description 18
- 230000009466 transformation Effects 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 9
- 238000012216 screening Methods 0.000 claims description 5
- 230000000306 recurrent effect Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims description 2
- 238000002592 echocardiography Methods 0.000 abstract 2
- 125000004122 cyclic group Chemical group 0.000 abstract 1
- 230000005684 electric field Effects 0.000 description 9
- 230000006870 function Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 230000006399 behavior Effects 0.000 description 4
- 238000013277 forecasting method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000033001 locomotion Effects 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000013135 deep learning Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000011158 quantitative evaluation Methods 0.000 description 2
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 1
- 208000025274 Lightning injury Diseases 0.000 description 1
- 101150114976 US21 gene Proteins 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
- G01S13/958—Theoretical aspects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/418—Theoretical aspects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/044—Recurrent networks, e.g. Hopfield networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Radar, Positioning & Navigation (AREA)
- Computing Systems (AREA)
- Remote Sensing (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Human Resources & Organizations (AREA)
- Biomedical Technology (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- Economics (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Geometry (AREA)
- Computer Networks & Wireless Communication (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Computer Hardware Design (AREA)
- Medical Informatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Entrepreneurship & Innovation (AREA)
- Electromagnetism (AREA)
- Educational Administration (AREA)
- Primary Health Care (AREA)
Abstract
本发明属于雷电预警技术领域,公开了一种预测雷电的方法。包括数据预处理,雷电预测和输出预测结果三个阶段,其中数据预处理阶段用图像表征待处理数据的属性并模糊处理;雷电预测阶段是运用卷积运算抽取所述经模糊处理图像的特征,通过循环神经网络编码后得到反映雷电属性变化的函数,再进行反卷积运算将该函数表征的雷电属性未来变化状况解码为动态灰度图像;该方法消除了非对流性回波或非闪电位置回波导致的冗余和干扰,构建了与雷电发生直接相关的雷暴云团的局部空间特征和时序变化规律,从而提高了雷电预测的准确率。
Description
技术领域
本发明属于雷电预警技术领域,具体公开了一种预测雷电的方法。
背景技术
雷电是强对流天气过程的产物。强对流天气过程演变迅速,云团生成、消散、合并和***规律复杂,导致雷电的预报预警难度大。目前雷电预报方法主要有如下几种:
利用大气电场仪或其组成的网络预报雷电。譬如,中国专利CN110297284A“一种基于三维大气电场仪的雷暴云移动路径跟踪方法”公开了利用大气电场仪探测雷暴云方位和跟踪其移动路径的方法。中国专利CN106353604B“一种基于大气电场背景噪声自动追踪的雷电预警方法”提出了一种剔除非雷暴时背景电场噪声,进而提高雷电预警准确度的方法。
基于雷电定位数据外推预测雷暴移动路径。譬如,论文《基于闪电数据的雷暴识别、追踪和外推方法》提出了先利用密度聚类算法识别雷暴,再利用卡尔曼滤波算法追踪和外推雷暴路径的方法。
基于气象雷达预报雷电。譬如,中国专利CN110927725“一种基于气象雷达的雷电预警和监测方法”提供了一种基于气象雷达回波功率和反射率因子的雷电预警方法。中国专利CN106019287A“一种基于小型雷达的输电线路雷电监测预警方法”公开了结合小型雷达反射率因子和雷电定位***数据预报输电杆塔雷击风险的方法。
基于卫星监测数据预报雷电。如中国专利CN10944990A“一种基于气象卫星数据的雷电预警方法”建立了云顶亮温变化率和雷电之间的关系,进而对雷电进行风险等级预警。
基于数值天气预报模式开展雷电预报。如中国专利CN109447315A“一种基于多时空尺度的电力气象数值天气预报方法和装置”提供了利用中尺度天气预报模式开展包括雷电在内的多种气象灾害多尺度预报方法。
除了以上回顾的单独利用大气电场仪、雷电定位***、雷达、卫星和数值天气预报模式开展雷电预报,还有综合多种手段开展雷电预报的方法,如中国专利CN109738970A“基于雷电数据挖掘实现雷电预警的方法、装置和存储介质”公开了综合利用大气电场强度、雷达回波、雷电定位数据和卫星云图,预测雷暴团移动速度和方向的方法。
总体而言,目前各种雷电预报方法提出的判据比较单一,通常利用某物理参量(如大气电场强度或其变化率、云顶亮温变化率、雷达回波dBZ)的阈值或相关关系公式对雷暴的发生或移动路径进行预测,难以准确刻画雷暴云团运动和生、消、聚、变的复杂规律。
发明内容
本发明的目的在于针对现有技术的不足,提出了一种预测雷电的方法,包括数据预处理,雷电预测和输出预测结果三个步骤,所述数据预处理是用图像表征待处理数据的属性并模糊处理;所述雷电预测是运用卷积运算抽取所述经模糊处理图像的特征,通过循环神经网络编码后得到反映雷电属性变化的函数,再进行反卷积运算将该函数表征的雷电属性未来变化解码为动态灰度图像。
进一步的,数据预处理包括以下步骤:
S1.对获取的雷达基数据进行解析并提取出雷达反射率因子数据,用最邻近法插值到以雷达站为中心的固定大小的等经纬度网格;
S2.取同一位置不同高度上的最大反射率因子值生成组合反射率,把所述组合反射率转换为灰度值,并保存为灰度图。
S3.在雷达基数据处理的范围内,以雷达扫描一次时间为周期,提取对应时间内雷电定位散点数据。
S4.将所述雷电定位散点数据网格化,采用与所述雷达基数据相同的网格数和网格分辨率,在每一个网格内对雷电进行计数,然后转换为灰度值,并保存为和所述雷达基数据图片格式相同的雷电定位散点数据的灰度图。
S5.对所述雷电定位散点数据的灰度图进行模糊处理。
优选的,所述编码和解码之间***数据筛选步骤;所述编码构成神经网络计算模块,所述神经网络模块包括:
将所述S2和S5步骤中灰度图经卷积网络变换和下采样后形成中间数据E0组成的输入层;
以所述中间数据E0作为输入,经过卷积网络变换后形成中间数据E1组成的第一层;
以所述中间数据E1作为输入,对其进行下采样,获得中间数据E12组成的下采样层;
以所述中间数据E12作为输入,对其进行卷积网络变换,获得中间数据E2组成的第二层;
接收上一个时次的E1和本时次的E0,输出本时次的E1以及接收上一个时次的E2和本时次的E12,输出本时次的E2的时间步进演算层;
数据筛选步骤包括:对所述编码构成神经网络计算模块的第一层和第二层每个时间步输出的中间数据进行编码,所述第一层得到中间数据A1,所述第二层得到中间数据A2;
所述解码包括以下步骤:
P1.以所述中间数据A2作为输入,经过卷积网络变换后得到中间数据D2;
P2.对所述中间数据D2进行上采样,获得中间数据D21;
P3.以所述中间数据A1作为输入,进行卷积网络变换后获得中间数据D1;
P4.对所述中间数据D21和D1在第三维进行拼接,将拼接后的数据进行卷积网络变换并进行上采样获得中间数据D10。
P5.接收上一个时次的所述中间数据D2和本时次的A2,输出本时次的A2;接收上一个时次的所述中间数据D1和本时次的A1,输出本时次的A1。
P6.对所述中间数据D10进行卷积网络变换,输出多个时次的预测数据,并把该预测数据保存为与输入相同格式的灰度图。
进一步的,输出预测结果包括如下步骤:
F1.将雷电预测结果的灰度图转换为彩色图;
F2.叠加地图作为底图;
F3.添加坐标轴和必要的说明文字。
相对于现有技术,本发明具有以下有益效果:
本发明提供的一种预测雷电的方法通过融合雷达反射率因子数据和雷电定位数据,在机器学习架构中进行自拟合和学习,刻画和雷电发生直接相关的雷暴云团的局部空间特征和时序变化规律,进而能够实现更加可靠而精准的雷电预报。
附图说明
图1为本发明实施例中步骤模块关系图
图2为本发明实施例中雷电定位数据灰度图及其高斯模糊图
图3为本发明实施例中雷电预测模块深度学习框架图
图4为本发明实施例中雷达反射率因子和雷电定位数据融合图
图5为本发明实施例中预测雷电和实测雷电分布图,图中第一行为预测结果,第二行为观测结果
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行描述,所描述的实施例是本发明一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域技术人员在没有做出创新劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
以下结合本说明书附图1至5和具体实施例对本发明作进一步的详细说明:
本发明公开的一种预测雷电的方法,该方法主要包含3个按顺序执行的步骤模块,其中模块1为模块2提供数据样本,模块2输出预测结果,模块3对预测结果进行加工处理,方便用户通过图片理解预测结果。
本发明采用的数据为雷达组合反射率因子图像和雷电定位散点图像,可将临近预报问题抽象为对应的图像时间序列的预测问题。
下面分别对本发明的3个模块进行详细说明。
模块1:数据预处理模块,包括对雷达和雷电数据的预处理。
雷达数据预处理是将雷达基数据解析后绘制为图片。
以中国气象局业务运行的S波段雷达为例,该雷达一次体扫时间6分钟,基数据有11个高度层,最大扫描半径460km。
在对雷达数据处理时,分两个步骤:
第一步:以雷达站为中心,划分出460km×460km的等经纬度网格,网格分辨率约为0.01°×0.01°,对每一个高度层的反射率因子数据用最邻近法插值。
第二步:取同一位置不同高度上的最大反射率因子数据生成组合反射率,按照如下公式将组合反射率转换为灰度值,
GS=dBZ/70×255 (1)
然后绘制灰度图,保存为png格式图片。
雷电定位数据预处理是将雷电定位散点数据绘制为图片。雷电定位数据主要信息包含雷电发生时间、经纬度位置等。处理步骤包括:
第一步:在雷达数据处理的范围内(以雷达为中心,460km×460km范围),以雷达一次体扫时间(6分钟)为周期,提取对应时间内雷电定位散点数据。
第二步:计算每一个1km×1km网格中的雷电数量,按照如下公式雷电数转换为灰度值,
GS=Flash/2×255 (2)
然后绘制灰度图,保存为png格式图片。
第三步:对雷电定位数据的灰度图进行高斯模糊处理。高斯模糊的窗口大小取15km×15km,标准差3。通过图2中模糊前后的图片对比看出高斯模糊效果。高斯模糊的作用是“模糊”掉雷电定位散点数据中比较分散的结果,保留团状、片状雷电定位结果。
模块2:雷电预测模块。该模块采用深度学习方法,包括数据编码、筛选和解码3个步骤,分别由对应的编码网络(E)、注意力网络(A)和解码网络(D)3个模块来执行。
编码网络包含:
输入层:接收10个连续时次的雷达组合反射率因子灰度图和雷电定位数据灰度图,两种数据叠加成460*460*2的数据后(叠加操作如图4所示),经过卷积网络变换和下采样后变成115*115*64的中间数据,该中间数据标记为E0。
第一层:以E0作为输入,经过卷积网络变换后变为115*115*64的中间数据:记为E1。
下采样层DS12:以E1作为输入,对E1进行下采样,获得24*24*128的中间数据,记为E12。
第二层:以E12作为输入,对E12进行卷积网络变换,获得中间数据E2。
时间信息传播模块:该模块采用循环神经网络中的GRU模型。其作用是接收上一个时次的E1和本时次的E0,输出本时次的E1;接收上一个时次的E2和本时次的E12,输出本时次的E2。
注意力网络包含:对编码网络的第一层和第二层每个时间步输出的中间数据进行编码,编码得到中间数据记为A1和A2。
解码网络包含:
第二层:以A2作为输入,经过卷积网络变换后变为24*24*128的中间数据:记为D2。
上采样层US21:以D2作为输入,对D2进行上采样,获得115*115*128的中间数据,记为D21。
第一层:以A1作为输入,进行卷积网络变换后获得中间数据D1,形状为115*115*64。
融合层:对D21和D1在第三维进行拼接,对拼接后的数据进行卷积网络变换并进行上采样获得460*460*32的中间数据D10。
时间信息传播模块:该模块采用循环神经网络中的GRU模型。其作用是接收上一个时次的D2和本时次的A2,输出本时次的A2;接收上一个时次的D1和本时次的A1,输出本时次的A1。
输出层:对D10进行卷积网络变换,输出连续10个时次的460*460*1的预测数据,并把预测数据保存为与输入相同格式的灰度图。
雷电预测模型所采用的优化算法为ADAM优化算法,学习率设为0.002,代价函数是绝对值误差和均方根误差之和:L1+L2,最大训练轮数为200。
模块3:预测结果输出模块。在本模块中,主要执行如下步骤:
1.将雷电预测结果的灰度图转换为彩色图;
2.叠加地图作为底图;
3.添加坐标轴和必要的说明文字。
获得预测结果后,采用如下方法对结果进行定量评估。对于一次预测,存在如下4种情况:
(1)预测有,实际有;将这种情况发生次数记为a。
(2)预测有,实际无;将这种情况发生次数记为b。
(3)预测无,实际有;将这种情况发生次数记为c。
(4)预测无,实际无;将这种情况发生次数记为d。
通过如下5个指标定量评价预测结果,分别为:
正确率(PC):PC=(a+d)/n,其中n=a+b+c+d
命中率(H):H=a/(a+c)
虚警率(FAR):FAR=b/(a+b)
关键成功指数(CSI):CSI=a/(a+b+c)
偏差(B):B=(a+b)/(a+c)
获取模块2得到的雷电预测结果灰度图图像数据,重新对应模块1的经纬度范围(以雷达为中心,460km×460km范围),根据给定色标绘制彩图并叠加地图。图5为2020年03月21日武汉市雷暴过程的预测结果,第一行为预测结果(北京时2020年03月21日22时18分到03月22日00时06分,每12分钟预测一次),第二行为预测结果对应的观测结果。模型正确预测了此次雷暴过程2小时内的雷电主要发生区域,雷电从无到有的趋势变化也得到很好的预测。
针对模块2得到的雷电预测结果图像数据进行定量评估。预测结果图像数据为460km×460km的网格数据,以5km×5km为统计单位对网格数据进行划分,若该5km×5km网格内存在不为0的值,则将新组成的网格记为1,否则记为0,最终得到92×92个网格,同理得到92×92的观测数据网格。将二者网格一一对应针对上文指出的4种情况进行匹配判断,得到评估结果。表1为2020年03月21日22时40分和2020年03月22日00时06分(即起报后的36分钟和2小时)的评估结果,两次雷电预测正确率均达到90%,随着预测时间增加虚警率有所上升,但都低于27%,CSI也有所下降,但高于0.4,可以认为虽然随着预测时间增加,模型预测能力有所下降,总体而言对于此次雷暴过程模型达到了较好的预测效果。
表1雷电预测定量评估结果
通过以上实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可通过软件编程加上通用计算机硬件设备方式来实现,也可以全部通过硬件实现。基于该理解,上述技术方案对现有技术贡献的内容可以通过软件产品的形式实现,该软件产品可以存储在计算设备可读介质中,例如硬盘、光盘等,包含若干指令用来让一台计算设备执行实施例或其中部分方法。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种预测雷电的方法,包括数据预处理,雷电预测和输出预测结果三个步骤,其特征为:所述数据预处理是用图像表征待处理数据的属性并模糊处理;所述雷电预测是运用卷积运算抽取所述经模糊处理图像的特征,通过循环神经网络编码后得到反映雷电属性变化的函数,再进行反卷积运算将该函数表征的雷电属性未来变化解码为动态灰度图像。
2.根据权利要求1所述的一种预测雷电的方法,其特征在于:所述数据预处理包括以下步骤:
S1.对获取的雷达基数据进行解析并提取出雷达反射率因子数据,用最邻近法插值到以雷达站为中心的固定大小的等经纬度网格;
S2.取同一位置不同高度上的最大反射率因子值生成组合反射率,把所述组合反射率转换为灰度值,并保存为灰度图。
S3.在雷达基数据处理的范围内,以雷达扫描一次时间为周期,提取对应时间内雷电定位散点数据。
S4.将所述雷电定位散点数据网格化,采用与所述雷达基数据相同的网格数和网格分辨率,在每一个网格内对雷电进行计数,然后转换为灰度值,并保存为和所述雷达基数据图片格式相同的雷电定位散点数据的灰度图。
S5.对所述雷电定位散点数据的灰度图进行模糊处理。
3.根据权利要求2的一种预测雷电的方法,其特征在于:所述编码和解码之间***数据筛选步骤;所述编码构成神经网络计算模块,所述神经网络模块包括:
将所述S2和S5步骤中灰度图经卷积网络变换和下采样后形成中间数据E0组成的输入层;
以所述中间数据E0作为输入,经过卷积网络变换后形成中间数据E1组成的第一层;
以所述中间数据E1作为输入,对其进行下采样,获得中间数据E12组成的下采样层;
以所述中间数据E12作为输入,对其进行卷积网络变换,获得中间数据E2组成的第二层;
接收上一个时次的E1和本时次的E0,输出本时次的E1以及接收上一个时次的E2和本时次的E12,输出本时次的E2的时间步进演算层;
数据筛选步骤包括:对所述编码构成神经网络计算模块的第一层和第二层每个时间步输出的中间数据进行编码,所述第一层得到中间数据A1,所述第二层得到中间数据A2;
所述解码包括以下步骤:
P1.以所述中间数据A2作为输入,经过卷积网络变换后得到中间数据D2;
P2.对所述中间数据D2进行上采样,获得中间数据D21;
P3.以所述中间数据A1作为输入,进行卷积网络变换后获得中间数据D1;
P4.对所述中间数据D21和D1在第三维进行拼接,将拼接后的数据进行卷积网络变换并进行上采样获得中间数据D10。
P5.接收上一个时次的所述中间数据D2和本时次的A2,输出本时次的A2;接收上一个时次的所述中间数据D1和本时次的A1,输出本时次的A1。
P6.对所述中间数据D10进行卷积网络变换,输出多个时次的预测数据,并把该预测数据保存为与输入相同格式的灰度图。
4.如权利要求1所述的一种预测雷电的方法,其特征在于所述输出预测结果包括如下步骤:
F1.将雷电预测结果的灰度图转换为彩色图;
F2.叠加地图作为底图;
F3.添加坐标轴和必要的说明文字。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110474699.5A CN113204903B (zh) | 2021-04-29 | 2021-04-29 | 一种预测雷电的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110474699.5A CN113204903B (zh) | 2021-04-29 | 2021-04-29 | 一种预测雷电的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113204903A true CN113204903A (zh) | 2021-08-03 |
CN113204903B CN113204903B (zh) | 2022-04-29 |
Family
ID=77029446
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110474699.5A Active CN113204903B (zh) | 2021-04-29 | 2021-04-29 | 一种预测雷电的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113204903B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114049545A (zh) * | 2021-11-05 | 2022-02-15 | 中山大学 | 一种基于点云体素的台风定强方法、***、设备及介质 |
CN115169431A (zh) * | 2022-05-16 | 2022-10-11 | 湖南师范大学 | 一种雷暴识别方法及*** |
CN116629142A (zh) * | 2023-07-24 | 2023-08-22 | 合肥工业大学 | 基于transformer机制的雷电定位轨迹预测方法、***和存储介质 |
CN117115177A (zh) * | 2023-08-22 | 2023-11-24 | 南通大学 | 基于动态通道图卷积与多尺度注意力的闪电通道分割方法 |
CN117809192A (zh) * | 2024-03-01 | 2024-04-02 | 南京信息工程大学 | 一种基于denclue聚类算法的雷暴识别方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6405134B1 (en) * | 2000-08-30 | 2002-06-11 | Weatherdata, Inc. | Method and apparatus for predicting lightning threats based on radar and temperature data |
JP2009192312A (ja) * | 2008-02-13 | 2009-08-27 | Nippon Telegr & Teleph Corp <Ntt> | 発雷予測装置およびその方法、プログラム |
CN109738970A (zh) * | 2018-12-07 | 2019-05-10 | 国网江苏省电力有限公司电力科学研究院 | 基于雷电数据挖掘实现雷电预警的方法、装置和存储介质 |
CN110579823A (zh) * | 2019-09-02 | 2019-12-17 | 中国电力科学研究院有限公司 | 一种强降雨短临预报方法及*** |
CN110796299A (zh) * | 2019-10-23 | 2020-02-14 | 国网电力科学研究院武汉南瑞有限责任公司 | 一种雷电预测方法 |
CN110989040A (zh) * | 2019-11-29 | 2020-04-10 | 张烨方 | 一种基于切片处理的人工智能雷电临近预警方法及*** |
US20200132884A1 (en) * | 2018-10-30 | 2020-04-30 | Climacell Inc. | Forecasting method with machine learning |
CN111666903A (zh) * | 2020-06-10 | 2020-09-15 | 中国电子科技集团公司第二十八研究所 | 一种卫星云图中雷暴云团的识别方法 |
CN111881723A (zh) * | 2020-06-09 | 2020-11-03 | 南瑞集团有限公司 | 雷电地闪回击波形的自动识别方法、***及电子设备 |
CN111897030A (zh) * | 2020-07-17 | 2020-11-06 | 国网电力科学研究院有限公司 | 一种雷暴预警***及方法 |
CN112363140A (zh) * | 2020-11-05 | 2021-02-12 | 南京叁云科技有限公司 | 一种基于循环神经网络的热力约束外推客观订正方法 |
CN112396116A (zh) * | 2020-11-24 | 2021-02-23 | 武汉三江中电科技有限责任公司 | 一种雷电检测方法、装置、计算机设备及可读介质 |
-
2021
- 2021-04-29 CN CN202110474699.5A patent/CN113204903B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6405134B1 (en) * | 2000-08-30 | 2002-06-11 | Weatherdata, Inc. | Method and apparatus for predicting lightning threats based on radar and temperature data |
JP2009192312A (ja) * | 2008-02-13 | 2009-08-27 | Nippon Telegr & Teleph Corp <Ntt> | 発雷予測装置およびその方法、プログラム |
US20200132884A1 (en) * | 2018-10-30 | 2020-04-30 | Climacell Inc. | Forecasting method with machine learning |
CN109738970A (zh) * | 2018-12-07 | 2019-05-10 | 国网江苏省电力有限公司电力科学研究院 | 基于雷电数据挖掘实现雷电预警的方法、装置和存储介质 |
CN110579823A (zh) * | 2019-09-02 | 2019-12-17 | 中国电力科学研究院有限公司 | 一种强降雨短临预报方法及*** |
CN110796299A (zh) * | 2019-10-23 | 2020-02-14 | 国网电力科学研究院武汉南瑞有限责任公司 | 一种雷电预测方法 |
CN110989040A (zh) * | 2019-11-29 | 2020-04-10 | 张烨方 | 一种基于切片处理的人工智能雷电临近预警方法及*** |
CN111881723A (zh) * | 2020-06-09 | 2020-11-03 | 南瑞集团有限公司 | 雷电地闪回击波形的自动识别方法、***及电子设备 |
CN111666903A (zh) * | 2020-06-10 | 2020-09-15 | 中国电子科技集团公司第二十八研究所 | 一种卫星云图中雷暴云团的识别方法 |
CN111897030A (zh) * | 2020-07-17 | 2020-11-06 | 国网电力科学研究院有限公司 | 一种雷暴预警***及方法 |
CN112363140A (zh) * | 2020-11-05 | 2021-02-12 | 南京叁云科技有限公司 | 一种基于循环神经网络的热力约束外推客观订正方法 |
CN112396116A (zh) * | 2020-11-24 | 2021-02-23 | 武汉三江中电科技有限责任公司 | 一种雷电检测方法、装置、计算机设备及可读介质 |
Non-Patent Citations (3)
Title |
---|
DIEU TIEN BUI等: "A novel deep learning neural network approach for predicting flash flood susceptibility A case study at a high frequency tropical storm area", 《SCIENCE OF THE TOTAL ENVIRONMENT》 * |
杨仲江等: "序列结构的RNN模型在闪电预警中的应用", 《灾害学》 * |
琚泽立等: "基于潜势预报和雷达回波特征的雷电预报方法", 《电瓷避雷器》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114049545A (zh) * | 2021-11-05 | 2022-02-15 | 中山大学 | 一种基于点云体素的台风定强方法、***、设备及介质 |
CN115169431A (zh) * | 2022-05-16 | 2022-10-11 | 湖南师范大学 | 一种雷暴识别方法及*** |
CN116629142A (zh) * | 2023-07-24 | 2023-08-22 | 合肥工业大学 | 基于transformer机制的雷电定位轨迹预测方法、***和存储介质 |
CN116629142B (zh) * | 2023-07-24 | 2023-09-29 | 合肥工业大学 | 基于transformer机制的雷电定位轨迹预测方法和*** |
CN117115177A (zh) * | 2023-08-22 | 2023-11-24 | 南通大学 | 基于动态通道图卷积与多尺度注意力的闪电通道分割方法 |
CN117809192A (zh) * | 2024-03-01 | 2024-04-02 | 南京信息工程大学 | 一种基于denclue聚类算法的雷暴识别方法 |
CN117809192B (zh) * | 2024-03-01 | 2024-04-26 | 南京信息工程大学 | 一种基于denclue聚类算法的雷暴识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113204903B (zh) | 2022-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113204903B (zh) | 一种预测雷电的方法 | |
CN113902897B (zh) | 目标检测模型的训练、目标检测方法、装置、设备和介质 | |
CN108732550B (zh) | 用于预测雷达回波的方法和装置 | |
CN111627015B (zh) | 小样本缺陷识别方法、装置、设备及存储介质 | |
Zhang et al. | Weather visibility prediction based on multimodal fusion | |
CN110705115A (zh) | 一种基于深度置信网络的气象预报方法和*** | |
CN112287983B (zh) | 一种基于深度学习的遥感图像目标提取***和方法 | |
CN108986142A (zh) | 基于置信图峰值旁瓣比优化的遮挡目标跟踪方法 | |
CN116310850B (zh) | 基于改进型RetinaNet的遥感图像目标检测方法 | |
CN115937659A (zh) | 基于Mask-RCNN的室内复杂环境下多目标检测方法 | |
CN115861756A (zh) | 基于级联组合网络的大地背景小目标识别方法 | |
CN117523420B (zh) | 一种基于雷达产品数据的闪电落区识别方法及*** | |
CN117671480A (zh) | 基于视觉大模型的滑坡自动识别方法、***和计算机设备 | |
CN109636194B (zh) | 一种输变电项目重大变动多源协同检测方法与*** | |
CN116580123A (zh) | 一种区域大气环境光化学网格模型图像生成*** | |
CN115456957A (zh) | 一种全尺度特征聚合的遥感影像变化检测的方法 | |
Socaci et al. | XNow: A deep learning technique for nowcasting based on radar products’ values prediction | |
CN115641498A (zh) | 基于空间多尺度卷积神经网络的中期降水预报后处理订正方法 | |
CN116263735A (zh) | 神经网络的鲁棒性评估方法、装置、设备及存储介质 | |
CN113989720A (zh) | 目标检测方法、训练方法、装置、电子设备及存储介质 | |
CN112418730A (zh) | 一种运输***的响应指标预估模型的构建方法 | |
Ma et al. | UAV Target Tracking and Detection Based on Faster R-CNN Improved Networks | |
CN117367435B (zh) | 疏散路径规划方法、装置、设备及存储介质 | |
CN116469013B (zh) | 道路积水预测方法、装置、计算机设备及存储介质 | |
CN117131783B (zh) | 基于多模态学习的输电线路风险预测模型、方法和*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |