一维垂直各向异性介质频率域海洋可控源电磁资料反演方法

本文提出了一维垂直各向异性(VTI)介质倾斜偶极源频率域海洋可控源电磁(CSEM)资料高斯-牛顿反演方法.在电阻率各向异性介质水平偶极源和垂直偶极源海洋CSEM正演算法的基础上,利用欧拉旋转方法,实现了各向异性介质倾斜偶极源海洋CSEM正演算法.海洋可控源电磁场关于地下介质横向电阻率(ρ_h)和垂向电阻率(ρ_v)的偏导数(即灵敏度矩阵)是解析计算的,结合垂直各向异性介质横向电阻率与垂向电阻率的关系,将各向异性率融入到正则化因子选择中,实现了正则化因子的自适应选择.理论模型合成数据和实测资料反演算例表明,我们提出的反演方法能够较准确的重构海底围岩和基岩的各向异性电阻率以及高阻薄层的埋藏深度、厚度和垂向电阻率.

Frequency-domain inversion of marine CSEM data in one-dimensional vertically anisotropic structures

The marine controlled-source electromagnetic (CSEM) method has proven to be an important addition to seismic imaging techniques in exploration of offshore hydrocarbon reservoirs and near-surface investigations. In inversion and interpretation of industrial CSEM data sets, this method can be used in a number of situations by assuming isotropy. However, the presence of electric anisotropy in the Earth's crust, due to thin layer interbedding or grain alignments in the sediments, can significantly alter the response measured by the EM receivers. Ignoring anisotropy in interpreting marine CSEM data may lead to distorted images of seabed conductivity structures, even misinterpretation.#br#In this paper, we present an inversion method for frequency domain marine controlled-source electromagnetic data generated by a titled dipole source in vertically anisotropic stratified media. This approach is based on the Gauss-Newton scheme. We extend and generalize the formulation of marine CSEM fields to calculation of the electromagnetic fields excited by arbitrarily oriented dipole sources. The partial derivatives of the electromagnetic fields with respect to both the horizontal and vertical resistivity are analytically calculated. Based on the relationship between the horizontal resistivity (ρ_h) and the vertical resistivity (ρ_v) of the inversion model, an adaptive selection method for regularization factors is proposed to balance the effects of the data misfit and the structural constraint.#br#The synthetic and real data inversion tests indicate that our inversion method can reconstruct the anisotropic resistivity of the overburden layer and the basement, and the burial depth, the thickness and the vertical resistivity of the reservoir layer can be well recovered. A better reconstruction can be obtained with multiple frequency and multiple component data sets. Combing inline and broadside geometry data sets can provide high resolution in reconstructing the burial depth of the reservoir layer and resistivity of the anisotropic basement.