CN106353604A - 一种基于大气电场背景噪声自动追踪的雷电预警方法 - Google Patents

一种基于大气电场背景噪声自动追踪的雷电预警方法 Download PDF

Info

Publication number
CN106353604A
CN106353604A CN201610775673.3A CN201610775673A CN106353604A CN 106353604 A CN106353604 A CN 106353604A CN 201610775673 A CN201610775673 A CN 201610775673A CN 106353604 A CN106353604 A CN 106353604A
Authority
CN
China
Prior art keywords
power
electric field
atmospheric electric
frame
background noise
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
Application number
CN201610775673.3A
Other languages
English (en)
Other versions
CN106353604B (zh
Inventor
陶汉涛
谷山强
吴大伟
王海涛
冯万兴
赵淳
郭钧天
卢恩泽
周自强
张磊
陈玥
姜志博
何君
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan NARI Ltd
Original Assignee
Wuhan NARI Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan NARI Ltd filed Critical Wuhan NARI Ltd
Priority to CN201610775673.3A priority Critical patent/CN106353604B/zh
Publication of CN106353604A publication Critical patent/CN106353604A/zh
Application granted granted Critical
Publication of CN106353604B publication Critical patent/CN106353604B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08Measuring electromagnetic field characteristics
    • G01R29/0807Measuring electromagnetic field characteristics characterised by the application
    • G01R29/0814Field measurements related to measuring influence on or from apparatus, components or humans, e.g. in ESD, EMI, EMC, EMP testing, measuring radiation leakage; detecting presence of micro- or radiowave emitters; dosimetry; testing shielding; measurements related to lightning
    • G01R29/0842Measurements related to lightning, e.g. measuring electric disturbances, warning systems

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Relating To Insulation (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于大气电场背景噪声自动追踪的雷电预警方法,主要包括:S1:输入关注区域实测的大气电场数据;S2:等时间间隔截取输入数据后进行分帧处理;S3:以帧为单位求取每一时间段大气电场的功率谱;S4:根据所求大气电场功率谱的特征估算背景噪声功率谱;S5:计算每一帧大气电场功率与背景噪声功率的比值进行雷电预警。该方法可自动追踪非雷暴时刻变化缓慢的背景电场,消除了噪声对雷暴前期电场突变判断的影响,为电场环境复杂、气候差异较大的区域提供了一种适应性强、实用性高的雷电预警方法。

Description

一种基于大气电场背景噪声自动追踪的雷电预警方法
技术领域
本发明属于雷电预警技术领域,尤其涉及一种基于大气电场背景噪声自动追踪的雷电预警方法。
背景技术
雷电是自然界中频繁发生的一种剧烈放电现象,给人类的劳动生产以及生命安全带来了极大的危害。据统计,我国每年雷击次数高达1000万余次,因雷击造成的人员伤亡3000~4000人,直接经济损失约100亿元人民币。现有的防雷措施有多种,如安装避雷器(针、线)、降低接地电阻等,但这些被动式的雷电防护措施已无法满足人类生产提出的防雷新要求。因此,有必要开拓以雷电预警为基础的主动性雷电防护新思路,在雷暴临近前采取有效的处理措施或制定合理的应急预案以减小户外设施的雷击故障率以及因雷击造成的经济损失和人员伤亡。
国内外研究学者对于雷电预警已开展大量的研究,主要是依据雷暴前大气电场的幅值特征,辅以雷电定位、气象雷达、卫星云图等多维雷电观测资料进行雷电活动预警。如《一种雷电预警装置及雷电预警的实现方法》(201210557928.0)是通过判断雷电波和电场是否均达到预警设定值来启动报警程序;《雷电预警方法及装置》(201010567357.X)公开了一种根据计算得到的雷电源距离及雷电发生的频数进行不同等级雷电报警的方法;《适用于输电线路施工作业的雷电预警系统及预警方法》(201310726107.X)综合了雷达回波、雷电定位以及大气电场信息实现不同等级的雷电预警。上述雷电预警方法通过设置大气电场预警阈值来完成雷暴天气的判断,简单、可行,在工程应用中具有一定的效果。但所设阈值需要前期对雷暴天气与非雷暴天气的大气电场数据进行大量的对比分析,而且预警效果受大气电场仪所处环境较大,导致实际应用中误报、漏报较多。
有鉴于此,有必要提供一种基于大气电场背景噪声自动追踪的雷电预警方法,以解决上述问题。
发明内容
针对上述现有技术存在的不足,本发明提供一种基于大气电场背景噪声自动追踪的雷电预警方法,依据雷暴天气与非雷暴天气大气电场的变化特征,在求取大气电场功率谱的基础上,采用自适应噪声谱估计方法计算大气电场背景噪声的功率谱,并以大气电场功率与背景噪声功率的比值作为雷电预警的依据。
为实现此目的,本发明所采用的技术方案是:一种基于大气电场背景噪声自动追踪的雷电预警方法,其特征在于,包括以下步骤:
S1:输入关注区域实测的大气电场数据;
S2:等时间间隔截取输入数据后进行分帧处理;
S3:以帧为单位求取每一时间段大气电场的功率谱;
S4:根据所求大气电场功率谱的特征估算背景噪声功率谱;
S5:计算每一帧大气电场功率与背景噪声功率的比值进行雷电预警。
作为进一步可选的技术方案,上述方案中,步骤S4中根据所求大气电场功率谱的特征,采用自适应噪声谱估计方法估算背景噪声的功率谱,包括:
a) 采用一阶递归公式计算大气电场的平滑功率谱;
b) 计算平滑功率谱中的功率基准值;
c) 计算每一帧平滑功率与功率基准值的比值;
d) 根据上述比值计算背景噪声的平滑系数;
e) 根据上述系数计算背景噪声功率谱。
作为进一步可选的技术方案,上述方案中,所述的自适应噪声谱估计方法,功率基准值计算方法为:若当前帧的平滑功率小于前一帧的功率基准值,则将当前帧的平滑功率作为当前帧的功率基准值;否则,对前一帧的功率基准值采用一阶函数进行自适应调整,比较调整后功率与当前帧平滑功率的大小,取较小值作为当前帧的功率基准值。
作为进一步可选的技术方案,上述方案中,所述的功率基准值计算方法,实测大气电场在当前时间段内第一帧的功率基准值为上一时间段尾部平滑功率的最小值,若为初次运行雷电预警程序,则将第一帧的平滑功率作为第一帧的功率基准值。
与现有技术相比,本发明具有如下的有益效果:该方法可自动追踪非雷暴时刻变化缓慢的背景电场,消除了噪声对雷暴前期电场突变判断的影响,为电场环境复杂、气候差异较大的区域提供了一种适应性强、实用性高的雷电预警方法。
附图说明
图1是本发明一种基于大气电场背景噪声自动追踪的雷电预警方法的流程图。
图2是本发明实施例的非雷暴天气雷电预警过程图。
图3是本发明实施例的雷暴天气雷电预警过程图。
具体实施方式
为了更好地理解本发明,下面结合实施例进一步阐明本发明的内容,但本发明的内容不仅仅局限于下面的实施例。本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样在本申请所列权利要求书限定范围之内。
通常情况下,非雷暴天气的大气电场变化缓慢,幅值也较小,其统计变化范围为0.019~0.310 kV/m;然而雷暴发生前,大气电场即开始剧烈抖动,幅值也增大,可达到8kV/m,甚至更高。根据这一特点,再考虑环境中噪声电场的影响,可提出一种基于大气电场背景噪声自动追踪的雷电预警方法。该方法在大气电场的实时监测过程中,依据非雷暴时刻大气电场变化缓慢的特点自动追踪大气环境中的背景电场,再根据雷暴前期大气电场幅值突变的特征进行雷电预警,为电场环境复杂、气候差异较大的区域提供了一种适应性强、实用性高的雷电预警方法。
如图1所示,在本发明的一个具体实施例中,依据关注区域实测的大气电场数据进行雷电预警,具体包括以下步骤:
S1:采用已标定好的大气电场仪实时探测关注区域的大气电场,并将探测到的大气电场数据作为雷电预警程序的输入;
S2:等时间间隔截取输入的实测大气电场数据,对截取后的大气电场数据进行分帧处理;
S3:对分帧处理后的大气电场数据采用傅里叶变换依次求取当前时间段内每一帧大气电场的幅值谱,进而获得当前时间段的大气电场功率谱;
S4:采用自适应噪声谱估计方法估算大气电场背景噪声的功率谱,主要包括:
a) 采用一阶递归公式计算大气电场平滑功率谱P(λ, k)为
b) 计算平滑功率谱中的功率基准值P ref(λ, k)有
c) 计算每一帧平滑功率与功率基准值的比值S r(λ, k)有
d) 根据上述比值计算背景噪声的平滑系数α s(λ, k)的方法为
① 判断是否为雷暴前期的大气电场
② 计算为雷暴前期大气电场的概率S p(λ, k)
③ 计算背景噪声的更新系数α s(λ, k)
e) 根据上述系数计算背景噪声功率谱N(λ, k)有
其中, Y(λ, k)为傅里叶变换后得到的大气电场幅值谱;λ为分帧后的帧序号;k为傅里叶变换后的频段序号;δ(k)为与频段有关的噪声更新判断阈值;η、γ、β、α pα d为常数。
S5:依次计算每一帧大气电场功率与背景噪声功率的比值,若大于预设阈值,则发出相应的雷电预警信息,否则进行下一帧比值的判断。完成当前时间段最后一帧比值的判断后,等待下一时间段大气电场数据的输入,重复步骤S1~S5。
本说明书未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (4)

1.一种基于大气电场背景噪声自动追踪的雷电预警方法,其特征在于,包括以下步骤:
S1:输入关注区域实测的大气电场数据;
S2:等时间间隔截取输入数据后进行分帧处理;
S3:以帧为单位求取每一时间段大气电场的功率谱;
S4:根据所求大气电场功率谱的特征估算背景噪声功率谱;
S5:计算每一帧大气电场功率与背景噪声功率的比值进行雷电预警。
2.根据权利要求1所述的基于大气电场背景噪声自动追踪的雷电预警方法,其特征在于,步骤S4中根据所求大气电场功率谱的特征,采用自适应噪声谱估计方法估算背景噪声的功率谱,包括:
a) 采用一阶递归公式计算大气电场的平滑功率谱;
b) 计算平滑功率谱中的功率基准值;
c) 计算每一帧平滑功率与功率基准值的比值;
d) 根据上述比值计算背景噪声的平滑系数;
e) 根据上述系数计算背景噪声功率谱。
3.根据权利要求2所述的基于小型雷达的输电线路雷电监测预警方法,其特征在于,所述的自适应噪声谱估计方法,功率基准值计算方法为:若当前帧的平滑功率小于前一帧的功率基准值,则将当前帧的平滑功率作为当前帧的功率基准值;否则,对前一帧的功率基准值采用一阶函数进行自适应调整,比较调整后功率与当前帧平滑功率的大小,取较小值作为当前帧的功率基准值。
4.根据权利要求3所述的基于小型雷达的输电线路雷电监测预警方法,其特征在于,所述的功率基准值计算方法,实测大气电场在当前时间段内第一帧的功率基准值为上一时间段尾部平滑功率的最小值,若为初次运行雷电预警程序,则将第一帧的平滑功率作为第一帧的功率基准值。
CN201610775673.3A 2016-08-31 2016-08-31 一种基于大气电场背景噪声自动追踪的雷电预警方法 Active CN106353604B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610775673.3A CN106353604B (zh) 2016-08-31 2016-08-31 一种基于大气电场背景噪声自动追踪的雷电预警方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610775673.3A CN106353604B (zh) 2016-08-31 2016-08-31 一种基于大气电场背景噪声自动追踪的雷电预警方法

Publications (2)

Publication Number Publication Date
CN106353604A true CN106353604A (zh) 2017-01-25
CN106353604B CN106353604B (zh) 2019-06-07

Family

ID=57857043

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610775673.3A Active CN106353604B (zh) 2016-08-31 2016-08-31 一种基于大气电场背景噪声自动追踪的雷电预警方法

Country Status (1)

Country Link
CN (1) CN106353604B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110297284A (zh) * 2019-07-30 2019-10-01 南京信息工程大学 一种基于三维大气电场仪的雷暴云移动路径跟踪方法
CN110927725A (zh) * 2019-11-30 2020-03-27 兰州大学 一种基于气象雷达的雷电预警和监测方法
CN113985148A (zh) * 2021-11-04 2022-01-28 湖南省气象灾害防御技术中心 一种雷电预报方法及雷电预报设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103472342A (zh) * 2013-09-25 2013-12-25 中国气象科学研究院 雷电监测预警方法及系统
CN103745580A (zh) * 2013-12-25 2014-04-23 国网电力科学研究院武汉南瑞有限责任公司 适用于输电线路施工作业的雷电预警系统及预警方法
CN105353225A (zh) * 2015-10-13 2016-02-24 中国石油化工股份有限公司 可预测雷电运动趋势的雷电预警方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103472342A (zh) * 2013-09-25 2013-12-25 中国气象科学研究院 雷电监测预警方法及系统
CN103745580A (zh) * 2013-12-25 2014-04-23 国网电力科学研究院武汉南瑞有限责任公司 适用于输电线路施工作业的雷电预警系统及预警方法
CN105353225A (zh) * 2015-10-13 2016-02-24 中国石油化工股份有限公司 可预测雷电运动趋势的雷电预警方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
李颖等: "大气电场的FFT频谱分析及雷暴预报研究", 《气象科学》 *
王强等: "电场时序差分在雷电预警中的有效性分析", 《气象科学》 *
高菊霞等: "利用底面大气电场监测实现雷电预警", 《陕西气象》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110297284A (zh) * 2019-07-30 2019-10-01 南京信息工程大学 一种基于三维大气电场仪的雷暴云移动路径跟踪方法
CN110927725A (zh) * 2019-11-30 2020-03-27 兰州大学 一种基于气象雷达的雷电预警和监测方法
CN110927725B (zh) * 2019-11-30 2023-03-14 兰州大学 一种基于气象雷达的雷电预警和监测方法
CN113985148A (zh) * 2021-11-04 2022-01-28 湖南省气象灾害防御技术中心 一种雷电预报方法及雷电预报设备

Also Published As

Publication number Publication date
CN106353604B (zh) 2019-06-07

Similar Documents

Publication Publication Date Title
Wang et al. Rapid rises in the magnitude and risk of extreme regional heat wave events in China
CN104122488A (zh) 一种输电线路故障校核诊断方法
CN103236666B (zh) 输电线路差异化防雷方法
Chen et al. Research on the characteristics of urban rainstorm pattern in the humid area of Southern China: A case study of Guangzhou City
CN103337133A (zh) 基于识别预报技术的电网雷暴灾害预警系统及方法
CN104202719B (zh) 基于位置可信度的人数测定和人群态势监测方法及系统
CN106353604A (zh) 一种基于大气电场背景噪声自动追踪的雷电预警方法
CN105137286A (zh) 输电线路雷击监测装置及防雷水平评估的方法
CN106019287A (zh) 一种基于小型雷达的输电线路雷电监测预警方法
CN102446390A (zh) 一种对输电线路附近监控区域进行安全检测、预警的方法和系统
CN106157540B (zh) 一种基于雷电预警系统的雷电预警及雷云轨迹预测方法
CN107870275A (zh) 基于大数据的避雷器运行状态评价方法
CN109374995A (zh) 低压配电系统、线路老化和反窃电监测方法及对应系统
CN110147741A (zh) 一种用于电力网络管理的遥感森林树高提取方法
CN104793062A (zh) 一种基于农业物联网的田间终端设备防雷方法
CN113222260A (zh) 输电线路的雷击跳闸率预测模型
CN204089886U (zh) 基于气象信息的电网动态监控预警系统
CN106680574B (zh) 一种变电站设备过电压感知和数据处理方法
Su-Ping et al. Quality control and analysis of global gauge-based daily precipitation dataset from 1980 to 2009
CN107221882B (zh) 一种线路避雷器安装选点方法
CN115356525A (zh) 一种基于闪电定位技术的雷击接闪系统及其使用方法
CN111222662A (zh) 一种电网台风洪涝灾害预警方法和装置
CN205864628U (zh) 输电线路异物和缺失智能视频在线监测系统
CN206075498U (zh) 基于数值模拟天气预报的电网风害预警系统
CN108629078A (zh) 基于监控信息的输电网雷击灾害可靠性评估方法

Legal Events

Date Code Title Description
C06 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