埋地目标体矢量电磁散射的一种快速正演算法

利用积分方程方法以及半空间并矢格林函数的快速算法对埋地目标体矢量电磁散射进行正演计算。首先,利用半空间电磁并矢格林函数建立起埋地目标体的体积分方程。然后通过将空间偏导转移至格林函数谱域积分的积分号之外,并采用离散复镜像方法来近似余下的零阶索末菲积分,进一步得到并矢格林函数各个分量的闭合形式。由于避免了对索末菲积分的繁琐数值计算,使得生成反应矩阵和计算散射场时由半空间并矢格林函数计算带来的瓶颈问题得到较好的克服,因而极大地提高了埋地目标体电磁响应正演计算效率,同时也能保证足够的精度。     

电磁波法探测地下异常体的数值模型和实验模型

本文研究均匀有耗媒质全空间及半空间中三维异常体散射的数值模型和实验模型。用矩量法解体积分方程,导出自小块积分的简单准确的计算公式,说明非自小块积分以及索末菲积分的简化计算方法。给出两个利用金属镜象板的水池测试装置,说明与测试方法有关的若干主要问题。计算结果与测量结果基本吻合。     

二维均匀介质中积分方程的快速精确算法

采用稳定型双共轭梯度快速Fourier变换(BCGS-FFT)算法精确计算二维均匀介质中的积分方程.采用一种新的插值函数作为基函数和试探函数对积分方程进行弱化离散,离散后的积分方程采用稳定型双共轭梯度迭代方法进行求解,从而得到异常体内电场的分布.计算时采用快速Fourier变换技术将积分方程内Green函数与电场的乘积表示成褶积形式以加快计算速度.数值计算举例说明了算法的精确性和有效性.     

单程波算子积分解的象征表示

单程波波场延拓算子在地震偏移成像中有重要应用.单程波波场延拓算子按其实现方式可分为Kirchhoff积分、空间隐式有限差分和Fourier变换方法，他们代表了算子的不同表示方法，当截断使用这些方法时会得到不同的精度.象征表示对这些方法的导出和精度分析有重要作用.算子作用于正弦波函数所得函数称为算子的象征.算子的象征是褶积算子Fourier变换的推广.Fourier变换方法则直接用象征函数的可分表示求出.空间隐式有限差分则可以用象征函数的Padè近似或部分分式导出.单程波算子在深度域的积分称为单程波算子积分解.本文推导了单程波算子积分解的象征表达式，给出了算子象征的代数运算的头几阶表达式，这些表达式还未在前人文献中发现.Kirchhoff积分所需格林函数可以通过象征函数和鞍点法导出.基于积分解的象征表达式给出了非对称走时公式，对改善Kirchhoff积分的聚焦性能有重要意义.     

格林函数的奇异性处理

格林函数在解决电磁场和直流电场问题的积分方程法和边界单元法等方法中有着广泛的应用.在应用格林函数时会遇到格林函数奇异性的问题,即其在源区域的内部,当场点和源点重合时,会碰到对格林函数为无穷大,积分为奇异积分.对于不同问题中遇到的不同形式的格林函数奇异性问题的各种处理方法进行了分析和评述.对于标量格林函数奇异性问题处理方法有：挖去法、级数展开法、绕开法和解析法等.对于并矢格林函数奇异性问题处理方法有：分量处理法、源并矢法和拟源并矢法等处理方法.通过对实际工作中所遇到的新方法新问题中的不同形式的格林函数的研究,提供了奇异性问题的处理的方法和途径.     

任意点震源τ-P域地震图的解及AVO分析

用严格的数学推导得出τ-P域地震图的解析解。分析表明,τ-P域中地震图的子波可表示为激发震源子波的积分形式,决定反射波振幅的因素是平面波反射系数、反射界面上覆介质的波速以及方向因子。这些结果说明:1.在τ-P域进行AVO分析是可行的,并具有更高的精度;2.从一种新的途径可以进行速度函数的反演计算。     

一种随机粗糙海面电磁散射的极化核函数

从磁场的积分方程出发,应用迭代方法推导了二次迭代条件下粗糙海面的近似电磁散射极化核函数,并对其进行了分析.这种函数的表示是在光学近似即KA近似的基础上的一种扩展形式,因而可称为扩展KA近似(EKA).Elfouhaily采用了与Holliday类似的方法分析了双站配置条件下散射的极化核函数,试图保留所有的高次项,但是由于处理的复杂性他只给出了斜率项在散射项中的表示形式,并称之为KA近似的"补充项".本文在此基础上进一步考虑了Elfouhaily近似的形式,给出了表示双站散射极化核函数的完整表达形式,并以此为基础推算了后向散射条件下极化核函数的理论结果,并与Elfouhaily的结果进行了对比,说明了极化核函数完整表达形式的合理性.近似散射极化核函数的分析对于深入理解粗糙表面(特别是海面)电磁散射特性以及散射近似模型的建立等问题的研究具有重要意义.     

沿球面大圆的积分逆变换及其应用

本文给出了从沿球面上大圆的积分值直接反推球面上函数分布的关系式,即得到了球面上大圆的积分变换和反变换对。如利用的是整个大圆的积分,得到的是球面上任何点及其对蹠点上函数的平均值,或者说得到的是球面上函数的偶对称部分。如果进行积分的球面域小于半球面,反演将是唯一的。所得结果可应用于从地震表面波资料反演三维地球内部构造的地球物理问题。     

利用高斯求积和连分式展开计算电磁张量格林函数积分

针对于利用快速汉克尔变换求解电磁张量格林函数中含有贝塞尔函数积分项时，积分核函数不总是随贝塞尔函数（Jv(λr)）参数λ的增大而趋于收敛，从而产生错误结果的问题，我们利用一种直接数值积分——结合连分式展开的高斯求积方法进行计算，通过对不同模型的试算表明这种方法总能够保证电磁张量格林函数积分的正确计算.     

地震道重建和重采样的分形插值方法研究

为了对地震数据进行高保真重建,本文给出了分形插值函数的一种替代形式,即利用某一函数的逐段积分表达形式构建相应的分形插值函数,这种形式在以往的地震数据处理文献中未曾提及.这类分形插值函数在分形和线性插值之间建立了一种简单的关系.实际上,利用这种简单关系,许多分形插值函数问题可以简化为一种简单的向量/矩阵表示.此外,利用这种替代形式,借助于压缩因子的显式表达,我们可以进行高精度的地震数据重建.本文借助于这种分形插值函数的替代形式,发展了一种新的地震数据重建和重采样的方法.为了检验本文方法的有效性,我们采用济阳凹陷10个在一条线上的地震台站的数据对该方法进行了检验,实际地震图重建结果表明了本文提出方法的有效性.     

各向异性半空间上方垂直磁偶极子电磁场的球面波展开式

文中推导出单轴各向异性媒质半空间上方垂直磁偶极源电磁场的精确解析表达式. 应用圆柱波函数的球面波展开式和超几何函数理论，场分量中的索末菲尔德型积分被表示成快速、绝对收敛的球面波函数系展开式；展开系数是以物性参数为复宗量的勒让德多项式. 该展开式数学物理意义明显，并且不受场点和源点的位置、媒质的物性参数和频率等条件的限制. 利用本文的结果可十分方便地计算和分析任意场点处的电磁场分布.     

Analytical Solutions of Three‐Dimensional Contaminant Transport Models with Exponential Source Decay

Despite the availability of numerical models, interest in analytical solutions of multidimensional advection‐dispersion systems remains high. Such models are commonly used for performing Tier I risk analysis and are embedded in many regulatory frameworks dealing with groundwater contamination. In this work, we develop a closed‐form solution of the three‐dimensional advection‐dispersion equation with exponential source decay, first‐order reaction, and retardation, and present an approach based on some ease of use diagrams to compare it with the integral open form solution and with earlier versions of the closed‐form solution. The comparison approach focuses on the relative differences associated with source decay and the effect of simulation time. The analysis of concentration contours, longitudinal sections, and transverse sections confirms that the closed‐form solutions studied can be used with acceptable approximation in the central area of a plume bound transversely within the source width, both behind and beyond the advective front and for concentration values up to two orders of magnitude less than the initial source concentration. As the proposed closed‐form model can be evaluated without nested numerical computations and with simple mathematical functions, it can be very useful in risk assessment procedures.     

Integral inversion of SST data of type GRACE

Satellite missions CHAMP and GRACE dedicated to global mapping of the Earth’s gravity field yield accurate satellite-to-satellite tracking (SST) data used for recovery of global geopotential models usually in a form of a finite set of Stokes’s coefficients. The US-German Gravity Recovery And Climate Experiment (GRACE) yields SST data in both the high-low and low-low mode. Observed satellite positions and changes in the intersatellite range can be inverted through the Newtonian equation of motion into values of the unknown geopotential. The geopotential is usually approximated in observation equations by a truncated harmonic series with unknown coefficients. An alternative approach based on integral inversion of the SST data of type GRACE into discrete values of the geopotential at a geocentric sphere is discussed in this article. In this approach, observation equations have a form of Green’s surface integrals with scalar-valued integral kernels. Despite their higher complexity, the kernel functions exhibit features typical for other integral kernels used in geodesy for inversion of gravity field data. The two approaches are discussed and compared based on their relative advantages and intended applications. The combination of heterogeneous gravity data through integral equations is also outlined in the article. panovak@kma.zcu.cz     

The disturbance due to a point source in a homogeneous liquid layer over a heterogeneous liquid half-space

Summary The problem of a periodic point source in a homogeneous liquid layer overlying a heterogeneous liquid half-space is discussed. After obtaining the formal solution, the path of integration for the displacement potential of the layer is transformed from the positive real axis to the positive imaginary axis and the Sommerfeld contour and the latter is further distorted to the modified Sommerfeld contour. The residues of the integrand at the poles contained within the Sommerfeld loop constitute the normal mode solution to the problem. The integrands in the expressions for the integrals along the imaginary axis are expanded in a series of negative powers of exponentials and then some of the terms in these expansions are evaluated approximately. This gives various waves reflected from the interface and the integral along the Sommerfeld loop vanishes. The frequency equation is obtained, also by the principle of constructive interference. An expression for the reflection coefficient at an interface of two liquid media, the upper medium being homogeneous and the lower one inhomogeneous, is obtained.     

利用积分方程法计算双侧向测井仪的井眼校正曲线

对电测井积分方程中的各积分项进行了理论推导，利用任意闭合曲面对曲面外一点所张立体角为零的特点将电势积分项由对柱体侧面的积分转化为对柱体底面的积分，从而使在纵向采用线性插值条件下给出了电势积分项的解析解．对电流积分项进行了化简，使之由曲面积分化为一维数值积分．利用所得结果计算了双侧向测井仪的井眼校正曲线．     

ON THE CONSTRUCTION OF THE 3D BAND-LIMITED EXTRAPOLATION OPERATOR IN THE SPACE-FREQUENCY DOMAIN1

Based on an expansion of the band-limited 3D extrapolation operator in terms of orthogonal Chebyshev polynomials, a closed form expression of the space-frequency response is presented. A key step is an evaluation of the (inverse) 2D Fourier transform of circularly symmetric functions, which is related to the (zero-order) Hankel transform. Hankel transforms of individual members of the orthogonal set of polynomials are available from tables and summation of series; hence, the real and the imaginary parts of the space-frequency response can be found in terms of cylindrical and spherical Bessel functions, respectively. The procedure permits an efficient and accurate evaluation of the space-frequency response.     

Exact Elastodynamic Green Functions for Simple Types of Anisotropy Derived from Higher-Order Ray Theory

Using higher-order ray theory, we derived exact elastodynamic Green functions for three simple types of homogeneous anisotropy. The first type displays an orthorhombic symmetry, the other two types display transverse isotropy. In all cases, the slowness surfaces of waves are either ellipsoids, spheroids or spheres. All three Green functions are expressed by a ray series with a finite number of terms. The Green functions can be written in explicit and elementary form similar to the Stokes solution for isotropy. In two Green functions, the higher-order ray approximations form a near-singularity term, which is significant near a kiss singularity. In the third Green function, the higher-order ray approximations also form a near-field term, which is significant near the point source. No effect connected with the line singularity was observed.