CN108562896A - 一种基于三维正压浅海大陆架模型的深层海流反演方法 - Google Patents
一种基于三维正压浅海大陆架模型的深层海流反演方法 Download PDFInfo
- Publication number
- CN108562896A CN108562896A CN201810366515.1A CN201810366515A CN108562896A CN 108562896 A CN108562896 A CN 108562896A CN 201810366515 A CN201810366515 A CN 201810366515A CN 108562896 A CN108562896 A CN 108562896A
- Authority
- CN
- China
- Prior art keywords
- ocean
- follows
- coordinate system
- equation group
- sea
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
Abstract
一种基于三维正压浅海大陆架模型的深层海流反演方法,属高频地波雷达遥感探测技术领域,采用归一化的坐标系的三维非线性正压浅海大陆架模式,利用高频地波雷达探测到的海洋表面流数据,考虑复杂地形的边界问题,参数化物理过程,简化海洋控制方程组,再建立三维网格,采用差分法求解离散方程组进而反演得到相应的深层流数据信息,包括步骤:选择坐标系;导出归一化坐标系下的控制方程组;参数化物理过程;迭代法求解方程组得出深层流信息;实测数据与反演结果对比分析。本发明对比试验结果发现,实测数据曲线图和反演结果曲线图的整体趋势接近,反演曲线可以反映出海洋深层流变化的整体趋势,符合海洋动力学运动规律。
Description
技术领域
本发明涉及高频地波雷达遥感探测技术领域,尤其涉及一种基于三维正压浅海大陆架模型的深层海流反演方法。
背景技术
为了获得海洋表面以下海流的运动信息,现有技术采用直接接触测量水体的方式,所以获得大面积水域的数据则需要投放相当数量的直接测量仪器。该方法对周围测量环境的要求高,成本高,探测效率低。高频地波雷达遥感技术是近三十年来发展起来的探测海洋表面动力学参数的技术,覆盖范围广,探测精度和效率都较高。但是由于电磁波在水体中的衰减严重,使得高频地波雷达遥感探测技术仅限于探测海洋水体表面,不能直接获得海洋内部水体运动的相关信息是高频地波雷达遥感技术在测量海洋流参数方面的不足。
发明内容
针对现有技术的不足,本发明所解决的问题是利用高频地波雷达探测得到的海洋表面流参数信息即可得到深层流参数信息。
为解决上述技术问题,本发明采用的技术方案是一种基于三维正压浅海大陆架模型的深层海流反演方法,采用归一化的坐标系的三维非线性正压浅海大陆架模式,利用高频地波雷达探测到的海洋表面流数据,考虑复杂地形的边界问题,参数化物理过程,简化海洋控制方程组,再建立三维网格,采用差分法求解离散方程组进而反演得到相应的深层流数据信息,包括如下步骤:
(1)选择σ坐标系;
假设在高频地波雷达所探测的每场数据的累积时间内,海水在垂直方向上处于静力学平衡状态;海水处于正压模式,海水不可压缩,海水的运动速度不随时间变化而变化;为了便于分层,采用归一化的σ坐标系,使水浅处与水深处海水都具有相同的垂向分层数,并且每一层的侧边界不再变化,从海底到海面的高度处处规范化为1;
在σ坐标系下x轴和y轴不变,用σ取代z轴,其中z和σ的转换关系如下:
坐标变换后,设总的海水深度为H=h+ζ,海洋表面z=ζ,σ=0,在海底z=-h,σ=-1;
(2)导出归一化坐标系下的控制方程组;
在σ坐标系中,利用dσ和dz之间的微分系数关系将海洋基本控制方程组中各个方程转化为如下形式:
1)连续方程表达式如下:
2)方向动量方程的非线性项采用能量二次守恒的半动量形式,表达式如下:
3)流体静力学方程如下:
表达式中的u、v、w分别代表流速在x、y、σ方向上的分量;H=h+ζ是总海深;ζ是静止海面状态下的浪高;h是静态水域深度;ρ是海水密度;R是地球半径;是维度;是科氏参数;Ω是地转角速度;AH是水平方向涡动粘滞系数,取常数5000m2/s;AV是垂直方向涡动粘滞系数;ξx,ξy是在水平方向涡动粘滞系数中由于海底地形和水位有关的量,接近无穷小,忽略不计;
由公式(02)得到垂直流速的表达式如下:
在动量方程中,水平压强梯度力方程如下:
在实际观测资料中发现,海水水平密度变化较小,所以在计算中可以忽略水平压强梯度力第二项,简化的水平压强梯度力方程如下:
(3)参数化物理过程;
海洋表面限制条件由高频地波雷达监测的表面流参数值决定,设定海洋表面垂直速度w=0,动力学参数限制条件如下:
其中τa表面风应力矢量,τb是海底摩擦力矢量,表达式如下所示:
τa=ρaCa|va|va (13)
τb=ρbCb|vb|vb (14)
其中ρa是空气密度;va是风速矢量;Ca=(0.73+0.069|va|)×10-3是海水对风的拖曳系数;ρb是海水密度;vb是靠近海底的流速矢量;是海底摩擦应力拖曳系数;k是Kalman常数;Δσ是最后一层的厚度;z0是海底粗糙度;根据Prandtl混合长度理论,得垂直方向涡动粘滞系数表达式如下:
其中代表分子粘性系数;代表混合长度;s是反映海底粗糙度的参数,z'=(1+σ)H代表考察点到海底的距离;
(4)建立三维网格,离散化方程组;
采用以下方式分别将u、v、w在x、y、σ方向上进行离散,以u为例:
u的值为同一网格中四个结点的加权平均值,θ是空间加权系数,当θ<0.5时,上述表达式为一阶,当θ=0.5时,上述表达式为二阶,采用二阶精度,其格式如下:
(5)迭代法求解方程组得出深层流信息;
在公式(18)中,后一层结点上的值不仅与上一层结点值有关,还和本层已计算出的结点值有关,利用计算机,运用迭代法求解离散化的海洋控制方程组;
(6)实测数据与反演结果对比分析;
分析数值反演角度图发现,随着海流深度的增加,实测海流的流向和反演得出的海流的流向都呈逐渐右偏的趋势;分析数值反演速度图发现,随着海流深度的增加,实测海流的流速和反演得出的海流的流速都呈逐渐减小的趋势,在靠近海底处,由于底部摩擦力的作用,流速都趋近于零。
本发明的有益效果:对比试验结果发现,实测数据曲线图和反演结果曲线图的整体趋势接近,反演曲线可以反映出海洋深层流变化的整体趋势,符合海洋动力学运动规律。
附图说明
图1是本发明流程示意图;
图2是坐标变换前的坐标系示意图;
图3是坐标变换后的坐标系示意图;
图4是探测点在14:10的角度变化图;
图5是探测点在14:30的角度变化图;
图6是探测点在14:10的速度变化图;
图7是探测点在14:30的速度变化图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的说明,但不是对本发明的限定。
图1示出了一种基于三维正压浅海大陆架模型的深层海流反演方法,采用归一化的坐标系的三维非线性正压浅海大陆架模式,利用高频地波雷达探测到的海洋表面流数据,考虑复杂地形的边界问题,参数化物理过程,简化海洋控制方程组,再建立三维网格,采用差分法求解离散方程组进而反演得到相应的深层流数据信息,包括如下步骤:
(1)选择σ坐标系;
假设在高频地波雷达所探测的每场数据的累积时间内,海水在垂直方向上处于静力学平衡状态;海水处于正压模式,海水不可压缩,海水的运动速度不随时间变化而变化;为了便于分层,把如图2所示的变换前的坐标系变换为归一化的σ坐标系,如图3所示,使水浅处与水深处海水都具有相同的垂向分层数,并且每一层的侧边界不再变化,从海底到海面的高度处处规范化为1;
在σ坐标系下x轴和y轴不变,用σ取代z轴,其中z和σ的转换关系如下:
坐标变换后,设总的海水深度为H=h+ζ,海洋表面z=ζ,σ=0,在海底z=-h,σ=-1;
(2)导出归一化坐标系下的控制方程组;
在σ坐标系中,利用dσ和dz之间的微分系数关系将海洋基本控制方程组中各个方程转化为如下形式:
1)连续方程表达式如下:
2)方向动量方程的非线性项采用能量二次守恒的半动量形式,表达式如下:
3)流体静力学方程如下:
表达式中的u、v、w分别代表流速在x、y、σ方向上的分量;H=h+ζ是总海深;ζ是静止海面状态下的浪高;h是静态水域深度;ρ是海水密度;R是地球半径;是维度;是科氏参数;Ω是地转角速度;AH是水平方向涡动粘滞系数,取常数5000m2/s;AV是垂直方向涡动粘滞系数;ξx,ξy是在水平方向涡动粘滞系数中由于海底地形和水位有关的量,接近无穷小,忽略不计;
由公式(02)得到垂直流速的表达式如下:
在动量方程中,水平压强梯度力方程如下:
在实际观测资料中发现,海水水平密度变化较小,所以在计算中可以忽略水平压强梯度力第二项,简化的水平压强梯度力方程如下:
(3)参数化物理过程;
海洋表面限制条件由高频地波雷达监测的表面流参数值决定,设定海洋表面垂直速度w=0,动力学参数限制条件如下:
其中τa表面风应力矢量,τb是海底摩擦力矢量,表达式如下所示:
τa=ρaCa|va|va (13)
τb=ρbCb|vb|vb (14)
其中ρa是空气密度;va是风速矢量;Ca=(0.73+0.069|va|)×10-3是海水对风的拖曳系数;ρb是海水密度;vb是靠近海底的流速矢量;是海底摩擦应力拖曳系数;k是Kalman常数;Δσ是最后一层的厚度;z0是海底粗糙度;根据Prandtl混合长度理论,得垂直方向涡动粘滞系数表达式如下:
其中代表分子粘性系数;代表混合长度;s是反映海底粗糙度的参数,z'=(1+σ)H代表考察点到海底的距离;
(4)建立三维网格,离散化方程组;
采用以下方式分别将u、v、w在x、y、σ方向上进行离散,以u为例:
u的值为同一网格中四个结点的加权平均值,θ是空间加权系数,当θ<0.5时,上述表达式为一阶,当θ=0.5时,上述表达式为二阶,采用二阶精度,其格式如下:
(5)迭代法求解方程组得出深层流信息;
在公式(18)中,后一层结点上的值不仅与上一层结点值有关,还和本层已计算出的结点值有关,利用计算机,运用迭代法求解离散化的海洋控制方程组;
(6)实测数据与反演结果对比分析;
实验结果图中实线代表海流计的实测数据,虚线代表利用本发明方法反演得出的结果;
图4和图5分别为探测点在14:10和14:30两个不同时刻数值反演的角度图。分析数值反演角度图发现,随着海流深度的增加,实测海流的流向和反演得出的海流的流向都呈逐渐右偏的趋势;图6和图7分别为探测点在14:10和14:30两个不同时刻数值反演的速度图。分析数值反演速度图发现,随着海流深度的增加,实测海流的流速和反演得出的海流的流速都呈逐渐减小的趋势,在靠近海底处,由于底部摩擦力的作用,流速都趋近于零;
随着深度的变化,反演海流的整个运动趋势是符合海洋动力学规律的;虽然实测海流数据的流向变化图和流速变化图在小幅度上都出现了反复波动,但并不影响随着深度增加,流向整体右偏和流速逐渐减小的整体趋势;
综上,实测数据曲线图和反演结果曲线图的整体趋势接近,反演曲线可以反映出海洋深层流变化的整体趋势,符合海洋动力学运动规律。
以上结合附图对本发明的实施方式做出了详细说明,但本发明不局限于所描述的实施方式。对于本领域技术人员而言,在不脱离本发明的原理和精神的情况下,对这些实施方式进行各种变化、修改、替换和变型仍落入本发明的保护范围内。
Claims (1)
1.一种基于三维正压浅海大陆架模型的深层海流反演方法,采用归一化的坐标系的三维非线性正压浅海大陆架模式,利用高频地波雷达探测到的海洋表面流数据,考虑复杂地形的边界问题,参数化物理过程,简化海洋控制方程组,再建立三维网格,采用差分法求解离散方程组进而反演得到相应的深层流数据信息,包括如下步骤:
(1)选择σ坐标系;
假设在高频地波雷达所探测的每场数据的累积时间内,海水在垂直方向上处于静力学平衡状态;海水处于正压模式,海水不可压缩,海水的运动速度不随时间变化而变化;为了便于分层,采用归一化的σ坐标系,使水浅处与水深处海水都具有相同的垂向分层数,并且每一层的侧边界不再变化,从海底到海面的高度处处规范化为1;
在σ坐标系下x轴和y轴不变,用σ取代z轴,其中z和σ的转换关系如下:
坐标变换后,设总的海水深度为H=h+ζ,海洋表面z=ζ,σ=0,在海底z=-h,σ=-1;
(2)导出归一化坐标系下的控制方程组;
在σ坐标系中,利用dσ和dz之间的微分系数关系将海洋基本控制方程组中各个方程转化为如下形式:
1)连续方程表达式如下:
2)方向动量方程的非线性项采用能量二次守恒的半动量形式,表达式如下:
3)流体静力学方程如下:
表达式中的u、v、w分别代表流速在x、y、σ方向上的分量;H=h+ζ是总海深;ζ是静止海面状态下的浪高;h是静态水域深度;ρ是海水密度;R是地球半径;是维度;是科氏参数;Ω是地转角速度;AH是水平方向涡动粘滞系数,取常数5000m2/s;AV是垂直方向涡动粘滞系数;ξx,ξy是在水平方向涡动粘滞系数中由于海底地形和水位有关的量,接近无穷小,忽略不计;
由公式(02)得到垂直流速的表达式如下:
在动量方程中,水平压强梯度力方程如下:
在实际观测资料中发现,海水水平密度变化较小,所以在计算中可以忽略水平压强梯度力第二项,简化的水平压强梯度力方程如下:
(3)参数化物理过程;
海洋表面限制条件由高频地波雷达监测的表面流参数值决定,设定海洋表面垂直速度w=0,动力学参数限制条件如下:
其中τa表面风应力矢量,τb是海底摩擦力矢量,表达式如下所示:
τa=ρaCa|va|va (13)
τb=ρbCb|vb|vb (14)
其中ρa是空气密度;va是风速矢量;Ca=(0.73+0.069|va|)×10-3是海水对风的拖曳系数;ρb是海水密度;vb是靠近海底的流速矢量;是海底摩擦应力拖曳系数;k是Kalman常数;Δσ是最后一层的厚度;z0是海底粗糙度;根据Prandtl混合长度理论,得垂直方向涡动粘滞系数表达式如下:
其中代表分子粘性系数;代表混合长度;s是反映海底粗糙度的参数,z'=(1+σ)H代表考察点到海底的距离;
(4)建立三维网格,离散化方程组;
采用以下方式分别将u、v、w在x、y、σ方向上进行离散,以u为例:
u的值为同一网格中四个结点的加权平均值,θ是空间加权系数,当θ<0.5时,上述表达式为一阶,当θ=0.5时,上述表达式为二阶,采用二阶精度,其格式如下:
(5)迭代法求解方程组得出深层流信息;
在公式(18)中,后一层结点上的值不仅与上一层结点值有关,还和本层已计算出的结点值有关,利用计算机,运用迭代法求解离散化的海洋控制方程组;
(6)实测数据与反演结果对比分析;
分析数值反演角度图发现,随着海流深度的增加,实测海流的流向和反演得出的海流的流向都呈逐渐右偏的趋势;分析数值反演速度图发现,随着海流深度的增加,实测海流的流速和反演得出的海流的流速都呈逐渐减小的趋势,在靠近海底处,由于底部摩擦力的作用,流速都趋近于零。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810366515.1A CN108562896A (zh) | 2018-04-23 | 2018-04-23 | 一种基于三维正压浅海大陆架模型的深层海流反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810366515.1A CN108562896A (zh) | 2018-04-23 | 2018-04-23 | 一种基于三维正压浅海大陆架模型的深层海流反演方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108562896A true CN108562896A (zh) | 2018-09-21 |
Family
ID=63536365
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810366515.1A Pending CN108562896A (zh) | 2018-04-23 | 2018-04-23 | 一种基于三维正压浅海大陆架模型的深层海流反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108562896A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111859748A (zh) * | 2020-07-15 | 2020-10-30 | 天津大学 | 一种基于垂向混合坐标的海洋内波模拟方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103091711A (zh) * | 2013-01-24 | 2013-05-08 | 中国石油天然气集团公司 | 全波形反演方法及装置 |
CN106019394A (zh) * | 2016-04-27 | 2016-10-12 | 中国地质科学院矿产资源研究所 | 海洋大地电磁场非线性共轭梯度三维并行反演方法 |
CN107607936A (zh) * | 2017-08-31 | 2018-01-19 | 武汉大学 | 一种高频天地波雷达海洋表面流反演方法 |
-
2018
- 2018-04-23 CN CN201810366515.1A patent/CN108562896A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103091711A (zh) * | 2013-01-24 | 2013-05-08 | 中国石油天然气集团公司 | 全波形反演方法及装置 |
CN106019394A (zh) * | 2016-04-27 | 2016-10-12 | 中国地质科学院矿产资源研究所 | 海洋大地电磁场非线性共轭梯度三维并行反演方法 |
CN107607936A (zh) * | 2017-08-31 | 2018-01-19 | 武汉大学 | 一种高频天地波雷达海洋表面流反演方法 |
Non-Patent Citations (2)
Title |
---|
李自立: "地波雷达表面流探测与深层流反演算法研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
李自立等: "反演深层海流运动的差分方法研究", 《华中科技大学学报(自然科学版)》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111859748A (zh) * | 2020-07-15 | 2020-10-30 | 天津大学 | 一种基于垂向混合坐标的海洋内波模拟方法 |
CN111859748B (zh) * | 2020-07-15 | 2024-04-02 | 天津大学 | 一种基于垂向混合坐标的海洋内波模拟方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Walker | Physical and logistical considerations of using ultrasonic anemometers in aeolian sediment transport research | |
Walker et al. | Simulation and measurement of surface shear stress over isolated and closely spaced transverse dunes in a wind tunnel | |
Feddersen | Observations of the surf-zone turbulent dissipation rate | |
Walstra et al. | Modeling nearshore morphological evolution at seasonal scale | |
CN113124941B (zh) | 一种河道流量非接触式测量及精确计算方法 | |
CN107451372B (zh) | 一种运动波与动力波相结合的山区洪水过程数值模拟方法 | |
De Marchis et al. | Wind-and tide-induced currents in the Stagnone lagoon (Sicily) | |
CN104268848A (zh) | 一种海洋内波波速监测的方法 | |
Pidcock et al. | The spatial variability of vertical velocity in an Iceland basin eddy dipole | |
Davies et al. | An intercomparison and validation of a range of turbulence closure schemes used in three dimensional tidal models | |
Bjørnestad et al. | Lagrangian measurements of orbital velocities in the surf zone | |
Dey et al. | Reynolds stress anisotropy in flow over two-dimensional rigid dunes | |
CN108562896A (zh) | 一种基于三维正压浅海大陆架模型的深层海流反演方法 | |
Blenkinsopp et al. | Validation of volume continuity method for estimation of cross-shore swash flow velocity | |
CN110412312A (zh) | 涌潮流速垂向分布的一种计算方法 | |
Oetjen et al. | Experimental models of coarse-clast transport by tsunamis | |
Bernardes et al. | The alignment of the mean wind and stress vectors in the unstable surface layer | |
Sarkar et al. | Formation and migration of ripple pattern due to pure wave | |
Marchenko et al. | Asymmetric tide in Lake Vallunden (Spitsbergen) | |
Babanin et al. | Wind measurements near the surface of waves | |
Hinata et al. | Dependence of wind-driven current on wind stress direction in a small semienclosed, homogeneous rotating basin | |
Mukaro et al. | Wave height and wave velocity measurements in the vicinity of the break point in laboratory plunging waves | |
Bakhoday Paskyabi et al. | The influence of surface gravity waves on the injection of turbulence in the upper ocean | |
Uchiyama | Modeling three-dimensional cohesive sediment transport and associated morphological variation in estuarine intertidal mudflats | |
Lewis | Cross-shelf variations of near-inertial current oscillations |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180921 |