海洋环流模式中不同近似假设下的海表高度

Boussinesq近似是现代海洋环流模式中经常采用的假设，但随着海洋模式的不断发展完善以及气候研究应用的需要，有必要估算Boussinesq近似造成的模式误差。分别利用一个非Boussinesq近似的海洋模式与另一个结构相同且采用Boussinesq近似的模式计算海表高度，并同时利用模式预报的温度、盐度资料计算了比容异常高度。分析结果显示，这3种不同定义的海表高度无论空间结构，还是时间演变，都基本类似，尤其在热带海区最接近，差值≤1cm。Boussinesq近似意味着在模式中以体积守恒代替质量守恒，通常的做法是对其进行简单的质量补偿来保持质量守恒。比较说明，以质量补偿方法进行的高度订正对减小Boussinesq近似带来的误差没有本质的意义。     

Response behaviour of triangular tension leg platforms under impact loading

Excellent station keeping characteristics and relative insensitivity with increasing water depth make triangular tension leg platforms (TLPs) a proven concept in deep water oil exploration. TLPs are often subjected to less probable forces which arise due to collision of ships, icebergs or any other huge sea creature. Dynamic analysis of two triangular TLP models at water depths 1200 and 527.8 m is performed under regular waves along with impulse load acting at an angle of 45 degrees at the TLP column. Hydrodynamic forces on these TLPs are evaluated using modified Morison equation, based on water particle kinematics arrived at using Stokes’ fifth order wave theory. Based on numerical studies conducted, it is seen that impulse loading acting on corner column of TLP significantly affect its response while that acting on pontoons dose not affect TLPs behaviour.     

A new approach for predicting the extreme response of structures

In this article a method to estimate the mean upcrossing intensity, μ(u), of a stochastic process is proposed. It is assumed that the stochastic process is a sum of a Gaussian process and a second-order correction term. The method is based on the two-dimensional Saddlepoint approximation. The accuracy of the method is tested on processes having analytical solutions for μ(u). Numerical examples are given where the stochastic process represents (i) the horizontal response of a floating offshore structure in a Gaussian sea, and (ii) the response of a structure subjected to a Gaussian wind velocity process. In addition, the estimates are compared to empirical upcrossing intensities of simulated responses. For case (ii), the obtained μ(u) estimates are compared to estimates obtained by numerical integration.     

Wave-induced porewater pressure in a seabed with inhomogeneous permeability

We have developed an analytical treatment for the wave-induced response of a seabed with inhomogeneous permeability, by employing the boundary layer approximation. By assuming an exponential function for vertically distributed permeability, the inhomogeneous consolidation equation was formulated as Bessel's differential equation. Subsequently, we proposed some representative quantities for examining the effects of inhomogeneity of permeability. The main effects of inhomogeneity are that the boundary layer becomes thinner and the phase delay in the direction of depth becomes larger as the inhomogeneity becomes significant. In addition, it was shown that the boundary layer is properly evaluated by using the averaged permeability.     

Structural monitoring of offshore platforms using impulse and relaxation response

Monitoring offshore platforms, long span bridges, high rise buildings, TV towers and other similar structures is essential for ensuring their safety in service. Continuous monitoring assumes even greater significance in the case of offshore platforms, which are highly susceptible to damage due to the corrosive environment and the continuous action of waves. Also, since a major part of the structure is under water and covered by marine growth, even a trained diver cannot easily detect damage in the structure. In the present work, vibration criterion is adopted for structural monitoring of jacket platforms. Artificial excitation of these structures is not always practicable and ambient excitation due to wind and waves may not be sufficient for collecting the required vibration data. Alternate methods can be adopted for the same purpose, for example, the application of an impact or a sudden relaxation of an applied force for exciting the structure. For jacket platforms, impact can be applied by gently pushing the structure at the fender while relaxation can be accomplished by pulling the structure and then suddenly releasing it using a tug or a supply vessel in both cases. The present study is an experimental investigation on a laboratory model of a jacket platform, for exploring the feasibility of adapting vibration responses due to impulse and relaxation, for structural monitoring. Effects of damage in six members of the platform as well as changes in deck masses were studied. A finite element model of the structure was used to analyze all the cases for comparison of the results as well as system identification. A data acquisition and analysis procedure for obtaining the response signatures of the platform due to the impulse and relaxation procedure was also developed for possible adoption in on-line monitoring of offshore platforms. From the study, it has been concluded that both impulse and relaxation responses are useful tools for monitoring offshore jacket platforms. The present work forms the basis for the development of an automated, on-line monitoring system for offshore platforms, using neural networks.     

Chebyshev逼近滤波器在位场分离中的应用

在对经典FIR数字滤波器的设计方法进行研究的基础上,提出了一种可以用于位场分离的基于Chebyshev最佳一致逼近原理的FIR滤波器的设计方法。在理论模型实验中,采用基于Hanning窗的低通滤波器计算出的区域异常最大误差为6.266×10^-6 m/s² ,均方差为2.115×10^-6 m/s² ,最大百分比误差为22.2%,而且计算点在±9 km以外的误差均大于10.1%。而利用最佳一致逼近原理分离出的区域场和局部场与理论异常值拟合得较好,曲线基本重合。分离出的区域异常最大误差为3.101×10^-6 m/s² ,均方差为0.989×10^-6 m/s² ,最大百分比误差仅在边部的几个数据上,为7.76%,其余各点的误差均小于4.1%。实例检验中将该方法用于孙吴—嘉荫剖面布格重力异常场的分离,分离出的区域场中局部场残留少,分离彻底,效果较为理想。     

Hierarchical decomposition of variance with applications in environmental mapping based on satellite images

A quadtree-based image segmentation procedure (HQ) is presented to map complex environmental conditions. It applies a hierarchical nested analysis of variance within the framework of multiresolution wavelet approximation. The procedure leads to an optimal solution for determining mapping units based on spatial variability with constraints on the arrangement and shape of the units. Linkages to geostatisiics are pointed out, but the HQ decomposition algorithm does not require any homogeneity criteria. The computer implementation can be parameterized by either the number of required mapping units or the maximum within-unit variance, or it can provide a spectrum of significances of nested ANOVA. The detailed mathematical background and methodology is illustrated by a salt-affected grassland mapping study (Hortobágy, Hungary), where heterogeneous environmental characteristics have been sampled and predicted based on remotely sensed images using these principles.     

Weak-contrast approximation of the elastic scattering matrix in anisotropic media

The plane-wave reflection and transmission coefficients at a plane interface between two anisotropic media constitute the elements of the elastic scattering matrix. For a 1-D anisotropic medium the eigenvector decomposition of the system matrix of the transformed elasto-dynamic equations is used to derive a general expression for the scattering matrix. Depending on the normalization of the eigenvectors, the expressions give scattering coefficients for amplitudes or for vertical energy flux.Computing the vertical slownesses and the corresponding polarizations, the eigenvector matrix and its inverse can be found. We give a simple formula for the inverse, regardless of the normalization of the eigenvectors. When the eigenvectors are normalized with respect to amplitudes of displacement (or velocity), the calculation of the scattering matrix for amplitudes is simplified.When the relative changes in all parameters are small, a weak-contrast approximation of the scattering matrix, based on the exactly determined polarization vectors in an average medium, is obtained. The same approximation is also derived directly from the transformed elasto-dynamic equations for a smooth vertically inhomogeneous medium, proving the consistency of the approximation.For monoclinic media, with the mirror symmetry plane parallel to the interface, the approximative scattering matrix is given in terms of analytic expressions for the non-normalized eigenvectors and vertical slownesses. For transversely isotropic media with a vertical axis of symmetry (VTI) and isotropic media, explicit solutions for the weak-contrast approximations of the scattering matrices have been obtained. The scattering matrix for amplitudes for isotropic media is well known. The scattering matrix for vertical energy flux may have applications in AVO analysis and inversion due to the reciprocity of the reflection coefficients for converted waves.Numerical examples for monoclinic and VTI media provide good agreement between the approximative and the exact reflection matrices. It is, however, expected that the approximations cannot be used when the symmetry properties of the two media are very different. This is because the approximation relies on a small relative contrast between the eigenvectors in the two media.Presented at the Workshop Meeting on Seismic Waves in Laterally Inhomogeneous Media, Castle of Trest, Czech Republic, May 22–27, 1995.