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.     

基于BP神经网络的地图数字化误差纠正

由于在数字化采集过程中不可避免地会引入系统误差和异常误差,因此消除和削弱这些误差的影响是提高空间数据质量的关键。然而由于图纸变形不均匀,扫描误差又极其复杂,用常规的多项式拟合技术只能消除部分有规律的系统误差,很难完全消除它们对地图数字化坐标的影响。BP神经网络是一个高度非线性映射系统,能以任意精度逼近。结合地图数字化坐标改正的特点,本文给出了基于BP神经网络地图数字化坐标误差纠正的方法,并通过实例验证了该方法的有效性。     

利用GPS水准点精化局部CQG2000

CQG2000模型的精度仍不能满足中、大比例尺的测图需要,如何充分利用现有成果,提高其价值和实用性,是目前需要迫切解决的问题。本文提出了通过加密GPS水准点来提高现有CQG2000精度的方法,并研究了GPS水准点数目和分布对提高结果的影响。     

一种改进的周跳多项式拟合方法

多项式拟合探测周跳受相位观测值所含误差的影响,随着采样数的增加,其探测周跳的能力随之降低,通常将观测数据分成若干段来进行处理。这样不可避免的增加工作强度,降低工作效率。本文根据多项式拟合探测周跳的基本原理,对多项式进行改进。经实验证明,利用改进后的多项式进行拟合能够有效地弥补传统方法的不足,使周跳的探测与修复更加稳定可靠。     

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%。实例检验中将该方法用于孙吴—嘉荫剖面布格重力异常场的分离,分离出的区域场中局部场残留少,分离彻底,效果较为理想。     

区域似大地水准面确定的最小二乘支持向量机方法

支持向量机(SVM)是近年来发展起来的机器学习的新方法,它较好地解决小样本、非线性、高维数、局部极小点等实际问题.文中研究支持向量机的拓展算法--最小二乘支持向量机(LSSVM),并将其应用于确定大面积复杂似大地水准面.通过工程实例并与神经网络模型和二次曲面多项式拟合模型相比较,验证确定区域似大地水准面的LSSVM方法的有效性.     

基于岭估计的有理多项式参数求解方法

在使用最小二乘法解算卫星遥感影像的RPC参数时,如果控制点非均匀分布或模型过度参数化,其法方程系数矩阵很容易产生病态,获得的解将偏离真值,甚至得到错误的解.使用岭估计可改善法方程的状态,保证解稳定.采用岭估计方法,通过所获取的不同岭参数对SPOT和QuickBird影像进行实验,证实L曲线法是一种稳定的、有效的岭参数确定方法,可显著提高RPC参数的解算精度.     

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