CN109902136A - 一种多个海洋学模型参数的同时率定方法 - Google Patents
一种多个海洋学模型参数的同时率定方法 Download PDFInfo
- Publication number
- CN109902136A CN109902136A CN201910159901.8A CN201910159901A CN109902136A CN 109902136 A CN109902136 A CN 109902136A CN 201910159901 A CN201910159901 A CN 201910159901A CN 109902136 A CN109902136 A CN 109902136A
- Authority
- CN
- China
- Prior art keywords
- coefficient
- empirical
- depth
- wind
- friction factor
- 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.)
- Granted
Links
Classifications
-
- 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
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
本发明提供了一种多个海洋学模型参数的同时率定方法,包括构建地波雷达观测系统,获取表层流场观测数据;获取研究海域海水层结和水深,由表层流场推导出深度平均的流场;构建球坐标系下考虑风应力和海表面气压的二维正压海洋伴随同化模型;给定二维正压海洋伴随同化模型中风应力拖曳系数和底摩擦系数的经验公式,给出经验公式中经验系数调整策略;以二维正压海洋伴随同化模型为基础,同化深度平均的流场数据,同时对风应力拖曳系数和底摩擦系数经验公式中的经验系数进行优化调整,从而实现对于风应力拖曳系数和底摩擦系数的同时率定。本发明既保证了模型参数具有一定的物理意义,也实现了多个模型参数的同时、高效率定。
Description
技术领域
本发明涉及海洋学模型参数率定领域,具体涉及一种多个海洋学模型参数的同时率定方法。
背景技术
海洋学模型通过求解Navier-Stokes控制方程,获得海洋环境变量的时空分布,具有定量、预报性强的特色,在基础学科研究和资源开发活动中起着越来越重要的作用。此外,精确的海洋数值预报结果是开展海上溢油漂移预测和海上遇险目标漂移预测的基础。
海洋学模型中,存在着众多需要参数化的过程。海洋上边界处的风场是海洋能量的重要来源,海面风场对于海洋的影响常用依赖风速的风应力计算公式进行参数化,其中风应力拖曳系数的选取对于风应力的计算至关重要。底部摩擦效应对于海洋能量耗散的影响,是一个经典的海洋学问题,底摩擦效应常用二次形式的依赖于流速的计算公式进行参数化,其中底摩擦系数的取值很大程度上影响能量耗散的计算结果。可以看出,风应力拖曳系数和底摩擦系数的合理取值及相互匹配,是决定海洋数值模拟结果准确与否的重要因素。
众多学者基于实验室实验或者海上观测实验,单独给出了风应力拖曳系数或者底摩擦系数的经验公式。一般情况下,经验公式是根据一定的物理意义给出的,但其中的经验系数需要根据实验数据进行率定。而经验系数在不同的实验条件、不同的海域、不同的海洋环境下,取值往往不同,因此,不假思索的直接采用文献中的经验系数,可能会造成数值模拟结果与实际情况相差较大。此外,在利用海洋学模型开展数值模拟时,传统的试错法很难实现对于多个模型参数进行高效、同步率定。
在此背景下,如何利用实验数据或者观测数据,对研究海域的风应力拖曳系数和底摩擦系数经验公式中的经验系数同时进行自动化率定,是进一步提高海洋数值模拟精度和预报精度的重要问题,也是本发明主要解决的技术问题。
发明内容
有鉴于此,本发明提供了一种多个海洋学模型参数的同时率定方法,解决了上述问题中的至少一个。
为了达到上述目的,本发明的技术方案是这样实现的:
一种多个海洋学模型参数的同时率定方法,包括:
构建地波雷达观测系统,获取表层流场观测数据;
获取研究海域海水层结和水深,由表层流场推导出深度平均的流场;
构建球坐标系下考虑风应力和海表面气压的二维正压海洋伴随同化模型;
给定二维正压海洋伴随同化模型中风应力拖曳系数和底摩擦系数的经验公式,给出经验公式中经验系数调整策略;
以二维正压海洋伴随同化模型为基础,同化深度平均的流场数据,根据给出的经验系数调整策略,同时对风应力拖曳系数和底摩擦系数经验公式中的经验系数进行优化调整,从而实现对于风应力拖曳系数和底摩擦系数的同时率定。
可选地,获取研究海域海水层结和水深,由表层流场推导出深度平均的流场;包括:
根据研究海域海水层结和水深,求解垂向本征方程得到垂向本征函数;
利用模态分解理论和最小二乘理论,由表层流场构造全水深的流场,进而得到深度平均的流场。
可选地,构建球坐标系下考虑风应力和海表面气压的二维正压海洋伴随同化模型;包括:
所述二维正压海洋伴随同化模型包含正向模块、反向模块和梯度模块。
可选地,所述正向模块的控制方程为:
其中,t代表时间,λ和φ分别是经度和纬度,z为深度,自静止海表面向上为正,u和v分别表示λ和φ方向上深度平均流速分量,p为大气压强,ρ为水体密度,a=Rcosφ,R是地球半径,g是重力加速度,Ω为地球自转角速度,f=2Ωsinφ是Coriolis参数,Ah是水平涡动粘性系数,ρa为空气密度,Cs为风应力拖曳系数,W为海面10米高度的风速,h为静水深,ζ为海表面起伏,为引潮势,Cd为底摩擦系数,Δ是拉普拉斯算子。
可选地,反向模块的控制方程是针对于同化流速数据而推导出来的,具体为:
其中,t代表时间,λ和φ分别是经度和纬度,z为深度,自静止海表面向上为正,u和v分别表示λ和φ方向上深度平均流速分量,p为大气压强,ρ为水体密度,a=Rcosφ,R是地球半径,g是重力加速度,Ω为地球自转角速度,f=2Ωsinφ是Coriolis参数,Ah是水平涡动粘性系数,ρa为空气密度,Cs为风应力拖曳系数,W为海面10米高度的风速,h为静水深,ζ为海表面起伏,为引潮势,Cd为底摩擦系数,Ku和Kv是权重矩阵,理论上是观测误差矩阵的逆矩阵,u和为深度平均流速u和v的观测值,ζa、ua和va分别是ζ、u和v的伴随变量,Δ是拉普拉斯算子。
可选地,给定二维正压海洋伴随同化模型中风应力拖曳系数和底摩擦系数的经验公式,给出经验公式中经验系数调整策略;包括:
给定的风应力拖曳系数的经验公式为:
其中,Cs为风应力拖曳系数,W为10米高度风速,C0为风速小于11米时的常数分布的风应力拖曳经验系数,k、b、A、B、C为经验系数。
可选地,给定二维正压海洋伴随同化模型中风应力拖曳系数和底摩擦系数的经验公式,给出经验公式中经验系数调整策略;包括:
给定的底摩擦系数的经验公式为:
其中,Cd为底摩擦系数,g为重力加速度,Cchezy为Chezy系数,n为曼宁经验系数,α为经验系数。
可选地,给定二维正压海洋伴随同化模型中风应力拖曳系数和底摩擦系数的经验公式,给出经验公式中经验系数调整策略;包括:
以所述梯度模块为基础,经验公式中经验系数调整策略如下:
其中,pc代表任一经验系数,pcold为优化调整前的pc值,pcnew为优化调整后的pc值,ss为调整步长,为代价函数关于经验系数的梯度。
可选地,所述代价函数关于经验系数的梯度是根据所述梯度模块和分步求导公式求得的。
可选地,以二维正压海洋伴随同化模型为基础,同化深度平均的流场数据,根据给出的经验系数调整策略,同时对风应力拖曳系数和底摩擦系数经验公式中的经验系数进行优化调整,从而实现对于风应力拖曳系数和底摩擦系数的同时率定,包括:
风应力拖曳系数和底摩擦系数的同时率定是指,对风应力拖曳系数和底摩擦系数经验公式中的经验系数进行同时率定,从而实现对具有一定物理意义的风应力拖曳系数和底摩擦系数的同时率定。
与现有技术相比,本发明的优点是:本发明利用伴随同化方法同化地波雷达观测数据,将参数率定问题转化为数值迭代过程,能够自适应地同时率定多个海洋学模型参数;此外,通过引入风应力拖曳系数和底摩擦系数的经验公式,对公式中的经验系数进行率定,而不是直接对风应力拖曳系数和底摩擦系数进行率定,避免了传统伴随同化方法率定参数时不考虑物理意义的弊端。基于地波雷达观测数据的多个海洋学模型参数的同时率定,为进一步提高海洋数值预报精度奠定了坚实的基础。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是本发明一种实施方式的流程图;
图2是本发明对风应力拖曳系数和底摩擦系数经验公式中的经验系数进行优化调整的流程示意图。
具体实施方式
下面结合附图对本发明的较佳实施例进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
如图1所示,一种多个海洋学模型参数的同时率定方法,包括以下步骤:
S100:构建地波雷达观测系统,获取表层流场观测数据;
具体地,地波雷达是一种新兴的海洋监测技术,具有超视距、大范围、全天候以及低成本等优点,能够从雷达回波中提取风场、浪场、流场等海况信息;本步骤中,可以构造多站地波雷达组网的观测系统,从而提高表层流场的观测精度。
S200:获取研究海域海水层结和水深,由表层流场推导出深度平均的流场。
具体地,包括:
S201:根据研究海域海水层结和水深,求解垂向本征方程得到垂向本征函数。
其中,垂向本征方程如下所述:
其中,ψ为垂向本征函数,z为垂向坐标(从海表面起,向上为正),N为表征海水层结的浮频率,h为水深,c为波速。由公式(1)可知,垂向本征函数具有正交性。
利用WOA数据集或者HYCOM模式实时预报系统提供的温度、盐度数据,可以计算出研究海域的浮频率。由浅水波的频散关系可求得正压模态的波速;利用WKB近似,可以得到斜压模态的波速。将浮频率和波速带入垂向本征方程(1),可以求解出不同模态的垂向本征函数。
S202:利用模态分解理论和最小二乘理论,由表层流场构造全水深的流场,进而得到深度平均的流场。
其中,所述模态分解理论为:
其中,u和v分别为东西向和南北向流速,t为时间,um和vm是第m模态(m=0为正压模态,m>0为斜压模态)的流速分量,ψm为第m模态的垂向本征函数(即式(1)的解),Um和Vm分别是um和vm对应的时间变化系数。
最小二乘理论是使得地波雷达观测系统所得表层流速观测值和利用前M个垂向模态构造的表层流速值之间的误差平方达到最小(以流速分量u为例):
其中,z1代表表层,M为所选取的模态数(包含第0模态的正压模态和前M个斜压模态)。
通过将已经求解得到的前M个垂向本征函数带入式(3),可在最小二乘意义下,求得Um的值;进而利用模态分解公式(2),可以得到全水深的流场。将全水深的流场进行垂向加权平均,可以得到垂向平均的流场。
S300:构建球坐标系下考虑风应力和海表面气压的二维正压海洋伴随同化模型;
具体地,二维正压海洋伴随同化模型包含正向模块、反向模块和梯度模块。
其中,正向模块的控制方程为:
其中,t代表时间,λ和φ分别是经度和纬度,z为深度,自静止海表面向上为正,u和v分别表示λ和φ方向上深度平均流速分量,p为大气压强,ρ为水体密度(1025kg·m-3),a=Rcosφ,R是地球半径,g是重力加速度,Ω为地球自转角速度,f=2Ωsinφ是Coriolis参数,Ah是水平涡动粘性系数,ρa为空气密度(1.27kg·m-3),Cs为风应力拖曳系数,W为海面10米高度的风速,h为静水深,ζ为海表面起伏,为引潮势,Cd为底摩擦系数,Δ是拉普拉斯算子:
正向模块的初始条件采用HYCOM模式实时预报系统提供的流场数据,上边界条件采用NCEP气候预测系统提供的再分析风场、气压场,开边界条件采用OSU全球潮汐模型提供的M2、S2、K1、O1、N2、K2、P1和Q1八个主要分潮的水位数据,闭边界条件为法向流速为0,水深由ETOPO5数据插值得到。
为构建反向模块和梯度模块,特构造表征深度平均流场模拟结果和观测结果之间误差的代价函数:
其中,J(u,v)为代价函数,Ku和Kv是权重矩阵,理论上是观测误差矩阵的逆矩阵,u和为深度平均流速u和v的观测值,Σ是正向模块的时空计算区域。
进而根据拉格朗日乘子法,构造拉格朗日函数:
其中,ζa、ua和va分别是ζ、u和v的伴随变量。
令拉格朗日函数关于伴随变量ζa、ua和va的偏导数为0,可以得到正向模块的控制方程(式(4))。令拉格朗日函数关于模型变量ζ、u和v的偏导数为0,可以得到反向模块(伴随模型)的控制方程,具体如下:
反向模块的初始条件为伴随变量ζa、ua和va的值为0,侧边界条件为伴随变量ua和va的法向分量为0。
令拉格朗日函数关于风应力拖曳系数和底摩擦系数的偏导数为0,可以求出梯度模块(代价函数关于模型参数的梯度表达式),具体如下:
S400:给定二维正压海洋伴随同化模型中风应力拖曳系数和底摩擦系数的经验公式,给出经验公式中经验系数调整策略;
本步骤中,所述风应力拖曳系数经验公式采用依赖于风速的分段函数形式,具体如下:
其中,Cs为风应力拖曳系数,W为10米高度风速,C0为风速小于11米时的常数分布的风应力拖曳经验系数,k、b、A、B、C为经验系数。上述经验系数C0、k、b、A、B、C均假定为常数、需要进行率定。
所述底摩擦系数的经验公式采用依赖于水深的形式,具体如下:
其中,Cd为底摩擦系数,g为重力加速度,Cchezy为Chezy系数,n为曼宁经验系数,α为经验系数,假定上述经验系数均为常数且需要进行率定。
本步骤中,根据梯度模块和分步求导公式,可知代价函数关于经验公式(式(10)和式(11))中经验系数的梯度表达式为:
根据最速下降法,给出经验公式中经验系数调整策略,如下:
其中,pc代表任一经验系数,pcold为优化调整前的pc值,pcnew为优化调整后的pc值,ss为调整步长,为代价函数关于经验系数的梯度(即:式(12))。
S500:以二维正压海洋伴随同化模型为基础,同化深度平均的流场数据,根据给出的经验系数调整策略,同时对风应力拖曳系数和底摩擦系数经验公式中的经验系数进行优化调整,从而实现对于风应力拖曳系数和底摩擦系数的同时率定。
具体地,风应力拖曳系数和底摩擦系数的同时率定是指,不是对其中的一个参数进行单独率定,也不是直接假定参数是常数或者空间分布的形式直接进行率定,而是对风应力拖曳系数和底摩擦系数经验公式中的经验系数进行同时率定,从而实现对具有一定物理意义的风应力拖曳系数和底摩擦系数的同时率定。
如图2所示,对风应力拖曳系数和底摩擦系数经验公式中的经验系数进行优化调整的流程可细分为如下步骤:
S501:根据参考文献,给定风应力拖曳系数经验公式(式(10))和底摩擦系数经验公式(式(11))中经验系数C0、k、b、A、B、C、n、α的初始猜测值;
S502:对正向模块(式(4))进行数值离散,在时间上正向积分正向模块,然后在时间上反向积分反向模块(式(8)),进而计算出步骤S300中代价函数值(式(6))和步骤S400中代价函数关于经验系数的梯度值(式(12));
S503:根据步骤S400中经验系数调整策略(式(13)),同时对风应力拖曳系数和底摩擦系数经验公式中的经验系数进行优化调整;
S504:重复步骤S502、S503,直至相邻两次的代价函数的差小于1×10-4或者重复次数达到100次;将所得经验系数C0、k、b、A、B、C、n、α的值带入步骤S400中风应力拖曳系数经验公式(式10)和底摩擦系数经验公式(式(11)),即可得到最终基于地波雷达观测数据所率定出的风应力拖曳系数和底摩擦系数。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (10)
1.一种多个海洋学模型参数的同时率定方法,其特征在于,包括:
构建地波雷达观测系统,获取表层流场观测数据;
获取研究海域海水层结和水深,由表层流场推导出深度平均的流场;
构建球坐标系下考虑风应力和海表面气压的二维正压海洋伴随同化模型;
给定二维正压海洋伴随同化模型中风应力拖曳系数和底摩擦系数的经验公式,给出经验公式中经验系数调整策略;
以二维正压海洋伴随同化模型为基础,同化深度平均的流场数据,根据给出的经验系数调整策略,同时对风应力拖曳系数和底摩擦系数经验公式中的经验系数进行优化调整,从而实现对于风应力拖曳系数和底摩擦系数的同时率定。
2.根据权利要求1所述的同时率定方法,其特征在于,获取研究海域海水层结和水深,由表层流场推导出深度平均的流场;包括:
根据研究海域海水层结和水深,求解垂向本征方程得到垂向本征函数;
利用模态分解理论和最小二乘理论,由表层流场构造全水深的流场,进而得到深度平均的流场。
3.根据权利要求1所述的的同时率定方法,其特征在于,构建球坐标系下考虑风应力和海表面气压的二维正压海洋伴随同化模型;包括:
所述二维正压海洋伴随同化模型包含正向模块、反向模块和梯度模块。
4.根据权利要求3所述的同时率定方法,其特征在于,所述正向模块的控制方程为:
其中,t代表时间,λ和φ分别是经度和纬度,z为深度,自静止海表面向上为正,u和v分别表示λ和φ方向上深度平均流速分量,p为大气压强,ρ为水体密度,a=R cosφ,R是地球半径,g是重力加速度,Ω为地球自转角速度,f=2Ωsinφ是Coriolis参数,Ah是水平涡动粘性系数,ρa为空气密度,Cs为风应力拖曳系数,W为海面10米高度的风速,h为静水深,ζ为海表面起伏,为引潮势,Cd为底摩擦系数,Δ是拉普拉斯算子。
5.根据权利要求3所述的同时率定方法,其特征在于,反向模块的控制方程是针对于同化流速数据而推导出来的,具体为:
其中,t代表时间,λ和φ分别是经度和纬度,z为深度,自静止海表面向上为正,u和v分别表示λ和φ方向上深度平均流速分量,p为大气压强,ρ为水体密度,a=R cosφ,R是地球半径,g是重力加速度,Ω为地球自转角速度,f=2Ωsinφ是Coriolis参数,Ah是水平涡动粘性系数,ρa为空气密度,Cs为风应力拖曳系数,W为海面10米高度的风速,h为静水深,ζ为海表面起伏,为引潮势,Cd为底摩擦系数,Ku和Kv是权重矩阵,理论上是观测误差矩阵的逆矩阵,u和为深度平均流速u和v的观测值,ζa、ua和va分别是ζ、u和v的伴随变量,Δ是拉普拉斯算子。
6.根据权利要求1所述的的同时率定方法,其特征在于,给定二维正压海洋伴随同化模型中风应力拖曳系数和底摩擦系数的经验公式,给出经验公式中经验系数调整策略;包括:
给定的风应力拖曳系数的经验公式为:
其中,Cs为风应力拖曳系数,W为10米高度风速,C0为风速小于11米时的常数分布的风应力拖曳经验系数,k、b、A、B、C为经验系数。
7.根据权利要求1所述的的同时率定方法,其特征在于,给定二维正压海洋伴随同化模型中风应力拖曳系数和底摩擦系数的经验公式,给出经验公式中经验系数调整策略;包括:
给定的底摩擦系数的经验公式为:
其中,Cd为底摩擦系数,g为重力加速度,Cchezy为Chezy系数,n为曼宁经验系数,α为经验系数。
8.根据权利要求3所述的的同时率定方法,其特征在于,给定二维正压海洋伴随同化模型中风应力拖曳系数和底摩擦系数的经验公式,给出经验公式中经验系数调整策略;包括:
以所述梯度模块为基础,经验公式中经验系数调整策略如下:
其中,pc代表任一经验系数,pcold为优化调整前的pc值,pcnew为优化调整后的pc值,ss为调整步长,为代价函数关于经验系数的梯度。
9.根据权利要求8所述的同时率定方法,其特征在于,所述代价函数关于经验系数的梯度是根据所述梯度模块和分步求导公式求得的。
10.根据权利要求1所述的同时率定方法,其特征在于,以二维正压海洋伴随同化模型为基础,同化深度平均的流场数据,根据给出的经验系数调整策略,同时对风应力拖曳系数和底摩擦系数经验公式中的经验系数进行优化调整,从而实现对于风应力拖曳系数和底摩擦系数的同时率定,包括:
风应力拖曳系数和底摩擦系数的同时率定是指,对风应力拖曳系数和底摩擦系数经验公式中的经验系数进行同时率定,从而实现对具有一定物理意义的风应力拖曳系数和底摩擦系数的同时率定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910159901.8A CN109902136B (zh) | 2019-03-04 | 2019-03-04 | 一种多个海洋学模型参数的同时率定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910159901.8A CN109902136B (zh) | 2019-03-04 | 2019-03-04 | 一种多个海洋学模型参数的同时率定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109902136A true CN109902136A (zh) | 2019-06-18 |
CN109902136B CN109902136B (zh) | 2020-02-21 |
Family
ID=66946199
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910159901.8A Active CN109902136B (zh) | 2019-03-04 | 2019-03-04 | 一种多个海洋学模型参数的同时率定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109902136B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111797578A (zh) * | 2020-06-30 | 2020-10-20 | 自然资源部第一海洋研究所 | 一种定量计算降雨致海面应力的方法 |
CN113076373A (zh) * | 2021-02-26 | 2021-07-06 | 广东科诺勘测工程有限公司 | 一种海域流场及冲淤深度实时水文监测大数据展示与空间查询方法及系统 |
CN113378443A (zh) * | 2021-08-12 | 2021-09-10 | 中国地质大学(武汉) | 一种地波雷达数据融合同化方法及计算机设备 |
CN113901741A (zh) * | 2021-10-20 | 2022-01-07 | 哈尔滨工业大学(威海) | 海洋-电离层耗散动力学仿真方法、仿真系统、设备及应用 |
CN114814276A (zh) * | 2022-03-21 | 2022-07-29 | 汕头大学 | 海上风电设备运行导致周边海水垂向运动流速的计算方法 |
CN115081254A (zh) * | 2022-08-19 | 2022-09-20 | 中交第四航务工程勘察设计院有限公司 | 全球波浪数学模型分区块高效率定方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140032191A1 (en) * | 2010-10-26 | 2014-01-30 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | Tool for rapid configuration of a river model using imagery-based information |
CN107729656A (zh) * | 2017-10-19 | 2018-02-23 | 福建四创软件有限公司 | 基于台风过程风场作用下的风暴潮漫堤风险预警方法 |
-
2019
- 2019-03-04 CN CN201910159901.8A patent/CN109902136B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140032191A1 (en) * | 2010-10-26 | 2014-01-30 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | Tool for rapid configuration of a river model using imagery-based information |
CN107729656A (zh) * | 2017-10-19 | 2018-02-23 | 福建四创软件有限公司 | 基于台风过程风场作用下的风暴潮漫堤风险预警方法 |
Non-Patent Citations (2)
Title |
---|
ALEX E. HAY 等: "Turbulence and Drag in a High Reynolds Number Tidal Passage Targetted for In-Stream Tidal Power", 《IEEE》 * |
黄士力 等: "上海黄浦江及沿海高分辨率非结构网格风暴潮数值模型", 《海洋预报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111797578A (zh) * | 2020-06-30 | 2020-10-20 | 自然资源部第一海洋研究所 | 一种定量计算降雨致海面应力的方法 |
CN111797578B (zh) * | 2020-06-30 | 2022-05-06 | 自然资源部第一海洋研究所 | 一种定量计算降雨致海面应力的方法 |
CN113076373A (zh) * | 2021-02-26 | 2021-07-06 | 广东科诺勘测工程有限公司 | 一种海域流场及冲淤深度实时水文监测大数据展示与空间查询方法及系统 |
CN113076373B (zh) * | 2021-02-26 | 2023-11-21 | 广东科诺勘测工程有限公司 | 一种海域流场及冲淤深度实时水文监测大数据展示与空间查询方法及系统 |
CN113378443A (zh) * | 2021-08-12 | 2021-09-10 | 中国地质大学(武汉) | 一种地波雷达数据融合同化方法及计算机设备 |
CN113378443B (zh) * | 2021-08-12 | 2021-11-05 | 中国地质大学(武汉) | 一种地波雷达数据融合同化方法及计算机设备 |
CN113901741A (zh) * | 2021-10-20 | 2022-01-07 | 哈尔滨工业大学(威海) | 海洋-电离层耗散动力学仿真方法、仿真系统、设备及应用 |
CN113901741B (zh) * | 2021-10-20 | 2023-08-18 | 哈尔滨工业大学(威海) | 海洋-电离层耗散动力学仿真方法、仿真系统、设备 |
CN114814276A (zh) * | 2022-03-21 | 2022-07-29 | 汕头大学 | 海上风电设备运行导致周边海水垂向运动流速的计算方法 |
CN114814276B (zh) * | 2022-03-21 | 2023-08-18 | 汕头大学 | 海上风电设备运行导致周边海水垂向运动流速的计算方法 |
CN115081254A (zh) * | 2022-08-19 | 2022-09-20 | 中交第四航务工程勘察设计院有限公司 | 全球波浪数学模型分区块高效率定方法及装置 |
CN115081254B (zh) * | 2022-08-19 | 2022-11-15 | 中交第四航务工程勘察设计院有限公司 | 全球波浪数学模型分区块高效率定方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109902136B (zh) | 2020-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109902136A (zh) | 一种多个海洋学模型参数的同时率定方法 | |
Yamazaki et al. | A physically based description of floodplain inundation dynamics in a global river routing model | |
CN109657418A (zh) | 一种基于mike21的湖泊水环境容量计算方法 | |
Haidvogel et al. | Ocean forecasting in terrain-following coordinates: Formulation and skill assessment of the Regional Ocean Modeling System | |
Westerink et al. | Tide and storm surge predictions using finite element model | |
Haas | Assessment of energy production potential from ocean currents along the United States coastline | |
CN105894106B (zh) | 一种海洋模式和气象模式的一体化耦合方法 | |
Lewis et al. | Assimilation of Doppler radar current data into numerical ocean models | |
Lima et al. | An investigation of ocean model uncertainties through ensemble forecast experiments in the Southwest Atlantic Ocean | |
Fan et al. | Innovative approach to design truncated mooring system based on static and damping equivalent | |
Wang et al. | Estimation of bottom friction coefficient in multi‐constituent tidal models using the adjoint method: temporal variations and spatial distributions | |
Belyaev et al. | Comparison of methods for ARGO drifters data assimilation into a hydrodynamical model of the ocean | |
Zou et al. | Tide-surge and wave interaction in the Gulf of Maine during an extratropical storm | |
Diansky et al. | Numerical simulation of the Caspian sea circulation using the marine and atmospheric research system | |
Berntsen | A perfectly balanced method for estimating the internal pressure gradients in σ-coordinate ocean models | |
Koçyigit et al. | Three-dimensional numerical modelling of wind-driven circulation in a homogeneous lake | |
Zhang et al. | A modelling study of ice effect on tidal damping in the Bohai Sea | |
Zhang et al. | Hydrodynamic model with wave–current interaction in coastal regions | |
Jiang et al. | Estimating smoothly varying open boundary conditions for a 3D internal tidal model with an improved independent point scheme | |
Kurapov et al. | Representer-based analyses in the coastal upwelling system | |
Kordzadze et al. | Coastal forecasting system for the easternmost part of the Black Sea | |
Yan et al. | Seasonal variability of the North Equatorial Current transport in the western Pacific Ocean | |
Parent et al. | Comparative assimilation of Topex/Poseidon and ERS altimeter data and of TAO temperature data in the Tropical Pacific Ocean during 1994–1998, and the mean sea-surface height issue | |
Ma et al. | A high-precision hydrodynamic model coupled with the hydrological habitat suitability model to reveal estuarine vegetation distribution | |
Liu et al. | Impact of assimilation of moored velocity data on low‐frequency current estimation in northwestern tropical Pacific |
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 |