Improved assessments of wave-induced fatigue for free spanning pipelines

For subsea pipeline projects, the costs related to seabed correction and free span intervention are often considerable. Development of reliable methods for fatigue analyses of pipelines in free spans contributes to minimize costs without compromising pipeline integrity. Assessment of wave-induced fatigue damage on multi-span pipelines is investigated, and improved analysis methods are suggested in this paper. A time-domain (TD) algorithm is developed, which accounts for non-linear hydrodynamic loading and dynamic interaction between adjacent spans. The proposed TD approach is employed to evaluate linearized frequency-domain (FD) solutions from recognized design standards and to study the dynamic response of multi-span pipelines to direct wave loading. Differences between multi- and single-span analyses are described for the first time, and the common assumption that the main fatigue damage contribution comes from the fundamental mode is demonstrated not to hold for multi-spans. An improved FD solution capable of predicting multi-mode response is derived and demonstrated to give accurate fatigue life estimates for multi-span pipelines.     

Estimation of mean and extreme waves in the East China Seas

Accurately estimating the mean and extreme wave statistics and better understanding their directional and seasonal variations are of great importance in the planning and designing of ocean and coastal engineering works. Due to the lack of long-term wave measurement data, the analysis of extreme waves is often based on the numerical wave hind-casting results. In this study, the wave climate in the East China Seas (including the Bohai Sea, the Yellow Sea and the East China Sea) for the past 35 years (1979–2013) is hind-casted using a third generation wave model – WAMC4 (Cycle 4 version of WAM model). Two sets of reanalysis wind data from NCEP (National Centers for Environmental Prediction, USA) and ECMWF (European Centre for Medium-range Weather Forecasts) are used to drive the wave model to generate the long-term wave climate. The hind-casted waves are then analysed to study the mean and extreme wave statistics in the study area. The results show that the mean wave heights decrease from south to north and from sea to land in general. The extreme wave heights with return periods of 50 and 100 years in the summer and autumn seasons are significantly higher than those in the other two seasons, mainly due to the effect of typhoon events. The mean wave heights in the winter season have the highest values, mainly due to the effect of winter monsoon winds. The comparison of extreme wave statistics from both wind fields with the field measurements at several nearshore wave observation stations shows that the extreme waves generated by the ECMWF winds are better than those generated by the NCEP winds. The comparison also shows the extreme waves in deep waters are better reproduced than those in shallow waters, which is partly attributed to the limitations of the wave model used. The results presented in this paper provide useful insight into the wave climate in the area of the East China Seas, as well as the effect of wind data resolution on the simulation of long-term waves.     

Interaction of cnoidal waves with an array of vertical concentric porous cylinders

An array of large concentric porous cylinder arrays is mounted in shallow water exposed to cnoidal waves. The interactions between waves and cylinders are studied theoretically using an eigenfunction expansion approach. Semi-analytical solutions of hydrodynamic loads and wave run-up on each cylinder are obtained using first approximation to cnoidal waves. The square array configuration of four-legged identical concentric porous cylinder is investigated in present study. Numerical results reveal the variation of dimensionless wave force and wave run-up on individual cylinder with angle of incidence, porosity parameter, spacing between outer and inner cylinders, spacing between concentric porous cylinders and wave parameter. Different mechanism of wave force is found under different range of scattering parameter.     

Experimental and theoretical study of a cylindrical oscillating water column device with a quadratic power take-off model

The wave power extraction by a cylindrical oscillating water column (OWC) device with a quadratic power take-off (PTO) model was studied experimentally and theoretically. In the experiment, a scaled model OWC was tested in a wave flume, with an orifice being used to simulate a quadratic PTO mechanism. In the theoretical analysis, the quadratic PTO model was linearized based on Lorenz's principle of equivalent work, which allows us to perform a frequency domain analysis using an eigen-function matching method. The effects of higher harmonic components and the spatial non-uniformity of the surface velocity inside the chamber were discussed. A semi-analytical model was proposed to understand the viscous loss affecting the measured capture length. Our treatment of the quadratic PTO model was validated by comparing quasi-linear theoretical capture length and the laboratory measurement. Our results also showed that the effects of spatial non-uniformity and viscous loss could be noticeable for shorter waves.     

Stability and capsizing analysis of nonlinear ship rolling in wind and stochastic beam seas

Considering the actual seaway condition, stability and capsizing of nonlinear ship rolling system in stochastic beam seas is of significant importance for voyage safety. Safe zone are defined in the phase space plan of the unperturbed Hamilton system to qualitatively distinguish ship motions as capsize and noncapsize. Capsize events are defined by solutions passing out of the safe zone. The probability of such an occurrence is studied by virtue of the random Melnikov function and the concept of phase space flux. In this paper, besides conventional wave excitation, the effect of wind load is also taken into account. The introduction of wind load will lead to asymmetry, in other words, it transforms the symmetric heteroclinic orbits into asymmetric homoclinic orbits. For asymmetric dynamical system, the orbital analytic solutions and its power spectrum are not readily available, and the technique of discrete time Fourier transformation (DTFT) is used. In the end, as verification of theoretical critical significant wave height, capsizing probability contour diagram is generated by means of numerical simulation. The contour diagram shows that these analytical methods provide reliable and predictive results about the likelihood of a vessel capsizing in a given seaway condition.     

CFD investigations of OXYFLUX device,an innovative wave pump technology for artificial downwelling of surface water

No other environmental variable of such ecological importance to estuarine and coastal marine ecosystems around the world has changed so drastically, in such a short period of time, as dissolved oxygen in coastal waters. The prevalent methods for counteracting anoxic sea events are indirect measures which aim to cut-down anthropic loads introduced in river and marine environments. To date, no direct approaches, like artificial devices have been investigated except the WEBAP and OXYFLUX devices. The present paper adopts a numerical approach to the analysis of the pumped surface water as well as the analysis of the dynamic response of the OXYFLUX device. By means of a CFD-RANS code and through the application of overset grid method, the 1/16 OXYFLUX model’s dynamic response and pumping performance are evaluated. The appropriate grid is selected after an extensive sensitivity analysis carried out on 9 different grids. The CFD model is validated by comparing numerical and physical results of heave decay test, heave response, and surface water discharge under the action of regular waves. The extensive comparison with experimental results shows consistently accurate predictions. The main findings of the study show that nonlinear effects remarkable reduce the dynamic behaviour of the OXYFLUX and generate an unexpected second harmonic for pitch response intensifying the overtopping discharge also for small waves caused by the summer's low intensity winds.     

基于OS071X测波雷达的波浪探测分析

文章利用舟山海洋环境监测站2014年第16号台风"凤凰"过境影响期间,对OS071X测波雷达和3m波浪浮标在同一海域内的波浪探测资料进行了对比验证。通过分析,OS071X测波雷达探测反演的有效波高、周期与3m波浪浮标探测结果对比,二者按先后顺序的均方根误差分别为0.27m和0.29m,相关系数分别为0.93和0.90。结果表明,OS071X测波雷达系统的探测精度得到了验证,能够满足实时探测近岸海域波浪的需求。     

锂光卤石X-ray衍射标准图谱

盐湖卤水提锂过程中需了解与掌握锂元素结晶析出与赋存状态,锂光卤石(LiCl·MgCl_2·7H_2O)作为重要的中间产物,其结构与种类均需鉴别,但由于缺乏锂光卤石的X-ray衍射标准图谱,使得常见的Xray衍射方法无法应用。利用相图原理,通过湿固相法制备高纯锂光卤石,然后对其进行X-ray衍射分析,从而提供锂光卤石X-ray衍射的标准图谱,为高镁锂比盐湖卤水提锂过程分析提供可靠的理论支撑。     

FLOOD MODELLING WITH GIS-DERIVED DISTRIBUTED UNIT HYDROGRAPHS

The concept of a spatially distributed unit hydrograph is based on the fact that the unit hydrograph can be derived from the time–area curve of a watershed by the S-curve method. The time–area diagram is a graph of cumulative drainage area contributing to discharge at the watershed outlet within a specified time of travel. Accurate determination of the time–area diagram is made possible by using a GIS. The GIS is used to describe the connectivity of the links in the watershed flow network and to calculate distances and travel times to the watershed outlet for various points within the watershed. Overland flow travel times are calculated by the kinematic wave equation for time to equilibrium; channel flow times are based on the Manning and continuity equations. To account for channel storage, travel times for channel reaches are increased by a percentage depending on the channel reach length and geometry. With GIS capability for rainfall mapping, the assumption of a uniform spatial rainfall distribution is no longer necessary; hence the term, spatially distributed unit hydrograph. An example of the application for the Waiparous Creek in the Alberta Foothills is given. IDRISI is used to develop a simple digital elevation model of the 229 km² watershed, using 1 km × 1 km grid cells. A grid of flow directions is developed and used to create an equivalent channel network. Excess rainfall for each 1 km × 1 km cell is individually computed by the Soil Conservation Service (SCS) runoff curve method and routed through the equivalent channel network to obtain the time–area curve. The derived unit hydrograph gave excellent results in simulating an observed flood hydrograph. The distributed unit hydrograph is no longer a lumped model, since it accounts for internal distribution of rainfall and runoff. It is derived for a watershed without the need for observed rainfall and discharge data, because it is essentially a geomorphoclimatic approach. As such, it allows the derivation of watershed responses (hydrographs) to inputs of various magnitudes, thus eliminating the assumption of proportionality of input and output if needed. The superposition of outputs is retained in simulating flood hydrographs by convolution, since it has been shown that some non-linear systems satisfy the principle of superposition. The distributed unit hydrograph appears to be a very promising rainfall runoff model based on GIS technology.     

柴达木低品位固体钾矿溶解转化率与品位关系

我国是缺钾大国,探明的可溶性钾盐十分有限,柴达木盆地是我国盐湖钾盐的聚集地,含水层骨架中含有数亿吨低品位固体钾矿,开发出这些低品位固体钾盐对增加国内钾肥供给,保障粮食安全具有重大意义。青海盐湖工业股份有限公司历经十年开发出低品位固体钾盐的浸泡式溶解转化方法,试验表明,固体钾盐的溶解转化率随固体钾品位的增加而增加,累计溶解转化率为56%～98%。该技术在柴达木多个盐湖得到推广应用,察尔汗盐湖保有固体钾盐达2.96×10~8 t,平均品位KCl为1.24%,根据试验数据回归分析溶解转化率77%,可溶解转化出钾盐(KCl)约2.28×10~8 t,柴达木盆地盐湖保有3.43×10~8 t固体钾盐,平均品位KCl为1.25%,溶解转化率78%,可溶解转化出钾盐(KCl)2.68×10~8 t,大幅度增加了钾盐储备量。