任意各向异性介质三维有限元航空电磁响应模拟

在航空电磁探测方法应用中,地下探测环境复杂,地下介质物性参数存在各向异性特征。如果采用常规各向同性介质模型,在资料解释过程中会产生严重的偏差。本文基于矢量有限元法开展三维任意各向异性频率域航空电磁响应模拟计算研究。首先,将总场分解成一次场和二次场,对空气介质的均匀全空间进行一次场解析计算,同时利用矢量有限元法对二次电场的双旋度方程进行求解。为提高求解速度,采用共享内存直接求解器PARDISO对大规模稀疏矩阵并行计算,大大提高了三维模型的计算速度。之后开展了4种典型目标模型的三维各向异性介质模拟：围岩各向同性-目标体各向异性（绕z轴旋转）模型;围岩各向同性-目标体各向异性（绕x轴旋转）模型;围岩各向异性-目标体各向异性（绕z轴旋转）模型;围岩各向异性-目标体各向异性（绕x轴旋转）模型。对比分析了不同模型在不同旋转角度情况下,磁场实虚分量的变化特征,进而总结了各向异性参数对航空电磁响应的影响规律和识别方法。     

Numerical Calculation for Nonlinear Waves in Water of Arbitrarily Varying Depth with Boussinesq Equations

Based on the high order nonlinear and dispersive wave equation with a dissipalive term, a numerical model for nonlinear waves is developed. It is suitable to calculate wave propagation in water areas with an arbitrarily varying bottom slope and a relative depth h/ L0≤ 1. By the application of the completely implicit slagger grid and central difference algorithm, discrete governing equations are obtained. Although the central difference algorithm of second-order accuracy both in time and space domains is used to yield the difference equations, the order of truncation error in the difference equation is the same as that of the third-order derivatives of the Boussinesq equation. In this paper, the correction to the first-order derivative is made, and the accuracy of the difference equation is improved. The verifications of accuracy show that the results of the numerical model are in good agreement with those of analytical solutions and physical models.     

SURFACE DISPLACEMENTS OF A NON-HOMOGENEOUS ELASTIC HALF-SPACE SUBJECTED TO UNIFORM SURFACE TRACTIONS. PART I: LOADING ON ARBITRARILY SHAPED AREAS

A numerical technique is presented for the analysis of surface displacements of a non-homogeneous elastic half-space subjected to vertical and/or horizontal surface loads uniformly distributed over an arbitrarily shaped area. The non-homogeneity considered is a particular form of power variation of Young's modulus with depth. Since the exponent which determines the degree of non-homogeneity may vary from zero to unity, both the homogeneous half-space and the Gibson soil may be included as limiting cases in a single numerical scheme. In order to account for the arbitrary shape of the loading, the boundary of the loaded area is linearized piecemeal. This enables the modeling of any load pattern according to the desired degree of accuracy. Special attention is focused on the integration scheme, since the singularity associated with the Green's function becomes progressively more pronounced the greater the non-homogeneity parameter gets. The performance of the numerical procedure is studied using analytical solutions for rectangular shaped areas. Further comparisons with well-known solutions based on integral transform techniques for a uniformly distributed load acting on a circular area of the non-homogeneous soil mass show excellent agreement as well. © 1997 John Wiley & Sons, Ltd.     

双排弹性桩隔振屏障对平面SH波的多重散射

提出了一种新的求解任意排列、任意半径的弹性桩对平面SH波的多重散射的理论方法,以解决以往假设单重散射的计算方法中不考虑桩列作为整体屏障从而忽略桩间相互干涉关系的不足,并且可用于分析多排弹性桩对平面SH波的散射性状。在数值计算分析中讨论了散射重数,排间距,桩土剪切模量比,桩数等对双排弹性桩屏障隔离效果的影响,从而对实际工程中利用排桩进行振动污染的治理和屏障隔振设计提出了重要的理论依据。     

等效源法三维随机点位场数据处理和转换

为了实现曲面随机点位场数据的曲面延拓和转换,以磁异常位场数据为例,采用一组磁偶极子作为等效源,置于观测面下方的一个曲面上,把观测磁异常作为这组磁偶极子所产生磁异常的边界条件,通过求解线性方程组的方法反演磁偶极子磁矩的大小,再根据反演结果正演所要计算的磁异常数据,实现了曲面随机点磁异常位场数据的向上延拓、向下延拓、求导以及化极处理.在数据量较大时,为了提高反演计算的速度,把磁异常数据和磁偶极子分成若干小块,再利用各块磁异常数据分别反演该块数据下方磁偶极子的磁矩,并通过迭代计算来逐步取得更准确的反演结果.模型试验表明,磁异常位场数据向上延拓的均方根误差小于±2nT,向下延拓和化极也可以取得较高的精度,所提出的分块处理方法提高了延拓和转换的速度,实际资料处理给出了曲面随机点航磁异常数据向下延拓和化极的一个例子.     

Prestack reverse-time depth migration of arbitrarily wide-angle wave equations

Based on arbitrarily wide-angle wave equations,a reverse-time propagation scheme is developed by substituting the partial derivatives of depth and time with central differences. The partial derivative of horizontal direction is replaced with high order difference. The imaging condition is computed by solving the eikonal equations. On the basis of above techniques,a prestack reverse-time depth migration algorithm is developed. The processing exam-ples of synthetic data show that the method can remove unwanted internal reflections and decrease the migration noise. The method also has the advantage of fidelity and is applicable of dip angle reflector imaging.