基于深度加权的重力梯度张量数据的3D聚焦反演

针对重力梯度数据聚焦反演结果中存在的"上漂"现象,在经典Tikhonov正则化理论框架下,引入最小支撑泛函数对反演模型进行约束以避免反问题解的不稳定,并针对重力梯度数据聚焦反演中存在的趋肤效应,在模型目标函数中引入指数深度加权函数。通过理论模型,对部分重力梯度张量分量进行了单独以及联合聚焦反演,验证了基于深度加权的聚焦反演方法的有效性,并将该反演方法运用到涩北一号气田区的实际数据中,反演结果较好地反映出气田位置。     

基于粒子群算法的重力张量反演研究

重力梯度张量的定义是对重力位求二阶导数,相比于传统的布格重力异常,它在反映地下密度异常分布上有着更高的灵敏度,且能够进一步准确直接地反映目标体的边界.但对单个张量分量作反演时,可能会丢失一部份有用的信息,从而造成反演结果的误差.而全张量反演是将重力的五个梯度张量联合起来进行反演,这样做可以综合更丰富的场源信息.相比于传统做法中的布格重力异常反演和单重力张量分量反演,得到的反演结果不仅有了更高的分辨率,在识别目标体特征上也有更好的效果.粒子群算法是一种基于群体智能的优化迭代算法,这里利用粒子群算法对重力张量单分量、布格重力异常和全张量分别进行反演,并对结果进行简要分析.     

三维重力梯度局部光滑约束反演

针对重力勘探中光滑反演存在的分辨率较低的问题,本文提出一种基于地质体埋深、地层倾向等一定先验信息的局部光滑约束三维反演算法,并提出了一种光滑反演中粗糙度矩阵的存储方式,该存储方式可以将M×N维粗糙度矩阵存储为M×2维,减少了计算机计算内存,且详细介绍了该存储方式下粗糙度矩阵与其他矩阵相乘时,粗糙度矩阵所存储的位置信息的读取方式以及与其他矩阵逐列逐步相乘最终得到计算结果的过程。最后,利用文中提出的算法对理论模型和实测数据进行反演试算,结果表明局部光滑反演算法相比于全局光滑反演结果更加准确,且该算法在一定噪声水平下依然稳定,在实际生产中有效可行。     

基于断层擦痕侧伏角公式的构造应力张量反演

利用Aleksandrowski推导的断层擦痕侧伏角公式,提出了一种改进的构造应力张量反演法。根据不同应力比值下构造应力与断层擦痕侧伏角的关系,利用最小二乘法,详细推导了构造应力张量的反演方法,最后利用Matlab软件编制程序实现了该计算过程。在此基础上,现场调查了江西武山铜矿北矿带的断层滑动数据,利用上述改进方法反演了该地区的构造应力张量方向特征。将上述反演结果与由震源机制解获得的最新构造应力场结果进行对比分析,并结合研究区断层的展布特征及所在的地质构造背景,获得了该地区自晚侏罗世以来的构造应力场特征,表明该地区构造应力场具有一定的稳定性。改进后的构造应力张量求解过程及工程应用表明,将Aleksandrowski的图示法和Etchecopar等的反演法结合起来反演构造应力张量,是对已有方法的改进,具有更加简单、实现容易等特点,计算结果比较符合实际。     

基于非结构化网格的瞬变电磁2．5维有限元正演模拟

利用Delaunay三角化这种网格非结构化方法。通过编程实现了二维模型的非结构化三角形网格剖分,并编写了中心回线法瞬变电磁2．5维有限元正演程序。与前人计算结果对比,在取得相同计算精度的情况下,与结构化网格相比,非结构化网格所需网格和节点数量大大减少,计算效率更高。通过将非结构化网格法引入到瞬变电磁2．5维正演模拟中,实现了对复杂二维地电模型的有限元数值模拟,提高了现有有限元算法的应用范围。     

油气模型的重力梯度张量研究

以位场数据处理中的傅里叶变换为出发点,推导出由重力异常计算梯度张量的余弦谱变换公式,进而通过离散余弦变换实现梯度张量的数值计算。对不同贮油厚度的三度体球冠模型梯度张量的研究表明:当油气含量较小时,其梯度张量与均匀三度体球冠模型所表现出的特征相似,V_xx、V_yy呈“两高夹一低”的对称特征;V_xy总体上呈现出向中心变小的趋势,在4个边缘的中部分别有4个呈两两对称的相对较高值;V_zz四周低,中间高;V_zx、V_zy呈现出西低东高和南低北高的特征。而当贮油厚度逐渐增大时,V_xx、V_yy的形态则渐变为“两低夹一高”;V_xy在四周有相对低值,中部有极大值;V_zz四周高,中间低,变化更加尖锐;V_zx及V_zy呈西高东低和南高北低的特征。对三度体球冠模型重力梯度张量的研究,可以为寻找储油背斜构造提供理论依据。     

基于非结构三角网格的海洋CSEM和MT二维联合反演研究

采用非结构三角形网格结合快速Occam算法对海洋可控源电磁(CSEM)与大地电磁(MT)数据进行二维联合反演试验研究.非结构三角形网格能够准确地模拟起伏地形和复杂地质构造,对反演目标区域采用精细网格剖分,其他区域采用粗网格剖分,在满足精度的前提下减少了不必要的计算量.将CSEM与MT数据加入到同一反演数据集中,通过对联合反演数据权重公式进行推导,构建CSEM及MT数据相关权重因子,控制不同数据的拟合权重来实现联合反演.最后对不同模型进行了反演计算,结果表明联合反演比单一反演对海底构造与异常体的还原度更高,验证了联合反演算法的可靠性.     

基于非结构化网格的2．5D直流电阻率模拟

结合大型通用有限元软件ANSYS所提供的强大网格剖分功能,实现了2D地质模型的非结构化三角形网格剖分,并编写了2.5D直流电阻率有限元法FORTRAN语言计算程序.经计算对比表明,在节点数基本相同时,非结构化三角形网格比传统的矩形网格、矩形～三角形网格具有更高的计算精度,可以很好地拟合复杂地形和地质体边界.在此基础上,采用非结构化三角形网格,计算了起伏地形下2.5D地电模型的视电阻率异常,并利用比较法进行了地形改正.与水平地形时的结果对比表明,比较法可以较好地消除地形影响,突出局部地质体的异常.     

基于非结构网格的带地形MT二维Occam反演及应用

本文对基于自适应非结构三角形网格的带地形MT数据二维Occam反演进行应用研究.自适应非结构三角形网格能够准确地模拟起伏地形和复杂地质构造,正演网格由有限元解的后验误差估计指导自动细化,保证了模型响应的准确性,反演网格在反演目标区域采用精细网格剖分,在模型边界区域采用粗网格剖分,在满足反演精度的前提下减少了不必要的反演...     

基于非结构化网格的大规模磁测数据正则化反演研究

磁测数据解析正演的内存需求与模型规模、观测点数成正比,以其为依托开展大规模磁测数据反演时,计算成本较大。文章利用非结构化网格有限单元法进行磁测正演,分析了地表观测网格、边界节点数和局部网格密度等因素对正演计算精度的影响;在此基础上,基于灵敏度矩阵隐式存储的高斯-牛顿法求解正则化反演目标函数,实现了大规模磁测数据三维反演。理论模型试算和铀矿区域实测数据案例表明,文章所述方法可保证正演计算的精确度,提高了计算效率、减少了反演的内存消耗,可用于复杂条件下大规模磁测数据的三维反演。     

重力梯度数据协克里金三维反演确定岩脉倾向

采用协克里金法处理重力或重力梯度数据反演地下3D密度时,岩脉倾角的先验信息是否合理决定了反演结果的好坏。为了使协克里金法在先验信息不充足时也能得到理想的反演效果,我们提出角度扫描的方法来获取岩脉倾角信息:在保持其他反演参数不变的前提下,假设模型的倾斜角度未知,每次引入不同角度进行协克里金反演,再根据先验信息剔除不合理密度值,得到密度分布;统计预测数据与测量数据拟合残差的标准差,根据极小值对应的角度估计岩脉的倾角。本文基于4种不同倾斜角度(0°,45°,90°,135°)岩脉模型的重力梯度垂直分量Tzz验证方法有效性,都能够准确地估计出岩脉倾斜角度。这使得协克里金法在岩脉倾角未知的情况下,也能得到较好的反演结果。     

利用重力梯度张量研究黑龙江省虎林盆地的断裂构造

以重力异常为依据,利用快速傅里叶变换计算了虎林盆地的重力全梯度张量。根据Vxx和Vyy的零值线能突出场源 体的边界以及走向近似于x或y轴方向的地质体、Vzz能够圈定异常体的位置和范围以及Vzx,Vzy和Vxy能够突出显示NS, EW和NE（NW）向的断裂构造线,分析虎林盆地重力异常及其梯度张量展布特征。推断研究区存在近NE,NW,NEE,近 EW向和NS向断裂,对其中5条深大断裂进行了描述。根据重力异常及其全梯度张量和断裂展布特征,可以将虎林盆地划 分为3个隆起和2个凹陷,与迎春隆起、虎林河凹陷、虎林中央隆起、东林子凹陷和兴凯湖隆起5个构造单元存在很好的对 应关系。     

基于Zoeppritz方程的纵横波模量反演

叠前地震反演的主要目的是根据地震数据精确地求取地下介质的弹性参数,进而进行流体识别和油藏描述.常规叠前反演方法使用Zoeppritz方程近似式计算反射系数,但是含有丰富信息的大偏移距地震数据无法被有效利用,针对这个问题,将纵横波模量引入Zoeppritz方程,建立了基于Zoeppritz方程的纵横波模量广义线性反演方法,适用于大角度地震数据.模型测试和实际地震资料应用表明,此方法能够可靠的得到地下介质的纵横波模量,相对于间接反演减小了间接反演带来的累积误差,为后续流体识别和含油气检测提供了很好的数据基础.     

重力异常及其梯度张量DEXP定量解释方法的影响因素分析

快速成像反演方法是近几年重磁勘探定量解释的一个发展热点,由于其在计算过程中不需要加入先验信息,故而可以较为快速地估算场源的深度与密度等相关参量。DEXP(depth from extreme point)成像法由于在深度加权函数中考虑了构造指数,使得成像结果更为准确。笔者基于DEXP快速成像法的基本理论,首先,推导了重力异常及重力梯度张量的DEXP变换形式,并将其应用到场源解释之中;然后,采用理论模型试验,分别分析了数据的点距、误差、计算范围以及背景场对成像结果的影响;最后,将该成像方法应用于实测的全张量重力梯度数据,并与前人研究结果进行了对比与分析。理论模型试验与实际应用均表明:DEXP成像法能够有效压制数据噪声的影响,具有计算稳定性和准确性特点;数据的点距、计算范围和背景场均对DEXP成像结果具有一定的影响,因此在实际数据处理时,应该综合考虑它们对成像结果的影响,并且需要进行相关数据预处理以提高定量解释的精度。     

张量CSAMT数据的二维反演

在一单行本中,我们提出了一种ＣＳＡＭＴ正演模拟算法并说明了这种算法在解释野外数据中的应用效果,在本文中,提出了一种用于二维电阻率模型的ＣＳＡＭＴ数据反演算法,这种算法包含一个三维源并可以反演（ＴＥ）模式及（ＴＭ）模式下的数据,用这种算法对假设的数据进行试算,结果表明,不用做平面波校正就能对远区场,过渡区场及近区场进行反演,最终的模型结果显示,原始电阻率模型可被恢复,而不存在任何由上源于效应而产生的     

基于OpenMP的重力张量并行正演

重力张量是重力位的二阶空间导数,对密度体的变化和细节部份反映更为灵敏.对于复杂的重力密度模型来说,张量的解析公式是很难推导的.为了模拟这种情况,将复杂模型进行有限元剖分,计算每个单元的重力异常对测点张量的影响,最后叠加得到整个复杂模型在测点处的重力张量.采用Delaunay四面体非结构化网格对密度体进行剖分,经分析表明,这种剖分方式具有较高的精度.有限元剖分的方式拟合复杂形体缺点之一是计算量比较大,而并行化能有效解决这个问题.这里基于OpenMP并行模型编写了重力张量并行正演程序,分析了不同情况下并行的执行性能,为大规模复杂模型的重力张量正演提供了一种并行策略和思路.     

基于OpenMP的重力张量并行正演CSCD

重力张量是重力位的二阶空间导数,对密度体的变化和细节部份反映更为灵敏。对于复杂的重力密度模型来说,张量的解析公式是很难推导的。为了模拟这种情况,将复杂模型进行有限元剖分,计算每个单元的重力异常对测点张量的影响,最后叠加得到整个复杂模型在测点处的重力张量。采用Delaunay四面体非结构化网格对密度体进行剖分,经分析表明,这种剖分方式具有较高的精度。有限元剖分的方式拟合复杂形体缺点之一是计算量比较大,而并行化能有效解决这个问题。这里基于OpenMP并行模型编写了重力张量并行正演程序,分析了不同情况下并行的执行性能,为大规模复杂模型的重力张量正演提供了一种并行策略和思路。     

网格剖分对大地电磁测深反演精度的影响

网格剖分的合适与否,影响着反演结果的精度。用不同模型、不同反演方法,对比验证粗细两种不同的网格对反演精度的影响。研究结果表明,从整体上看两种不同网格剖分方式,粗网格在深部比细网格模拟精度高,但是近地表开始阶段,粗网格与细网格波动幅度都较大,都脱离正常值,粗网格模拟精度不如细网格;粗网格的整体变化幅度比细网格缓和,其中对低阻异常体反演精度较好,但对高阻异常体反演精度较差,与真实值相去甚远。     

基于深度学习的重力梯度事后误差补偿方法

旋转式重力梯度仪动态测量的误差补偿是航空重力梯度测量的关键技术.因深度学习直接提取数据特征,有望补偿以往通过建模难以补偿的误差,从而进一步提高误差补偿效果,笔者提出了基于深度学习的航空重力梯度测量事后误差补偿方法,利用航空试验实测数据建立数据集,搭建并训练由梯度仪运动参数映射到梯度仪输出噪声的神经网络,利用该网络预测梯...     

基于非结构化有限元的三维井地电阻率法约束反演

电磁探测反演是典型的不适定问题,易造成反演结果的多解性,不适定性是反演自身固有的特征,没有求解的附加信息这一本质困难是很难克服的,解决该问题的有效方法是研究约束反演。本文采用目前较为主流的高斯牛顿—共轭梯度法(GN-CG),在反演目标函数中直接施加约束条件,将介质电阻率的取值范围作为先验信息和约束条件以外点罚函数法的方式引入到反演目标函数中,与常规三维电阻率反演目标函数相比,增加了不等式约束项的目标函数,理论上可以压制反演的多解性。通过多种理论模型的测试结果表明,本文基于不等式约束的三维井地电阻率反演算法有效地改善了反演结果的精度,以惩罚函数法施加不等式约束条件的方式是现实可行及有效的。