一种二维和三维嵌套海洋流体动力学数值模式及其在北部湾潮汐和潮流数值模拟中的应用

研究设计了一种二维和三维嵌套、外模态和内模态分离的水动力学数值模式,既可用于潮汐,也可用于风暴潮的数值计算。该模型对全部海区进行二维计算,对其中重点关心的海区同时进行三维计算。在三维计算区域,采用了内、外模态既分离又耦合的计算技术。数值格式采用全部交错的网格结构,三维模型中垂直方向采用σ-坐标代替通常的z-坐标,垂直涡动粘性系数由混合长度理论确定,垂直粘性项采用隐式差分格式。作为算例,本文对南海北部湾潮波进行了细网格的数值计算。     

A quasi-three-dimensional numerical prediction model of salinity structure in Bohai Sea and Huanghai Sea

Ａｑｕａｓｉ－ｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌｎｕｍｅｒｉｃａｌｐｒｅｄｉｃｔｉｏｎｍｏｄｅｌｏｆｓａｌｉｎｉｔｙｓｔｒｕｃｔｕｒｅｉｎＢｏｈａｉＳｅａａｎｄＨｕａｎｇｈａｉＳｅａ￥ＳｕｎＷｅｉｙａｎｇａｎｄＷａｎｇＺｏｎｇｓｈａｎ（Ｒｅｃｅｉｖｅ...     

Three-dimensional structure of tidal current in the East China Sea and the Yellow Sea

A three-dimensional tidal current model is developed and applied to the East China Sea (ECS), the Yellow Sea and the Bohai Sea. The model well reproduces the major four tides, namely M₂, S₂, K₁ and O₁ tides, and their currents. The horizontal distributions of the major four tidal currents are the same as those calculated by the horizontal two-dimensional models. With its high resolutions in the horizontal (12.5 km) and the vertical (20 layers), the model is used to investigate the vertical distributions of tidal current. Four vertical eddy viscosity models are used in the numerical experiments. As the tidal current becomes strong, its vertical shear becomes large and its vertical profile becomes sensitive to the vertical eddy viscosity. As a conclusion, the HU (a) model (Davieset al., 1997), which relates the vertical eddy viscosity to the water depth and depth mean velocity, gives the closest results to the observed data. The reproduction of the amphidromic point of M₂ tide in Liaodong Bay is discussed and it is concluded that it depends on the bottom friction stress. The model reproduces a unique vertical profile of tidal current in the Yellow Sea, which is also found in the observed data. The reason for the reproduction of such a unique profile in the model is investigated.     

Numerical Study on Tidal Mixing in the Bohai Sea

The spatial and temporal variability of tidal mixing in Bohai Sea is studied using a numerical approach. In calculating tidal mixing, accurate barotropic tidal current is obtained via a harmonic analysis package utilizing the simulated current output from a high-resolution regional ocean model. And a “small-scale” roughness map is adopted to describe the detailed topographic features of Bohai Sea. It is shown that the tidal mixing estimated in Bohai Sea is much higher than the level of global background, and fluctuates considerably at some regions within a single day. In Liaodong Bay, Bohai Bay and Bohai Strait, the mixing varies greatly, with the peak value of O (10^?2) m² s^?1. The order of magnitude of mixing in Laizhou Bay is about O (10^?5～10^?3) m² s^?1. Mixing with background level of O (10^?5) m² s^?1 only appears in central area. Result also shows that rough topography plays relatively a more important role than tidal current in enhancing diapycnal mixing in Bohai Sea. The distributions of tidal mixing in selected sections reveal that the vertical stratification in Bohai Sea is not obvious, generally renders a barotropic structure.     

渤海潮流及湖余流的数值计算

渤海是我国重点保护水系之一。为了控制水域污染,制定废水排海标准,首要步骤须查明入海污染物质的扩散过程。在渤海,潮流是形成污染扩散的主要营力之一,为此首先计算渤海的潮流及潮余流场。近二十年来,随着电子计算机的迅速发展和计算技术的进步,潮汐动力学进入了数值研究的新阶段。这一新方向的特点,是以求解自然条件下长波方程的初值——边值问题为出发点,研究以往所不能解决的许多复杂现象。     

Vertical structure of the M<Subscript>2</Subscript> tidal current in the Yellow Sea

The vertical structure of the M₂ tidal current in the Yellow Sea is analyzed from data acquired using an acoustic Doppler current profiler. The observed vertical profiles of the M₂ tidal current are decomposed into two rotating components of counter-clockwise and clockwise, and restructured using a simple one-point model with a constant vertical eddy viscosity. The analyzed results show that the internal fictional effect dominates the vertical structure of the tidal current in the bottom boundary layer. In the Yellow Sea, the effect of the bottom friction reduces the current speed by about 20–40% and induces the bottom phase advance by about 15–50 minutes. In the shallower coastal regions, the effects of bottom topography are more prominent on the vertical structure of tidal currents. The vertical profile of the tidal current in summer, when the water column is strongly stratified, is disturbed near the pycnocline layer. The stratification significantly influences the vertical shear and distinct seasonal variation of the tidal current.     

A Three-dimensional Tidal Model in Boundary-fitted Curvilinear Grids

In hydrodynamic models of marine and estuarine currents the use of boundary-fitted curvilinear grids not only makes the model grids fit to the coastline and bathymetry well, but also makes the kinetic boundary conditions simple and more accurate. Because of these advantages, a three-dimensional tidal model with boundary-fitted curvilinear grids has been developed to simulate both tide and current in estuarine and shelf water. The basic idea is to use a set of coupled σ-stretched and elliptic transformations to map the physical space into a corresponding transformed space such that all boundaries are coincident with co-ordinate lines and the transformed grids are rectangular. The hydrodynamic equations in the transformed space are solved in a rectangular mesh of the transformed grid system. The application of the model in the simulation of the M₂,S₂ , K₁and O₁tidal waves in the Bohai Sea, China, shows that the numerical results are in good agreement with the observations. The method of boundary-fitted curvilinear grids is effective for improving simulation accuracy of current in the estuarine and shallow seas, especially at the coastal regions where the current was usually impractical by a uniform Cartesian grid system with the shoreline and bathymetry represented by numerous stair-steps.     

Three-dimensional Layer-integrated Modelling of Estuarine Flows with Flooding and Drying

Details are given of the refinement and application of a thee-dimensional (3-D) layer-integrated numerical model of tidal circulation, with the aim of simulating severe tidal conditions for practical engineering applications. The mode splitting strategy has been used in the model. A set of depth-integrated 2-D equations are first solved to give the pressure gradient, and the layer-integrated 3-D equations are then solved to obtain the vertical distributions of the flow velocities. Attention has been given to maintaining consistency of the physical quantities derived from the 2-D and 3-D equations. A TWO=layer mixing length turbulence model for the vertical shear stress distribution has been included in the model. Emphasis has been focused on applying the model to a real estuary, which is geometrically complicated and has large tidal ranges giving rise to extensive flooding and drying. The model has been applied to three examples, including: wind-driven flow in a rectangular lake, tidal circulation in a model rectangular harbour, and tidal circulation in a large estuary. Favourable results have been obtained for both the simple and complex flow beds.     

A three-dimentional model of the northeast East China Sea, with application to the M2 current field

A simple three-dimensional tidal model is used to examine the M2 tidal current distribution in a northeastern part of the East China Sea, especially the vertical variation of the current in the region. Computed M2 current is compared with observations available and found to be in good agreement.Main features of the calculating method in this study are: (1) Vertical variation of the tidal current is taken as a funetion of the depth-mean velocity: (2) the method is applicable to a variety of the vertical eddy viscosities; (3) it has a fine vertical resolution, especially near the sea bootom. So, this method not only enables us to get a steady state solution easily but also depicts effects of the friction on the vertical variation of the current much better.     

Numerical modelling of Lagrangian residual current in the Bohai Strait

The Bohai Strait is the only passage-way of the semienclosed Bohai Sea into the outer ocean. The various dynamical processes have great effects on both sides of the strait, especially on the Bohai Sea. Tidal current dominates in the strait, and among the tida! components the M, component tide is prevailing.The seawater in the strait is homogeneous vertically almost the whole year because of strong tidal mixing. The M, component tide was simulated numerically using two dimensional barotropical model, and the Eulerian tide-induced residual currents were also deduced. The results showed that the tidal currents were dominated by the coastal line and topography. There exist many eddies in the tide-induced residual current fields, particularly near the northern coast.On the basis of the simulation of the tidal currents, the movements of the water particles in the strait were tracked by means of Lagrangian techniques. Consequently the Lagrangian residual currents were obtained. Somewhat similar to the Eulerian     

A three-dimensional semi-implicit unstructured grid finite volume ocean model

A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C -unstructured non-orthogonal grid in the horizontal direction and z -level grid in the vertical direction. The three-dimensional model is discretized by the semi-implicit finite-volume method, in that the free-surface and the vertical diffusion are semi-implicit, thereby removing stability limitations associated with the surface gravity wave and vertical diffusion terms. The remaining terms in the momentum equations are discretized explicitly by an integral method. The partial cell method is used for resolving topography, which enables the model to better represent irregular topography. The model has been tested against analytical cases for wind and tidal oscillation circulation, and is applied to simulating the tidal flow in the Bohai Sea. The results are in good agreement both with the analytical solutions and measurement results.     

秦皇岛海域海流特征及规模化养殖对其影响的观测研究

秦皇岛海域是辽东湾与渤海中部及渤海湾进行物质和能量交换的重要通道。本文基于海床基观测平台获取的夏秋季海流连续观测资料，运用调和分析和滤波等方法对该海域的海流特征及其对规模化养殖的响应进行了研究。结果表明:秦皇岛海域最显著的潮流是M₂分潮流，其最大流速介于20.0~36.9 cm/s之间，远小于辽东湾东部海域M₂分潮流最大流速；秋季秦皇岛海域余流流速介于0.2~2.5 cm/s之间，整体上较辽东湾东侧海域余流弱，辽东湾底层可能存在逆时针的弱环流系统；夏季秦皇岛海域M₂和K₁分潮流的最大流速均大于秋季；养殖活动对余流影响较大，养殖区中部A7、A8站余流的垂向平均流速比养殖区边缘A6站分别减小76%和18%左右。     

Three-dimensional modeling of tidal circulation and mixing over the Java Sea

A combination of a three-dimensional hydrodynamic model and in-situ measurements provides the structures of barotropic tides, tidal circulation and their relationship with turbulent mixing in the Java Sea, which allow us to understand the impact of the tides on material distribution. The model retains high horizontal and vertical resolutions and is forced by the boundary conditions taken from a global model. The measurements are composed of the sea level at coastal stations and currents at moorings embedded in Seawatch buoys, in addition to hydrographic data. The simulated tidal elevations are in good agreement with the data for the K₁ and M₂ constituents. The K₁ tide clearly shows the lowest mode resonance in the Java Sea with intensification around the nodal point in the central region. The M₂ tide is secondary and propagates westward from the eastern open boundary, along with a counterclockwise amphidromic point in the western part. The K₁ tide produces a major component of tidal energy, which flows westward and dissipates through the node region near the Karimata Strait. Meanwhile, the M₂ tide dissipates in the entire Java Sea. However, the residual currents are mainly induced by the M₂ tide, which flows westward following the M₂ tidal wave propagation. The tidal mixing is mainly caused by K₁ tide which peaks at the central region and is consistent with the uniform temperature and salinity along the vertical dimension. This mixing is expected to play an important role in the vertical exchange of nutrients and control of biological productivity.     

渤海西南部海域潮流数值计算

采用二维有限差分解法对渤海西南部海域进行潮流数值模拟,得出M_2分潮的潮波图,椭圆长短轴图,最大流速分布和不同时刻的潮流场图,并根据实测资料对其进行了验证。从而可系统地了解这一海区M_2分潮的潮波系统及逐时潮流场等海洋要素的分布情况。     

中国陆架潮流沉积体系和模式

在1991年中法合作渤海潮流沉积研究的基础上,查阅了国内外有关研究成果,分析了中国陆架的水深地形、沉积地貌与潮流动力的关系,认为潮流对中国陆架的海底地貌和沉积的形成发育起了主导作用。当潮流流速大于3节时,潮流的侵蚀作用是主要的,往复潮流多形成冲刷深槽,大大刷深了海峡或水道。当潮流流速1－3节时,潮流的沉积作用是主要的,多形成浅滩,即潮流沙脊和潮流沙席。以M2分潮椭率绝对值0．4为界,大于0．久者意味着潮流旋转性强,多形成潮流沙席;小于0．4者意味着潮流往复性强,多形成潮流沙脊。提出了我国邻近陆架发育了5个现代潮流沉积地貌体系：（1）黄海东部潮流沉积体系,它由西朝鲜湾潮流沙脊和其南部的沙席两者组成;（2）渤海东部潮流沉积体系,它由老铁山水道冲刷槽,辽东浅滩沙脊和渤中浅滩沙席三者组成;（3）长江口外潮流沉积体系,它由江苏滨外潮流沙脊和长江口浅滩潮流沙席组成;（4）台湾滨外潮流沉积体系,它由台湾海峡冲刷槽、台湾浅滩沙脊、澎湖水道冲刷槽和台中浅滩沙席四者组成;（5）琼州海峡潮流沉积体系,它由琼州海峡冲刷槽、东浅滩沙脊和西浅滩沙脊三者组成。此外,在东海陆架上还有冰后期海侵早期形成的残留潮流沉积体系。全新世陆架浅海潮流沉积模式可分海峡一浅     

风暴潮预报模式在渤海海域中的应用研究

首先运用COHEREN S三维多功能大陆架水动力学模型,对渤海四个主要分潮进行了调和分析,建立了渤海天文潮预报模型。结合9712号热带风暴的路径记录资料,运用风场、气压场估算公式,估算了风暴所经地带风场、气压场的变化过程。以此为输入条件,运用COHEREN S模型模拟了9712号热带风暴作用下渤海的增水过程,还模拟讨论了底摩阻系数对水位的影响。     

渤海开边界潮汐的伴随法反演

潮汐潮流数值模拟中的一个主要难点在于开边界条件的确定。本文采用伴随法 ,由渤海沿岸 1 9个验潮站的潮汐调和常数来反演渤海海域的开边界条件 ,以实现渤海潮波的数值模拟。计算所得调和常数与实测值之差的绝对平均值 :m1 潮波振幅差为 1 4cm ,迟角差为5 0°;M2 潮波振幅差为 2 4cm ,迟角差为 5 0°。数值模拟结果正确地反映了渤海m1 和M2 潮波的基本特征     

渤海潮波系统的数值模拟

利用二维非线性潮波方程组,讨论了渤黄海主要分潮(全日潮、半日潮及浅水分潮) 数值模拟中的有关问题。数值模拟中同时考虑了4个主要分潮(M2,S2,K1,O1)和两个浅水分潮(M4,MS4)。分析表明,在渤黄海潮波系统数值模拟中,稳定后选取14 d的数值模拟结果进行调和分析能够取得最佳(最合理)的调和分析结果。计算出调和常数的模拟值与实测值之差的绝对平均值:M2分潮的振幅差为4cm,迟角差为3.3°,S2分潮的振幅差为2cm,迟角差为4.2°,K1 分潮的振幅差为1cm,迟角差为3.7°,O1分潮的振幅差为2 cm,迟角差为5.5°。实验结果较好地体现了渤黄海潮波系统的特征。     

ADCIRC模式在渤海M2分潮模拟中的应用研究

利用有限元方法的ADCIRC (Advanced Circulation Model)海洋模式,建立了渤海高分辨率的二维潮汐潮流模型,模式结果与实测资料吻合良好.模式成功模拟出了M2分潮在渤海的2个无潮点和3个圆流点,位置与前人的研究结果基本一致.M2分潮流在渤海中央为顺时针旋转的旋转流,在辽东湾、渤海湾和莱州湾基本为...     

南海北部陆架陆坡区海流观测研究

针对2006-2009年期间,南海北部陆架陆坡区3个站ADCP海流连续观测资料,采用功率谱分析、潮流调和分析方法,重点分析了陆架陆坡区100 m,200 m和1 200 m水深海域海流的垂向结构,探讨了环流的季节变化和空间分布特征,特别讨论了南海暖流和北陆坡流的时空变化特征。结果表明,陆架陆坡区潮流类型属于不规则日潮,深水站点中层表现为正规全日潮类型,垂向为"三层结构",甚至更加复杂。O₁,K₁,M₂,S₂等分潮总体上为顺时针旋转,在深水站点,基本表现为西北－东南走向的往复流形态。从能量角度看,表层和底层海流中,潮流所占份额较大,分别占30%~40%和40%~50%,中层较小,约为20%。对东沙群岛西南陆架陆坡区环流,观测计算结果证实了西向强流的存在,且垂向结构具有显著的季节变化,在200 m水深处没有明显的南海暖流,只是10~30 m以上层次存在逆风海流。海南岛以东海域连续15个月表层环流的结果表明,冬季明显受到南海暖流的影响,存在东北向的逆风海流,夏秋季的环流表现为西南向,流速较强,夏季也存在逆风情况,造成上述情形的原因可能是该地南海暖流的流轴具有季节性变化——冬季偏南,夏季偏北。