首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 125 毫秒
1.
海底表面磁源瞬变响应建模及海水影响分析   总被引:11,自引:4,他引:7       下载免费PDF全文
刘长胜  林君 《地球物理学报》2006,49(6):1891-1898
根据电磁场理论,推导了磁偶源和接收点均位于海水中时层状海底模型的频域电磁场响应一般表达式,并通过此式,得到了海水为均匀半空间和有限海水深度两种情况下,垂直磁偶极装置、中心回线和重叠回线分别置于均匀半空间海底表面时的瞬变电磁响应(磁场和感应电压)表达式. 这些表达式将瞬变响应和海底的电导率等参数有机联系在一起,为海底瞬变电磁法的正演计算和反演解释提供了理论基础. 仿真计算表明,海水的存在不仅使得瞬变响应曲线形态发生变化,而且影响其对海底电导率的分辨能力.  相似文献   

2.
矩形大定源层状模型瞬变电磁响应计算   总被引:4,自引:2,他引:2       下载免费PDF全文
针对矩形大定源介绍了一种层状导电模型瞬变电磁响应的快速计算方法.其基本思想是,将大定源回线源用较少的有限大小方形回线叠加.在频率域,该叠加过程导出矩形大定源回线模型响应可以表示为一个标准的Hankel积分,叠加效应仅仅是相应Bessel函数的空间积分过程.瞬变电磁响应可以利用余弦变换由频率域响应获得.均匀半空间模型计算结果和解析解较为吻合,层状模型模拟结果符合物理规律,不同测点位置感应电动势的视电阻率具有很好的一致性.新方法误差只是来源于圆回线对小方形回线的等效误差.与均匀半空间解析结果比较,只用16个方形回线等效,相对误差就可以小于10-3.  相似文献   

3.
瞬变电磁法中最常用矩形或圆形回线,但有时由于地形限制,只能使用不规则回线发射和接收.如果仍用现有的常规回线理论进行处理解释,得到的结果会存在很大的偏差.本文首先进行了不规则回线瞬变电磁法一维正演理论研究,基于电偶极子源的频率域响应公式,通过沿回线积分和时频转换,推导出不规则回线源在水平层状介质中的时间域响应公式.采用欧拉算法、高斯积分和快速汉克尔变换,计算不规则回线内任意一点处的磁场响应.利用改进的二分搜索法计算全区视电阻率,并在此基础上利用烟圈法反演电阻率和深度.通过四个典型地电模型的正反演计算,表明烟圈法反演能够有效反映地电模型的大致形态,可以用于不规则回线瞬变电磁数据的快速反演解释.  相似文献   

4.
针对我国西南印度洋合同区热液硫化物矿快速、有效以及便捷地质勘探装备的要求,研发了由甲板控制系统、万米光电复合缆、仪器舱拖体和天线拖体组成的深海6000 m拖曳式瞬变电磁系统.为了便于拖体布放和快速发现异常,选择重叠回线收发装置类型;采用理论数值计算确定了在拖曳高度不大于50 m的前提下,观测1~100ms窗范围内的二次场响应,可以发现近海底深海热液硫化物矿堆;另外,通过不同拖延深度海上试验,研究了拖曳深度对瞬变电磁的影响规律及仪器性能.大洋第30航次第二航段在西南印度洋脊热液区,应用该系统发现了明显的瞬变电磁异常.印证了系统的有效性.  相似文献   

5.
层状介质中地下瞬变电磁场全空间效应   总被引:2,自引:1,他引:1  
全空间效应是矿井瞬变电磁方法固有的特殊理论问题之一.建立了3种典型的层状地质模型,采用时域有限差分法求解Maxwell方程,并引入电磁波衰减因子EA(dB),来分析全空间效应的影响特征及影响因素.研究结果表明,地下瞬变电磁法工作装置中,分离回线装置对低阻目标层的探测效果优于重叠回线装置;相对于半空间低阻层瞬变响应,全空间低阻层瞬变响应的时间范围更宽,因而需要更长的观测时间;全空间效应影响因素为顶、底板内低阻层相对电阻率差异大小和目标层深度与厚度,当目标层电阻率大于顶板内低阻层电阻率时,所观测响应主要为低阻层响应,此时将会增大资料解释的难度.此外,研究还显示,瞬变电磁法对低阻层的纵向分辨率与低阻层厚度相关,厚度越小分辨率越高,当厚度较大时,其纵向分辨能力较差.研究成果为实际工作参数选择及资料处理与解释提供了参考.  相似文献   

6.
大回线源瞬变电磁响应理论研究回顾及展望   总被引:4,自引:3,他引:1  
解析公式的求解是电磁探测方法基础理论研究的重要内容.以往用于大回线瞬变电磁法的理论公式,有圆形重叠回线表达式、中心回线表达式.为了更接近勘探工程中使用的矩形回线源,和为提高效率在回线中心1/3区域观测的实际情况,采用了大定源回线的解析公式.由此定义的视电阻率,克服了偏离中心点观测引起的边缘效应,提高了解释精度,也将大定...  相似文献   

7.
多孔径瞬变电磁场物理模拟   总被引:1,自引:0,他引:1       下载免费PDF全文
为了提高瞬变电磁法对地探测精度,提出了一种新颖的瞬变电磁装置一多孔径发射装置,并进行了多孔径瞬变电磁场物理模拟观测.分别使用四孔径发射装置和常规回线发射装置,进行不同收一发距情况下的一次场观测,并对铜板在不同埋藏深度情况下的二次响应进行观测,通过比较两种不同装置下的响应曲线和计算观测数据误差,最终发现使用四孔径发射装置可使一次场强度提高34%,可使二次场强度提高31%.结论表明,多孔径发射源可以改变电磁场的空间分布,认为瞬变电磁场存在类相干现象.研究成果对提高瞬变电磁的探测精度具有很重要的意义.  相似文献   

8.
瞬变电磁接收装置对浅层探测的畸变分析与数值剔除   总被引:13,自引:11,他引:2       下载免费PDF全文
从理论上对瞬变电磁系统的接收装置频率特性进行分析,得出了瞬变信号在早期发生畸变的根本原因,揭示出接收装置在发射电流关断期间和电流关断后对一次场和二次场的影响关系,给出基于导电球体模型下接收装置的谐振频率与最小取样时间的关系图,进而讨论了接收装置对浅层目标体探测的影响.当瞬变电磁测量系统在电流关断开始时刻记录全程瞬变响应和发射电流波形时,如果接收装置的位置和谐振频率已知,就可以通过数值计算方法从根本上剔除接收装置对早期瞬变信号的影响,即使接收装置存在过渡过程,也同样可以实现近地表浅层的探测,从而缩短瞬变电磁法的浅层勘探盲区,提高近地表的探测分辨率和精度.  相似文献   

9.
对同一回线和"8"字形回线进行实验对比研究.采用控制变量法对两种装置进行了多方面实验对比,并在同一实际工作区域运用同一回线和"8"字形回线进行实地测量,分析数据质量,对比解释成果.比较二者在地面瞬变电磁勘探应用中的优劣.得到结论:1)"8"字形回线比同一回线有更好的信噪比.2)"8"字形回线装置在抗干扰能力优于同一回线装置,晚期信号的稳定性优于同一回线装置.3)关断时间相同的情况下,同一回线比"8"字形回线有更好的分辨能力,能获得更加丰富的浅层地质信息.4)在信号的灵敏度反应方面,"8"字形回线不如同一回线装置.5)在深大断裂构造或者目标体跟围岩电阻率差距较大的地区可以采用"8"字形回线.  相似文献   

10.
瞬变电磁(TEM)视电阻率早、晚期计算公式在求取时进行了简化近似,造成部分地电信息的失真.瞬变电磁响应信号衰减过程与介质的电性特征有关,由此引入瞬变电磁响应微分衰减速度和二阶微分衰减速度两个解释参数,力图避开视电阻率计算,以衰减速度表示瞬变电磁探测结果的电性特征.首先研究了不同导电均匀半空间与低阻目标体地电模型瞬变电磁响应微分衰减速度在时间上的变化特征,在此基础上对瞬变电磁响应二阶微分衰减速度进行了分析,通过绘制TEM二阶微分衰减拟断面图实现了低阻目标体地电模型的电性层划分.研究结果表明:电性分界面与二阶微分衰减的"0"等值线相对应,可以用"0"值等值线作为电性分界面的划分依据,即瞬变电磁微分解释方法能够区分不同地电断面,可以用于瞬变电磁探测成果解释.  相似文献   

11.
A three-dimensional finite-element time-domain forward-modelling algorithm is developed to simulate transient electromagnetics excited by grounded-wire sources. The main advantage of this finite-element time-domain algorithm is that full transmitting-current waveforms and complex-shaped sources resulting from topography can be directly dealt with in this algorithm. The models used to test this algorithm include a homogeneous half-space model, a stratified-medium model, the model of a complex conductor at a vertical contact and the Ovoid Zone massive sulfide deposit at Voisey's Bay, Canada. The homogeneous half-space model is used to determine the truncation boundary for a computational domain, and to compare with the electromagnetic responses excited by step-off, step-on and direct current waveforms. For the stratified-medium model, results demonstrate that full transmitting waveforms have strong effects on the observed electromagnetic responses. The model of a complex conductor at a vertical contact is designed for the grounded electrical source airborne transient electromagnetic method and is also used to examine the effectiveness of the broadside and inline configurations for such a vertical, thin plate embedded in the subsurface. The area of the Ovoid Zone massive sulfide deposit possesses non-negligible topography, the effects of which on the shapes of the grounded-wire sources must be taken into account when implementing the finite-element time-domain solver. The results show that both the broadside and inline electromagnetic responses are strongly affected by the massive conductive ore body.  相似文献   

12.
This paper presents a computational method for the interpretation of electromagnetic (EM) profile data in the frequency domain using a thin plate model within a two-layer earth. The modelling method is based on an integral equation formulation, where the conductor is represented by a lattice structure composed of two-dimensional surface elements. Several approximations are used to simplify the theoretical basis and to decrease the computation time. The simple parametric model allows efficient use of optimization methods. We employ a linearized inversion scheme based on singular value decomposition and adaptive damping. The new forward computation method and the parameter optimization are combined in the computer program, emplates . The modelling examples demonstrate that the approximate method is capable of describing the characteristic behaviour of the EM response of a thin plate-like conductor in conductive surroundings. The efficacy of the inversion is demonstrated using both synthetic and field data. An optional depth compensation method is used to improve the interpreted values of the depth of burial. The results show that the method is cost effective and suitable for interactive interpretation of EM data.  相似文献   

13.
A survey using Crone pulse electromagnetic (PEM) equipment was conducted near the village of Gani in Andhra Pradesh, India. The anomalies obtained had a shape typical for an inclined sheet conductor. From the multichannel TEM response, the multifrequency Slingram response was calculated. A comparison of the Slingram response of the Gani conductors with the free-air model response of a thin sheet conductor shows phase rotation and attenuation of the anomaly vector, thereby confirming the presence of a conductive overburden. Interpretation of the Gani results using free-air phasor diagrams causes an overestimation of depth and conductance of the target conductor unless data obtained at a sufficiently low frequency are used for the interpretation. For areas like Gani this frequency can be as high as 336 Hz.  相似文献   

14.
瞬变电磁隧道超前预报成像技术   总被引:17,自引:3,他引:14       下载免费PDF全文
为解决隧道掌子面前含水带病害快速有效探测问题,尝试把对含水带结构反映敏感的瞬变电磁法引入到隧道掌子面前进行工作.文中介绍了瞬变电磁超前预报的工作装置形式;通过对掌子面特定环境的分析,提出引用“浮动薄板”理论,以二次电导微分参数为特征量建立隧道超前预报成像系统.推导出以等效导电薄板为虚拟像源的磁场响应与电导之间的非线性关系式. 通过引入辅助函数,采用遗传算法求得电导参数.最终以二次电导微分参数绘制成像剖面.对地电模型进行数值模拟和对应用实例进行了成像,结果表明成像方法对隧道掌子面前方水体病害预报效果明显.  相似文献   

15.
地磁Z分量日变"低点位移"异常成因分析   总被引:2,自引:0,他引:2  
胡久常 《地震》2011,31(2):79-86
本文基于对地磁Z分量日变"低点位移"异常和太阳静日变化内外源场等效涡旋电流体系时空演化特征的分析,提出了"电磁炉效应"假说,以解释地磁Z分量日变"低点位移"异常的成因.本文认为,起源于地球内部和外部涡旋电流体系所产生的太阳静日变化Sq,作为一个持续的交变磁场,切割地壳深部的高导低速层,形成涡旋电流并产生热能,致使高导低...  相似文献   

16.
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.  相似文献   

17.
Theoretically exact type curves for a semi-infinite thin conductor are presented for various dip, angles, depth of burial and conductance. The study shows that the common-loop response shape is sensitive to small changes in conductor dip, but is affected more subtly by comparable strike variations. For large sheet conductors a decrease in the strike angle results in a broadening but unlike that for a finite plate there is no reduction in peak amplitude. For dipping conductors, response asymmetry and the direction and magnitude of peak amplitude displacement can be used to assess the disposition and quality of the conductor. A generalized interpretation scheme is proposed, based on dimensionless response characteristics and normalized decay curves, to facilitate the rapid in-field determination of conductor dip, conductance and depth of burial, for any time regime.  相似文献   

18.
To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional (2D) timedomain transient electromagnetic secondary field of the line source based on the DuFort-Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force (Emf), and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the flat layered model with high precision—the maximum relative error is less than 0.01% between our method and the analytical method—and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.  相似文献   

19.
The magnetometric resistivity (MMR) method uses a sensitive magnetometer to measure the low-level, low-frequency magnetic fields associated with the galvanic current flow between a pair of electrodes. While the MMR anomalies of simple structures such as dikes and vertical contacts have been determined analytically, there is a lack of systematic information on the expected responses from simple three-dimensional bodies. We determine the characteristic anomalies associated with square, plate-like conductors, which are excellent models of many base metal mineral deposits. The anomalies of plates of finite size are determined numerically using an integral equation method. A plate is subdivided into many sections and the current flow within each section is solved by equating the electrical field within each section to the tangential electrical field just outside it. When the plate size is small in relation to either the depth or the transmitter spacing, the shape and amplitude of the anomaly produced is closely approximated by a current dipole model of the same length and depth. At the other extreme, a large plate is represented by a half-plane. The dipole and half-plane models are used to bracket the behaviour of plates of finite size. The form of a plate anomaly is principally dependent on the shape, depth and orientation of the plate. A large, dipping plate near the surface produces a skewed anomaly highly indicative of its dip, but the amount of skew rapidly diminishes with increased depth or decreased size. Changes in plate conductivity affect the amplitude of the anomaly, but have little effect on anomaly shape. A current channelling parameter, determined from the conductivity contrast, can thus be used to scale the amplitude of an anomaly whose basic shape has been determined from geometrical considerations. The separation into geometrical and electrical factors greatly simplifies both the interpretation and modelling of MMR anomalies, particularly in situations with multiple plates. An empirical formula, using this separation, predicts the anomaly of two or more parallel plates with different conductances. In addition, the relation between the resolution of two vertical, parallel plates of equal conductance and their separation is determined. The ability of the integral equation method to model plate-like structures is demonstrated with the interpretation of an MMR anomaly in a survey conducted at Cork Tree Well in Western Australia. The buried conductor, a mineralized graphitic zone, is modelled with a vertical, bent plate. The depth to the top of the plate, and the plate conductance, is adjusted to fit the anomaly amplitude as closely as possible. From the modelling it would appear that this zone is not solely responsible for the observed anomaly.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号