Evaluation of mine tailings effects on a benthic marine infaunal community over 29 years

Benthic infaunal and sediment data collected over 29 years were used to examine the extent and geographic range of effects and recovery from submarine tailings deposition. Empirical data were used to determine extreme and moderate effects, and identify near-, mid- and far-field zones. A simple probability test using overlap in frequency distributions was then used to determine less obvious effects, identify "normal" or "reference" conditions, and verify the geographic range of effect zones. Tailings thickness and sediment particulate copper were elevated from the outfall to 16-20 km distant, particularly below discharge depth (50 m). Changes in tailings thickness and copper levels before, during and after mining showed three distinct impact zones below discharge depth: near-field (<5 km from outfall); mid-field (5-16 km); and far-field (20 + km). Consistent faunal declines during mining were noted at sediment particulate copper levels >300 microg g(-1), and sediment tailings thickness > 15-20 cm. Extreme impoverishment at these sediment levels were; < 100 animals/m2, <2.5 g/m2 biomass and <7 taxa per station. No stations with >300 microg g(-1) copper and 20 cm tailings had more than 40 taxa. Amphipods were particularly affected by tailings, with distributions reflecting the instability in sediments from heavy deposition and/or shifting of tailings due to resuspension. The probability testing method assumed that reference conditions existed in far-field stations. There was no difference in abundance distributions between near-, mid- and far-field zones. Total taxa per station during mining was significantly reduced (<45) in both near- and mid-field stations compared with far-field (less than 10% overlap in distributions). One midfield station in Quatsino Sound (23) was in reference condition throughout the study period. Post-closure stations in the near- and mid-field had total taxa numbers approaching far-field conditions (overlap of 62% in distributions). Within 3 years following closure, total taxa values were within the far-field range for all stations sampled. Biomass showed a distinct decline in the near-field during mining, but with more overlap in near-field and far-field distributions (27%) than for total taxa per station, and was returning to far-field conditions in the post-closure years. A multivariate (Bray-Curtis) measure of dissimilarity indicated that the near- and mid-field infaunal compositions were distinct from the far-field during mining (< 10% overlap in distributions). The mid-field composition overlapped more with the far-field in the post-closure years, but the near-field composition did not, suggesting that the mine tailings are still affecting the fauna. The most abundant taxon in each of the reference and near-field station groups both showed significantly distinct relative abundance distributions between near- and far-field stations. In summary, the probability method showed that species richness, biomass, species composition and indicator taxa were useful for differentiating affected and non-affected stations. Polychaetes recolonize stable tailings most quickly, and have dominated the tailings stations in the post-closure years. Amphipods have recolonized sporadically, but seem to be highly sensitive to tailings stability. There was no evidence of recovery in dominant bivalves or echinoderms within three years after mine closure.     

A B-spline based higher order panel method for analysis of steady flow around marine propellers

A higher order panel method based on B-spline representation for both the geometry and the solution is developed for the analysis of steady flow around marine propellers. The self-influence functions due to the normal dipole and the source are desingularized through the quadratic transformation, and then shown to be evaluated using conventional numerical quadrature. By selecting a proper order for numerical quadrature, the accuracy of the present method can be increased to the machine limit. The far- and near-field influences are shown to be evaluated based on the same far-field approximation, but the near-field solution requires subdividing the panels into smaller subpanels continuously, which can be effectively implemented due to the B-spline representation of the geometry. A null pressure jump Kutta condition at the trailing edge is found to be effective in stabilizing the solution process and in predicting the correct solution. Numerical experiments indicate that the present method is robust and predicts the pressure distribution on the blade surface, including very close to the tip and trailing edge regions, with far fewer panels than existing low-order panel methods.     

Modeling of coastal effluent transport: an approach to linking far-field and near-field models

One of the challenges in effluent transport modeling in coastal tidal environments is the proper specification of initial dilution in connection with the far-field transport phenomena. An approach of external linkage of far-field and near-field effluent transport models is presented, and applied to simulating the effluent transport in the Port Angeles Harbor, Washington in the Strait of Juan de Fuea. A near-field plume model was used to calculate the effluent initial dilution and a three-dimensional (3-D) hydrodynamic model was developed to simulate the tidal circulation and far-field effluent transport in the Port Angeles Harbor. The hydrodynamic model was driven by tides and surface winds. Observed water surface elevation and velocity data were used to calibrate the model over a period covering the neap-spring tidal cycle. The model was also validated with observed surface drogue trajectory data. The model successfully reproduced the tidal dynamics in the study area and good agreements between model results and observed data were obtained. It is demonstrated that the linkage between the near-field and far-field models in effluent transport modeling can be achieved through iteratively adjusting the model grid sizes such that the dilution ratio and effluent concentration in the circulation model grid cell match the concentration calculated by the near-field plume model.     

NURBS高阶面元法在多体水动力分析中应用

当位于海上的多个浮体距离相近时,它们间将产生显著的水动力相互作用。使用NURBS(非均匀有理B样条)高阶面元法进行多体水动力分析,速度势用任意阶的B样条函数近似,物面用NURBS曲面表示,由于水动力与几何表示相分离,就可以方便地实现计算机辅助设计与水动力分析的集成。对截断圆柱与椭球间的水动力相互作用和半潜式平台与箱型驳船间的水动力相互作用的计算显示此种方法可以对多体间的水动力相互作用给出准确的分析,计算时间和占用内存大幅减少。     

Accurate Acoustic Signal Parameter Estimation for Marine Geodesy Surveys

Acoustic signal parameter estimation is important for diverse marine geodesy surveys and several other applications. However, the received signal from a far-field target characterized by planar wavefront propagation is frequently affected by strong nearby interfering signals. Their presence deteriorates the performance of direction-of-arrival (DOA) estimation for far-field target. In order to enhance the reception of signal from far-field target, the near-field/far-field (NFFF) beamformer is proposed. Such a beamformer optimizes beam pattern for far-field detection by maximizing beamformer output in the direction of the far-field target with the imposed condition to eliminate interfering signals generated in near-field locations. As the interference suppression only occurs at the position of near-field interference, a possible blind zone for far-field detection in conventional methods is not created. The NFFF beamformer is applicable for coherent signals and the scenario with multi interferences. For stationary situation where interferers locations are fixed, the NFFF beamformer computations do not require time updates with associated computational load. Furthermore the proposed method can be extended to several new situations such as acoustic monitoring performed from a stationary platform subjected to water currents, waves, winds and other variables, all of them generating nearby interferences and also to different array configurations including 2D and 3D arrays.     

南海北部海域油气资源调查技术及其应用研究——震源子波检测及模拟分析

在野外采集中考察气枪的稳定性以及室内进行确定性子波反褶积处理,都要求能够获得气枪阵列的实时远场子波。从震源子波检测方法出发,阐述了近场子波采集、幅值转换及记录流程,并列举出实测近场子波时需要考虑解决的问题;接着介绍了震源远场子波的获取方法和主要途径;震源近场子波及远场子波的实测方案,检测震源子波所需的仪器设备和采集技术方法。文中分析了气枪子波的气泡干扰问题,虚反射与后续气泡震荡间的关系以及吸收衰减对子波的影响,并根据实测远场子波,对齐氏和范氏两种气体模型进行了模拟效果对比,结果表明:采用范氏气体模型所模拟的子波与实际子波更为吻合,模拟出的气枪的子波波形、振幅大小、气泡周期上,与业界通用的Nucleus软件的模拟结果吻合较好,在后续的使用中将继续优化。     

Dissolved and particulate aluminum in the Columbia River and coastal waters of Oregon and Washington: Behavior in near-field and far-field plumes

The distribution of dissolved (soluble and total) and particulate (leachable and total) aluminum was examined in the Columbia River and estuary, in near-field and far-field river plumes, and in adjacent coastal waters of Washington and Oregon during the River Influence on Shelf Ecosystems (RISE) cruise of May/June 2006. Dissolved and particulate aluminum (Al) concentrations were significantly greater in the river than in the coastal waters that mixed to form the plume. Dissolved Al concentrations in the Columbia River (∼80 nM) were low relative to other major rivers. Leachable and total particulate Al concentrations within the river reached concentrations greater than 1000 nM and 18,000 nM, respectively. Dissolved Al within the Columbia River estuary showed a significant removal (∼60%) at salinities between 0 and 10 with salt-induced flocculation of colloidal Al complexes and enhanced particle scavenging being probable explanations for aluminum removal. Dissolved and particulate Al concentrations were significantly greater in near-field plumes relative to surrounding coastal waters. As the plume advected from near-field to far-field away from the river mouth, dilution of the plume with lower dissolved Al surface waters as well as particle scavenging along the flow path appeared to be controlling dissolved Al distributions. Particle settling as well as dilution with lower particle-load waters led to observed decreases in particulate Al as the plume moved from near-field to far-field. However, the percent-leachable particulate aluminum in both near-field and far-field plumes was remarkably constant at ∼7%. Dissolved and particulate Al in a far-field plume over 100 km southwest of the Columbia River mouth were over an order-of-magnitude greater than surrounding waters, illustrating the importance of the Columbia River plume as a mechanism for transporting Al offshore. Aluminum could be used to trace the input of biologically-required elements such as iron into waters off the shelf.     

Scour and Burial Mechanics of Objects in the Nearshore

A process-based, numerical, hydrodynamic vortex lattice mine scour/burial model (VORTEX) is presented that simulates scour and burial of objects of arbitrary shape resting on a granular bed in the nearshore. There are two domains in the model formulation: a far-field where burial and exposure occur due to changes in the elevation of the seabed and a near-field involving scour and transport of sediment by the vortices shed from the object. The far-field burial mechanisms are based on changes in the equilibrium bottom profiles in response to seasonal changes in wave climate and accretion/erosion waves spawned by fluxes of sediment into the littoral cell. The near-field domain consists of one grid cell extracted from the far-field that is subdivided into a rectangular lattice of panels having sufficient resolution to define the shape of the object. The vortex field induced by the object is constructed from an assemblage of horseshoe vortices excited by local pressure gradients and shear over the lattice panels. The horseshoe vortices of each lattice panel release a pair of vortex filaments into the neighboring flow. The induced velocity of these trailing vortex filaments causes scour of the neighboring seabed and induces hydrodynamic forces on the object. Scour around the object and its subsequent movement into the scour depression contribute to burial, while far-field changes in local sand level may increase burial depth or expose the object. Scour and burial predictions of mines and mine-like objects were tested in field experiments conducted in the nearshore waters off the Pacific coast of California at Scripps Pier, the Gulf Coast of Florida at Indian Rocks, and off the Atlantic coast of Massachusetts at Martha's Vineyard. Model predictions of mine scour and burial are in reasonable agreement with field measurements and underwater photographs.     

A comparison of complete and approximate solutions for second-order diffraction loads on arrays of vertical ci cular cylinders

This paper presents a comparison between two theoretical methods for computing the second-order diffraction loads on arrays of bottom-mounted, surface-piercing vertical circular cylinders in regular waves. One method presents a complete solution for the second-order hydrodynamic loads on the cylinder array via a numerical integration over the mean fluid free-surface. The other method is based on a large spacing approximation between the array members and involves the solution of a set of equivalent isolated body problems to obtain estimates for the second-order hydrodynamic loads. Numerical results for a pair of cylinders indicate very good agreement between the two methods at center-to-center spacing of both three and five radii, indicating that the approximate method may be sufficient to compute hydrodynamic interference effects to the second-order in many practical engineering situations.     

On two-dimensional moonpool resonance for twin bodies in a two-layer fluid

This paper studies the moonpool resonance of two heaving rectangular bodies in a two-layer fluid. A mathematical model is proposed based on an eigenfunction matching approach. The motion of the two-dimensional bodies is assumed to be vertical and harmonic. Heave added mass and damping coefficients are computed to examine the hydrodynamic behavior of the twin bodies. The free surface and internal wave elevations are obtained near the resonant frequencies. The presented results and analyses reveal that there exist both Helmholtz and higher-order resonances in the two-layer fluid system, which is similar to the single-layer fluid case. It is also found that the resonances are closely associated with the free surface elevation inside moonpool gap, not the wave elevation at the interfacial surface. In addition, parametric studies have been performed to identify the dependencies of hydrodynamic behavior on geometry and density stratification.     

一种污水排海的近远区计算与控制方法──以大亚湾水容量计算及污水排海为例

针对点容量计算中混合区限制的要求，采用了动量积分法、经验公式法和有限差分法，分别对污水排放后可能引起的近区、过渡区和远区的污水浓度变化进行了模拟预测，并就一定的限制性条件对容量点近区、高浓度混合区及容量点远区影响范围进行了分析。在大亚湾水容量计算及污水排海规划中的应用表明，该法具有计算机时较省、精度相对较高等优点。     

Numerical study of enhanced energy dissipation near a seamount

Using a two-dimensional primitive equation model, we examine nonlinear responses of a semidiurnal tidal flow impinging on a seamount with a background Garrett-Munk-like (GM-like) internal wavefield. It is found that horizontally elongated pancake-like structures of high vertical wavenumber near-inertial current shear are created both in the near-field (the region over the slope of the seamount) and far-field (the region over the flat bottom of the ocean). An important distinction is that the high vertical wavenumber near-inertial current shear is amplified only at mid-latitudes in the far-field (owing to a parametric subharmonic instability (PSI)), whereas it is amplified both at mid-and high-latitudes (above the latitude where PSI can occur) in the near-field. In order to clarify the generating mechanism for the strong shear in the near-field, additional numerical experiments are carried out with the GM-like background internal waves removed. The experiments show that the strong shear is also created, indicating that it is not caused by the interaction between the background GM-like internal waves and the semidiurnal internal tides. One possible explanation is proposed for the amplification of high vertical wavenumber near-inertial current shear in the near-field where tide residual flow resulting from tide-topography interaction plays an important role in transferring energy from high-mode internal tides to near-inertial internal waves.     

应用比例边界有限元法求解狭缝对双箱水动力的影响

比例边界有限元法(SBFEM)是一种半解析数值分析的新方法,既融合了有限元法和边界元法的优点,又有其特有的优点。用该方法可求解有限水深下狭缝对双箱水动力作用的影响,为波浪与多浮体超大型结构的相互作用探索一些规律。整个计算域划分成2个无限子域和4个有限子域,并利用加权余量法在各个子域上推导了SBFEM的积分方程;计算了4个数值算例并与边界元等其它数值方法进行了比较,验证了该方法是一种用很少单元便能得到精确结果的高效方法。应用SBFEM对不同箱体宽度、不同狭缝宽度、不同吃水深度条件的双箱作了计算,得出了狭缝对双箱水动力干涉影响的一些规律,对超大型浮体水动力分析和结构设计具有一定的参考价值。     

A numerical study of second-order springing of an elastic body using higher-order boundary element method (HOBEM)

Second-order springing on an elastic body with forward speed is analyzed by numerical simulations. The boundary-value problem for the velocity potential is solved by means of the direct time-domain higher-order boundary element method (HOBEM). The free-surface boundary condition in the boundary-value problem is approximated on the mean surface up to second order by use of perturbation and Taylor-series expansion methods. The body boundary condition for an elastic body is derived with various quantities which are redefined in the generalized mode. These variables such as mode shape, normal vector, etc. are obtained by using directional derivative and continuum mechanics, and the same mathematical expressions are used to obtain several second-order generalized forces. To validate the numerical results, the second-order hydrodynamic force on the bottom-mounted rigid/elastic cylinders without forward speed is compared with other semi-analytic results. The property of second-order forces on an elastic ship is studied by changing the flexural rigidity and forward speed with elastic response. It is confirmed that the second-order velocity potential is important for a body with forward speed and investigation should be made more on numerical methods for accurate computation of the second-order velocity-potential force with forward speed.     

内波远距离传播演化的实验研究

为了研究内波远距离传播过程中的演化规律,本文采用图像测速法(PIV)分别对内波近场和远场的速度场进行测量。实验中同时采用两台CCD相机对实验区域进行拍摄,根据实验结果对内波能量和垂向模态结构进行计算分析。实验结果表明,在近场区域生成的内波主要表现为内波射线结构。内波射线在经过海表面反射后,其能量在空间上出现非对称结构,能量在加强区域较减弱区增加约15%。在远场,内波射线结构不再清晰,内波主要表现为低模态内波结构。内波射线在反射时能量衰减显著,损失约为50%;低模态内波可以离开内波生成源地远距离传播,传播过程中能量损失较小,在远场传播过程中(第一模态内波半波长的距离)能量损失约20%。低模态内波的传播相速度介于垂向第一模态和第二模态相速度之间。     

Numerical investigation of the hydrodynamic interaction between two underwater bodies in relative motion

The hydrodynamic interaction between an Autonomous Underwater Vehicle (AUV) manoeuvring in close proximity to a larger underwater vehicle can cause rapid changes in the motion of the AUV. This interaction can lead to mission failure and possible vehicle collision. Being self-piloted and comparatively small, an AUV is more susceptible to these interaction effects than the larger body. In an aim to predict the manoeuvring performance of an AUV under the effects of the interaction, the Australian Maritime College (AMC) has conducted a series of computer simulations and captive model experiments. A numerical model was developed to simulate pure sway motion of an AUV at different lateral and longitudinal positions relative to a larger underwater vehicle using Computational Fluid Dynamics (CFDs). The variables investigated include the surge force, sway force and the yaw moment coefficients acting on the AUV due to interaction effects, which were in turn validated against experimental results. A simplified method is presented to obtain the hydrodynamic coefficients of an AUV when operating close to a larger underwater body by transforming the single body hydrodynamic coefficients of the AUV using the steady-state interaction forces. This method is considerably less time consuming than traditional methods. Furthermore, the inverse of this method (i.e. to obtain the steady state interaction force) is also presented to obtain the steady-state interaction force at multiple lateral separations efficiently. Both the CFD model and the simplified methods have been validated against the experimental data and are capable of providing adequate interaction predictions. Such methods are critical for accurate prediction of vehicle performance under varying conditions present in real life.     

柔性仿金枪鱼月牙形尾鳍水动力分析

为研究鱼类高速游动机理,以金枪鱼月牙形尾鳍为对象,采用面元法计算分析了其非定常水动力性能。假设尾鳍在横移和摇摆的同时,以某一匀速向前运动,并假设其在弦向和展向以某一给定规律发生变形,以模仿柔性变形。探讨了前进速度、横移和摇摆的幅度、频率及其相位差对柔性尾鳍推进性能的影响,并与做相应运动的刚性尾鳍进行了对比。     

波浪与结构物作用分析的一种高阶边界元方法——自由项和柯西主值积分的直接数值计算

采用直接数值计算方法计算了势流问题高阶边界元方法中的自由项系数和柯西主值积分,建立了波浪与结构物作用的一种高阶边界元方法.通过算例研究了物体表面上固角系数的计算精度和不同网格剖分、不同阶高斯积分点对柯西主值积分的影响.对截断圆柱上的波浪作用力与解析解做了对比,发现本方法具有很高的计算精度,随网格的加密迅速收敛于解析解.     

Extension of Haskind's relations to cylindrical wave fields in the context of an interaction theory

The Direct Matrix Method Interaction Theory (IT) proposed by Kagemoto and Yue [1] speeds up the computation of hydrodynamic coefficients for large arrays of bodies when compared to direct calculations using standard Boundary Element Method (BEM) solvers. One of the most computationally expensive parts of the matrix method is the calculation of two hydrodynamic operators, known as Diffraction Transfer Matrix (DTM) and Radiation Characteristics (RC), which describe the way an isolated geometry scatters and radiates waves, respectively. A third operator, called Force Transfer Matrix (FTM), was introduced by McNatt et al. [2] to facilitate the calculation of the forces exerted on the bodies. In this paper, a novel set of relations between the FTM and RC components is obtained using the Kochin functions specific to the cylindrical basis solutions. They extend the classical Haskind's relations, valid with incident plane waves, to the cylindrical components of the scattered and radiated fields. Moreover, an alternative demonstration of the identities is given, which does not rely on the far-field asymptotic representation of the potential. Additional expressions are provided that relate the hydrodynamic coefficients and the RC for isolated bodies as well as for arrays, and numerical checking of the derived mathematical expressions is presented. These new relations can be used to speed up calculation of the hydrodynamic operators required for the use of the IT and to test its accuracy.     

A Numerical Model for Mixing Areas of Pollutant Discharging into Tidal Flows

- A combined numerical model for computing mixing areas of pollution vertical jet discharging into tidal flows has been developed. This numerical model is composed of a 2-D depth-averaged dynamic far-field numerical model and a 2-D vertical dynamic near-field numerical model. The former uses finite node method to compute velocity, and improved finite node method to compute pollutant concentration distribution; the latter is a k-e turbulence model, and uses SIMPLE (Semi-Implicit Method for Pressure Linked Equations) method to compute velocity. At the same time, the former provides boundary conditions for the latter. This model can simulate both far- field pollutant concentration distribution and near-field vertical recirculation quickly and precisely. This model has been verified by measured data of pollutant F of the Dachang reach of the Changjing River and test data presented by Chen el al. (1992). On the basis of verification, the authors use a designed hydrograph to compute this mixing area for a cer