一种抑制虚假对流的方法

1.本发明涉及一种气象数据处理方法，尤其涉及一种抑制虚假对流的方法。


背景技术：

2.数值天气预报模型（nwps）中的初始条件对短期天气预报的技能起着重要作用。目前，雷达观测、闪电数据、卫星数据等已经通过各种方法进行同化，以改善nwps的初始场，这些研究在提高强对流过程的命中率方面取得了良好的效果。但是，在实际天气预报中，经常会出现虚假预报的情况。对于快速移动的对流系统，模拟和观测之间可能存在明显的位置偏差，此时如果同化高时空分辨率的观测数据，很可能会出现两个雨带（一个是实际的，一个是假的）。虽然降水的命中率提高了，但空报率仍然很高。还有其他一些情况，数据同化可以改善强对流过程的预报，但同时也会产生一些虚假的降水。vendrasco等人（2016）表明，使用三维变异（3dvar）技术同化雷达数据可能会产生虚假降水，以及降水位置和数量的巨大误差。这些问题在快速循环同化中可能会恶化。
3.目前已有一些关于如何抑制假性对流的研究。其中一个原因可能是在引入高时空分辨率数据的情况下，初始分析中的动力和微物理场缺乏适当的平衡。vendrasco等人（2016）试图通过在代价函数中加入大尺度分析约束来最小化高估的预报，lin等人（2021年）通过谱nudging方法对区域模型施加了一个大尺度约束以抑制虚假降水。除此之外，reen（2017）认为虚假对流可能是由于噪声造成的，可以使用数字滤波方法（lynch，1993）来试图减少它，然而，数字滤波也可能去除真正的大气特征。上述方法都能显著减少数据同化过程中引入的高估降水，但它们对不准确的驱动场导致的虚假降水带的抑制作用很小。
4.此外，一些学者试图减少虚假对流区域的水汽或水物质。fierro等人（2019）将闪电密度数据同化为雷电发生区的水汽质量混合比，然后在整个领域内去除等量的水汽质量。但结果显示，这些在闪电区域外的水汽质量的增量负调整对虚假对流的影响可以忽略不计。主要原因是只有初始时刻的水汽被轻微调整，但对动力场和热力场没有足够的调整，新的对流仍然会在虚假区域被触发，周围区域的水汽在动力场的作用下会迅速补充到水汽减少的区域，用物理初始化方法去除初始场的水汽，可能只影响初始时刻的降水，但对流会在随后的时刻迅速重新发展。
5.晴空回波或无雨回波也可用来抑制虚假的对流，目前也取得了一些比较明显的效果，但仍然存在一些问题。首先，晴空回波的观测量很大，同化的工作量也会很大。第二，晴空回波的探测精度比降水回波低，观测误差不易给出（rennie等，2018）。第三，以往的研究表明，在通过3dvar将水物质混合比作为反射率数据同化的控制变量时，由于背景反射率小，代价函数梯度极大，可能导致同化效率低下（sun等，1997），而在同化晴空回波时，这个问题会更加明显，因为晴空回波区域的背景反射率趋于小值（kong，2017）。
6.一般来说，有效抑制假对流需要同时考虑动态和热力调整。一方面，需要从热力和动力方面抑制假对流区域的对流发展。另一方面，需要削弱假对流区与周围其他区域的热和动量交换，从而限制假对流区对流的再生。
7.集合平方根滤波器（ensrf）使用集合预报来估计流依赖的背景误差协方差，它已被用于同化常规观测、雷达、卫星和闪电数据等（wang等，2015；gao和min，2018）。ensrf具有更大的灵活性，gan等人(2021)提出了一种通过ensrf方法同化闪电数据转换的二维w
max
的方案，结果表明，同化二维w
max
可以有效改善热力场，为对流的发展提供温暖湿润的环境，并可以调整动力场，使低层出现强辐合，高层出现强辐散，由于w
max
与模式的动力和热力变量之间有很强的关联性，较大的w
max
可以用来把中尺度和小尺度的信息很好地加入到背景中，并能在很大程度上减少工作量。


技术实现要素：

8.本发明所要解决的技术问题是提供一种有效降低预报空报率的抑制虚假对流的方法。
9.为解决上述问题，本发明所述的一种抑制虚假对流的方法，包括以下步骤：
⑴
确定虚假对流的范围：将观测的最大雷达回波小于10 dbz且背景场中的最大回波与观测回波的差大于10 dbz的范围确定为虚假对流的范围；
⑵
将0≤雷达回波≤10 dbz的区域确定为非降水回波区域，定义0w
max
为背景场非降水回波区最大垂直速度的平均值；
⑶
利用集合均方根滤波的方法，首先同化三维的雷达回波，然后在虚假对流区同化所述0w
max
值，循环同化四次得到分析场；
⑷
利用分析场继续预报得到预报场；
⑸
进行预报检验。
10.所述步骤
⑶
集合均方根滤波的方法中同化三维的雷达回波和同化二维的0w
max
时分析变量包括风场（u,v,w）、扰动位温（prt）、扰动位势（ph）、水汽混合比（qv）、雨水混合比（qr）、冰混合比（qi）、霰混合比（qg）、云混合比（qc）、雪混合比（qs）。
11.所述集合均方根滤波的方法中同化三维的雷达回波时采用的水平局地化方案和垂直局地化方案是gc函数。
12.所述集合均方根滤波的方法中同化二维的0w
max
采用gc函数作为水平局地化函数。
13.本发明与现有技术相比具有以下优点：1、本发明基于中尺度数值模式wrf和ensrf同化方法，通过在虚假对流区域同化0w
max
值，从动力场和热力场两方面对背景场进行调整，可以减小虚假对流区域的背景场中的水汽和上升运动，有效抑制该区域对流的发展，还可以减少周边区域向虚假对流区域的水汽输送，并且削弱地面冷池，不利于虚假对流的再生和维持，对于虚假降水和虚假回波的预报具有显著的抑制效果，也可以有效降低预报空报率。
14.2、相比于已有的同化晴空回波、或者同化非降水回波的方法，本发明提出的在虚假对流区域同化0w
max
值的方法抑制虚假降水预报的效果更加显著。
15.3、相比于已有的同化三维的晴空回波、或者同化非降水回波的方法，本发明提出的在虚假对流区域同化二维0w
max
值的方法很大程度上减小了计算量。
附图说明
16.下面结合附图对本发明的具体实施方式作进一步详细的说明。
17.图1为本发明中试验设计的流程图，其中rf、da、w
max
分别表示雷达回波、资料同化、最大垂直速度，ctl、rda、rda_0wmax分别表示控制试验、同化雷达回波的试验、以及同时同化雷达和0w
max
值的试验，spin
‑
up可以理解为：因为中小尺度模式都是由无云初始场开始的，所以开始需要一定的时间产生相应的云水信息，冷启动之后会在一段时间内产生实际的云水信息，但是这几个小时内的预报是不准确的，将这几个小时称为“spin
‑
up”时间。
18.图2为本发明同化的雷达回波和同化的0w
max
范围。其中（a1）是2018年6月27日12 utc观测的最大回波，（a2）是2018年6月27日13 utc观测的最大回波，（a3）是2018年6月27日14 utc观测的最大回波，（a4）是2018年6月27日15 utc观测的最大回波；（b1）是同化2018年6月27日12 utc的0w
max
的范围，（b2）是同化的2018年6月27日13 utc的0w
max
的范围，（b3）是同化的2018年6月27日14 utc的0w
max
的范围，（b4）是同化的2018年6月27日15 utc的0w
max
的范围。
19.图3为本发明预报的最大雷达回波和最大回波的预报性能图。（a）为观测的2018年6月27日16 utc的最大雷达回波，（b）为ctl试验预报的2018年6月27日16 utc的最大雷达回波，（c）为rda试验预报的2018年6月27日16 utc的最大雷达回波，（d）为rda_0wmax试验预报的2018年6月27日16 utc的最大雷达回波，（e）为阈值取10 dbz的预报性能图，（f）为阈值取20 dbz的预报性能图，（g）为阈值取30 dbz的预报性能图。
20.图4为本发明累积降水分布和累积降水的预报性能图。（a）为观测的2018年6月27日15 utc至2018年6月27日18 utc的累积降水，（b）为ctl试验预报的2018年6月27日15 utc至2018年6月27日18 utc的累积降水，（c）为rda试验预报的2018年6月27日15 utc至2018年6月27日18 utc的累积降水，（d）为rda_0wmax试验预报的2018年6月27日15 utc至2018年6月27日18 utc的累积降水，（e）为阈值取1 mm的预报性能图，（f）为阈值取5 mm的预报性能图，（g）为阈值取10 mm的预报性能图。
具体实施方式
21.一种抑制虚假对流的方法，包括以下步骤：
⑴
确定虚假对流的范围：10 dbz被认为是区分降水回波和非降水回波的阈值。将观测的最大雷达回波小于10 dbz（包括观测缺省的区域）且背景场中的最大回波与观测回波的差大于10 dbz的范围确定为虚假对流的范围。
22.⑵
将0≤雷达回波≤10 dbz的区域确定为非降水回波区域，定义0w
max
为背景场非降水回波区最大垂直速度的平均值；非降水回波区对流发展很弱，所以该区域的最大垂直速度是接近于0 m/s的小值。本发明中根据2018年6月27日个例背景场计算得到的0w
max
为0.2 m/s。
23.⑶
利用集合均方根滤波（ensrf）的方法，首先同化三维的雷达回波（>10 dbz），然后在虚假对流区同化0w
max
值，循环同化四次得到分析场。
24.由于最大垂直速度是一个特殊的二维的变量，ensrf具有很强的灵活性，易于建立最大垂直速度和其他变量之间建立平衡约束，此外ensrf方法避免了产生扰动观测时生成
的采样误差，所以本发明中选用ensrf同化闪电转换的最大垂直速度。
25.集合均方根滤波的方法中同化三维的雷达回波和同化二维的0w
max
时分析变量包括风场（u,v,w）、扰动位温（prt）、扰动位势（ph）、水汽混合比（qv）、雨水混合比（qr）、冰混合比（qi）、霰混合比（qg）、云混合比（qc）、雪混合比（qs）。其中：集合均方根滤波的方法中同化三维的雷达回波时采用的水平局地化方案和垂直局地化方案是gc函数（gaspari和cohn，1999）。因为弱对流区垂直速度的垂直变化不明显，集合均方根滤波的方法中同化二维的0w
max
采用gc函数作为水平局地化函数，没采用垂直局地化函数。
26.⑷
利用分析场继续预报得到预报场。
27.⑸
进行预报检验。
28.为了验证本发明对强对流预报的改善效果，选取了2018年6月27日发生在我国江苏和安徽地区的一次强对流过程进行个例分析，本次试验总共设计一组控制实验以及两组循环同化试验，试验设计的流程图如图1所示。图1中rf、da、w
max
分别表示雷达回波、资料同化、最大垂直速度，ctl、rda、rda_0wmax分别表示控制试验、同化雷达回波的试验、以及同时同化雷达和0w
max
的试验。spin
‑
up可以理解为：因为中小尺度模式都是由无云初始场开始的，所以开始需要一定的时间产生相应的云水信息，冷启动之后会在一段时间内产生实际的云水信息，但是这几个小时内的预报是不准确的，将这几个小时称为“spin
‑
up”时间。
29.详细的试验设计如表1所示：表1 试验设计表在上述试验中，试验ⅰ（ctl）是控制试验，试验ⅱ（rda）是利用ensrf方法同化雷达回波的试验，同化时间是2018年6月27日12、13、14、15 utc，试验ⅲ（rda_0wmax）是利用ensrf方法同化雷达回波和虚假对流区的0w
max
的试验。本研究选用wrf 4.0版本，21个全球集合成员预报系统的集合预报被用为wrf模式提供初始场和边界场，为了使样本离散度增强，选用rrtmg和goddard辐射参数化方案使21个样本扰动成为42个样本，除此之外，选用的参数化方案有thompson微物理方案、monin
‑
obukhov近地层方案、mellor
‑
yamada
‑
janjic tke行星边界层方案、unified noah陆面方案，tiedtke积云参数化方案（只在d01中使用）。本研究采用两重嵌套网格，d01的水平分辨率为9 km，模式网格为171
×
191，d02水平分辨率为3 km，模式网格为211
×
271，垂直方向有50层。
30.本发明采用ensrf方案同化雷达回波和虚假对流区的二维0w
max
值，图2所示是四个同化时刻的雷达回波和0w
max
值的分布情况。仅同化雷达回波的试验（rda）存在比较明显的虚假对流，主要位于真实对流的周围。
31.图3（a
‑
d）和图4（a
‑
d）显示了雷达最大反射率和累积降水的分布。2018年6月27日16 utc，有两个强回波区域：一个在安徽和江苏两省交界处，另一个在海洋上空。rda试验很好地模拟了强反射率，但反射率的范围比观测回波要大。rda_0wmax试验预报的反射率的范围小于rda试验。在降水方面，rda实验在江苏中部出现了明显的虚假降水，但rda_0wmax试验模拟的降水范围与观测值接近。
32.图3（e
‑
g）和图4（e
‑
g）分别是不同阈值的雷达回波和累积降水的预报性能图，预报性能图可以反映很多信息，包括预报偏差（fr，图中虚线所示，越接近1表明预报偏差越小，大于1表示高估，小于1表示低估）、命中率（pod，越高表明预报越好）、临界成功指数（csi，图中曲线所示，范围0
‑
1，约接近于1越好）、成功率（sr，1减去空报率，越接近于1越好），综合来说在预报性能图中约接近于右上角，表明预报效果越准确。反射率的预报性能图显示rda_0wmax实验的空报率小于rda实验，降水的预报性能图也可以得到同样的结论。综上所述，通过ensrf在虚假对流区域同化0w
max
可以有效地抑制虚假对流。