On the exploration of a marine aquifer offshore Israel by long-offset transient electromagnetics

The existence of aquifers extending from land beneath the sea floor up to a distance of several kilometres has been observed and examined all over the world. The coastal aquifer of Israel is a heavily used groundwater reservoir which has to be constantly monitored to ensure the drinking water supply. Former land-based electromagnetic measurements show that it is, in several places, blocked to seawater intrusion and is consequently a candidate for submarine extension. Multicomponent long-offset transient electromagnetic measurements were carried out offshore on the coast of Israel. We deployed a 400-m-long grounded dipole as transmitter and several electric and magnetic receivers on the sea floor up to a distance of 4.8 km from the coast. Altogether, we deployed 8 transmitter positions and received data sets at 14 receiver stations onshore and offshore, with offsets of mostly 400 and 800 m. In this paper, we present the survey and 1D Occam and Marquardt inversions of the offshore horizontal electric components in the broadside and inline configuration. In addition, the vertical magnetic component in the broadside position is also considered. Only single inversions, both single offset and single component, were used to detect the aquifer under sea sediments. We prove the submarine existence of the Israeli coastal aquifer up to a distance to the coast of approximately 3.2 to 3.6 km using all measured long-offset transient electromagnetic components. In addition, we present modelling studies with synthetic data derived from a subsurface model adjusted to our survey area with very shallow water from 10 to 50 m. As part of the planning before the survey, a parameter study of the expected subsurface, the examination of the airwave phenomenon and the justification for our 1D inversion strategy are shown. More detailed eigen parameter analyses are added to explain the measured data.     

Investigating the effect of data quality on time domain electromagnetic discrimination

Using field data and numerical simulations we investigate the effect of data quality on time domain electromagnetic discrimination. Data quality decreases when measurements contain responses not accounted for by our mathematical modelling. This can include instrument noise, inaccurately reported position and orientation information, geologic contributions to the signal, and loss of validity of the forward modelling. Survey design is critical to data quality in order to have sufficient sampling of data anomalies, and also to ensure that each target is illuminated such that both the axial and transverse components of the polarization can be excited and measured. For dipole model based discrimination algorithms, success is contingent upon the accuracy with which the components of the polarization tensor can be estimated. Field data from different survey modes are analysed to identify noise sources and provide quantitative estimates of the noise in each survey. Inversion results show that increased noise levels lead to greater spread in recovered parameters. Monte Carlo simulations are performed in order to investigate the importance of other data quality factors. Analysis of inversion results from the simulations show that anomaly size, signal to noise ratio, positioning error, line spacing and station spacing all play a role in the spread of recovered parameters. Through the analysis of our simulation results we propose a figure of merit as a means of quantifying different data quality factors with a single number and relate this number to the accuracy with which parameters can be estimated.     

Application of a library based method to time domain electromagnetic data for the identification of unexploded ordnance

A simple library based algorithm for the identification of unexploded ordnance is tested on time domain electromagnetic data. A library of polarization tensors is generated from data acquired on a test stand over a collection of different UXO. The objective of the algorithm is to determine which target of the library is most likely to have produced an observed data anomaly. For each target in our library a non-linear inverse problem is solved for the position and orientation that minimizes the least-squares difference between the observed data anomaly and the data predicted from each target. This technique avoids direct inversion for polarization tensor, making it feasible for cases where sensor data quality may not be sufficient to support confident estimation of model parameters. For cases where the background soil response is significant, we also invert for the t⁻¹ decay characteristic of viscous remnant magnetic soil. We present results of preliminary tests of the library technique to Geonics EM63 time domain electromagnetic data collected on a test plot seeded with UXO. These tests demonstrate an excellent ability to accurately identify isolated targets. Misidentification of single targets occur when data anomalies have low signal to noise ratios or when targets within the library have similar polarization tensors.     

海底油气藏及天然气水合物的时频电磁辨识

提出了伪随机多频海洋电磁法观测方案.采用伪随机多频信号作为激励场源,多偏移距同线偶极-偶极同时观测,相关辨识海底地电系统的频率特性和冲激响应,可以在时间域和频率域同时辨识海底高阻薄层.在时间域,瞬变冲激时刻可以直接指示海底地层电导率的变化;在频率域,利用多个频率的电场响应计算的频散率及其道闻变化量,相对相位道间变化量对高阻薄层有很好的反映.从而实现对海底油气及天然气水合物的多参数辨识.     

时间域航空电磁2.5维非线性共轭梯度反演

对于时间域航空电磁法二维和三维反演来说,最大的困难在于有效的算法和大的计算量需求.本文利用非线性共轭梯度法实现了时间域航空电磁法2.5维反演方法,着重解决了迭代反演过程中灵敏度矩阵计算、最佳迭代步长计算、初始模型选取等问题.在正演计算中,我们采用有限元法求解拉式傅氏域中的电磁场偏微分方程,再通过逆拉氏和逆傅氏变换高精度数值算法得到时间域电磁响应.在灵敏度矩阵计算中,采用了基于拉式傅氏双变换的伴随方程法,时间消耗只需计算两次正演,从而节约了大量计算时间.对于最佳步长计算,二次插值向后追踪法能够保证反演迭代的稳定性.设计两个理论模型,检验反演算法的有效性,并讨论了选择不同初始模型对反演结果的影响.模型算例表明：非线性共轭梯度方法应用于时间域航空电磁2.5维反演中稳定可靠,反演结果能够有效地反映地下真实电性结构.当选择的初始模型电阻率值与真实背景电阻率值接近时,能得到较好的反演结果,当初始模型电阻率远大于或远小于真实背景电阻率值时反演效果就会变差.     

Petrophysics and mineral exploration: a workflow for data analysis and a new interpretation framework

As mineral exploration seeks deeper targets, there will be a greater reliance on geophysical data and a better understanding of the geological meaning of the responses will be required, and this must be achieved with less geological control from drilling. Also, exploring based on the mineral system concept requires particular understanding of geophysical responses associated with altered rocks. Where petrophysical datasets of adequate sample size and measurement quality are available, physical properties show complex variations, reflecting the combined effects of various geological processes. Large datasets, analysed as populations, are required to understand the variations. We recommend the display of petrophysical data as frequency histograms because the nature of the data distribution is easily seen with this form of display. A petrophysical dataset commonly contains a combination of overlapping sub-populations, influenced by different geological factors. To understand the geological controls on physical properties in hard rock environments, it is necessary to analyse the petrophysical data not only in terms of the properties of different rock types. It is also necessary to consider the effects of processes such as alteration, weathering, metamorphism and strain, and variables such as porosity and stratigraphy. To address this complexity requires that much more supporting geological information be acquired than in current practice. The widespread availability of field portable instruments means quantitative geochemical and mineralogical data can now be readily acquired, making it unnecessary to rely primarily on categorical rock classification schemes. The petrophysical data can be combined with geochemical, petrological and mineralogical data to derive explanations for observed physical property variations based not only on rigorous rock classification methods, but also in combination with quantitative estimates of alteration and weathering. To understand how geological processes will affect different physical properties, it is useful to define three end-member forms of behaviour. Bulk behaviour depends on the physical properties of the dominant mineral components. Density and, to a lesser extent, seismic velocity show such behaviour. Grain and texture behaviour occur when minor components of the rock are the dominate controls on its physical properties. Grain size and shape control grain properties, and for texture properties the relative positions of these grains are also important. Magnetic and electrical properties behave in this fashion. Thinking in terms of how geological processes change the key characteristics of the major and minor mineralogical components allows the resulting changes in physical properties to be understood and anticipated.     

Joint inversion of long‐offset and central‐loop transient electromagnetic data: Application to a mud volcano exploration in Perekishkul,Azerbaijan

Mud volcanism is commonly observed in Azerbaijan and the surrounding South Caspian Basin. This natural phenomenon is very similar to magmatic volcanoes but differs in one considerable aspect: Magmatic volcanoes are generally the result of ascending molten rock within the Earth's crust, whereas mud volcanoes are characterised by expelling mixtures of water, mud, and gas. The majority of mud volcanoes have been observed on ocean floors or in deep sedimentary basins, such as those found in Azerbaijan. Furthermore, their occurrences in Azerbaijan are generally closely associated with hydrocarbon reservoirs and are therefore of immense economic and geological interest. The broadside long‐offset transient electromagnetic method and the central‐loop transient electromagnetic method were applied to study the inner structure of such mud volcanoes and to determine the depth of a resistive geological formation that is predicted to contain the majority of the hydrocarbon reservoirs in the survey area. One‐dimensional joint inversion of central‐loop and long‐offset transient electromagnetic data was performed using the inversion schemes of Occam and Marquardt. By using the joint inversion models, a subsurface resistivity structure ranging from the surface to a depth of approximately 7 km was determined. Along a profile running perpendicular to the assumed strike direction, lateral resistivity variations could only be determined in the shallow depth range using the transient electromagnetic data. An attempt to resolve further two‐dimensional/three‐dimensional resistivity structures, representing possible mud migration paths at large depths using the long‐offset transient electromagnetic data, failed. Moreover, the joint inversion models led to ambiguous results regarding the depth and resistivity of the hydrocarbon target formation due to poor resolution at great depths (>5 km). Thus, 1D/2D modelling studies were subsequently performed to investigate the influence of the resistive terminating half‐space on the measured long‐offset transient electromagnetic data. The 1D joint inversion models were utilised as starting models for both the 1D and 2D modelling studies. The results tend to show that a resistive terminating half‐space, implying the presence of the target formation, is the favourable geological setting. Furthermore, the 2D modelling study aimed to fit all measured long‐offset transient electromagnetic Ex transients along the profile simultaneously. Consequently, 3125 2D forward calculations were necessary to determine the best‐fit resistivity model. The results are consistent with the 1D inversion, indicating that the data are best described by a resistive terminating half‐space, although the resistivity and depth cannot be determined clearly.     

平面波域反数据处理压制多次波方法研究

在地震勘探领域,尤其是海洋地震勘探中,多次波一直是影响地震处理与解释的主要因素之一.本文基于x-t域的反馈模型及多次波衰减的理论,详尽推导了x-t域反数据处理的方法;基于平面波域多次波的产生机制,借鉴x-t域反数据处理的方法,推导出在平面波域进行反数据处理的原理.文中给出一个有限差分的模拟数据进行测试,处理结果表明,本方法可以有效地衰减表层相关多次波,提高地震数据分析的精确性,在保护了一次波能量的同时,可以更加有效快捷地去除多次波.     

利用相关域小波变换进行SWD资料预处理

随钻地震(SWD)的波场十分复杂,对钻头有效信号和地表机械干扰成分的分析是SWD重要的资料预处理步骤.本文利用有效信号和噪声带有周期性或时延差异等时间结构特征,引入相关域小波变换进行SWD信号检测和分析.有效信号在钻柱内往复多次传播,因而带有严格的周期性,泥浆泵等机械发出的噪声也是周期性的,这些成分在自相关域内可以得到很好的凸显.SWD波场的各种成分,由于到达各个接收道的时延不同,在互相关域的特定时延处也能够得到凸显.利用小波变换对这些在相关域内得到凸显的成分进行多分辨分析,能够获得优势频率范围、周期、衰减等主要特征.根据这些信息,设计出合理的SWD处理方法,初步得到了有效信号的直达波.数据试处理结果表明,相关域小波变换是随钻地震的一个有效的预处理方法.     

The application of spatial derivatives to non‐potential field data interpretation

Source/body edge detection is a common feature in the processing and interpretation of potential field data sets. A wide range of spatial derivatives is available to enhance the information contained in the basic data. Here the ability of these procedures to assist with the mapping interpretation of non‐potential field data is considered. The study uses airborne electromagnetic (conductivity) data but also provides a general context for other conductivity/resistivity data, provided the non‐potential field nature of active and thus spatially‐focused, measurements is acknowledged. The study discusses and demonstrates the application of a range of common spatial derivative procedures, including the analytic signal and upward continuation, to both magnetic and conductivity data. The ability of the tilt derivative to provide enhanced mapping of conductivity data is considered in detail. Tilt and its associated functions are formed by taking combinations of vertical and horizontal derivatives of the data set. Theoretical forward modelling studies are first carried out to assess the performance of the tilt derivative in relation to the detection and definition of concealed conductivity structure. The tilt derivative embodies automatic gain control that normalizes the detection and definition of both weak and strong conductivity gradients across an appropriate subsurface depth range. The use of high‐order spatial derivatives inevitably results in a degree of noise (cultural perturbation) amplification that is survey and technique specific. Both of these aspects are considered using practical case studies of jointly obtained magnetic and conductivity data at a variety of spatial scales.     

A hybrid analysis method for linear dynamic soil-structure interaction in time and frequency domain

IntroductionThe analysis of dynamic soil-structure interaction for important engineering project is still based on linear model (including equivalent linear model) with complex damping, and traditional frequency domain method (Lysmer, et al, 1975, 1981; DING, et al, 1999). Namely, first calculating frequency domain solution by Fourier transform, and then calculating time domain solution by Fourier inverse transform. The motion equation of a system in frequency domain is usually written as (…     

Artificial neural networks for removal of couplings in airborne transient electromagnetic data

Modern airborne transient electromagnetic surveys typically produce datasets of thousands of line kilometres, requiring careful data processing in order to extract as much and as reliable information as possible. When surveys are flown in populated areas, data processing becomes particularly time consuming since the acquired data are contaminated by couplings to man‐made conductors (power lines, fences, pipes, etc.). Coupled soundings must be removed from the dataset prior to inversion, and this is a process that is difficult to automate. The signature of couplings can be both subtle and difficult to describe in mathematical terms, rendering removal of couplings mostly an expensive manual task for an experienced geophysicist. Here, we try to automate the process of removing couplings by means of an artificial neural network. We train an artificial neural network to recognize coupled soundings in manually processed reference data, and we use this network to identify couplings in other data. The approach provides a significant reduction in the time required for data processing since one can directly apply the network to the raw data. We describe the neural network put to use and present the inputs and normalizations required for maximizing its effectiveness. We further demonstrate and assess the training state and performance of the network before finally comparing inversions based on unprocessed data, manually processed data, and artificial neural network automatically processed data. The results show that a well‐trained network can produce high‐quality processing of airborne transient electromagnetic data, which is either ready for inversion or in need of minimal manual processing. We conclude that the use of artificial neural network scan significantly reduce the processing time and its costs by as much as 50%.     

基于神经网络的时间域直升机电磁数据电导率深度成像

采用G-S变换以及高斯数值积分法,形成了时间域直升机的航空电磁响应正演样本集,分析了飞机测量过程中吊舱高度变化对电磁响应的影响,并将吊舱高度的变化等效成电导率为零的假层厚度的变化,以去除高度计等的影响.以假层半空间模型为基础,研究了基于人工神经网络的电导率深度成像算法,通过分析两个三层模型的电导率深度成像结果得出,神经网络方法计算时间域航空电磁探测的视电导率精度较高,特别是对高阻层的视电导率计算.     

复阻尼振动方程频域解法中有关问题的研究

以地震反应谱的形式，阐述了复阻尼系统在频域求解过程中容易出现的问题。通过理论与计算分析得知，这些问题对反应谱计算有很大的影响。另外，还初步分析了时域精确法及频域分析方法的计算效率，认为时域精确法比频域分析方法的计算效率高。     

An operational modal analysis method in frequency and spatial domain

A frequency and spatial domain decomposition method （FSDD） for operational modal analysis （OMA） is presented in this paper, which is an extension of the complex mode indicator function （CMIF） method for experimental modal analysis （EMA）. The theoretical background of the FSDD method is clarified, Singular value decomposition is adopted to separate the signal space from the noise space. Finally, an enhanced power spectrum density （PSD） is proposed to obtain more accurate modal parameters by curve fitting in the frequency domain. Moreover, a simulation case and an application case are used to validate this method.     

OBC与拖缆数据的子波域匹配滤波方法研究

海底电缆（OBC）采集方式的出现,导致OBC采集与拖缆采集数据进行拼接成为三维连片地震资料处理中的一个重要环节.由于采集方式不同,数据在振幅、频率、相位等方面均存在较大差异.为了消除和减小这些差异,为后续的地质解释和油气预测提供便利,本文提出了一种在子波域的1范数匹配滤波方法来计算匹配滤波因子,以用来提高OBC数据和拖缆数据的匹配精度,改善连片处理的数据品质,为后续的地质解释和油气预测提供便利.通过实际资料检验,子波域匹配滤波方法在OBC与拖缆数据拼接应用中取得了较好的效果.     

2D共炮时间域高斯波束偏移

针对传统射线方法在奇异区成像精度不高,而2D频率域高斯波束叠前深度偏移需要计算成像点处每个频率的格林函数,影响计算效率的问题,本文通过使用复走时代替实走时,改变频率域下成像公式的积分顺序,给出了在时间域下进行高斯波束偏移的方法和计算公式.本文使用复杂数值模型验证了2D时间域高斯波束叠前偏移方法的正确性,并同传统射线偏移成像结果做了对比.对比结果表明时间域高斯波束偏移在成像精度上优于传统射线偏移.     

Monitoring and volumetric estimation of injected CO2 using 4D seismic,petrophysical data,core measurements and well logging: a case study at Ketzin,Germany

More than 50 000 tons of CO₂ have been injected at Ketzin into the Stuttgart Formation, a saline aquifer, at approximately 620 m depth, as of summer 2011. We present here results from the 1^st repeat 3D seismic survey that was performed at the site in autumn 2009, after about 22 000 tons of CO₂ had been injected. We show here that rather complex time‐lapse signatures of this CO₂ can be clearly observed within a radius of about 300 m from the injection well. The highly irregular amplitude response within this radius is attributed to the heterogeneity of the injection reservoir. Time delays to a reflection below the injection level are also observed. Petrophysical measurements on core samples and geophysical logging of CO₂ saturation levels allow an estimate of the total amount of CO₂ visible in the seismic data to be made. These estimates are somewhat lower than the actual amount of CO₂ injected at the time of the survey and they are dependent upon the choice of a number of parameters. In spite of some uncertainty, the close agreement between the amount injected and the amount observed is encouraging for quantitative monitoring of a CO₂ storage site using seismic methods.     

三维地形频率域人工源电磁场的边界元模拟方法

提出了一种用边界元法计算频率域人工源三维地形电磁场的数值模拟方法.首先用矢量积分理论和电磁场边界条件,将上半空间（空气）和下半空间（地下介质）两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中,假设单元中电磁场为无限大气空间电磁场与地形影响的叠加,并假设地形影响为常项,这样既保证了计算精度又使得计算方法简便.通过分解和计算,每一个矢量面积分方程分解为对应三个坐标方向的三个常量线性方程,这些线性方程组成了对角占优的线性方程组,可用SSOR方法求解.文中给出了垂直磁偶源的垂直磁场地形影响的例子.