舟山绿色石化三维潮流特性数值研究

本文基于MIKE3平台建立宁波、舟山及其附近海域三维水动力模型,通过与实测数据进行潮流验证的基础上,对比了三维模型与二维模型的区别,并主要对舟山绿色石化三期围垦工程实施前后周边潮流场的变化影响进行了预测分析。结果表明:（1）同样的初始条件下,二维模型和三维模型的潮位和潮流验证结果都是良好的,而在某些情况下二维数值结果更接近于实测值;（2）三维垂向上的流速分布呈现出随水深增加而减小的趋势,总体上潮流拟合效果较好,其中0.6H水深以下的流速相对误差结果在10%～15%之间,考虑模型可能受未考虑风场和波浪影响的缘故,局部测站表层计算流速较实测值稍大;而垂向上的流向分布近乎一致;（3）围垦工程前后并未对工程周边潮流场有明显影响,仅在大鱼山岛东北侧、东南侧和西北角受“截弯取直”和小岛屿间水道围堵的影响,有局部潮流形态的改变。整体上看,以顺应岸线方向上的围垦工程布置方案会对局部近岸潮流产生影响,对外海大范围潮流场的影响几乎不计。     

The morphodynamic modelling of tidal sand waves on the shoreface

An artificial sand wave on the Dutch shoreface of the North Sea has been studied in conditions with relatively strong tidal currents in the range of 0.5 to 1 m/s and sediments in the medium sand size range of 0.2 to 0.5 mm. The sand wave is perpendicular to the tidal current and has a maximum height and length of the order of 5 m and 1 km, respectively. The sand wave is dynamically active and shows migration rates of the order of a few metres per year. A numerical morphodynamic model (DELFT3D model) has been used to simulate the morphological behaviour of the sand wave in the North Sea. This model approach is based on the numerical solution of the three-dimensional shallow water equations in combination with a surface wave propagation model (wind waves) and the advection–diffusion equation for the sediment particles with online bed updating after each time step. The model results show that the sand wave grows in the case of dominant bed-load transport (weak tidal currents; relatively coarse sediment; small roughness height; low waves) and that the sand wave decays in the case of dominant suspended transport (strong currents, relatively fine sediment, large roughness height; storm waves).     

Numerical modeling of long waves in shallow water using Incremental Differential Quadrature Method

Incremental Differential Quadrature Method (IDQM) as a rapid and accurate method for numerical simulation of Nonlinear Shallow Water (NLSW) waves is employed. To the best of authors’ knowledge, this is the first endeavor to exploit DQM in coastal hydraulics. The one-dimensional NLSW equations and related boundary conditions are discretized in space and temporal directions by DQM rules and the resulting system of equations are used to compute the state variables in the entire computational domain. It was found that the splitting of total simulation time into a number of smaller time increments, could significantly enhance the performance of the proposed method. Furthermore, results of this study show two main advantages for IDQM compared with other conventional methods, namely; unconditional stability and minimal computational effort. Indeed, using IDQM, one can use a few grid points (in spatial or time direction) without imposing any stability condition on the time step to obtain an accurate convergent solution.     

Shifts in the timing of spawning in sole linked to warming sea temperatures

Phenotypic traits such as peak spawning time may vary within and differ between populations in relation to environmental factors, such as temperature. Sole (Solea solea) is a valuable, commercially exploited species that spawns in late winter or spring. The date of peak spawning was estimated for each year for seven stocks from monthly fish samples collected from commercial fisheries since 1970. Four out of seven stocks showed a significant long-term trend towards earlier spawning (Irish Sea, east-central North Sea, southern North Sea, eastern English Channel) at a rate of 1.5 weeks per decade. The other three stocks (Bristol Channel, western English Channel and western-central North Sea) failed to show a relationship, but the available time series were limited for these stocks (< 10 years). Sea surface temperature during winter significantly affected the date of peak spawning, although the effect differed between stocks. The implications of the effect of winter temperature on the timing of spawning for the population dynamics are discussed.     

Comparison of two three-dimensional hydrodynamic modeling systems for coastal tidal motion

A comparison of two three-dimensional numerical modeling systems for tidal elevations and velocities in the coastal waters is presented. The two modeling systems are: (1) the Princeton Ocean Model (POM) and (2) the MIKE 3 flow model. The model performance results for Singapore's coastal waters show that the predicted tidal elevations from the two hydrodynamic modeling systems are almost identical and are in very good agreement with field measurement data. The simulated tidal current velocities match well with field measurement data at the selected stations, but it seems that the POM provides the slightly better simulation, compared to the MIKE 3 flow model. The depth profiles of the velocities obtained from the two modeling systems may be greatly different at some time, due to the vertical diffusion coefficient calculated from different turbulent sub-models in the two modeling systems. The POM generally predicts larger peak tidal velocities. The maximum speed differences for the model results from the two modeling systems occur in the top and differ from time to time and from location to location, reaching up to 20%.     

Modelling 3D hydrodynamics governing island-associated sandbanks in a proposed tidal stream energy site

A 3D numerical modelling study to investigate the existing hydrodynamic regime of the Pentland Firth Inner Sound Channel, Scotland, UK is presented. Hydrodynamics that govern some sensitive sedimentary deposits in the Inner Sound Channel are discussed. A 3D hydrodynamic model Delft3D is set up for Pentland Firth, located between Orkney Islands and mainland Scotland and a full sensitivity analysis of the numerical model is carried out. The current model set up sufficiently captures the existing hydrodynamics during a full spring-neap tidal cycle inside Pentland Firth. Using model results, the 3D structure of the dynamics of the tidal flows in the Inner Sound Channel is investigated. The temporal variability of tidal flows, the residual tidal flows in the channel and local flow interactions with the Island of Stroma are described. It is proved that the tidally dominant flows drive the sediment transport gradient model to explain the principle maintenance mechanisms of two island-associated sandbanks present in the Inner Sound. The present study provides detailed information on the physics of the tidal regime in the Inner Sound and explains the presence of sandbanks in an area of high tidal flows. Due to extremely high tidal flows, Inner Sound is considered as one of the most favourable sites for tidal energy extraction in the UK. The findings of this study will be very useful in assessing the significance of impacts of future tidal energy extraction on natural hydrodynamics and sediment dynamics of the area.     

围填海工程对半封闭海湾水动力环境的影响

为研究围填海工程对海湾水动力环境的影响,本文基于二维数值模型MIKE21,建立了东山湾附近海域的潮流模型。对比观测数据发现,大潮期间的最高、最低潮位模拟误差在6 cm以内,小潮期间的误差相对较大;流速和流向的模拟误差在9%左右,最大误差出现在转流时刻;总体来看,模拟结果与观测数据吻合良好。在此模型基础上,研究了东山湾围填海前、后潮流动力、水体半交换时间和纳潮量等水动力要素的变化。结果表明:围填海后大范围的涨落流态与围填海前保持一致,大体上仍然呈现S-N走向,涨潮流偏N向,落潮流偏S向,往复流特征较为明显;从局部流场来看,涨落潮流场发生了一定的变化,围填海区域南、北两侧的流矢变化较为明显,涨潮流矢由偏N向改为偏E向,而落潮流矢由偏S向改为偏W向,同时受到围填海区域岸线的遮蔽效应,围填海南、北两侧水域的流速也有一定减弱,而西侧的流速则有一定增强;围填海实施前的水体半交换时间为220.5h,实施后时间为239.4 h;纳潮量变化为–2.5%左右。研究表明,围填海工程对东山湾水动力环境的影响主要集中在工程区域附近,其对泥沙冲淤、生态环境等的影响将在后续研究中进一步探讨。     

渤、黄、东海8个主要分潮的数值模拟研究

应用MIKE数值模拟软件,采用无结构三角形网格,建立一套计算区域包括整个渤海、黄海、东海以及东海大陆架和琉球群岛的高分辨率数值模型,考虑了实际水深和岸线,外海开边界采用西北太平洋大模型结果的潮位提供,模拟了东中国海潮波的波动过程,对潮波垂直运动过程进行调和分析,得到了渤海、黄海、东海的M2,S2,K1,O1以及N2,K2,P1,Q1八个主要分潮的传播和分布特征。利用中国沿海14个潮位站的调和常数对模型结果进行了验证,验证结果显示模型较为准确可靠。研究结果表明:4个主要半日潮(全日潮)在渤、黄、东海的传播情形基本相似,即潮波在渤海、黄海、东海沿岸的传播性质上类似沿岸开尔文波的传播形态,并且成功再现了计算海域的4个半日分潮无潮点和2个全日分潮无潮点。全日潮振幅各无潮点附近振幅最小,而海湾的波腹区振幅最大,东海潮差呈现近岸方向振幅大、离岸方向振幅小,浙闽沿海振幅也较大,黄海振幅相对较小,渤海振幅在辽东湾和渤海湾顶最大,两个无潮点周边振幅较小。     

基于动力条件改变的笔架山连岛坝侵蚀原因分析

利用MIKE21软件的波浪谱模型及潮流模型对辽宁锦州大笔架山岛连岛坝海域的波浪、潮流等水动力条件进行了数值计算,并对锦州港建设前后的动力条件变化进行了分析.结果表明：锦州港的建设完全阻断了偏西向波浪对连岛坝的塑造作用,坝体的侵蚀确系单侧（偏东向）波浪作用的结果,坝体局部侵蚀严重区域形成四处明显的潮沟,潮沟处流速增大,加剧了对坝体的侵蚀.     

长江口新桥水道潮汐不对称性研究

涨潮槽是河口重要的地貌类型,其潮汐动力受地形摩擦、径流及风浪等因素影响而发生不对称现象,直接影响槽内的泥沙动力过程。本文基于实测数据与数值模型计算潮汐不对称偏度,研究长江口南支最大的涨潮槽——新桥水道的潮汐不对称性及其影响因素。结果表明：（1）受余流与各分潮的相互作用控制,新桥水道可分为三个区域,上段为涨潮优势区域,中段为涨落潮转换区域,下段是落潮优势区域。（2）新桥水道欧拉余流表现出明显的分段差异,上段指向陆而在中下段指向海并受到径流影响。斯托克斯余流则均指向陆并向陆增强。（3）新桥水道内中下段潮汐不对称受地形的影响较大,其中扁担沙的北移增强了新桥水道中段的涨潮优势特性,新的新桥通道的形成促进了新桥水道下段的落潮优势特性。     

Numerical and Experimental Study of the 3D Effect on Connecting Arm of Vertical Axis Tidal Current Turbine

Vertical axis tidal current turbine is a promising device to extract energy from ocean current. One of the important components of the turbine is the connecting arm, which can bring about a significant effect on the pressure distribution along the span of the turbine blade, herein we call it 3D effect. However, so far the effect is rarely reported in the research, moreover, in numerical simulation. In the present study, a 3D numerical model of the turbine with the connecting arm was developed by using FLUENT software compiling the UDF (User Defined Function) command. The simulation results show that the pressure distribution along the span of blade with the connecting arm model is significantly different from those without the connecting arm. To facilitate the validation of numerical model, the laboratory experiment has been carried out by using three different types of NACA aerofoil connecting arm and circle section connecting arm. And results show that the turbine with NACA0012 connecting arm has the best start-up performance which is 0.346 m/s and the peak point of power conversion coefficient is around 0.33. A further study has been performed and a conclusion is drawn that the aerofoil and thickness of connecting arm are the most important factors on the power conversion coefficient of the vertical axis tidal current turbine.     

Numerical study of power production from tidal energy in the Khuran Channel and its feedback on background hydrodynamics

This study focuses on the development of a farm of tidal turbines in the Khuran Channel. The important factors include the location of turbines and their hydrodynamic effects on the environment. A three-dimensional circulation model for the Persian Gulf is employed for the comprehensive evaluation of the tidal energy potential throughout the study area. The model is validated by using in situ observations of water level and current data.The appropriate potential points for extracting the tidal e...     

沙埕港湾口断面潮流及余流特征分析

基于对沙埕港湾口断面的连续走航观测资料,成功构建了沿走航断面的10个站点的连续海流序列,并分析了潮流、余流、潮通量等水文要素。分析结果表明,沙埕港湾口水道潮流类型为正规半日潮流,涨潮最先出现在中下层而落潮最先出现在上层,涨(落)潮转流相差约为30min。水道内潮流为往复流,M2和S2分潮流流速较大,倾角基本沿水道主轴方向。沙埕港湾口断面余流呈2层结构,10m以浅基本为东南向余流流出湾口,核心位于湾口断面南侧。10m以深多为西北向流入湾内,入流核心位于湾口断面中部的底层区域。对潮通量的计算表明,通过湾口进入沙埕港的潮通量约为1.63×108m3。     

环台湾岛海域全日分潮的特征和潮汐、潮流的综合性质

用97版POC海洋模式，对环台湾岛海域的全日分潮和整个潮汐，潮流综合特征进行三维数值研究。研究海域全日分潮是由太平洋传入的，且来自台湾岛北部海区传入的潮波穿越海峡。由吕宋海峡传入的全日分潮对维持南海的潮运动起着重要作用。全日分潮最大流同时线分布表层有5个圆流点，其中4个本文首次得到。台湾海峡及其以北海区和台湾东部洋区为不规则半日潮区，台湾东南为全日分潮为主的混合潮区。台湾岛北部为气旋式余流涡旋区，环绕台浅滩为反气旋余流涡旋区，澎湖水道开始的转向流预示着海峡及其邻近海区的涨潮流或落潮流的来临。     

基于FVCOM的南海北部海域潮汐潮流数值模拟

基于非结构三角形网格的FVCOM(finite-volume coastal ocean model )数值模型, 对南海北部海域的潮汐、潮流进行了精细化数值模拟研究, 并根据模拟结果详细分析了M2, S2, K1, O1 分潮的潮汐和潮流特征。研究结果表明: 神泉港到甲子港海域表现为正规全日潮性质, 珠江口附近海区潮汐以不正规半日潮为主, 其他海域主要表现为不规则全日潮; 陆架海域和深水海域主要表现为往复流, 陆架坡折区存在较强的旋转流, 陆架坡折区为不规则半日潮流和不规则全日潮流的分界线; 东沙群岛附近海域以不规则全日潮流为主, 旋转方向为顺时针; 整个海域的最大流速分布与等深线基本平行, 东沙群岛附近速度明显变大, 最大值出现在台湾浅滩附近, 最大值超过70 cm/s; 南海潮波系统以巴士海峡传入的大洋潮波为主, 分为三支潮流, 以不同的形式进出南海北部海域; 余流在台湾浅滩附近达到最大, 超过6 cm/s, 自南向北进入台湾海峡, 近岸余流自东向西沿岸流动。本研究在东沙群岛周边的模拟结果与前人基于实测资料的分析吻合较好, 并且由于采用了高精度的三角网格, 本文对东沙群岛周边海域的潮汐潮流结构和性质的刻画和分析是迄今为止较为精细的, 同时本研究还提高了对沿岸验潮站调和常数的模拟精度。     

Favorable conditions for cold-water intrusion from the Kuroshio intermediate layer into Osaka Bay

From 1980 to 1995, in August, the bottom layer of Osaka Bay was occupied by cold, nutrient-rich water compared with that observed during both previous and subsequent decades. To investigate the mechanisms for the intrusion of bottom-layer cold water into Osaka Bay, the intrusion into Osaka Bay via the Kii Channel is simulated using a finite-volume coastal ocean model with unstructured triangular cell grids. The initial conditions, boundary conditions, and surface temperature given to the model are obtained from daily reanalysis data provided by the Japan Coastal Ocean Predictability Experiment. The model shows that cold water uplifted on the eastern side of the Kii Peninsula is propagated westward at 1.0 m/s as a coastal boundary current; it reaches the Kii Channel mouth when the Kuroshio axis is located around 74 km south of Cape Shionomisaki. However, the modeled cold water mass at the Kii Channel mouth does not intrude further to the north of the Kii Channel; therefore, another mechanism is required to explain the cold-water intrusion into the bottom layer of Osaka Bay. A plausible explanation is the estuarine circulation established by the freshwater supply at the bay head. When the river runoff is included in the model without forced vertical mixing, the temperature in Osaka Bay decreases 6.6 days later than the temperature decreases at the Kii Channel mouth. Furthermore, the shoreward current speed in the bottom layer of the modeled estuarine circulation is 15 cm/s, which provides the mechanism required for the cold water mass to pass the Kii Channel.     

The influence of tides,wind and waves on the redistribution of nourished sediment at Terschelling,The Netherlands

A calibrated morphodynamic model of the barrier island of Terschelling, The Netherlands [Grunnet, N.M., Walstra, D.J.R., Ruessink, B.G., 2004. Process-based modelling of a shoreface nourishment. Coastal Eng. 51/7, 581–607], comprising the island and its two adjacent tidal inlet systems, is applied to identify the relative contribution of tides, wind and waves to the cross-shore and alongshore redistribution of a 2 Mm³ nourishment supplied to the nearshore zone along the island. Several model simulations with varying combinations of horizontal and vertical tide, wind and wave forcing were designed to investigate the effect of each individual forcing on a large spatio-temporal scale (order of kilometres and months, respectively). As expected, stirring and transport by waves and wave-induced currents are predicted to be by far the dominant contributor to the net sediment transport along the coast of Terschelling. Because of the strong obliquity of the winds and the relatively small tidal currents in front of the island (≈ 0.5 m/s), alongshore wind-driven currents increase sediment transport rates and horizontal tides virtually have no net transport capacity. This motivated a local model of the study area along the closed coast of Terschelling, not including tidal inlets and further simplifying tidal boundary definitions by omitting the horizontal tides: morphodynamic simulations of the local model show virtually identical results as the larger model predictions. The reduction in complexity in setting up a local model instead of a regional model coupled with the corresponding significant reduction in computational time points to an increasing applicability of complex process-based models.     

基于船载ADCP观测对罗源湾湾口断面潮流及余流的分析

基于对罗源湾可门水道的25 h连续走航ADCP观测,成功构建了沿走航断面共12个站位的连续海流时间序列,并对这些站位的潮流、余流以及潮通量等进行了分析。结果表明可门水道内的潮流为正规半日潮流,驻波性质明显,涨潮首先出现在水道中下层而退潮则首先发生在水道上层。水道内的潮流为往复流,水道南部M₂分潮流流速较大,并且其倾角自北向南逐渐增加。此外,水道两端的浅水区域内浅水分潮M₄振幅较显著。可门水道内余流呈现出两层结构,20 m以浅余流沿东北向流出海湾,并且出流的核心位置偏南,而20 m以深的余流沿西南向流入湾内,入流的流核位于偏北的近底层区域。对潮通量的积分计算表明通过可门水道进入罗源湾的潮通量约为4.81×108 m3。     

Nearshore directional wave measurements by surface-following buoy and acoustic Doppler current profiler

Directional energy spectra of nearshore surface waves were measured for a 3-year period (2004–2007) at a site with mean depth 14 m and mean tidal range 2.1 m. Triaxys surface-following wave buoys reported hourly directional wave energy spectra and wave parameters near the offshore end of the Savannah River Entrance Channel, Georgia, USA. An acoustic Doppler current profiler (ADCP) was located beside the wave buoy for 3 months. Directional and non-directional surface wave energy spectra and the corresponding bulk wave parameters (height, period, and direction) are compared for the two systems. Most parameters derived from the spectra agree closely; the most significant differences were found at the upper and lower frequency measurement limits, where signal-to-noise ratios were lower. The wave buoy consistently reports a small amount of energy below 0.05 Hz that does not appear in the ADCP-derived spectra and does not appear to be related to the mooring system. This leads to larger mean and peak periods reported by the buoy. All directional spectra were computed using the Maximum Entropy Method for both instruments, but the buoy, with spectra derived from six independent time series, provides lower directional resolving power than the ADCP, which utilizes twelve time series. Both systems gave similar results defining mean and peak wave directions, with the primary difference being that the ADCP indicates energy to be more tightly concentrated around the peak direction.     

Numerical simulation of tidally induced eddies in the Bungo Channel: A possible role for sporadic Kuroshio-water intrusion (kyucho)

It is well known that the sudden intrusion of Kuroshio warm water into the Bungo Channel (kyucho) is regulated by spring–neap tidal forcing. In order to clarify the physical background behind this regulation, numerical experiments are carried out using a high-resolution non-hydrostatic three-dimensional model. We first reproduce the strong mixing region off the east coast of the Bungo Channel resulting from tidal flow interaction with complicated land configurations during spring tides; behind islands and headlands, small-scale eddies satisfying an approximate cyclostrophic balance are generated. As a result, averaged over the whole model domain, the tidal-mean energy dissipation rate reaches ≈1.6?×?10^?6?W?kg^?1. The model predicted energy dissipation rates at the location and times of direct microstructure measurements in the Bungo Channel are comparable to the observed values. We next examine whether or not strong tidal mixing thus reproduced can inhibit the northward intrusion of Kuroshio warm water in the Bungo Channel. It is shown that the Kuroshio warm water can (or cannot) pass through the tidal mixing regions off the east coast of the Bungo Channel during periods of weakened (or enhanced) tidal mixing at neap (or spring) tides. This indicates that taking into account the realistic spring–neap modulation of tidal mixing intensity is indispensable to further increase the ability of the existing forecast system for kyucho in the Bungo Channel.