面向目标自适应有限元法的带地形三维大地电磁各向异性正演模拟

传统三维大地电磁各向异性模拟均是基于规则六面体网格,计算精度有限且较难拟合复杂地质条件.本文采用面向目标自适应非结构矢量有限元法,对三维大地电磁各向异性介质进行模拟.首先从电场双旋度方程出发,利用伽辽金方法建立变分方程;然后利用电流密度连续性条件构建适合大地电磁各向异性问题的加权后验误差估计方法,实现面向目标的网格自适应正演;最后通过典型算例分析各向异性对网格自适应和大地电磁响应的影响特征以及各向异性的识别方法.本文算法能够高精度地拟合起伏地表和任意各向异性介质,适用于分析复杂地电条件大地电磁响应特征,为提高大地电磁资料解释水平提供了理论基础.     

电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟

现有海洋可控源电磁三维数值模拟方法大多基于电导率各向同性介质理论,不能模拟海底地层电导率各向异性的实际情况.本文给出了电导率各向异性三维介质中电性源海洋可控源电磁二次电场的边值问题以及相应的变分问题,采用长方体单元对研究区域剖分,将场分量定义在剖分单元的边上,利用矢量有限单元法求解变分问题,实现了电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟.这个新的正演方法可以计算电导率任意各向异性三维地电模型的海洋可控源电磁响应,基于二次场矢量有限元法直接求解电磁场,避免了传统有限元方法可能遇到的伪解问题和难于处理电场法向分量不连续的问题,提高了数值模拟计算精度.一维电导率各向异性模型电磁场数值解与解析解吻合得相当好,无论在源附近还是远离源处相对误差均不超过1%.电导率各向异性二维模型的计算结果与已有文献采用的非结构有限元模拟结果十分吻合.三维地电模型数值模拟结果显示,电导率各向异性张量电导率主轴分量和欧拉角对不同装置海洋可控源电磁响应均有着明显的影响.     

Anisotropy of magnetic remanence: A brief review of mineralogical sources,physical origins,and geological applications,and comparison with susceptibility anisotropy

The magnetic fabric of rocks and sediments is most commonly characterized in terms of the anisotropy of low-field magnetic susceptibility (AMS). However, alternative methods based on remanent magnetization (measured in the absence of a magnetic field) rather than induced magnetization (measured in the applied field) have distinct advantages for certain geological applications. This is particularly true for; (1) adjunct studies in paleomagnetism, in order to assess the fidelity with which a natural remanence records the paleofield orientation; (2) studies of weakly magnetic or weakly deformed rocks, for which susceptibility anisotropy is very difficult to measure precisely; and (3) quantitative applications such as strain estimation. The fundamental differences between susceptibility and remanence (and their respective anisotropies) are due to several factors: (1) susceptibility arises from all of the minerals present in a sample, whereas remanence is carried exclusively by a relatively small number of ferromagnetic minerals; (2) ferromagnetic minerals are generally more anisotropic than para- and diamagnetic minerals; (3) for ferromagnetic minerals, remanence is inevitably more anisotropic than susceptibility; and (4) a number of common minerals, including single-domain magnetites, possess an inverse anisotropy of susceptibility, i.e., they tend to have minimum susceptibility parallel to the long axis of an individual particle; remanence is immune to this phenomenon. As a consequence of all these factors, remanence anisotropy may generally provide a better quantitative estimate of the actual distribution of particle orientations in a rock sample.Contribution number 9102 of the Institute for Rock Magnetism, University of Minnesota.     

带地形的大地电磁各向异性二维模拟及实例对比分析

带地形的大地电磁二维正演数值模拟多数基于电性各向同性理论,由于地球内部电性各向异性现象的普遍存在,基于电性各向异性理论研究地形起伏情况下大地电磁二维正演数值模拟就显得非常迫切.本文首先由麦克斯韦方程出发,引入张量电导率,求得一组关于平行走向的电场分量Ex和磁场分量Hx的二阶偏微分方程,使用有限差分法求解出Ex和Hx的近似解,并以此求得其他场分量;其次,引入地形因素,改变变量在网格节点中的排列方式,选择交错排列方式从而给有限差分系数矩阵的最大带宽分配合理的存储空间;最后,使用Weaver的方法解决TM模式下,在地-空分界面垂直于构造走向的一些区域存在不同电导率的问题.通过对带地形的二维电性各向异性结构做正演模拟,研究地形因素对大地电磁响应的影响;以电性各向异性理论为基础,将地形因素引入对实测大地电磁资料的处理中,通过做二维正演拟合和未引入地形因素的结果做对比,说明电性各向异性现象的普遍存在,认识地形因素对观测大地电磁场的影响,为今后分析解释实测大地电磁资料包含地形因素和电性各向异性情况提供理论基础和技术指导.     

安宁河断裂带测区古构造残余应力随深度分布

以岩体正交异性弹性理论为基础，用Ｘ射线法，在安宁河断裂带测区，选取了９个钻孔的岩芯，测量了岩体中宏观残余应力场的水平和铅直三维主应力及最大剪和随深度的分布。分析了其沿深度的分布规律，求得了其铅直梯度值。     

回线源瞬变电磁法有限体积三维任意各向异性正演及分析

目前,瞬变电磁法(TEM)数据基本都是基于各向同性模型进行反演解释,这对于存在明显电性各向异性的勘探区域会产生较大的反演解释误差.为分析电各向异性对回线源瞬变电磁信号的影响方式与程度,本文通过求解离散化的全张量电导率时间域Helmholtz方程,实现了基于有限体积法的TEM任意各向异性的三维正演算法.该算法采用基于交错网格的拟态有限体积法(MFV)对时域Maxwell方程组进行空间域离散,并利用后退欧拉算法(Backward Euler Method)进行时间域离散.为提高时域电磁场的求解精度与效率,该算法将时间分段等步长算法与方程直接求解法相结合.通过对一维各向异性模型以及三维复杂各向同性模型进行测试,验证了本算法对于回线源瞬变电磁响应计算的正确性及有效性.最后,通过对几类典型电各向异性介质中大回线源瞬变电磁信号响应的分析,总结了不同电各向异性类型对TEM电磁信号的影响模式,结果表明,主轴各向异性情况下TEM信号主要受水平方向电导率的影响,倾斜各向异性对TEM信号的影响程度远大于水平各向异性,而通过水平各向异性信号能较清晰判断出各向异性主轴方向.     

基于牛顿迭代法求解群速度的地震波射线追踪模拟

地下介质中普遍存在着各向异性,当前基于各向异性的地震波射线追踪多是在弱各向异性介质中进行且采用群速度近似表示方法,这些近似方法在强各项异性介质中会导致很大误差而无法真正模拟地震波的传播规律。根据地下普遍存在各向异性的事实和地震波基本传播规律,提出利用牛顿迭代法高效求解群速度,基于Paraview平台自动化构建三维地质模型,采用最短路径法进行地震波射线追踪模拟及可视化,实现对复杂三维地质的速度不均匀性和各向异性的表达,为三维地质模型的构建和地震波射线追踪模拟及可视化提供一种新思路,并以华北克拉通山西断陷带北部局部区域为例进行研究。结果表明,该方法能够减少由各向异性对地震波传播模拟造成的影响,清晰表达了研究区地质结构和各向异性特点,在对复杂三维地质结构的解读中能够较好应用。     

轴对称各向异性介质波动方程深度偏移的对称非平稳相移方法

由所建立的三维qP波相速度表示式出发,导出并解析求解各向异性介质中的频散方程,得到三维各向异性介质中的相移算子,进而将以相移算子为基础的对称非平稳相移方法推广到各向异性介质,发展了一个三维各向异性介质的深度偏移方法. 文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致,将各向同性、弱各向异性及强各向异性统一在一个模型中. 所建立的各向异性介质对称非平稳相移波场延拓算子可以同时适应速度及各向异性参数横向变化;文中给出的算例虽然是针对二维VTI介质的,但所提出的算法同样适用于三维TI介质.     

Influence of magnetic interactions on the anisotropy of magnetic susceptibility: the case of single domain particles packed in globule aggregates

The influence of magnetic interactions on the anisotropy of magnetic susceptibility (AMS) have been largely studied by several theoretical models or experiments. Numerical models have shown that when magnetostatic interactions occur, the distributions of particles over the volume rather than their individual orientations control the AMS. We have shown recently from a comprehensive rock magnetic study and from a theoretical 2-dimensional (2-D) model that single domain particles closely packed in globule aggregates could produce strong local random interaction magnetic fields which could influence the magnetic susceptibility and decrease the degree of anisotropy. In this paper, we first present in detail this 2-D theoretical model and then we extend it to the 3-D case. The possible distribution function of the magnetostatic interaction fields comprises two extreme states: it is either isotropic or ordered. The former case corresponds to the thermal-demagnetized state while the second case corresponds to the alternating field (AF) demagnetized state. We show that when easy axes of magnetization are not uniformly distributed, the degree of anisotropy decreases as the interaction field increases in both AF- and thermal-demagnetized states in 2-D and 3-D geometry. Thus we conclude that random magnetic fields generated by a random arrangement of magnetic particles over the sample volume decrease the degree of anisotropy of AMS and may alter the magnetic fabric.     

Magnetotelluric inversion for azimuthally anisotropic resistivities employing artificial neural networks

An extension of an artificial neural network (ANN) approach to solve the magnetotelluric (MT) inverse problem for azimuthally anisotropic resistivities is presented and applied for a real dataset. Three different model classes, containing general 1-D and 2-D azimuthally anisotropic features, have been considered. For each model class, characteristics of three-layer feed forward ANNs trained through an error back propagation algorithm have been adjusted to approximate the inverse modeling function. It appears that, at least for synthetic models, reasonable results would be obtained by applying the amplitudes of the complex impedance tensor elements as inputs. Furthermore, the Levenberg-Marquart algorithm possesses optimal performance as a learning paradigm for this problem. The evaluation of applicability of the trained ANNs for unknown data sets excluded from the learning procedure reveals that the trained ANNs possess acceptable interpolation and extrapolation abilities to estimate model parameters accurately. This method was also successfully used for a field dataset wherein anisotropy had been previously recognized.     

频率域海洋可控源电磁垂直各向异性三维反演

地层宏观电性各向异性会对可控源电磁响应产生重要影响.由于海底地层电性结构常表现为电导率各向异性,若仅对海洋可控源电磁(MCSEM)数据进行常规各向同性反演,有可能无法获得准确的反演解释结果,从而削弱MCSEM技术的可靠性.本文实现了电导率垂直各向异性(VTI)条件下频率域海洋可控源电磁数据三维反演算法.其中,三维正演采用基于二次场控制方程的交错网格有限体积法,并利用直接矩阵分解技术来求解离散所得的大型线性方程组,有利于快速计算多场源的响应.反演采用具有近似二次收敛性的高斯牛顿算法对目标函数进行最优化.最后,对具有VTI电性各向异性特征的盐丘构造模型的MCSEM合成数据分别进行了电导率各向同性和垂直各向异性三维反演,结果表明:各向同性三维反演算法无法对受VTI介质影响的MCSEM数据进行正确的反演解释,而垂直各向异性三维反演能够获得更为可靠的地下电阻率结构和异常体分布,展现出对海底电性各向异性结构更为优良的反演解释能力.     

起伏地形下各向异性的2D大地电磁无网格法数值模拟

作为一种基于节点的计算方法,无网格法具有构造高阶导数方便,自适应分析便利的优点,特别适合复杂地质构造的数值模拟.本文针对实际地球物理勘探中存在的起伏地形和各向异性的地电结构,提出用无网格法来模拟大地电磁响应,采用复合二次径向基函数构造了形函数,推导了大地电磁无网格法等价线性方程组,研究了系数矩阵的压缩存储方法以及大型稀疏复线性方程组快速求解算法,实现了起伏地形下各向异性的2D大地电磁高精度数值模拟.基于层状模型验证了算法的正确性,计算结果表明:无论是TM模式还是TE模式,计算相对误差均小于1%;通过对地垒和地堑模型的模拟,得出了起伏地形对视电阻率和相位的影响规律;对起伏地形下含有不同各向异性系数异常体的模型进行了数值模拟,为开展复杂地质结构的电磁场特征研究以及地形校正奠定了理论基础.随着计算科学的快速发展,无网格法必将成为新的高精度电磁场数值模拟方法.     

用有限差分法计算各向异性介质中多分量感应测井的响应

感应测井电阻率是砂泥岩储层中定量评价含油气性的主要参数之一,然而由于储层电各向异性的存在,传统的感应测井仪得到的电阻率主要响应地层的水平电阻率,比地层的垂直电阻率小,由此计算的含水饱和度偏高．多分量感应测井则可以通过不同方向的发射和接收,测量多个分量,弥补上述缺陷,为储层的含油气性评价提供准确信息．本文用交错网格有限差分法导出了主轴各向异性介质中Maxwell方程的离散化关系式,计算了三维各向异性介质中多分量感应测井的电磁响应．利用三层模型的数值模拟,分析了不同磁场分量对各向异性介质的响应特性,用多层模型和倾斜各向异性地层的响应计算检验了所编制的交错网格有限差分法程序．     

三维各向异性介质中的波动方程叠前深度偏移方法

基于三维VTI各向异性介质的频散关系,构建波数项和空间项分离的单程波算子表达式,以优化算法,确定算子的待定系数,实现广角逼近三维VTI介质的广义相移算子,发展了可灵活处理强或弱各向异性介质的波动方程叠前深度偏移方法.文中同时也针对其工业应用建议了三维VTI各向异性介质中可提高计算效率的频率相关变步长波场深度延拓算法及稀...     

电导率各向异性的海洋电磁三维有限单元法正演

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.     

3D anisotropic modeling and identification for airborne EM systems based on the spectral-element method

The airborne electromagnetic (AEM) method has a high sampling rate and survey flexibility. However, traditional numerical modeling approaches must use high-resolution physical grids to guarantee modeling accuracy, especially for complex geological structures such as anisotropic earth. This can lead to huge computational costs. To solve this problem, we propose a spectral-element (SE) method for 3D AEM anisotropic modeling, which combines the advantages of spectral and finite-element methods. Thus, the SE method has accuracy as high as that of the spectral method and the ability to model complex geology inherited from the finite-element method. The SE method can improve the modeling accuracy within discrete grids and reduce the dependence of modeling results on the grids. This helps achieve high-accuracy anisotropic AEM modeling. We first introduced a rotating tensor of anisotropic conductivity to Maxwell’s equations and described the electrical field via SE basis functions based on GLL interpolation polynomials. We used the Galerkin weighted residual method to establish the linear equation system for the SE method, and we took a vertical magnetic dipole as the transmission source for our AEM modeling. We then applied fourth-order SE calculations with coarse physical grids to check the accuracy of our modeling results against a 1D semi-analytical solution for an anisotropic half-space model and verified the high accuracy of the SE. Moreover, we conducted AEM modeling for different anisotropic 3D abnormal bodies using two physical grid scales and three orders of SE to obtain the convergence conditions for different anisotropic abnormal bodies. Finally, we studied the identification of anisotropy for single anisotropic abnormal bodies, anisotropic surrounding rock, and single anisotropic abnormal body embedded in an anisotropic surrounding rock. This approach will play a key role in the inversion and interpretation of AEM data collected in regions with anisotropic geology.     

Quantitative measures of textural anisotropy resulting from magmatic compaction illustrated by a sample from the Palisades sill, New Jersey

Early in the crystallization of many tholeiitic basaltic magmas, plagioclase crystals cluster together into a 3-D cellular network, which forms a passive marker capable of recording the deformation that accompanies compaction of crystal mush. Although irregular in detail, the overall network is initially isotropic and only becomes anisotropic as a result of compaction. We have developed four independent methods to quantify the 3-D textural anisotropy of a basalt sample using at least three non-parallel thin sections. Three of the methods are based on the geometrical properties of digitized maps of the feldspar chain networks. One approach focuses on the angular variation of the mean intercept along parallel traverses through the network, another examines the orientation and size distribution of individual links, and the third considers the average shape of interstitial regions outlined by the plagioclase network. The fourth technique approximates the textural anisotropy by the variogram anisotropy of a scanned thin section image. We illustrate the methods using five oriented non-parallel thin sections from a sample of diabase 146 m above the base of the 300-m-thick Palisades sill of New Jersey. Compaction of crystal mush in this sill has previously been postulated on the basis of chemical evidence. The 3-D feldspar network anisotropy based on the first three approaches suggests nearly uniaxial compaction on the order of 8.6% in a direction within 3° of the intersection of the columnar joints at the sample site. A rigorous statistical test based on the statistics of elliptically contoured non-normal multivariate distributions documents that the link-vector distribution in vertical sections are statistically anisotropic at a 95% confidence level and that the overall compaction is 7.9±2.6%. The orientation and magnitude of the 3-D textural anisotropy determined by the image variogram of the non-opaque minerals is almost identical to the mean feldspar network anisotropy; 8.5% compaction in a direction 10° from the columnar joint intersections. The major silicate textural and feldspar network anisotropy axes both plunge almost directly down dip of the sill. On the other hand, the major axis of the variogram anisotropy of the opaque minerals is approximately parallel to the strike of the sill and to the major axis of the anisotropic magnetic susceptibility. The anisotropy of the silicate mineral fabric may reflect down-dip flow of a deformable melt-rich crystal mush, whereas the AMS and opaque textural anisotropy reflects the influence of gravitational stresses during the growth of magnetite in the final stages of melt crystallization. Evidently the Palisades sill was not originally horizontal but was intruded in an orientation close to its present attitude.     

各向异性标准化方差计算重磁源边界

在重磁源边界定位方法中,传统的梯度方法易受干扰的影响使计算的边界混乱,而且在弱异常处由于叠加异常的影响很难识别场源边界.本文首先利用坐标旋转构造了各向异性高斯函数,提出了各向异性标准化方差计算重磁源边界的方法.理论分析与模型实验详细阐明了该方法的数学含义,并通过干扰分析验证了方法的稳定性与有效性,结合中扬子地区航磁异常...     

Seismoelectric numerical simulation in 2D vertical transverse isotropic poroelastic medium

Seismoelectric coupling in an electric isotropic and elastic anisotropic medium is developed using a primary–secondary formulation. The anisotropy is of vertical transverse isotropic type and concerns only the poroelastic parameters. Based on our finite difference time domain algorithm, we solve the seismoelectric response to an explosive source. The seismic wavefields are computed as the primary field. The electric field is then obtained as a secondary field by solving the Poisson equation for the electric potential. To test our numerical algorithm, we compared our seismoelectric numerical results with analytical results obtained from Pride's equation. The comparison shows that the numerical solution gives a good approximation to the analytical solution. We then simulate the seismoelectric wavefields in different models. Simulated results show that four types of seismic waves are generated in anisotropic poroelastic medium. These are the fast and slow longitudinal waves and two separable transverse waves. All of these seismic waves generate coseismic electric fields in a homogenous anisotropic poroelastic medium. The tortuosity has an effect on the propagation of the slow longitudinal wave. The snapshot of the slow longitudinal wave has an oval shape when the tortuosity is anisotropic, whereas it has a circular shape when the tortuosity is isotropic. In terms of the Thomsen parameters, the radiation anisotropy of the fast longitudinal wave is more sensitive to the value of ε, while the radiation anisotropy of the transverse wave is more sensitive to the value of δ.     

基于有限体积法的二维大地电磁各向异性数值模拟

为了计算带任意地形的各向异性介质中二维大地电磁响应,本文在非结构化网格的基础上,采用有限体积法,开发了二维大地电磁各向异性正演模拟的新算法.首先,从Maxwell方程出发,推导二维各向异性介质中大地电磁场的边值问题;然后,采用三角网格自动生成技术对求解区域进行非结构化网格剖分,进而构建节点中心控制体积单元,利用有限体积方法,得到求解边值问题的大型稀疏线性方程组;最后,利用Pardiso精确地计算了大地电磁响应值.三个各向异性模型的计算结果表明,本文开发的有限体积算法,不仅能够高精度求解带任意地形的大地电磁电导率各向异性问题,而且对于同一模型,该方法的计算消耗和精度都与有限单元法相当.因此,有限体积法是处理电磁法各向异性问题的一种有效方法.