CN115078848A - 基于闪电信号的电离层无源被动探测方法 - Google Patents
基于闪电信号的电离层无源被动探测方法 Download PDFInfo
- Publication number
- CN115078848A CN115078848A CN202210815772.5A CN202210815772A CN115078848A CN 115078848 A CN115078848 A CN 115078848A CN 202210815772 A CN202210815772 A CN 202210815772A CN 115078848 A CN115078848 A CN 115078848A
- Authority
- CN
- China
- Prior art keywords
- lightning
- value
- satellite
- model
- propagation
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 19
- 230000001902 propagating effect Effects 0.000 claims abstract description 9
- 239000006185 dispersion Substances 0.000 claims abstract description 4
- 239000005433 ionosphere Substances 0.000 claims description 23
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000000034 method Methods 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 6
- 230000000052 comparative effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/0807—Measuring electromagnetic field characteristics characterised by the application
- G01R29/0814—Field 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/0842—Measurements related to lightning, e.g. measuring electric disturbances, warning systems
-
- 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
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了基于闪电信号的电离层无源被动探测方法,具体涉及闪电信号探测技术领域,包括如下步骤,S1.根据IGRF模型确定卫星所在位置的磁倾角与磁场强度,S2.由磁离子理论中的A‑H(阿普顿‑哈特利)公式,即电磁波在等离子体中传播的色散关系式,计算路径L上的折射率分布,S3.IRI模型和IGRF模型共同给出了A‑H公式中的X,Y,Z的值,S4.电离层进行垂直方向上的网格化,S5.利用观测值对状态预测值进行更新,得到更新后的状态分析值,S6.高斯误差协方差模型假设背景场数据之间高斯相关,其协方差正比于观测值的平方,S7.对IRI模型给出的电子密度剖面进行了调整修正,利用闪电信号这一自然源发射的VLF电磁波,通过卫星接收闪电信号的时间来对电离层进行被动探测,可降低成本。
Description
技术领域
本发明涉及闪电信号探测领域,具体为基于闪电信号的电离层无源被动探测方法。
背景技术
闪电是大气中一种超长距离的强放电现象,是对流性天气系统发展到一定阶段的产物,有效监测闪电信号,对于研究闪电对强对流天气的指示作用,以及闪电对大气环境和气候的影响等具有指导意义。
专利号为CN110837006B的申请文件公开了一种基于星地同步观测对比的卫星闪电探测评估方法,其特征在于,包括如下步骤,S1:对一个典型区域在特定时间段内星地同步观测的闪电进行对比,以说明卫星与地基闪电探测结果的相关性和差异性;S2:为解决卫星与地基探测闪电信号不同表现形式的问题,根据地闪回击数据的时空粒度特征,从卫星探测的逐帧像元数据中提取一种新的卫星探测闪电数据单元spot,使其在时空粒度上与地基stroke具有可比性;S3:为解决卫星与地基观测闪电种类不完全一致的问题,找到卫星与地基可以同时观测闪电的特定位置;S4:对特定区域的卫星与地基闪电探测结果进行比对,定义二者的匹配度。
但是现有技术中闪电信号电离层的探测需要借助卫星发射电磁波,地面接收站接收电磁波,成本较高。
发明内容
本发明的目的在于提供基于闪电信号的电离层无源被动探测方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:
基于闪电信号的电离层无源被动探测方法,包括如下步骤:
S1.根据IGRF模型确定卫星所在位置的磁倾角与磁场强度,由此确定穿过卫星的磁力线路径,即卫星接收到的闪电VLF波在电离层中的传播路径L,假设传播路径中任一点的磁倾角为θ,磁场强度为B;
S2.由磁离子理论中的A-H(阿普顿-哈特利)公式,即电磁波在等离子体中传播的色散关系式,计算路径L上的折射率分布:
S3.IRI模型和IGRF模型共同给出了A-H公式中的X,Y,Z的值,求闪电VLF波在电离层中传播路径上任一点的折射率n,由折射率与波的传播速度的关系表达式:v=nrC,其中v为电磁波的传播速度,nr为折射率n的实部,C为光速;则波在路径上某一位置处的传播时间为:闪电VLF波在电离层中传播的总时间为:其中h1为电离层底部的高度,h为卫星高度,闪电VLF波在电离层中的传播时间与闪电信号的频率有关,可得闪电信号在电离层中传播的时频曲线;
S4.对电离层进行垂直方向上的网格化,则在一个网格内两个不同频率闪电信号的传播时间差为:其中dh为该网格垂直方向的长度,则两个不同频率闪电信号在电离层中传播至卫星处的总时间差为:其中h1为电离层底部的高度,h为卫星高度;
其中Nen为电离层剖面第n个网格上IRI模型给出的电子密度值,式中,Hk+1是k+1时刻的观测算子,表示观测值和状态值之间的函数关系:Hk+1=(dh1,dh2,...,dhn),其中dhn表示第n个网格的长度,则表示两个不同频率闪电信号在电离层中总的传播时间差;是k+1时刻卫星接收频率为f1和f2的闪电信号时间差的观测值,上标0表示观测,式中,Kk+1是K+1时刻的增益矩阵,Pk+1是k+1时刻状态预测值的误差协方差矩阵:
S6.高斯误差协方差模型假设背景场数据之间高斯相关,其协方差正比于观测值的平方,计算方法如下式所示:Bij=σiσjρij,其中Bij为在电离层剖面两个不同网格中频率f1和f2频率的闪电信号传播时间差的误差协方差,ρij为相关系数,ηb为比例系数,l(H)为高度方向的相关长度,Hij是网格点间在高度方向上的距离;
S7.利用k+1时刻观测值调制背景场,得到k+1时刻的两个不同频率闪电信号在网格中传播时间差的分析值,进而得到各网格点上的电子密度分析值,即对IRI模型给出的电子密度剖面进行了调整修正。
优选的,S4中采用卡尔曼滤波算法综合利用观测信息和预测模型与观测数据的误差统计特征对状态量进行估计,使状态估计值的误差最小,进而调整IRI模型给出的电离层电子密度剖面。
与现有技术相比,本发明的有益效果是:利用闪电信号这一自然源发射的VLF电磁波,通过卫星接收闪电信号的时间来对电离层进行被动探测,可降低成本。
附图说明
图1为对比例的探测原理图;
图2为本发明的原理图;
图3为本发明闪电信号在电离层中传播的时频曲线图;
图4为本发明实施例中电子密度剖面图;
图5为本发明卫星接收闪电信号的时频图;
图6为本发明实施例中修正后的电子密度剖面图;
图7为本发明更新前的状态预测值与更新后的分析值的误差分析图。
具体实施方式
对比例
如图1所示,基于Chapman函数的电离层TEC同化模型计算流程,主要包括4个步骤,
S1.利用背景场模型和计算时刻前两天的历史TEC数据,获取电离层特性参数η和α的值,输入参数包括年份、年积日、世界时、地理经纬度、太阳辐射强度I、电离层温度T;
式中,E(βc,λc)表示(βc,λc)处的TEC残差值,(m),m(cosθc)为完全正则化的非整阶缔合Legendre函数,Ckm和Skm为完全正则化的球冠谐函数系数,Kmax和M分别为球冠谐模型的最大阶数和最大次数
S3.利用当前η、α、I、T计算区域内所有格网点上的电离层TEC值,并对计算值利用S2中函数式进行残差改正得到新的TEC值,最后利用新的TEC值对η、α、I、T进行纠正;
S4.将纠正后的特性参数值输入背景场模型对电离层TEC值进行预测。
实施例
如图2-7所示,根据IGRF模型确定卫星所在位置(165km处)的磁倾角与磁场强度,由此确定了穿过卫星的磁力线路径,即卫星接收到的闪电VLF波在电离层中的传播路径L,并得到路径上任意位置处的磁倾角θ和磁场强度B;
根据IRI模型确定了闪电信号在电离层中传播路径L上K+1时刻任意位置的电子密度,电子密度剖面图如图4所示;
由磁离子理论中的A-H(阿普顿-哈特利)公式,即电磁波在等离子体中传播的色散关系式,计算路径L上的折射率分布;IRI模型和IGRF模型共同给出了A-H公式中的X,Y,Z的值,由此可求得闪电VLF波在电离层中传播路径上任一点的折射率n,由折射率与波的传播速度的关系表达式,v=nrC,其中v为电磁波的传播速度,nr为折射率n的实部,C为光速,则波在路径上某一位置处的传播时间为:对电离层进行垂直方向上的网格化,将电离层底部(65km)至卫星(165km)的纵向空间上等分为1000个网格,每个网格高度100m,并假设每一个网格内部电子密度,磁倾角,磁场强度均不变。则在一个网格内两个不同频率闪电信号(2kHz和4kHz)的传播时间差为:其中dh为该网格垂直方向的长度;则两个不同频率闪电信号在电离层中传播至卫星处的总时间差为:其中h1为电离层底部的高度,h为卫星高度;公式给出的是两个不同频率的闪电VLF波在电离层中传播时间差的预测模型,采用卡尔曼滤波算法综合利用观测信息和预测模型与观测数据的误差统计特征对状态量进行估计,使状态估计值的误差最小,进而调整IRI模型给出的电离层电子密度剖面。
其中Nen为电离层剖面第n个网格上IRI模型给出的电子密度值,式中,Hk+1是k+1时刻的观测算子,表示观测值和状态值之间的函数关系:Hk+,=(dh1,dh2,...,dhn)=(100,100...100),其中dhn表示第n个网格的长度,则表示两个不同频率闪电信号在电离层中总的传播时间差,卫星接收闪电信号的时频图图5所示;是k+1时刻卫星接收频率为f1和f2的闪电信号时间差的观测值,上标0表示观测,该观测值由卫星实际接收闪电信号的时频图得到:式中,Kk+1是K+1时刻的增益矩阵,Pk+1是k+1时刻状态预测值的误差协方差矩阵:高斯误差协方差模型假设背景场数据之间高斯相关,其协方差正比于观测值的平方,计算方法如下式所示:Bij=σiσjρij, 其中Bij为在电离层剖面两个不同网格中频率f1和f2频率的闪电信号传播时间差的误差协方差,ρij为相关系数,ηb为比例系数,l(H)为高度方向的相关长度,Hij是网格点间在高度方向上的距离;
利用k+1时刻观测值调制背景场,得到k+1时刻的两个不同频率闪电信号在网格中传播时间差的分析值,进而得到各网格点上的电子密度分析值,即对IRI模型给出的电子密度剖面进行了调整修正。修正后的电子密度剖面分析值如图6所示;更新前的状态预测值与更新后的分析值的误差分析如图7所示。
Claims (2)
1.基于闪电信号的电离层无源被动探测方法,其特征在于:包括如下步骤:
S1.根据IGRF模型确定卫星所在位置的磁倾角与磁场强度,由此确定穿过卫星的磁力线路径,即卫星接收到的闪电VLF波在电离层中的传播路径L,假设传播路径中任一点的磁倾角为θ,磁场强度为B;
S2.由磁离子理论中的A-H公式,即电磁波在等离子体中传播的色散关系式,计算路径L上的折射率分布:
S3.IRI模型和IGRF模型共同给出了A-H公式中的X,Y,Z的值,求闪电VLF波在电离层中传播路径上任一点的折射率n,由折射率与波的传播速度的关系表达式:v=nrC,其中v为电磁波的传播速度,nr为折射率n的实部,C为光速;则波在路径上某一位置处的传播时间为:闪电VLF波在电离层中传播的总时间为:其中h1为电离层底部的高度,h为卫星高度,闪电VLF波在电离层中的传播时间与闪电信号的频率有关,可得闪电信号在电离层中传播的时频曲线;
S4.对电离层进行垂直方向上的网格化,则在一个网格内两个不同频率闪电信号的传播时间差为:其中dh为该网格垂直方向的长度,则两个不同频率闪电信号在电离层中传播至卫星处的总时间差为:其中h1为电离层底部的高度,h为卫星高度;
其中Nen为电离层剖面第n个网格上IRI模型给出的电子密度值,式中,Hk+1是k+1时刻的观测算子,表示观测值和状态值之间的函数关系:Hk+1=(dh1,dh2,...,dhn),其中dhn表示第n个网格的长度,则表示两个不同频率闪电信号在电离层中总的传播时间差;是k+1时刻卫星接收频率为f1和f2的闪电信号时间差的观测值,上标0表示观测,式中,Kk+1是K+1时刻的增益矩阵,Pk+1是k+1时刻状态预测值的误差协方差矩阵:
S6.高斯误差协方差模型假设背景场数据之间高斯相关,其协方差正比于观测值的平方,计算方法如下式所示:Bij=σiσjρij,其中Bij为在电离层剖面两个不同网格中频率f1和f2频率的闪电信号传播时间差的误差协方差,ρij为相关系数,ηb为比例系数,l(H)为高度方向的相关长度,Hij是网格点间在高度方向上的距离;
S7.利用k+1时刻观测值调制背景场,得到k+1时刻的两个不同频率闪电信号在网格中传播时间差的分析值,进而得到各网格点上的电子密度分析值,即对IRI模型给出的电子密度剖面进行了调整修正。
2.根据权利要求1所述的基于闪电信号的电离层无源被动探测方法,其特征在于:S4中采用卡尔曼滤波算法综合利用观测信息和预测模型与观测数据的误差统计特征对状态量进行估计,使状态估计值的误差最小,进而调整IRI模型给出的电离层电子密度剖面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210815772.5A CN115078848B (zh) | 2022-07-12 | 2022-07-12 | 基于闪电信号的电离层无源被动探测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210815772.5A CN115078848B (zh) | 2022-07-12 | 2022-07-12 | 基于闪电信号的电离层无源被动探测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115078848A true CN115078848A (zh) | 2022-09-20 |
CN115078848B CN115078848B (zh) | 2024-04-09 |
Family
ID=83259172
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210815772.5A Active CN115078848B (zh) | 2022-07-12 | 2022-07-12 | 基于闪电信号的电离层无源被动探测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115078848B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006079041A1 (en) * | 2005-01-24 | 2006-07-27 | Quasar Federal Systems, Inc. | Integrated sensor system monitoring and characterizing lightning events |
JP2008122100A (ja) * | 2006-11-08 | 2008-05-29 | Toshiba Corp | 電離層電子密度算出装置 |
RU2011128236A (ru) * | 2011-07-07 | 2013-02-10 | Александр Васильевич Тертышников | Способ зондирования ионосферы, тропосферы, геодвижений и комплекс для его реализации |
CN103278704A (zh) * | 2013-05-30 | 2013-09-04 | 重庆绿色智能技术研究院 | 一种基于北斗的三维闪电探测系统和方法 |
CN104793204A (zh) * | 2015-04-21 | 2015-07-22 | 武汉大学 | 一种中间层-电离层综合探测系统 |
CN111007489A (zh) * | 2019-12-05 | 2020-04-14 | 武汉大学 | 一种浮标式电离层双程斜向返回探测系统及方法 |
CN111007490A (zh) * | 2019-12-05 | 2020-04-14 | 武汉大学 | 一种基于浮标地理信息的天波超视距雷达坐标配准方法 |
CN112649899A (zh) * | 2020-11-19 | 2021-04-13 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种全球电离层数据同化和预报方法 |
CN114167505A (zh) * | 2021-11-22 | 2022-03-11 | 武汉大学 | 基于LoranC甚低频信号的低电离层D层探测系统及方法 |
CN114384564A (zh) * | 2021-12-31 | 2022-04-22 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种基于多源数据驱动的电离层层析成像方法 |
-
2022
- 2022-07-12 CN CN202210815772.5A patent/CN115078848B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006079041A1 (en) * | 2005-01-24 | 2006-07-27 | Quasar Federal Systems, Inc. | Integrated sensor system monitoring and characterizing lightning events |
JP2008122100A (ja) * | 2006-11-08 | 2008-05-29 | Toshiba Corp | 電離層電子密度算出装置 |
RU2011128236A (ru) * | 2011-07-07 | 2013-02-10 | Александр Васильевич Тертышников | Способ зондирования ионосферы, тропосферы, геодвижений и комплекс для его реализации |
CN103278704A (zh) * | 2013-05-30 | 2013-09-04 | 重庆绿色智能技术研究院 | 一种基于北斗的三维闪电探测系统和方法 |
CN104793204A (zh) * | 2015-04-21 | 2015-07-22 | 武汉大学 | 一种中间层-电离层综合探测系统 |
CN111007489A (zh) * | 2019-12-05 | 2020-04-14 | 武汉大学 | 一种浮标式电离层双程斜向返回探测系统及方法 |
CN111007490A (zh) * | 2019-12-05 | 2020-04-14 | 武汉大学 | 一种基于浮标地理信息的天波超视距雷达坐标配准方法 |
CN112649899A (zh) * | 2020-11-19 | 2021-04-13 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种全球电离层数据同化和预报方法 |
CN114167505A (zh) * | 2021-11-22 | 2022-03-11 | 武汉大学 | 基于LoranC甚低频信号的低电离层D层探测系统及方法 |
CN114384564A (zh) * | 2021-12-31 | 2022-04-22 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种基于多源数据驱动的电离层层析成像方法 |
Non-Patent Citations (3)
Title |
---|
TIAN XIANG, MORAN LIU , SHIMIN HE AND CHEN ZHOU: "Assimilation and Inversion of Ionospheric Electron Density Data Using Lightning Whistlers", REMOTE SENSING, vol. 15, no. 12, 10 July 2023 (2023-07-10) * |
周晨;赵正予;杨国斌;邓峰;: "电离层高频信道互易性研究:中纬度电离层斜向探测实验", 电子与信息学报, no. 01, 15 January 2011 (2011-01-15) * |
胡小希,周晨,赵军,刘祎,刘默然,赵正予: "基于优化神经网络算法的电离层foF2预测", 电波科学学报, vol. 33, no. 06, 15 December 2018 (2018-12-15) * |
Also Published As
Publication number | Publication date |
---|---|
CN115078848B (zh) | 2024-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109946762B (zh) | 一种基于概率分布短时预报降水的方法及系统 | |
Kazumori | Satellite radiance assimilation in the JMA operational mesoscale 4DVAR system | |
CN109344812B (zh) | 一种改进的基于聚类的单光子点云数据去噪方法 | |
CN113281742A (zh) | 一种基于滑坡形变信息和气象数据的sar滑坡预警方法 | |
KR101994387B1 (ko) | 수증기 가중치 병합을 이용한 산악지역 강수실황 예측을 위한 장치 및 방법 | |
KR101291980B1 (ko) | 기상레이더 반사도 자료의 통합품질지수 생성 방법 | |
CN113251914A (zh) | InSAR技术与长短时记忆神经网络结合的地表形变预测方法 | |
CN114966233A (zh) | 基于深度神经网络的雷电预报系统及方法 | |
Otsuka et al. | GSMaP RIKEN Nowcast: Global precipitation nowcasting with data assimilation | |
CN111666896A (zh) | 一种基于线性融合模型的遥感影像时空融合方法 | |
CN114186413A (zh) | 一种基于地表形变和孕灾环境条件的滑坡易发性评价方法 | |
CN116029162B (zh) | 利用星载gnss-r数据的洪涝灾害淹没范围监测方法和系统 | |
CN115078848A (zh) | 基于闪电信号的电离层无源被动探测方法 | |
Bouilloud et al. | Radar rainfall estimation for the post-event analysis of a Slovenian flash-flood case: application of the Mountain Reference Technique at C-band frequency | |
Kida et al. | Improvement of rain/no-rain classification methods for microwave radiometer observations over the ocean using a 37 GHz emission signature | |
Islam et al. | Using S-band dual polarized radar for convective/stratiform rain indexing and the correspondence with AMSR-E GSFC profiling algorithm | |
CN116361621A (zh) | 基于风云四号卫星的大气波导监测诊断方法和系统 | |
CN110727719B (zh) | 一种基于动力松弛逼近的闪电定位资料同化方法 | |
Tonboe et al. | The EUMETSAT OSI SAF Sea Ice Concentration Algorithm Algorithm Theoretical Basis Document | |
Behrangi et al. | Summertime evaluation of REFAME over the Unites States for near real-time high resolution precipitation estimation | |
KR102670033B1 (ko) | 다중 우적계 관측장비를 이용한 복잡지형 강우장 생성을 위한 시스템 및 방법 | |
Combarnous et al. | An observation operator for geostationary lightning imager data assimilation in the French storm-scale numerical weather prediction system AROME | |
Overeem et al. | Retrieval algorithm for rainfall mapping from microwave links in a cellular communication network. | |
CN114019255B (zh) | 一种雷暴中心的定位和预测方法及系统 | |
Ishii et al. | Study on measurement performance of future space-based Doppler wind lidar in Japan |
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 |