[发明专利]一种基于三维正压浅海大陆架模型的深层海流反演方法

摘要

一种基于三维正压浅海大陆架模型的深层海流反演方法，属高频地波雷达遥感探测技术领域，采用归一化的坐标系的三维非线性正压浅海大陆架模式，利用高频地波雷达探测到的海洋表面流数据，考虑复杂地形的边界问题，参数化物理过程，简化海洋控制方程组，再建立三维网格，采用差分法求解离散方程组进而反演得到相应的深层流数据信息，包括步骤：选择坐标系；导出归一化坐标系下的控制方程组；参数化物理过程；迭代法求解方程组得出深层流信息；实测数据与反演结果对比分析。本发明对比试验结果发现，实测数据曲线图和反演结果曲线图的整体趋势接近，反演曲线可以反映出海洋深层流变化的整体趋势，符合海洋动力学运动规律。

主权项

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是地球半径；是维度；是科氏参数；Ω是地转角速度；A_H是水平方向涡动粘滞系数，取常数5000m²/s；A_V是垂直方向涡动粘滞系数；ξ_x,ξ_y是在水平方向涡动粘滞系数中由于海底地形和水位有关的量，接近无穷小，忽略不计；由公式(02)得到垂直流速的表达式如下：在动量方程中，水平压强梯度力方程如下：在实际观测资料中发现，海水水平密度变化较小，所以在计算中可以忽略水平压强梯度力第二项，简化的水平压强梯度力方程如下：(3)参数化物理过程；海洋表面限制条件由高频地波雷达监测的表面流参数值决定，设定海洋表面垂直速度w＝0，动力学参数限制条件如下：其中τa表面风应力矢量，τb是海底摩擦力矢量，表达式如下所示：τa＝ρaCa|va|va                      (13)τb＝ρbCb|vb|vb                       (14)其中ρ_a是空气密度；v_a是风速矢量；C_a＝(0.73+0.069|v_a|)×10^‑3是海水对风的拖曳系数；ρ_b是海水密度；v_b是靠近海底的流速矢量；是海底摩擦应力拖曳系数；k是Kalman常数；Δσ是最后一层的厚度；z₀是海底粗糙度；根据Prandtl混合长度理论，得垂直方向涡动粘滞系数表达式如下：其中代表分子粘性系数；代表混合长度；s是反映海底粗糙度的参数，z'＝(1+σ)H代表考察点到海底的距离；(4)建立三维网格，离散化方程组；采用以下方式分别将u、v、w在x、y、σ方向上进行离散，以u为例：u的值为同一网格中四个结点的加权平均值，θ是空间加权系数，当θ＜0.5时，上述表达式为一阶，当θ＝0.5时，上述表达式为二阶，采用二阶精度，其格式如下：(5)迭代法求解方程组得出深层流信息；在公式(18)中，后一层结点上的值不仅与上一层结点值有关，还和本层已计算出的结点值有关，利用计算机，运用迭代法求解离散化的海洋控制方程组；(6)实测数据与反演结果对比分析；分析数值反演角度图发现，随着海流深度的增加，实测海流的流向和反演得出的海流的流向都呈逐渐右偏的趋势；分析数值反演速度图发现，随着海流深度的增加，实测海流的流速和反演得出的海流的流速都呈逐渐减小的趋势，在靠近海底处，由于底部摩擦力的作用，流速都趋近于零。