共查询到19条相似文献，搜索用时 448 毫秒
1.
油气藏自身会产生激电效应,明确海洋可控源电磁场激电效应特征非常重要。本文基于自适应有限元方法实现了带激电效应的2.5D海洋可控源频率域电磁场正演。采用非结构化三角网格单元剖分地电模型,可以模拟任意起伏地形和复杂地电构造。利用Cole-Cole模型表述油气藏的激电效应,实现了考虑激电效应的海洋电磁复电阻率正演算法。将本文的计算结果与已发表的一维模型结果对比,检验了本文算法的正确性。最后,设计了不同的极化模型,针对不同的场源方位角、激电参数和地形,计算了的电磁场响应结果,分析了激电效应对海洋可控源电磁场的影响规律,结果表明油气藏的激电效应和海底地形对海洋可控源电磁场有一定的影响。 相似文献
2.
全局弱式无单元法是在有限单元法基础上发展起来的一种数值模拟方法,它采用局部支持域内的节点信息来构造形函数实现局部精确逼近,摆脱了单元,仅依赖于节点信息,具有预处理简单、模拟精度高、灵活性强的特点,适用于复杂地电条件下直流电阻率正演模拟.本文采用RPIM构造直流电阻率全局弱式无单元法形函数,利用RPIM形函数推导了直流电阻率全局弱式无单元法方程.然后,编制了直流电阻率全局弱式无单元法正演模拟Fortran程序,利用该程序对典型的地电模型进行了正演模拟,并将正演结果与基于线性插值的FEM正演结果及解析解进行对比,结果表明采用RPIM形函数的全局弱式无单元法用于直流电阻率正演模拟的正确性及有效性,且在同等条件下,全局弱式无单元法模拟精度高于矩形剖分的FEM,更有利于指导电法勘探的数据解译;利用该程序对复杂地电模型进行了正演模拟,结果表明全局弱式无单元法对复杂地电模型模拟效果良好,适应性强,灵活性高,可任意加密节点提高模拟精度. 相似文献
3.
4.
CSAMT在山区金属矿勘查中,采用各种滤波和相位积分之类的处理方法,校正因地形起伏和局部电性不均匀引起的静态效应,往往难保奏效,开发消除静态效应的新方法是提高CSAMT资料处理与解释水平和方法应用效果的重要研究课题. 本文以如何消除地形影响为重点,对起伏地形下CSAMT二维大地三维源地电模型,采用加权余弦数值积分法,进行波数域电磁场二维有限单元法正演. 为模拟复杂地形地电模型,选取交叉对称网格三角形剖分法,实现了在国内常用赤道电偶极装置的CSAMT二维正演计算;在二维正演的基础上,开发了基于奥克姆反演法的CSAMT二维反演技术,研制出一套起伏地形下CSAMT二维正反演处理与解释方法技术系统. 通过理论模型试算和实测数据处理证实,本系统能有效地削减起伏地形影响. 在找矿应用中,该系统反演的电阻率断面,极大地消除了起伏地形影响和静态效应,突显出清晰的控矿构造和矿体的异常,取得了重要成效. 相似文献
5.
基于非结构网格有限元方法开展了三维复杂地电模型的线源井地电法的高效正演模拟研究,探讨了通过求取电场响应导数来刻画目标体边界范围、采用差异场地形校正技术来消除地形影响等措施对井地电法成像的效果和精度的影响。并通过对比与解析解,验证了本文数值解算法的有效性。模型计算结果表明：积水巷道的空间位置和走向均会引起视电阻率的显著变化,视电阻率变化率的极值准确且清晰地指示了巷道边界的位置;电位的归一化总水平导数极大地提高了井地电法对目标体复杂边界位置的识别能力;地形对井地电场分布的影响也很大,其视电阻率响应与地形形状近似呈对称关系,利用差异场技术能有效地削弱地形对井地电法高精度成像的影响。 相似文献
6.
7.
作为一种基于节点的计算方法,无网格法具有构造高阶导数方便,自适应分析便利的优点,特别适合复杂地质构造的数值模拟.本文针对实际地球物理勘探中存在的起伏地形和各向异性的地电结构,提出用无网格法来模拟大地电磁响应,采用复合二次径向基函数构造了形函数,推导了大地电磁无网格法等价线性方程组,研究了系数矩阵的压缩存储方法以及大型稀疏复线性方程组快速求解算法,实现了起伏地形下各向异性的2D大地电磁高精度数值模拟.基于层状模型验证了算法的正确性,计算结果表明:无论是TM模式还是TE模式,计算相对误差均小于1%;通过对地垒和地堑模型的模拟,得出了起伏地形对视电阻率和相位的影响规律;对起伏地形下含有不同各向异性系数异常体的模型进行了数值模拟,为开展复杂地质结构的电磁场特征研究以及地形校正奠定了理论基础.随着计算科学的快速发展,无网格法必将成为新的高精度电磁场数值模拟方法. 相似文献
8.
利用CSAMT电磁场响应提取激电效应的方法初探 总被引：1，自引：0，他引：1
在传统的CSAMT电磁场一雏正演中加入柯尔-柯尔模型,模拟可极化地电断面的CSAMT电磁场响应.通过对均匀半空间、H型地电断面的理论计算,分析了柯尔-柯尔模型参数对CSAMT的Ex、Hy影响特征.基于电场相对磁场对激电效应的反映更加灵敏,提出分别利用Hy与Ex分量反演电阻率、厚度以及激电参数的方法.反演计算表明,该方法克服了电阻率、厚度、激电参数一同反演时受初始模型影响太大,容易陷入局部极小值的问题,可以在一定程度上提取出合适的激电参数. 相似文献
9.
地电体对频率域电磁波激发源的响应为电磁感应和激电效应的综合响应.传统CSAMT法进行数据正反演时认为大地介质电阻率是与频率无关的实数,而实际上因为激电效应,地下可极化体的电阻率是一个与频率相关的复数.为推进二者总体响应研究,并扩展激电法的应用范围,同时提高电磁法勘探的精度,本文基于Dias模型,以复电阻率代替不考虑地电体极化效应的直流电阻率,对CSAMT场源一维层状模型进行了正演模拟,为提取CSAMT信号中所含激电信息提供理论基础.结果表明,考虑激电参数后,视电阻率及相位响应曲线出现明显异常（包括远场、过渡场、近场）;极化前后振幅比值异常峰值、相位差值异常峰值可直观体现激电异常;异常峰值与极化层层厚、埋深以及电阻率变化有连续的对应关系.认为从频率域电磁法信号中提取激电信息有乐观的前景. 相似文献
10.
11.
本文实现了2.5维电导率正交各向异性海洋可控源电磁等参有限元数值模拟.利用傅里叶变换导出了电导率正交各向异性2.5维海洋可控源电磁法波数域电磁场耦合方程,采用伽里金加权余量法推导了相应的有限元方程;采用任意四边形单元对研究区域进行剖分,在单元中进行双二次插值,将有限元方程化为线性代数方程组;最后,求解线性方程组并进行反傅里叶变换获得空间域电磁场值.这个方法可以模拟海底起伏地形条件下地下任意形状电导率正交各向异性的复杂模型.与一维模型的数值模拟结果对比表明,电磁场数值解与解析解吻合.二维模型的计算结果与二维自适应非结构有限元模拟结果也吻合.水平海底二维地电模型考察了不同各向异性系数对海洋可控源电磁响应的影响特征.海底起伏地形地电模型的数值结果表明,电导率各向异性对海洋可控源电磁响应影响明显,有可能淹没海底地形和高阻油气藏引起的异常. 相似文献
12.
Complex terrain causes great MT noise.This paper puts forward a FEM model method using adaptive topography and quadratic elements based on studies by previous researchers.This method can model all kinds of complicated terrain and geoelectric bodies preferably.The numeric modeling,calculation of the auxiliary field and definition of resistivity are deduced by electromagnetic equations.Lastly,several examples are presented,which show the method is rapid,effective and has highly accurate. 相似文献
13.
We conducted a study on the numerical calculation and response analysis of a transient electromagnetic field generated by a ground source in geological media. One solution method, the traditional discrete image method, involves complex operation, and its digital filtering algorithm requires a large number of calculations. To solve these problems, we proposed an improved discrete image method, where the following are realized: the real number of the electromagnetic field solution based on the Gaver–Stehfest algorithm for approximate inversion, the exponential approximation of the objective kernel function using the Prony method, the transient electromagnetic field according to discrete image theory, and closed-form solution of the approximate coefficients. To verify the method, we tentatively calculated the transient electromagnetic field in a homogeneous model and compared it with the results obtained from the Hankel transform digital filtering method. The results show that the method has considerable accuracy and good applicability. We then used this method to calculate the transient electromagnetic field generated by a ground magnetic dipole source in a typical geoelectric model and analyzed the horizontal component response of the induced magnetic field obtained from the “ground excitation–stratum measurement” method. We reached the conclusion that the horizontal component response of a transient field is related to the geoelectric structure, observation time, spatial location, and others. The horizontal component response of the induced magnetic field reflects the eddy current field distribution and its vertical gradient variation. During the detection of abnormal objects, positions with a zero or comparatively large offset were selected for the drillhole measurements or a comparatively long observation delay was adopted to reduce the influence of the ambient field on the survey results. The discrete image method and forward calculation results in this paper can be used as references for relevant research. 相似文献
14.
可控源电磁法具有分辨率高及抗干扰能力强等特点,是一种重要的地电磁勘探方法.目前,可控源电磁法的高精度正演计算一直是其核心研究问题之一.传统积分方程法一般采用近似积分公式、简单矩形网格和近似的奇异性体积分计算技术,制约了体积分方程法处理复杂地下异常体的能力,降低了计算精度.针对上述问题,本文基于完全积分公式、四面体非结构化网格和奇异体积分的精确解析解来高精度求解复杂可控源电磁模型的正演响应.首先,从电场积分公式出发,推导了可控源电磁问题满足的积分方程;其次,借助于非结构化四面体网格离散技术,实现了地下复杂异常体的有效模拟.最后,利用散度定理把强奇异值体积分转换为一系列弱奇异性的面积分公式,并通过推导获得了这些弱奇异性的面积分公式的解析解,从而最终实现三维可控源电磁问题的高精度积分求解.以块状低阻体地电模型为测试模型,采用本文提出的积分方程方法获得的数值解与其他公开数值算法解进行对比分析,其对比结果具有高度的吻合性,验证了算法的正确性;同时,设计了球状及复杂地电模型进行算法收敛性测试,进一步验证算法的正确性以及能够处理地下复杂模型的能力. 相似文献
15.
从电偶源三维地电断面可控源电磁法的二次电场边值问题及其变分问题出发,采用任意六面体单元对研究区域进行剖分,并且在单元分析中同时对电导率及二次电场进行三线性插值,实现电导率分块连续变化情况下,基于二次场的可控源电磁三维有限元数值模拟.这个新的可控源电磁三维正演方法可以模拟实际勘探中地下任意形状及电性参数连续变化的复杂模型.理论模型的计算结果表明,均匀大地计算的视电阻率误差和相位误差分别为0.002%和0.0005°.分层连续变化模型的有限元计算结果表明,其与对应的分层均匀模型解析结果有明显差异.三维异常体组合模型以及倾斜异常体等复杂模型的有限元计算结果也有效地反映了异常形态. 相似文献
16.
Three-dimensional numerical modeling of full-space transient electromagnetic responses of water in goaf 总被引：2，自引：0，他引：2
The full-space transient electromagnetic response of water-filled goaves in coal mines were numerically modeled. Traditional numerical modeling methods cannot be used to simulate the underground full-space transient electromagnetic field. We used multiple transmitting loops instead of the traditional single transmitting loop to load the transmitting loop into Cartesian grids. We improved the method for calculating the z-component of the magnetic field based on the characteristics of full space. Then, we established the fullspace 3D geoelectrical model using geological data for coalmines. In addition, the transient electromagnetic responses of water-filled goaves of variable shape at different locations were simulated by using the finite-difference time-domain (FDTD) method. Moreover, we evaluated the apparent resistivity results. The numerical modeling results suggested that the resistivity differences between the coal seam and its roof and floor greatly affect the distribution of apparent resistivity, resulting in nearly circular contours with the roadway head at the center. The actual distribution of apparent resistivity for different geoelectrical models of water in goaves was consistent with the models. However, when the goaf water was located in one side, a false low-resistivity anomaly would appear on the other side owing to the full-space effect but the response was much weaker. Finally, the modeling results were subsequently confirmed by drilling, suggesting that the proposed method was effective. 相似文献
17.
We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets. 相似文献
18.
The induced polarization (IP) in rocks and minerals is of significance to the marine controlled-source electromagnetic (CSEM) field. We propose an adaptive finite-element algorithm for the 2.5D frequency-domain forward modeling of marine CSEM that considers the induced polarization. The geoelectrical model is discretized using an unstructured triangular elemental grid that accommodates the complex topography and geoelectrical structures. We use the Cole–Cole model to describe the IP and develop a complex resistivity forward modeling algorithm. We compare the simulation results with published 1D model results and subsequently calculate the electromagnetic field for variable azimuth sources, IP parameters, and topography. Finally, we analyze the IP effect on the marine CSEM field and show that IP of oil reservoirs and topography affects the marine CSEM electromagnetic field. 相似文献
19.
在吸收并修正部分前人工作的基础上，给出了电导率分块均匀的瞬变电磁25维有限单元法正演算法. 在该算法中，（1）采用在矩形网格中进一步三角细化的剖分方式，同时考虑到矩形单元中心节点的性质，求解线性方程组前用高斯消元法消除矩形网格中心节点上的待求未知量，这样一方面较准确地模拟出任意二维复杂地电断面，另一方面又大大地节省了计算量；（2） 采用计算二次场的算法，事先仅需获得电场和磁场沿地电构造走向方向两个分量的一次场，故而不但计算精度较高，而且不明显增加计算量. 最后，对几例典型模型进行正演，并与解析解对比，验证了本文方法的有效性. 相似文献