共查询到20条相似文献，搜索用时 46 毫秒
1.
Zhi-Yong Zhang Han-Dong Tan Kun-Peng Wang Chang-Hong Lin Bin Zhang Mao-Bi Xie 《应用地球物理》2016,13(1):13-24
Traditional two-dimensional (2D) complex resistivity forward modeling is based on Poisson’s equation but spectral induced polarization (SIP) data are the coproducts of the induced polarization (IP) and the electromagnetic induction (EMI) effects. This is especially true under high frequencies, where the EMI effect can exceed the IP effect. 2D inversion that only considers the IP effect reduces the reliability of the inversion data. In this paper, we derive differential equations using Maxwell’s equations. With the introduction of the Cole–Cole model, we use the finite-element method to conduct 2D SIP forward modeling that considers the EMI and IP effects simultaneously. The data-space Occam method, in which different constraints to the model smoothness and parametric boundaries are introduced, is then used to simultaneously obtain the four parameters of the Cole—Cole model using multi-array electric field data. This approach not only improves the stability of the inversion but also significantly reduces the solution ambiguity. To improve the computational efficiency, message passing interface programming was used to accelerate the 2D SIP forward modeling and inversion. Synthetic datasets were tested using both serial and parallel algorithms, and the tests suggest that the proposed parallel algorithm is robust and efficient. 相似文献
2.
本文基于非结构网格实现了海洋可控源电磁法三维有限元正演模拟.该算法采用完全非结构网格剖分,可以模拟任意起伏地形和复杂地电模型.为了避免场源的奇异性,采用一次场/二次场分解算法,一次场由基于Schelkunoff势函数的一维解析公式得到.为了提高算法的精度和效率,采用对测点附近单元和异常体区域进行体积约束加密的方法,实现了非结构网格的局部加密.一、二维模型计算和分析表明,本文采用的局部加密方法能够明显地改善算法的精度,最大相对误差基本在1%以内.对三维模型计算及对比分析,说明了该算法对三维可控源电磁正演的实用性.复杂海底地形模型的正演模拟表明,海底地形对电磁场的影响很大,在进行海洋可控源电磁资料解释时,地形的影响有必要考虑在内. 相似文献
3.
4.
Seafloor massive sulfide (SMS) deposits are generated by high-temperature hydrothermal systems. Their precious resources have attracted global interest. A number of investigations with controlled-source electromagnetic (CSEM) methods have been implemented in recent years. There are three major problems with SMS surveying using EM methods. First, SMS imaging techniques for hydrothermal systems have a limited range. Simulations and applications have validated only simple layered models. Second, their inversion efficiencies must be improved further. Laterally constrained inversions and spatially constrained inversions are usually used to map geological structures. However, choosing their suitable weighting parameters is inefficient. Third, the effects of induced polarization (IP) on ore deposits are not considered in such inversions. A non-polarizable model is unable to accurately depict a polarizable model. To resolve these problems, an advanced strategy is used to improve the efficiency of the pseudo-3D inversion process. The proposed imaging method has the ability to map complex 3D geoelectrical structures, and therefore, it can both obtain information regarding surface ore deposits and distinguish between active and inactive hydrothermal systems. However, this method can also be used to depict the distributions of alteration zones and buried deposits. Furthermore, the influences of IP on the inversion are discussed with respect to the Cole- Cole model, and it is shown that the effects of IP on polarizable deposits cannot be ignored during the inversion. 相似文献
5.
本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究. 相似文献
6.
本文采用有理函数Krylov子空间模型降阶算法实现了同时求解多频可控源电磁法三维正演响应的快速计算.首先采用基于Yee氏交错网格的拟态有限体积法实现控制方程的空间离散,将任意频率的电场响应表示为关于频率参数的传递函数.采用有理函数Krylov子空间算法求解该传递函数.针对构建m维有理函数Krylov子空间需要求解m次（几十到上百）关于有理函数极点和离散控制方程系数矩阵的线性方程组的问题,本文提出采用单个重复极点的有理函数Krylov子空间模型降阶算法,结合直接法求解器PARDISO,采用Gram-Schmidt方法,只需要1次系数矩阵分解和m次矩阵回代即可实现有理函数Krylov子空间的构建,极大地减少了计算量.针对最优化有理函数极点选取问题,本文根据传递函数的有理函数Krylov子空间投影算法的误差分析理论,引入关于单个重复极点的收敛率函数,通过求解有理函数的最大收敛率直接给出最优化的单个重复极点公式.最终实现了不同发射频率的可控源电磁法三维正演响应的快速计算.分别计算了典型层状模型多发射频率的CSAMT和海洋CSEM的正演响应,通过与解析解的对比验证了本文算法在多发射频率正演的计算精度和计算效率;并通过一个三维海洋CSEM勘探设计最优化发射频率和接收区域选取的例子进一步说明本文算法的优点. 相似文献
7.
复电阻率法在矿产、油气勘探调查中发挥着重要作用.为了认识复杂构造的复电阻率法电磁场的变化规律,本文基于自适应有限元方法,采用非结构化网格,引入Cole-Cole模型,实现了电偶源2.5D复电阻率法电磁场正演,可以模拟复杂地形和地电结构,正演结果更符合野外实际地质情况.通过将本文的计算结果与半空间模型解析解、层状介质和起伏模型结果进行对比,验证了本文算法的正确性.最后,基于复杂地电模型,通过正演模拟,系统分析了地形、激电参数、复杂构造对复电阻率法电磁场的影响特征. 相似文献
8.
从电偶源三维地电断面可控源电磁法的二次电场边值问题及其变分问题出发,采用任意六面体单元对研究区域进行剖分,并且在单元分析中同时对电导率及二次电场进行三线性插值,实现电导率分块连续变化情况下,基于二次场的可控源电磁三维有限元数值模拟.这个新的可控源电磁三维正演方法可以模拟实际勘探中地下任意形状及电性参数连续变化的复杂模型.理论模型的计算结果表明,均匀大地计算的视电阻率误差和相位误差分别为0.002%和0.0005°.分层连续变化模型的有限元计算结果表明,其与对应的分层均匀模型解析结果有明显差异.三维异常体组合模型以及倾斜异常体等复杂模型的有限元计算结果也有效地反映了异常形态. 相似文献
9.
可控源电磁法具有分辨率高及抗干扰能力强等特点,是一种重要的地电磁勘探方法.目前,可控源电磁法的高精度正演计算一直是其核心研究问题之一.传统积分方程法一般采用近似积分公式、简单矩形网格和近似的奇异性体积分计算技术,制约了体积分方程法处理复杂地下异常体的能力,降低了计算精度.针对上述问题,本文基于完全积分公式、四面体非结构化网格和奇异体积分的精确解析解来高精度求解复杂可控源电磁模型的正演响应.首先,从电场积分公式出发,推导了可控源电磁问题满足的积分方程;其次,借助于非结构化四面体网格离散技术,实现了地下复杂异常体的有效模拟.最后,利用散度定理把强奇异值体积分转换为一系列弱奇异性的面积分公式,并通过推导获得了这些弱奇异性的面积分公式的解析解,从而最终实现三维可控源电磁问题的高精度积分求解.以块状低阻体地电模型为测试模型,采用本文提出的积分方程方法获得的数值解与其他公开数值算法解进行对比分析,其对比结果具有高度的吻合性,验证了算法的正确性;同时,设计了球状及复杂地电模型进行算法收敛性测试,进一步验证算法的正确性以及能够处理地下复杂模型的能力. 相似文献
10.
大规模地球物理电磁数据的定量解释需要发展高效、稳定的三维正反演算法.本文通过求解离散化的三维电场矢量Helmholtz方程,实现了基于有限体积法的频率域可控源电磁(CSEM)三维正演算法.为模拟具有强电性差异的三维电性介质,该算法采用拟态有限体积法(MFV)对Maxwell方程组进行离散化;另外,为获得稳定、高精度的正演数值结果,采用直接矩阵分解技术来求解离散所得到的大型稀疏线性方程组.对于具有多个发射源的CSEM测量来说,一次矩阵分解结果能够用于同频率下所有场源的正演计算.为降低场源奇异性及边界条件对数值精度的影响,采用虚拟场源校正技术,避免了散射场公式中在构建场源项时所需的大量时间.对于具有多个频率的CSEM的模拟计算,采用分频并行策略来加快三维正演计算.最后,通过与一维层状模型及三维模型的数值结果的对比验证了本文所开发的正演算法对频率域CSEM模拟计算的准确性及有效性,表明该正演算法能够有效应用于三维介质的数值计算.另外,对于多频率CSEM的并行测试结果表明基于分频并行策略的并行计算能够显著地降低正演计算时间. 相似文献
11.
12.
为了分析海洋中含金属矿的储油气砂层的成分、结构对于海洋可控源电磁场的影响特征,引入多相微增等效介质模型,结合非结构化网格实现了海洋可控源电磁三维矢量有限元正演.首先对多相微增模型进行了介绍,并分析不同参数影响下的等效电导率变化特征.而后,详细推导了频率域可控源电磁法非结构化矢量有限元方程.接着,利用预处理的IDR （s）迭代算法求解线性方程组.最后设计典型的海洋地电模型,通过正演计算,验证了本文算法的正确性和有效性,同时研究了海洋储油砂层中的金属矿的含量、孔隙度、含水饱和度等参数对于可控源电磁场响应特征的影响. 相似文献
13.
一种激发极化法2.5维正演的自适应有限元方法（英文） 总被引：2，自引：2，他引：0
传统的基于结构化网格有限元法采用的单元比较规则如矩形等,且网格剖分和加密要靠手动实现,所以传统的基于结构化网格有限元法不能准确和灵活地模拟复杂介质。本文采用易于模拟复杂介质模型的非结构化三角形网格进行剖分,且利用对偶加权后验误差估计指导网格自动细化过程,然后在电位模拟的基础上计算雅可比偏导矩阵,并依据Seigel（1959）理论实现激发极化法2.5维自适应有限元正演模拟算法。通过对垂直接触面模型进行正演分析,接收点附近网格得到了明显加密,电位数值解平均相对误差收敛到0.4%,视极化率平均相对误差收敛到1.2%,表明经自适应网格细化后,该算法数值解最终能收敛到精确解附近。最后对两个较复杂模型进行了正演计算与分析,进一步验证了该算法的准确性和灵活性。 相似文献
14.
海底采集到的电磁数据按照其主要包含的信息及研究目的大致可分为海洋可控源电磁场(CSEM)信号、天然场源大地电磁场(MT)信号、海洋环境电磁场信号以及其他随机干扰信号.常常通过计算功率谱密度、时频分析和极化分析的方法研究海洋电磁场特征.本文介绍一种新方法——时频方向谱分析法及其在实测海洋电磁数据处理中的应用,该方法能够在一定的时间-频率尺度上有效分辨场源信号的运动方向.对于海洋CSEM数据,利用该方法可以估算发射源的运动方向,进而在发射源或采集站方位信息缺失情况下,实现海洋CSEM数据的旋转电性轴处理.对于海洋电磁数据,利用该方法可以详细分析海水运动感应电磁场的信号特征. 相似文献
15.
Jianping Xiang N.B. Jones Daizhan Cheng F.S. Schlindwein 《Geophysical Prospecting》2002,50(6):565-576
The problem of discrimination between a valid induced polarization (IP) response and electromagnetic (EM) coupling effects is considered and an effective solution is provided. First, a finite dimensional approximation to the Cole‐Cole model is investigated. Using the least‐squares approach, the parameters of the approximate model are obtained. Next, based on the analysis of overvoltage, a finite dimensional structure of the IP model is produced. Using this overvoltage‐based structure, a specific finite dimensional approximation of the Cole‐Cole model is proposed. Summarizing the analysis of the finite dimensional IP model, it is concluded that the proposed IP model, which fits the field data much better than the traditional Cole‐Cole model, is essentially an RC‐circuit. From a circuit‐analysis point of view, it is well known that an electromagnetic effect can be described by an RL‐circuit. The simulation results on experimental data support this conception. According to this observation, a new method to discriminate between a valid IP response and EM coupling effects is proposed as follows: (i) use a special finite dimensional model for IP–EM systems; (ii) obtain the parameters for the model using a least‐squares approach; (iii) separate RC‐type terms and RL‐type terms – the first models the IP behaviour, the latter represents the EM part. Simulation on experimental data shows that the method is very simple and effective. 相似文献
16.
陆地可控源电磁法的观测资料可依据频段范围近似地划分为近区场、中间区场及远区场,但采用测量相互正交电、磁分量,并计算视电阻率的资料处理方式只适用于远区场数据.为更有效地利用陆地可控源电磁法不同区间场的观测资料,本文结合三维数值模拟技术并采用电场分量直接进行反演的策略,对不同区间电场的响应特征与探测效果进行了分析.数值模拟结果表明:近区电场的异常响应最明显,异常响应不随频率发生显著变化,但纵向分辨能力差;远区电场异常响应随频率发生显著变化,其探测深度取决于频率的高低;中间区场较为复杂,地表电场异常响应的等值线中心并不是位于异常体中心正上方,而是在沿场传播方向上向异常体与围岩的分界面处偏移,并且发现中间区场资料的加入会影响反演结果的准确性.综合合成数据和野外实测资料的反演结果,发现结合近区场和远区场资料而舍弃中间区场资料的反演效果更佳,这为陆地可控源电磁法资料的反演解释提供了一种有效途径. 相似文献
17.
本文以非结构化网格剖分为基础,推导了电导率呈任意各向异性条件下的海洋可控源电磁法二次场磁矢量位、标量位所满足的有限元方程.将不完全LU分解预条件因子(ILU)与Induced dimension reduction(IDR(s))迭代算法相结合对有限元方程进行求解.之后,采用加权移动最小二乘法对二次场矢量位、标量位进行求导得到相应的电磁场各分量.最后,开展了两个地电模型的电磁场计算与分析.结果表明:本文所提出的电导率呈任意各向异性条件下的海洋可控源电磁法三维正演算法正确;不完全LU分解预处理的IDR(s)算法计算效率明显优于常用的ILU-QMR、ILU-BICGSTAB算法;算法具有良好的通用性,可用于陆地电磁、航空电磁、井中电磁等地球物理勘探方法在电导率任意各向异性方面的研究. 相似文献
18.
为了提高海洋可控源电磁的地形响应模拟精度,探讨了三维交错网格剖分有限体积法求解频率域电磁场解的数值方法.在似稳场近似条件下,推导了有限差分法和有限体积法低频电磁场控制方程的离散表达式.通过对离散表达式的分析比较表明,有限体积法较之有限差分具有更高的计算精度,同时也具有与有限差分相当的计算效率.对含地形界面网格单元电导率采用加权平均处理,通过与二维有限元程序的计算结果对比,验证了该方法可以有效提高有限体积法对地形变化的模拟精度.应用有限体积法计算了一个三维带地形储层模型的可控源电磁响应,分析表明地形变化对电场分量影响明显,磁场分量对地形变化不敏感. 相似文献
19.
地层宏观电性各向异性会对可控源电磁响应产生重要影响.由于海底地层电性结构常表现为电导率各向异性,若仅对海洋可控源电磁(MCSEM)数据进行常规各向同性反演,有可能无法获得准确的反演解释结果,从而削弱MCSEM技术的可靠性.本文实现了电导率垂直各向异性(VTI)条件下频率域海洋可控源电磁数据三维反演算法.其中,三维正演采用基于二次场控制方程的交错网格有限体积法,并利用直接矩阵分解技术来求解离散所得的大型线性方程组,有利于快速计算多场源的响应.反演采用具有近似二次收敛性的高斯牛顿算法对目标函数进行最优化.最后,对具有VTI电性各向异性特征的盐丘构造模型的MCSEM合成数据分别进行了电导率各向同性和垂直各向异性三维反演,结果表明:各向同性三维反演算法无法对受VTI介质影响的MCSEM数据进行正确的反演解释,而垂直各向异性三维反演能够获得更为可靠的地下电阻率结构和异常体分布,展现出对海底电性各向异性结构更为优良的反演解释能力. 相似文献
20.
在使用电偶极发射源的可控源电磁法(CSEM)勘探中,发射源的方位、长度、形状等对观测数据有重要的影响,然而现有的大部分三维数值模拟方法没有全面地将这些因素考虑进来,很多都只能应对非常简单的场源形态,例如单一方位的点电偶极子,这有可能显著降低模拟结果的准确性.本文实现了基于交错网格有限体积(FV)离散的海洋CSEM三维正演算法,能够模拟形态相对复杂的场源,包括任意方位的有限长直导线和弯曲导线发射源.该算法使用一次场/二次场方法,只需对二次场使用FV法求解,避免了场源的奇异性问题;一次场的计算为一维正演问题,使用准解析法求解,并且只要在计算一次场时考虑复杂的场源形态便可以实现同样场源的三维正演.通过与一维理论模型的解析解对比验证了三维程序的准确性,并针对三维理论模型进行了一系列正演测试,初步考察了场源形态对三维正演结果的影响. 相似文献