Elongated horizontal and vertical receivers in three-dimensional electromagnetic modelling and inversion

Electromagnetic geophysical methods often rely on measurements of naturally occurring or artificially impressed electric fields. It is technically impossible, however, to measure the electric field directly. Instead, the electric field is approximated by recording the voltage difference between two electrodes and dividing the obtained voltage by the distance between the electrodes. Typically, modelling and inversion algorithms assume that the electric fields are obtained over infinitely short point-dipoles and thus measured fields are assigned to a single point between the electrodes. Such procedures imply several assumptions: (1) The electric field between the two electrodes is regarded as constant or being a potential field and (2) the receiver dimensions are negligible compared to the dimensions of the underlying modelling grid. While these conditions are often fulfilled for horizontal electric fields, borehole sensors for recordings of the vertical electric field have dimensions in the order of ≈100 m and span several modelling grid cells. Observations from such elongated borehole sensors can therefore only be interpreted properly if true receiver dimensions and variations of electrical conductivity along the receiver are considered. Here, we introduce a numerical solution to include the true receiver geometry into electromagnetic modelling schemes, which does not rely on such simplifying assumptions. The algorithm is flexible, independent of the chosen numerical method to solve Maxwell's equations and can easily be implemented in other electromagnetic modelling and inversion codes. We present conceptual modelling results for land-based controlled source electromagnetic scenarios and discuss consideration of true receiver geometries for a series of examples of horizontal and vertical electric field measurements. Comparison with Ez data measured in an observation borehole in a producing oil field shows the importance of both considering the true length of the receiver and also its orientation. We show that misalignment from the vertical axis as small as 0.1° may seriously distort the measured signal, as horizontal electric field components are mapped into the desired vertical component. Adequate inclusion of elongated receivers in modelling and inversion can also help reducing effects of static shift when interpreting (natural source) magnetotelluric data.     

Application of Hilbert‐like transforms for enhanced processing of full tensor magnetic gradient data

Commonly, geomagnetic prospection is performed via scalar magnetometers that measure values of the total magnetic intensity. Recent developments of superconducting quantum interference devices have led to their integration in full tensor magnetic gradiometry systems consisting of planar‐type first‐order gradiometers and magnetometers fabricated in thin‐film technology. With these systems measuring directly the magnetic gradient tensor and field vector, a significantly higher magnetic and spatial resolution of the magnetic maps is yield than those produced via conventional magnetometers. In order to preserve the high data quality in this work, we develop a workflow containing all the necessary steps for generating the gradient tensor and field vector quantities from the raw measurement data up to their integration into high­resolution, low­noise, and artefactless two‐dimensional maps of the magnetic field vector. The gradient tensor components are processed by superposition of the balanced gradiometer signals and rotation into an Earth‐centred Earth‐fixed coordinate frame. As the magnetometers have sensitivity lower than that of gradiometers and the total magnetic intensity is not directly recorded, we employ Hilbert‐like transforms, e.g., integration of the gradient tensor components or the conversion of the total magnetic intensity derived by calibrated magnetometer readings to obtain these values. This can lead to a better interpretation of the measured magnetic anomalies of the Earth's magnetic field that is possible from scalar total magnetic intensity measurements. Our conclusions are drawn from the application of these algorithms on a survey acquired in South Africa containing full tensor magnetic gradiometry data.     

用模式匹配算法研究层状各向异性倾斜地层中多分量感应测井响应

利用数值模式匹配技术,研究并建立层状各向异性倾斜地层中多分量感应测井响应的快速数值模拟算法.首先将位于井轴上的三个相互正交磁偶极子转化成关于极角θ的三个谐变分量的叠加,以便将三个正交磁偶极子电磁场的正演问题完全简化成三个谐变分量电磁场的轴对称定解问题,并给出电磁场各个谐变分量在井轴上满足的边界条件,保证电磁场在井轴附近仍然可解.然后利用数值模式匹配技术建立电磁场各个谐变分量的正演过程,得到层状各向异性地层中磁流源并矢Green函数的半解析表达式,给出计算层状各向异性倾斜地层中多分量感应测井响应的具体方法,最后通过数值计算结果证明该算法的有效性并考察几种不同情况下多分量感应测井响应特征.     

Magnetic data analysis at low latitudes using magnitude transforms

Magnitude transforms include magnitude magnetic anomalies (MMA), their gradients and Laplacians. They can be calculated from the total magnetic anomalies or other component anomalies. Magnitude magnetic anomalies have a space distribution different from that of the component anomalies. Their values are non‐negative and their respective patterns are similar to the positive gravity anomalies. Magnitude transforms are an effective tool for magnetic data analysis due to their simplified pattern and direct correlation with the space location of the source. They have advantages over the traditionally used reduction‐to‐the‐pole (RTP) transform, especially at low magnetic latitudes. The calculation of magnitudes of the anomalous field requires the total field data to be transformed into the component anomalies, while the reduction‐to‐the‐pole transform also includes a rotation of the magnetization vector, the orientation of which is usually assumed. For equal latitudes, the transfer functions of component‐component transforms in the frequency domain show better stability than the component‐component‐rotation transfer function. This is illustrated by a comparison of analytical expressions, and synthetic models of magnetic fields. The Dixon seamount case shows the possibilities for an improved data analysis and more confident source recognition at low latitudes using magnetic transforms.     

Two-Dimensional Magnetotelluric Modelling of Ore Deposits: Improvements in Model Constraints by Inclusion of Borehole Measurements

A combination of magnetotelluric (MT) measurements on the surface and in boreholes (without metal casing) can be expected to enhance resolution and reduce the ambiguity in models of electrical resistivity derived from MT surface measurements alone. In order to quantify potential improvement in inversion models and to aid design of electromagnetic (EM) borehole sensors, we considered two synthetic 2D models containing ore bodies down to 3000 m depth (the first with two dipping conductors in resistive crystalline host rock and the second with three mineralisation zones in a sedimentary succession exhibiting only moderate resistivity contrasts). We computed 2D inversion models from the forward responses based on combinations of surface impedance measurements and borehole measurements such as (1) skin-effect transfer functions relating horizontal magnetic fields at depth to those on the surface, (2) vertical magnetic transfer functions relating vertical magnetic fields at depth to horizontal magnetic fields on the surface and (3) vertical electric transfer functions relating vertical electric fields at depth to horizontal magnetic fields on the surface. Whereas skin-effect transfer functions are sensitive to the resistivity of the background medium and 2D anomalies, the vertical magnetic and electric field transfer functions have the disadvantage that they are comparatively insensitive to the resistivity of the layered background medium. This insensitivity introduces convergence problems in the inversion of data from structures with strong 2D resistivity contrasts. Hence, we adjusted the inversion approach to a three-step procedure, where (1) an initial inversion model is computed from surface impedance measurements, (2) this inversion model from surface impedances is used as the initial model for a joint inversion of surface impedances and skin-effect transfer functions and (3) the joint inversion model derived from the surface impedances and skin-effect transfer functions is used as the initial model for the inversion of the surface impedances, skin-effect transfer functions and vertical magnetic and electric transfer functions. For both synthetic examples, the inversion models resulting from surface and borehole measurements have higher similarity to the true models than models computed exclusively from surface measurements. However, the most prominent improvements were obtained for the first example, in which a deep small-sized ore body is more easily distinguished from a shallow main ore body penetrated by a borehole and the extent of the shadow zone (a conductive artefact) underneath the main conductor is strongly reduced. Formal model error and resolution analysis demonstrated that predominantly the skin-effect transfer functions improve model resolution at depth below the sensors and at distance of \(\sim \) 300–1000 m laterally off a borehole, whereas the vertical electric and magnetic transfer functions improve resolution along the borehole and in its immediate vicinity. Furthermore, we studied the signal levels at depth and provided specifications of borehole magnetic and electric field sensors to be developed in a future project. Our results suggest that three-component SQUID and fluxgate magnetometers should be developed to facilitate borehole MT measurements at signal frequencies above and below 1 Hz, respectively.     

Analysis of a guided wave along a conducting structure in a borehole

A guided wave along a borehole is often observed in borehole radar measurements. These guided waves deform the antenna pattern and can cause artefacts in radar measurements. A water-filled borehole or a conducting logging cable can function as a waveguide for electromagnetic waves under some conditions. We describe the theoretical characteristics of such a guided wave in a borehole and compare them with our experiments. The measured signal discussed was obtained with a directional borehole radar. This radar uses a cylindrical conformal array antenna as receiver and is a model of a conducting structure in a borehole. The induced field around the borehole was compared with the theory. The most fundamental symmetrical and asymmetrical modes were TM₀₁ and HE₁₁, and they were identified in the measured signals using time–frequency distribution analysis and by observation of the azimuthal field distribution of the magnetic field.     

An alternative approach in establishing relation between vertical and horizontal gradients of 2D potential field

The relation in which the vertical and horizontal gradients of potential field data measured along a profile across a two‐dimensional source are a Hilbert transform pair is re‐established using complex domain mathematics. In addition, a relation between the measured field and its vertical gradient in terms of a closed‐form formula is also established. The formula is based on hypersingular or Hadamard's finite‐part integral. To estimate the vertical gradient directly from the field data, Linz's algorithm of computing Hadamard's finite‐part integral is implemented. Numerical experiments are conducted on synthetically generated total magnetic intensity data with a mild level of noise contamination. A model of a magnetically polarised vertical thin sheet buried at a finite depth within a non‐magnetic half‐space was considered in generating the synthetic response. The results from numerical experiments on the mildly noise‐contaminated synthetic response are compared with those from using classical Fourier and robust regularised Hilbert transform‐based techniques.     

Magnetic mapping for the detection and characterization of UXO: Use of multi-sensor fluxgate 3-axis magnetometers and methods of interpretation

Fluxgate vector magnetometers can be powerful instruments in magnetic mapping and in the detection of unexploded ordnance (UXO). Being lightweight and having low power requirements, they allow the development of multi-sensor use. The main problems in using fluxgate magnetometers arise from calibration errors and drift but these can be overcome using a quick and simple method of calibration in the field. This method also has the advantage of compensating permanent and induced magnetic fields generated by magnetized objects carried with the sensors. Measurement accuracy is similar to that obtained with scalar magnetometers. Multi-magnetometer profiling allows direct inversion of raw magnetic data along profiles in order to locate and characterize dipoles typically generated by UXO. The method used is a non-linear inverse procedure that estimates the three coordinates and magnetization vector of the dipole. Application of the technique to both a synthetic case and a field example illustrates the advantages of the method compared with the use of the analytic signal.     

Electric field data in inductive source electromagnetic surveys

Measurement of the electric field data due to an inductive loop source in a controlled source electromagnetic survey is not common, because electric field data, usually involving grounded electrodes, are expensive to acquire and difficult to interpret. With the recently developed capability of versatile three‐dimensional inversion, we revisit the idea of measuring electric field in a large ground loop survey for mineral exploration. The three‐dimensional modelling and inversion approach helps us quantitatively understand the detectability and recoverability of the proposed survey configuration. Our detectability study using forward modelling shows that the relative anomaly (percentage difference) in electric field does not decay with a lower induction number, but the conventional magnetic field data (dB/dt) does. Our recoverability study examines how much and what kind of information can be extracted from electric field data for the reconstruction of a three‐dimensional model. Synthetic inversions show the following observations. (i) Electric field data are good at locating lateral discontinuity, whereas dB/dt has better depth resolution. (ii) Electric field is less sensitive to the background conductivity and, thus, is prone to misinterpretation because of a bad initial model in inversion. We recommend warm‐starting the electric field inversion with an initial model from a separate dB/dt inversion. (iii) Electric field data may be severely contaminated by near‐surface heterogeneity, but an inversion can recover the deep target concealed by the geologic noise. (iv) Even one line of single‐component electric field data can greatly improve the horizontal resolution in a dB/dt inversion. Finally, we investigate a field dataset of both electric field and dB/dt measurements at a uranium deposit. The field example confirms that the electric field and magnetic field data contain independent information that is crucial in the accurate recovery of subsurface conductivity. Our synthetic and field examples demonstrate the benefit of acquiring electric field data along with magnetic field data in an inductive source survey.     

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.     

Field instrumentation for high‐resolution parallel monitoring of bedrock erosion and bedload transport

River incision is fundamental in shaping the Earth's surface. In mountainous regions with steep river beds, fluvial bedrock erosion by bedload transport is an important mechanism forming channels. However, there are only a few complete field datasets that can be used to improve process understanding and evaluate erosion models, especially at the process scale. To provide a simultaneous dataset of hydraulics, bedload transport and bedrock erosion at high temporal and spatial resolution, a new measuring device has been installed in the Erlenbach, a gauged stream in the Swiss Pre‐Alps. In this stream, bedload transport rates can be calculated from surveying deposits and from geophone plate sensors and bedload transport samples can be taken directly by an automated moving basket system. To measure bedrock erosion rates simultaneously, two natural stone slabs were mounted flush with the channel bed in a steel frame hosting various measurement devices. Force sensors below the slabs record normal stress and shear stress. At‐a‐point erosion rates on the slab surfaces are continuously measured at sub‐millimetre precision at three locations on each slab. In addition, the slab topography is monitored following erosive flood events. In this article (i) the ‘erosion scale’ device is described, (ii) data resolution and data quality is assessed by means of tests and event data, and (iii) the first transport event is discussed. The erosion scales are confirmed to provide data at high spatio‐temporal resolution for process analysis. The preliminary data show evidence for the tools effect in bedrock erosion. The bedrock slabs can be exchanged to obtain measurements for catchments with different lithologies for comparison. Copyright © 2014 John Wiley & Sons, Ltd.     

强剩磁强退磁条件下的二维井中磁测反演

强剩磁、强退磁改变了总磁化强度的大小和方向,给磁测资料解释带来困难.为此,本文利用二维井中磁测数据反演磁化强度矢量的二维分布.首先利用井中磁测的磁异常模量反演磁化强度大小的分布.然后,在已知磁化强度大小分布的前提下,拟合磁场分量,反演磁化强度方向的分布.其中,磁化强度大小和方向均用共轭梯度法求解,并通过预优矩阵改善磁化强度大小的反演效果.理论模拟说明,该方法能准确获得磁化强度矢量分布.磁化强度矢量反演结果包括感磁、剩磁及退磁的影响,这为研究强剩磁、高磁化率矿床提供了一种有效方法.     

An effective data processing flow for the acoustic reflection image logging

Recent studies have revealed the great potential of acoustic reflection logging in detecting near borehole fractures and vugs. The new design of acoustic reflection imaging tool with a closest spacing of 10.6m and a certain degree of phase steering makes it easier to extract the reflection signals from the borehole mode waves. For field applications of the tool, we had developed the corresponding processing software: Acoustic Reflection Imaging. In this paper, we have further developed an effective data processing flow by employing multi‐scale slowness‐time‐coherence for reflection wave extraction and incorporating reverse time migration for imaging complicated subtle structures with the strong effects of borehole environment. Applications of the processing flow to synthetic data of acoustic reflection logging in a fractured formation model and interface model with fluid filled borehole generated by 2D finite difference method, and to the physical modelling data from a laboratory water tank, as well as to the field data from two wells in a western Chinese oil field, demonstrate the validity and capability of our multi‐scale slowness‐time‐coherence and reverse time migration algorithms.     

ASPI experiment: measurements of fields and waves on board the INTERBALL-1 spacecraft

The plasma-wave experiment ASPI (analysis of spectra of plasma waves and instabilities) on board the INTERBALL spacecraft is a combined wave diagnostics experiment. It performs measurements of the DC and AC magnetic field vector by flux-gate and searchcoil sensors, the DC and AC electric field vector by Langmuir double probes and the plasma current by Langmuir split probe. Preliminary data analysis shows the low noise levels of the sensors and the compatibility of new data with the results of previous missions. During several months of in-orbit operation a rich collection of data was acquired, examples of which at the magnetopause and plasma sheet are presented in second part of the paper.     

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.     

根据钻孔崩落资料确定剑川地区应力场方向

本文论述了钻孔崩落的力学机制、产生环境及其与现今构造应力场的关系,介绍了超声波井下电视的测井原理、操作技术和钻孔崩落测井资料分析方法。根据一口800米深钻孔的超声波井下电视测井获得的崩落资料,确定剑川地震活动区最大水平主应力方向为N20°E。这与近两年来在该区采用水压致裂法测定的最大水平主压应力方向N15°E是十分一致的     

利用2.5维数值模式匹配算法高效高精度建立垂直井眼中多分量阵列感应井眼校正库

本文将应用2.5维数值模式匹配算法研究建立柱状横向同性地层中偏心条件下磁流源并矢Green函数的高效算法,并利用该算法高效计算多分量阵列感应井眼校正库.首先,通过Fourier展开技术将偏心条件下磁流源并矢Green函数的数值模拟转化为一系列的轴对称问题;然后利用模式匹配算法求解轴对称问题,得到柱状介质中偏心条件下磁流源并矢Green函数的半解析解;最后,结合中国石油集团测井有限公司开发的三维感应测井仪器(TDIT)的结构参数,推导出复合线圈系视电导率张量的计算公式.在此基础上,根据井场实际地层资料归纳出理论地层模型相应参数的变化范围,高效高精度建立一套直接应用于工业生产的多分量阵列感应井眼校正库,并利用井眼校正库详细地分析考察各模型参数对TDIT响应的影响.     

Interaction between plasma structures in an unstable ionospheric plasma

A vortex structure renders additional stability to plasma irregularities stretched along magnetic field lines. Plasma irregularities extended over several tens of kilometers are registered with rocket and satellite equipment in the topside ionosphere. The registered scale of irregularities depends on the spatial and time resolution of the equipment used. Irregular structures were registered in the ionosphere during experiments with barium clouds and jets, when a plasma irregularity separated into strata extended over several meters and several kilometers across the geomagnetic field. It has been indicated that plasma vortices can be generated in an unstable plasma in a situation when its quasi-neutrality is disturbed. Local geomagnetic field disturbances will be caused by the appearance of a proper vortex magnetic field. Plasma vortices can interact in an inhomogeneous plasma with an unstable electron component. Such interactions are related to the transformation of the phase volume of free electrostatic oscillations in the frequency-wave vector space.     

Optimization approach to automatic first arrival picking for three‐component three‐dimensional vertical seismic profiling data

A new, adaptive multi‐criteria method for accurate estimation of three‐component three‐dimensional vertical seismic profiling of first breaks is proposed. Initially, we manually pick first breaks for the first gather of the three‐dimensional borehole set and adjust several coefficients to approximate the first breaks wave‐shape parameters. We then predict the first breaks for the next source point using the previous one, assuming the same average velocity. We follow this by calculating an objective function for a moving trace window to minimize it with respect to time shift and slope. This function combines four main properties that characterize first breaks on three‐component borehole data: linear polarization, signal/noise ratio, similarity in wave shapes for close shots and their stability in the time interval after the first break. We then adjust the coefficients by combining current and previous values. This approach uses adaptive parameters to follow smooth wave‐shape changes. Finally, we average the first breaks after they are determined in the overlapping windows. The method utilizes three components to calculate the objective function for the direct compressional wave projection. An adaptive multi‐criteria optimization approach with multi three‐component traces makes this method very robust, even for data contaminated with high noise. An example using actual data demonstrates the stability of this method.     

Estimating magnetic dike parameters using a non‐linear constrained inversion technique: an example from the Särna area,west central Sweden

In this paper, we describe a non‐linear constrained inversion technique for 2D interpretation of high resolution magnetic field data along flight lines using a simple dike model. We first estimate the strike direction of a quasi 2D structure based on the eigenvector corresponding to the minimum eigenvalue of the pseudogravity gradient tensor derived from gridded, low‐pass filtered magnetic field anomalies, assuming that the magnetization direction is known. Then the measured magnetic field can be transformed into the strike coordinate system and all magnetic dike parameters – horizontal position, depth to the top, dip angle, width and susceptibility contrast – can be estimated by non‐linear least squares inversion of the high resolution magnetic field data along the flight lines. We use the Levenberg‐Marquardt algorithm together with the trust‐region‐reflective method enabling users to define inequality constraints on model parameters such that the estimated parameters are always in a trust region. Assuming that the maximum of the calculated g_zz (vertical gradient of the pseudogravity field) is approximately located above the causative body, data points enclosed by a window, along the profile, centred at the maximum of g_zz are used in the inversion scheme for estimating the dike parameters. The size of the window is increased until it exceeds a predefined limit. Then the solution corresponding to the minimum data fit error is chosen as the most reliable one. Using synthetic data we study the effect of random noise and interfering sources on the estimated models and we apply our method to a new aeromagnetic data set from the Särna area, west central Sweden including constraints from laboratory measurements on rock samples from the area.