三维不规则地形河谷场地地震响应分析方法研究

本文提出了不规则地形地震响应分析的自由场计算方法,并以小湾拱坝场址作为三维场地模型,比较了三维不规则场地侧边界自由场的精确解和近似解;同时,采用大型有限元软件Ansys/Ls-Dyna计算了自由场对场地内部场点地震响应的影响.计算结果表明,三维场地模型的二维自由场近似解与精确解有很大差异,其结果对三维场地内部场点地震响应影响显著,因此,要准确求解三维不规则地形的地震响应必须采用真实的二维侧边界模型.     

瞬变电磁场时域格林函数解

近源时间域电磁场具有信号强、探测深度大、精度高等优点,但传统勘探电磁场理论中偶极子近似在近源会引起较大误差,导致这一优势的发挥受到了制约.开展直接时间域电磁场解析式研究,是解决这一问题的途径之一.本文提出在点电荷微元假设下,引入时域格林函数,求取瞬变电磁场时间域解析解.采用积分运算法,把电磁场阻尼波动方程的求解问题转化为求其格林函数积分形式解的问题;建立辅助路径解决奇点问题,利用复分析中的约当引理、留数定理和广义函数等理论和方法,推导计算出时间域格林函数的时空四重广义积分.得到达朗贝尔方程的直接时域格林函数精确解析式,与传统方法“比拟”出的公式具有相同的形式,验证了本文推导的时域格林函数解析公式的正确性;推导出扩散方程的直接时间域解析解.通过与时变点电荷源时间域的电磁响应近似表达式进行对比,得出本文所推导的公式计算精度较高的结论;建立了全空间回线源瞬变电磁场问题的直接时间域求解公式.为解决全场区瞬变电磁场精细探测直接时域解析问题提供了基础理论.     

均匀半空间瞬变电磁场直接时域响应数值分析

近源时域电磁场具有信号强、探测深度大和精度高等优点,但传统瞬变电磁场理论中偶极子近似在近源区会引起较大误差,推导瞬变电磁场直接时域解析式是解决这一问题的关键.本文在点电荷微元假设下通过时域格林函数,采用分离变量等方法推导出了上半空间一次有源波动场和反射波的时域解析式和下半空间二次无源波动场的时域解析式,结合均匀半空间瞬变电磁场的边界条件给出了均匀半空间瞬变电磁场的直接时域解析式,进而利用第一型曲线积分,通过沿回线源叠加推导出圆回线源在瞬变电磁场中的直接时域解析式.然后在半空间表面上,与传统的电偶极源假设下的表达式作了比较.数值结果表明两者在远源区的计算结果相差甚微,而近源区则存在很大误差.本文利用真正点元（点电荷）严密推导给出的均匀半空间表面上瞬变电磁场的直接时域解析式适用于全场区探测,克服了偶极子假设下只适用远场区的不足,为瞬变电磁法的进一步发展和实际勘探提供了新的理论基础.     

地磁场的电流环模式

地磁场的模式除了球谐模式和偶极子模式以外,还有电流环模式。地磁场的球谐模式使用广泛,但它不能反应地磁场场源的实际状况;地磁场的偶极子模式只能很近似地反应地磁场场源的状况;为此,有些学者提出用若干个电流环表示地磁场的分布,偶极子是电流环的特殊情况,当电流环的半径为零时,电流环就变为偶极子,电流环模式比偶极子模式能更好地表示地磁场场源的实际状况。     

点电荷微元与偶极子源的时域电磁场响应对比(英文)

电偶极子依靠两个等量异号电荷的加速运动实现电磁场辐射,同样地,时变点电荷微元可以依靠单个电荷的加速运动产生电磁场。基于这种电磁场理论实质,本文开展了点电荷微元与偶极子源的时域电磁场响应对比研究。首先,给出三维导电空间点源的时间域位函数;然后,利用辅助位函数与电场、磁场响应之间的关系,给出均匀无耗全空间的点电荷载流微元时间域电场、磁场的表达式。并从阶跃关断时刻、近源区场、远源区场等方面对比分析全空间中两种源的异同。研究结果表明:在近源区,两种源产生的场差别较大;在远源区,两种源产生的场几乎相同。通过实际观测的数据和使用电偶极子源公式对简化模型计算结果的对比,认为点电荷微元场的空间分布更符合实际场源产生的电场、磁场分布特征。     

ATI介质中四次时差系数解析近似

非双曲（远偏移距）时差为各向异性介质中正、反演研究,特别是各向异性参数估计提供了重要信息.本文在任意空间取向TI（ATI）介质水平界面同类反射波四次时差系数(A₄)精确解的基础上,进一步讨论我们推导得出的ATI介质中四次时差系数解析近似解,比较随CMP测线方位变化的近似解与精确解之间的差别,为利用近似解来解析研究ATI介质中非双曲时距以及参数反演提供有价值的信息．结合实际岩性资料的数值研究表明,ATI条件下四次时差系数近似解与精确解之间存在差别,不仅表现在A₄数的大小及符号特征上,更突出地表现在A₄系数随方位的变化特征上;在强各向异性条件下,近似解相比精确解存在较大误差．但在各向异性参数满足0<ε-δ<0.15、｜δ｜<0.20的情况下,对于TI对称轴的特殊倾角范围(75°～80°),近似解与精确解的差别很小,可用近似解进行各向异性观测解释及参数反演．     

CHAMP卫星时代全球岩石圈磁场模型的近似误差估计

岩石圈磁场模型的误差估计是地磁导航、定量地球物理与地质解释的重要基础.本文基于统计学观点,将由CHAMP卫星磁测数据构建的系列具有代表性的全球岩石圈磁场模型视为随机观测信号,对其在球谐域与空间域分别进行平均值与标准差的统计分析,得到CHAMP卫星时代全球岩石圈磁场模型的球谐域与空间域近似误差水平、误差分布特征、不同类型磁场分量之间的误差比例关系以及全球岩石圈磁场能量与误差水平随高度的衰减特征.本文所提的误差估计方法与结果,尽管存在一定的近似,但是可以为估计与了解岩石圈磁场模型在球谐域与空间域的误差水平与分布特征提供参考资料与定量依据.     

考虑雷暴活动的全球近地稳态大气电场的精确数学解

全球大气电模式是指大气和地表构成的闭合回路,其中考虑雷暴为主要充电电源,并将雷暴作为一个垂直电偶极子来考虑.本文从电流连续方程出发,采用球谐函数展开法,边界条件采用下边界为地表大气电势,上边界为电离层底部电势,非齐次电流源函数项采用垂直偶极子假设.通过所做的非线性函数变换,得出了考虑雷暴活动的全球电模式的近地稳态解析解,并证明其能展开为有限多项式级数之和.所得结论与过去的理论分析结果和实验观测是一致的.该文的工作是对以前相关工作仅能得到近似解的一种改进,同时本文所得公式可用于今后对大气电参量进行相关的数值计算及模拟之中.     

梯度法在球谐综合中的适用性分析

针对地球表面球谐综合过程中要求考虑高程信息的问题,在球谐综合过程中,基于地形模型DTM2006.0和超高阶的重力场模型EGM2008,分别在椭球坐标系与球坐标系下,利用梯度法对球谐综合过程中的重力场泛函进行一定阶次扩展,分析梯度法在椭球坐标系与球坐标系下的适用性.在实验区域,计算结果表明,在球谐综合过程中,利用梯度法计算重力场泛函,它能够将高程信息带来的近似误差减少到对实际应用无关紧要的程度,可计算任意密集网格点上的表面重力场泛函.     

简谐横荷载作用下Winkler地基上有限长梁1/3次亚谐共振分析

基于Winkler地基模型及Euler - Bernoulli梁理论,建立了Winkler地基上有限长梁的非线性运动方程.运用Galerkin方法对运动方程进行一阶模态截断,得到离散的非线性振动方程,利用多尺度法求得该系统1/3次亚谐共振的一阶近似解.分析了长细比、弹性模量、地基刚度、阻尼、密度等参数对其亚谐共振幅频响...     

Comparison of the time-domain electromagnetic field from an infinitesimal point charge and dipole source

An electromagnetic field is generated through the accelerating movement of two equal but opposite charges of a single dipole. An electromagnetic field can also be generated by a time-varying infinitesimal point charge. In this study, a comparison between the electromagnetic fields of an infinitesimal point charge and a dipole has been presented. First, the time-domain potential function of a point source in a 3D conductive medium is derived. Then the electric and magnetic fields in a 3D homogeneous lossless space are derived via the relation between the potential and field. The field differences between the infinitesimal point charge and the dipole in the step-off time, far-source, and near-source zones are analyzed, and the accuracy of the solutions from these sources is investigated. It is also shown that the field of the infinitesimal point charge in the near-source zone is different from that of the dipole, whereas the far-source zone fields of these two sources are identical. The comparison of real and simulated data shows that the infinitesimal point charge represents the real source better than the dipole source.     

可控源电磁三维频率域有限元模拟

本文采用电磁场的磁矢量位和电标量势,将Maxwell方程组化为位势的类似于Helmholtz型方程,并引入罚项及稳定化方法克服了电磁三维有限元算法中的伪解及数值不稳定性;采用人工边界把计算区域局域化,将均匀半空间中水平电偶极子源产生的位势值作为人工边界上的第一类边界条件以表示源的作用,减少了实际的计算区域.理论模型和复杂模型的计算结果均表明,可控源电磁三维有限元数值模拟给出了稳定、可靠的电磁场分布.     

On the electrical reflectivity tensor in d.c. electric field modelling

We use the integral equation for a d.c. electric field, published in the literature, to introduce the concept of the electrical reflectivity tensor into d.c. electric field modelling. It is shown that in d.c. electric field modelling, the electric reflectivity tensor can be obtained in exactly the same way as in electromagnetic modelling. As a result, for a d.c. electric field, the quasi‐linear and the quasi‐analytical approximations, as well as the quasi‐analytical series, can be constructed in exactly the same way as in electromagnetic modelling. If the primary field is uniform, and if the anomalous body is a uniform circular cylinder or a uniform sphere, the reflectivity tensor is zero order (constant), relating to the free surface charge density. Thus, for some homogeneous bodies that have simple shapes and are embedded in a uniform primary field, the electrical reflectivity tensor is not only a mathematical mechanism for obtaining approximate solutions, but also a physical reality. Indeed, the free surface charge density is defined as the change of the electric displacement vector across the boundary surface under consideration. If the primary field is caused by a point source, and if the anomalous body is a uniform sphere, the reflectivity tensor is second order, varying slowly within the sphere. The relationship to the free surface charge density can be established only when both the reflectivity tensor and the free surface charge density are approximated by the first terms of their series solutions. If the point source is far from the centre of the sphere, the corresponding reflectivity tensor reduces to zero order, and is independent of the observation position within the sphere, i.e. it is a constant. Therefore, the basic idea of the quasi‐analytical approximation, i.e. taking the reflectivity tensor outside the integral operator, is justified in the case considered here.     

基于电场矢量波动方程的三维可控源电磁法有限单元法数值模拟

从可控源电磁法的基本原理出发,推导了基于电场矢量波动方程的三维边值问题,利用广义变分原理,把边值问题转换为变分问题,并引入散度条件,避免了伪解的出现,使有限元计算在理论上更加完备.在准静态近似条件下,把水平电偶极子在空中和大地的远区电场闭合表达式作为有限元计算中的区域外边界条件,解决了边界条件加载的困难;把应用于地震模拟中的伪delta函数引入到可控源电磁法中的三维有限元模拟中,消除了源点的奇异性,提高了方程组的稳定性.通过对均匀大地和层状介质模型的模拟,检验了程序的正确性,并对典型的地质体模型进行了数值模拟,分析了其变化规律.     

单一层面高频电偶极子场的近似计算

本文对一个倾斜平面界面的高频电磁波场近似计算问题做了解答。这是钻井高频电磁波勘探方法中的基本问题。文中写出了问题的严格解,并用鞍点法对φ=0的条件,得出其近似计算式和数值结果.这个结果具有便于计算和明确的物理意义。根据计算结果,指出了特别在高导介质中,界面附近的干涉异常,必须予以充分重视。     

On small-scale near-surface distortion in controlled source tensor electromagnetics1

Based on a Born approximation of a thin sheet integral equation, it is shown that small-scale surficial heterogeneity significantly distorts the electromagnetic field excited by electric dipoles only when either the source or the receiver are located on the heterogeneities. When a surface heterogeneity is beneath the source, the associated distortion of the electromagnetic field is manifest as a change in the effective electric dipole moment. Hence the magnetotelluric transfer functions and impedance relations remain undistorted in this case. When a surface heterogeneity is beneath the receiver, the electric field is severely distorted, but the magnetic field is only slightly distorted. The impedance tensor is therefore strongly distorted, but the tipper vector is almost unaltered. Since the controlled source tipper is a function of 1D earth conductivity, it is proposed that tipper data should be used in the first stage of 1D interpretation. For a 1D earth, the tipper vector must always point towards the source and, in the near-field limit, should have unit length. These two necessary conditions must be met by the measured tipper before it is interpreted one dimensionally.     

Contribution of the relative velocity between source and receivers to electromagnetic fields in the sea

The analysis of the primary electromagnetic fields caused by steady state or transient electric current flowing along a current loop moving with a constant velocity below the sea surface has several applications. It supports the analysis of submarine physical data and it is useful for protecting ships from the threat of sea mines. The usual approach to the solution for the primary field starts from a hertz vector potential in the frequency domain due to a magnetic dipole. Subsequently it employs Fourier, Laplace, and Hankel transforms to describe the time variation of the primary electromagnetic induced field due to a loop. The result is applicable to both shallow and deep sea water environments. Because of the difference in velocity between source and receiver, a careful application of the convolution integral is necessary in order to adapt the source pulse solution to any type of transmitting current waveform. Furthermore, since the scattered field represents a fraction of the primary field, even minor differences in it caused by the differential velocity renders inadequate interpretation of EMI data.     

虚拟波动域三维海洋可控源电磁场正演模拟

本文基于对应原理将似稳态条件下频率域电磁场扩散方程转换成虚拟波动域电磁场波动方程,采用高阶时域有限差分进行求解,引入复频移完全匹配层吸收边界条件,降低了内存需求,提高了计算效率,并在虚拟波动域用伪δ函数离散电偶极源,实现了虚拟波动域任意取向电偶极源三维海洋可控源电磁场高阶时域有限差分正演算法.通过与拟解析解和频率域三维可控源电磁场数值模拟结果的对比,验证了本文算法的正确性和高效性,且探讨了网格参数和边界条件对不同频率电磁场模拟结果的影响.