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

2.
The theory behind transient electromagnetic surveys can be well described in terms of transverse magnetic and transverse electric modes. Soundings using transverse magnetic and transverse electric modes require different source configurations. In this study, we consider an alternating transverse magnetic field excitation by a circular electric dipole. The circular electric dipole transmitter is a horizontal analogue of the vertical electric dipole. Offshore surveys using circular electric dipole might represent an alternative to the conventional marine controlled‐source electromagnetic method at shallow sea and/or for exploring relatively small targets. Field acquisition is carried out by recording either electric or magnetic responses. Electric responses bear information on the 1D structure of a layered earth and successfully resolve high‐resistivity targets in marine surveys. Land‐based circular electric dipole soundings are affected by induced polarisation. On the contrary, magnetic responses are absent on the surface of a 1D earth, and as a result, they are very sensitive to any and even very small 3D conductivity perturbations. In addition, they are sensitive to induced polarisation or some other polarisation effects in the subsurface. At present, circular electric dipole transmitters and magnetic receivers are successfully used in on‐land mineral and petroleum exploration.  相似文献   

3.
罗鸣  裴建新  徐震寰 《地球物理学报》2017,60(12):4901-4915
众所周知,电阻率各向异性对海洋可控源电磁(CSEM)响应有显著影响.在模拟和解释海洋可控源电磁资料时,通常假设接收站布放于平稳的海底,发射源以理想的水平电偶极源形式按预设路径拖曳前进.然而,在实际海洋作业中,自由下沉的接收站可能沉降在局部倾斜的海底面上,在海底的方位可能是任意取向的.同时,由于洋流等因素的影响,发射源偶极子会发生旋转、倾斜,其位置也会偏离预设的拖曳路径,所以发射源和接收站姿态的变化都将对电磁资料产生影响.本文计算了电阻率各向异性介质中,发射源任意姿态任意位置激发、接收站任意姿态任意位置接收的海洋CSEM响应,分析观测系统引起数据误差的机制,分别在电阻率各向同性、围岩电阻率各向异性和高阻储层电阻率各向异性三种模型中,讨论发射源和接收站的姿态及位置变化对响应的影响.计算结果表明,电磁场振幅误差主要来源于发射源和接收站的倾斜角,而相位的变化为发射源位置和接收站倾斜角、位置共同影响的结果.  相似文献   

4.
This paper presents a new explicit method for the estimation of layered vertical transverse isotropic (VTI) anisotropic parameters from walkaway VSP data. This method is based on Dix‐type normal moveout (NMO) inversion. To estimate interval anisotropic parameters above a receiver array, the method uses time arrivals of surface‐related double‐reflected downgoing waves. A three‐term NMO approximation function is used to estimate NMO velocity and a non‐hyperbolic parameter. Assuming the vertical velocity is known from zero‐offset VSP data, Dix‐type inversion is applied to estimate the layered Thomsen anisotropic parameters ?, δ above the receivers array. Model results show reasonable accuracy for estimates through Dix‐type inversion. Results also show that in many cases we can neglect the influence of the velocity gradient on anisotropy estimates. First breaks are used to estimate anisotropic parameters within the walkaway receiver interval. Analytical uncertainty analysis is performed to NMO parameter estimates. Its conclusions are confirmed by modelling.  相似文献   

5.
Recently, new on‐shore acquisition designs have been presented with multi‐component sensors deployed in the shallow sub‐surface (20 m–60 m). Virtual source redatuming has been proposed for these data to compensate for surface statics and to enhance survey repeatability. In this paper, we investigate the feasibility of replacing the correlation‐based formalism that undergirds virtual source redatuming with multi‐dimensional deconvolution, offering various advantages such as the elimination of free‐surface multiples and the potential to improve virtual source repeatability. To allow for data‐driven calibration of the sensors and to improve robustness in cases with poor sensor spacing in the shallow sub‐surface (resulting in a relatively high wavenumber content), we propose a new workflow for this configuration. We assume a dense source sampling and target signals that arrive at near‐vertical propagation angles. First, the data are preconditioned by applying synthetic‐aperture‐source filters in the common receiver domain. Virtual source redatuming is carried out for the multi‐component recordings individually, followed by an intermediate deconvolution step. After this specific pre‐processing, we show that the downgoing and upgoing constituents of the wavefields can be separated without knowledge of the medium parameters, the source wavelet, or sensor characteristics. As a final step, free‐surface multiples can be eliminated by multi‐dimensional deconvolution of the upgoing fields with the downgoing fields.  相似文献   

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

7.
The time-domain controlled source electromagnetic method is a geophysical prospecting tool applied to image the subsurface resistivity distribution on land and in the marine environment. In its most general set-up, a square-wave current is fed into a grounded horizontal electric dipole, and several electric and magnetic field receivers at defined offsets to the imposed current measure the electromagnetic response of the Earth. In the marine environment, the application often uses only inline electric field receivers that, for a 50% duty-cycle current waveform, include both step-on and step-off signals. Here, forward and inverse 1D modelling is used to demonstrate limited sensitivity towards shallow resistive layers in the step-off electric field when transmitter and receivers are surrounded by conductive seawater. This observation is explained by a masking effect of the direct current signal that flows through the seawater and primarily affects step-off data. During a step-off measurement, this direct current is orders of magnitude larger than the inductive response at early and intermediate times, limiting the step-off sensitivity towards shallow resistive layers in the seafloor. Step-on data measure the resistive layer at times preceding the arrival of the direct current signal leading to higher sensitivity compared to step-off data. Such dichotomous behaviour between step-on and step-off data is less obvious in onshore experiments due to the lack of a strong overlying conductive zone and corresponding masking effect from direct current flow. Supported by synthetic 1D inversion studies, we conclude that time-domain controlled source electromagnetic measurements on land should apply both step-on and step-off data in a combined inversion approach to maximize signal-to-noise ratios and utilize the sensitivity characteristics of each signal. In an isotropic marine environment, step-off electric fields have inferior sensitivity towards shallow resistive layers compared to step-on data, resulting in an increase of non-uniqueness when interpreting step-off data in a single or combined inversion.  相似文献   

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

9.
The receiver function method was originally developed to analyse earthquake data recorded by multicomponent (3C) sensors and consists in deconvolving the horizontal component by the vertical component. The deconvolution process removes travel path effects from the source to the base of the target as well as the earthquake source signature. In addition, it provides the possibility of separating the emergent P and PS waves based on adaptive subtraction between recorded components if plane waves of constant ray parameters are considered. The resulting receiver function signal is the local PS wave's impulse response generated at impedance contrasts below the 3C receiver.We propose to adapt this technique to the wide‐angle multi‐component reflection acquisition geometry. We focus on the simplest case of land data reflection acquisition. Our adapted version of the receiver function approach consists in a multi‐step procedure that first removes the P wavefield recorded on the horizontal component and next removes the source signature. The separation step is performed in the τ?p domain while the source designature can be achieved in either the τ?p or the t?x domain. Our technique does not require any a priori knowledge of the subsurface. The resulting receiver function is a pure PS‐wave reflectivity response, which can be used for amplitude versus slowness or offset analysis. Stack of the receiver function leads to a high‐quality S wave image.  相似文献   

10.
罗鸣  李予国 《地球物理学报》2015,58(8):2851-2861
地下介质的电阻率常常表现为各向异性,海底裂隙地层和层状沉积序列可能形成宏观电阻率各向异性.在解释海洋电磁资料时,电阻率各向异性的影响不应该被忽略,否则可能会得到错误的海底地电模型.作者编写了电阻率任意各向异性一维层状介质海洋可控源电磁场计算程序,计算了电阻率各向异性层状模型的海洋可控源电磁响应,讨论了覆盖层和高阻储层分别具有电阻率各向异性时的电磁场响应特征.  相似文献   

11.
We developed a new marine controlled‐source electromagnetic receiver for detecting methane hydrate zones and oil and gas reservoirs on the seafloor, which is not imaged well by seismic reflection surveys. To determine the seafloor structure, the electromagnetic receiver should have low noise, power consumption, clock drift error, and operating costs while being highly reliable. Because no suitable receiver was available in our laboratory, we developed a new marine controlled‐source electromagnetic receiver with these characteristics; the receiver is equipped with acoustic telemetry modem and an arm‐folding mechanism to facilitate deployment and recovering operations. To demonstrate the applicability of our new receiver, we carried out a field experiment offshore of Guangzhou in the South China Sea, where methane hydrates have been discovered. We successfully obtained controlled‐source electromagnetic data along a profile about 13 km long. All six new receivers were recovered, and high‐quality electromagnetic data were obtained. Relatively high apparent resistivity values were detected. The results of the offshore field experiment support the claim that the electromagnetic data obtained using the new receiver are of sufficient quality for the survey target.  相似文献   

12.
A comprehensive controlled source electromagnetic (CSEM) modelling study, based on complex resistivity structures in a deep marine geological setting, is conducted. The study demonstrates the effects of acquisition parameters and multi‐layered resistors on CSEM responses. Three‐dimensional (3D) finite difference time domain (FDTD) grid‐modelling is used for CSEM sensitivity analysis. Interpolation of normalized CSEM responses provides attributes representing relative sensitivity of the modelled structures. Modelling results show that fine grid, 1 × 1 km receiver spacing, provides good correlations between CSEM responses and the modelled structures, irrespective of source orientation. The resolution of CSEM attributes decreases for receiver spacing >2 × 2 km, when using only in‐line data. Broadside data in the grid geometry increase data density by 100 – approximately 200% by filling in in‐line responses and improve the resolution of CSEM attributes. Optimized source orientation (i.e., oblique to the strike of an elongated resistor) improves the structural definition of the CSEM anomalies for coarse‐grid geometries (receiver spacing ≥3 × 3 km). The study also shows that a multi‐resistor anomaly is not simply the summation but a cumulative response with mutual interference between constituent resistors. The combined response of constituent resistors is approximately 50% higher than the cumulative response of the multi‐resistor for 0.5 Hz at 4000 m offset. A gradual inverse variation of offset and frequency allows differentiation of CSEM anomalies for multi‐layered resistors. Similar frequency‐offset variations for laterally persistent high‐resistivity facies show visual continuity with varying geometric expressions. 3D grid‐modelling is an effective and adequate tool for CSEM survey design and sensitivity analysis.  相似文献   

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

14.
本文针对多辐射场源地空瞬变电磁法理论,建立了相应的多分量全域视电阻率定义方法.分析了利用磁场强度进行多辐射场源地空系统全域视电阻率定义的优点,针对磁场强度的各个分量,提出了各自的全域视电阻率算法,实现了多分量、全时域、全空域视电阻率计算,并分析了偏移距对全域视电阻率的影响.通过调整源的相对位置及电流方向等参数,多辐射源瞬变电磁地空系统不仅可以加强不同分量信号强度,削弱随机干扰,还可更好地分辨地下异常体的位置.通过对所设计模型的处理,证实了多辐射场源地空系统多分量全域视电阻率算法的有效性,也验证了多辐射场源地空系统的优势.  相似文献   

15.
To investigate the characteristics of the anisotropic stratum, a multi‐azimuth seismic refraction technique is proposed in this study since the travel time anomaly of the refraction wave induced by this anisotropic stratum will be large for a far offset receiver. To simplify the problem, a two‐layer (isotropy–horizontal transverse isotropy) model is considered. A new travel time equation of the refracted P‐wave propagation in this two‐layer model is derived, which is the function of the phase and group velocities of the horizontal transverse isotropic stratum. In addition, the measured refraction wave velocity in the physical model experiment is the group velocity. The isotropic intercept time equation of a refraction wave can be directly used to estimate the thickness of the top (isotropic) layer of the two‐layer model because the contrast between the phase and group velocities of the horizontal transverse isotropic medium is seldom greater than 10% in the Earth. If the contrast between the phase and group velocities of an anisotropic medium is small, the approximated travel time equation of a refraction wave is obtained. This equation is only dependent on the group velocity of the horizontal transverse isotropic stratum. The elastic constants A11, A13, and A33 and the Thomsen anisotropic parameter ε of the horizontal transverse isotropic stratum can be estimated using this multi‐azimuth seismic refraction technique. Furthermore, under a condition of weak anisotropy, the Thomsen anisotropic parameter δ of the horizontal transverse isotropic stratum can be estimated by this technique as well.  相似文献   

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

17.
Using approximate boundary conditions, expressions for electromagnetic fields have been derived for a thin, highly resistive layer lying between two homogeneous layers excited by an electric dipole grounded on the surface of the earth. The variations of the fields with the parameter T/T1 (ratio of the transverse resistance of the thin layer to the transverse resistance of the first layer) were studied in relation to frequency, time, the normalized separation source—receiver, and the angle between the source and the radius to the observation point. For a value of h2/h1 (ratio of thickness of second layer to the thickness of the first layer) approximately equal to 0.2, the general three-layer medium case gives the same results as this approach. It was found that the electric fields have a very strong dependence on the parameter T (transverse resistance) which characterizes the thin, highly resistive layer. However, the magnetic fields depend only very weakly on this parameter.  相似文献   

18.
Elastic imaging from ocean bottom cable (OBC) data can be challenging because it requires the prior estimation of both compressional‐wave (P‐wave) and shear‐wave (S‐wave) velocity fields. Seismic interferometry is an attractive technique for processing OBC data because it performs model‐independent redatuming; retrieving ‘pseudo‐sources’ at positions of the receivers. The purpose of this study is to investigate multicomponent applications of interferometry for processing OBC data. This translates into using interferometry to retrieve pseudo‐source data on the sea‐bed not only for multiple suppression but for obtaining P‐, converted P to S‐wave (PS‐wave) and possibly pure mode S‐waves. We discuss scattering‐based, elastic interferometry with synthetic and field OBC datasets. Conventional and scattering‐based interferometry integrands computed from a synthetic are compared to show that the latter yields little anti‐causal response. A four‐component (4C) pseudo‐source response retrieves pure‐mode S‐reflections as well at P‐ and PS‐reflections. Pseudo‐source responses observed in OBC data are related to P‐wave conversions at the seabed rather than to true horizontal or vertical point forces. From a Gulf of Mexico OBC data set, diagonal components from a nine‐component pseudo‐source response demonstrate that the P‐wave to S‐wave velocity ratio (VP/VS) at the sea‐bed is an important factor in the conversion of P to S for obtaining the pure‐mode S‐wave reflections.  相似文献   

19.
In the field of seismic interferometry, researchers have retrieved surface waves and body waves by cross‐correlating recordings of uncorrelated noise sources to extract useful subsurface information. The retrieved wavefields in most applications are between receivers. When the positions of the noise sources are known, inter‐source interferometry can be applied to retrieve the wavefields between sources, thus turning sources into virtual receivers. Previous applications of this form of interferometry assume impulsive point sources or transient sources with similar signatures. We investigate the requirements of applying inter‐source seismic interferometry using non‐transient noise sources with known positions to retrieve reflection responses at those positions and show the results using synthetic drilling noise as source. We show that, if pilot signals (estimates of the drill‐bit signals) are not available, it is required that the drill‐bit signals are the same and that the phases of the virtual reflections at drill‐bit positions can be retrieved by deconvolution interferometry or by cross‐coherence interferometry. Further, for this case, classic interferometry by cross‐correlation can be used if the source power spectrum can be estimated. If pilot signals are available, virtual reflection responses can be obtained by first using standard seismic‐while‐drilling processing techniques such as pilot cross‐correlation and pilot deconvolution to remove the drill‐bit signatures in the data and then applying cross‐correlation interferometry. Therefore, provided that pilot signals are reliable, drill‐bit data can be redatumed from surface to borehole depths using this inter‐source interferometry approach without any velocity information of the medium, and we show that a well‐positioned image below the borehole can be obtained using interferometrically redatumed reflection responses with just a simple velocity model. We discuss some of the practical hurdles that restrict the application of the proposed method offshore.  相似文献   

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

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

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