技术领域
[0001] 本
发明属于雷达回波
信号仿真技术领域，特别是涉及一种机载气象雷达低空风切变三维回波仿真方法。
背景技术
[0002] 低空风切变是飞机在
起飞和着陆阶段遇到的危害最大的天气现象之一。当飞机穿过风切变区域的风场中心时，受急剧变化的风速和下沉气流的影响，会使飞机接近
失速。由于在飞机起降阶段可操纵的高度空间受限，若风切变强度较大，则极易造成严重的坠机事故。由于真实的风切变回
波数据较难获取，因此开展高保真的低空风切变三维回波仿真方法的研究具有十分重要的现实意义。
[0003] 气象目标模型的建立不同于一般的目标建模。常见的目标建模，大多是基于目标的起伏特性在时间上的相关性来建立相应的统计模型，而对于气象目标模型来说明，不仅在时间维具有相关性，在空间维也具有相关性，因此不能简单地采用统计模型的方法。
[0004] 低空分切变是一种分布式气象目标，常用的低空风切变建模方法主要分为3类：
[0005] (1)类似于一般目标的外场实测法，利用以由机场的气象雷达网和塔台组成的监测网实际测量数据为依据建立的模型，如美国的联合机场天气研究(Joint Airport Weather Studies，JAWS)计划等。该类模型基于探测到的真实数据进行建模，模型数据真实可靠，但其无法反映出风切变的本质特征及动态发展过程，不能随气象条件的不同而变化。
[0006] (2)工程化的简化模型，该模型通常根据物理概念用简单的数学形式拟合而成。虽然该类模型实现较为简单，但其忽略了其它因素对雷达回波的影响，因而并不能体现出风切变风场的全部特点。
[0007] (3)以描述小尺度气象的大
气动力学和微物理方程为
基础的小尺度气象模型。该类模型有NASA的Langley研究中心建立的TASS(Terminal Area Simulation System)模型、简化解析模型以及ARPS(Advanced Regional Prediction System)模型。TASS模型能够模拟整个风场从发生到发展的动态全过程，不但能够反映风场的速度特征，还能提供大气中的相应参数，如
温度、
水分和反射率因子等。
[0008] 机载气象雷达回波信号仿真越来越得到重视，一方面，随着雷达的不断发展，功能愈加强大，种类愈加繁多，利用外场实验对雷达性能和指标的测试已经不能满足需求；另一方面，对于低空风切变的研究，其具有发生突然、持续时间短、影响强度大、危险性高和不易探测等特点，因此实测数据不易获取。那么，建立合理的
算法模型，利用
软件或
硬件对雷达回波信号仿真就成为重要的手段。
发明内容
[0009] 为了解决上述问题，本发明的目的在于真实低空风切变回波数据难以获取的情况下，提供一种机载气象雷达低空风切变三维回波仿真方法。
[0010] 为了达到上述目的，本发明提供的机载气象雷达低空风切变三维回波仿真方法包括以下几个步骤：
[0011] 1)利用包括场景范围、边界条件在内的场景参数，建立
三维网格化的计算域；
[0012] 2)对上述建立的计算域进行
流体仿真，生成原始风场数据；
[0013] 3)采用反距离加权插值法对上述步骤2)中生成的原始风场数据进行修正，以使原始风场数据均匀化；
[0014] 4)建立仿真场景，初始化包括雷达参数、载机参数、扫描参数和场景参数在内的场景仿真参数，并按照场景参数的设置读取步骤3)中的修正风场数据；
[0015] 5)根据步骤4)中建立的仿真场景，模拟机载气象雷达扫描低空风切变，对雷达波束内的散射点进行相干
叠加，形成一个扫描方位的雷达回波信号。
[0016] 6)利用步骤4)中的载机参数、场景参数和扫描参数分别更新飞机
位置信息、修正风场数据和雷达波束扫描
角度信息；
[0017] 7)依据雷达参数中的扫描范围判断扫描是否完成，若完成，输出一个扫描周期的回波信号；若未完成，返回步骤5)直至扫描结束。
[0018] 在步骤1)中，所述的场景范围为中心对称的圆柱形，边界条件为一个标准
大气压的开放性出口，不限制流体流出方向，三维网格化的计算域由非均匀网格构成。
[0019] 在步骤2)中，所述的对上述建立的计算域进行流体仿真，生成原始风场数据的方法是：设定由流体特性和边界条件构造的流体仿真的初始条件，利用
迭代连续方程与动量方程计算三维网格内的原始风场数据，进行流体仿真。
[0020] 所述的初始条件中以参考
密度为1.225kg/m3的理想气体设定仿真的流体特性，初始条件中依据的边界条件为一个标准大气压的开放性出口，不限制流体流出方向。
[0021] 在步骤3)中，所述的反距离加权插值法是将插值点与已知点间距离的倒数作为权值进行加权求和，其表达式如下：
[0022]
[0023] 式中，Z为插值点位置的待估参数，比如密度、速度，Zi表示第i个已知点的参数，di为第i个已知点与插值点的距离，p为幂指数，控制权值的大小，N表示加权计算用到的已知点的总数。
[0024] 在步骤5)中，所述的对雷达波束内的散射点进行相干叠加，形成一个扫描方位的雷达回波信号的方法是：利用逐脉冲撒点法实现，通过计算雷达扫描方向上各散射点的回波，叠加后形成该方位的雷达回波信号。
[0025] 在步骤6)中，所述的利用步骤4)中的载机参数、场景参数和扫描参数分别更新飞机位置信息、修正风场数据和雷达波束扫描角度信息的方法是：利用飞行速度更新飞机位置状态，利用风场风速更新修正风场数据，并按照机载气象雷达的扫描参数更新波束扫描角度信息。
[0026] 在步骤7)中，所述的扫描范围指的是依据机载气象雷达的工作模式，工作于风切变模式下的雷达扫描范围为±90°。
[0027] 本发明提供的机载气象雷达低空风切变三维回波仿真方法具有如下优点和积极效果：与
现有技术相比，本发明针对低空风切变回波数据难以直接获取的问题，基于TASS模型的设计思想，结合低空风切变的气象特性，利用计算
流体力学仿真的方法建立较为真实的三维风场模型，从而实现高保真的机载气象雷达低空风切变三维回波仿真，可为后续
数据处理和低空风切变检测技术的研究等提供可靠的数据支持。本发明方法可真实地反映出风切变回波信号的气象特征，仿真数据有效可靠。
附图说明
[0028] 图1为计算域模型示意图；
[0029] 图2为三维风场的速度分布图；
[0030] 图3为三维风场垂直切面的速度分布图；
[0031] 图4为三维风场垂直切面的密度分布图；
[0032] 图5为提供的机载气象雷达低空风切变三维回波仿真方法
流程图；
[0033] 图6为雷达回波信号的功率分布图；
[0034] 图7为雷达回波信号的速度分布图；
[0035] 图8为方位角为0°的雷达回波信号速度分布图。
具体实施方式
[0036] 下面结合附图和具体
实施例对本发明提供的机载气象雷达低空风切变三维回波仿真方法进行详细说明。
[0037] 如图5所示，本发明提供的机载气象雷达低空风切变三维回波仿真方法包括按顺序进行的下列步骤：
[0038] 1)利用包括场景范围、边界条件在内的场景参数，建立三维网格化的计算域；
[0039] 微下击暴流是一种特殊的低空风切变，本发明方法采用微下击暴流模型建模，如图1所示，构造的场景范围为中心对称的圆柱形，圆柱半径为3000米，圆柱高度为3000米。边界条件为开放式出口，圆柱顶端设置流体入口，入口为圆形，半径为500m，其余部分为流体出口区域，不限制流体流出方向，场景压强为1个标准大气压，入口流体为0℃，仿真场景为25℃。三维网格化的计算域采用非均匀网格拟合，生成计算
节点。
[0040] 2)对上述建立的计算域进行流体仿真，生成原始风场数据；
[0041] 流体仿真的初始条件由流体特性和边界条件构造，初始条件中以参考密度为1.225kg/m3的理想气体设定仿真的流体特性，初始条件中依据的边界条件即为步骤1中所述的边界条件。利用迭代连续方程与动量方程计算三维网格内的风场模型数据，进行流体仿真。设置迭代收敛因子为0.7，建立风切变的下沉气流的迭代连续方程、动量方程，通过该方程构造仿真模型。根据连续性条件(
质量守恒定律)，流出量与流入量的差应等于单位体积内流体质量的变化，从而得到连续方程：
[0042]
[0043] 式中，ρ为大气密度，t为时间，为哈密顿运算符，V为速度矢量，Sm是源项。
[0044] 流体内微元体上满足动量守恒，动量方程如下：
[0045]
[0046] 式中，p为微元体上的流体压力，V为速度矢量，t为时间，为哈密顿运算符，g为作用在微元体上的重力，F为其他外部体积力，τ为分子粘性作用下产生的作用在微体元表面上的粘性
应力张量。
[0047] 流体仿真后，获得的三维风场速度分布如图2所示，三维风场垂直切面的速度分布如图3所示。由图2-图3可见，生成的三维风场速率明显呈现出水平涡旋的特性，在水平涡旋的下方，速率出现最大值，周围的速度出现梯度递减，水平方向上风速呈现出反“S”特性，符合低空风切变的速度分布特点。
[0048] 在25℃时，一个标准大气压条件下生成的三维风场垂直切面的密度分布，如图4所示。由图4可以看出，低空风切变撞击地面时形成高密度区域，在涡旋中心处形成低密度区，随着风场的扩散，密度缓慢回到平均水平，达到与周围密度一致，呈现出高-低-高的变化趋势，本发明的流体仿真结果符合低空风切变的密度分布特点。
[0049] 3)采用反距离加权插值法对上述步骤2)中生成的原始风场数据进行修正，以使原始风场数据均匀化；
[0050] 由于步骤1)中的网格采用非均匀网格拟合，因此产生的原始风场数据不便于后续处理，本发明方法采用插值预处理方式将生成的原始风场数据均匀化。均匀化的方式采用反距离加权法，反距离加权插值法是将插值点与已知点间距离的倒数作为权值进行加权求和，其表达式如下：
[0051]
[0052] 式中，Z为插值点位置的待估参数(如密度、速度等)，Zi表示第i个已知点的参数，di为第i个已知点与插值点的距离，p为幂指数，控制权值的大小，N表示加权计算用到的已知点的总数。
[0053] 4)建立仿真场景，初始化包括雷达参数、载机参数、扫描参数和场景参数在内的场景仿真参数，并按照场景参数的设置读取步骤3)中的修正风场数据；
[0054] 雷达参数包括脉冲重复
频率为7000Hz，发射脉宽为1μs，天线增益为30dB，雷达发射功率为100W，
信噪比为5dB；载机参数包括飞机高度为600米，飞行速度为75m/s；扫描参数按照现有机载气象雷达工作在风切变模式下设置，天线扫描的角度范围为±90°，天线下视角为2.5°，扫描
角速度为 场景参数假设飞机处于(0，0)位置，风场中心处于飞机正前方6公里处，并按照该设置读取步骤3)中的修正风场数据；
[0055] 5)根据步骤4)中建立的仿真场景，模拟机载气象雷达扫描低空风切变，对雷达波束内的散射点进行相干叠加，形成一个扫描方位的雷达回波信号。
[0056] 低空风切变是分布式目标，具有分布式目标的特性，每个距离单元中由大量的散射点组成。计算每个散射点的雷达回波，通过相干叠加生成每个距离单元的回波信号。将一个波束范围内所有散射点的回波信号按上述方法进行叠加，即得到一个扫描方位的雷达回波信号，其表达式为：
[0057]
[0058] 式中，Ai(t)为第i个散射点的回波幅度， 为第i个散射点的回波
相位，N0为接收机噪声。
[0059] 根据气象雷达方程推导出回波幅度公式：
[0060]
[0061] 式中，Ai为第i个散射点的幅度，P为雷达发射功率，G为天线增益，Ze为风场反射率因子，Ri(t)为t时刻第i个散射点与雷达之间的距离，λ为发射信号
波长。
[0062] 单个散射点回波相位可以表示为：
[0063]
[0064] 式中，λ为发射信号波长，Ri为飞机与第i个散射点间的距离。回波信号的相位中包含散射点的速度信息。
[0065] 6)利用步骤4)中的载机参数、场景参数和扫描参数分别更新飞机位置信息、修正风场数据和雷达波束扫描角度信息；
[0066] 利用飞行速度更新飞机位置状态，利用风场风速更新修正风场数据，并按照机载气象雷达的扫描参数更新波束扫描角度信息。
[0067] 7)依据雷达参数中的扫描范围判断扫描是否完成，若完成，输出一个扫描周期的回波信号；若未完成，返回步骤5)直至扫描结束。
[0068] 依据机载气象雷达的工作模式，工作于风切变模式下的雷达扫描范围为±90°。图6为雷达回波信号的功率分布图，如图6所示，风切变区域位于扫描中心的正前方，扫描点距离风场中心位置约为6km，存在的风切变回波区域与构建的风场模型相匹配。图7为完整扫描范围内雷达回波信号的速度分布图。以一个方位下的速度估计结果为例，飞机正前方(方位角为0°)回波信号的风速分布如图8所示。由图7-图8可见，风切变在飞机正前方约3km到
9km的范围。在风场中心位置附近，风速估计结果呈对称分布，且符合真实风切变信号沿距离方向呈反“S”特性的分布特点，在风切变区域之外，杂乱的速度分布是由于噪声引起的。
[0069] 本发明以TASS模型的设计思想为基础，结合气象目标的微物理特性，利用计算流体力学的仿真方法构造低空风切变的三维风场模型，在对非均匀的原始场景数据进行预处理之后，利用逐脉冲撒点仿真法，模拟机载气象雷达扫描的实际过程，实现了高保真的低空风切变三维回波仿真。
