一种陆表温度数据同化处理方法、装置、设备及存储介质与流程

1.本技术涉及数据同化处理技术领域，具体而言，涉及一种陆表温度数据同化处理方法、装置、设备及存储介质。


背景技术：

2.准确的气候预测具有重要的社会经济价值，目前主要的气候预测方式是基于耦合气候模式的预测。耦合气候模式能合理模拟大气，海洋，陆地以及海冰等地球系统中各个子系统之间的相互作用过程，从而更全面模拟全球气候的演变过程。为与耦合模式相对应，传统的适用于单独陆面模式的陆表温度同化方法也需要发展到能适应多圈层耦合模式的耦合同化技术，以适应气候预测系统发展的需要。
3.陆表温度是全球天气和气候系统中一个重要变量，具有很强的日变化特性和空间差异性。现有的耦合同化系统中，针对陆表温度的耦合同化方案较少，已有的针对陆表温度变量的耦合同化算法也仅是针对时空变率不大的月尺度数据，时间跨度较大，预测精度较低。为了能够提高预测精度，需要使用时空变率更高的日尺度和小时尺度等高频陆表温度数据。但是现有的耦合同化系统不能够适应时空变率较大的观测资料。


技术实现要素：

4.有鉴于此，本技术的目的在于提供一种陆表温度数据同化处理方法、装置、设备及存储介质，以克服现有技术中的问题。
5.第一方面，本技术实施例提供了一种陆表温度数据同化处理方法，应用于耦合同化系统，所述耦合同化系统包括应用观测卫星数据，所述观测卫星数据提供待预测区域的陆表温度全球观测信息，并在对应的观测时刻生成观测数据；所述耦合同化系统中包含有与所述陆表温度具有关联关系的陆面其他相关关键变量数据；所述方法包括：基于各个格点上观测数据的温度均值数据和预设时间位相中包含的观测数据，计算得到该时间位相中各个格点的观测数据相较于所述温度均值数据的观测距平数据；针对所述观测距平数据进行傅式变换，得到所述观测距平数据对应的空间数据，并对所述空间数据进行随机扰动处理，得到处理后的随机集合样本；根据上个时刻的模式误差集合样本和当前时刻的随机集合样本，得到当前时刻三维空间的模式误差集合样本，并对当前时刻三维空间的模式误差集合样本进行傅式反变换，得到当前时刻二维空间的模式误差集合样本；根据当前时刻二维空间的模式数据的模式误差集合样本和所述陆面其他相关关键变量数据的模式误差集合样本之间的关联关系，计算得到预报时刻的集合样本，并针对预报时刻的集合样本使用目标集合kalman滤波同化算法处理，得到预报时刻的分析场。
6.在本技术一些技术方案中，上述方法通过以下方式确定预设时间位相中包含的观测数据，包括：根据卫星陆表温度产品时间频次特征，确定与所述时间频次特征对应的同化窗
口；所述同化窗口用于将观测值周期划分为多个时间位相；根据各个所述时间位相对应的时间数据和生成所述观测数据的观测时刻，将各个所述观测数据归整到各个所述时间位相中。
7.在本技术一些技术方案中，上述方法对所述空间数据进行随机扰动处理，得到处理后的随机集合样本，包括：根据所述空间数据，确定所述空间数据的误差方差；根据所述空间数据的误差方差和生成的随机数，得到所述随机集合样本。
8.在本技术一些技术方案中，上述目标集合kalman滤波同化算法包括原始集合kalman滤波同化算法和基于svd的背景场误差矩阵分解算法。
9.第二方面，本技术实施例提供了一种陆表温度数据同化处理装置，驻留于耦合同化系统，所述耦合同化系统包括应用观测卫星数据，所述观测卫星数据提供待预测区域的陆表温度全球观测信息，并在对应的观测时刻生成观测数据；所述耦合同化系统中包含有与所述陆表温度具有关联关系的陆面其他相关关键变量数据；所述装置包括：预处理模块，用于基于各个格点上观测数据的温度均值数据和预设时间位相中包含的观测数据，计算得到该时间位相中各个格点的观测数据相较于所述温度均值数据的观测距平数据；变换模块，用于针对所述观测距平数据进行傅式变换，得到所述观测距平数据对应的空间数据，并对所述空间数据进行随机扰动处理，得到处理后的随机集合样本；扰动模块，用于根据上个时刻的模式误差集合样本和当前时刻的随机集合样本，得到当前时刻三维空间的模式误差集合样本，并对当前时刻三维空间的模式误差集合样本进行傅式反变换，得到当前时刻二维空间的模式误差集合样本；同化模块，用于根据当前时刻二维空间的模式数据的模式误差集合样本和所述陆面其他相关关键变量数据的模式误差集合样本之间的关联关系，计算得到预报时刻的集合样本，并针对预报时刻的集合样本使用目标集合kalman滤波同化算法处理，得到预报时刻的分析场。
10.第三方面，本技术实施例提供了一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述的陆表温度数据同化处理方法的步骤。
11.第四方面，本技术实施例提供了一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器运行时执行上述的陆表温度数据同化处理方法的步骤。
12.本技术的实施例提供的技术方案可以包括以下有益效果：本技术方法包括基于各个格点上观测数据的温度均值数据和预设时间位相中包含的观测数据，计算得到该时间位相中各个格点的观测数据相较于所述温度均值数据的观测距平数据；针对所述观测距平数据进行傅式变换，得到所述观测距平数据对应的空间数据，并对所述空间数据进行随机扰动处理，得到处理后的随机集合样本；根据上个时刻的模式误差集合样本和当前时刻的随机集合样本，得到当前时刻三维空间的模式误差集合样本，并对当前时刻三维空间的模式误差集合样本进行傅式反变换，得到当前时刻二维空间的模式误差集合样本；根据当前时刻二维空间的模式数据的模式误差集合样本和所述陆面
其他相关关键变量数据的模式误差集合样本之间的关联关系，计算得到预报时刻的集合样本，并针对预报时刻的集合样本使用目标集合kalman滤波同化算法处理，得到预报时刻的分析场。本技术的方法针对日尺度以上频率陆表温度资料特点，在时间上通过日变化距平订正算法来满足日内尺度的陆表温度变量的时间变率要求，在空间上进行随机协调扰动，保证陆表温度空间过渡的连续性，缓解不同格点耦合同化不一致时产生的模式冲击，保持耦合同化过程中的稳定性和物理协调性。
13.为使本技术的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
14.为了更清楚地说明本技术实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍， 应当理解，以下附图仅示出了本技术的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
15.图1示出了本技术实施例所提供的一种陆表温度数据同化处理方法的流程示意图；图2示出了本技术实施例所提供的一种观测数据进行了划分归类的示意图；图3示出了本技术实施例所提供的一种陆表温度数据同化处理装置示意图；图4为本技术实施例提供的一种电子设备的结构示意图。
具体实施方式
16.为使本技术实施例的目的、技术方案和优点更加清楚，下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，应当理解，本技术中附图仅起到说明和描述的目的，并不用于限定本技术的保护范围。另外，应当理解，示意性的附图并未按实物比例绘制。本技术中使用的流程图示出了根据本技术的一些实施例实现的操作。 应该理解，流程图的操作可以不按顺序实现，没有逻辑的上下文关系的步骤可以反转顺序或者同时实施。 此外，本领域技术人员在本技术内容的指引下，可以向流程图添加一个或多个其他操作，也可以从流程图中移除一个或多个操作。
17.另外，所描述的实施例仅仅是本技术一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本技术实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本技术的实施例的详细描述并非旨在限制要求保护的本技术的范围，而是仅仅表示本技术的选定实施例。基于本技术的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本技术保护的范围。
18.需要说明的是，本技术实施例中将会用到术语“包括”，用于指出其后所声明的特征的存在，但并不排除增加其它的特征。
19.数值天气预报和气候预测研究中，由于全球范围内台站观测条件的限制，覆盖全球的陆面台站观测资料十分稀少，一般很难真正应用于气候预测业务。相对而言，卫星遥感资料则具有覆盖范围广、同步性好、时空分辨率高和观测要素全面等优点，可以较好弥补全球台站观测不足的缺点，因此卫星遥感资料的应用对数值预报和气候预测水平的提高具有
重要作用。
20.陆面温度是全球气候系统一个重要组成部分，是影响全球陆地水分和能量循环的重要因素，通过影响潜热和显热通量控制土壤—大气之间的能量和水分收支平衡，进而对全球天气和气候预测及水文预报产生显著影响。由于数值模式模拟的陆表温度在量值和时空分布上存在很大不确定性，而卫星遥感反演陆表温度又由于受扫描时间和地表粗糙度等参数限制，在量值上也存在较大误差，陆表温度同化则是实现卫星遥感观测和模式模拟最优融合从而获取准确陆表温度信息的一种有效手段。
21.准确的气候预测具有重要的社会经济价值，目前主要的气候预测方式为耦合气候模式预测。耦合气候模式能合理模拟大气，海洋，陆地以及海冰等地球系统中各个子系统之间的相互作用过程，从而更全面模拟全球气候的演变过程。为与耦合模式相对应，传统的适用于单独陆面模式的陆表温度同化方法也需要发展到能适应多圈层耦合模式的耦合同化技术，以适应气候预测系统发展的需要。
22.现有的针对陆面温度变量的耦合同化算法仅是针对时空变率不大的月尺度数据，时间跨度较大，预测精度较低。为了能够提高预测精度，需要使用时空变率更高的日尺度和小时尺度等高频陆表温度数据。但是现有的耦合同化系统不能够适应时空变率较大的观测资料。
23.基于此，本技术实施例提供了一种陆表温度数据同化处理方法、装置、设备及存储介质，下面通过实施例进行描述。
24.图1示出了本技术实施例所提供的一种陆表温度数据同化处理方法的流程示意图，其中，该方法应用于耦合同化系统，所述耦合同化系统包括应用观测卫星数据，所述观测卫星数据提供待预测区域的陆表温度全球观测信息，并在对应的观测时刻生成观测数据；所述耦合同化系统中包含有与所述陆表温度具有关联关系的陆面其他相关关键变量数据；该方法包括步骤s101-s104；具体的：s101、基于各个格点上观测数据的温度均值数据和预设时间位相中包含的观测数据，计算得到该时间位相中各个格点的观测数据相较于所述温度均值数据的观测距平数据；s102、针对所述观测距平数据进行傅式变换，得到所述观测距平数据对应的空间数据，并对所述空间数据进行随机扰动处理，得到处理后的随机集合样本；s103、根据上个时刻的模式误差集合样本和当前时刻的随机集合样本，得到当前时刻三维空间的模式误差集合样本，并对当前时刻三维空间的模式误差集合样本进行傅式反变换，得到当前时刻二维空间的模式误差集合样本；s104、根据当前时刻二维空间的模式数据的模式误差集合样本和所述陆面其他相关关键变量数据的模式误差集合样本之间的关联关系，计算得到预报时刻的集合样本，并针对预报时刻的集合样本使用目标集合kalman滤波同化算法处理，得到预报时刻的分析场。
25.本技术的方法针对日尺度以上频率陆表温度资料特点，在时间上通过日变化距平订正算法来满足日内尺度的陆表温度变量的时间变率要求，在空间上进行随机协调扰动，保证陆表温度空间过渡的连续性，缓解不同格点耦合同化不一致时产生的模式冲击，保持耦合同化过程中的稳定性和物理协调性。
26.下面对本技术的一些实施方式作详细说明。在不冲突的情况下，下述的实施例及实施例中的特征可以相互组合。
27.s101、基于各个格点上观测数据的温度均值数据和预设时间位相中包含的观测数据，计算得到该时间位相中各个格点的观测数据相较于所述温度均值数据的观测距平数据。
28.本技术中的陆表温度数据同化处理方法应用于耦合同化系统，耦合同化系统是指在考虑数据时空分布以及观测场和背景场误差的基础上,在多圈层全耦合数值预报模型的动态运行过程中融合新的观测数据的同化系统。它是在过程模型的动态框架内,通过数据同化算法不断融合时空上离散分布的不同来源和不同分辨率的直接或间接观测信息来自动调整模型轨迹,以改善动态模型状态的估计精度，提高模型预测能力，具体应用在天气预测等场景中。现有的耦合同化系统的对于输入数据的要求具有一定的限制，如果输入数据比较频繁（即通过更高的频率获取数据）的话，对于耦合同化系统的冲击较大，可以导致耦合同化系统不能够正常工作。一般来说，在具体实施的时候，现有技术中较多使用逐月数据，相较于逐月数据来说，逐日数据和逐小时数据就是高频率数据。相较于现有技术针对逐月数据处理，本技术可以理解为对于逐日数据和逐小时数据的处理。为了方便描述，后续本技术以逐日数据为例进行介绍本方案。本技术针对通过更高频率采集到的较多数据进行处理，以使耦合同化系统能够接受处理之后的数据。
29.本技术的耦合同化系统包括应用观测卫星数据，所述观测卫星数据提供待预测区域的陆表温度全球观测信息。这里待预测区域可以是某个地区也可以是整个地球。观测卫星在进行观测的时候，会在预设的一些观测时刻生成观测数据。例如，观测卫星对a城市进行观测，并在19.10、19.30、19.50分别获取了a城市的陆表温度20℃、19℃、18℃。为了避免直接使用观测数据给耦合同化系统带来的影响，本技术对于各个格点上得到的观测数据进行了日距平数据处理。
30.本技术中观测数据是遥感卫星分别从对应的格点上观测得到的，本技术以格点为依据对各个观测数据进行处理。统计记得到各个格点上观测得到的观测数据的个数和各个观测数据的具体温度数值，然后根据该格点上观测数据的个数和各个观测数据的具体温度数值，计算得到在该格点上观测得到的所有观测数据的温度均值数据。例如，在某个格点上进行了三次观测，分别得到28℃、28.1℃和28.2℃，则该格点的温度均值数据就是28.1℃。
31.需要说明的是，卫星进行观测时候，由于某些原因，可能会存在未采集到观测数据的问题，导致某些格点上会有观测数据的缺测。在计算温度均值数据的时候，本技术仅考虑有值的格点。
32.在另一方面，为了避免较多数据的处理，本技术将观测数据进行了划分归类。具体的如图2所示包括步骤s201~s202：s201、根据卫星陆表温度产品时间频次特征，确定与所述时间频次特征对应的同化窗口；所述同化窗口用于将观测值周期划分为多个时间位相；s202、根据各个所述时间位相对应的时间数据和生成所述观测数据的观测时刻，将各个所述观测数据归整到各个所述时间位相中。
33.本技术考虑了卫星陆表温度产品时间频次特征，然后基于卫星陆表温度产品时间频次特征，确定与所述时间频次特征对应的同化窗口。在确定了同化窗口之后，本技术通过
同化窗口将观测周期划分为多个时间位相。例如，同化窗口为一小时，如果观测时间周期为二十四小时的话，那么就可以得到二十四个时间位相分别为：一点到两点、两点到三点
……
十一点到十二点。然后本技术中各个观测数据都有对应的观测时刻，通过观测时刻与各个时间位相进行对比，确定各个观测时刻被哪个时间位相所包含，就将与该观测时刻对应的观测数据归在该时间位相下。
34.在得了各个格点上观测数据的温度均值数据和预设时间位相中包含的观测数据之后，通过二者作差即可得到该时间位相中各个格点的观测数据相较于所述温度均值数据的观测距平数据。本技术基于观测距平数据对耦合模式中的数据进行订正，降低了同化对耦合模式带来的冲击和导致的不协调性。
35.在具体实施的时候，本技术是基于集合kalman滤波同化算法的框架实现的。在考虑陆表温度日变化距平信息后，在集合kalman滤波同化算法中，各时步中的样本时间演变预报方程可表示为：其中，为模式积分时步，为集合扰动样本，n为集合扰动样本总数，x为陆表温度的日均值距平，上标“+”和
“‑”
分别表示分析场和背景场，表示耦合模式中包括陆表温度距平变量的非线性积分算子。各时步中状态变量的分析更新方程可表示为：其中，为陆表温度日距平观测，为实现陆表温度日距平观测向背景场转换的观测算子，为kalman增益矩阵，可表示为：其中，上标
ꢀ“‑”
表示背景场预报时步， t为矩阵转置，为陆表温度日距平观测误差协方差矩阵。其中b矩阵通过集合样本统计得到：误差协方差矩阵。其中b矩阵通过集合样本统计得到：为陆表温度日变化距平信息的背景场误差协方差矩阵，表示陆表温度日距平变量的集合背景场误差向量，其中为陆表温度日距平变量的集合样本均值。
36.在进行同化分析时，为避免直接更新陆表温度变量可能会对模式状态造成的过大冲击而导致不稳定现象，此时的耦合同化中的更新变量为陆表温度日均值距平信息。通过陆表温度观测信息对耦合模式中陆表温度日均值距平信息的调整，可以更好地保持耦合模式陆表温度变量的连续性，同时又合理引入观测的陆表温度日距平变化信息，从而既在模式中充分保留了卫星遥感陆表温度的日变化信息，又合理保证了耦合模式连续积分的稳定性。
37.s102、针对所述观测距平数据进行傅式变换，得到所述观测距平数据对应的空间数据，并对所述空间数据进行随机扰动处理，得到处理后的随机集合样本。
38.所述对所述空间数据进行随机扰动处理，得到处理后的随机集合样本，包括：根据
所述空间数据，确定所述空间数据的误差方差；根据所述空间数据的误差方差和生成的随机数，得到所述随机集合样本。
39.s103、根据上个时刻的模式误差集合样本和当前时刻的随机集合样本，得到当前时刻三维空间的模式误差集合样本，并对当前时刻三维空间的模式误差集合样本进行傅式反变换，得到当前时刻二维空间的模式误差集合样本。
40.本技术在得到观测距平数据之后，同样的对于耦合模式而言，也需要确定出耦合模式中模式数据的耦合距平信息。具体的，本技术的耦合同化系统包括耦合距平分析子模块，通过该耦合距平分析子模块确定出模式距平数据。例如，耦合距平分析子模块，根据确认的卫星观测陆表温度日变化位相信息和日变化距平信息，得到各个位相上观测记录相比日均值的距平信息。同样，对于耦合模式而言，在其中加入陆表温度日变化位相和距平分析子模块，通过该子模块实现在对应同化窗口中给出耦合模式中陆表温度的日变化位相和日变化信息。然后，根据确认的卫星观测陆表温度日变化位相信息和日变化距平信息，得到各个位相上观测记录相比日均值的距平信息。同样，耦合模式中的距平分析子模块也会在同化时刻对应生成模式陆表温度的日均值距平信息。在进行同化分析时，为避免直接更新陆表温度变量可能会对模式状态造成的过大冲击而导致不稳定现象，此时的耦合同化中的更新变量为陆表温度日均值距平信息。通过陆表温度观测信息对耦合模式中陆表温度日均值距平信息的调整，可以更好地保持耦合模式陆表温度变量的连续性，同时又合理引入观测的陆表温度日距平变化信息，从而既在模式中充分保留了卫星遥感陆表温度的日变化信息，又合理保证了耦合模式连续积分的稳定性。
41.由于卫星遥感陆表温度观测资料在空间上存在很多缺测，如果直接使用传统的集合kalman滤波同化算法，在水平空间对陆表温度变量进行直接更新，由于相邻格点上因为缺测状态的不同而出现并不都进行同化分析更新的情况，从而导致同化后的陆表温度变量出现空间上的不连续，这一不连续状态经过耦合模式积分的放大效应极易导致误差增长，从而产生很大不确定性并产生溢出导致耦合模式运行中断，引起同化失败。为了避免这一现象，本技术进一步的对观测距平数据进行了随机扰动处理。
42.也就是说，本技术并没有直接使用观测距平数据，而是对观测距平数据进行调整。在对观测距平数据进行调整时候，本技术没有直接对观测距平数据进行操作，而是对观测距平数据对应的空间数据进行处理，能够进一步的降低对于耦合模式的影响。具体的，本技术通过正向傅式变换，得到观测距平数据对应的空间数据，之后，本技术对空间数据进行了扰动处理。另一方面，集合kalman滤波同化算法自身由于使用集合样本求解背景场误差协方差矩阵，会由于模式误差项的缺失在集合同化运行过程中导致滤波发散现象。为避免这一现象，本发明在集合kalman滤波同化算法中引入基于fourier变换的二维monte carlo随机扰动算法，一方面可以更好地实现水平方向上陆表温度日距平模式误差的二维协调扰动，另一方面也可以克服集合kalman滤波同化方法中由于集合样本求背景场误差协方差矩阵而缺失的模式误差项。该算法步骤如下：针对水平二维空间上n
×
m网格内的陆表温度日距平变量场，对其进行离散形式的fourier变换：
其中， l和p为计数，和分别为n和m方向的波数。假设的形式如下：其中是一组用来引入随机位相平移的随机数。此时，变量场在空间上的协方差通过变换可表示为：通过采用该算法，在集合同化中叠加到集合样本上的随机扰动场是作用于fourier空间而不是直接作用于前面得到的陆表温度日距平变量，能有效避免直接在空间不同点的陆表温度日距平上加随机扰动导致的空间各网格点扰动间的不协调现象，并保证各模式格点引入的陆表温度距平模式误差在随模式积分时间变化时也同时保持良好的空间协调性，从而保证陆表温度集合预报样本能综合体现陆表温度空间变化的协调相关关系。
43.此时，在考虑了水平方向上陆表温度日距平变量的协调扰动后，在每个集合kalman滤波同化时步中的集合样本时间演变预报方程可进一步改写为：其中就为二维空间场上添加的陆表温度日距平模式误差项，用来克服集合kalman滤波同化方法中由于集合样本求背景场误差协方差矩阵而缺失的模式误差项。中体现的是陆表温度日距平模式误差（进而也是保证陆表温度变量误差自身）能保持二维空间上的协调关系，从而保证了整个空间场上高频陆表温度观测资料耦合同化的效果。
44.s104、根据当前时刻二维空间的模式数据的模式误差集合样本和所述陆面其他相关关键变量数据的模式误差集合样本之间的关联关系，计算得到预报时刻的集合样本，并针对预报时刻的集合样本使用目标集合kalman滤波同化算法处理，得到预报时刻的分析场。
45.本技术针对卫星遥感陆表温度资料的耦合同化整体是基于集合kalman滤波同化框架。耦合同化时状态变量的更新方程前面已给出。集合kalman滤波同化的另一个优势就是每个时步能对背景场误差协方差矩阵进行实时更新，从而及时调整整体误差结构而不断优化同化效果。集合kalman滤波同化的主要计算量也来自对背景场误差协方差矩阵的计算和更新。为耦合同化时对背景场误差协方差矩阵求解时的计算量，并保持其稳定性和协调性，在其中引入了背景场误差协方差矩阵svd分解算法，该算法主要针对求解耦合同化分析矩阵时涉及到的背景场误差协方差矩阵的求逆计算量过大的问题：
其中，上标“+”和
“‑”
分别表示分析场和背景场，为观测误差向量，=为陆表温度日距平变量的集合观测误差向量，表示陆表温度日距平变量的集合背景场误差向量，为预报时刻陆表温度日变化距平的背景场误差协方差矩阵，为分析时刻陆表温度日变化距平的背景场误差协方差矩阵。上式为陆表温度背景场误差矩阵的更新方程，其中括号中的背景场误差协方差矩阵求逆运算往往计算量很大，尤其是当观测向量的维度m很大时，其求解计算量正比于m2，计算量呈几何级数增加。为减小耦合同化时的计算量并保持计算精度，引入一种基于svd分解的误差协方差求逆算法，针对上式括号中的求逆运算，采用一种基于svd的背景场误差矩阵分解算法的计算策略：其中，的维度为m
×
n，n为集合样本数，相比m2有明显减小。此时对进行svd分解运算，得到，然后就有：其中n
×
n的上对角矩阵，包含的是svd分解后的非零特征值，而向量中则储存着svd分解后的特征向量，并且主要的特征信息都储存在中的前n个特征向量中。通过这种svd变换策略，可以将背景场误差协方差的求逆运算的计算量从m
×
m量级减小到n
×
n量级，从而大大减小了背景场误差协方差矩阵求逆运算量，并增加了背景场集合样本的代表性和稳定性，从而能稳定实现在耦合模式中同化时陆表温度变量的耦合同化分析更新。
46.为进一步增加耦合同化的稳定性，在集合kalman滤波同化框架中，为避免直接在陆表温度误差距平变量集合预报样本上加随机扰动后可能对陆表温度的陆面其他相关关键变量数据（土壤湿度、土壤温度、积雪）等变量之间的物理协调性的破坏，采用多参数平衡扰动算法改进了集合kalman滤波同化方案，该方案首先通过耦合模式集合样本统计陆表温度与土壤湿度、土壤温度、积雪等关键陆表变量之间的物理约束关系，然后将这些建立的物理约束关系作为构建背景场误差矩阵的物理约束条件，从而保证构建耦合模式集合误差时，能保持陆表温度与其他关键陆面变量的随机扰动之间的平衡关系，在同化时保留背景场集合预报样本中关键陆面变量之间的物理协调性，增加耦合同化系统的稳定性，并改善耦合同化效果。
47.我们将陆表温度距平与土壤湿度距平、土壤温度距平、积雪距平之间的物理约束关系分别定义为，，，这些物理约束关系都是通过耦合模式内部的非线性关系计算得到，因此在同化中的每个集合样本中都是能得到保持的。在考虑陆表温度日距平与其他关键陆面变量之间约束关系后，在每个集合kalman滤波同化时步中的集合样本预报方程可扩展为：其中，为模式积分时步，为集合扰动样本，n为集合扰动样本总数，x为陆表温度距平，上标“+”和
“‑”
分别表示分析场和背景场，，，分别为陆
表温度距平与土壤湿度距平、土壤温度距平、积雪距平之间的物理约束关系，表示耦合模式中包括陆表温度距平、土壤湿度距平及其约束关系、土壤温度距平及其约束关系、积雪距平及其约束关系的非线性积分算子。此时集合kalman滤波同化中背景场误差协方差矩阵可表示为：其中，由集合样本统计得到的背景场误差向量中增加了包含陆表温度距平与其他变量距平之间约束关系而引入的模式误差部分，其中，为土壤湿度距平约束关系集合样本均值，为土壤温度距平约束关系集合样本均值，为积雪距平约束关系集合样本均值。此时同化中的kalman增益矩阵可扩展为：通过观测时刻将陆表温度距平和其他变量距平的物理约束关系作用于增益矩阵和背景场误差协方差矩阵，一方面对集合预报样本统计缺失的模式误差项进行补偿，另一方面对背景场误差协方差矩阵中缺失的陆表温度与其他关键陆面变量的物理协调关系进行补偿。
48.本技术提供了一种完全适用于高频次（日尺度和小时尺度）卫星遥感陆表温度产品在耦合模式中的协调耦合同化算法，能有效地提升数值预报模式和气候模式系统对包括我国自主风云卫星在内的国内外多种气象卫星陆表温度产品的耦合同化能力。该方法还真正实现了耦合气候模式对卫星遥感陆表温度信息的协调同化应用，从而使得卫星遥感陆表温度产品信息能真正应用于全耦合数值预报和气候模式，并通过耦合器对整个耦合模式各分量产生物理协调的影响，从而能更好地为进行数值天气预报和气候预测提供更准确和更协调的全球陆表状态信息，可以为现有耦合模式和耦合预测系统提供高质量的、准确性和精度都更高、物理协调的陆面初值，更好地发挥多种卫星遥感陆表温度产品在耦合模式系统中的优势和价值。
49.本技术可用于制作各行业所需的气象陆面数据。例如，国家防总需要未来数月内全国干旱、高温热浪、寒潮、农业墒情、农作物产量等预测和规划，这些都可以用高频次的卫星陆表温度资料耦合同化技术提供精细化的合理初值来进行预测。
50.图3示出了本技术实施例所提供的一种陆表温度数据同化处理装置的结构示意图，驻留于耦合同化系统，驻留于耦合同化系统，所述耦合同化系统包括应用观测卫星数据，所述观测卫星数据提供待预测区域的陆表温度全球观测信息，并在对应的观测时刻生
成观测数据；所述耦合同化系统中包含有与所述陆表温度具有关联关系的陆面其他相关关键变量数据；所述装置包括：预处理模块，用于基于各个格点上观测数据的温度均值数据和预设时间位相中包含的观测数据，计算得到该时间位相中各个格点的观测数据相较于所述温度均值数据的观测距平数据；变换模块，用于针对所述观测距平数据进行傅式变换，得到所述观测距平数据对应的空间数据，并对所述空间数据进行随机扰动处理，得到处理后的随机集合样本；扰动模块，用于根据上个时刻的模式误差集合样本和当前时刻的随机集合样本，得到当前时刻三维空间的模式误差集合样本，并对当前时刻三维空间的模式误差集合样本进行傅式反变换，得到当前时刻二维空间的模式误差集合样本；同化模块，用于根据当前时刻二维空间的模式数据的模式误差集合样本和所述陆面其他相关关键变量数据的模式误差集合样本之间的关联关系，计算得到预报时刻的集合样本，并针对预报时刻的集合样本使用目标集合kalman滤波同化算法处理，得到预报时刻的分析场。
51.通过以下方式确定预设时间位相中包含的观测数据，包括：根据卫星陆表温度产品时间频次特征，确定与所述时间频次特征对应的同化窗口；所述同化窗口用于将观测值周期划分为多个时间位相；根据各个所述时间位相对应的时间数据和生成所述观测数据的观测时刻，将各个所述观测数据归整到各个所述时间位相中。
52.所述对所述空间数据进行随机扰动处理，得到处理后的随机集合样本，包括：根据所述空间数据，确定所述空间数据的误差方差；根据所述空间数据的误差方差和生成的随机数，得到所述随机集合样本。
53.所述目标集合kalman滤波同化算法包括原始集合kalman滤波同化算法和基于svd的背景场误差矩阵分解算法。
54.如图4所示，本技术实施例提供了一种电子设备，用于执行本技术中的陆表温度数据同化处理方法，该设备包括存储器、处理器、总线及存储在该存储器上并可在该处理器上运行的计算机程序，其中，上述处理器执行上述计算机程序时实现上述的陆表温度数据同化处理方法的步骤。
55.具体地，上述存储器和处理器可以为通用的存储器和处理器，这里不做具体限定，当处理器运行存储器存储的计算机程序时，能够执行上述的陆表温度数据同化处理方法。
56.对应于本技术中的陆表温度数据同化处理方法，本技术实施例还提供了一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，该计算机程序被处理器运行时执行上述的陆表温度数据同化处理方法的步骤。
57.具体地，该存储介质能够为通用的存储介质，如移动磁盘、硬盘等，该存储介质上的计算机程序被运行时，能够执行上述的陆表温度数据同化处理方法。
58.在本技术所提供的实施例中，应该理解到，所揭露系统和方法，可以通过其它的方式实现。以上所描述的系统实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间
的耦合或直接耦合或通信连接可以是通过一些通信接口，系统或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
59.所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
60.另外，在本技术提供的实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
61.所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本技术的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行本技术各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器（read-only memory， rom）、随机存取存储器（random access memory ，ram）、磁碟或者光盘等各种可以存储程序代码的介质。
62.应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释，此外，术语“第一”、“第二”、“第三”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
63.最后应说明的是：以上所述实施例，仅为本技术的具体实施方式，用以说明本技术的技术方案，而非对其限制，本技术的保护范围并不局限于此，尽管参照前述实施例对本技术进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本技术揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本技术实施例技术方案的精神和范围。都应涵盖在本技术的保护范围之内。因此，本技术的保护范围应所述以权利要求的保护范围为准。