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.     

海底可控源电磁接收机及其水合物勘查应用

海洋可控源电磁法在国外已成为海底天然气水合物调查的有效手段之一.为实现我国海域深水条件下水合物的海洋可控源电磁探测,本文从方法原理出发,采用低功耗嵌入式控制、前端低噪声斩波放大、高精度时间同步和水声通讯等技术,设计并开发了由承压舱、玻璃浮球、采集电路、电场与磁场传感器、姿态测量装置、声学释放器、USBL定位信标、测量臂、水泥块等部件组成的海底可控源电磁接收机,实现了海洋微弱电磁场信号的高精度采集.海底可控源电磁接收机具有高可靠性、低噪声、低功耗和低时漂的特点.利用研制的海底可控源电磁接收机,在琼东南海域进行水合物勘查,采集得到了可靠的人工源电磁场数据.通过数据处理及反演,获得了研究区海底的电阻率模型,结合地震资料,对高阻异常体进行推断解释,其结果为天然气水合物钻探井位布置提供了电性依据.     

Trans‐dimensional Bayesian inversion of controlled‐source electromagnetic data in the German North Sea

This paper presents the first controlled‐source electromagnetic survey carried out in the German North Sea with a recently developed seafloor‐towed electrical dipole–dipole system, i.e., HYDRA II. Controlled‐source electromagnetic data are measured, processed, and inverted in the time domain to estimate an electrical resistivity model of the sub‐seafloor. The controlled‐source electromagnetic survey targeted a shallow, phase‐reversed, seismic reflector, which potentially indicates free gas. To compare the resistivity model to reflection seismic data and draw a combined interpretation, we apply a trans‐dimensional Bayesian inversion that estimates model parameters and uncertainties, and samples probabilistically over the number of layers of the resistivity model. The controlled‐source electromagnetic data errors show time‐varying correlations, and we therefore apply a non‐Toeplitz data covariance matrix in the inversion that is estimated from residual analysis. The geological interpretation drawn from controlled‐source electromagnetic inversion results and borehole and reflection seismic data yield resistivities of ～1 Ωm at the seafloor, which are typical for fine‐grained marine deposits, whereas resistivities below ～20 mbsf increase to 2–4 Ωm and can be related to a transition from fine‐grained (Holocene age) to unsorted, coarse‐grained, and compacted glacial sediments (Pleistocene age). Interface depths from controlled‐source electromagnetic inversion generally match the seismic reflector related to the contrast between the different depositional environments. Resistivities decrease again at greater depths to ～1 Ωm with a minimum resistivity at ～300 mbsf where a seismic reflector (that marks a major flooding surface of late Miocene age) correlates with an increased gamma‐ray count, indicating an increased amount of fine‐grained sediments. We suggest that the grain size may have a major impact on the electrical resistivity of the sediment with lower resistivities for fine‐grained sediments. Concerning the phase‐reversed seismic reflector that was targeted by the survey, controlled‐source electromagnetic inversion results yield no indication for free gas below it as resistivities are generally elevated above the reflector. We suggest that the elevated resistivities are caused by an overall decrease in porosity in the glacial sediments and that the seismic reflector could be caused by an impedance contrast at a thin low‐velocity layer. Controlled‐source electromagnetic interface depths near the reflector are quite uncertain and variable. We conclude that the seismic interface cannot be resolved with the controlled‐source electromagnetic data, but the thickness of the corresponding resistive layer follows the trend of the reflector that is inclined towards the west.     

A land controlled‐source electromagnetic experiment using a deep vertical electric dipole: experimental settings,processing, and first data interpretation

A multichannel borehole‐to‐surface controlled‐source electromagnetic experiment was carried out at the onshore CO₂ storage site of Hontomín (Spain). The electromagnetic source consisted of a vertical electric dipole located 1.5 km deep, and the electric field was measured at the surface. The subsurface response has been obtained by calculating the transfer function between the transmitted signal and the electric field at the receiver positions. The dataset has been processed using a fast processing methodology, appropriate to be applied on controlled‐source electromagnetics (CSEM) data with a large signal‐to‐noise ratio. The dataset has been analysed in terms of data quality and repeatability errors, showing data with low experimental errors and good repeatability. We evaluate if the induction of current along the casing of the injection well can reproduce the behaviour of the experimental data.     

Feasibility of signal enhancement with multiple grounded‐wire sources for a frequency‐domain electromagnetic survey

Frequency‐domain electromagnetic methods with a grounded‐wire source are powerful tools in geophysical exploration. However, the signal may be too weak to guarantee the quality of survey data in complex electromagnetic environments, especially when the receiver is located in the air for the newly developed grounded‐source airborne frequency‐domain electromagnetic method. In this paper, a signal enhancement method with multiple sources is proposed to solve this problem. To evaluate the signal enhancement effect, we compared the signals generated by a single source and multiple sources with equal electric moment. The signal differences caused by synchronisation error and separation distance between source elements were analysed, and the methods to achieve maximum signal were introduced. Besides, we discussed the interaction between adjacent source elements to ensure the system safety, including the changes in output current and the safe distance between two sources using a dual‐source model. Lastly, a comprehensive field experiment was designed and conducted to test the multiple‐source method. The data processing results are comparable for single and dual sources, and the signal‐to‐noise ratio of dual source is higher in the field test. The subsurface resistivity structure at the test site is consistent with the previous controlled‐source audio‐frequency magnetotellurics method. These results show that signal enhancement with multiple sources is feasible. This study provides guidance to the application of multiple sources in field surveys when the survey environment is complex and rigorous.     

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.     

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.     

1D joint multi‐offset inversion of time‐domain marine controlled source electromagnetic data

The accurate estimation of sub‐seafloor resistivity features from marine controlled source electromagnetic data using inverse modelling is hindered due to the limitations of the inversion routines. The most commonly used one‐dimensional inversion techniques for resolving subsurface resistivity structures are gradient‐based methods, namely Occam and Marquardt. The first approach relies on the smoothness of the model and is recommended when there are no sharp resistivity boundaries. The Marquardt routine is relevant for many electromagnetic applications with sharp resistivity contrasts but subject to the appropriate choice of a starting model. In this paper, we explore the ability of different 1D inversion schemes to derive sub‐seafloor resistivity structures from time domain marine controlled source electromagnetic data measured along an 8‐km‐long profile in the German North Sea. Seismic reflection data reveal a dipping shallow amplitude anomaly that was the target of the controleld source electromagnetic survey. We tested four inversion schemes to find suitable starting models for the final Marquardt inversion. In this respect, as a first scenario, Occam inversion results are considered a starting model for the subsequent Marquardt inversion (Occam–Marquardt). As a second scenario, we employ a global method called Differential Evolution Adaptive Metropolis and sequentially incorporate it with Marquardt inversion. The third approach corresponds to Marquardt inversion introducing lateral constraints. Finally, we include the lateral constraints in Differential Evolution Adaptive Metropolis optimization, and the results are sequentially utilized by Marquardt inversion. Occam–Marquardt may provide accurate estimation of the subsurface features, but it is dependent on the appropriate conversion of different multi‐layered Occam model to an acceptable starting model for Marquardt inversion, which is not straightforward. Employing parameter spaces, the Differential Evolution Adaptive Metropolis approach can be pertinent to determine Marquardt a priori information; nevertheless, the uncertainties in Differential Evolution Adaptive Metropolis optimization will introduce some inaccuracies in Marquardt inversion results. Laterally constrained Marquardt may be promising to resolve sub‐seafloor features, but it is not stable if there are significant lateral changes of the sub‐seafloor structure due to the dependence of the method to the starting model. Including the lateral constraints in Differential Evolution Adaptive Metropolis approach allows for faster convergence of the routine with consistent results, furnishing more accurate estimation of a priori models for the subsequent Marquardt inversion.     

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.     

Controlled‐source electromagnetic monitoring of reservoir oil saturation using a novel borehole‐to‐surface configuration

To advance and optimize secondary and tertiary oil recovery techniques, it is essential to know the areal propagation and distribution of the injected fluids in the subsurface. We investigate the applicability of controlled‐source electromagnetic methods to monitor fluid movements in a German oilfield (Bockstedt, onshore Northwest Germany) as injected brines (highly saline formation water) have much lower electrical resistivity than the oil within the reservoir. The main focus of this study is on controlled‐source electromagnetic simulations to test the sensitivity of various source–receiver configurations. The background model for the simulations is based on two‐dimensional inversion of magnetotelluric data gathered across the oil field and calibrated with resistivity logs. Three‐dimensional modelling results suggest that controlled‐source electromagnetic methods are sensitive to resistivity changes at reservoir depths, but the effect is difficult to resolve with surface measurements only. Resolution increases significantly if sensors or transmitters can be placed in observation wells closer to the reservoir. In particular, observation of the vertical electric field component in shallow boreholes and/or use of source configurations consisting of combinations of vertical and horizontal dipoles are promising. Preliminary results from a borehole‐to‐surface controlled‐source electromagnetic field survey carried out in spring 2014 are in good agreement with the modelling studies.     

Three‐dimensional inversion of towed streamer electromagnetic data

A towed streamer electromagnetic system capable of simultaneous seismic and electromagnetic data acquisition has recently been developed and tested in the North Sea. We introduce a 3D inversion methodology for towed streamer electromagnetic data that includes a moving sensitivity domain. Our implementation is based on the 3D integral equation method for computing responses and Fréchet derivatives and uses the re‐weighted regularized conjugate gradient method for minimizing the objective functional with focusing regularization. We present two model studies relevant to hydrocarbon exploration in the North Sea. First, we demonstrate the ability of a towed electromagnetic system to detect and characterize the Harding field, a medium‐sized North Sea hydrocarbon target. We compare our 3D inversion of towed streamer electromagnetic data with 3D inversion of conventional marine controlled‐source electromagnetic data and observe few differences between the recovered models. Second, we demonstrate the ability of a towed streamer electromagnetic system to detect and characterize the Peon discovery, which is representative of an infrastructure‐led shallow gas play in the North Sea. We also present an actual case study for the 3D inversion of towed streamer electromagnetic data from the Troll field in the North Sea and demonstrate our ability to image all the Troll West Oil and Gas Provinces and the Troll East Gas Province. We conclude that 3D inversion of data from the current generation of towed streamer electromagnetic systems can adequately recover hydrocarbon‐bearing formations to depths of approximately 2 km. We note that by obviating the need for ocean‐bottom receivers, the towed streamer electromagnetic system enables electromagnetic data to be acquired over very large areas in frontier and mature basins for higher acquisition rates and relatively lower cost than conventional marine controlled‐source electromagnetic methods.     

海洋电磁低时钟漂移及自动增益采集技术研究

作为地震探测的有效补充,电磁法在海洋油气资源勘探中发挥着越来越重要的作用.本文针对海洋电磁法中极低时钟误差、大动态范围采集等要求,提出使用GPS和高精度原子钟解决低时钟漂移问题;采集电路中设计增益可调的放大电路,MCU通过对采集到的信号进行实时分析,实现增益动态调整,解决海洋电磁信号大动态范围采集问题;研制的海洋电磁信号数据记录单元每天时间误差小于0.3 ms、电场通道动态范围可以达到160 dB、磁场通道动态范围可以达到134 dB.通过在室内指标测试、室外试验,表明设计的数据记录单元能够稳定可靠工作.为证明数据纪录单元的有效性,将数据记录单元和国外商用MT仪器进行野外一致性对比实验,数据记录单元与商用仪器结果基本一致;在水深100 m海域进行了海底MT信号采集,得到了有效的结果.该研究为海洋电磁信号低时钟漂移、大动态范围采集提供了一种有效解决方案.     

海洋拖曳式水平电偶源数值模拟与电场接收机研制

海洋天然气水合物资源调查中,电磁法已成为地震勘探最有效的补充手段之一,而我国海洋水合物电磁探测中的拖曳式电场接收机此前是空白.本文从方法原理、仪器设备和海洋试验等方面,对海洋拖曳式水平电偶极-偶极方法进行研讨.建立三维电阻率模型,进行数值模拟,讨论收发距、发射频率、拖曳装置离海底高度等参数与电场异常响应的关系,针对特定模型给出了最佳观测参数以指导仪器硬件设计与海洋测试.在此基础上,采用高级嵌入式计算机的智能化控制、前端低噪声斩波放大、高精度同步采集等技术,研制了由拖体、承压舱、采集电路、低噪声电场传感器等组成的海洋拖曳式水平轴向电场接收机,其系统噪声达到2 nV/√HZ @(1~100Hz).将自主研制的电场接收机用于国内首次拖曳式电偶极-偶极方法的海洋试验与评估,并对海试资料进行了处理与分析.研究结果表明,所研发的仪器达到了设计指标,能够用于海洋天然气水合物调查.     

Multi‐dimensional analyses of the SEAM controlled source electromagnetic data—the story of a blind test of interpretation workflows

Using a subset of the SEG Advanced Modeling Program Phase I controlled‐source electromagnetic data, we apply our standard controlled‐source electromagnetic interpretation workflows to delineate a simulated hydrocarbon reservoir. Experience learned from characterizing such a complicated model offers us an opportunity to refine our workflows to achieve better interpretation quality. The exercise proceeded in a blind test style, where the interpreting geophysicists did not know the true resistivity model until the end of the project. Rather, the interpreters were provided a traditional controlled‐source electromagnetic data package, including electric field measurements, interpreted seismic horizons, and well log data. Based on petrophysical analysis, a background resistivity model was established first. Then, the interpreters started with feasibility studies to establish the recoverability of the prospect and carefully stepped through 1D, 2.5D, and 3D inversions with seismic and well log data integrated at each stage. A high‐resistivity zone is identified with 1D analysis and further characterized with 2.5D inversions. Its lateral distribution is confirmed with a 3D anisotropic inversion. The importance of integrating all available geophysical and petrophysical data to derive more accurate interpretation is demonstrated.     

The marine controlled source electromagnetic response of a steel borehole casing: applications for the NEPTUNE Canada gas hydrate observatory

Gas hydrates are a potential energy resource, a possible factor in climate change and an exploration geohazard. The University of Toronto has deployed a permanent seafloor time‐domain controlled source electromagnetic (CSEM) system offshore Vancouver Island, within the framework of the NEPTUNE Canada underwater cabled observatory. Hydrates are known to be present in the area and due to their electrically resistive nature can be monitored by 5 permanent electric field receivers. However, two cased boreholes may be drilled near the CSEM site in the near future. To understand any potential distortions of the electric fields due to the metal, we model the marine electromagnetic response of a conductive steel borehole casing. First, we consider the commonly used canonical model consisting of a 100 Ωm, 100 m thick resistive hydrocarbon layer embedded at a depth of 1000 m in a 1 Ωm conductive host medium, with the addition of a typical steel production casing extending from the seafloor to the resistive zone. Results show that in both the frequency and time domains the distortion produced by the casing occurs at smaller transmitter‐receiver offsets than the offsets required to detect the resistive layer. Second, we consider the experimentally determined model of the offshore Vancouver Island hydrate zone, consisting of a 5.5 Ωm, 36 m thick hydrate layer overlying a 0.7 Ωm sedimentary half‐space, with the addition of two borehole casings extending 300 m into the seafloor. In this case, results show that the distortion produced by casings located within a 100 m safety zone of the CSEM system will be measured at 4 of the 5 receivers. We conclude that the boreholes must be positioned at least 200 m away from the CSEM array so as to minimize the effects of the casings.     

Two‐dimensional controlled‐source electromagnetic interferometry by multidimensional deconvolution: spatial sampling aspects

We use numerically modelled data sets to investigate the sensitivity of electromagnetic interferometry by multidimensional deconvolution to spatial receiver sampling. Interferometry by multidimensional deconvolution retrieves the reflection response below the receivers after decomposition of the fields into upward and downward decaying fields and deconvolving the upward decaying field by the downward decaying field. Thereby the medium above the receiver level is replaced with a homogeneous half‐space, the sources are redatumed to the receiver level and the direct field is removed. Consequently, in a marine setting the retrieved reflection response is independent of any effect of the water layer and the air above. A drawback of interferometry by multidimensional deconvolution is a possibly unstable matrix inversion, which is necessary to retrieve the reflection response. Additionally, in order to correctly separate the upward and the downward decaying fields, the electromagnetic fields need to be sampled properly. We show that the largest possible receiver spacing depends on two parameters: the vertical distance between the source and the receivers and the length of the source. The receiver spacing should not exceed the larger of these two parameters. Besides these two parameters, the presence of inhomogeneities close to the receivers may also require a dense receiver sampling. We show that by using the synthetic aperture concept, an elongated source can be created from conventionally acquired data in order to overcome these strict sampling criteria. Finally, we show that interferometry may work under real‐world conditions with random noise and receiver orientation and positioning errors.     

Insight into the marine controlled‐source electromagnetic signal propagation

In marine controlled‐source electromagnetic (CSEM) surveys the subsurface is explored by emitting low‐frequency signals from an electric dipole source close to the sea‐bed. The main goal is often to detect and describe possible thin resistive layers beneath the sea‐bed. To gain insight into how CSEM signals propagate, it is informative to study a stratified model. The electromagnetic field is then given in terms of integrals over TE‐ and TM‐polarized plane‐wave constituents. An asymptotic evaluation of the field integrals for large propagation distances results in explicit spatial expressions for the field components and the derived expressions can be used to analyse how the CSEM signals propagate. There are two major signal pathways in a standard CSEM model. One of these pathways is via the thin resistive layer and the resulting response is accounted for by a pole in the reflection response for the TM mode. The signal is propagating nearly vertically down to the resistor from the source, then guided while attenuated along the reservoir, before propagating nearly vertically up to the receiver. The response is slightly altered by the sea‐bed interface and further modified in shallow water due to multiple reflections between the sea‐surface and sea‐bed at both the source and receiver sides. The other major signal pathway is via the resistive air half‐space, the so‐called airwave. The airwave is generated by the TE mode and interacts with the subsurface via vertically propagating signals reflected between the sea‐surface and subsurface at both the source and receiver sides.     

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.     

An improved interpolation scheme at receiver positions for 2.5D frequency-domain marine controlled-source electromagnetic forward modelling

An improved interpolation scheme is presented for 2.5-dimensional marine controlled-source electromagnetic forward modelling. For the marine controlled-source electromagnetic method, due to the resistivity contrast between the seawater and seafloor sedimentary layers, it is difficult to compute the electromagnetic fields accurately at receivers, which are usually located at the seafloor. In this study, the 2.5-dimensional controlled-source electromagnetic responses are simulated by the staggered finite-difference method. The secondary-field approach is used to avoid the source singularities, and the one-dimensional layered background model is used for calculating the primary fields excited by the source quasi-analytically. The interpolation of electromagnetic fields at the cell nodes for the whole computational domain to the receiver locations is discussed in detail. Numerical tests indicate that the improved interpolation developed is more accurate for simulating the electromagnetic responses at receivers located at the seafloor, compared with the linear or rigorous interpolation.     

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.