一种适用于小雨量级的gnss极化信号正演方法及系统
阅读说明：本技术 一种适用于小雨量级的gnss极化信号正演方法及系统 (GNSS polarized signal forward modeling method and system suitable for small rainfall magnitude ) 是由 白伟华 孙越强 苏豆豆 杜起飞 刘黎军 李伟 王先毅 蔡跃荣 曹光伟 夏俊明 孟 于 2021-07-06 设计创作，主要内容包括：本发明公开了一种适用于小雨量级的GNSS极化信号正演方法及系统,所述小雨的降雨率小于等于2.5毫米/小时,所述方法包括：获取待测小雨雨区内的降雨信息；所述小雨雨区有GNSS极化无线电掩星事件发生,并且掩星的射线路径穿过雨区；分别建立雨滴形状模型和雨滴谱分布模型；计算雨滴的相对复介电常数；结合雨滴形状模型计算单一雨滴水平极化前向散射幅度和垂直极化前向散射幅度；根据雨滴谱分布模型和单一雨滴水平极化前向散射幅度及垂直极化前向散射幅度,计算差分相位延迟；根据差分相位延迟对整个雨区内射线总路径长度进行积分,得到正演后的极化相移。(The invention discloses a GNSS polarized signal forward modeling method and system suitable for light rain magnitude, wherein the rainfall rate of the light rain is less than or equal to 2.5 mm/h, and the method comprises the following steps: acquiring rainfall information in a small rain area to be detected; the small rain area has a GNSS polarized radio occultation event, and a ray path of occultation passes through the rain area; respectively establishing a raindrop shape model and a raindrop spectrum distribution model; calculating the relative complex dielectric constant of the raindrops; calculating the horizontal polarization forward scattering amplitude and the vertical polarization forward scattering amplitude of the single raindrop by combining the raindrop shape model; calculating differential phase delay according to a raindrop spectrum distribution model and the single raindrop horizontal polarization forward scattering amplitude and the single raindrop vertical polarization forward scattering amplitude; and integrating the total path length of rays in the whole rain area according to the differential phase delay to obtain the forward polarization phase shift.)
技术领域
本发明涉及降雨遥感探测研究领域，具体涉及到适用于小雨量级降雨事件的GNSS极化信号的正演方法及系统。
背景技术
降雨是研究气候变化的基本变量之一，其对农业生产、生态环境以及人们的日常生活有着重要的作用。目前国际上探测降雨的手段多种多样，有天基降水雷达系统如带有降水雷达的热带降雨测量任务TRMM、带有双频降水雷达的全球降水测量任务GPM，此外还有高光谱红外探测仪、微波探测仪、地基降水雷达等。然而这些手段均无法从观测值中直接或间接提取到影响水蒸汽或降水的因素，均不适合观测强降水云的内部结构。全球卫星导航系统(global navigation satellite system,GNSS)的无线电掩星(radio occultation,RO)技术根据无线电链路传播时的附加时延和弯曲信息反演得到折射率、温度、密度等大气物理参量廓线，是进行行星电离层和大气层探测的一种重要手段。与探空仪或气象卫星等相比，该技术具有全天候、精度高、垂直分辨率高、系统误差低等优点，且目前有研究表明GNSS RO技术可以利用极化无线电掩星(polarimetric radio occultation,PRO)的潜在能力，共同探测和量化强降水事件中的大气热力学廓线和降雨信息。
为有效利用GNSS RO这一丰富可靠的数据源，提高对降雨事件预报的准确性，基于信号通过大尺寸的非球形水凝物会引起极化差分相位延迟的事实，目前国内还未发展此类气象业务，对GNSS极化无线电掩星探测和量化降雨的理论和方法尚在探索当中，对其反演降水的理论和方法还未发展成熟。
发明内容
本发明的目的在于克服现有技术缺陷，提出了一种适用于小雨量级的GNSS极化信号正演方法及系统。
为了实现上述目的，本发明提出了一种适用于小雨量级的GNSS极化信号正演方法，所述小雨的降雨率小于等于2.5毫米/小时，所述方法包括：
获取待测小雨雨区内的降雨信息；所述小雨雨区有GNSS极化无线电掩星事件发生，并且掩星的射线路径穿过雨区；
根据降雨信息分别建立雨滴形状模型和雨滴谱分布模型；
根据降雨信息计算雨滴的相对复介电常数；
结合雨滴形状模型计算单一雨滴水平极化前向散射幅度和垂直极化前向散射幅度；
根据雨滴谱分布模型和单一雨滴水平极化前向散射幅度及垂直极化前向散射幅度，计算差分相位延迟；
根据差分相位延迟对整个雨区内射线总路径长度进行积分，得到正演后的极化相移。
作为上述方法的一种改进，所述待测小雨雨区内的降雨信息包括：雨区范围、沿射线路径的降雨率信息、频率、等体积直径、总路径长度、空气温度和雨滴倾斜角度。
作为上述方法的一种改进，所述雨滴形状模型采用SC雨滴形状模型；具体包括：
设置雨滴形状为扁椭球体，轴比和等体积直径D满足下式：
其中，b为短半径，a为长半径，D为等体积直径。
作为上述方法的一种改进，所述雨滴谱分布模型采用MP分布雨滴谱分布模型或JD雨滴谱分布模型；当采用MP分布雨滴谱分布模型时，雨滴谱分布N(D)满足下式：
N(D)＝N0e-ΛD N0＝8000,Λ＝4.1R-0.21
当采用JD雨滴谱分布模型时，雨滴谱分布N(D)满足下式：
N(D)＝N0e-ΛD N0＝30000,Λ＝5.7R-0.21；
其中，N0为常数，Λ为与降雨率有关的参数，R为降雨率，单位为mm/h，D为雨区内雨滴的等体积直径。
作为上述方法的一种改进，所述计算雨滴的相对复介电常数ε，满足下式：
ε＝(ε0-ε1)/[1-i(f/f1)]+(ε1-ε2)/[1-i(f/f2)]+ε2
其中，ε的实数部分ε′，虚数部分ε″满足下式：
其中，
ε0＝77.66-103.3·T1
ε1＝0.0671·ε0
ε2＝3.52
f1＝20.20+146.4·T1+316·T1 2
f2＝39.8·f1
T1＝1-300/T
其中，f为频率，f1为弛豫频率，f2为与弛豫频率相关的系数，ε0为静介电常数，ε1为与静介电常数相关的系数，ε2为常数，T为大气温度，T1为与大气温度相关的系数。
作为上述方法的一种改进，所述根据雨滴谱分布模型和单一雨滴水平极化前向散射幅度及垂直极化前向散射幅度，计算差分相位延迟；具体包括：
根据所选择的雨滴谱分布模型和计算得到的水平极化前向散射幅度fh(D)及垂直极化前向散射幅度fv(D)，采用下式计算差分相位延迟Kdp为：
其中，表示取fh(D)-fv(D)的实数部分，N(D)为雨滴谱分布，λ为GNSS极化信号波段波长，σ倾斜角分布标准差，D为雨区内雨滴的等体积直径，θ为雨滴倾斜角度。
作为上述方法的一种改进，所述根据差分相位延迟对整个雨区内射线总路径长度进行积分，得到正演后的极化相移；具体包括：
在雨区内沿RO射线路径上每个点的差分相位延迟在整个雨区内射线路径总长度L上进行积分，得到正演后的极化相移ΔΦ为：
ΔΦ＝∫LKdp(l)dl
其中，Kdp(l)为单个点l的差分相位延迟，ΔΦ的单位为mm。
一种用于小雨量级的GNSS极化信号正演系统，所述小雨的降雨率小于等于2.5毫米/小时，其特征在于，所述系统包括：数据获取模块、模型建立模块、相对复介电常数计算模块、前向散射幅度计算模块、差分相位延迟计算模块和极化相移计算模块；其中，
所述数据获取模块，用于获取待测小雨雨区内的降雨信息；所述小雨雨区有GNSS极化无线电掩星事件发生，并且掩星的射线路径穿过雨区；
所述模型建立模块，用于根据降雨信息分别建立雨滴形状模型和雨滴谱分布模型；
所述相对复介电常数计算模块，用于根据降雨信息计算雨滴的相对复介电常数；
所述前向散射幅度计算模块，用于结合雨滴形状模型计算单一雨滴水平极化前向散射幅度和垂直极化前向散射幅度；
所述差分相位延迟计算模块，用于根据雨滴谱分布模型和单一雨滴水平极化前向散射幅度及垂直极化前向散射幅度，计算差分相位延迟；
所述极化相移计算模块，用于根据差分相位延迟对整个雨区内射线总路径长度进行积分，得到正演后的极化相移。
与现有技术相比，本发明的优势在于：
1、本发明具备操作简单、计算快速、准确性高、精度高等优点，能够满足小雨量级降雨事件正演极化相移的需求；
2、本发明采用T-matrix法可以计算任何形状的降雨粒子，因此具备较高的普适性；
3、本发明可以方便进行大规模降雨正演模拟的计算，可以满足对大量不同形状的雨滴粒子进行前向散射特性计算的需求；
4、通过本发明的方法，可以利用现有的降雨产品数据集进行正演模拟，在不增加任何成本和观测的条件下，完成对降雨区的正演，对降雨气象研究领域具有重要贡献。
附图说明
图1是本发明适用于小雨量级的GNSS极化信号正演方法流程图；
图2是本发明所需准备的小雨雨区相关信息图；
图3本发明实施例中2018/12/19小雨雨区模拟图，其中，图3(a)为各个降雨层的降雨率分布，Layer2所在高度为4km，Layer1所在高度为2km，Surface为地表面，图3(b)为RO事件穿过各个高度层的雨区；
图4是本发明实施例中2018/12/19 16:54发生的RO事件射线切点轨迹变化图。
具体实施方式
本发明的目的在于提供一种适用于小雨量级的GNSS极化信号正演方法，其具备简单、快速、适合大规模计算、模拟效果好、精度高、普适性高等优点，能够满足利用GNSS极化信号对小雨量级降雨事件进行正演模拟的需求。
为完善GNSS极化测降雨的理论，探索GNSS PRO技术量化降雨的方法，本发明提出了一种适用于小雨量级降雨事件的GNSS极化信号正演方法，该方法的模型和算法主要使用于小雨量级条件下的GNSS极化信号正演，具体指降雨率小于等于2.5毫米/小时。该方法计算简单快速且准确性高，具备较高的普适性，并方便进行大规模计算，有助于开展GNSS PRO技术量化探测降雨参数机理深入研究和指标误差传递量化分析，指导后续载荷方案设计、测试和研发，以及为后续数据反演奠定基础。
为实现上述目的，本发明提出了一种适用于小雨量级的GNSS极化信号正演方法，主要包括以下7个步骤：
数据准备，雨滴形状模型选择，雨滴谱模型选择，雨滴相对复介电常数计算，前向散射幅度计算，差分相位延迟计算和极化相移计算。
所述的数据准备，是指获取所需研究的小雨雨区内的降雨信息，包括沿射线路径的降雨率信息、等体积直径、总路径长度、空气温度、雨滴倾斜角度等信息。需注意的是，所要研究的小雨雨区须有GNSS极化无线电掩星事件发生，且掩星的射线路径要穿过雨区。
所述的雨滴形状模型选择，是指选择Steinert-Chandra(SC)雨滴形状模型作为小雨雨区的雨滴模型。将雨滴形状看作扁椭球体，其模型主要由轴比和等体积直径的关系来描述，表达式为
其中，b为短半径，a为长半径，D为等体积直径。
所述的雨滴谱模型选择，是指选择Marshall–Palme(MP)分布或Joss Drizzle(JD)雨滴谱分布模型作为小雨雨区的雨滴谱模型，其表达式及经验参数如下表所示：
所述的雨滴相对复介电常数计算，是指将雨滴近似为纯水来计算其相对复介电常数。本发明采用Liebe1993模型公式来计算雨滴的相对复介电常数，公式如下：
ε＝(ε0-ε1)/[1-i(f/f1)]+(ε1-ε2)/[1-i(f/f2)]+ε2
式中，
ε0＝77.66-103.3·T1
ε1＝0.0671·ε0
ε2＝3.52
f1＝20.20+146.4·T1+316·T1 2
f2＝39.8·f1
T1＝1-300/T
其中，f为频率，f1为弛豫频率，f2为与弛豫频率相关的系数，ε0为静介电常数，ε1为与静介电常数相关的系数，ε2为常数，T为大气温度，T1为与大气温度相关的系数。
所述的前向散射幅度计算，是指采用T-matrix法结合所选择的雨滴形状模型计算单一雨滴水平极化前向散射幅度fh(D)和垂直极化前向散射幅度fv(D)。采用T-matrix法不仅计算快速且结果精确，适合大规模计算任何形状的非球形雨滴。T-matrix法不属于本专利范畴，不予详述。
所述的差分相位延迟计算，是指根据所选择的小雨雨滴谱分布模型和计算得到的水平极化前向散射幅度fh(D)与垂直极化前向散射幅度fv(D)采用下式进行计算，
其中，代表取fh(D)-fv(D)的实数部分，N(D)为雨滴谱分布，λ为GNSS极化信号波段波长，D为雨区内雨滴的等体积直径，θ为雨滴倾角。
所述的极化相移计算，主要是指雨区内沿RO射线路径上每个点的差分相位延迟在整个雨区内射线路径总长度L上的积分，公式如下：
ΔΦ＝∫LKdp(l)dl
其中，Kdp为单个点的差分相位延迟，ΔΦ为极化相移，单位为mm。
经过以上步骤处理即可得到所研究小雨雨区正演模拟的极化相移。
下面结合附图和实施例对本发明的技术方案进行详细的说明。
实施例1
如图1所示，本发明实施例的一种适用于小雨量级的GNSS极化信号正演方法，主要由数据准备、雨滴形状模型选择、雨滴谱模型选择、雨滴相对复介电常数计算、前向散射幅度计算、差分相位延迟计算、极化相移计算七大步骤组成。
本实施例中所用降雨产品为GPM DPR，所用极化掩星产品为PAZ PRO数据。将两种产品进行匹配后，筛选得到发生过RO事件的小雨雨区。以下以2018/12/1916:54掩星事件发生时所穿过的小雨雨区进行正演模拟为例实施该发明中的适用于小雨量级的GNSS极化信号正演方法。
本实施例中步骤一，数据准备。根据图2获取小雨雨区的相关信息，包括降雨率、空气温度、频率、雨区范围等信息。其中，降雨率、空气温度信息来自GPM DPR日产品，等体积直径来自DPR 2A级轨道产品，频率为L1波段频率，大约为1.57542GHz，雨滴倾角设置为0°，雨区范围由掩星事件发生的开始切点位置和结束切点位置所决定，切点位置信息由PAZ PRO数据产品中GPS位置、LEO位置及时间信息进行计算。该小雨雨区地理范围广，覆盖北纬34°-38°，垂直分布可达4km高度，降雨带以及延伸范围内降雨单元的最高降雨率为5.0298mm/h，整个区域内各个降雨层的平均降雨率分别为1.0070、1.1456和0.8432mm/h，沿射线路径平均降雨率为1.4317mm/h，整体降雨率偏小。GPM卫星扫描雨区时刻为16:16，时间相差不到一小时。经过可视化后，得到该小雨雨区模拟图及RO事件穿越雨区模拟图(见图3)。可知每层降雨信息丰富且均与RO事件相交，RO射线所经过雨区得降雨率均在3mm/h以下。图4为本发明实施例中2018/12/19 16:54发生的RO事件射线切点轨迹变化图，加粗线段代表20km以下高度的RO射线切点轨迹，可见RO发生时各条射线的切点轨迹变化大致可呈现出一条曲线。
本实施例中步骤二，雨滴形状模型选择。主要是选择SC模型作为实施例雨区的雨滴形状模型。
本实施例中步骤三，雨滴谱模型选择。分别选择MP分布、JD分布作为实施例雨区的雨滴谱分布模型进行正演。
本实施例中步骤四，雨滴相对复介电常数计算。采用Liebe1993模型公式计算雨滴的相对复介电常数。
本实施例中步骤五，前向散射幅度计算。采用T-matrix法结合SC雨滴形状进行单一雨滴的水平极化和垂直极化前向散射幅度计算。
本实施例中步骤六，差分相位延迟计算。结合MP雨滴谱分布模型和步骤五计算得到的水平极化和垂直极化前向散射幅度来计算群雨滴的差分相位延迟。结合JD雨滴谱分布模型和步骤五计算得到的水平极化和垂直极化前向散射幅度来计算群雨滴的差分相位延迟。
本实施例中步骤七，极化相移计算。根据得到的差分相位延迟对整个雨区内射线总路径长度进行积分，得到正演后的极化相移。
采用本方法对该小雨雨区的GNSS极化信号进行正演，得到MP雨滴谱分布正演得到的极化相移为5.6636mm，JD雨滴谱分布正演得到的极化相移为5.8649。PAZ卫星实测数据天线相位校正的极化相移为6.0154mm。两个正演模拟结果与实测数据之间相差不大，分别为0.3518mm、0.1505mm，较为接近，说明了本发明适用于小雨量级的GNSS极化信号正演方法的准确性和可靠性。
实施例2
一种用于小雨量级的GNSS极化信号正演系统，所述小雨的降雨率小于等于2.5毫米/小时，所述系统包括：数据获取模块、模型建立模块、相对复介电常数计算模块、前向散射幅度计算模块、差分相位延迟计算模块和极化相移计算模块；其中，
所述数据获取模块，用于获取待测小雨雨区内的降雨信息；所述小雨雨区有GNSS极化无线电掩星事件发生，并且掩星的射线路径穿过雨区；
所述模型建立模块，用于根据降雨信息分别建立雨滴形状模型和雨滴谱分布模型；
所述相对复介电常数计算模块，用于根据降雨信息计算雨滴的相对复介电常数；
所述前向散射幅度计算模块，用于结合雨滴形状模型计算单一雨滴水平极化前向散射幅度和垂直极化前向散射幅度；
所述差分相位延迟计算模块，用于根据雨滴谱分布模型和单一雨滴水平极化前向散射幅度及垂直极化前向散射幅度，计算差分相位延迟；
所述极化相移计算模块，用于根据差分相位延迟对整个雨区内射线总路径长度进行积分，得到正演后的极化相移。
最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。
- 上一篇：一种医用注射器针头装配设备
- 下一篇：基于周期性圆环-抛物线型复合孔阵列的数字超表面