潮汐和流影响下长江口波浪场数值计算

采用SWAN模型和REF/DIF模型进行嵌套计算的方法来获取长江口海域实际波浪场.其中设计一种根据入射波向即时生成计算网格的方法,解决REF/DIF模型对于波浪入射角的限制从而实现两种模型的嵌套.为考虑水流和潮位分布的空间差异对波浪传播变形的作用,利用二维流场模型计算长江口的水位和流场过程,在实际波浪计算中引入了水位和水流作用.计算结果与观测资料的对比表明:1)SWAN模型和REF/DIF模型的嵌套计算方法可以作为提高浅水区域波浪计算精度的一种有效途径;2)水位和水流对长江口波浪计算的影响显著,考虑了水位和水流条件后,尤其是在大潮期间,能比较显著地提高计算精度.     

渤海海峡及附近水域水位控制的组织与实施

针对渤海海峡及附近水域开展的检查测量与航路扫海测量面临的现实困难,综合运用多种技术方法实施了水位控制:采用多站分区的最小二乘拟合法和基于潮汐模型与余水位监控法等两套水位控制方案以降低作业风险;采用虚拟站技术实现海上定点站的水位推算以提高作业效率。并对技术方案设计至水位改正数计算中的每个关键步骤都进行严密的论证分析,评估对应的量化精度指标,体现了全过程质量控制的思想。评估结果表明,水位控制的精度满足规范与行业标准的要求,达到了方案设计的目标。     

基于余水位精化潮高模型的水位生成研究

研究了基于余水位精化潮高模型的水位生成方法。所提方法首先构建网格潮汐场模型,其次采用余水位提高内插模型精度。通过网格搜索策略实现了水深点所属网格的快速定位及水位生成,探索了该模型在水位改正及预报中的应用。实验表明：所提方法既可应用于水位改正又可用于水位预报,且较现有方法在精度上有所提高。     

珠江口及邻近海域实时水位监控的可行性分析

采用基于POM模式与blending同化法构建了珠江口及邻近海域的13个主要分潮的精密潮汐模型,空间分辨率为1'×1'。精度评估表明,模型的总体综合预报误差小于9. 3cm。在广州海事测绘中心管辖的长期验潮站控制下,以余水位传递的方式,珠江口海域的实时水位计算精度能满足水位改正的相关指标要求。对于大鹏湾海域,余水位的空间一致性较强,后续若建立长期站,则也能实现水位改正精度指标下的实时水位监控。     

海南岛博鳌港洪水过程模拟

采用一个包括陆地淡水径流和风暴潮影响的模型计算了海南省博鳌港的洪水水位。模型中的外海潮汐考虑K^1和O^1分潮，风暴增水和河流流量的变化过程采用峰值函数的形式模拟，取得了与实测结果相近的计算结果。模拟实验表明，洪水水位峰值与河流洪峰流量，风暴增水峰值，系统的过水断面形状以及洪水历时的长短等因素有关。对于博鳌港口门通航条件整治，采用复合的梯形过水断面最为合适，因为相对于天然的口门过水断面，复合梯形断面能使湾内的最高洪水水位下降约20％。     

粤东船舶定线制测量中水位控制的实施

针对2014年粤东水域船舶定线制扫测项目海域范围大、潮汐类型变化明显、布站困难的现实,经分析论证和数据试验验证,采用基于潮汐模型和余水位监控法实施水位控制。验证表明,该方法所提供的测区水位数据不仅能满足相关技术标准的精度要求,而且可以提高作业效率、降低作业成本与风险。水位控制的实施过程都贯穿着"分析论证、质量控制"的作业思路,可推广应用至其他测量任务。     

珠江河口一维河网、三维河口湾水动力连接计算

对珠江河口一维河网和三维河口湾斜压模型的水动力连接计算进行了研究。通过一维河网求解公式递推出口门连接处水位、流量关系式,然后依据三维模型计算得到口门连接处流量值,根据口门连接处的水位及流量连接条件计算出口门连接处的水位值,以此水位作为边界条件分别计算一维河网和三维河口湾斜压模型,实现了一维和三维斜压模型的水动力连接计算。对模型的连接计算进行验证对比,结果基本令人满意。     

胶州湾全域潮流物理模型的相似条件及模型比尺和总体设计

胶州湾全域的变态物理模型以团岛至薛家岛的狭窄口门处为单边界.依据相似理论,确定胶州湾全域物理模型的水平比尺和垂直比尺分别为1∶700和1∶70.由几何比尺给出模型的流速比尺、糙度比尺、水流时间比尺、流量比尺及潮量比尺.设计胶州湾整体模型的试验水池为宽40 m、长50m、深1 m.试验采用翻板尾门加双向泵以水位控制为主的...     

基于ARMA模型船摇预报的船用稳定平台PID控制算法研究

传统船用稳定平台的PID控制器在闭环调节过程中存在固有的时滞,这会导致稳定平台的稳定精度无法进一步提高,针对此问题提出采用ARMA模型对船的姿态进行预测,经过处理后得到下一时刻稳定平台将要产生的误差,并将此误差作为输入施加到PID控制器中,以实现事先调节,最终达到提高稳定平台的稳定精度的目的。在ARMA模型中,采用最小二乘法对模型的参数进行实时估计,并提出使用预测值的均方误差最小的原则确定模型阶数和单次预测所使用样本数。试验结果表明,ARMA模型的预测结果具有较高的准确度,相对于传统PID控制器而言改进后的PID控制器一定程度上提高了船用稳定平台的稳定精度。     

基于潮汐模型与余水位监控法的海图测量水位改正及应用

面对近年海图测量逐渐由近海向外延伸的发展现状,其对水位改正的要求越来越高,水位改正是海图测量中的一个尤为重要的技术环节,对保证测量精度起着极其重要作用。实际测量作业中,在一些潮汐性质变化复杂的海域,采用传统水位改正方法面临设站数量多、难度大、易丢失、精度难保证等难点。针对测区潮波分布的复杂性与特殊性,基于POM模式,采用"混合"同化法同化由T/P卫星测高数据反演的沿迹潮汐参数构建了高分辨率高精度的潮汐模型,利用余水位的空间相关性,称之为基于潮汐模型与余水位监控的水位改正法。通过海图测量实例中基于潮汐模型与余水位监控法的应用,对测量成果精度进行了分析和验证,为海图测量在水位改正方法的选择提供了技术支持,进而对保证测量精度有着重要的现实意义。     

The application of the self-tuning neural network PID controller on the ship roll reduction in random waves

In this paper, we present a mathematical model including seakeeping and maneuvering characteristics to analyze the roll reduction for a ship traveling with the stabilizer fin in random waves. The self-tuning PID controller based on the neural network theory is applied to adjust optimal stabilizer fin angles to reduce the ship roll motion in waves. Two multilayer neural networks, including the system identification neural network (NN1) and the parameter self-tuning neural network (NN2), are adopted in the study. The present control technique can save the time for searching the optimal PID gains in any sea states. The simulation results show that the present developed self-tuning PID control scheme based on the neural network theory is indeed quite practical and sufficient for the ship roll reduction in the realistic sea.     

基于强化学习的锚泊辅助动力定位系统智能决策研究

针对半潜平台锚泊辅助动力定位系统的最优定位点问题,设计了基于强化学习中深度神经网络的Q学习(DQN)控制策略的锚泊辅助动力定位的智能决策系统。该决策系统中DQN方法与比例—积分—微分(PID)控制方法相结合使用,实现系统优化。在基于机器人操作系统(ROS)平台的动力定位时域模拟程序中进行数值仿真,仿真结果验证了该系统在定位点决策问题上的可靠性和有效性,从而使半潜平台在面对未知海况时,均能寻找到最优定位点,在保证锚泊辅助动力定位系统可靠性的同时降低功率消耗,提高经济性。     

Boussinesq水波方程新型数值解法

为建立高效的Boussinesq类水波数值模型,提出了一种新型的、基于有限差分和有限体积方法的混合数值格式。针对守恒形式的一维控制方程,在等间距矩形控制体内对其进行积分并离散,采用有限体积方法计算界面数值通量,剩余源项采用有限差分方法计算。其中,采用MUSTA格式并结合高精度状态插值方法计算控制体界面数值通量。时间积分则采用具有TVD性质的三阶龙格-库塔多步积分法进行。除验证模型外,重点对MUSTA格式和广泛使用的HLL格式进行了比较。结果表明,MUSTA格式可用于Boussinesq类水波方程数值求解,综合考虑数值精度、计算效率、程序编制和实际应用这几个方面,其较HLL格式更具有优势。     

An efficient shock capturing algorithm to the extended Boussinesq wave equations

A hybrid finite-volume and finite-difference method is proposed for numerically solving the two-dimensional (2D) extended Boussinesq equations. The governing equations are written in such a way that the convective flux is approximated using finite volume (FV) method while the remaining terms are discretized using finite difference (FD) method. Multi-stage (MUSTA) scheme, instead of commonly used HLL or Roe schemes, is adopted to evaluate the convective flux as it has the simplicity of centred scheme and accuracy of upwind scheme. The third order Runge–Kutta method is used for time marching. Wave breaking and wet–dry interface are also treated in the model. In addition to model validation, the emphasis is given to compare the merits and limitations of using MUSTA scheme and HLL scheme in the model. The analytical and experimental data available in the literature have been used for the assessment. Numerical tests demonstrate that the developed model has the advantages of stability preserving, shock-capturing and numerical efficiency when applied in the complex nearshore region. Compared with that using HLL scheme, the proposed model has comparable numerical accuracy, but requires slightly less computation time and is much simpler to code.     

A Hybrid Finite-Volume and Finite Difference Scheme for Depth-Integrated Non-Hydrostatic Model

A depth-integrated, non-hydrostatic model with hybrid finite difference and finite volume numerical algorithm is proposed in this paper. By utilizing a fraction step method, the governing equations are decomposed into hydrostatic and non-hydrostatic parts. The first part is solved by using the finite volume conservative discretization method, whilst the latter is considered by solving discretized Poisson-type equations with the finite difference method. The second-order accuracy, both in time and space, of the finite volume scheme is achieved by using an explicit predictor-correction step and linear construction of variable state in cells. The fluxes across the cell faces are computed in a Godunov-based manner by using MUSTA scheme. Slope and flux limiting technique is used to equip the algorithm with total variation dimensioning property for shock capturing purpose. Wave breaking is treated as a shock by switching off the non-hydrostatic pressure in the steep wave front locally. The model deals with moving wet/dry front in a simple way. Numerical experiments are conducted to verify the proposed model.     

南海表面浮力通量、水型变性及南海"暖水"的消长分析

利用英国南安普敦海洋中心(SOC)海-气界面的热量和淡水通量资料以及世界大洋图集(WOA01)的海面温度和盐度资料计算了南海表面各月的有效浮力通量及水型变性矢量场,分析了南海表层暖水形成和发展的季节性特征。研究结果表明,南海表层暖水的发展、维持以及消亡在很大程度上受到海洋表层浮力通量的影响;此外,南海上层水体密度的垂直结构和变化也深受表面浮力通量和表层水型变性的影响。在中、北部,南海的垂向结构季节性变化较为明显,其中冬季表层水体的下沉强度最大,深度最深,而夏季表层水体却无下沉趋势。     

Fòlaga: A low-cost autonomous underwater vehicle combining glider and AUV capabilities

The paper describes the current developments of a class of low-cost, light-weight autonomous underwater vehicles for coastal oceanographic applications; the vehicle class is named Fòlaga, the Italian name of an aquatic bird that swims on the water surface and dives to catch fish. The main design characteristics of the most recent vehicle of the class, the Fòlaga III, are reviewed. Navigation and control system design are discussed, with particular attention to the diving phase, which is accomplished as in oceanographic gliders by varying the vehicle buoyancy and attitude. Experimental results show that the PID robust controllers implemented are effective in the diving control phase. Finally, a distributed cooperation algorithm to be applied by a team of Fòlaga-like vehicles in adaptive oceanographic sampling applications is described. The algorithm optimizes area coverage while taking into account the accuracy in the reconstruction of the oceanographic field and inter-vehicle communication through a range constraint. The resulting dynamic programming algorithm can be implemented in a distributed fashion among the team components.     

Weak constraint data assimilation for tides in the Arctic Ocean

Estimation of the open-boundary inputs solving a weak constraint variational formulation for an Arctic tide model is considered as an ill-posed problem in the sense that the solution is very sensitive to the data noise and to grid size. Mathematically, spatial discretization of a cost function to be minimized and penalization of normal flow through the open boundary act as regularization of the problem. An heuristic choosing rule for the regularization parameter is applied to assess a suitable spatial resolution and the weight referred to the open boundary penalty. It is shown that these provide a better fit of the solution to a control data set compared with a finer grid, the value of the energy flux through the open boundary being in agreement with other model estimates. The M₂ solution obtained is much closer to the control data than other modern solutions while the accuracy of the simulated K₁ constituent is within the same error level. The tidal maps for these waves exhibit certain distinctions in comparison with other charts.     

Tidal dynamics in the strong tidal current environment of the Uldolmok waterway, southwestern tip off the Korean peninsula

Uldolmok waterway, located between an island off the southwestern tip of Korean peninsula and mainland, is famous for its strong tidal current that has a maximum current of about 6.0m/s. A series of field observations along with numerical modeling have been carried out in order to understand the tidal dynamics in terms of the force balance along the whole waterway and the energy balance in the narrowest part of the waterway. First, analysis of the ADCP current and the tide level variation data reveals that the tidal dynamics along the total waterway (channel) is balanced dominantly between the pressure gradient and linear bottom frictional forces, with the phase lag of sea level difference for the semi-diurnal constituents leading the current phase about by 10°. Secondly, the result of the numerical modeling reveals that the tidal energy flux vector flows toward the narrowest section, indicating that there should be related nonlinear processes. Through the numerical model experiment with multi-components, the convergence of (M₂ + S₂) tidal energy flux of 6.68 × 10⁷ Joule/s in the narrow area of the Uldolmok waterway is explained mainly by the energy consumption of 73% through the nonlinear generation of shallow water components and by the bottom frictional energy dissipation of 27%. This reveals that the remarkably strong nonlinear process dominates in the narrowest section of the Uldolmok waterway, compared with other areas, such as Yellow and East China Seas where the total M2 energy flux through the open boundary is balanced in terms of the bottom dissipation (Kang et al. 2003; Choi 1980).     

渤海湾水位、流场的一个数值预报模式

本文利用普林斯顿海洋模式(POM),建立了渤、黄海风增水预报模式和渤海湾水位海流预报模式;利用正交曲线网格提高重点区域的分辨率;采用美国海洋大气局(NOAA)全球预报风场和气压场作为模式表面强迫场,将计算域边界上的天文潮预报值与风增水模式预报的余水位相叠加构建模式的边界条件,在正压条件下,模拟了渤海湾2002年的水位流场过程。结果表明,模式能够较好地再现计算域内天文潮和综合水位的预报,域内10个潮位站模式与实测分析的m1和M2分潮的振幅与迟角差均不超过5.1cm和6.3°,15个潮流站模式与实测分析的m1和M2分潮流的振幅与迟角差均不超过7.5cm/s和15.8°,模式预报的水位值与塘沽站实测值非常接近,预报精度较单纯的天文潮预报有明显提高。