基于非结构网格剖分的海洋可控源电磁三维正则化反演研究

为更好地处理与解释复杂海底地形条件下测得的海洋可控源电磁数据,本文提出了一种基于非结构网格剖分的频率域海洋可控源电磁数据三维正则化反演方法.该方法首先对海洋地电模型以非结构四面体单元进行离散,然后基于矢量有限元方法获得海洋可控源电磁响应和灵敏度信息,最后采用共轭梯度法求解高斯-牛顿反演方程计算模型修正量.为提高反演的稳定性,通过在反演过程中采用对数转换方法实现反演模型参数的上下限约束.本文分别测试了单测线水平海底地形反演算例和面积性测量的起伏海底地形反演算例.反演结果表明,本文提出的频率域海洋可控源电磁三维反演能够准确地恢复高阻储油层的位置和电阻率信息,且计算效率较高,可用于实测海洋电磁资料的处理与解释.     

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

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

2.5维电导率正交各向异性海洋可控源电磁等参有限元数值模拟

本文实现了2.5维电导率正交各向异性海洋可控源电磁等参有限元数值模拟.利用傅里叶变换导出了电导率正交各向异性2.5维海洋可控源电磁法波数域电磁场耦合方程,采用伽里金加权余量法推导了相应的有限元方程;采用任意四边形单元对研究区域进行剖分,在单元中进行双二次插值,将有限元方程化为线性代数方程组;最后,求解线性方程组并进行反傅里叶变换获得空间域电磁场值.这个方法可以模拟海底起伏地形条件下地下任意形状电导率正交各向异性的复杂模型.与一维模型的数值模拟结果对比表明,电磁场数值解与解析解吻合.二维模型的计算结果与二维自适应非结构有限元模拟结果也吻合.水平海底二维地电模型考察了不同各向异性系数对海洋可控源电磁响应的影响特征.海底起伏地形地电模型的数值结果表明,电导率各向异性对海洋可控源电磁响应影响明显,有可能淹没海底地形和高阻油气藏引起的异常.     

利用可探测度探讨海洋可控源电磁响应的等值性

储层特征参数的计算方法简单,能够直观地显示地下电阻率结构的变化,对海洋可控源电磁(CSEM)资料的定性解释非常重要.本文对比分析了不同水深环境下常规归一化振幅、改进的归一化振幅、有效异常和可探测度对于地下油气储层的识别能力.结果表明相较于其他三个参数,可探测度适用于不同水深环境且其分辨能力较高.另外,电磁场响应的等值性是造成反演解释多解性的主要原因之一,研究海洋CSEM法的等值现象有助于资料解释.本文计算了含中间薄层1D模型的可控源电磁场响应,结果表明海洋CSEM法存在不完全的S等值和T等值现象,但不存在H等值现象.在储层特征参数研究的基础上,采用可探测度对海洋地电模型的中间薄层进行识别,有利于解决海洋CSEM法的等值问题,从而提高电磁资料解释的质量.     

一维电阻率各向异性对海洋可控源电磁响应的影响研究

地下介质的电阻率常常表现为各向异性,海底裂隙地层和层状沉积序列可能形成宏观电阻率各向异性.在解释海洋电磁资料时,电阻率各向异性的影响不应该被忽略,否则可能会得到错误的海底地电模型.作者编写了电阻率任意各向异性一维层状介质海洋可控源电磁场计算程序,计算了电阻率各向异性层状模型的海洋可控源电磁响应,讨论了覆盖层和高阻储层分别具有电阻率各向异性时的电磁场响应特征.     

三维任意各向异性介质中海洋可控源电磁法正演研究

由于海底介质受沉积环境的影响,层理发育呈现明显各向异性特征.对于海洋可控源电磁法各向异性的研究以往主要局限于一维和二维模型,为更深入了解复杂情况下海底各向异性对海洋可控源电磁响应的影响规律,本文开展三维任意各向异性介质中海洋可控源电磁法正演研究.采用交错网格有限差分技术,通过对任意各向异性介质电导率张量实行体积和空间电流密度平均,完成海洋可控源电磁二次散射电场的离散化,成功实现任意各向异性介质中海洋可控源电磁正演模拟.通过对几种典型各向异性电性模型条件下海洋电磁电场多分量响应及分布特征和各向同性情况的对比分析,总结电各向异性对海洋电磁响应的影响规律和识别方法.本文算法研究及算例可为海洋可控源电磁数据精细化处理解释提供技术支撑.     

基于三维模拟的海洋CSEM资料处理

海洋可控源电磁法已经成为海洋油气勘探一个重要工具,但是其资料处理和解释还处于定性和一维模拟阶段.在积分方程三维模拟的基础上对Troll油田实测数据进行了处理,采用人机交互三维模拟寻找背景模型和异常体初始模型,最后对异常体电阻率采用准线性近似快速反演,取得了定量的结果.同时,说明对于二维测线和二维模型依然可以用三维来模拟,其结果优于二维反演.在电子计算机技术快速发展的今天,可以预计三维反演将成为资料处理解释的主流.     

层状垂直各向异性介质海洋CSEM发射源及接收站姿态和位置对电磁响应影响研究

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

电磁传播电阻率测井的二维全参数反演方法研究

Baker Hughes INTEQ的MPR是应用较广的传播电阻率测井仪器，它有4条相位差和4条幅度比测井曲线，本文以它在垂直井中的测井响应为例开展电磁传播电阻率资料的二维反演方法研究.由于其测井响应对地层纵向边界的探测灵敏度比对其他地层参数的大得多，在其资料反演中必须重视地层纵向边界位置的提取和反演精度.本文首先应用特征识别方法由分辨率最好的测井曲线提取地层纵向边界，接着对浅、深探测的两条相位差曲线分别做纵向一维反演，利用此结果来建立二维反演的初始地层模型，采用改进的阻尼型正则化高斯 牛顿优化算法，实现了MPR资料的二维全参数反演.对含23个纵向分层的两种典型地层模型的数值模拟测井资料的二维反演得到了令人满意的结果，它可显著提高MPR的资料解释水平.发现只用它的5条测井曲线就可得到与8条测井曲线几乎相同的反演效果.     

陆地可控源电磁法三维勘探的激电效应影响研究

可控源电磁法(CSEM)的正演模拟与反演解释中常忽略激电介质的影响.本文基于电磁三维正反演开源平台ModEM,分别开展了陆地CSEM三维正反演以及考虑激电介质的陆地CSEM三维正演研究,并与一维模型解析解对比验证了三维计算的准确性.数值模拟结果显示,模型含激电效应与不含激电效应的电场E_x频率响应差异明显,激电效应会导致电场E_x幅值减小、相位增加,且对E_x的影响程度与激电异常体的频谱特性及观测的收发距均相关.合成数据的陆地CSEM反演结果显示,数据中含激电效应会使得反演出的异常体在深度、规模和电阻率值上都与真实异常体存在较大偏差,忽略观测数据中的激电效应会明显影响观测的准确性.本文研究表明当陆地CSEM三维观测资料中包含激电效应时,有必要考虑其影响.     

基于高斯牛顿法的频率域可控源电磁三维反演研究

三维反演解释是电磁法勘探发展的重要趋势,而如何提高三维反演的可靠性、稳定性和计算效率是算法开发者们目前的研究重点.本文实现了一种频率域可控源电磁(CSEM)三维反演算法.其中正演基于拟态有限体积法离散化,利用直接矩阵分解技术来求解大型线性系统方程,不仅准确、稳定,而且特别有利于含有大量发射场源位置的CSEM勘探情况;对目标函数的最优化采用高斯牛顿法(GN),具有近似二次的收敛性;使用预条件共轭梯度法(PCG)求解每次GN迭代所得到的法方程,避免了显式求解和存储灵敏度矩阵,减小了计算量.以上这些方法的结合应用,使得本文的三维反演算法准确、稳定且高效.通过陆地和海洋CSEM勘探场景中的典型理论模型的反演测试,验证了本文算法的有效性.     

陆地可控源电磁法探测效果的频率响应

陆地可控源电磁法的观测资料可依据频段范围近似地划分为近区场、中间区场及远区场,但采用测量相互正交电、磁分量,并计算视电阻率的资料处理方式只适用于远区场数据.为更有效地利用陆地可控源电磁法不同区间场的观测资料,本文结合三维数值模拟技术并采用电场分量直接进行反演的策略,对不同区间电场的响应特征与探测效果进行了分析.数值模拟结果表明:近区电场的异常响应最明显,异常响应不随频率发生显著变化,但纵向分辨能力差;远区电场异常响应随频率发生显著变化,其探测深度取决于频率的高低;中间区场较为复杂,地表电场异常响应的等值线中心并不是位于异常体中心正上方,而是在沿场传播方向上向异常体与围岩的分界面处偏移,并且发现中间区场资料的加入会影响反演结果的准确性.综合合成数据和野外实测资料的反演结果,发现结合近区场和远区场资料而舍弃中间区场资料的反演效果更佳,这为陆地可控源电磁法资料的反演解释提供了一种有效途径.     

多道瞬变电磁法共中心点道集数据联合反演

多道瞬变电磁法是目前地面资源勘查领域的研究热点,开展该方法观测数据的一维反演对完成装备研发及技术推广具有积极的推动作用.针对多道瞬变电磁法特殊的源波形和数据采集方式,本文对多道瞬变电磁数据的响应提取和一维反演进行了研究.首先,采用反卷积方法从伪随机二进制序列全波形响应数据中提取大地脉冲响应.为了避免脉冲响应中空气波的干扰,提出采用反向积分算法,把脉冲响应转换成下降沿阶跃响应,并将其整理成共中点道集.最后,为充分利用各个偏移距下的响应数据对不同深度范围目标体的分辨能力,采用基于OCCAM算法的联合反演方法对共中心点道集数据进行联合反演.结果表明在不同目标体深度以及有噪声干扰的情况下,共中心点道集数据联合反演计算均可获得较好的结果.     

基于多震相联合反演江苏地区一维速度模型研究

利用地震走时数据,采用联合反演方法获取了江苏地区的一维P波速度模型。与仅采用初至波走时的传统天然地震走时获取方法相比,该方法充分利用了大量存在的续至波参与反演,能有效改进中下地壳的反演能力。针对地震震相目录中常存在震相标识错误的问题,采用的自动判别筛选震相方法能最大限度提高数据走时的精度,可以对不同震相进行有效区分。与其他常用一维速度模型相比,本文反演的模型对Pg、Pn震相走时拟合效果最佳,残差最小。当所用走时数据拥有较高定位精度时,该反演方法能为研究区三维速度结构成像和地震定位提供较可靠的一维速度模型。     

Large-scale three-dimensional inversion of EarthScope MT data using the integral equation method

In this paper we apply 3D inversion to MT data collected in the Northwestern United States as a part of the EarthScope project. By the end of 2009 MT data had been collected from 262 stations located throughout Oregon, Washington, Idaho, and most of Montana and Wyoming. We used data from 139 MT stations in this analysis. We developed fully parallelized rigorous 3D MT inversion software based on the integral equation method with variable background conductivity. We also implemented a receiver footprint approach which considerably reduced the computational resources needed to invert the large volumes of data covering vast areas. The data set used in the inversion was obtained through the Incorporated Research Institutions for Seismology (IRIS). The inversion domain was divided into 2.7 M cells. The inverted electrical conductivity distribution agrees reasonably well with geological features of the region.     

频率域海洋可控源电磁垂直各向异性三维反演

地层宏观电性各向异性会对可控源电磁响应产生重要影响.由于海底地层电性结构常表现为电导率各向异性,若仅对海洋可控源电磁(MCSEM)数据进行常规各向同性反演,有可能无法获得准确的反演解释结果,从而削弱MCSEM技术的可靠性.本文实现了电导率垂直各向异性(VTI)条件下频率域海洋可控源电磁数据三维反演算法.其中,三维正演采用基于二次场控制方程的交错网格有限体积法,并利用直接矩阵分解技术来求解离散所得的大型线性方程组,有利于快速计算多场源的响应.反演采用具有近似二次收敛性的高斯牛顿算法对目标函数进行最优化.最后,对具有VTI电性各向异性特征的盐丘构造模型的MCSEM合成数据分别进行了电导率各向同性和垂直各向异性三维反演,结果表明:各向同性三维反演算法无法对受VTI介质影响的MCSEM数据进行正确的反演解释,而垂直各向异性三维反演能够获得更为可靠的地下电阻率结构和异常体分布,展现出对海底电性各向异性结构更为优良的反演解释能力.     

Parameter optimization of gravity density inversion based on correlation searching and the golden section algorithm

For density inversion of gravity anomaly data, once the inversion method is determined, the main factors affecting the inversion result are the inversion parameters and subdivision scheme. A set of reasonable inversion parameters and subdivision scheme can, not only improve the inversion process efficiency, but also ensure inversion result accuracy. The gravity inversion method based on correlation searching and the golden section algorithm is an effective potential field inversion method. It can be used to invert 2D and 3D physical properties with potential data observed on flat or rough surfaces. In this paper, we introduce in detail the density inversion principles based on correlation searching and the golden section algorithm. Considering that the gold section algorithm is not globally optimized, we present a heuristic method to ensure the inversion result is globally optimized. With a series of model tests, we systematically compare and analyze the inversion result efficiency and accuracy with different parameters. Based on the model test results, we conclude the selection principles for each inversion parameter with which the inversion accuracy can be obviously improved.     

基于拟态有限体积法的频率域可控源三维正演计算

大规模地球物理电磁数据的定量解释需要发展高效、稳定的三维正反演算法.本文通过求解离散化的三维电场矢量Helmholtz方程,实现了基于有限体积法的频率域可控源电磁(CSEM)三维正演算法.为模拟具有强电性差异的三维电性介质,该算法采用拟态有限体积法(MFV)对Maxwell方程组进行离散化;另外,为获得稳定、高精度的正演数值结果,采用直接矩阵分解技术来求解离散所得到的大型稀疏线性方程组.对于具有多个发射源的CSEM测量来说,一次矩阵分解结果能够用于同频率下所有场源的正演计算.为降低场源奇异性及边界条件对数值精度的影响,采用虚拟场源校正技术,避免了散射场公式中在构建场源项时所需的大量时间.对于具有多个频率的CSEM的模拟计算,采用分频并行策略来加快三维正演计算.最后,通过与一维层状模型及三维模型的数值结果的对比验证了本文所开发的正演算法对频率域CSEM模拟计算的准确性及有效性,表明该正演算法能够有效应用于三维介质的数值计算.另外,对于多频率CSEM的并行测试结果表明基于分频并行策略的并行计算能够显著地降低正演计算时间.