共查询到19条相似文献，搜索用时 859 毫秒
1.
在地质情况复杂化以及地质勘探精细化的要求下,有必要开展张量CSAMT正反演研究.考虑到各向异性对正反演结果影响较大,本文开展基于主轴各向异性介质的张量CSAMT正演研究.根据含源的电场控制方程,引入主轴各向异性介质的张量电阻率,将电场总场分解为一次场和二次场,一次场由Key的算法得到,二次场利用有限差分法计算,即实现了主轴各向异性介质的张量CSAMT三维正演.通过与二维有限元算法结果比较来验证算法的正确性;通过设计几组主轴各向异性三维模型,分析卡尼亚电阻率和相位的响应特征.结果表明卡尼亚电阻率ρxy和ρyx分别对异常体X方向和Y方向的电阻率变化敏感,说明张量CSAMT方法对识别各向异性具有优势.本研究为开展任意各向异性介质的张量CSAMT正反演研究提供基础. 相似文献
2.
有限元法2.5维CSAMT数值模拟 总被引：9，自引：7，他引：2
目前国内大多数CSAMT资料处理都用一维反演方法,即假定地球介质是水平层状地电结构,这对于横向变化较大的介质结构的反演则是不真实的．为了合理的模拟可控信号源的三维特征应该使用三维有限元方法．受计算机内存和计算速度的限制,三维方法的实用化受到约束．在许多情况下,地电结构沿走向变化很小,只沿倾向发生变化．这种地电结构是二维的,而人工源是三维的,因此CSAMT资料的观测可用2．5维有限元方法进行数值模拟．本文从麦克斯韦方程组出发,建立了2．5维有限元CSAMT数值模拟方法,其核心是把地电参数变化小的走向方向转化成波数域,用一系列波数模拟三维源的特征．用一个横向均匀的三层地电结构模型展示了2．5维数值模拟的特征,并和一维模拟结果进行了比较,证实了有限元方法2．5维CSAMT数值模拟的可靠性．在此基础上对一个地电结构已知的实际模型进行了2．5维数值模拟并和该剖面的野外实测剖面数据进行了比较,进一步有效地说明了本文介绍的2．5维CSAMT数值模拟方法是仿真和可靠的,为在此基础上的2．5维反演打下了良好的基础． 相似文献
3.
利用三维（3D）全空间积分方程准线性解作为数值模拟手段,对包含电离层、大气层和地球介质层的“地-电离层”典型异常目标体多层介质模型进行了数值模拟,得到了偶极源长度100 km、电流200 A、收发距离远达1600 km的合理的异常电性目标体的电阻率-频率响应结果.通过对不同埋深异常目标体的电阻率-频率特征的讨论,认为当考虑电离层和大气层的“地-电离层”大尺度深层横向不均匀复杂介质模拟时,电磁场对深部目标体仍有很好的异常响应,但当异常体电阻率及同一电阻率的岩石埋深不同时,其电阻率-频率的响应特征有很大不同,这正反映了不同岩石物性参数的电阻率-频率响应特征.文中结果表明从长偶极大功率源（WEM）激励的电磁场观测资料区分岩石的电阻率参数是可能的,这为利用WEM观测资料建立电性参数和岩石/地层参数经验关系奠定了基础. 相似文献
4.
《地球物理学进展》2017,(5)
可控源音频大地电磁测深(CSAMT)在实际生产中具有广泛的应用.本文选取一组精度较高的波数,利用基于双二次插值的有限元算法对A型地电模型进行模拟计算,将计算结果与解析解对比,验证本文算法的正确性.再根据生产中常见的地球物理异常体特征,设置了不同的二维地电模型进行有限元模拟,分别得到TE及TM模式下的视电阻率及相位等值线断面图,通过对比分析可知模拟结果均能在不同程度上反映出异常体特征,进一步说明算法的正确性.后文基于沁水盆地北部煤层气富集与地下水、地质构造之间的关系,设计了不同地球物理模型并开展典型薄层模型正演模拟.模拟结果表明:CSAMT能有效获得地层的电性结构特征,为CSAMT在煤层气勘探可行性提供了理论基础,CSAMT可以作为煤层气非地震勘探的技术手段. 相似文献
5.
可控源电磁法(CSEM)的正演模拟与反演解释中常忽略激电介质的影响.本文基于电磁三维正反演开源平台ModEM,分别开展了陆地CSEM三维正反演以及考虑激电介质的陆地CSEM三维正演研究,并与一维模型解析解对比验证了三维计算的准确性.数值模拟结果显示,模型含激电效应与不含激电效应的电场E_x频率响应差异明显,激电效应会导致电场E_x幅值减小、相位增加,且对E_x的影响程度与激电异常体的频谱特性及观测的收发距均相关.合成数据的陆地CSEM反演结果显示,数据中含激电效应会使得反演出的异常体在深度、规模和电阻率值上都与真实异常体存在较大偏差,忽略观测数据中的激电效应会明显影响观测的准确性.本文研究表明当陆地CSEM三维观测资料中包含激电效应时,有必要考虑其影响. 相似文献
6.
在一些地层层理发育的地区,地下介质存在显著的电各向异性,此时基于各向同性模型解释含各向异性效应的可控源音频大地电磁（CSAMT）测深观测数据会导致错误的结果.本文通过引入3×3的对称正定张量表征电导率各向异性,采用非结构四面体网格和矢量有限元方法离散电场满足的矢量Helmholtz方程,并将电磁场源等效为系列电偶极子,实现任意各向异性介质中CSAMT高效数值模拟.本文首先通过层状各向异性模型检验三维有限元算法的精度和有效性,进一步建立三维地电模型研究异常体各向异性和围岩各向异性对CSAMT响应的影响,最后使用视电阻率极性图来识别各向异性电导率主轴方向.数值模拟结果表明,各向异性电导率对CSAMT视电阻率幅值及分布规律都有很大影响,视电阻率极性图能够很好地识别各向异性主轴方向. 相似文献
7.
8.
地电体对频率域电磁波激发源的响应为电磁感应和激电效应的综合响应.传统CSAMT法进行数据正反演时认为大地介质电阻率是与频率无关的实数,而实际上因为激电效应,地下可极化体的电阻率是一个与频率相关的复数.为推进二者总体响应研究,并扩展激电法的应用范围,同时提高电磁法勘探的精度,本文基于Dias模型,以复电阻率代替不考虑地电体极化效应的直流电阻率,对CSAMT场源一维层状模型进行了正演模拟,为提取CSAMT信号中所含激电信息提供理论基础.结果表明,考虑激电参数后,视电阻率及相位响应曲线出现明显异常（包括远场、过渡场、近场）;极化前后振幅比值异常峰值、相位差值异常峰值可直观体现激电异常;异常峰值与极化层层厚、埋深以及电阻率变化有连续的对应关系.认为从频率域电磁法信号中提取激电信息有乐观的前景. 相似文献
9.
地球物理学中关于电磁波勘探研究通常采用的是地球半空间模型.然而,对于几十公里的有限长电缆源（长偶极源）,远距离电磁波场探测必须要考虑电离层的影响,它是一个全空间问题.关于包含电离层、空气层和地球介质（我们称“地-电离层”模式）的电磁波场特征的研究在国外较少,国内几乎是空白.本文采用全空间积分方程法首先对小尺度的可控源电磁波场特征进行了研究,由于此时电离层的影响可忽略,它应该和半空间成熟的CSAMT模拟结果一致,对比结果表明,二者是一致的,验证了全空间模拟方法的可靠性和有效性.随后进行了50 km长电缆电离层和空气层高度都为100 km的“地-电离层”模式大尺度电磁波场模拟,以探讨大尺度可控源电磁波场的特征.给定频率的“地-电离层”模式电磁场的衰减曲线表明长电缆远距离电磁波场由于受电离层的作用存在衰减逐渐变小的过渡场和衰减变小的波导场.为了探讨复杂介质“地-电离层”模式电磁波特征,对“地-电离层”模式的典型地盾和地台多层介质模型进行了数值模拟,得到了偶极源长度50 km、电流200A、收发距离远达1600 km和2500 km的合理的电磁场结果.最后,对一简单含油储层结构模型进行了长偶极、大功率、远距离电磁波场响应计算.储层横向不均匀复杂结构模拟的结果表明,考虑电离层和大气层的“地-电离层”模式大尺度深层复杂介质模拟时,电磁场对深部目标体仍有很好的异常响应. 相似文献
10.
《地球物理学进展》2016,(2)
电性源瞬变电磁具有快捷、探测深度大、对异常体探测效果好的特点,然而对其观测分量、装置、电流波形的选择还没有进行过详细的阐述,为此本文做了详细的研究.通过对电性源赤道装置与轴向装置的对比研究,表明赤道装置可以对地下结构取得较好的探测效果,在探测相同深度时,轴向装置收发距较小,信噪比较高.通过对电流波形的研究表明,阶跃激励源激发的场视电阻率从早期到晚期很好的反映了地质结构从浅到深的变化,可以很好的勾勒地电结构的变化,而脉冲波激发的场对目标体反映明显,然而其视电阻率在过渡场阶段存在畸变现象.通过对Hz场和Ex场的研究表明,Hz场在小收发距时即可进行大深度探测,在小收发距下信号强度较高,其应具有较高的信噪比,Ex场需要在大收发距观测,对高阻体探测效果比Hz好. 相似文献
11.
与可控源音频大地电磁（CSAMT）相比,广域电磁法通过采用全区视电阻率定义,突破了卡尼亚视电阻率所需的远区条件限制,极大拓展了可控源电磁观测区域和探测深度.考虑到电偶源激发场的三维特征以及地下复杂三维结构,为提高广域电磁数据解释精度,本文实现了基于二次耦合势的广域电磁法三维正演计算.该算法利用Helmholtz定理将麦克斯韦方程转化为库伦规范下的磁矢势和电标势耦合方程,有效改善了离散所得大型线性方程组的谱性质,并通过强加散度条件来消除电场伪解的影响.此外,采用散射场方法,其中一次场使用准解析法求解,二次场使用有限体积法求解,克服了局部激发场源奇异性问题.通过与一维层状模型下电偶源产生的电磁场准解析解对比,验证了本文算法的正确性.在此基础上,利用本文的正演算法对比分析了广域电磁法与CSAMT对典型三维目标体的探测能力,结果表明在相同的观测条件下,广域电磁法能够更准确地反映地下目标体信息,拥有更优的分辨能力. 相似文献
12.
研讨了频率域电磁法中不同源装置的大地电磁测深、线源频率电磁测深和偶极源频率电磁测深阻抗视电阻率的源效应影响特征。在唯象分析的基础上，提出了几种电磁测深法阻抗视电阻率的相互换算法──源效应校正法（大地电磁测深二维ＴＥ极化视电阻率和其它两种电磁法的赤道装置二维阻抗视电阻车）。模型试验表明，利用这一源效应校正法可以由大地电磁二维视电阻率近似地计算出线源频率电磁二维阻抗视电阻率。这一方法被尝试应用于由线源频率电磁二维阻抗视电阻率估算偶极源频率电磁二维阻抗视电阻率。 相似文献
13.
THE DISPERSION RELATION, OF EARTH''S ELECTROMAGNETIC RESPONSE FUNCTIONS AND THE JOINT INVERSION OF DATA 
下载免费PDF全文
下载免费PDF全文 On the basis of the dispersion relations of MT field, the necessity and applied prospects of the joint inversions using a pair of MT response functions which are correlative with the dispersion relations, are infered. A filter coefficient algorithm is made, with which the corresponding impedance phase data can be estimated using a set of apparent resistivities. The tests for the observed MT data show that when comparing the impedance phase estimated using the dispersion relation with the ob served phase, it can be checked whether the dispersion relation between observed apparent resistivity and phase data is satisfied or not, and that the use of the phase data corrected using the dispersion relation in the joint inversion is advantageous to obtain more confident results. It is shown that joint inversions are more advantageous than single parameter inversions, and that in the most case the joint inversion using the apparent resistivities of impedance real and imaginary parts is more advantageous than the jointinversion using the normal apparent resistivity and impedance phase. The existence of the dipersion relations between the ratio apparent resistivity and corresponding impedance phase of the orthogonal electric and magnetic field horizontal Components in the frequency EM sounding with horizontal electric dipole(FEMS) are discussed, the better effect of the joint inversion using the pair of EM response functions is obtained. The problems on the one-dimensional joint inversion for the MT and FEMS apparent resistivities, for which the observed frequency bands partly overlape each other, are studied. It is shown that this joint inversion is applicable and effective:the joint inversions of the practical data for two kinds of EM methods at two sites give the results well corresponding to the drilling data. The simulated MT inversions for the data of two kinds of EM methods are made, and more confident results also are obtained. 相似文献
14.
本文提出了能提高异常体分辨能力,同时得到绝对电导率的地面磁电阻率数据三维反演方法.磁电阻率响应用准直流的低频磁场代替;数值模拟由频率域电场满足的Helmholtz方程出发,采用三维交错网格有限差分法;长直导线源作为发射源,其中源的计算包含在背景场中;结合地面磁电阻率数据各分量的特点,选择y分量进行反演研究;反演采用三维非线性共轭梯度反演技术,为了提高异常体的深度分辨能力,进行迭代重构反演;用印模法对初始模型进行重构,采用的是辅模型在浅部,元模型在深部的组合方式.从合成数据和实际数据的反演结果可以得到以下的认识:(1)由频率域麦克斯韦方程组出发,低频磁场数据反演可以直接得到电导率,而不是相对电导率之比;(2)采用印模法组合初始模型,进行迭代重构反演,可以提高地面磁电阻率数据反演对异常体的分辨能力,确定埋深位置,同时不会丧失对于浅部异常体的分辨能力;(3)在结合印模法的地面磁电阻率数据三维反演中,深部异常体的分辨能力受地表不均匀导电体影响较小;(4)确定印模深度可以采用上一次重构反演结束时的模型变化量,通过相邻两次重构反演结束时的模型变化量之差来确定迭代重构是否终止.因为静磁场与重力场在数学上的相似性,本文的反演方法可以被运用到重力场等位场的地面数据的反演中. 相似文献
15.
为了易于模拟野外复杂地形和地下任意形状地电体模型,将有限元单元网格设计为三角单元;并考虑到野外实际勘探中,地球介质的电性参数均是连续变化的情况,单元内的场值和电性参数被设计为双线性变化;推导出二维起伏地形条件下大地电磁法有限元数值模拟算法;根据单元节点主场值和线性插值形函数间的关系,计算出单元节点的辅助场值;在二维起伏地形情况下,定义TE、TM模式视电阻率和阻抗相位.4个模型的计算的结果与解析法的均方根误差小于1%,地形模拟与前人的计算结果相符,模拟倾斜界面异常体,能有效的反映出其异常形态. 相似文献
16.
In order to interpret field data from small-loop electromagnetic (EM) instruments with fixed source–receiver separation, 1D inversion method is commonly used due to its efficiency with regard to computation costs. This application of 1D inversion is based on the assumption that small-offset broadband EM signals are insensitive to lateral resistivity variation. However, this assumption can be false when isolated conductive bodies such as man-made objects are embedded in the earth. Thus, we need to clarify the applicability of the 1D inversion method for small-loop EM data. In order to systematically analyze this conventional inversion approach, we developed a 2D EM inversion algorithm and verified this algorithm with a synthetic EM data set. 1D and 2D inversions were applied to synthetic and field EM data sets. The comparison of these inversion results shows that the resistivity distribution of the subsurface constructed by the 1D inversion approach can be distorted when the earth contains man-made objects, because they induce drastic variation of the resistivity distribution. By analyzing the integrated sensitivity of the small-loop EM method, we found that this pitfall of 1D inversion may be caused by the considerable sensitivity of the small-loop EM responses to lateral resistivity variation. However, the application of our 2D inversion algorithm to synthetic and field EM data sets demonstrate that the pitfall of 1D inversion due to man-made objects can be successfully alleviated. Thus, 2D EM inversion is strongly recommended for detecting conductive isolated bodies, such as man-made objects, whereas this approach may not always be essential for interpreting the EM field data. 相似文献
17.
从电偶源三维地电断面可控源电磁法的二次电场边值问题及其变分问题出发,采用任意六面体单元对研究区域进行剖分,并且在单元分析中同时对电导率及二次电场进行三线性插值,实现电导率分块连续变化情况下,基于二次场的可控源电磁三维有限元数值模拟.这个新的可控源电磁三维正演方法可以模拟实际勘探中地下任意形状及电性参数连续变化的复杂模型.理论模型的计算结果表明,均匀大地计算的视电阻率误差和相位误差分别为0.002%和0.0005°.分层连续变化模型的有限元计算结果表明,其与对应的分层均匀模型解析结果有明显差异.三维异常体组合模型以及倾斜异常体等复杂模型的有限元计算结果也有效地反映了异常形态. 相似文献
18.
局部畸变问题曾经困扰大地电磁资料反演解释几十年,大地电磁三维数值模拟技术的发展为剖析局部畸变特点和得到可靠的反演成像结果提供了技术基础。本文采用三维数值模拟成像方法对典型三维局部畸变模型进行模拟分析。三维数值模拟结果显示:电场分量垂直电性分界面的极化模式视电阻率曲线(对应二维情况下TM模式)在穿越低阻异常体界面时,曲线会先上移后下移,而在穿越高阻异常体界面时,曲线会先下移后上移,这与电性分界面处积累面电荷产生的二次电场有关。三维模型中XY模式、YX模式视电阻率和相位在三维异常体附近的水平变化是呈现近似垂向对称的,该现象与电场垂直跨越电性界面时视电阻率的变化规律是吻合的,当测线分别沿X方向和Y方向展布时,三维情况下的XY和YX模式分别对应二维情况下的TM模式。低阻小异常体对区域构造响应的畸变影响比高阻小异常体要严重。低阻小异常体对二维区域响应的两种极化模式视电阻率和相位都有非常明显的畸变影响,相比较而言对TE模式的畸变要大于TM模式,因此我们在做二维反演解释时,可优先考虑拟合TM模式数据。位于小异常体中心上方测点的三维畸变响应虽然与对应真实二维区域响应的差异比较大,但可以等效于某种二维模型响应,这种由局部畸变造成的假二维响应在实际野外数据的解释中是需要注意的。 相似文献
19.
电性源地空瞬变电磁法具有工作效率高、勘探深度大、采集信号信噪比高、适用于地形地质条件复杂地区等优点.但是,到目前为止,由于尚未建立起该方法的解释系统,大大制约了该方法的发展.本文旨在建立起完整的地空电磁探测系统,丰富整个探测系统的理论.本文围绕地空瞬变电磁法全域视电阻率定义、瞬变电磁虚拟波场的克希霍夫偏移成像、逆合成孔径成像方法三个科学问题进行了系统研究.提出了用磁场强度定义全域视电阻率的迭代算法,理论模型试验结果表明计算出的视电阻率曲线首支趋于第一层电阻率,尾支趋于最后一层电阻率,实现了全空域、全时域视电阻率的计算;在先前研究的基础上,实现了适合电性源地空装置的瞬变电磁虚拟波场的克希霍夫偏移成像;采用相关迭加技术,实现孔径内多测点数据合成,将传统的单点处理方式发展成为逐点推移多次覆盖的逆合成孔径处理方法.层状模型试验表明:(1)全域视电阻率能够光滑、完整、渐变地反映出模型的电性信息变化;(2)当改变三层模型中间层电阻率时,全域视电阻率曲线随着参数的改变分异明显,对电性层的识别容易且直观;(3)由于在电阻率计算中同时考虑了接收机高度、偏移距、时间等各参数的影响,全域视电阻率可实现全空域、全时域的视电阻率计算.含水采空区的复杂模型算例表明:(1)根据不同测线的全域视电阻率结果可以看出,在靠近采空区的位置,全域视电阻率断面可以清晰地反映出采空区的空间位置,随着测线离采空区越来越远,采空区异常越来越弱直至消失;(2)波场变换和偏移成像的结果显示存在两个电性差异较大的界面,上界面指示地表,由于空气和大地之间的电性差异较大,故该界面波场信号反映强烈,遍布整个区域,下界面异常信号则主要集中在中部,向外逐渐减弱,指示采空区;(3)逆合成孔径成像结果表明地表界面在合成前后没有变化,而采空区异常合成后范围明显变小,且异常边界清晰,指示的采空区位置与模型吻合很好.本文借助于逆合成孔径雷达成像的基本思想,建立了一套电性源瞬变电磁地空逆合成孔径成像方法.基于反函数思想结合迭代算法提出的电性源地空瞬变电磁法的全域视电阻率定义方法,实现了全空域、全时域的视电阻率计算;借鉴瞬变电磁拟地震偏移成像算法,实现了瞬变场的三维成像;借鉴逆合成孔径雷达的思想,提出电性源瞬变电磁地空逆合成孔径算法,进一步提高了成像的分辨率.采空区模型算例表明相关叠加合成确实具有增强有用信号、提高信噪比、提高分辨率的诸多优点,证实了瞬变电磁地空逆合成孔径成像方法的有效性. 相似文献