CN106134480B - F1层不完全发育电离层模型及反演方法 - Google Patents
F1层不完全发育电离层模型及反演方法Info
- Publication number
- CN106134480B CN106134480B CN200910122586.8A CN200910122586A CN106134480B CN 106134480 B CN106134480 B CN 106134480B CN 200910122586 A CN200910122586 A CN 200910122586A CN 106134480 B CN106134480 B CN 106134480B
- Authority
- CN
- China
- Prior art keywords
- layer
- inversion
- height
- virtual height
- frequency
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 239000005433 ionosphere Substances 0.000 claims abstract description 38
- 238000005259 measurement Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000011161 development Methods 0.000 claims description 21
- 238000001514 detection method Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000009191 jumping Effects 0.000 claims 1
- 238000013507 mapping Methods 0.000 abstract 1
- 230000014509 gene expression Effects 0.000 description 15
- 238000004422 calculation algorithm Methods 0.000 description 10
- 238000011160 research Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000010845 search algorithm Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供一种F1层不完全发育电离层模型及反演方法,包括下列步骤:步骤A:进行数据预处理,根据实测垂测图,确定当前电离层存在的层数,以及频率-虚高描迹信息,并根据描迹信息进行各层频点抽取;步骤B:进行分层处理,并根据各层描迹信息,首先对E层进行反演,在E层反演的基础上,反演F1层和F2层;步骤C:判断F1层是否充分发育,在反演计算时根据F2层实测描迹前端的频率-虚高数据,推测F1层的发育情况;根据F1层的发育情况按照相应的模型进行反演计算。本发明能够有效的反演被遮蔽的F1层参数,同样能够满足实测虚高与计算虚高均方根误差最小。
Description
技术领域
本发明涉及电离层垂直探测技术领域,特别涉及一种F1层不完全发育电离层模型及反演方法。
背景技术
在众多的电离层探测技术中,迄今为止,基于地面的垂直探测仍然是短波通信、导航、雷达电波修正以及空间物理研究等方面最重要、最经济的探测手段。数十年来,世界各地许多地方都布置了电离层垂直探测站,而垂直探测获得的垂测电离图的反演技术也受到十分广泛的重视。垂测电离图是描述电离层虚高对探测频率关系的曲线图,由这种电离图,通过某种反演技术可以获得电离层电子浓度随高度分布的剖面信息,由电离图反演电离层剖面,国内外已经做过许多研究工作,并形成了若干种反演算法,Kelso系数法和Budden的矩阵法不需要知道电离层剖面模型,直接由垂测电离图反演电离层剖面,但这些方法需要很高质量的电离图,描迹完整,分辨率高。而且当出现电离层谷时,处理十分困难,计算结果的误差很大。当前反演算法主要采用的是一种模式法,即事先假设一个电离层电子浓度剖面模型的解析形式,并与观测的电离图描迹拟合,从而反演这个模型的参数,文献[Huang Xueqing,Bodo W.Reinisch.Automatic calculation of electron density profiles from digital ionograms 2.Trueheight inversion of topside ionograms with the profile-fitting method[J].Radio Sci.Vol.17,No.4:837-844.Bodo W.Reinisch,Huang Xueqin,Automatic calculation ofelectron density profile from digital ionograms 3.Processing of bottomsideionograms[J].Radio Sci.1983,Vol.18,No.3:477-492.]的反演算法是一种比较新的算法,它基于质量很高的数字电离图,对于底部电离层,E层的反演采用一种搜索算法,由于E层参数变化范围不大,这种算法是很有效的,对于谷区采用反抛物线模型,F层采用切比雪夫多项式拟合。这种反演方法也要求电离图的质量比较高。
为了保证剖面的光滑性,Dyson建立了一种准抛物段模型电离层[P.L.Dyson,J.A.Bennett,A model of the vertical distribution of the electronconcentration in the ionosphere and its application to oblique propagation studies[J].J.atmos.terr.Phys,1988,Vol.50,No.3:251-262],在各层之间引入一个连接层,将参数的反演转化为一个有约束条件的最优化问题求解。
为了便于射线追踪和电离层重构的应用,需要能够解析表达电离层剖面,同时满足电子浓度连续和导数连续的要求,目前通用的垂测电离图反演算法是通过对频高图中每一点的虚高和频率进行逐点反演计算,最后给出的是该垂测图对应的电子浓度点迹分布图,然后再利用QPS电离层模型,拟合反演的点迹,这样增加了处理步骤,更是增加了拟合带来的误差。
然而,在实际工程中,比如天波超视距雷达系统,需要获得当前电离层的特征参数foE(E层临界频率),hmE(E层最大电离高度)、foF1(F1层临界频率)、hmF1(F1最大电离高度)、foF2(F2层临界频率)、hmF2(F2层最大电离高度),从而获得实时大区域电离层的状态,提高雷达目标的定位精度。但是在实际情况中,垂测图的变化并不是规则的,并且对垂测图的判读是有误差的,这就需要一种能够有效克服判读误差,电离图描迹不规则,甚至是某些层特别是F1层被Es遮蔽时的反演算法。
发明内容
本发明的目的在于,提供一种F1层不完全发育电离层模型及反演方法,能够对正规的典型的垂测图进行反演,特别能够对F1层不完全发育,F1层被Es层遮蔽的垂测图进行反演。
本发明的有益效果是:依照本发明的基于QPS模型的垂测电离图反演方法,直接利用准抛物段模型(QPS模型)来表征电离层电子浓度剖面,并通过调整模型参数,比较实测虚高与计算虚高,使实测虚高与计算虚高均方根误差最小时对应的那组参数即为所需要的电离层特征参数,该方法一定程度上克服了判读误差,在虚高最小均方根误差的意义上是最优的。对Es层遮蔽F1层,甚至部分F2层的情况,提出了利用F2层信息和E层的信息进行反推的算法,有效地反演了被遮蔽的F1层参数,同样能够满足实测虚高与计算虚高均方根误差最小。
附图说明
图1为垂测电离层反演流程图。
具体实施方式
以下,参考附图1详细描述本发明的F1层不完全发育电离层模型及反演方法。
本发明的F1层不完全发育电离层模型及反演方法,包括下列步骤:
步骤100:进行数据预处理操作。
对于频高图中变化比较缓慢的地方,抽取有限的数据参与反演,就能满足上面的要求,使运算时间降低60%,甚至更多。其中,各层频点抽取的原则如下:
如果各层的频点少于10个,则不抽取。否则,
E层保证频点数大于10,小于20个,且最后3个频点保留。
F1层保证频点数大于10,小于25个,且最后5个频点保留。
F2层保证频点数大于10,小于25个,且最后8个频点保留。
其中,保留最后几个频点的目的是因为它们决定了临界频率、最大电离高度的反演最终结果,如果保留的少,就会影响精度。
表1 数据预处理反演参数对比表
步骤200:各层(E、F1、F2)反演处理过程:
每一层(E、F1、F2)参数的反演最终确定层的三个参数临界频率fc、最大电离高度rm和半厚度ym。采用一种搜索的算法来实现每层电离图的反演。下面针对E层来说,假设在电离图E描迹上取K个频率点,与其相应的K个虚高观测值为h′k,计算的虚高为h″k,计算值与观测值的均方误差为:
搜索的结果是找到一组参数fcE、rmE和ymE,使得ε最小。
E层参数反演过程如下:
①利用国际参考电离层模型(IRI)预测fcE、rmE和ymE,确定三个参数的搜索范围。
②选择E描迹的K个频率点,确定与之对应的K个虚高观测值。
③设置fcE、rmE和ymE的初始值。
④根据(22)式计算均方根误差ε。
⑤fcE增加一个步长,如果fcE超出搜索范围,跳转步骤⑨。
⑥rmE增加一个步长,如果rmE超出搜索范围,返回至步骤⑤。
⑦ymE增加一个步长,如果ymE超出搜索范围,返回至步骤⑥。
⑧根据(22)式计算均方根误差ε,返回至步骤⑦。
⑨比较计算的均方根误差ε,得到最小的ε,从而确定最佳的一组参数fcE、rmE和ymE。
确定E层参数以后,根据同样的步骤确定F1层的参数和F2层的参数。
步骤300:判断F1层是否被Es全遮蔽,如果是,则F1层按充分发育处理,否则,判断F1层是否正常发育;
其中,可以按照如下方法判断F1层是否正常发育:
根据判别准则,即对F2层频率-虚高数据前端的频点和虚高进行多项式拟合,并求拟合后的斜率最小值,若斜率的最小值不大于一定的数值(该数值需要多幅数据图形统计确定),则认为F1层完全发育,按照完全发育模型进行反演计算;否则说明F1层可能未充分发育,此时在计算F1层的虚高时,应用F1层未充分发育的反演方法。
其中,F1层被Es全遮蔽时,根据F2层参数推F1层参数,并应用F1层完全发育算法进行。
另外,各层发育完全时,根据步骤200中的反演步骤可以得到电子浓度剖面和反演虚高。
对于F1层不完全发育时,依照未充分发育的QPS模型进行垂直电离层反演;对于F1层完全发育时,依照基准QPS模型进行垂直电离层反演。
具体地,基准QPS模型如下:
QP电离层模型电子浓度随高度的分布可描述为如下的函数形式:
其中rm=最大电子浓度对应的高度(km),rb=电离层底部高度(km),ym=rm-rb为电离层半厚度(km),fc=临界频率(MHz)。
(1)式可以改写为:
其中:
如果电离层由E层、F1层、F2层组成(QPS模型),则各层的电子浓度分布为(这里表示为等离子体频率的剖面分布):
E层
第一连接层,连接E层和F1层
F1层
第二连接层,连接F1和F2层
F2层 (3)
各段等离子密度的梯度为:
各层的参数分别为:E层由底高rbE,半厚度ymE,临界频率fcE确定,bE=aE(rbE/ymE)2。F1层由层的底高rbF1,半厚度ymF1,临界频率fcF1确定,bF1=aF1(rbF1/ymF1)2。F2层由层的底高rbF2,半厚度ymF2,临界频率fcF2确定,bF2=aF2(rbF2/ymF2)2。aj1、bj1、rj1、和aj2、bj2、rj2的表达式可根据各层连接处电子浓度和梯度相等的关系推出有:
aj1=aE,rj1=rmE,
aj2=aF1,rj2=rmF1,
其中:rc1为第一个连接层与F1层交点处的高度,rc2为第二个连接层与F2层交点处的高度。
对于F1层不完全发育时连接层的表达式:
若F1层未充分发育,则F1层与F2层之间的连接层起始高度应该从F1层测到的最后一个工作频点对应的电离层真高连起而不再是从F1最大电子浓度对应的真高处连起,因此该连接层参数对应的公式需要重新推导。电离层其它各分段层的参数表达式形式不变。推导过程如下:
对于连接层与F1层连接点处,依据电离层的连续性,F1层最后一个工作频点对应的等离子体频率应用连接层参数表达为:
式中,foF1由垂测自动判图给出,rp为F1层最后一个工作频点对应的电离层真高,对应表达式为:
由电离层的连续可导性可得:
令:
联立(5)、(7)、(8)解方程组可得
对于连接层与F2层的连接点处,依据电离层的连续可导性可以得到下面两个表达式:
令:
把表达式(13)、(14)分别带入(11)、(12)并联立方程组解得:
联立(10)、(16)解关于rj2的方程组可得
联立(9)、(15)可得下面的表达式:
解该方程得到
把表达式(17)带入表达式(19)消掉rj2可以得到rc2关于B、rp等参数的表达式:
确定rc2的表达式后,就可以得到rj2、aj2、bj2的具体表达式,即可以得到F1层未充分发育时F1层与F2层之间连接层参数的表达式:
至此,F1层不完全发育的电离层模型已经建立。
综上所述,依照本发明的F1层不完全发育电离层模型及反演方法,能够对正规的典型的垂测图进行反演,特别能够对F1层不完全发育,F1层被Es层遮蔽的垂测图进行反演,能够为研究及工程需求提供大量数据支持。
以上是为了使本领域普通技术人员理解本发明,而对本发明所进行的详细描述,但可以想到,在不脱离本发明的权利要求所涵盖的范围内还可以做出其它的变化和修改,这些变化和修改均在本发明的保护范围内。
Claims (5)
1.一种F1层不完全发育电离层反演方法,其特征在于,包括下列步骤:
步骤A:进行分层处理,依据垂直探测电离图探测到的频高图分为E-F1-F2三层模式,进行数据预处理,根据实测垂测图,确定当前电离层存在的层数,以及频率-虚高描迹信息,并根据描迹信息进行各层频点抽取;
步骤B:根据各层描迹信息,首先对E层进行反演,在E层反演的基础上,反演F1层和F2层;
步骤C:判断F1层是否充分发育,在反演计算时根据F2层实测描迹前端的频率-虚高数据,推测F1层的发育情况;根据F1层的发育情况按照相应的模型参数进行反演计算;其中,
在所述步骤A中,根据描迹信息进行各层频点抽取的原则为:如果各层的频点少于10个,则不抽取;否则进行抽取,并且,抽取后E层反演用的频点数大于10、小于20个,且最后3个频点保留;抽取后F1层反演用的频点数大于10、小于25个,且最后5个频点保留;抽取后F2层反演用的频点数大于10、小于25个,且最后8个频点保留;
在所述步骤C中,在判断F1层是否充分发育之前,进一步判断Es层是否完全遮蔽F1层,如果完全遮蔽,按F1层完全发育模型进行反演;
在所述步骤C中,判断F1层是否充分发育包括下列步骤:对F2层频率-虚高数据前端的频点和虚高进行多项式拟合,并求拟合后的斜率最小值,若斜率的最小值不大于设置的门限值,则认为F1层完全发育;否则说明F1层未充分发育。
2.如权利要求1所述的F1层不完全发育电离层反演方法,其特征在于,通过F1层不完全发育电离层模型来表征电离层电子浓度剖面,通过调整模型参数,使实测虚高与计算虚高均方根误差ε最小时对应的那组电离层的特征参数作为反演结果。
3.如权利要求2所述的F1层不完全发育电离层反演方法,其特征在于,电离层的特征参数包括临界频率fc、最大电离高度rm和半厚度ym。
4.如权利要求3所述的F1层不完全发育电离层反演方法,其特征在于,在所述步骤B中,根据各层描迹信息,首先对E层进行反演包括下列步骤:
步骤11:利用国际参考电离层模型IRI预测临界频率fcE、最大电离高度rmE和半厚度ymE,确定三个参数的搜索范围;
步骤12:选择E描迹的K个频率点,确定与之对应的K个虚高观测值;
步骤13:设置fcE、rmE和ymE的初始值;
步骤14:根据初始值计算虚高值,根据计算出的虚高值与所述虚高观测值的均方误差公式计算均方根误差ε;
步骤15:fcE增加一个步长,如果fcE超出搜索范围,跳转步骤19;否则执行步骤16;
步骤16:rmE增加一个步长,如果rmE超出搜索范围,返回至步骤15;否则执行步骤17;
步骤17:ymE增加一个步长,如果ymE超出搜索范围,返回至步骤16;否则执行步骤18;
步骤18:根据当前计算出的虚高值与所述虚高观测值的均方误差公式计算均方根误差ε,返回至步骤17;
步骤19:比较计算的均方根误差ε,得到最小的ε,从而确定最佳的一组参数fcE、rmE和ymE;
确定E层参数以后,根据同样的步骤确定F1层的参数和F2层的参数;
其中,步骤14和步骤18中的计算值与观测值的均方误差公式为:
其中,在电离图E描迹上取K个频率点,与其相应的K个虚高观测值为h′k,计算的虚高为h″k。
5.如权利要求4所述的F1层不完全发育电离层反演方法,其特征在于,在所述步骤B中,在E层反演的基础上,反演F1层和F2层时,如果电波在F1层反射,则虚高由E、E与F1层的连接层1和F1层3部分的虚高组成;如果电波在F2层反射,虚高由E、连接层1、F1层、F1与F2层的连接层2和F2层5部分的虚高组成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910122586.8A CN106134480B (zh) | 2009-09-18 | 2009-09-18 | F1层不完全发育电离层模型及反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910122586.8A CN106134480B (zh) | 2009-09-18 | 2009-09-18 | F1层不完全发育电离层模型及反演方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106134480B true CN106134480B (zh) | 2013-04-17 |
Family
ID=57251423
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910122586.8A Expired - Fee Related CN106134480B (zh) | 2009-09-18 | 2009-09-18 | F1层不完全发育电离层模型及反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106134480B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105160156A (zh) * | 2015-08-10 | 2015-12-16 | 中国电子科技集团公司第二十二研究所 | 一种融合数据预处理的垂测电离图反演方法 |
CN105184039A (zh) * | 2015-06-17 | 2015-12-23 | 中国电子科技集团公司第二十二研究所 | 一种电离层垂直剖面建模及参数反演方法 |
CN105549007A (zh) * | 2016-01-05 | 2016-05-04 | 中国电子科技集团公司第二十二研究所 | 一种基于交叠多项式模型的垂测电离图反演方法 |
CN105701276A (zh) * | 2016-01-05 | 2016-06-22 | 中国电子科技集团公司第二十二研究所 | 一种改进的电离层垂直剖面建模方法 |
-
2009
- 2009-09-18 CN CN200910122586.8A patent/CN106134480B/zh not_active Expired - Fee Related
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105184039A (zh) * | 2015-06-17 | 2015-12-23 | 中国电子科技集团公司第二十二研究所 | 一种电离层垂直剖面建模及参数反演方法 |
CN105184039B (zh) * | 2015-06-17 | 2018-07-27 | 中国电子科技集团公司第二十二研究所 | 一种电离层垂直剖面建模及参数反演方法 |
CN105160156A (zh) * | 2015-08-10 | 2015-12-16 | 中国电子科技集团公司第二十二研究所 | 一种融合数据预处理的垂测电离图反演方法 |
CN105160156B (zh) * | 2015-08-10 | 2017-12-26 | 中国电子科技集团公司第二十二研究所 | 一种融合数据预处理的垂测电离图反演方法 |
CN105549007A (zh) * | 2016-01-05 | 2016-05-04 | 中国电子科技集团公司第二十二研究所 | 一种基于交叠多项式模型的垂测电离图反演方法 |
CN105701276A (zh) * | 2016-01-05 | 2016-06-22 | 中国电子科技集团公司第二十二研究所 | 一种改进的电离层垂直剖面建模方法 |
CN105549007B (zh) * | 2016-01-05 | 2018-05-22 | 中国电子科技集团公司第二十二研究所 | 一种基于交叠多项式模型的垂测电离图反演方法 |
CN105701276B (zh) * | 2016-01-05 | 2018-09-25 | 中国电子科技集团公司第二十二研究所 | 一种改进的电离层垂直剖面建模方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103870845B (zh) | 点云聚类去噪过程中新的k值优化方法 | |
CN109143196B (zh) | 基于k分布海杂波幅度模型的三分位点参数估计方法 | |
CN110659693B (zh) | 基于k近邻分类的配电网快速拓扑识别方法、系统及介质 | |
CN106134480B (zh) | F1层不完全发育电离层模型及反演方法 | |
CN109448087B (zh) | 一种盾构隧道表面点云数据分段方法 | |
CN108280852B (zh) | 一种基于激光点云数据的门窗点云形状检测方法及系统 | |
CN106548484A (zh) | 基于二维凸包的产品模型散乱点云边界特征提取方法 | |
CN100573557C (zh) | 一种快速退火的基于mrf的sar图像分割方法 | |
CN110503721B (zh) | 基于加权径向基函数插值的断裂地形保持方法 | |
CN104899448B (zh) | 一种集合卡尔曼滤波静态局地化方案的自适应补偿方法 | |
CN102542556A (zh) | 超声图像乳腺肿瘤自动提取方法 | |
CN105320809A (zh) | 一种针对风电场空间相关性的风速预测方法 | |
CN113435630B (zh) | 一种产流模式自适应的流域水文预报方法及系统 | |
CN106296680A (zh) | 一种基于区域的多特征融合高分辨率遥感影像分割方法 | |
CN104155694B (zh) | 一种反射转换横波共检波点叠加剖面的剩余静校正方法 | |
CN103049653B (zh) | 基于em算法的g0分布参数最大似然估计方法 | |
CN105389821A (zh) | 一种基于云模型和图割相结合的医学图像分割方法 | |
CN109389553A (zh) | 基于t样条的气象传真图等值线插值方法 | |
CN102005049B (zh) | 基于单边广义高斯模型的sar图像变化检测阈值方法 | |
CN111104746A (zh) | 一种基于小波分析的河漫滩高程确定方法 | |
CN105426605A (zh) | 多屏三维航迹及噪声等值线实时显示方法 | |
CN105572741B (zh) | 一种计算3d高频静校正量的方法 | |
CN112129989B (zh) | 一种基于自适应聚类的电压暂降分段刻画方法 | |
CN112381029B (zh) | 一种基于欧氏距离的机载LiDAR数据建筑物提取方法 | |
CN115184998A (zh) | 一种基于改进U-net神经网络的瑞利波频散曲线自动提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
GR03 | Grant of secret patent right | ||
GRSP | Grant of secret patent right | ||
DC01 | Secret patent status has been lifted | ||
DCSP | Declassification of secret patent | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130417 Termination date: 20180918 |
|
CF01 | Termination of patent right due to non-payment of annual fee |