一种大气温度廓线预测方法

1.本发明涉及大气温度预测技术领域，具体涉及一种大气温度廓线预测方法。


背景技术：

2.大气廓线数据除了直接应用在天气预报，气候变化监测等方面，在定量遥感和光学工程中亦有广泛应用，因为大气廓线从根本上决定了大气的光学特性，可用来预测大气透过率和大气背景辐射。大大气参数空间分布结构复杂，根据不同的参数特性从近地面到太空可划分为对流层，平流层，中间层，热层和逸散层。
3.大气温度空间分布数据的预测模型可以看作是时间序列的预测模型。虽然目前已经有部分温度预测的相关研究，但由于现有的温度预测方法大都基于季节性差分自回归滑动平均模型(seasonal autoregressive integrated moving average,sarima)实施的，而sarima模型是线性模型，在处理线性相关关系较弱的平流层数据时，预测精度较低，平流层大气温度数据的预测成为了整个大气温度廓线预测模型的短板。


技术实现要素：

4.本发明的目的在于提供一种大气温度廓线预测方法，以解决现有的大气温度预测方法在处理线性相关关系较弱的平流层数据时预测精度低的问题。
5.本发明解决上述技术问题的技术方案如下：
6.本发明提供一种大气温度廓线预测方法，所述大气温度廓线预测方法包括：
7.s1：获取多个测量高度层；
8.s2：分别获取多个所述测量高度层的中每层的测量数据；
9.s3：将每层所述测量数据分别输入sarima模型中，得到与每层所述测量数据相对应的各层预测数据；
10.s4：根据所有的预测层数据，得到一条大气廓线；
11.s5：将所述大气廓线输入降噪自编码器中，得到预测大气温度廓线。
12.可选择地，所述步骤s1包括：
13.s11：获取目标测量点的大气温度廓线数据；
14.s12：对所述目标测量点的大气温度廓线数据进行编号，得到多个测量高度层。
15.可选择地，所述步骤s3中，所述sarima模型为：
[0016][0017]
其中，φ(
·
)和φ(
·
)表示自回归多项式，且φ(z)＝1-φ1z-φ2z
2-...-φ
pzp
，φ(z)＝1-φ1z-φ2z
2-...-φ
pzp
；b和bs是延迟算子，且bx
t
＝x
t-1
，bsx
t
＝x
t-s
，s表示延迟的时间长度；和表示差分算子，且d和d表示差分进行地次数；x
t
是时间序列在t时刻的观测值；θ(
·
)和θ(
·
)表示滑动平均多项式，且θ(z)＝1-θ1z-θ2z
2-...-θ
qzq
，θ(z)＝1-θ1z-θ2z
2-...-θ
qzq
，q和q限制滑动平均多项式地结构；ε
t
是在t
时刻的估计误差；z为变量。
[0018]
可选择地，所述步骤s5中，所述降噪自编码器的表达式为：
[0019]
d＝f(x)＝f(d+noise)
[0020]
其中，x表示大气廓线，d表示预测大气廓线，且x＝d+noise，其中noise表示噪音，f(x)表示映射函数。
[0021]
本发明具有以下有益效果：
[0022]
本发明所提供的大气温度廓线预测方法充分发挥了sarima模型和降噪自编码器(denoising autoencoder,dae)的优势。sarima模型小数据量的预测任务中有着明显的优势，它能充分挖掘数据之间的相关关系。但是sarima模型是线性模型，在处理线性相关关系较弱的平流层数据时，预测的精度较低，平流层大气温度数据的预测成为了整个大气温度模型的短板。本发明使用dae在垂直空间上对sarima模型的预测结果进行降噪优化。每个高度层的sarima模型的建模过程相互独立，在每个高度层上的预测误差也相互独立。在每个高度层上sarima模型的预测误差近似服从均值为0，方差为σ2(n)的正太分布，该误差分布只与高度层z有关。从垂直空间上来看，sarima模型的预测结果是受到噪声影响的大气温度廓线，噪声在每个高度层上稳定，服从近似的正态分布，通过引入降噪模型便可以对sarima模型进行优化。降噪自编码器又为模型引入了新的知识，使模型能够更好的对大气温度数据进行建模。根据实验，相比于单纯的sarima模型，sarima-dae混合模型在大气温度数据的预测中表现优异，特别是在线性相关关系较弱的平流层大气温度预测中预测精度得到了极大的提升。
附图说明
[0023]
图1为本发明所提供的大气温度廓线预测方法流程图；
[0024]
图2为图1中步骤s1的分步骤流程图；
[0025]
图3为本发明所提供的大气温度廓线预测方法的数据空间结构示意图；
[0026]
图4为不同高度层大气温度随时间的变化曲线图；
[0027]
图5为不同时间的大气温度廓线图。
具体实施方式
[0028]
以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
[0029]
实施例1
[0030]
本发明提供一种大气温度廓线预测方法，参考图1所示，所述大气温度廓线预测方法包括：
[0031]
s1：获取多个测量高度层；
[0032]
具体地，参考图2所示，所述步骤s1包括：
[0033]
s11：获取目标测量点的大气温度数据；
[0034]
s12：对所述目标测量点的大气温度数据进行编号，得到多个测量高度层。
[0035]
s2：分别获取多个所述测量高度层的中每层的测量数据；
[0036]
s3：将每层所述测量数据分别输入sarima模型中，得到与每层所述测量数据相对
应的各层预测数据；
[0037]
其中，sarima模型的具体结构为：
[0038][0039]
其中，φ(
·
)和φ(
·
)表示自回归多项式，且φ(z)＝1-φ1z-φ2z
2-...-φ
pzp
，φ(z)＝1-φ1z-φ2z
2-...-φ
pzp
；b和bs是延迟算子，且bx
t
＝x
t-1
，bsx
t
＝x
t-s
，s表示延迟的时间长度；和表示差分算子，且d和d表示差分进行地次数；x
t
是时间序列在t时刻的观测值；θ(
·
)和θ(
·
)表示滑动平均多项式，且θ(z)＝1-θ1z-θ2z
2-...-θ
qzq
，θ(z)＝1-θ1z-θ2z
2-...-θ
qzq
，q和q限制滑动平均多项式地结构；ε
t
是在t时刻的估计误差；z为变量。
[0040]
s4：根据所有的预测数据，得到一条大气廓线；
[0041]
s5：将所述大气廓线输入降噪自编码器(denoising autoencoder,dae)中，得到预测大气温度廓线。
[0042]
其中，所述降噪自编码器的表达式为：
[0043]
d＝f(x)＝f(d+noise)
[0044]
其中，x表示大气廓线，d表示预测大气廓线，且x＝d+noise，其中noise表示噪音，f(x)表示映射函数。
[0045]
实施例2
[0046]
实验所使用的数据是一组卫星探测的大气温度空间分布数据，该大气温度数据具有三维的空间分布结构，如图3所示，在空间上一共记录了96个高度层的大气温度数据。该大气温度数据是对中国广州省茂名市的一组卫星测量数据。
[0047]
在每个空间测量点上的大气温度数据可以从时间和空间两个视角来进行分析。从时间上来看，如图4所示，大气温度数据由一组时间序列数据组成，观测时间从2011年持续到2014年。其中，每个月的1号进行数据记录，每个空间点一共记录了48个时间点的观测数据。随着时间的变化而变化，大气温度有明显的季节性变化趋势，其中对流层和中间层的大气对流运动明显，大气温度的周期性变化更加明显，平流层的大气对流运动较弱，大气温度的变化的周期性特征较弱。从垂直空间上来看，如图5所示，大气温度数据由一组大气温度廓线组成。从图中可以明显的看出，大气温度廓线具有相似的空间结构，在不同高度层之间的大气温度数据具有明显的相关性。
[0048]
如图4所示，随着季节的更替，大气温度呈现出季节性变化趋势。大气温度时间序列相邻数据之间以及季节性相邻数据之间具有很强的相关性，使用sarima模型能够很好的对数据中的相关关系进行建模。sarima模型的具体结构为：
[0049][0050]
其中，φ(
·
)和φ(
·
)表示自回归多项式，且φ(z)＝1-φ1z-φ2z
2-...-φ
pzp
，φ(z)＝1-φ1z-φ2z
2-...-φ
pzp
；b和bs是延迟算子，且bx
t
＝x
t-1
，bsx
t
＝x
t-s
，s表示延迟的时间长度；和表示差分算子，且d和d表示差分进行地次数；x
t
是时间序列在t时刻的观测值；θ(
·
)和θ(
·
)表示滑动平均多项式，且θ(z)＝1-θ1z-θ2z
2-...-θ
qzq
，θ(z)＝1-θ1z-θ2z
2-...-θ
qzq
，q和q限制滑动平均多项式地结构；ε
t
是在t
时刻的估计误差；z为变量。
[0051]
降噪自编码器是自编码器的一种扩展，不同于自编码器用来完成自身到自身的映射，降噪自编码器用来进行噪声数据和干净数据之间的映射。理想的降噪自编码器为：
[0052]
d＝f(x)＝f(d+noise)
[0053]
其中，x表示大气廓线，d表示预测大气廓线，且x＝d+noise，其中noise表示噪音，f(x)表示映射函数。
[0054]
经过降噪映射，能够从受到噪声污染的数据x中重建干净的数据d。
[0055]
受到噪声污染的数据在自编码器的作用下重构成数据o，在进行重构的过程中先进行数据的编码，使数据从高维空间映射到低维空间，然后再从低维空间映射到高维空间。在降维的过程中维度较低的大气温度廓线数据的信息损失较小，而维度较高的噪声数据的损失较大，从而完成对数据的降噪。使用降噪后的数据o和纯净的数据d之间的误差用来训练降噪自编码器模型，在不断的迭代训练过程中使模型的误差不断地降低。
[0056][0057]
其中，d表示预测大气廓线；o表示降噪后的数据；di表示向量d中第i个元素；oi表示向量o中第i个元素。
[0058]
本发明具有以下有益效果：
[0059]
本发明所提供的大气温度廓线预测方法充分发挥了sarima模型和dae的优势。sarima模型小数据量的预测任务中有着明显的优势，它能充分挖掘数据之间的相关关系。但是sarima模型是线性模型，在处理线性相关关系较弱的平流层数据时，预测的精度较低，平流层大气温度数据的预测成为了整个大气温度模型的短板。本发明使用dae在垂直空间上对sarima模型的预测结果进行降噪优化。每个高度层的sarima模型的建模过程相互独立，在每个高度层上的预测误差也相互独立。在每个高度层上sarima模型的预测误差近似服从均值为0，方差为σ2(n)的正太分布，该误差分布只与高度层z有关。从垂直空间上来看，sarima模型的预测结果是受到噪声影响的大气温度廓线，噪声在每个高度层上稳定，服从近似的正态分布，通过引入降噪模型便可以对sarima模型进行优化。降噪自编码器又为模型引入了新的知识，使模型能够更好的对大气温度数据进行建模。根据实验，相比于单纯的sarima模型，sarima-dae混合模型在大气温度数据的预测中表现优异，特别是在线性相关关系较弱的平流层大气温度预测中预测精度得到了极大的提升。
[0060]
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

技术特征：
1.一种大气温度廓线预测方法，其特征在于，所述大气温度廓线预测方法包括：s1：获取多个测量高度层；s2：分别获取多个所述测量高度层中每层的测量数据；s3：将每层所述测量数据分别输入sarima模型中，得到与每层所述测量数据相对应的各层预测数据；s4：根据所有的预测层数据，得到一条大气廓线；s5：将所述大气廓线输入降噪自编码器中，得到预测大气温度廓线。2.根据权利要求1所述的多层大气温度预测方法，其特征在于，所述步骤s1包括：s11：获取目标测量点的历史大气温度廓线数据；s12：对测量点的大气温度廓线数据进行编号，得到多个测量高度层。3.根据权利要求1所述的多层大气温度预测方法，其特征在于，所述步骤s3中，所述sarima模型为：其中，φ(
·
)和φ(
·
)表示自回归多项式，且φ(z)＝1-φ1z-φ2z
2-...-φ
p
z
p
，φ(z)＝1-φ1z-φ2z
2-...-φ
p
z
p
；b和b
s
是延迟算子，且bx
t
＝x
t-1
，b
s
x
t
＝x
t-s
，s表示延迟的时间长度；和表示差分算子，且d和d表示差分进行地次数；x
t
是时间序列在t时刻的观测值；θ(
·
)和θ(
·
)表示滑动平均多项式，且θ(z)＝1-θ1z-θ2z
2-...-θ
q
z
q
，θ(z)＝1-θ1z-θ2z
2-...-θ
q
z
q
，q和q限制滑动平均多项式地结构；ε
t
是在t时刻的估计误差；z为变量。4.根据权利要求1所述的多层大气温度预测方法，其特征在于，所述步骤s5中，所述降噪自编码器的表达式为：d＝f(x)＝f(d+noise)其中，x表示大气廓线，d表示预测大气廓线，且x＝d+noise，其中noise表示噪音，f(x)表示映射函数。

技术总结
本发明公开了一种大气温度廓线预测方法，所述预测方法包括：S1：获取多个测量高度层；S2：分别获取多个所述测量高度层中每层的测量数据；S3：将每层所述测量数据分别输入SARIMA模型中，得到与每层所述测量数据相对应的各层预测数据；S4：根据所有的预测数据，得到一条大气廓线；S5：将所述大气廓线输入降噪自编码器中，得到预测大气温度廓线。得到预测大气温度廓线。


技术研发人员：郭兴 吴家骥 曹运华 朱松龄
受保护的技术使用者：西安电子科技大学
技术研发日：2021.12.06
技术公布日：2022/3/8