Signal detectability of marine electromagnetic methods in the exploration of resistive targets

We compare selected marine electromagnetic methods for sensitivity to the presence of relatively thin resistive targets (e.g., hydrocarbons, gas hydrates, fresh groundwater, etc.). The study includes the conventional controlled‐source electromagnetic method, the recently introduced transient electromagnetic prospecting with vertical electric lines method, and the novel marine circular electric dipole method, which is still in the stage of theoretical development. The comparison is based on general physical considerations, analytical (mainly asymptotic) analysis, and rigorous one‐dimensional and multidimensional forward modelling. It is shown that transient electromagnetic prospecting with vertical electric lines and marine circular electric dipole methods represent an alternative to the conventional controlled‐source electromagnetic method at shallow sea, where the latter becomes less efficient due to the air‐wave phenomenon. Since both former methods are essentially short‐offset time‐domain techniques, they exhibit a much better lateral resolution than the controlled‐source electromagnetic method in both shallow sea and deep sea. The greatest shortcoming of the transient electromagnetic prospecting with vertical electric lines and marine circular electric dipole methods comes from the difficulties in accurately assembling the transmitter antenna within the marine environment. This makes these methods significantly less practical than the controlled‐source electromagnetic method. Consequently, the controlled‐source electromagnetic method remains the leading marine electromagnetic technique in the exploration of large resistive targets in deep sea. However, exploring laterally small targets in deep sea and both small and large targets in shallow sea might require the use of the less practical transient electromagnetic prospecting with vertical electric lines and/or marine circular electric dipole method as a desirable alternative to the controlled‐source electromagnetic method.     

Time domain electromagnetic‐induced polarisation: extracting more induced polarisation information from grounded source time domain electromagnetic data

Electrical induced polarisation surveys are used to detect chargeable materials in the earth. For interpretation of time domain electrical‐induced polarisation data a common procedure is to first invert the direct current data (electric current on time) to recover conductivity and then invert the induced polarisation data (current off‐time) to recover chargeability. This direct current‐induced polarisation inversion procedure assumes that the off time data are free of secondary electromagnetic induction effects. To comply with this, early time data are often discarded or not recorded. For mid‐time data, an electromagnetic decoupling technique, which removes electromagnetic induction in the observations, needs to be implemented. Usually, responses from a half‐space or a layered earth are subtracted. Recent capability in three‐dimensional time domain electromagnetic forward modelling and inversion allows to revisit these procedures. In a Time domain electromagnetic‐induced polarisation survey, a high sampling rate allows early time channels of the electromagnetic data to be recorded. The recovery of chargeability then follows a three‐step workflow: (i) invert early time channel time domain electromagnetic data to recover the three‐dimensional conductivity; (ii) use that conductivity to compute the time domain electromagnetic response at later time channels and subtract this fundamental response from the observations to extract the induced polarisation responses, and (iii) invert the induced polarisation responses to recover a three‐dimensional chargeability. This workflow effectively removes electromagnetic induction effects in the observations and produces better chargeability and conductivity models compared with conventional approaches. In a synthetic example involving a gradient array, we show that the conductivity structure obtained from the early time channel data, which are usually discarded, is superior to that obtained from the steady state direct current voltages. This adds a further reason to collect these electromagnetic data.     

Two types of marine controlled source electromagnetic transmitters

Marine controlled source electromagnetic methods are used to derive the electrical properties of a wide range of sub‐seafloor targets, including gas hydrate reservoirs. In most marine controlled source electromagnetic surveys, the deep‐tow transmitter is used with a long horizontal electric dipole being towed above the seafloor, which is capable of transmitting dipole moments in the order of up to several thousand ampere‐metres. The newly developed deployed transmitter uses two horizontal orthogonal electrical dipoles and can land on the seafloor. It can transmit higher frequency electromagnetic signals, can provide accurate transmission orientation, and can obtain higher signal stacking, which compensates for the shorter source dipole length. In this paper, we present the study, key technologies, and implementation details of two new marine controlled source electromagnetic transmitters (the deep‐tow transmitter and the deployed transmitter). We also present the results of a marine controlled source electromagnetic experiment conducted from April to May 2014 in the South China Sea using both the deep‐tow transmitter and the deployed transmitter, which show that the two types of marine transmitters can be used as effective source for gas hydrate exploration.     

Electromagnetic fields generated by finite‐length wire sources: comparison with point dipole solutions

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.     

Elimination of the water-layer response from multi-component source and receiver marine electromagnetic data

This paper presents the theory to eliminate from the recorded multi‐component source, multi‐component receiver marine electromagnetic measurements the effect of the physical source radiation pattern and the scattering response of the water‐layer. The multi‐component sources are assumed to be orthogonally aligned above the receivers at the seabottom. Other than the position of the sources, no source characteristics are required. The integral equation method, which for short is denoted by Lorentz water‐layer elimination, follows from Lorentz' reciprocity theorem. It requires information only of the electromagnetic parameters at the receiver level to decompose the electromagnetic measurements into upgoing and downgoing constituents. Lorentz water‐layer elimination replaces the water layer with a homogeneous half‐space with properties equal to those of the sea‐bed. The source is redatumed to the receiver depth. When the subsurface is arbitrary anisotropic but horizontally layered, the Lorentz water‐layer elimination scheme greatly simplifies and can be implemented as deterministic multi‐component source, multi‐component receiver multidimensional deconvolution of common source gathers. The Lorentz deconvolved data can be further decomposed into scattering responses that would be recorded from idealized transverse electric and transverse magnetic mode sources and receivers. This combined electromagnetic field decomposition on the source and receiver side gives data equivalent to data from a hypothetical survey with the water‐layer absent, with idealized single component transverse electric and transverse magnetic mode sources and idealized single component transverse electric and transverse magnetic mode receivers. When the subsurface is isotropic or transverse isotropic and horizontally layered, the Lorentz deconvolution decouples into pure transverse electric and transverse magnetic mode data processing problems, where a scalar field formulation of the multidimensional Lorentz deconvolution is sufficient. In this case single‐component source data are sufficient to eliminate the water‐layer effect. We demonstrate the Lorentz deconvolution by using numerically modeled data over a simple isotropic layered model illustrating controlled‐source electromagnetic hydrocarbon exploration. In shallow water there is a decrease in controlled‐source electromagnetic sensitivity to thin resistors at depth. The Lorentz deconvolution scheme is designed to overcome this effect by eliminating the water‐layer scattering, including the field's interaction with air.     

Generalized Tikhonov's algorithm for accurate calculation of one-dimensional transient responses directly in time domain

In the last two decades, forward modelling for the time domain (transient) electromagnetic method has concentrated almost entirely on multi-dimensional models and algorithms. At the same time, the interpretation of real field data is still mainly one dimensional. This is caused by the lack of an efficient multi-dimensional acquisition procedure supported by sufficiently fast and reliable inversion software, on the one hand, and by the great efficiency of one-dimensional field set up and interpretation of the data on the other hand. The latter is particularly true for the short offset transient electromagnetic method, which is much less sensitive to multi-dimensional effects, compared to long offset methods. The most commonly used one-dimensional forward modelling algorithms are based on the spectral method, which requires calculating rapidly oscillating Fourier–Bessel (Hankel) integrals. Due to the very fast decay of short offset responses, the integrals become computationally unstable at late times of the transient process. Although this problem has been successfully solved for practically feasible measurement times of conventional short offset systems using transverse electric and mixed transverse electric and transverse magnetic fields, it turned out crucial for novel methods based on the use of unimodal transverse magnetic fields. These methods are much more sensitive to geoelectric parameters of the Earth in general and those of resistive targets, in particular, but they generate responses, which drop at late times significantly faster than those of conventional methods. Such behaviour of transverse magnetic fields represents severe computational problem for the spectral method, but is successfully solved by direct time domain algorithms. This article describes a generalization of the well-known Tikhonov's solution to a boundary value problem directly in time domain, which is applied to an arbitrary one-dimensional earth model excited by an arbitrary source. Contrary to existing spectral algorithms, the described method allows accurate calculations of both transverse electric and transverse magnetic transient responses at arbitrarily late times. On the other hand, it is more time efficient than finite-difference/finite element direct time domain algorithms and provides analytical late-stage asymptotic solutions.     

3D magnetotelluric modelling including surface topography

An edge finite‐element method has been applied to compute magnetotelluric (MT) responses to three‐dimensional (3D) earth topography. The finite‐element algorithm uses a single edge shape function at each edge of hexahedral elements, guaranteeing the continuity of the tangential electric field while conserving the continuity of magnetic flux at boundaries. We solve the resulting system of equations using the biconjugate gradient method with a Jacobian preconditioner. The solution gives electric fields parallel to the slope of a surface relief that is often encountered in MT surveys. The algorithm is successfully verified by comparison with other numerical solutions for a 3D‐2 model for comparison of modelling methods for EM induction and a ridge model. We use a 3D trapezoidal‐hill model to investigate 3D topographic effects, which are caused mainly by galvanic effects, not only in the Zxy mode but also in the Zyx mode. If a 3D topography were approximated by a two‐dimensional topography therefore errors occurring in the transverse electric mode would be more serious than those in the transverse magnetic mode.     

Effects of the sea floor topography on the 1D inversion of time‐domain marine controlled source electromagnetic data

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.     

海底电性源频率域CSEM勘探建模及水深影响分析

为了探索我国海域油气和水合物等高阻目标体CSEM勘探的可行性和方法技术,本文研究了在海水中水平电性源激励下有限水深海洋地电模型的频率域电磁响应,为进一步的1D和3D仿真计算奠定了理论基础.在推导电磁响应公式时,首先给出了各层介质的Lorentz势,然后根据Coulomb势与Lorentz势的关系,得到了各层介质的Coulomb势.各层介质中的电磁场均可以由Lorentz势或者Coulomb势计算得到,但在有限元计算时Coulomb势具有优势.长导线源的电磁场和势函数可以由电偶源的电磁场和势函数沿导线长度积分得到.文中具体给出了海水中水平电偶源和长导线源在海水层的电磁场公式,并根据该公式计算了不同水深环境下海底表面的电磁场分布,分析了海水深度对海底油气储层电磁异常的影响.结果表明,随着水深减小,异常幅度和形态特征发生明显变化.当水深很浅时(如50 m),只有同线方向的E_x和E_z两个电场分量存在明显异常.最后,以两个已知海底油田为例,计算了不同水深环境下可观测到的电场异常,展示了电性源频率域CSEM在海底勘探中(包括浅海环境)的良好应用前景.对于该方法实用化过程中还需进一步解决的问题,文中结尾部分也进行了初步探讨.     

The electromagnetic response of a horizontal electric dipole buried in a multi‐layered earth

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.     

Navigating marine electromagnetic transmitters using dipole field geometry

The marine controlled source electromagnetic (CSEM) technique has been adopted by the hydrocarbon industry to characterize the resistivity of targets identified from seismic data prior to drilling. Over the years, marine controlled source electromagnetic has matured to the point that four‐dimensional or time lapse surveys and monitoring could be applied to hydrocarbon reservoirs in production, or to monitor the sequestration of carbon dioxide. Marine controlled source electromagnetic surveys have also been used to target shallow resistors such as gas hydrates. These novel uses of the technique require very well constrained transmitter and receiver geometry in order to make meaningful and accurate geologic interpretations of the data. Current navigation in marine controlled source electromagnetic surveys utilize a long base line, or a short base line, acoustic navigation system to locate the transmitter and seafloor receivers. If these systems fail, then rudimentary navigation is possible by assuming the transmitter follows in the ship's track. However, these navigational assumptions are insufficient to capture the detailed orientation and position of the transmitter required for both shallow targets and repeat surveys. In circumstances when acoustic navigation systems fail we propose the use of an inversion algorithm that solves for transmitter geometry. This algorithm utilizes the transmitter's electromagnetic dipole radiation pattern as recorded by stationary, close range (<1000 m), receivers in order to model the geometry of the transmitter. We test the code with a synthetic model and validate it with data from a well navigated controlled source electromagnetic survey over the Scarborough gas field in Australia.     

On small-scale near-surface distortion in controlled source tensor electromagnetics1

Based on a Born approximation of a thin sheet integral equation, it is shown that small-scale surficial heterogeneity significantly distorts the electromagnetic field excited by electric dipoles only when either the source or the receiver are located on the heterogeneities. When a surface heterogeneity is beneath the source, the associated distortion of the electromagnetic field is manifest as a change in the effective electric dipole moment. Hence the magnetotelluric transfer functions and impedance relations remain undistorted in this case. When a surface heterogeneity is beneath the receiver, the electric field is severely distorted, but the magnetic field is only slightly distorted. The impedance tensor is therefore strongly distorted, but the tipper vector is almost unaltered. Since the controlled source tipper is a function of 1D earth conductivity, it is proposed that tipper data should be used in the first stage of 1D interpretation. For a 1D earth, the tipper vector must always point towards the source and, in the near-field limit, should have unit length. These two necessary conditions must be met by the measured tipper before it is interpreted one dimensionally.     

Field test of sub‐basalt hydrocarbon exploration with marine controlled source electromagnetic and magnetotelluric data

The recent use of marine electromagnetic technology for exploration geophysics has primarily focused on applying the controlled source electromagnetic method for hydrocarbon mapping. However, this technology also has potential for structural mapping applications, particularly when the relative higher frequency controlled source electromagnetic data are combined with the lower frequencies of naturally occurring magnetotelluric data. This paper reports on an extensive test using data from 84 marine controlled source electromagnetic and magnetotelluric stations for imaging volcanic sections and underlying sediments on a 128‐km‐long profile. The profile extends across the trough between the Faroe and Shetland Islands in the North Sea. Here, we focus on how 2.5D inversion can best recover the volcanic and sedimentary sections. A synthetic test carried out with 3D anisotropic model responses shows that vertically transverse isotropy 2.5D inversion using controlled source electromagnetic and magnetotelluric data provides the most accurate prediction of the resistivity in both volcanic and sedimentary sections. We find the 2.5D inversion works well despite moderate 3D structure in the synthetic model. Triaxial inversion using the combination of controlled source electromagnetic and magnetotelluric data provided a constant resistivity contour that most closely matched the true base of the volcanic flows. For the field survey data, triaxial inversion of controlled source electromagnetic and magnetotelluric data provides the best overall tie to well logs with vertically transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data a close second. Vertical transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data provided the best interpreted base of the volcanic horizon when compared with our best seismic interpretation. The structural boundaries estimated by the 20‐Ω·m contour of the vertical resistivity obtained by vertical transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data gives a maximum geometric location error of 11% with a mean error of 1.2% compared with the interpreted base of the volcanic horizon. Both the model study and field data interpretation indicate that marine electromagnetic technology has the potential to discriminate between low‐resistivity prospective siliciclastic sediments and higher resistivity non‐prospective volcaniclastic sediments beneath the volcanic section.     

基于改进的接收点插值算法的频率域海洋可控源电磁法2.5维正演

在海洋可控源电磁法勘探中,接收站常置于海底.在进行海洋电磁场模拟时,由于海水和海底介质存在显著电性差异,这给海底接收点处场值的求取带来困难.本文提出一种新的接收点插值算法,该算法考虑到海底电场法向分量不连续性问题,用法向电流分量进行插值以准确求取海底任意接收点处电磁场值.本文利用交错网格有限差分法实现了二维介质中频率域海洋可控源法(CSEM)正演.对构造走向做傅里叶变换,将三维电磁模拟问题转换为波数域2.5维问题,即三维场源激励下针对二维地电模型的电磁模拟问题.使用交错网格有限差分法,基于一次场/二次场分离方法导出波数域二次电场离散形式,并进一步求得波数域电磁场.采用本文提出的改进的插值算法可求得海底任意接收点处波数域电磁场,采用傅里叶逆变换对波数域电磁场进行积分可得到接收点处空间域电磁场.模型算例表明,与常规的线性插值和严格插值算法相比,本文提出的改进的插值算法具有更高的精度.     

多源频率域地空系统三维电磁响应分析

地空电磁法已经成为深部资源勘探的重要地球物理方法,但对频率域地空系统的三维多源电磁响应特征研究较少.本文设计了多种激励源组合方式,采用非结构化有限元数值模拟方法,对三维地电模型的空中垂直磁场的响应特征进行了研究.首先推导了基于电场的双旋度公式及其变分形式,加入罚项以减少伪解的影响.接着把有限元稀疏矩阵方程转换为频率的函数,采用Krylov子空间投影方法,通过模型降阶算法降低稀疏矩阵的阶数,实现多频点的快速计算.建立了三维低阻体模型、高阻体模型以及两个相邻低阻体模型,分别采用单源、双源、三源和四源激励模式,从垂直磁场的总场、二次场响应和全域视电阻率等方面进行分析比较.结果表明:多源地空电磁法不仅可以增加总场的强度,而且可以改变异常体的二次电磁响应分布规律.各电偶源延长线呈正三角形分布的三源和矩形分布的四源激励模式在增强信号强度以及削弱异常体的边界效应方面具有一定的优势,是一种优化的多源激励方式.     

复杂场源形态的海洋可控源电磁三维正演

在使用电偶极发射源的可控源电磁法(CSEM)勘探中,发射源的方位、长度、形状等对观测数据有重要的影响,然而现有的大部分三维数值模拟方法没有全面地将这些因素考虑进来,很多都只能应对非常简单的场源形态,例如单一方位的点电偶极子,这有可能显著降低模拟结果的准确性.本文实现了基于交错网格有限体积(FV)离散的海洋CSEM三维正演算法,能够模拟形态相对复杂的场源,包括任意方位的有限长直导线和弯曲导线发射源.该算法使用一次场/二次场方法,只需对二次场使用FV法求解,避免了场源的奇异性问题;一次场的计算为一维正演问题,使用准解析法求解,并且只要在计算一次场时考虑复杂的场源形态便可以实现同样场源的三维正演.通过与一维理论模型的解析解对比验证了三维程序的准确性,并针对三维理论模型进行了一系列正演测试,初步考察了场源形态对三维正演结果的影响.     

一维电阻率各向异性对海洋可控源电磁响应的影响研究

地下介质的电阻率常常表现为各向异性,海底裂隙地层和层状沉积序列可能形成宏观电阻率各向异性.在解释海洋电磁资料时,电阻率各向异性的影响不应该被忽略,否则可能会得到错误的海底地电模型.作者编写了电阻率任意各向异性一维层状介质海洋可控源电磁场计算程序,计算了电阻率各向异性层状模型的海洋可控源电磁响应,讨论了覆盖层和高阻储层分别具有电阻率各向异性时的电磁场响应特征.     

Sensitivity and inversion of marine electromagnetic data in a vertically anisotropic stratified earth

The controlled‐source electromagnetic (CSEM) and magnetotelluric method (MT) are two techniques that can be jointly used to explore the resistivity structure of the earth. Such methods have, in recent years, been applied in marine environments to the exploration and appraisal of hydrocarbons. In many situations the electric properties of the earth are anisotropic, with differences between resistivity in the vertical direction typically much higher than those in the horizontal direction. In cases such as this, the two modes of the time‐harmonic electromagnetic field are altered in different ways, implying that the sensitivity to the earth resistivity may vary significantly from one particular resistivity component (scalar, horizontal or vertical) to another, depending on the measurement configuration (range, azimuth, frequency or water depth). In this paper, we examine the sensitivity of the electromagnetic field to a vertically anisotropic earth for a typical set of configurations, compare inversion results of synthetic data characterizing a vertically anisotropic earth obtained using the isotropic and anisotropic assumptions and show that correctly accounting for anisotropy can prevent artefacts in inversion results.     

海底油气藏地质模型的冲激响应

海洋可控源电磁法(mCSEM)的时间域冲激响应特征可以反映海底油气高阻薄层.本文计算了水平电偶极子源均匀大地半空间,海洋均匀双半空间和海洋四层模型的阶跃响应和冲激响应,提出了瞬变冲激时刻的概念.分析了水平电偶源瞬变冲激时刻与介质电导率的指示关系.对于海底油气高阻薄层宜采用多偏移距同时测量方式,由于在低电导率介质中电磁能量传播得要快,在适当的收发距瞬变冲激时刻会提前到达,提出的瞬变冲激时刻道间变化量可以明确指示高阻薄层的存在及埋深.文中还分析了海水深度对瞬变冲激时刻的影响.由于“天波”干扰,瞬变冲激响应受到一定收发距观测的限制.消除 “天波”影响是时间域和频率域mCSEM数据处理的研究热点.     

利用数值积分提高一维模型电偶源电磁测深响应计算精度

讨论了直接利用数值积分提高电偶源电磁测深响应计算精度的方法．具体为对Hankel积分进行直接积分，结合连分式展开方法以提高积分求和的收敛速度．利用该方法对均匀半空间和层状(两层)模型的电磁测深响应进行了模拟，结果表明与常规的快速Hankel滤波方法相比，采用直接数值积分能明显提高电偶源频率测深响应计算精度．从而为获得高精度瞬变测深晚期响应提供算法基础．