首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
Electromagnetic fields in a steel-cased borehole   总被引：1，自引：0，他引：1  
The development of an electromagnetic numerical modelling scheme for a magnetic dipole in an arbitrary casing segment in an inhomogeneous conductivity background has been difficult, due to the very high electrical conductivity and magnetic permeability contrasts between the steel casing and the background medium. To investigate the effect of steel casing efficiently, we have developed an accurate but simple finite‐element modelling scheme to simulate electromagnetic fields in a medium of cylindrically symmetric conductivity structures. In order to preserve the cylindrical symmetry in the resulting electromagnetic fields, a horizontal loop current source is used throughout. One of the main advantages of the approach is that the problem is scalar when formulated using the azimuthal electric field, even if the casing is both electrically conductive and magnetically permeable. Field calculations have been made inside the cased borehole as well as in another borehole which is not cased. Careful analyses of the numerical modelling results indicate that the anomaly observed in a cross‐borehole configuration is sensitive enough to be used for tomographic imaging.  相似文献   

2.
Very early times in the order of 2–3 μs from the end of the turn‐off ramp for time‐domain electromagnetic systems are crucial for obtaining a detailed resolution of the near‐surface geology in the depth interval 0–20 m. For transient electromagnetic systems working in the off time, an electric current is abruptly turned off in a large transmitter loop causing a secondary electromagnetic field to be generated by the eddy currents induced in the ground. Often, however, there is still a residual primary field generated by remaining slowly decaying currents in the transmitter loop. The decay disturbs or biases the earth response data at the very early times. These biased data must be culled, or some specific processing must be applied in order to compensate or remove the residual primary field. As the bias response can be attributed to decaying currents with its time constantly controlled by the geometry of the transmitter loop, we denote it the ‘Coil Response’. The modelling of a helicopter‐borne time‐domain system by an equivalent electronic circuit shows that the time decay of the coil response remains identical whatever the position of the receiver loop, which is confirmed by field measurements. The modelling also shows that the coil response has a theoretical zero location and positioning the receiver coil at the zero location eliminates the coil response completely. However, spatial variations of the coil response around the zero location are not insignificant and even a few cm deformation of the carrier frame will introduce a small coil response. Here we present an approach for subtracting the coil response from the data by measuring it at high altitudes and then including an extra shift factor into the inversion scheme. The scheme is successfully applied to data from the SkyTEM system and enables the use of very early time gates, as early as 2–3 μs from the end of the ramp, or 5–6 μs from the beginning of the ramp. Applied to a large‐scale airborne electromagnetic survey, the coil response compensation provides airborne electromagnetic methods with a hitherto unseen good resolution of shallow geological layers in the depth interval 0–20 m. This is proved by comparing results from the airborne electromagnetic survey to more than 100 km of Electrical Resistivity Tomography measured with 5 m electrode spacing.  相似文献   

3.
In present‐day land and marine controlled‐source electromagnetic (CSEM) surveys, electromagnetic fields are commonly generated using wires that are hundreds of metres long. Nevertheless, simulations of CSEM data often approximate these sources as point dipoles. Although this is justified for sufficiently large source‐receiver distances, many real surveys include frequencies and distances at which the dipole approximation is inaccurate. For 1D layered media, electromagnetic (EM) fields for point dipole sources can be computed using well‐known quasi‐analytical solutions and fields for sources of finite length can be synthesized by superposing point dipole fields. However, the calculation of numerous point dipole fields is computationally expensive, requiring a large number of numerical integral evaluations. We combine a more efficient representation of finite‐length sources in terms of components related to the wire and its end points with very general expressions for EM fields in 1D layered media. We thus obtain a formulation that requires fewer numerical integrations than the superposition of dipole fields, permits source and receiver placement at any depth within the layer stack and can also easily be integrated into 3D modelling algorithms. Complex source geometries, such as wires bent due to surface obstructions, can be simulated by segmenting the wire and computing the responses for each segment separately. We first describe our finite‐length wire expressions and then present 1D and 3D examples of EM fields due to finite‐length sources for typical land and marine survey geometries and discuss differences to point dipole fields.  相似文献   

4.
We study a new marine electromagnetic configuration that consists of a ship‐towed inductive source transmitter and a series of remote electric dipole receivers placed on the seafloor. The approach was tested at the Palinuro Seamount in the southern Tyrrhenian Sea, at a site where massive sulphide mineralization has been previously identified by shallow drilling. A 3D model of the Palinuro study area was created using bathymetry data, and forward modelling of the electric field diffusion was carried out using a finite volume method. These numerical results suggest that the remote receivers can theoretically detect a block of shallowly buried conductive material at up to ～100 m away when the transmitter is located directly above the target. We also compared the sensitivity of the method using either a horizontal loop transmitter or a vertical loop transmitter and found that when either transmitter is located directly above the mineralized zone, the vertical loop transmitter has sensitivity to the target at a farther distance than the horizontal loop transmitter in the broadside direction by a few tens of metres. Furthermore, the vertical loop transmitter is more effective at distinguishing the seafloor conductivity structure when the vertical separation between transmitter and receiver is large due to the bathymetry. As a horizontal transmitter is logistically easier to deploy, we conducted a first test of the method with a horizontal transmitter. Apparent conductivities are calculated from the electric field transients recorded at the remote receivers. The analysis indicates higher apparent seafloor conductivities when the transmitter is located near the mineralized zone. Forward modelling suggests that the best match to the apparent conductivity data is obtained when the mineralized zone is extended southward by 40 m beyond the zone of previous drilling. Our results demonstrate that the method adds value to the exploration and characterization of seafloor massive sulphide deposits.  相似文献   

5.
Lake sediments may serve as archives on paleoclimatic fluctuations, geomagnetic field variations and volcanic activities. Lake Holzmaar in Eifel/Germany is a maar lake and its lacustrine sediments provide paleoclimatic proxy data. Therefore, knowledge about the geometry and, especially, about the thickness of the sediments is very important for determining an optimum drilling location for paleoclimatic studies. We have developed a floating in‐loop transient electromagnetic method field set up (Float‐transient electromagnetic method) with a transmitter and receiver size of 18 × 18 m2 and 6 × 6 m2 respectively. This special set up enables in‐loop transient electromagnetic method measurements on the surface of freshwater lakes that define the geometry and the thickness of sediments beneath such lakes thus helping to determine optimum drilling locations. Due to the modular design of the new Float‐transient electromagnetic method field set up, this system can be handled by two operators and can easily be transported. Sixteen in‐loop soundings were carried out on the surface of Lake Holzmaar. The transient electromagnetic method data could not be interpreted by conventional 1D inversions because of the 3D distribution of subsurface conductivity caused by the lake's geometry. Three‐dimensional finite element modelling was applied to explain the observed transients and the 3D conductivity distribution beneath the lake was recovered by taking its geometry into account. The 3D interpretation revealed approximately 55 m thick sediments beneath 20 m deep water in the central part of the lake.  相似文献   

6.
The electromagnetic response of a horizontal electric dipole transmitter in the presence of a conductive, layered earth is important in a number of geophysical applications, ranging from controlled‐source audio‐frequency magnetotellurics to borehole geophysics to marine electromagnetics. The problem has been thoroughly studied for more than a century, starting from a dipole resting on the surface of a half‐space and subsequently advancing all the way to a transmitter buried within a stack of anisotropic layers. The solution is still relevant today. For example, it is useful for one‐dimensional modelling and interpretation, as well as to provide background fields for two‐ and three‐dimensional modelling methods such as integral equation or primary–secondary field formulations. This tutorial borrows elements from the many texts and papers on the topic and combines them into what we believe is a helpful guide to performing layered earth electromagnetic field calculations. It is not intended to replace any of the existing work on the subject. However, we have found that this combination of elements is particularly effective in teaching electromagnetic theory and providing a basis for algorithmic development. Readers will be able to calculate electric and magnetic fields at any point in or above the earth, produced by a transmitter at any location. As an illustrative example, we calculate the fields of a dipole buried in a multi‐layered anisotropic earth to demonstrate how the theory that developed in this tutorial can be implemented in practice; we then use the example to examine the diffusion of volume charge density within anisotropic media—a rarely visualised process. The algorithm is internally validated by comparing the response of many thin layers with alternating high and low conductivity values to the theoretically equivalent (yet algorithmically simpler) anisotropic solution, as well as externally validated against an independent algorithm.  相似文献   

7.
Steel well casings in or near a hydrocarbon reservoir can be used as source electrodes in time‐lapse monitoring using grounded line electromagnetic methods. A requisite component of carrying out such monitoring is the capability to numerically model the electromagnetic response of a set of source electrodes of finite length. We present a modelling algorithm using the finite‐element method for calculating the electromagnetic response of a three‐dimensional conductivity model excited using a vertical steel‐cased borehole as a source. The method is based on a combination of the method of moments and the Coulomb‐gauged primary–secondary potential formulation. Using the method of moments, we obtain the primary field in a half‐space due to an energized vertical steel casing by dividing the casing into a set of segments, each assumed to carry a piecewise constant alternating current density. The primary field is then substituted into the primary–secondary potential finite‐element formulation of the three‐dimensional problem to obtain the secondary field. To validate the algorithm, we compare our numerical results with: (i) the analytical solution for an infinite length casing in a whole space, excited by a line source, and (ii) a three‐layered Earth model without a casing. The agreement between the numerical and analytical solutions demonstrates the effectiveness of our algorithm. As an illustration, we also present the time‐lapse electromagnetic response of a synthetic model representing a gas reservoir undergoing water flooding.  相似文献   

8.
Over the past decade the typical size of airborne electromagnetic data sets has been growing rapidly, along with an emerging need for highly accurate modelling. One‐dimensional approximate inversions or data transform techniques have previously been employed for very large‐scale studies of quasi‐layered settings but these techniques fail to provide the consistent accuracy needed by many modern applications such as aquifer and geological mapping, uranium exploration, oil sands and integrated modelling. In these cases the use of more time‐consuming 1D forward and inverse modelling provide the only acceptable solution that is also computationally feasible. When target structures are known to be quasi layered and spatially coherent it is beneficial to incorporate this assumption directly into the inversion. This implies inverting multiple soundings at a time in larger constrained problems, which allows for resolving geological layers that are undetectable using simple independent inversions. Ideally, entire surveys should be inverted at a time in huge constrained problems but poor scaling properties of the underlying algorithms typically make this challenging. Here, we document how we optimized an inversion code for very large‐scale constrained airborne electromagnetic problems. Most importantly, we describe how we solve linear systems using an iterative method that scales linearly with the size of the data set in terms of both solution time and memory consumption. We also describe how we parallelized the core region of the code, in order to obtain almost ideal strong parallel scaling on current 4‐socket shared memory computers. We further show how model parameter uncertainty estimates can be efficiently obtained in linear time and we demonstrate the capabilities of the full implementation by inverting a 3327 line km SkyTEM survey overnight. Performance and scaling properties are discussed based on the timings of the field example and we describe the criteria that must be fulfilled in order to adapt our methodology for similar type problems.  相似文献   

9.
电磁场数值模拟的背景场/异常场算法是三维正演的有效策略之一,优点为采用解析法计算电磁场背景场代替场源项、克服了场源奇异性,缺点为不适用于发射源布置于起伏地表或背景模型复杂的情形.总场算法是直接对电磁场总场开展数值模拟,其难点是有效加载场源、保证近区与过渡区数值解精度.本文以水平电偶源形式分段加载接地长导线源,并以电场总场Helmholtz方程为矢量有限元法控制方程,实现了基于非结构化四面体网格剖分的接地长导线源频率域电磁法三维正演.通过与均匀全空间中水平电偶源产生的电场解析解对比,验证了本文算法的正确性,并分析了四面体外接圆半径与其最短棱边的最大比值和四面体二面角最小值对数值解精度的影响规律.通过与块状高导体地电模型的积分方程法、有限体积法和基于磁矢量势Helmholtz方程的有限元法数值解对比,进一步验证了本文算法正确性,同时说明了非结构化四面体网格能够更加精细地剖分电性异常体,利于获得精确数值解.  相似文献   

10.
Calculations of emf induced in a vertical axis coil by a current step in a straight wire segment of finite length were carried out using an equivalent dipole placed along the source wire. An approximation is valid for homogeneous and layered earth models. The location of the equivalent dipole is calculated by forcing the early- and late-time asymptotes of the transient sounding curve of the equivalent dipole to match those of the finite wire. This approximation works because the early-time asymptote of the emf depends on the component of the receiver position which is parallel to the wire while the late-time asymptote does not. Analytical integration of the early-time asymptote yields an expression for the equivalent dipole location in closed form. The coil can be placed anywhere in the vicinity of the finite wire. Square or rectangular loop sources can be simulated by one or more finite wire segments depending on the source-receiver geometry. The equivalent radius calculated for central loop soundings agrees well with the value derived using a circular loop with the same area as the square loop. Results show that acceptable sounding curves can be generated by the equivalent dipole for coils placed as close as 0.5 source lengths from the finite wire segment. Higher accuracy can be obtained by splitting the finite wire into two or more subsegments. Results for layered models are slightly better than homogeneous earth models when the resistivity increases with depth and slightly worse for models with resistivity decreasing with depth. Approximate calculations are about 10 times less expensive than exact calculations depending on the method used for the numerical integration.  相似文献   

11.
固定翼时间域航空电磁探测系统在实际飞行测量过程中,发射线圈、接收线圈姿态和吊舱摆动状态不断变化,在测量数据中引入如发射磁矩方向、接收分量方向以及系统收发距等参数的误差,影响数据反演成像效果.本文基于固定翼时间域航空电磁正演理论,利用姿态变换,引入发射线圈、接收线圈双旋转矩阵;根据发射、接收线圈相对位置的几何关系,求得摆动格林张量;推导了任意姿态角度以及任意摆动角度情况下的固定翼航空电磁响应三分量计算表达式.通过层状大地模型的仿真计算,分别研究了发射、接收线圈各姿态以及吊舱摆动状态对航空电磁响应的影响,得出发射线圈、接收线圈俯仰旋转和吊舱同向摆动对系统电磁响应影响最强;仿真分析了实际测量中,三种角度同时存在情况下,航空电磁响应的定量变化规律.在此基础上,讨论了响应系数与大地电导率的关系,同时给出基于响应系数的固定翼航空电磁系统线圈姿态和摆动状态校正方法,准二维层状大地模型反演结果表明,校正后数据的反演精度提高了33.1%.  相似文献   

12.
In near-surface transient electromagnetic studies, it is desirable to measure the transient response starting from the earliest possible time. This requires the current in the transmitter loop to be switched off quickly, which necessitates working with a low transmitter current. As for deep-target transient electromagnetic studies, the transmitter current is as high as possible. The transmitter current's turn-off waveform and total duration affect the transient voltage response, especially at early times, which is to be accounted for when interpreting transient electromagnetic data. This article discusses the difference in switching off low and high current in a horizontal loop used as the source of the primary magnetic field in the transient electromagnetic method. Low and high currents are turned off in fundamentally different ways. When the current to be switched off is low, the loop can be represented as a symmetric combination of two transmission lines grounded at the middle of the loop perimeter. Such a representation of a loop allows calculating the current turn-off waveform at any point of the loop. The waveform and total duration of switching off a low current does not depend on its magnitude, but is determined by the period of natural oscillations of the current in the loop and the resistance of a shunting resistor. Switching off a low current in a loop can be represented as the sum of stepped current waves travelling along the loop wire. As a consequence, the current at different points of the loop perimeter is turned off at different times. In contrast to a low current, a high current is switched off linearly in time and synchronously at all points of the loop perimeter. The wave phenomena appear only at the very beginning of the current shutdown for a time interval that is much less than the total current turn-off duration. Presentation of the loop using a simple lumped-circuit model predicts the waveform and duration of the high current turn-off that coincide with the measured ones. There are two reasons why the article may be of interest to those engaged in the theory and/or practice of electromagnetic geophysical methods. First, it contributes to a general understanding of how the current in the transmitter loop is turned off. Second, the article shows how the parameters of a transmitter loop determine the current turn-off duration and thus the minimum depth of the transient electromagnetic sounding method.  相似文献   

13.
Time‐domain marine controlled source electromagnetic methods have been used successfully for the detection of resistive targets such as hydrocarbons, gas hydrate, or marine groundwater aquifers. As the application of time‐domain marine controlled source electromagnetic methods increases, surveys in areas with a strong seabed topography are inevitable. In these cases, an important question is whether bathymetry information should be included in the interpretation of the measured electromagnetic field or not. Since multi‐dimensional inversion is still not common in time‐domain marine controlled source electromagnetic methods, bathymetry effects on the 1D inversion of single‐offset and multi‐offset joint inversions of time‐domain controlled source electromagnetic methods data are investigated. We firstly used an adaptive finite element algorithm to calculate the time‐domain controlled source electromagnetic methods responses of 2D resistivity models with seafloor topography. Then, 1D inversions are applied on the synthetic data derived from marine resistivity models, including the topography in order to study the possible topography effects on the 1D interpretation. To evaluate the effects of topography with various steepness, the slope angle of the seabed topography is varied in the synthetic modelling studies for deep water (air interaction is absent or very weak) and shallow water (air interaction is dominant), respectively. Several different patterns of measuring configurations are considered, such as the systems adopting nodal receivers and the bottom‐towed system. According to the modelling results for deep water when air interaction is absent, the 2D topography can distort the measured electric field. The distortion of the data increases gradually with the enlarging of the topography's slope angle. In our test, depending on the configuration, the seabed topography does not affect the 1D interpretation significantly if the slope angle is less or around 10°. However, if the slope angle increases to 30° or more, it is possible that significant artificial layers occur in inversion results and lead to a wrong interpretation. In a shallow water environment with seabed topography, where the air interaction dominates, it is possible to uncover the true subsurface resistivity structure if the water depth for the 1D inversion is properly chosen. In our synthetic modelling, this scheme can always present a satisfactory data fit in the 1D inversion if only one offset is used in the inversion process. However, the determination of the optimal water depth for a multi‐offset joint inversion is challenging due to the various air interaction for different offsets.  相似文献   

14.
Artificial neural networks were used to implement an automatic inversion of frequency‐domain airborne electromagnetic (AEM) data that do not require a priori information about the survey area. Two classes of model, i.e. homogeneous half‐space models and horizontally layered half‐space models with two layers, are used in this 1D inversion, and for each data point the selection of the class of 1D model is performed prior to the inversion, also using an artificial neural network. The proposed inversion method was tested in a survey area situated in Austria, northwest of Vienna in the Bohemian Massif. The results of the inversion were compared with the geological setting, logging results, and seismic and gravimetric measurements. This comparison shows a good correlation between the AEM models and the known geological and geophysical data.  相似文献   

15.
A fully three‐dimensional finite‐element algorithm has been developed for simulating controlled‐source electromagnetic surveys. To exploit the advantages of geometric flexibility, frequency‐domain Maxwell's equations of the secondary electric field were discretised using edge‐based finite elements while the primary field was calculated analytically for a horizontally layered‐earth model. The resulting system of equations for the secondary field was solved using a parallel version of direct solvers. The accuracy of the algorithm was successfully verified by comparisons with integral‐equations and iterative solutions, and the applicability to models containing large conductivity contrasts was verified against published data. The advantages of geometry‐conforming meshes have been demonstrated by comparing different mesh systems to simulate an inclined sheet model. A comparison of the performance between direct and iterative solvers demonstrated the superior efficiency of direct solvers, particularly for multisource problems.  相似文献   

16.
Accurate modelling of the conductivity structure of mineralisations can often be difficult. In order to remedy this, a parametric approach is often used. We have developed a parametric thin‐sheet code, with a variable overburden. The code is capable of performing inversions of time‐domain airborne electromagnetic data, and it has been tested successfully on both synthetic data and field data. The code implements an integral solution containing one or more conductive sheets, buried in a half‐space with a laterally varying conductive overburden. This implementation increases the area of applicability compared to, for example, codes operating in free space, but it comes with a significant increase in computational cost. To minimise the cost, the code is parallelised using OpenMP and heavily optimised, which means that inversions of field data can be performed in hours on multiprocessor desktop computers. The code models the full system transfer function of the electromagnetic system, including variable flight height. The code is demonstrated with a synthetic example imitating a mineralisation buried underneath a conductive meadow. As a field example, the Valen mineral deposit, which is a graphite mineral deposit located in a variable overburden, is successfully inverted. Our results match well with previous models of the deposit; however, our predicted sheet remains inconclusive. These examples collectively demonstrate the effectiveness of our thin‐sheet code.  相似文献   

17.
全波形时间域航空电磁响应三维有限差分数值计算   总被引：11，自引：8，他引：3       下载免费PDF全文
全波形时域航空电磁测量具有解释精度高、分辨率高等优点,已成为新一代时域航空电磁探测发展的趋势.为此,本文基于无源有限差分法,将发射电流离散为若干梯形脉冲元,将脉冲元在时间上产生的电磁场转换为空间初始场,在差分迭代过程中分时引入差分方程,逐步将有源有限差分方程转换为无源有限差分方程,最终实现了发射电流为任意波形时,全波形时间域航空电磁响应的三维数值计算.计算结果表明:在梯形波发射时,对于典型地电模型的航空时间域电磁响应,计算稳定时间大于10 ms;与积分方程法相比,发射线圈中心点二次场响应平均误差小于1%.任意源三维全波形数值计算的实现为全波形三维反演和仪器设计奠定了基础.  相似文献   

18.
Theoretical and physical scale model responses for moving source and fixed loop electromagnetic exploration systems are presented for the case of a vertical, thin, plate conductor located first in free-space and then in a conductive host. The results are presented in the form of anomaly index diagrams plotted in complex space. The two modeling approaches show general agreement on the behavior of the current gathering effect that operates when the conductive host is present. Both approaches show an onset of a strong positive effect at higher frequencies in the responses provided by both horizontal coplanar coil systems and by fixed loop systems. Agreement is also shown on an unanticipated negative effect that a conductive host causes in the responses provided by the vertical coplanar coincident coil system. Both modeling approaches demonstrate that the responses provided by fixed loop systems are predominantly caused by the current gathering effect. The extent of this effect for fixed loop systems depends on the size of the transmitter loop with respect to the target and on the depth of the target. The theoretical modeling demonstrates that both moving source and fixed loop systems operated over a target located in a conductive host, provide responses that become almost identical and independent of the conductance of the target at high frequency.  相似文献   

19.
We recently proposed an efficient hybrid scheme to absorb boundary reflections for acoustic wave modelling that could attain nearly perfect absorptions. This scheme uses weighted averaging of wavefields in a transition area, between the inner area and the model boundaries. In this paper we report on the extension of this scheme to 2D elastic wave modelling with displacement‐stress formulations on staggered grids using explicit finite‐difference, pseudo‐implicit finite‐difference and pseudo‐spectral methods. Numerical modelling results of elastic wave equations with hybrid absorbing boundary conditions show great improvement for modelling stability and significant absorption for boundary reflections, compared with the conventional Higdon absorbing boundary conditions, demonstrating the effectiveness of this scheme for elastic wave modelling. The modelling results also show that the hybrid scheme works well in 2D rotated staggered‐grid modelling for isotropic medium, 2D staggered‐grid modelling for vertically transversely isotropic medium and 2D rotated staggered‐grid modelling for tilted transversely isotropic medium.  相似文献   

20.
本文实现了一种面向目标自适应海洋可控源电磁三维矢量有限元方法.为满足三维复杂电性结构模拟的需求,网格剖分采用非结构化六面体.在组装刚度矩阵之后,形成的大型复数线性方程组分解为等价的实数形式,利用带预条件的广义最小残差法进行求解.在获得微分方程的解之后,为提高解的准确性,通过面向目标的自适应误差估计来指示网格细化,重点加密能使观测点数值模拟精度提高的网格.对于大规模三维数据,为了使模型空间的并行计算达到均衡负载的效果,我们使用METIS函数库来进行网格计算任务量的划分.最后,通过对比一维解析解与三维自适应矢量有限元计算结果,验证了程序的正确性;通过自适应过程中误差指示子的分布,验证了面向目标自适应的有效性;通过对三维复杂模型进行均衡负载下的并行计算,测试了程序的可扩展性.  相似文献   

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

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