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

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

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.     

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.     

Fast computation of Hankel Transform using orthonormal exponential approximation of complex kernel function

The computation of electromagnetic (EM) fields, for 1-D layered earth model, requires evaluation of Hankel Transform (HT) of the EM kernel function. The digital filtering is the most widely used technique to evaluate HT integrals. However, it has some obvious shortcomings. We present an alternative scheme, based on an orthonormal exponential approximation of the kernel function, for evaluating HT integrals. This approximation of the kernel function was chosen because the analytical solution of HT of an exponential function is readily available in literature. This expansion reduces the integral to a simple algebraic sum. The implementation of such a scheme requires that the weights and the exponents of the exponential function be estimated. The exponents were estimated through a guided search algorithm while the weights were obtained using Marquardt matrix inversion method. The algorithm was tested on analytical HT pairs available in literature. The results are compared with those obtained using the digital filtering technique with Anderson filters. The field curves for four types (A-, K-, H-and Q-type) of 3-layer earth models are generated using the present scheme and compared with the corresponding curves obtained using the Anderson sc heme. It is concluded that the present scheme is more accurate than the Anderson scheme     

碳酸盐岩溶洞型储层反演方法

常规声波反演的方法原理和反演技术以层状介质为基础,其研究目标多是层状储层。碳酸盐岩溶洞型储层具有非规则形态、非均匀散布的特征,常规声波反演技术有其不适应之处。笔者研究的测井约束多重反演技术,解决了非层状、非均匀储层的地震反演问题,得出反映碳酸盐岩溶洞储层信息的波阻抗数据体,提取了突出溶洞型储集体低速特征的差异波阻抗,为寻找碳酸盐岩溶洞型油气藏提供了可靠的依据。