基于商用微波回传链路的近地面大气折射率廓线监测方法

摘要

本申请公开了基于商用微波回传链路的近地面大气折射率廓线监测方法，包括如下步骤：基于预设区域内商用微波的回传链路，获取每条所述回传链路的晴空实测衰减；设定并初始化近地面大气折射率廓线表征参数；根据所述商用微波的天线参数、链路接收端和发射端空间位置，基于所述表征参数和电磁传播计算模型，得到每条链路的晴空理论衰减；根据所述晴空实测衰减和所述晴空理论衰减，计算晴空衰减目标函数；基于所述晴空衰减目标函数，使用寻优算法，得到所述近地面大气折射率廓线表征参数得到最优解，完成近地面大气折射率廓线监测。本申请实现了区域范围内近地面大气折射率廓线的有效监测，具有成本低、时空分辨率高、覆盖范围广等优点。

权利要求

1.基于商用微波回传链路的近地面大气折射率廓线监测方法，其特征在于，包括如下步骤：基于预设区域内商用微波回传链路，获取每条所述商用微波回传链路的晴空实测衰减；

设定并初始化近地面大气折射率廓线的表征参数；

根据商用微波的天线参数、链路接收端和发射端空间位置，基于所述表征参数和电磁传播计算模型，得到每条商用微波回传链路的晴空理论衰减；

根据所述晴空实测衰减和所述晴空理论衰减，计算晴空衰减目标函数；

基于所述晴空衰减目标函数，使用寻优算法，得到所述近地面大气折射率廓线表征参数得到最优解，完成近地面大气折射率廓线监测。

2.根据权利要求1所述的基于商用微波回传链路的近地面大气折射率廓线监测方法，其特征在于，获取每条所述商用微波回传链路的所述晴空实测衰减的方法包括：选择所述预设区域内N条稳定运行的商用微波回传链路，并获取所述商用微波回传链路的发射端空间位置、接收端空间位置和天线参数；

按照预设时间间隔，分别计算每条商用微波回传链路的发射端平均功率和接收端平均功率；

根据所述发射端平均功率和所述接收端平均功率，得到每条所述商用微波回传链路的所述晴空实测衰减。

3.根据权利要求1所述的基于商用微波回传链路的近地面大气折射率廓线监测方法，其特征在于，所述电磁传播计算模型包括抛物方程法。

4.根据权利要求1所述的基于商用微波回传链路的近地面大气折射率廓线监测方法，其特征在于，根据商用微波的所述天线参数、所述链路接收端和发射端空间位置，基于所述表征参数和所述电磁传播计算模型，利用分步广角抛物方程模型，得到每条商用微波回传链路的所述晴空理论衰减。

5.根据权利要求1所述的基于商用微波回传链路的近地面大气折射率廓线监测方法，其特征在于，所述晴空衰减目标函数为：其中，M表示近地面大气折射率廓线表征参数，Regularizer(M)为正则化项。

6.根据权利要求5所述的基于商用微波回传链路的近地面大气折射率廓线监测方法，其特征在于，所述正则化项包括L1正则化项。

7.根据权利要求6所述的基于商用微波回传链路的近地面大气折射率廓线监测方法，其特征在于，所述近地面大气折射率廓线表征参数的最优解为使所述目标函数最小的解，表达式为其中，M表示近地面大气折射率廓线表征参数。

8.根据权利要求7所述的基于商用微波回传链路的近地面大气折射率廓线监测方法，其特征在于，所述寻优算法包括蚁群算法。

说明书

基于商用微波回传链路的近地面大气折射率廓线监测方法

技术领域

[0001] 本申请属于近地面大气环境参数获取技术领域，具体涉及基于商用微波回传链路的近地面大气折射率廓线监测方法。

背景技术

[0002] 近地面大气折射率廓线对电磁信号传播路径具有重要影响，异常的大气折射率垂直梯度往往会造成商用通信设备无法正常工作，同时也可能实现地基雷达的超视距探测。目前，近地面大气折射率廓线信息的监测手段较为有限。通过球载的无线电探空仪遥测大气的温、压、湿参数廓线，可以通过理论公式计算出相应的大气折射率廓线信息。然而，常规业务中无线电探空仪只每日定点放飞两次，导致无法获取精细的大气折射率廓线变化信息。此外，无线电探空仪在近地面的垂直分辨率也很低。折射计可以直接精确的测量大气折射率，但仪器笨重且价格昂贵，无法大范围推广。

[0003] 除了上述专业大气折射率探测方法，全球广泛存在的非合作电磁波源也可以用于大气折射率廓线的探测。最典型的例子是利用GPS信号穿过对流层所产生的大气折射量来计算大气折射率。该方法可以实现低成本、高垂直分辨率的全球大气折射率廓线监测，但是难以获取近地面范围的有效信息。除了GPS非合作电磁波源，近地面还广泛存在着大量商用微波回传链路，这些链路的运行频率在6～42GHz范围之间，链路长度在几千米～几十千米不等。商用微波回传链路电磁传播过程同样也会受到大气折射廓线影响进而造成接收端信号发生波动，因此这些信号之中也蕴含着大气折射率廓线信息，理论上可以用于大气折射率廓线监测，但如何利用商用微波的这种信号波动，却没有一套成熟可行的技术。

发明内容

[0004] 本申请提出了基于商用微波回传链路的近地面大气折射率廓线监测方法，根据区域内微波回传链路发射端和接收端功率，计算多链路晴空实测衰减；定义近地面大气折射率廓线表征参数并设置初始值，实现对大气折射率廓线的描述；基于链路接收端和发射端空间位置，结合电磁传播模型，计算多链路晴空理论衰减；根据实测衰减和理论衰减差异计算目标函数，通过一定寻优算法求解近地面大气折射率廓线表征参数最优解，实现区域范围内近地面大气折射率廓线有效反演。

[0005] 为实现上述目的，本申请提供了如下方案：

[0006] 基于商用微波回传链路的近地面大气折射率廓线监测方法，包括如下步骤：

[0007] 基于预设区域内商用微波回传链路，获取每条所述商用微波回传链路的晴空实测衰减；

[0008] 设定并初始化近地面大气折射率廓线的表征参数；

[0009] 根据商用微波的天线参数、链路接收端和发射端空间位置，基于所述表征参数和电磁传播计算模型，得到每条商用微波回传链路的晴空理论衰减；

[0010] 根据所述晴空实测衰减和所述晴空理论衰减，计算晴空衰减目标函数；

[0011] 基于所述晴空衰减目标函数，使用寻优算法，得到所述近地面大气折射率廓线表征参数得到最优解，完成近地面大气折射率廓线监测。

[0012] 可选的，获取每条所述商用微波回传链路的所述晴空实测衰减的方法包括：

[0013] 选择所述预设区域内N条稳定运行的商用微波回传链路，并获取所述商用微波回传链路的发射端空间位置、接收端空间位置和天线参数；

[0014] 按照预设时间间隔，分别计算每条商用微波回传链路的发射端平均功率和接收端平均功率；

[0015] 根据所述发射端平均功率和所述接收端平均功率，得到每条所述商用微波回传链路的所述晴空实测衰减。

[0016] 可选的，所述电磁传播计算模型包括但不限于抛物方程法。

[0017] 可选的，根据商用微波的所述天线参数、所述链路接收端和发射端空间位置，基于所述表征参数和所述电磁传播计算模型，利用分步广角抛物方程模型，得到每条商用微波回传链路的所述晴空理论衰减。

[0018] 可选的，所述晴空衰减目标函数为：

[0019]

[0020] 其中，Regularizer(M)为正则化项。

[0021] 可选的，所述正则化项包括但不限于L1正则化项。

[0022] 可选的，所述近地面大气折射率廓线表征参数的最优解为使所述目标函数最小的解，表达式为

[0023] 可选的，所述寻优算法包括但不限于蚁群算法。

[0024] 本申请的有益效果为：

[0025] 本申请公开了基于商用微波回传链路的近地面大气折射率廓线监测方法，在提取商用微波回传链路衰减信息的基础上，结合电磁传播计算模型，搜寻最优大气折射率廓线表征参数，实现了区域范围内近地面大气折射率廓线的有效监测；同时，本申请方法基于地面广泛现存的商用微波回传链路，无需投入额外仪器设备即可实现近地面大气折射率廓线信息获取，具有成本低、时空分辨率高、覆盖范围广等优点，应用价值极高。本申请方法可以作为近地面大气折射率廓线监测新方法应用到实际业务当中。

附图说明

[0026] 为了更清楚地说明本申请的技术方案，下面对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

[0027] 图1为本申请实施例的基于商用微波回传链路的近地面大气折射率廓线监测方法流程示意图；

[0028] 图2为本申请实施例的商用微波回传链路分布示意图；

[0029] 图3为本申请实施例的大气折射率廓线表征参数示意图；

[0030] 图4为本申请实施例的大气折射率廓线反演结果示意图。

具体实施方式

[0031] 下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。

[0032] 为使本申请的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本申请作进一步详细的说明。

[0033] 如图1所示，为本申请实施例的基于商用微波回传链路的近地面大气折射率廓线监测方法流程示意图，主要包括以下几个步骤：

[0034] 首先，基于预设区域内商用微波的回传链路，获取每条所述回传链路的晴空实测衰减。

[0035] 在本实施例中，该步骤的具体实现方法包括下面3个步骤：

[0036] 1.选择预设区域内4条稳定运行的商用微波回传链路MBL1,MBL2,MBL3,MBL4，其发射端空间位置分别为S1＝(xs1,ys1,zs1),S2＝(xs2,ys2,zs2),S3＝(xs3,ys3,zs3),S4＝(xs4,ys4,zs4)，接收端空间位置分别为E1＝(xe1,ye1,ze1),E2＝(xe2,ye2,ze2),E3＝(xe3,ye3,ze3),E4＝(xe4,ye4,ze4)，天线参数分别为A1,A2,A3,A4；

[0037] 2.按照预设时间间隔，计算每条链路的发射端平均功率TP1,TP2,TP3,TP4，接收端平均功率RP1,RP2,RP3,RP4；在本实施例中，以60min为时间间隔。

[0038] 3.根据接收端与发射端之间的功率差值，计算每条链路的晴空实测衰减PLi＝TPi‑RPi(i＝1,2,…,4)。

[0039] 在实际应用过程中，根据多条商用微波回传链路所记录的接收端和发射端功率，可获得一定区域内近地面平均大气折射率廓线信息，在本实施例中，商用微波回传链路分布示意如图2所示。

[0040] 然后，如图3所示，设置近地面大气折射率廓线表征参数M＝[k1,k2,k3,H1,H2]，并初‑1 ‑1 ‑1始化为k1＝‑10N‑unitkm ,k2＝10N‑unitkm ,k3＝‑10N‑unitkm ,H1＝30m,H2＝60m。

[0041] 再次，根据初始阶段得到的天线参数、链路接收端和发射端空间位置，以及上一步设置的近地面大气折射率廓线表征参数，结合电磁传播计算模型，并基于分步广角抛物方程模型，计算每条链路的晴空理论衰减PL01,PL02,PL03,PL04。在本实施例中，电磁传播计算模型包括但不限于抛物方程法。

[0042] 随后，根据晴空实测衰减和所述晴空理论衰减，计算晴空衰减目标函数。

[0043] 在本实施例中，目标函数J(M)为：

[0044]

[0045] 其中，Regularizer(M)为正则化项，在本实施例中，正则化项包括但不限于L1正则化项。

[0046] 最后，基于晴空衰减目标函数，使用寻优算法，得到近地面大气折射率廓线表征参数得到最优解，完成近地面大气折射率廓线监测，大气折射率廓线反演结果如图4所示。

[0047] 在本实施例中，近地面大气折射率廓线表征参数的最优解为上述目标函数最小的解，表达式为 寻优算法包括但不限于蚁群算法。

[0048] 以上所述的实施例仅是对本申请优选方式进行的描述，并非对本申请的范围进行限定，在不脱离本申请设计精神的前提下，本领域普通技术人员对本申请的技术方案做出的各种变形和改进，均应落入本申请权利要求书确定的保护范围内。