一种电离层水平不均匀结构重构方法技术

本发明专利技术公开一种电离层水平不均匀结构重构方法，包括如下步骤：（一）建立区域电离层电子浓度网格，在返回散射探测方位上，按照地面距离和高度分别划分网格；（二）准备反演算法的输入，（三）利用频率范围的实测前沿数据进行反演，采用求解非线性问题的Newton-Kontorovich方法进行迭代求解，（四）接着利用频率进行反演，（五）经过n次反演，最终得到二维电子浓度剖面，即我们的最优反演结果。本发明专利技术所公开的电离层水平不均匀结构重构方法，克服了现有技术中的缺点，基于解空间约束的思想，提出了频段递增逐步逼近反演算法，对解空间做出了合理限制。

【技术实现步骤摘要】

本专利技术涉及电离层研究及应用领域，尤其涉及一种电离层水平不均匀结构重构方 法。
技术介绍
高频返回散射探测作为重要的电离层探测手段，能够实现遥远区域电离层空间上 的连续监测，探测获取的高频返回散射扫频电离图显示了探测频率-群路径-回波能量三者 之间的关系。电离图包含了探测路径上的电离层状态信息，通过对其反演可以实时获取大 面积范围的电离层参数。由于电离层的时间聚焦和球形聚焦等效应，电离图上有着较清晰、 陡峭的前沿(也称最小群时延），一般能够准确判读。返回散射前沿除了受电离层电子浓度 分布影响之外，几乎不受任何其他因素（如天线波束、地面特性等）的影响。因此，返回散射 前沿被广泛用来进行电离层的反演研究。 不少学者已经致力于利用返回散射前沿反演电离层参数的研究，目前反演方法大 致可以归纳为如下几类：（1)基于迭代的拟合算法。首先是Rao(1974)利用返回散射电离图 前沿任意三组(P人f)层数据反演准抛物模型的三个参数，随后，Rao(1975)提出了一种利用 返回散射电离图离散散射源回波描迹反演电离层参数和散射源地面距离的方法。这两种反 演算法都假设电离层球形对称，不考虑地磁场的影响，并且都是反演单层准抛物模型参数。 1979年DuBroff将这种反演算法推广到电离层非球形对称的情况，假设了一个简单的梯度 电离层模型，反演包括准抛物模型的3个参数及其梯度共6个参数。（2)引入求解不适定问题 的理论与方法以解决返回散射电离图反演不稳定性问题。Chuang等（1977)利用地球物理反 演中发展起来的BG理论建立了斜向探测和返回散射探测反演模型，并用简单模型的垂测反 演结果验证了该方法的有效性。Fridman等（1994)以返回散射发射站处垂直探测得到的电 离层电子浓度剖面和返回散射前沿作为输入，利用求解非线性问题的Newton-Kontorovich 方法和线性不适定问题的Tikhonov正则化方法，反演电离层电子浓度的二维分布，获得电 离层的水平不均匀结构。遗传算法是近些年来发展起来的一种高效的非线性全局优化算 法，也有学者利用遗传算法来进行返回散射电离图的反演(谢树果，2005)。宋君等(2011)使 用模拟退火方法对返回散射电离图前沿进行了反演，但是该方法只适用于电离层较平稳时 期。现有的反演方法一般采用全频段反演的方式，在电离层存在较复杂水平不均匀特性时 不能很好地反映电离层电子浓度的变化趋势。
技术实现思路
本专利技术所要解决的技术问题就是提供一种可以如实反映电离层水平方向的不均 匀结构的电离层水平不均匀结构重构方法。 本专利技术采用如下技术方案： -种电离层水平不均匀结构重构方法，包括如下步骤： ( -）建立区域电离层电子浓度网格，在返回散射探测方位上，按照地面距离和高 度分别划分网格； (二)准备反演算法的输入，包括： ①实测的返回散射前沿以测(〇; ②利用发射站的垂测电离图反演得到返回散射发射站上空的电子浓度剖面No (h)； ③将实测前沿G^i(f)对应的探测频率范围划分为η个频段，每一频段的长 度记为八乜（1 = 1，2，~，]1)，获取频率范围，，~，=对应的实测前沿的测量误差先验值6142，~4 11，一般有61<62 < ··· < δη；(三)利用频率范围的实测前沿数据(}_()进行反演，采 用求解非线性问题的Newton-Kontorovich方法进行迭代求解，具体步骤如下： ①初始二维电子浓度剖面认为是水平均匀的，即步骤(一）中的电子浓度网格每一 地面距离上空的电子浓度剖面都等同于返回散射发射站处的电子浓度剖面No(h)，则有N (h，x)=No(h)，在此电离层模型下，利用电离层短波数字三维射线追踪技术合成频率范围 ); ②使用符号| | | |表示均方根，定义理论前沿值与实测前沿值之间的均方根误差 为： 其中Μ表示参与反演的前沿点个数； 如果下列条件满足 iGs^tfa^fa+AfiD-G^Ktf^fa+Afi]) I I <δχ (9) 则认为步骤(三)的①中的初始电子浓度剖面就是反演问题的最优解，否则进行下 面的迭代过程； ③求解频率范围的前沿数据对应 反演问题，上式中K(f，x)称为核函数，f表示工作频率，x表示地面距离，:表示均匀电离 层模型下返回散射前沿对应的地面散射点到发射站的地面距离，Gl(f)=G(f)-G〇(f)，其中G (f)表示实测的返回散射前沿，G〇(f)表示均匀电离层模型下的返回散射前沿，u(x)是待求 函数;根据Tikhonov正则化方法，求解反演问题：等价于求的极小值问题，后者又等价于求解下列 由（1 〇)式积分微分方程和(11)式边界条件组成的问题： 其中 Gi(f)=G^i(f)_G理仑⑴ 这里，;^ = ;1^，1；'2 = ;1^+八；|；'1，0^1](；〇=6獅（)，6理仑(;〇=6理仑（)； 上述问题中，u是我们待求的未知函数；K是已知核函数，按照 描述的方法计算，其中，上式中，fp表示等离子体频率;S = sinP，P表示射线的出射角，即射线传播方向与 垂直方向的夹角，有Sm=sini3m，其中妒表示返回散射前沿对应的射线出射角;Re表示地球半 径;下标"0"表示均匀电离层模型下的参量，参数q-般需要满足条件q 2 L2,其中L表示预期 的电子浓度水平变化尺度典型值，考虑到中炜地区电离层水平不均匀性的典型尺度通常为 上千公里量级或者更大，因此这里取值q=l〇 6km2,由于我们引进了电子浓度网格，可以将上 述问题表达式离散化，采用数值方法进行求解，离散化形式如下：满足边界条件： u〇 = ui ,Uv+i = uv (13) 其中假设网格是均匀划分的，v代表地面距离网格个数，w代表探测频率个数，d表示地 面距离网格的宽度，S表示探测频率步进；固定正则化参数α的值，利用Cholesky分解法求解上述v个方程可以得到v个未知 数的解= 1,…J)，这是Ua(x)的离散化形式；1的α中的最大值称为最优正则化参数ad，它对应的解W ad (X)>称为正则化解，也就是我们 要求的解；表示实测前沿与利用射线追踪合成的 理论前沿的差值：，f#Pf2分别表示起始工作频率和终止 工作频率，δ称为前沿测量误差的上界，将代入到（1)式，便得到反演后的电子浓度 剖面 N(h，x); ④在新的电子浓度剖面N(h，x)下，运用射线追踪技术合成频率范围 的理论前沿G舰()，将其代入(9)式，如果满足条件，则1^，1)就是反演问题的 最优解，记为N^ha)，反演过程结束，否则，转入步骤(三)的③； (四）接着利用频率范围的实测前沿数据6_( )进行反演，其反演方法与步骤(三)相同，只不过需将步骤(三）的①当前第1页1 2 3 本文档来自技高网...

【技术保护点】
一种电离层水平不均匀结构重构方法，其特征在于，包括如下步骤：(一)建立区域电离层电子浓度网格，在返回散射探测方位上，按照地面距离和高度分别划分网格；(二)准备反演算法的输入，包括：①实测的返回散射前沿G实测(f)；②利用发射站的垂测电离图反演得到返回散射发射站上空的电子浓度剖面N0(h)；③将实测前沿G实测(f)对应的探测频率范围[fa,fb]划分为n个频段，每一频段的长度记为Δfi(i＝1,2,…,n)，获取频率范围[fa,fa+Δf1]，[fa,fa+Δf1+Δf2]，…，[fa,fa+Δf1+Δf2+…+Δfn]＝[fa,fb]对应的实测前沿的测量误差先验值δ1，δ2，…，δn，一般有δ1≤δ2≤…≤δn；(三)利用频率范围[fa,fa+Δf1]的实测前沿数据G实测([fa,fa+Δf1])进行反演，采用求解非线性问题的Newton‑Kontorovich方法进行迭代求解，具体步骤如下：①初始二维电子浓度剖面认为是水平均匀的，即步骤(一)中的电子浓度网格每一地面距离上空的电子浓度剖面都等同于返回散射发射站处的电子浓度剖面N0(h)，则有N(h,x)＝N0(h)，在此电离层模型下，利用电离层短波数字三维射线追踪技术合成频率范围[fa,fa+Δf1]的理论前沿G理论([fa,fa+Δf1])；②使用符号|| ||表示均方根，定义理论前沿值与实测前沿值之间的均方根误差为：其中M表示参与反演的前沿点个数；如果下列条件满足||G理论([fa,fa+Δf1])‑G实测([fa,fa+Δf1])||≤δ1   (9)则认为步骤(三)的①中的初始电子浓度剖面就是反演问题的最优解，否则进行下面的迭代过程；③求解频率范围[fa,fa+Δf1]的前沿数据对应的反演问题，上式中K(f,x)称为核函数，f表示工作频率，x表示地面距离，表示均匀电离层模型下返回散射前沿对应的地面散射点到发射站的地面距离，G1(f)＝G(f)‑G0(f)，其中G(f)表示实测的返回散射前沿，G0(f)表示均匀电离层模型下的返回散射前沿，u(x)是待求函数；根据Tikhonov正则化方法，求解反演问题等价于求解Φα[u,G~1]=∫f1f2[Au-G~1(f)]2df+αΩ[u]]]>的极小值问题，后者又等价于求解下列由(10)式积分微分方程和(11)式边界条件组成的问题：∫0x0mK‾(ξ,t)u(t)dt-α{ddξ[q×du(ξ)dξ]-u(ξ)}=b(ξ)---(10)]]>满足边界条件du(ξ)dξ|ξ=0=du(ξ)dξ|ξ=x0m=0---(11)]]>其中K‾(ξ,t)=∫f1f2K(f,ξ)K(f,t)df]]>b(ξ)=∫f1f2K(f,ξ)G1(f)df]]>G1(f)＝G实测(f)‑G理论(f)这里，f1＝fa，f2＝fa+Δf1，G实测(F)＝G实测([fa,fa+Δf1])，G理论(f)＝G理论([fa,fa+Δf1])；上述问题中，u是我们待求的未知函数；K是已知核函数，按照K(f,x)=(S0m)32ρ2-(S0m)2B(x,x0m)-ρ22S0m(ρ2-ϵ~0)ϵ~0(S0m)22∫xx0mB(x,x′)ρ2ϵ~0(1ϵ~0∂ϵ~0∂z-]]>2RE)dx′]]>描述的方法计算，其中，B(x,x0m)=z0s(x)z0x(x0m)-z0x(x)z0s(x0m)(S0m)2×(∂(ρ2-ϵ~0)∂z-ρ2-ϵ~0ϵ~0∂ϵ~0∂z)+ρ2-ϵ~0ϵ~0×[z0x(x0m)(z0xx(x)z0s(x)-z0x(x)z0sx(x))-2z0xx(x)z0x(x)z0s(x0m)]]]>z0x=∂z0∂x]]>z0xx=∂2z0∂x2]]>z0s=∂z0∂S]]>z0sx=∂2z0∂S∂x]]>ρ＝1+h/REz=∫0hdh′/ρ]]>ϵ...

【技术特征摘要】

【专利技术属性】
技术研发人员：冯静，柳文，蔚娜，娄鹏，鲁转侠，杨龙泉，郭文玲，
申请(专利权)人：中国电子科技集团公司第二十二研究所，
类型：发明
国别省市：山东;37