首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
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.  相似文献   

2.
Airborne EM skin depths   总被引：1，自引：0，他引：1  
Skin depth is an electromagnetic (EM) scale length that provides a measure of the degree of attenuation experienced by a particular frequency of an EM system. As has been discussed in the literature, skin depth is not a complete measure of the depth of investigation, but the two may be related. Frequency‐domain airborne EM systems employ pairs of transmitter and receiver coils that use a frequency range from several hundred hertz to over 100 kHz. For elevated dipoles, both geometrical and frequency‐dependent attenuation of the induced fields must be considered. For airborne EM systems it is possible to define a skin depth based only on the electric field induced by the transmitter. A vertical decay scale length, here defined from the at‐surface position of maximum electric field, enables the same skin‐depth estimate to be obtained for both cases of vertical and horizontal dipolar excitation. Such dipolar skin depths associated with towed‐bird and fixed‐wing airborne systems are studied in relation to frequency, conductivity and sensor elevation. Dipolar skin depths are found to be much smaller than their plane‐wave counterparts except at high frequency (>50 kHz) and in combination with high conductivity. For the majority of airborne systems the influence of altitude on skin depth is highly significant. Dipolar skin depths increase with increasing sensor elevation. Low frequencies display the greatest sensitivity. At low elevation (<40 m), geometrical attenuation dominates the behaviour of the skin depth. The study indicates that typical low‐altitude airborne surveys provide vertically compact assessments of subsurface conductivity, well suited to near‐surface, environmental applications.  相似文献   

3.
为研究双轴各向异性介质多分量感应测井响应特征,本文基于三重傅里叶变换,推导任意方向偶极子源的谱域电磁场解析式;采用围线积分方法,自适应截断积分区间,结合谱域电磁场周期特性,实现三重傅里叶变换的精确快速积分;进而,针对双轴各向异性倾斜地层,模拟研究不同纵横向各向异性条件多分量感应测井响应特征.结果表明:利用谱域内电磁场在周向的周期特性简化解析式,可将计算速度提高4倍;自适应截断积分区间方法保证了计算精度,并极大地减少了积分节点数.对于倾斜双轴各向异性介质,倾角较大时,共面分量可反映地层横向各向异性,同轴分量可反映地层纵向各向异性;倾角较小时,同轴分量可反映地层横向各向异性,共面分量可反映地层纵向各向异性.  相似文献   

4.
The Millennium uranium deposit is located within the Athabasca Basin of northern Saskatchewan. The basement rocks, comprised primarily of paleo‐Proterozoic gneisses, are electrically resistive. However, the deposit is associated with highly conductive graphitic metasediments that are intercalated with the gneisses. An unconformity separates the basement rocks from the overlying, horizontally stratified, Proterozoic sandstones of the Athabasca Group (which are also highly resistive). The strike extents of the graphitic metasedimenary packages are extensive and therefore electromagnetic (EM) survey techniques are successful at identifying these zones but do not identify the specific locations where they are enriched in uranium. Through drilling it has been noted that hydrothermal processes associated with mineralization has altered the rocks in the vicinity of the deposits, which should in theory result in a resistivity low. A significant resistivity low has been mapped coincident with the Millennium deposit using ground resistivity survey techniques. However, a comparison of the airborne EM and ground resistivity results reveals that the two data sets have imaged different features. The resistive‐limit (on‐time) windows of the MEGATEM data show conductive features corresponding to lakes located to the west and south of the deposit. The late‐time windows show a feature to the east of the deposit, interpreted as being associated with the east‐dipping graphitic basement conductors (similar to that observed in historical ground EM data collected in this area). The early‐time TEMPEST windows (delay times less than 0.2 ms) show a broad resistivity low located at approximately the same location as where the alteration has been identified through drilling. Modelling the data is not easy but a response that decays prior to 0.3 ms is consistent with 500 Ωm material in the sandstone, a resistivity value close to the lower limit with respect to the hydrothermally altered Athabasca group sediments in this area. The MEGATEM system does not see a conductive zone over the alteration as clearly but the high signal‐to‐noise ratio in the late‐time MEGATEM data means that the conductive material at a greater depth is more coherently imaged.  相似文献   

5.
In this paper, the noise sources of an airborne electromagnetic frequency domain instrument used to measure sea-ice thickness are studied. The antennas are mounted on the wings of an aircraft. The paper presents real data examples showing that strong noise limited the accuracy of the thickness measurement to ± 0.5 m in the best case. Even drift cor­rection and frequency ?ltering did not reduce the noise to a level necessary for sea ice thickness measurements with an accuracy of 0.1 m. We show results of 3D ?nite element modeling of the coupling between transmitter and receiver coils and the aircraft, which indicate that wing ?exure is the primary cause of the strong noise. Wing de?ection angles below 5° relative to the fuselage are large enough to cause changes higher than the wanted signal from the seawater under the ice. Wing ?exure noise can be divided into an inductive and geometric contribution, both of the same order. Most of the wing ?exure signal appears on the inphase component only, hence the quadrature component should be taken for sea ice thickness retrievals when wing ?exure is present even when the inphase produces a larger ocean sig­nal. Results also show that pitch and roll movements of the aircraft and electromagnetic coupling between seawater and aircraft can contribute signi?cantly to the total noise. For ?ight heights of 30 m over the ocean these effects can change the sig­nal by about 10% or more. For highly quantitative measurements like sea-ice thickness all these effects must be taken into account. We conclude that a ?xed wing electromagnetic instrument for the purpose of measure­ments in a centimeter scale must include instrumentation to measure the relative position of the antenna coils with an accuracy of 1/10 mm. Furthermore the antenna separation distance should be as large as possible in order to increase the measured ratio of secondary to primary magnetic field strength.  相似文献   

6.
The diffusion of electromagnetic fields is dependent not only on conductivity, but also on magnetic permeability, dielectric permittivity and polarizability, i.e. dispersive conductivity. The long‐offset transient electromagnetic (LOTEM) method is mainly used to determine the spatial distribution of conductivity in the subsurface. However, earlier work on loop‐loop TEM suggests that transient EM methods can also be affected by induced polarization (IP). Numerous 1D forward calculations were carried out to study the IP effect on LOTEM data, using the Cole‐Cole relaxation model to simulate the polarizability of the ground. Besides the polarizability of each layer, the IP effect depends on the LOTEM field set‐up and the spatial distribution of conductivity in the ground. In particular, near‐surface layers with high chargeabilities can significantly distort the late time transients of the electric field components in the vicinity of the transmitter. The influence of polarizable layers on the magnetic field components can be neglected under normal circumstances. In 1997 and 1999, LOTEM measurements were carried out at Mt. Vesuvius in Italy to explore the geological structure of the volcano. Sensitivity studies on the effect of polarizable layers suggest that high chargeabilities in connection with conductive layers at greater depths would result in a detectable distortion of the electric field transients. Although the simultaneous IP measurements revealed high chargeabilities in a near‐surface layer, no evidence of IP effects could be found in the measured LOTEM data. We conclude that the observed chargeabilities are local and that 3D effects are probably present in the data. Another aspect is the measurement of the system response, which is usually measured by placing a receiver very close to the transmitter. Therefore, large distortions can be expected if near‐surface polarizable layers exist. This was verified in practice by field measurements in an area with high chargeabilities in Longerich, Cologne.  相似文献   

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.
The nth-order moments of the electromagnetic impulse response are useful for interpreting electromagnetic data. We have derived an analytic expression for the half-order moment of a conductive half-space. By inverting this expression, the measured half-order moment can be used to estimate an apparent conductivity of the ground. The first-order moment can also be used to estimate the half-space conductivity. A sensitivity analysis indicates that for an airborne EM configuration, the half-order moment will be most sensitive to material in the top 26–48 m, while the first-order moment will be sensitive to deeper material (down to depths between 66 and 127 m).  相似文献   

9.
利用三维时域有限差分技术（FDTD）研究建立倾斜各向异性地层中多分量电磁测井响应的数值模拟算法.首先,利用Yee氏交错差分格式和线性内插公式对各向异性介质中时域Maxwell方程进行离散化,得到各向异性介质中各个交错节点上电磁场的时域递推公式,同时将面积加权平均与旋转矩阵技术相结合确定各个网格节点的等效电导率张量,此外,单轴各向异性吸收边界条件（UPML）被用于降低求解区域外边界反射.在此基础上,利用单频正弦磁偶源计算时域电磁场空间分布以及特殊的双方程算法（2E2U）从数值结果中快速提取电磁场幅度和相位,确定多分量电磁测井仪器的三维响应.最后,通过数值结果证明该算法的有效性,并系统考察不同层厚、地层倾角以及侵入带半径等变化对多分量电磁测井响应的影响.  相似文献   

10.
We have derived asymptotic expansions that can be used to estimate the impulse response of a conductive halfspace when excited by an airborne electromagnetic system. Two expressions are required to calculate the response, one is applicable in the late-time regime and a second is used for the early-time regime. Fortunately, there is a small overlap in the late- and early-time regimes. Asymptotic expansions are commonly used for interpreting ground EM data, but our expressions are the first ones applicable to the airborne situation. The expansions are simple to program, can be calculated extremely rapidly and will be useful for understanding and interpreting airborne electromagnetic data. As an example of the use of the expressions, we have calculated the response of seawater and compared this with the airborne electromagnetic response measured off the coast of Perth, Western Australia.  相似文献   

11.
A large closed wire loop is generally used in field experiments for testing airborne electrical exploration equipment. Thus, methods are required for the precise calculation of an electromagnetic response in the presence of a closed wire loop. We develop a fast and precise scheme for calculating the transient response for such a closed loop laid out at the surface of a horizontally layered conductive ground. Our scheme is based on the relationship between the magnetic flux flowing through a closed loop and the current induced in it. The developed scheme is compared with 2D and 3D finite‐element modelling for several positions of an airborne electromagnetic system flying over a closed loop. We also study the coupling effect between the current flowing in the closed loop and the current flowing in the horizontally layered conductive medium. The result shows that for the central position of the transmitter, the difference between axisymmetrical finite‐element modelling and our scheme is less than 1%. Moreover, for the non‐coaxial transmitter–receiver–loop system, the solution obtained by our scheme is in good agreement with full 3D finite‐element modelling, and our total simulation time is substantially lower: 1 minute versus 120 hours.  相似文献   

12.
To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss–Newton optimization. Forward modeling is performed in the frequency domain based on the scattered secondary electrical field. Then, the inverse Fourier transform and convolution of the transmitting waveform are used to calculate the EM responses and the sensitivity matrix in the time domain for arbitrary transmitting waves. To optimize the computational time and memory requirements, we use the EM “footprint” concept to reduce the model size and obtain the sparse sensitivity matrix. To improve the 3D inversion, we use the OpenMP library and parallel computing. We test the proposed 3D parallel inversion code using two synthetic datasets and a field dataset. The time-domain airborne EM inversion results suggest that the proposed algorithm is effective, efficient, and practical.  相似文献   

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

14.
Electromagnetic fields in a non-uniform steel-cased borehole   总被引：1，自引：0，他引：1  
Since most oil wells are cased in steel, electromagnetic (EM) signals undergo severe attenuation as they diffuse across the casing. This paper examines an effect of non‐uniform casing properties on EM fields measured in a steel‐cased well embedded in a layered formation. We use a finite‐element method for computing secondary azimuthal electric fields in a cylindrically symmetric model, and analytically obtain primary fields for a homogeneous casing in a homogeneous whole space. Although steel casing largely masks EM signals induced into a layered formation, the induced signal is more pronounced in phase than in amplitude. The effect of casing non‐uniformity is quite large in measured fields but is highly localized. When electrical conductivity varies rapidly in the casing wall, the resulting EM fields also vary rapidly. A cross‐correlation function of these variations has strong peaks at two points, the interval between them being equal to the source–receiver distance. The high‐frequency coherent noise event caused by the non‐uniform casing can be greatly suppressed by low‐pass filtering to enhance EM signals indicating formation conductivity.  相似文献   

15.
万伟  唐新功  黄清华 《地球物理学报》2019,62(12):4846-4859
陆地可控源电磁法的观测资料可依据频段范围近似地划分为近区场、中间区场及远区场,但采用测量相互正交电、磁分量,并计算视电阻率的资料处理方式只适用于远区场数据.为更有效地利用陆地可控源电磁法不同区间场的观测资料,本文结合三维数值模拟技术并采用电场分量直接进行反演的策略,对不同区间电场的响应特征与探测效果进行了分析.数值模拟结果表明:近区电场的异常响应最明显,异常响应不随频率发生显著变化,但纵向分辨能力差;远区电场异常响应随频率发生显著变化,其探测深度取决于频率的高低;中间区场较为复杂,地表电场异常响应的等值线中心并不是位于异常体中心正上方,而是在沿场传播方向上向异常体与围岩的分界面处偏移,并且发现中间区场资料的加入会影响反演结果的准确性.综合合成数据和野外实测资料的反演结果,发现结合近区场和远区场资料而舍弃中间区场资料的反演效果更佳,这为陆地可控源电磁法资料的反演解释提供了一种有效途径.  相似文献   

16.
在直升机航空电磁法中,吊舱的姿态变化（如摆动、倾斜、偏航等）会使磁偶极子相对正常飞行时在方向和位置上发生变化,使得接收线圈处的电磁响应产生误差.传统的基于重叠偶极模型的姿态校正方法只考虑了磁偶极子方向上的变化,而忽略了磁偶极子位置上的变化,因此校正效果受到了一定的限制.本文分析了五种不同电性参数的均匀半空间模型中线圈系在姿态变化前后的电磁响应比,进而提出了一种新的姿态校正方法,该方法综合考虑了磁偶极子方向和位置变化对姿态变化前后的电磁响应比的影响,并推导出总姿态变化电磁响应比与方向和位置单独变化时电磁响应比的乘积近似成固定比例关系.为了验证本文提出的总姿态校正方法,本文分别使用本方案及基于重叠偶极模型的姿态校正方法对五组具有代表性的合成数据进行了校正.对比结果表明,本文的总姿态校正方法较之重叠偶极姿态校正方法能较大地提高校正精度.  相似文献   

17.
The relatively new subdiscipline of environmental geophysics has grown enormously in the last five years. The size and diversity of the field, and the associated literature, is such that it is extremely difficult to keep up with even a small portion of the field. Electrical and electromagnetic (E & EM) methods, including ground penetrating radar and time-domain reflectometry, play a central role in environmental geophysics. One reason for the utility of E & EM methods in groundwater studies is the similarity in the way that current flow and fluid flow depend on the connectivity and geometry of the pores in soils and rocks. Another reason is the influence of the pore water quality on the geophysical response. More than any other geophysical technique, E & EM methods are directly affected by the presence of conductive pore fluids in the subsurface, such as leachates from landfill sites and sea water invading a coastal groundwater supply that has been placed under stress because of population expansion. The chloride ion is one of the most electrically active of the naturally-occurring ions, and allows us to detect sea water incursion; leachates from landfill sites contain the by-products of organic decay, such as acetic acid, which are generally less conductive than chloride, but nonetheless enhance the pore water and formation electrical conductivities. Landfill leachate plumes are thus easily mapped. The shallow subsurface electrical and dielectric properties exhibit hysteresis due to seasonal changes in water content; the physical properties will be different for the same degree of saturation, depending on whether the water level is rising or falling. Topographic effects are also important; an empirical correction method works well to remove a background trend in the conductivity due to changes in elevation. Heterogeneity and anisotropy of the electric properties may be related to similar effects in the hydraulic properties. New technology and the adaptation of existing technology has lead to the development of fresh instruments, such as electrode arrays towed across the ground, resistivity logging while drilling, fast-rise time TEM, NMR combined with TEM, electric quadripole, et cetera. The applications of E & EM methods cover a wide range of geographic areas and groundwater problems, but have had particularly wide use for groundwater exploration in arid and semi-arid regions, for mapping and monitoring salt-water incursion in susceptible aquifers, and for mapping and monitoring contaminants.  相似文献   

18.
Very low frequency (VLF) military communications systems provide a primary field that can be used for shallow geophysical surveys to locate ground water contamination and vertical geologic contacts. Useful properties that can be easily obtained from the interaction of the earth and the primary field are the magnitude of the vertical secondary magnetic field, the surface impedence, and the phase angle between the electrical and magnetic horizontal components. The variations in the secondary magnetic field can be related to vertical geologic contacts, such as the edges of landfill trenches. The surface impedence yields an apparent terrain conductivity, which can be used to locate low-resistivity anomalies often associated with contaminated ground water. The phase angle gives information on vertical variations in resistivity, phase angles less than 45° indicating increasing resistivity with depth. The depth of penetration of the VLF field is about one skin depth. For a frequency of 20 kHz, the skin depth in meters is approximately equal to 3.67 where p is terrain resistivity in ohmmeters.  相似文献   

19.
Abstract

In a nonrotating system, the shear Reynolds stresses exerted by surface or internal gravity waves vanish on account of the exact quadrature between the horizontal and vertical orbital velocities. It is shown that a rotation of the system induces small in‐phase perturbations, resulting in a mean Reynolds stress which can generate low frequency currents. If both the wave field and the ocean are homogeneous with respect to the horizontal coordinates, the low‐frequency response is an undamped inertial oscillation. If either the wave field or the ocean are weakly inhomogeneous, the oscillation disperses in the vertical and horizontal directions due to phase‐mixing of modes with closely neighboring frequencies. Other effects which produce small frequency shifts also contribute to phase‐mixing, for example the horizontal component of the Coriolis vector and nonlinear interactions with geo‐strophic currents. The analysis is based on operator representations which avoid normal mode decomposition and yield simple integro‐differential operators for each phase‐mixing process. Numerical results are presented for a continuously stratified model typical for a shallow sea (Baltic). The orders of magnitude and qualitative features are in reasonable agreement with observations.  相似文献   

20.
A new approach to self-potential (SP) data interpretation for the recognition of a buried causative SP source system is presented. The general model considered is characterized by the presence of primary electric sources or sinks, located within any complex resistivity structure with a flat air-earth boundary. First, using physical considerations of the nature of the electric potential generated by any arbitrary distribution of primary source charges and the related secondary induced charges over the buried resistivity discontinuity planes, a general formula is derived for the potential and the electric field component along any fixed direction on the ground surface. The total effect is written as a sum of elementary contributions, all of the same simple mathematical form. It is then demonstrated that the total electric power associated with the standing natural electric field component can be written in the space domain as a sum of cross-correlation integrals between the observed component of the total electric field and the component of the field due to each single constitutive elementary charge. By means of the cross-correlation bounding inequality, the concept of a scanning function is introduced as the key to the new interpretation procedure. In the space domain, the scanning function is the unit strength electric field component generated by an elementary positive charge. Next, the concept of charge occurrence probability is introduced as a suitable function for the tomographic imaging of the charge distribution geometry underground. This function is defined as the cross-correlation product of the total observed electric field component and the scanning function, divided by the square root of the product of the respective variances. Using this physical scheme, the tomographic procedure is described. It consists of scanning the section, through any SP survey profile, by the unit strength elementary charge, which is given a regular grid of space coordinates within the section, at each point of which the charge occurrence probability function is calculated. The complete set of calculated grid values can be used to draw contour lines in order to single out the zones of highest probability of concentrations of polarized, primary and secondary electric charges. An extension to the wavenumber domain and to three-dimensional tomography is also presented and discussed. A few simple synthetic examples are given to demonstrate the resolution power of the new SP inversion procedure.  相似文献   

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

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