CN105184039A - 一种电离层垂直剖面建模及参数反演方法 - Google Patents
一种电离层垂直剖面建模及参数反演方法 Download PDFInfo
- Publication number
- CN105184039A CN105184039A CN201510335045.9A CN201510335045A CN105184039A CN 105184039 A CN105184039 A CN 105184039A CN 201510335045 A CN201510335045 A CN 201510335045A CN 105184039 A CN105184039 A CN 105184039A
- Authority
- CN
- China
- Prior art keywords
- layer
- paddy
- parameter
- represent
- delta
- 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
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种电离层垂直剖面建模及参数反演方法,所述方法包括以下步骤,步骤1:将电离层建模为包含E层、谷层、层、层的四层模型;步骤2:基于建立的电离层模型,推导各层回波虚高的计算公式;步骤3:利用实测电离图数据,结合各层回波虚高计算结果,进行各层电离层参数的反演。本发明所公开的电离层垂直剖面建模及参数反演方法,基于移位切比雪夫多项式模型的约束优化层参数、层参数的垂测电离图反演方法,即得到谷参数后,选取层较高区域回波描迹数据,在保证剖面连续光滑的约束条件下,计算层剖面多项式系数,该方法适用于层未充分发育的情况,可以有效提高该种情况下反演精度和稳定性。
Description
技术领域
本发明涉及电离层研究及应用领域,尤其涉及一种电离层垂直剖面建模及参数反演方法。
背景技术
电离层垂直探测(简称垂测)能够获得反映电离层虚高与频率关系的垂测电离图,所谓虚高,其实质是电磁波传播时间的度量,电磁波从垂测发射机发出,经电离层反射后到达垂测接收机所经历的时间乘以光速就称为虚高。垂测获得的虚高并不是电磁波在电离层中的真实反射高度,真实反射高度获得需要对垂测电离图进行反演,即利用垂测电离图频率-虚高描迹反演电离层剖面(电离层高度与等离子体频率或电子浓度的对应关系)。垂测电离图的反演对研究电离层结构和电离层波传播问题具有重要意义,一直以来受到十分广泛的重视,当然,反演也具有相当大的难度。
目前,应用较为普遍的垂测电离图反演方法是基于模式法思想发展的电离层参数反演方法,其中,黄雪钦等公开了一种基于移位切比雪夫多项式模型反演电离层剖面的方法,该方法中,将F层建模为移位切比雪夫多项式,求解满足计算虚高与实测虚高在最小二乘意义上最佳吻合的多项式系数,从而确定电离层剖面,后来,针对回波描迹中明显有F1层描迹的情况,做了模型上的一些改进,即将F1层和F2层分别用不同的移位切比雪夫多项式表示,但该方法更适用于F1层充分发育的情况,对于F1层未充分发育的情况,基于反演结果合成的垂测电离图与实测电离图在F1层临频附近差别较大,另外,未考虑层与层连接处的光滑性问题也是该方法的不足之处。
发明内容
本发明所要解决的技术问题就是提供一种电离层垂直剖面建模及参数反演方法。
本发明采用如下技术方案:
一种电离层垂直剖面建模及参数反演方法,其改进之处在于,所述方法包括以下步骤:
步骤1:将电离层建模为包含E层、谷层、F1层、F2层的四层模型,其中,E层和谷层剖面用抛物模型表示,F1层和F2层剖面用移位切比雪夫多项式模型表示;
步骤2:基于建立的电离层模型,推导各层回波虚高的计算公式;
步骤3:利用实测电离图数据,结合各层回波虚高计算结果,进行各层电离层参数的反演。
进一步的,所述步骤1具体包括:
步骤11:电离层电子浓度剖面具体形式如下式所示:
上式中各符号的具体含义如下:
E层:fNE表示E层等离子体频率;fCE表示E层临频;hmE表示E层峰高;ymE表示E层半厚;hbE=hmE-ymE表示E层底高;
谷层:fNV表示谷层等离子体频率;fCV表示谷层最小等离子体频率;hmV表示谷层等离子体频率为fCV时对应的电离层高度;ymV表示谷层半厚;h2=hmE+W,W定义为谷层宽度;
F1层:Ti(g)为移位切比雪夫多项式,具有下式所示形式:
fNF1表示F1层等离子体频率;fCF1表示F1层临频;Ai(i=0~I+1)为移位切比雪夫多项式系数,
且:
hmF1表示F1层峰高,且:
ΔfC为垂测电离图智能判读软件自动给出的fCF1相对于F1层模型设定临频的偏差;
F2层:移位切比雪夫多项式Ti(l)中的l具有下式所示形式:
fNF2表示F2层等离子体频率;fCF2表示F2层临频;Ci(i=0~N+1)为移位切比雪夫多项式系数,且:
hmF2表示F2层峰高,且:hmF2=CN+1,
E层和谷层的连接点位于E层峰高hmE,谷层与F1层的连接点位于高度h2处,并且在高度h2处的等离子体频率等于E层临频fCE,谷层包括两个部分:与E层的连接部分和与F1层的连接部分,这两部分的连接点位于高度h1处,F1层与F2层连接点位于F1层峰高hmF1处;
步骤12:为了使建立的电子浓度剖面满足连续光滑特性,在层与层的连接点处,基于连接点以上及以下电离层模型分别计算的等离子体频率的平方值以及剖面梯度应该相等,根据这一条件,限定相关参数之间的内在关系。
进一步的,所述步骤2具体包括:
步骤21:对于频率小于等于fCE的电波将在E层反射,进行E层回波虚高计算公式的推导;
步骤22:对于频率大于fCE、小于等于fCF1的电波将在F1层反射,进行F1层回波虚高计算公式的推导,主要包括在E层传播的群距离ΔhE(f)、在谷层中与E层连接部分传播的群距离Δh′J(f)、在谷层中与F1层连接部分传播的群距离Δh′V(f)和在F1层传播的群距离Δh′F1(f)的计算公式;
步骤23、对于频率大于fCF1、小于等于fCF2的电波将在F2层反射,进行F2层回波虚高计算公式的推导,主要包括Δh′E(f)、Δh′J(f)、Δh′V(f)、Δh′F1(f)和在F2层传播的群距离Δh′F2(f)的计算公式。
进一步的,所述步骤3具体包括:
步骤31:选取实测E层描迹数据,采用区域搜索的方法实现E层fCE、hmE(或hbE)、ymE三个参数的反演;
步骤32:选取F1层描迹中大于E层临频和F1层描迹最小虚高对应的频率之间的数据,基于搜索、迭代方法实现谷层参数的反演;
步骤33:选取F1层描迹中,F1层描迹最小虚高对应的频率到fCF1之间的数据,利用约束寻优法实现F1层参数的反演;
步骤34:选取F2层描迹数据,利用约束寻优法实现F2层参数的反演;
步骤35、对根据不同ΔfC得到的谷层、F1层和F2层参数,寻找使基于这些参数计算的所有数据点的计算虚高和实测虚高的误差和最小的一组参数,该组参数确定为最终的谷层、F1层和F2层参数。
本发明的有益效果在于:
本发明所公开的电离层垂直剖面建模及参数反演方法,基于移位切比雪夫多项式模型的约束优化F1层参数、F2层参数的垂测电离图反演方法,即得到谷参数后,选取F1层较高区域回波描迹数据,在保证剖面连续光滑的约束条件下,计算F1层剖面多项式系数,同样,在保证剖面连续光滑的约束条件下,选取F2层回波描迹数据,计算F2层剖面多项式系数,最后,基于所有数据点计算虚高和实测虚高误差和最小准则,选取对应初始设置下得到的剖面参数最终确定电离层剖面。该方法适用于F1层未充分发育的情况,可以有效提高该种情况下反演精度和稳定性。
附图说明
图1为本发明所公开的电离层垂直剖面建模及参数反演方法的流程图;
图2为三层电离层反演实例。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例1,如图1所示,本实施例公开了一种电离层垂直剖面建模及参数反演方法,包括以下步骤:
(1)建立电离层剖面数学模型:
本发明基于模式法的思想,将电离层建模为包含E层、谷层、F1层、F2层的四层模型,E层和谷层剖面用抛物模型表示,F1层和F2层剖面用移位切比雪夫多项式模型表示,电离层电子浓度剖面具有式(1)所示形式:
E层和谷层的连接点位于E层峰高hmE,谷层与F1层的连接点位于高度h2处,并且在高度h2处的等离子体频率等于E层临频fCE,谷层包括两个部分:与E层的连接部分和与F1层的连接部分,这两部分的连接点位于高度h1处,F1层与F2层连接点位于F1层峰高hmF1处,式(1)中各符号的具体含义如下:
E层:
fNE表示E层等离子体频率;fCE表示E层临频;hmE表示E层峰高;ymE表示E层半厚;hbE=hmE-ymE表示E层底高;
谷层:
fNV表示谷层等离子体频率;fCV表示谷层最小等离子体频率;hmV表示谷层等离子体频率为fCV时对应的电离层高度;ymV表示谷层半厚;h2=hmE+W,W定义为谷层宽度;
F1层:
Ti(g)为移位切比雪夫多项式,具有式(2)所示形式:
fNF1表示F1层等离子体频率;fCF1表示F1层临频;Ai(i=0~I+1)为移位切比雪夫多项式系数,且:
hmF1表示F1层峰高,且:
ΔfC为垂测电离图智能判读软件自动给出的fCF1相对于F1层模型设定临频的偏差;
F2层:
移位切比雪夫多项式Ti(l)中的l具有式(6)所示形式:
fNF2表示F2层等离子体频率;fCF2表示F2层临频;Ci(i=0~N+1)为移位切比雪夫多项式系数,且:
hmF2表示F2层峰高,且:
hmF2=CN+1(8)
为了使建立的电子浓度剖面满足连续光滑特性,在层与层的连接点处,基于连接点以上及以下电离层模型分别计算的等离子体频率的平方值以及剖面梯度应该相等,根据这一条件,限定相关参数之间的内在关系,即:
1)在h=hmF1处,有:
利用式(9)中的上式,可以由Ai(i=0~I+1)计算得到CN+1,如式(5)和式(7)所示。令则由式(9)中的下式可得:
式(10)为后续反演F2层系数Ci(i=0~N)的一个约束条件。
2)在h=h2处,有:
考虑到 并且令 以及 即:
式(12)也是后续反演F1层系数Ai(i=0~I)的一个约束条件,则由式(11)得到:
3)在h=h1处,有:
由式(14)中的上式可以得到:
hmV=hmE+[(4B2D+Q)D]1/2(15)
其中,式(15)结合式(13)中的上式以及W可以进一步得到:
由式(14)中的下式可以给出h1的计算结果为:
由以上推导可以看出,得到E层参数后,谷层只要确定B和W两个参数,则谷层剖面就可以确定,因此,在后续谷层参数反演中,我们只需确定B和W即可。
(2)各层反射回波虚高的计算:
通过对各层反射回波虚高的计算,得到计算虚高与实测虚高的误差量,从而可用于实现后续基于最小化误差量的各层参数反演。为了简化计算过程,同时又不会引入很大的误差,在计算电波在E层和谷层传播的群距离时没有考虑地磁场的影响,在计算电波在F1层和F2层传播的群距离时假设地磁场为一定值,即任意位置处的地磁场与垂测站上空300km高度处的地磁场一致。
1)E层回波虚高的计算
对于频率小于等于fCE的电波将在E层反射,回波虚高计算公式为:
其中,f为电波频率,hr为电波反射点处高度,μ′为群折射指数,不考虑地磁场时,具有如下形式:
式中fN表示对应位置处等离子体频率。
基于建立的E层电离层模型,则式(18)可以进一步计算得到:
2)F1层回波虚高的计算
对于频率大于fCE、小于等于fCF1的电波将在F1层反射,回波虚高计算公式为:
其中,式(21)中第二项为电波在E层传播的群距离,记为Δh′E(f),第三项为电波在谷层中与E层连接部分传播的群距离,记为Δh′J(f),第四项为电波在谷层中与F1层连接部分传播的群距离,记为Δh′V(f),第五项为电波在F1层传播的群距离,记为Δh′F1(f)。
计算Δh′E(f)、Δh′J(f)、Δh′V(f)时用到的μ′仍然具有式(19)的形式,计算结果分别为:
计算Δh′F1(f)时用到的μ′具有以下形式:
其中,
Yo=fH/f(30)
式中,fH为垂测站上空300km处磁旋频率,θ为垂测站上空300km处磁倾角。此时,无法直接给出Δh′F1(f)的解析表达式,只能采用数值积分方法计算Δh′F1(f),考虑到在反射点处,μ′为无穷大,为了能够进行数值计算,作以下变量替换,即:
则Δh′F1(f)可以写为:
根据建立的F1层电离层模型(式(1)中的第三个式子),可得到:
将式(35)代入到式(34),则Δh′F1(f)可进一步写为:
式中,
此时,在临近反射点(t→0)处,由式(26)~式(32),可以得到μo→0,M→1,从而由式(25)得到所以,只要磁倾角Si(f)就可以由式(37)计算得到(其中fCF1可由垂测电离图智能判读软件自动给出),从而可根据式(36)得到Δh′F1(f)。
3)F2层回波虚高的计算
对于频率大于fCF1、小于等于fCF2的电波将在F2层反射,回波虚高计算公式为:
其中,式(38)中第五项为电波在F1层传播的群距离,记为Δh′F1(f),第六项为电波在F2层传播的群距离,记为Δh′F2(f)。
此时,Δh′E(f)、Δh′J(f)、Δh′V(f)仍然可分别用式(22)、式(23)、式(24)计算,Δh′F1(f)仍然可以采用式(36)计算,但由于此时电波在F2层反射,其中的Si(f)的计算公式变为:
其中,
Δh′F2(f)采用与F1层回波虚高中Δh′F1(f)相同的计算方法得到:
式中,
此时,在临近反射点(s→0)处,由式(25)~式(32),得到fCF2可由垂测电离图智能判读软件自动给出。
(3)各层参数的反演:
基于计算得到的虚高和实测虚高的差最小化准则,实现各层参数的反演。
1)E层参数的反演
由式(1)中上式可知,决定E层剖面的三个参数主要是fCE、hmE(或hbE)、ymE,其中fCE可由垂测电离图智能判读软件自动给出,误差小于0.2MHz,本发明中采用一种区域搜索的方法实现E层参数的反演,具体为:
假设垂测电离图智能判读得到的E层描迹有K个点,其对应的工作频率和虚高分别为fk和h″(fk),读出的E层临频和最小虚高分别记为为和h″minE,则对参数fCE、hbE、ymE分别在[h″minE-δ1,h″minE+δ2]、[0,δ3](其中δ1、δ2和δ3是搜索范围控制量)以一定步进取值得到不同组参数,每一组参数根据式(20)计算得到K个点的h′(fk),然后计算实测虚高和计算虚高的误差平方和:
使ε达到最小的那组参数即确定为E层参数。
2)谷层参数的反演
本发明中基于搜索、迭代方法实现谷层B和W两个参数的反演。
在F1层描迹中,大于E层临频和F1层描迹最小虚高(记为h″minF1)对应的频率之间的数据对谷层参数比较敏感,因此,在谷层参数的反演过程中,选择这部分描迹点用于确定谷层参数,假设共有K个点,其对应的工作频率和虚高分别为fk和h″(fk)。
谷层参数反演的基本步骤为:
①设定W=0;
②设定B=0,I=7;
③基于最小二乘法计算F1层剖面系数Ai(i=0~I+1);
④检查计算的系数Ai(i=0~I)是否满足F1层剖面单调递增特性,(a)如果不满足,则I=I-1,如果I<0,则执行⑤,否则,执行③;(b)如果满足,对AI+1和前一次迭代记录的值进行比较,如果二者之差小于某一较小值(例如0.5km),则计算K个点实测虚高和计算虚高的误差平方和,记录当前B、W以及计算的误差平方和的值,否则,在没有超过限定的最大迭代次数情况下,根据Ai(i=0~I)自动调整I,按照式(12)更新B的值,执行③,若超过限定的最大迭代次数,则执行⑤;
⑤W=W+1(单位为km),如果W小于设定的搜索范围(例如0.7h″minF1-hmE),则执行②,否则,执行⑥;
⑥找出记录的误差平方和的最小值,该最小值对应的B、W即确定为谷层参数,如果没有记录下有效的B和W,则W=0,B=0。
其中,上述步骤③的具体方法为:
根据式(21)给出K个点实测虚高h″(fk)和计算虚高h′(fk)的误差平方和:
其中, 要使ε达到最小,即求解满足式(46)的系数Ai(i=0~I):
式(46)进一步化简为:
求解上述方程组即可得到系数Ai(i=0~I),然后按照式(4)可计算AI+1。
3)F1层参数的反演
选取F1层描迹中,F1层描迹最小虚高对应的频率到fCF1之间的数据用于确定F1层参数,假设共有K个数据点,其对应的工作频率和虚高分别为fk和h″(fk)。在读取了fCF1的情况下,F1层剖面由系数Ai(i=0~I+1)完全确定,可采用类似反演谷层参数时用到的方法来计算这些系数,这里需要注意的是,在谷层参数确定后,谷层与F1层交点处的剖面梯度也已经确定,那么,根据当前数据计算出的系数Ai(i=0~I)必须满足式(12),因此,F1层参数的反演实际上是一个约束优化问题,即:
可采用拉格朗日法求解上述问题进行F1层参数的反演,具体的为:
①根据式(48)建立一个新的函数:
②对式(49)中各个自变量求偏导数,建立方程组:
③求解上述方程组,得到系数Ai(i=0~I),检查计算的系数Ai(i=0~I)是否满足F1层剖面单调递增特性,(a)如果不满足,则I=I-1,如果I<0,执行④,否则,执行①;(b)如果满足,执行④;
④按照式(4)计算AI+1。
4)F2层参数的反演
F2层描迹中的数据用于确定F2层参数,假设共有K个数据点,其对应的工作频率和虚高分别为fk和h″(fk)。在读取了fCF2的情况下,F2层剖面由系数Ci(i=0~N+1)完全确定,在F1层参数确定后,F1层与F2层交点处的剖面梯度也已经确定,那么,根据当前数据计算出的系数Ci(i=0~N)必须满足式(10),因此,F2层参数的反演同样是一个约束优化问题,可采用类似反演F1层参数时用到的方法来计算这些系数:
式中, 同样,采用拉格朗日法求解上述问题进行F2层参数的反演,具体的为:
①根据式(48)建立一个新的函数:
②对式(49)中各个自变量求偏导数,建立方程组:
③求解上述方程组,得到系数Ci(i=0~N),检查计算的系数Ci(i=0~N)是否满足F2层剖面单调递增特性,(a)如果不满足,则N=N-1,如果N<0,执行④,否则,执行①;(b)如果满足,执行④;
④按照式(7)计算CN+1。
5)谷层、F1层、F2层参数的最终确定
对于未充分发育F1层,垂测电离图智能判读软件自动给出的fCF1相对于F1层模型设定临频的偏差ΔfC是一个未知量,理论上,ΔfC的取值介于0~fCF2-fCF1,因此,在反演谷层、F1层、F2层参数时,将ΔfC在0~fCF2-fCF1内遍历,选取使所有数据点计算虚高与实测虚高误差和最小的ΔfC对应的谷层、F1层、F2层参数作为最终的谷层、F1层、F2层参数。
按照上述方法对三层电离层反演后的结果如图2所示。
Claims (4)
1.一种电离层垂直剖面建模及参数反演方法,其特征在于,所述方法包括以下步骤:
步骤1:将电离层建模为包含E层、谷层、F1层、F2层的四层模型,其中,E层和谷层剖面用抛物模型表示,F1层和F2层剖面用移位切比雪夫多项式模型表示;
步骤2:基于建立的电离层模型,推导各层回波虚高的计算公式;
步骤3:利用实测电离图数据,结合各层回波虚高计算结果,进行各层电离层参数的反演。
2.根据权利要求1所述的电离层垂直剖面建模及参数反演方法,其特征在于,所述步骤1具体包括:
步骤11:电离层电子浓度剖面具体形式如下式所示:
上式中各符号的具体含义如下:
E层:fNE表示E层等离子体频率;fCE表示E层临频;hmE表示E层峰高;ymE表示E层半厚;hbE=hmE-ymE表示E层底高;
谷层:fNV表示谷层等离子体频率;fCV表示谷层最小等离子体频率;hmV表示谷层等离子体频率为fCV时对应的电离层高度;ymV表示谷层半厚;h2=hmE+W,W定义为谷层宽度;
F1层:Ti(g)为移位切比雪夫多项式,具有下式所示形式:
fNE1表示F1层等离子体频率;fCF1表示F1层临频;Ai(i=0~I+1)为移位切比雪夫多项式系数,且:
hmF1表示F1层峰高,且: ΔfC为垂测电离图智能判读软件自动给出的fCF1相对于F1层模型设定临频的偏差;
F2层:移位切比雪夫多项式Ti(l)中的l具有下式所示形式:
fNF2表示F2层等离子体频率;fCF2表示F2层临频;Ci(i=0~N+1)为移位切比雪夫多项式系数,且:
hmF2表示F2层峰高,且:hmF2=CN+1,
E层和谷层的连接点位于E层峰高hmE,谷层与F1层的连接点位于高度h2处,并且在高度h2处的等离子体频率等于E层临频fCE,谷层包括两个部分:与E层的连接部分和与F1层的连接部分,这两部分的连接点位于高度h1处,F1层与F2层连接点位于F1层峰高hmF1处;
步骤12:为了使建立的电子浓度剖面满足连续光滑特性,在层与层的连接点处,基于连接点以上及以下电离层模型分别计算的等离子体频率的平方值以及剖面梯度应该相等,根据这一条件,限定相关参数之间的内在关系。
3.根据权利要求1所述的电离层垂直剖面建模及参数反演方法,其特征在于,所述步骤2具体包括:
步骤21:对于频率小于等于fCE的电波将在E层反射,进行E层回波虚高计算公式的推导;
步骤22:对于频率大于fCE、小于等于fCF1的电波将在F1层反射,进行F1层回波虚高计算公式的推导,主要包括在E层传播的群距离Δh′E(f)、在谷层中与E层连接部分传播的群距离Δh′J(f)、在谷层中与F1层连接部分传播的群距离Δh′V(f)和在F1层传播的群距离Δh′F1(f)的计算公式;
步骤23、对于频率大于fCF1、小于等于fCF2的电波将在F2层反射,进行F2层回波虚高计算公式的推导,主要包括Δh′E(f)、Δh′J(f)、Δh′V(f)、Δh′F1(f)和在F2层传播的群距离Δh′F2(f)的计算公式。
4.根据权利要求1所述的电离层垂直剖面建模及参数反演方法,其特征在于,所述步骤3具体包括:
步骤31:选取实测E层描迹数据,采用区域搜索的方法实现E层fCE、hmE(或hbE)、ymE三个参数的反演;
步骤32:选取F1层描迹中大于E层临频和F1层描迹最小虚高对应的频率之间的数据,基于搜索、迭代方法实现谷层参数的反演;
步骤33:选取F1层描迹中,F1层描迹最小虚高对应的频率到fCF1之间的数据,利用约束寻优法实现F1层参数的反演;
步骤34:选取F2层描迹数据,利用约束寻优法实现F2层参数的反演;
步骤35、对根据不同ΔfC得到的谷层、F1层和F2层参数,寻找使基于这些参数计算的所有数据点的计算虚高和实测虚高的误差和最小的一组参数,该组参数确定为最终的谷层、F1层和F2层参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510335045.9A CN105184039B (zh) | 2015-06-17 | 2015-06-17 | 一种电离层垂直剖面建模及参数反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510335045.9A CN105184039B (zh) | 2015-06-17 | 2015-06-17 | 一种电离层垂直剖面建模及参数反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105184039A true CN105184039A (zh) | 2015-12-23 |
CN105184039B CN105184039B (zh) | 2018-07-27 |
Family
ID=54906117
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510335045.9A Expired - Fee Related CN105184039B (zh) | 2015-06-17 | 2015-06-17 | 一种电离层垂直剖面建模及参数反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105184039B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105549007A (zh) * | 2016-01-05 | 2016-05-04 | 中国电子科技集团公司第二十二研究所 | 一种基于交叠多项式模型的垂测电离图反演方法 |
CN105701276A (zh) * | 2016-01-05 | 2016-06-22 | 中国电子科技集团公司第二十二研究所 | 一种改进的电离层垂直剖面建模方法 |
CN109614585A (zh) * | 2018-12-03 | 2019-04-12 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种新的电离层区域重构方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102411664A (zh) * | 2010-09-21 | 2012-04-11 | 中国电子科技集团公司第二十二研究所 | 基于返回散射和斜测电离图联合反演电离层参数的方法 |
CN106134480B (zh) * | 2009-09-18 | 2013-04-17 | 中国电子科技集团公司第二十二研究所 | F1层不完全发育电离层模型及反演方法 |
CN106134481B (zh) * | 2009-09-18 | 2014-02-19 | 中国电子科技集团公司第二十二研究所 | 基于返回散射电离图反演Es层参数的方法 |
-
2015
- 2015-06-17 CN CN201510335045.9A patent/CN105184039B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106134480B (zh) * | 2009-09-18 | 2013-04-17 | 中国电子科技集团公司第二十二研究所 | F1层不完全发育电离层模型及反演方法 |
CN106134481B (zh) * | 2009-09-18 | 2014-02-19 | 中国电子科技集团公司第二十二研究所 | 基于返回散射电离图反演Es层参数的方法 |
CN102411664A (zh) * | 2010-09-21 | 2012-04-11 | 中国电子科技集团公司第二十二研究所 | 基于返回散射和斜测电离图联合反演电离层参数的方法 |
Non-Patent Citations (7)
Title |
---|
BODO W. REINISCH 等: "REAL TIME ELECTRON DENSITY PROFILES FROM IONOGRAMS", 《SPACE RESEARCH》 * |
BODO W.REINISCH 等: "Automatic calculation of electron density profiles from digital ionograms 3. Processing of bottomside ionograms", 《RADIO SCIENCE》 * |
HUANG XUEQIN 等: "Automatic calculation of electron density profiles from digital ionograms2.True height inversion of topside ionograms with the profile-fitting method", 《RADIO SCIENCE》 * |
J. E. TITHERIDGE: "The real height analysis of ionograms:A generalized formulation", 《RADIO SCIENCE》 * |
P. L. DYSON: "A model of the vertical distribution of the electron concentration in the ionosphere and its application to oblique propagation studies", 《JOURNAL OF ATMOSPHERIC AND TERRESTRIAL PHYSICS》 * |
王世凯 等: "垂测电离图F1层反演方法研究", 《电波科学学报》 * |
郑传青: "垂直探测频高图的反演", 《电波科学学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 | 中国电子科技集团公司第二十二研究所 | 一种改进的电离层垂直剖面建模方法 |
CN109614585A (zh) * | 2018-12-03 | 2019-04-12 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种新的电离层区域重构方法 |
CN109614585B (zh) * | 2018-12-03 | 2023-01-24 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种新的电离层区域重构方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105184039B (zh) | 2018-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102879820B (zh) | 基于三角网格的三维表层模型构建方法 | |
CN105676277B (zh) | 一种提高高陡构造速度反演效率的全波形联合反演方法 | |
US9753181B2 (en) | System and method for automatic local grid refinement in reservoir simulation systems | |
CN102759746A (zh) | 一种变偏移距垂直地震剖面数据反演各向异性参数方法 | |
CN108664705B (zh) | 一种基于OpenFOAM的模拟复杂地形地表粗糙度的方法 | |
CN106814391A (zh) | 基于菲涅尔体层析反演的地面微地震事件定位方法 | |
Karssenberg et al. | Conditioning a process-based model of sedimentary architecture to well data | |
CN105184039A (zh) | 一种电离层垂直剖面建模及参数反演方法 | |
CN105005682A (zh) | 一种垂测电离图反演方法 | |
CN105956271A (zh) | 一种隧道盾构施工引起地层位移的计算方法 | |
CN107817516A (zh) | 基于初至波信息的近地表建模方法及系统 | |
CN109444955A (zh) | 三维地震射线追踪的双线性走时扰动插值方法 | |
CN110515122B (zh) | 正演网格搜索定位及微地震信号识别方法及装置 | |
CN105093279A (zh) | 针对山前带的三维地震初至波菲涅尔体层析反演方法 | |
CN105160156A (zh) | 一种融合数据预处理的垂测电离图反演方法 | |
CN117331119A (zh) | 一种面向隧道探测的快速地震波走时计算方法 | |
CN107765306A (zh) | 一种vsp初始速度建模方法及装置 | |
CN111091235A (zh) | 变电站站区进出线路径的确定方法及装置 | |
CN105549007A (zh) | 一种基于交叠多项式模型的垂测电离图反演方法 | |
EP3186665B1 (en) | Method of characterising a subsurface volume | |
CN110909448B (zh) | 一种高频天波返回散射电离图反演方法 | |
CN107895076A (zh) | 一种基于目标的非贯穿型河道建模方法及系统 | |
KR20220100751A (ko) | 유입 풍속과 도시 형태 매개 변수를 이용한 도시 풍속 분석 시스템 및 이를 이용한 분석 방법 | |
CN111765870A (zh) | 基于倾斜摄影技术并分区域累加的土方量计算方法 | |
CN110602635A (zh) | 一种室内地图匹配增强定位方法、设备及存储设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180727 Termination date: 20210617 |