CN111275250B - 一种基于各向异性的强流区海表面温度预报方法 - Google Patents
一种基于各向异性的强流区海表面温度预报方法 Download PDFInfo
- Publication number
- CN111275250B CN111275250B CN202010044424.3A CN202010044424A CN111275250B CN 111275250 B CN111275250 B CN 111275250B CN 202010044424 A CN202010044424 A CN 202010044424A CN 111275250 B CN111275250 B CN 111275250B
- Authority
- CN
- China
- Prior art keywords
- omega
- diffusion
- model
- surface temperature
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000013277 forecasting method Methods 0.000 title claims description 7
- 238000009792 diffusion process Methods 0.000 claims abstract description 43
- 238000000034 method Methods 0.000 claims abstract description 36
- 239000004576 sand Substances 0.000 claims abstract description 3
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000011160 research Methods 0.000 claims description 4
- 239000012530 fluid Substances 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 238000001914 filtration Methods 0.000 abstract description 15
- 230000003292 diminished effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于各向异性的强流区海表面温度预报方法:(1)给定控制变量ω(x,y)的初始猜测值;(2)设T为研究变量,从“时间”t=0正向积分平流扩散模型到t=S得到TS;并运用参数化方法得到扩散系数a;(3)积分切线性伴随模型得到目标函数J(ω)的梯度g(ω);(4)最小化算法优化控制变量ω(x,y);(5)循环步骤(2)至步骤(4)直至控制变量ω(x,y)使目标函数J(ω)达最小值,输出此时的TS;(6)根据步骤(5)得到的TS建立分析场,所建立的分析场为强流区海表面温度预报提供初始场。本发明针对原始扩散滤波各向同性的特点,采用对扩散方程加平流项并参数化扩散系数的方法,将海流的影响加入到整个同化过程中使同化结果更加贴合真实值。
Description
技术领域
本发明涉及海洋环境工程,特别涉及一种基于各向异性的强流区海表面温度预报方法。
背景技术
海洋数值预报是指根据海洋实际情况,在一定的初值和边值条件下,通过大型计算机作数值计算,求解描写海洋演变过程的流体力学和热力学的方程组,预测未来一定时段的海洋运动状态和海洋现象的方法,即利用当前海洋状况作为输入数据而作出海洋预报的手段。如何得到更为准确的预报结果最重要的一方面就是初值问题,目前来说采用数据同化方法更够得到更为准确的初始场从而提高预报准确率,通常标准的三维变分数据同化方法被用来为海表面温度数值预报提供初始场。
一般来说,标准的三维变分数据同化方法能被写为下面目标函数最小化问题:
其中x是分析向量,xb是背景场,y是观测向量,H是从模型点插值到观测点的插值算子,R是观测误差协方差矩阵,(.)T是转置,(.)-1指的是逆矩阵。B是背景误差协方差矩阵,B矩阵直接决定观测值订正到背景场的程度.如果矩阵被高估了,背景项在目标函数中的作用就会变小,这样使得最小化主要针对观测项,而且会过分拟合观测场,忽略较好的背景信息。如果矩阵被低估了,观测项在目标函数中的作用就会变小,这使得最小化主要针对背景项.同样,如果过分拟合背景场,质量好的观测场的作用就会减弱.因此,变分同化方法中合理估计背景误差协方差对改善同化效果具有重要意义。
由于实际的海洋环境复杂并且真实情况未知,所以在任何数据同化方法中包括三维变分背景误差协方差矩阵B的确定都是一个挑战,如何解决这个问题成为数据同化研究者们的重点科研内容,目前有两种方法来确定B.第一种即是相关尺度法,但此方法无法求B-1并且需要大的存储空间来存储,具有一定的局限性。第二种方法是引进一个新的变量ω,被定义为:
ω=B-1(x-xb) (2)
用递归滤波来描述B:
Yi=αYi-1+(1-α)Xi (3)
Zi=αZi+1+(1-α)Yi (4)
其中Xi是在格点i处的初值,Yi是在从i=1到n滤波之后的值,Zi是在每个方向上滤波一遍之后的初值,α是滤波系数,它决定了观测信息传播到整个分析域的速度。对于多重滤波只需要一直重复(2)(3)的操作即可。多维滤波可以通过在每个方向上执行一维滤波来构建。
矩阵B被视为高斯型在递归滤波中可以被视为高斯滤波过程,而高斯滤波就是扩散滤波。通过求解扩散方程并结合最小化算法可以得到最后的分析场。但此过程扩散系数的确定经验性比较强,并且这种方法具有各向同性的性质,在强流区同化结果较不理想,进而导致在强流区的数值预报中预报结果较差。
发明内容
本发明的目的是克服现有技术中的不足,提供一种基于各向异性的强流区海表面温度预报方法,本发明针对原始扩散滤波各向同性的特点,采用对扩散方程加平流项并参数化扩散系数的方法,将海流的影响加入到整个同化过程中使同化结果更加贴合真实值。
本发明所采用的技术方案是:一种基于各向异性的强流区海表面温度预报方法,包括以下步骤:
步骤1,给定控制变量ω(x,y)的初始猜测值;
步骤2,设T为研究变量,从“时间”t=0正向积分平流扩散模型到t=S得到TS;并运用参数化方法得到扩散系数a;
步骤3,积分切线性伴随模型得到目标函数J(ω)的梯度g(ω);
步骤4,最小化算法优化控制变量ω(x,y);
步骤5,循环步骤2至步骤4直至控制变量ω(x,y)使目标函数J(ω)达最小值,输出此时的TS;
步骤6,根据步骤5得到的TS建立分析场,所建立的分析场为强流区海表面温度预报提供初始场。
进一步地,步骤2中,所述的平流扩散模型为:
进一步地,步骤2中,所述的运用参数化方法得到扩散系数a包括:
采用TS代表T|t=S,则目标函数J(ω)表示为:
式中,H表示从模型点插值到观测点的插值算子,R表示观测误差协方差矩阵,d是观测增量;
运用Smagorinsky扩散公式(6)参数化扩散系数a:
进一步地,步骤3中,所述的积分切线性伴随模型得到目标函数J(ω)的梯度g(ω)包括:
采用code2code形式,调用tapenade编译器自动生成切线性伴随模型代码;
从“时间”t=S到t=0积分切线性伴随模型得到切线性伴随模型在t=0时刻的解R0(ω);
目标函数J(ω)的梯度g(ω)即为-R0(ω)。
本发明的有益效果是:本发明一种基于各向异性的强流区海表面温度预报方法,采用基于传统扩散滤波构建流依赖的三维变分数据同化方法,用于为数值预报提供初始场。本发明将扩散方程加上平流项并按照同样的方法求解平流扩散方程,同时将扩散系数参数化来代替传统方法中的常数。相比于传统的扩散滤波数据同化方法,该方法体现出各向异性的特点,在同化过程中考虑了海流信息的影响,特别是在例如黑潮这种强流区,能够得到比传统方法误差更小的分析场,为数值预报的准确性提供了一定的保证。
附图说明
图1:本发明一种基于各向异性的强流区海表面温度预报方法流程示意图。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效,兹例举以下实施例,并配合附图详细说明如下:
传统扩散滤波扩散系数一般定为常数,存在各向同性的缺点,在流场比较强的海区同化结果误差较大。本发明在同化过程中加入了海流信息的影响,使用平流扩散方程进行滤波,并且使用Smagorinsky方程将扩散系数参数化,使得在强流区本发明方法比传统方法可以得到更加贴近真实场的分析场。
如附图1所示,一种基于各向异性的强流区海表面温度预报方法,采用基于流依赖的三维变分数据同化方法,通过加入海流信息来改进传统的扩散滤波数据同化方法是本发明的重要内容。具体包括以下步骤:
步骤1,给定控制变量ω(x,y)的初始猜测值,一般设为ω(x,y)=0。
步骤2,设T为研究变量,从“时间”t=0正向积分平流扩散模型到t=S得到TS:
式中,t表示时间;u表示东西下的流速;v表示南北下的流速;T(ω)表示研究变量;a表示扩散系数;S表示积分时长;D表示研究区域的宽度,L表示研究区域的长度;Ω是的中心域,Γ是的边界;n是Γ的外法向量。
如果用TS代表T|t=S,目标函数J(ω)则可以表示为:
式中,H表示从模型点插值到观测点的插值算子,R表示观测误差协方差矩阵,d是观测增量。
运用Smagorinsky扩散公式(6)参数化扩散系数a:
步骤3,积分切线性伴随模型得到目标函数J(ω)的梯度g(ω)。
步骤3.1,切线性伴随模型的获得:
由于扩散方程的引入,三维变分就变为特殊的四维变分,这里使用四维变分里常用的伴随法求目标函数的梯度,采用code2code形式,调用tapenade编译器自动生成切线性伴随模型代码。
步骤3.2,从“时间”t=S到t=0积分切线性伴随模型得到切线性伴随模型在t=0时刻的解R0(ω);目标函数J(ω)的梯度g(ω)即为-R0(ω)。
步骤4,最小化算法,即下降算法优化控制变量ω(x,y)。
步骤5,循环步骤2至步骤4直至控制变量ω(x,y)使目标函数J(ω)达最小值,输出此时的TS。
步骤6,根据步骤5得到的TS,则x=TS+xb即为我们所需要的分析场,其中xb为背景场,此分析场为强流区海表面温度预报提供初始场。
尽管上面结合附图对本发明的优选实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,并不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可以做出很多形式,这些均属于本发明的保护范围之内。
Claims (1)
1.一种基于各向异性的强流区海表面温度预报方法,其特征在于,包括以下步骤:
步骤1,给定控制变量ω(x,y)的初始猜测值;
步骤2,设T为研究变量,从“时间”t=0正向积分平流扩散模型到t=S得到TS;并运用参数化方法得到扩散系数a;
其中,所述的平流扩散模型为:
其中,所述的运用参数化方法得到扩散系数a包括:
采用TS代表T|t=S,则目标函数J(ω)表示为:
式中,H表示从模型点插值到观测点的插值算子,R表示观测误差协方差矩阵,d是观测增量;
运用Smagorinsky扩散公式(6)参数化扩散系数a:
步骤3,积分切线性伴随模型得到目标函数J(ω)的梯度g(ω),包括:
采用code2code形式,调用tapenade编译器自动生成切线性伴随模型代码;
从“时间”t=S到t=0积分切线性伴随模型得到切线性伴随模型在t=0时刻的解R0(ω);
目标函数J(ω)的梯度g(ω)即为-R0(ω);
步骤4,最小化算法优化控制变量ω(x,y);
步骤5,循环步骤2至步骤4直至控制变量ω(x,y)使目标函数J(ω)达最小值,输出此时的TS;
步骤6,根据步骤5得到的TS建立分析场,所建立的分析场为强流区海表面温度预报提供初始场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010044424.3A CN111275250B (zh) | 2020-01-15 | 2020-01-15 | 一种基于各向异性的强流区海表面温度预报方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010044424.3A CN111275250B (zh) | 2020-01-15 | 2020-01-15 | 一种基于各向异性的强流区海表面温度预报方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111275250A CN111275250A (zh) | 2020-06-12 |
CN111275250B true CN111275250B (zh) | 2022-05-17 |
Family
ID=71003132
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010044424.3A Expired - Fee Related CN111275250B (zh) | 2020-01-15 | 2020-01-15 | 一种基于各向异性的强流区海表面温度预报方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111275250B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111859249B (zh) * | 2020-06-08 | 2022-06-14 | 天津大学 | 一种基于解析四维集合变分的海洋数值预报方法 |
CN113945997B (zh) * | 2021-09-26 | 2023-10-10 | 天津大学 | 一种基于解析四维集合变分的提高海洋预报精度的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104412875A (zh) * | 2013-08-22 | 2015-03-18 | 杨琳 | 循环资源与循环能源的集成装置 |
CN105205484A (zh) * | 2014-11-26 | 2015-12-30 | 中国人民解放军第二炮兵工程大学 | 基于曲波变换与维纳滤波的合成孔径雷达目标检测方法 |
CN105981025A (zh) * | 2014-12-01 | 2016-09-28 | 哈尔滨工程大学 | 基于序列递归滤波三维变分的实测海洋环境数据同化方法 |
CN107247290A (zh) * | 2017-07-06 | 2017-10-13 | 西安交通大学 | 一种基于时空分数阶滤波的地震资料噪声压制方法 |
CN108469620A (zh) * | 2018-03-26 | 2018-08-31 | 江苏省有色金属华东地质勘查局地球化学勘查与海洋地质调查研究院 | 适用于辐射沙脊群浅水海域的水下地形测量方法 |
CN109754058A (zh) * | 2018-05-08 | 2019-05-14 | 天津大学 | 一种基于cgbp算法的深度基准面拟合方法 |
-
2020
- 2020-01-15 CN CN202010044424.3A patent/CN111275250B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104412875A (zh) * | 2013-08-22 | 2015-03-18 | 杨琳 | 循环资源与循环能源的集成装置 |
CN105205484A (zh) * | 2014-11-26 | 2015-12-30 | 中国人民解放军第二炮兵工程大学 | 基于曲波变换与维纳滤波的合成孔径雷达目标检测方法 |
CN105981025A (zh) * | 2014-12-01 | 2016-09-28 | 哈尔滨工程大学 | 基于序列递归滤波三维变分的实测海洋环境数据同化方法 |
CN107247290A (zh) * | 2017-07-06 | 2017-10-13 | 西安交通大学 | 一种基于时空分数阶滤波的地震资料噪声压制方法 |
CN108469620A (zh) * | 2018-03-26 | 2018-08-31 | 江苏省有色金属华东地质勘查局地球化学勘查与海洋地质调查研究院 | 适用于辐射沙脊群浅水海域的水下地形测量方法 |
CN109754058A (zh) * | 2018-05-08 | 2019-05-14 | 天津大学 | 一种基于cgbp算法的深度基准面拟合方法 |
Non-Patent Citations (10)
Title |
---|
A New Version of Regional Ocean Reanalysis for Coastal Waters of China and Adjacent Seas;Guijun Han;《ADVANCES IN ATMOSPHERIC SCIENCES》;20131231;全文 * |
Disjoint hypercyclic powers of weighted translations on groups;ZHANG Liang;《Czechoslovak Mathematical Journal》;20170810;全文 * |
Ensemble estimates of the wave state related parameters in a sea speay parameterizaion;Zhang Lianxin;《Springer Link》;20190508;全文 * |
Estimating Convection Parameters in the GFDL CM2.1 Model Using Ensemble Data Assimilation;Shan Li;《Journal of Advances in Modeling Earth Systems》;20180325;全文 * |
Water depth variation of southern Bohai Sea in recent 20 years;Yue Baojing;《Marine Geology & Quatemary Geology》;20100630;全文 * |
基于扩散滤波的多尺度三维变分研究;李冬;《海洋通报》;20110415;全文 * |
岛群河口开发利用的水环境效应研究;许婷;《中国优秀博硕士学位论文全文数据库(博士) 工程科技I辑》;20170501;全文 * |
水下热目标红外可探测性研究;郝果清;《数字通信世界》;20181201;全文 * |
海底地貌复杂度的表示与计算;翟京生;《中国海洋大学学报(自然科学版)》;20190815;全文 * |
集合滤波和三维变分混合数据同化方法研究术;吴新荣;《热带海洋学报》;20111115;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111275250A (zh) | 2020-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
GB2547816B (en) | Actually-measured marine environment data assimilation method based on sequence recursive filtering three-dimensional variation | |
CN111275250B (zh) | 一种基于各向异性的强流区海表面温度预报方法 | |
Mackman et al. | Investigation of an adaptive sampling method for data interpolation using radial basis functions | |
CN112861445A (zh) | 一种无网格数值模型的模拟方法 | |
Tey et al. | Moving least squares method and its improvement: A concise review | |
Liu et al. | Machine-learning-based prediction of regularization parameters for seismic inverse problems | |
CN110363713B (zh) | 基于递归样本缩放和双线性因子分解的高光谱图像降噪方法 | |
Hopkins | A simple generative model of trilobite segmentation and growth | |
Fuchs et al. | Scale-space approaches to FTLE ridges | |
Liu et al. | A novel framework for predicting active flow control by combining deep reinforcement learning and masked deep neural network | |
CN116306206A (zh) | 一种基于深度神经网络的翼型跨音速抖振流场快速预测方法 | |
CN114692529B (zh) | 一种cfd高维响应的不确定度量化方法、装置、计算机设备 | |
CN112541292B (zh) | 一种基于分布式光纤测温原理的海底电缆埋深估计算法 | |
Yusoff et al. | Flow visualization techniques: A review | |
Romano et al. | Convolutional-neural-network-based auto-encoder for synthetic upscaling of computational fluid dynamics simulations | |
CN113868853A (zh) | 一种梯度增强变保真代理模型建模方法 | |
CN113869350A (zh) | 基于空间特征差异的海流预测方法以及*** | |
He et al. | Aerodynamic data fusion with a multi-fidelity surrogate modeling method | |
de Leeuw et al. | Bm3d: motion estimation in time dependent volume data | |
Hurtado-de-Mendoza et al. | Non-modal analysis of linear multigrid schemes for the high-order Flux Reconstruction method | |
KR102364390B1 (ko) | 매트랩을 활용한 모델링 자동화 장치 및 방법 | |
CN115795696B (zh) | 容冰翼型设计过程中冰形的生成方法、装置、设备及介质 | |
CN117077551B (zh) | 基于梯度增强随机Co-Kriging模型的非嵌入式不确定性量化方法 | |
CN112084714B (zh) | 一种基于数据驱动的整车电磁干扰抑制方法 | |
CN109242937B (zh) | 将Mike流场结果文件转换成shapefile矢量文件的方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220517 |