(12)发明专利申请
(10)申请公布号 CN 110927725 A(43)申请公布日 2020.03.27
(21)申请号 201911198645.X(22)申请日 2019.11.30
(71)申请人 兰州大学
地址 730000 甘肃省兰州市城关区天水南
路222号(72)发明人 高学榜 谢莉 周俊 (74)专利代理机构 兰州智和专利代理事务所
(普通合伙) 62201
代理人 赵立权(51)Int.Cl.
G01S 13/95(2006.01)G01W 1/10(2006.01)
权利要求书2页 说明书5页 附图2页
CN 110927725 A(54)发明名称
一种基于气象雷达的雷电预警和监测方法(57)摘要
本发明提供了一种基于气象雷达的雷电预警和监测方法，主要包括如下步骤：S1.设雷电监测参数A为一个关于W波段雷达回波功率和S波段雷达反射率因子随雷电发展时间变化的函数；S2.计算得到雷暴云中颗粒系统携带电荷达到临界值时的折射率mwt'；S3.设雷电发生预警值为A'，雷电发生预警值A'可以由云颗粒折射率mwt'计算给出，当雷电监测参数A趋近于A'时，提示雷电即将发生。该方法基于气象雷达系统所提出，与传统监测雷电的技术相比，监测雷电的监测半径更大，可以增大到430km。此外，由于气象雷达系统的普及程度广泛，所以本方案实施更为方便、使用成本更低。
CN 110927725 A
权 利 要 求 书
1/2页
1.一种基于气象雷达的雷电预警和监测方法，其特征在于，包括如下步骤：S1.设雷电监测参数为A，且将雷电监测参数A设置为：
式中，为W波段雷达初始探测时刻雷暴云中颗粒的折射率且
mwt为W波段雷达此后每一时刻探测得到的雷暴云中颗粒的折射率，PWr(t0)为W波段雷达初始探测时刻接收到的回波功率，PWr(t)为W波段雷达此后每一时刻接收到的回波功率，z(λs,t0)为同时段S波段雷达初始探测时刻的反射率因子，z(λs,t)为同时段S波段雷达此后每一时刻的反射率因子；
S2.根据雷达测量到的当前雷暴云中颗粒系统的粒径分布和对应该粒径颗粒携带电荷的临界值，计算得到雷暴云中颗粒系统携带电荷达到临界值时的折射率mwt'；
S3.设雷电发生预警值为A'，将mwt＝mwt'带入步骤S1中函数中计算即可
得到雷电发生预警值A'；当雷电监测参数趋近于A'时，提示雷电即将
发生。
2.根据权利要求1所述的一种基于气象雷达的雷电预警和监测方法，其特征在于，步骤S1中，根据雷达方程可得：
式中，λ为雷达发射的电磁波的波长，R为雷达波探测距离，L为馈线损耗，Pr为雷达接收功率，Pt为雷达发射功率，h为雷达波有效照射深度，G是雷达天线增益，和θ分别为水平和垂直波束宽度，
为S波段雷达初始探测时刻雷暴云中颗粒的折射率，mst为S波段雷达此
后每一时刻雷暴云中颗粒的折射率，且
3.根据权利要求1所述的一种基于气象雷达的雷电预警和监测方法，其特征在于，步骤
2
CN 110927725 A
权 利 要 求 书
2/2页
S1中，t＝t0+△t，其中，t0为初始探测时刻，△t与雷达波的脉冲时间一致。
4.根据权利要求1所述的一种基于气象雷达的雷电预警和监测方法，其特征在于，步骤S2中，
式中，εεr为雷暴云中带电颗粒的介电常数，v为雷达初始探测时刻雷暴云中带电颗粒的介电常数，a为颗粒半径，η'为雷暴云中颗粒表面携带电荷的临界值，e是单个电子的电荷，
-12
me为电子的质量，εε，ω为W波段雷达电磁波的角0是真空介电常数且0＝8.854187817×10频率，γs为阻尼系数。
5.根据权利要求1所述的一种基于气象雷达的雷电预警和监测方法，其特征在于，步骤S2中，雷达初始探测时刻选择雷暴云初始形成不带电状态下的某一时刻。
3
CN 110927725 A
说 明 书
一种基于气象雷达的雷电预警和监测方法
1/5页
技术领域
[0001]本发明属于雷电预警与监测领域，具体涉及一种基于气象雷达的雷电预警和监测方法。
背景技术
[0002]雷电是一种常见的天气现象，由于雷电发生时会产生巨大的电脉冲辐射和回击电流，所以雷电被认为是对人类生产活动产生灾害的自然灾害之一。因此对雷电的准确预警直接影响到空中飞行器和微电子设备的安全应用，进而对雷电预警工作的研究越来越受到人们的重视。目前为止，一般用大气电场仪来监测和预警雷电，大气电场仪通过测量云雨颗粒起电过程中大气电场强度和极性的变化预警雷电的发展情况，它的优点在于可以全程监测雷电的发生、发展和消失，但是，大气电场仪的监测范围小，一般只能监测半径为8-9km范围内的雷电。
[0003]近年来，随着新一代天气雷达系统的建设，各个不同波段的气象雷达在全国范围内被布局。此外，与传统雷达相比，新布局的气象雷达监测准确度提高了很多。所以本发明基于广泛分布的气象雷达提出了一种新的预警和监测雷电的方法。
发明内容
[0004]本发明的目的在于针对现有技术存在的问题，提供一种监测半径更大、监测精度更高的雷电预警和监测方法。[0005]为此，本发明采用如下技术方案：
[0006]一种基于气象雷达的雷电预警和监测方法，包括如下步骤：[0007]S1.设雷电监测参数为A，且将雷电监测参数A设置为：
[0008]
[0009]式中，为W波段雷达初始探测时刻雷暴云中颗粒的折射率且mwt为W波段雷达此后每一时刻探测得到的雷暴云中颗粒的折射率，
PWr(t0)为W波段雷达初始探测时刻接收到的回波功率，PWr(t)为W波段雷达此后每一时刻接收到的回波功率，z(λz(λs,t0)为同时段S波段雷达初始探测时刻的反射率因子，s,t)为同时段S波段雷达此后每一时刻的反射率因子；
[0010]S2.根据雷达测量到的当前雷暴云中颗粒系统的粒径分布和对应该粒径颗粒携带电荷的临界值，计算得到雷暴云中颗粒系统携带电荷达到临界值时的折射率mwt'；
4
CN 110927725 A
说 明 书
2/5页
[0011]S3.设雷电发生预警值为A'，将mwt＝mwt'带入步骤S1中函数中计算
即可得到雷电发生预警值A'；当雷电监测参数即将发生。
[0012]进一步地，步骤S1中，根据雷达方程可得：
趋近于A'时，提示雷电
[0013]
[0014]
[0015]
式中，λ为雷达发射的电磁波的波长，R为雷达波探测距离，L为馈线损耗，Pr为雷达
为S波段雷达初始探测时刻雷暴云中颗粒的折射率，mst为S波段雷
接收功率，Pt为雷达发射功率，h为雷达波有效照射深度，G是雷达天线增益，和θ分别为水平和垂直波束宽度，
达此后每一时刻雷暴云中颗粒的折射率，且
进一步地，步骤S1中，t＝t0+△t，其中，t0为初始探测时刻，△t与雷达波的脉冲时间一致。
[0017]进一步地，步骤S2中，
[0018][0016]
式中，εεr为雷暴云中带电颗粒的介电常数，v为雷达初始探测时刻雷暴云中带电颗粒的介电常数，a为颗粒半径，η'为雷暴云中颗粒表面携带电荷的临界值，e是单个电子的电
-12
荷，me为电子的质量，εε，ω为W波段雷达电磁波0是真空介电常数且0＝8.854187817×10的角频率，γs为阻尼系数。[0020]进一步地，步骤S2中，雷达初始探测时刻选择雷暴云初始形成不带电状态下的某一时刻。
[0021]雷电是雷暴云中颗粒表面电荷聚集到一定程度的放电现象，基于大量的雷电统计数据，当雷暴云中颗粒表面电荷达到某一临界值时雷电发生。进一步研究发现，雷暴云中颗粒表面电荷的累积与雷暴云中颗粒的折射率有关，即
(式中，η为雷暴云中颗粒表面携带的电荷值，其余同上)，随着雷暴云中颗粒表面电荷的逐
5
[0019]
CN 110927725 A
说 明 书
3/5页
渐累积，颗粒的折射率也逐渐增加。因此，可以通过雷暴云中颗粒折射率的变化给定雷电监测参数A，并且随着雷暴云中颗粒表面携带电荷逐渐增大，监测参数A也逐渐增大，而当雷暴云中颗粒表面电荷达到某一临界值时，参数A达到最大值。故可以将该值作为雷电预警值A’，当雷电监测参数A趋近于A'时，提示雷电即将发生。
[0022]而气象雷达接收到的回波功率又与雷暴云中颗粒的折射率有关，因此，本发明通过气象雷达的相关参数进行计算即可得到雷电监测参数A。具体过程如下：[0023]由雷达方程可知，W波段雷达初始探测时刻接收到的回波功率PWr(t0)和此后每一时刻接收到的回波功率PWr(t)如下：
[0024]
[0025][0026]
从而，
[0027]
[0028]式中，为W波段雷达初始探测时刻获得的反射率因子和此后每一时刻反
射率因子的比值K，而根据不同波段雷达初始探测时刻获得的反射率因子和此后每一时刻反射率因子的比值相同的特性，以及雷暴云中颗粒表面电荷的累积容易影响W波段后向散射的强度，但不影响S波段后向散射的强度这一特点，可选择同时段S波段雷达初始探测时刻获得的反射率因子和此后每一时刻反射率因子的比值用于进行计算，即得：
[0029]
[0030]
另外，S波段和W波段雷达是目前监测云雨天气、雷暴天气的主要监测工具。因此，在本发明雷电监测和预警方法中选用这两种波段的雷达，可使监测更为方便和低成本。而且雷暴云中颗粒表面电荷对不同波段的雷达波的后向散射增强程度不同，选取的这两种波段雷达波长相差较大，能体现出带电影响的差异。
[0031]本发明雷电预警和监测方法的主要优点在于：[0032](1)该方法基于气象雷达系统，与传统的监测雷电技术相比，监测雷电的监测半径更大，可以增大到430km；[0033](2)由于气象雷达系统的普及程度广泛，所以本方案实施方便、使用成本低；
6
CN 110927725 A[0034]
说 明 书
4/5页
(3)由于气象雷达系统精度较高，所以本方案预警和监测雷电时间长，方位和距离
精确；
附图说明
[0035]图1为本发明雷电预警和监测的实施流程图；
[0036]图2为一次雷电发生过程中雷暴云颗粒系统的粒径分布图；[0037]图3为不同粒径下雷暴云中单颗粒携带的电荷临界值示意图。具体实施方式
[0038]下面结合具体实施例对本发明方案进行进一步说明。[0039]一种基于气象雷达的雷电预警和监测方法，包括如下步骤：[0040]S1.设雷电监测参数为A，且将雷电监测参数A设置为：
[0041]
[0042]式中，为W波段雷达初始探测时刻雷暴云中颗粒的折射率且mwt为W波段雷达此后每一时刻探测得到的雷暴云中颗粒的折射率，
PWr(t0)为W波段雷达初始探测时刻接收到的回波功率，PWr(t)为W波段雷达此后每一时刻接收到的回波功率，z(λz(λs,t0)为同时段S波段雷达初始探测时刻的反射率因子，s,t)为同时段S波段雷达此后每一时刻的反射率因子；且根据雷达方程可得：
[0043]
[0044]
[0045]
式中，λ为雷达发射的电磁波的波长，R为雷达波探测距离，L为馈线损耗，Pr为雷达
为S波段雷达初始探测时刻雷暴云中颗粒的折射率，mst为S波段雷
接收功率，Pt为雷达发射功率，h为雷达波有效照射深度，G是雷达天线增益，和θ分别为水平和垂直波束宽度，
达此后每一时刻雷暴云中颗粒的折射率，且
S2.根据雷达测量到的当前雷暴云中颗粒系统的粒径分布和对应该粒径颗粒携带电荷的临界值，计算得到雷暴云中颗粒系统携带电荷达到临界值时的折射率mwt'；且
[0046]
7
CN 110927725 A
说 明 书
5/5页
[0047][0048]
式中，εεr为雷暴云中带电颗粒的介电常数，v为雷达初始探测时刻雷暴云中带电颗
粒的介电常数，a为颗粒半径，η'为雷暴云中颗粒表面携带电荷的临界值，e是单个电子的电
-12
荷，me为电子的质量，εε，ω为W波段雷达电磁波0是真空介电常数且0＝8.854187817×10
成比例，为了计算方便，令
表示约化普朗克
的角频率，γs为阻尼系数。γs与弛豫时间相反，且与
kB＝1.38×10-23JK-1，为玻尔兹曼常数；
常数，T表示颗粒的温度，单位是开尔文(K)。
[0049]S3.设雷电发生预警值为A'，将mwt＝mwt'带入步骤S1中函数中计算
即可得到雷电发生预警值A'；当雷电监测参数趋近于A'时，提示雷电
即将发生。
[0050]进一步地，步骤S1中，t＝t0+△t，其中，t0为初始探测时刻，△t与雷达波的脉冲时间一致，本实施例取值为△t＝2μs。
[0051]雷电监测参数A选择上述函数形式，是因为与雷暴云中颗粒系统的粒
径分布和对应该粒径云颗粒携带电荷的临界值有关，并可以通过气象雷达的回波信息以及图3给出的相关数据进行计算。图3是Klett等人通过大量的研究发现，在雷暴天气中不同粒径云颗粒表面的累积电荷随着时间的变化会达到一个极限值，即不同粒径的云颗粒达到极限值时累积的电荷量数据(Pruppacher HR,Klett JD.In:Microphysics of clouds and precipitation.the Netherlands:Kluwer Academic Publishers；2003.)。本发明计算过程引用了上述研究成果。[0052]以一次雷电发生为例，通过气象雷达观测到当前雷暴云颗粒系统的粒径分布如图2所示，从该粒径分布图可以得出当前云颗粒系统的等效粒径为1.1μm，查阅图3可知对应该粒径的单颗粒携带电荷的临界数值为1.5×106。将上述数值带入
中可得出临界折射率mwt'＝1.89+125.73i；进一步
计算可得雷电发生预警值A'＝5.396。当雷电监测参数A趋近于该值时，提示雷电即将发生。
8
CN 110927725 A
说 明 书 附 图
1/2页
图1
9
CN 110927725 A
说 明 书 附 图
2/2页
图2
图3
10
因篇幅问题不能全部显示，请点此查看更多更全内容