联合反演

摘要： 如何利用观测到的地震图上尽可能多的信息来约束地下结构以及地震震源本身一直是地震学研究的前沿课题.近年来,随着计算机计算能力的提高,使用基于全波形反演的方法已被用于不同尺度结构成像中,并取得了良好的效果.但如何减小全波形反演对计算资源的巨大需求以及其反演的高度非线性仍是目前急需解决的问题.此外,对于区域以及全球尺度成像,全波形反演的波形的拟合仅限于相对较低的频率.目前,基于波形层析成像在区域尺度最高能拟合的频率大约为0.5 Hz,在全球尺度能拟合的频率更低,所以获得的波速模型的分辨率还有一定的改进空间.地震学体波和面波联合反演是另一种可以综合利用更多信息的成像方法.该种方法主要利用高频体波的走时信息以及面波的频散信息来约束地下结构.由于只需要求解高频近似下的波动方程,其效率较全波形反演有较大提高.相比于体波和面波数据单独反演,联合反演能利用体波和面波对地下结构约束的互补性来获得能同时拟合不同数据的波速结构模型.此外,体波和面波数据联合反演能获得更为准确的泊松比模型,因此可以更好地约束岩性、孔隙度、熔融程度等.鉴于目前海量的基于机器学习获得的不同震相的走时数据以及越来越多的密集流动地震观测,联合反演方法将在区域尺度地壳及上地幔结构的成像中发挥重要作用,进而为区域构造演化、地震灾害评估以及依赖波速结构模型的其他研究奠定基础.本文将回顾区域尺度常用的地震成像方法,介绍联合反演方法的基本原理以及川滇地区的应用,最后探讨一些未来的发展方向.

摘要： 接收函数与重力联合估计技术是一种估计地壳厚度和波速比参数的方法,自2011年由Lowry和Pérez-Gussinyé提出以来,经过国内外学者的发展和改进,该方法已日渐成熟,应用领域渐广.与单一接收函数H-κ叠加法相比,该方法由于引入了重力约束,可以在一定程度上降低估计的不确定性,获取台站下方更为准确的地壳厚度与波速比参数,为研究地壳结构与物质组成提供更为可靠的深部依据.本文首先简述了单一接收函数H-κ叠加法和重力似然函数估计的方法原理,并详细介绍了接收函数与重力联合估计技术的算法与流程;其次,剖析了该方法目前存在的一些技术问题,重点阐述了其相应的技术发展和各项优化;然后,归纳总结了该方法在国内外地区的应用进展;最后,展望了该方法未来的研究方向和发展趋势.

摘要： 作为青藏高原东北缘的祁连山,其东段与西段的地貌以及新生代地表变形样式差异显著.在新生代青藏高原向北东扩展这一共同构造背景之下,引起祁连山东段和西段这些差异的深部原因尚不清晰.本文基于一条宽频密集台站观测剖面,利用P波速度结构约束下接收函数和面波频散的联合反演方法给出祁连山东段的壳内及岩石圈地幔结构.结合前人在青藏高原东北缘的研究成果,对比祁连山东段和西段的岩石圈结构,探讨东、西段壳内变形方式的差异,从而揭示引起祁连山东西段这些差异背后的可能原因.东西段岩石圈结构的差异主要有:1)壳内低速层差异.在20~40 km的深度处,祁连西段(测线BB′)中地壳存在壳内低速物质(3.2~3.4 km·s^(-1)),而祁连东段(测线AA′)的中地壳并不存在低速层;2)岩石圈地幔差异.祁连西段的岩石圈地幔低速异常表现为分块的特征,而祁连东段的岩石圈地幔表现为比较统一的高速体.通过对比研究发现祁连西段壳内强度相对软弱,地壳的缩短变形程度更大,表现为高海拔的山脉以及更厚的地壳,而祁连东段壳内变形及缩短程度相对较弱,表现为地势相对较平坦的临夏盆地和陇中盆地以及较薄的地壳厚度.东、西段壳内强度以及变形方式的差异可能与深部岩石圈地幔的结构相关.

摘要： 射频大地电磁(RMT)是新兴的浅地表勘探法,其工作频段(10 kHz~1 MHz)处于音频大地电磁(AMT)与探地雷达(GPR)之间,随频率的增加,位移电流对RMT观测的影响逐渐增大,单一电阻率参数反演会影响数据解释的准确性.本文在同时考虑传导电流与位移电流影响的条件下,开展二维RMT数据电阻率和介电常数联合反演研究.为提高电阻率与介电常数之间的联合约束效果,在经典最小结构模型正则化约束的基础上,引入模糊C均值(FCM)聚类算法进行联合约束.模型试算表明,基于FCM聚类的二维RMT电阻率与介电常数联合反演可以在丰富解释参数的同时提高反演效果.

摘要： 鄂尔多斯盆地三叠系延长组致密砂岩储层具有地层水矿化度变化大、孔隙结构复杂、泥质含量较高等特点,含油饱和度评价困难。然而,致密储层钻井液侵入较浅,且多频介电测井蕴含地层水矿化度等信息,为饱和度评价提供了有利条件。基于阵列感应电阻率曲线以及介电测井4种频率的介电常数曲线和电导率曲线,提出一种联合电阻率测井和介电测井信息反演储层含油饱和度的方法。根据储层性质选择模型构建目标函数,电阻率测井响应采用西门杜方程,介电测井1 GHz频率下测量的介电常数和电导率使用复折射率模型,20~500 MHz频率测量的介电常数和电导率采用泥质砂岩响应模型。利用粒子群优化算法对目标函数进行寻优反演,并构建实例检验算法。将该方法应用于目标储层,获得饱和度、水相曲折度指数、冲洗带矿化度及泥质频散项体积含量。结果显示理论测井响应曲线和实际测井曲线的相对误差小于10%,计算的含油饱和度与岩心分析、试油数据吻合较好,表明该方法反演的参数结果有效、可靠,具有推广应用价值。

摘要： 单一类型地球物理数据的反演会存在多解性,联合反演方法能够有效克服多解性,其在过程的稳定性、模型的成像质量以及后续资料解释等方面都发挥出独特的优势。交叉梯度联合反演方法以其不依赖岩石物性关系的优点,成为研究的热点。本文以重震联合反演为例,详细研究了基于交叉梯度函数的联合反演方法以及反演中存在的相关问题。针对联合反演中不同数据灵敏度矩阵数值和反演系统差异较大的问题,利用更加灵活的反演策略,即分别交替进行数据拟合和交叉梯度结构约束。同时,还提出了正演网格和反演网格相分离的网格划分策略,平衡了计算效率和精度。采用合适的联合反演策略实现了重力数据和地震波初至数据之间的联合反演。通过模型试验表明,算法能够有效减小反演的多解性,提高反演的分辨率。

摘要： 地球物理反演问题通常是不适定的,反演结果往往具有多解性和不稳定性。为了改善单一地球物理反演中面临的这些问题,本文在单一重力/磁法反演的基础上,基于结构耦合的思想,引入交叉梯度算子,实现了密度和磁化率两种物性参数的相互耦合。通过设计不同理论模型组合,应用预条件共轭梯度算法对理论模型的重力和磁法数据开展交叉梯度联合反演计算,反演结果表明,对于单一棱柱体模型和组合棱柱体模型而言,重磁联合反演算法相比于单一反演结果都具有良好的稳定性,对于异常体的空间分布的刻画具有提升效果。

摘要： 地震动旋转分量可以提供传统平动分量所欠缺的空间梯度信息,理论上可以提升对震源参数的约束.本文对旋转分量用于矩张量反演进行了较为全面的数值试验.首先分析了均匀无限空间中位错点源产生的旋转波场与平动波场在辐射花样上的互补性,由此推测旋转分量用于矩张量反演能提供额外的约束.随后基于层状介质模型合成的平动和旋转记录,开展了地震矩张量反演测试研究.我们通过对比反演得到的矩张量、矩张量差异角ΔΦ和地震矩,定量评估在矩张量反演中引入旋转分量的作用,发现在合适的条件下,在原平动记录基础上加入旋转记录来联合反演会显著提升反演结果的可靠性,可减小矩张量差异角ΔΦ30%以上,极端条件下70%以上;如保持总记录数量不变(即在部分台站添加旋转分量替换其他台站的平动分量),某些矩张量元素反演效果会变好,但某些元素可能会变差.整体而言,在多数台站分布的情况下加入旋转都能起到强于增加同等数目平动台站的作用.同时,我们还发现在低频情况下,旋转分量对于深度非常敏感,引入旋转分量对于约束矩心深度有非常显著的效果.此外,还尝试了不同的试验条件,发现在低频、震中距较小、信噪比较低、速度结构模型较为准确的情况下,加入旋转分量对结果的提升更明显.整体而言,尽管旋转分量的使用需要满足一定的条件,将旋转用于矩张量反演是可行的,它能为矩张量反演提供与平动分量互补的信息,在稀疏台站和密集平动的台站的情况下加入旋转观测将显著增强对地震的监测能力.

摘要： 华南陆块属于比较稳定的陆块,地震活动不论是强度还是频度都比西部和华北地区弱很多,2019年在广西北流与广东化州交界处发生的M_(S)5.2中强地震为探讨华南板内中强地震成因提供了理想场所。本文计算了北流地震震源区布设的7个宽频带流动地震台站、120个短周期流动台站和4个固定台站的远震P波接收函数,结合前人在该区域的面波频散曲线,采用线性阻尼最小二乘反演算法联合反演了研究区地壳三维S波速度结构。得到如下结果:①研究区S波速度具有明显横向不均匀性,上下地壳速度在断裂附近存在明显差异性,新丰断裂南侧及东侧上地壳整体S波速度较高;②石窝断裂附近存在一条NW向连续性较好的低速带,震中区位置10 km深处存在明显局部高速异常,并随着深度减小异常逐渐增加,在浅处其与焦林断裂和新丰断裂附近的高速异常汇聚成一条NNE向连续异常带;③结合研究区地质调查结果和余震分布,震中区的局部高速异常可能是NW向石窝断裂南西盘发生右行走滑,并受到新丰断裂高速体的限制和互相作用发生了应力积累,从而引发此次中强地震;④地壳厚度泊松比反演结果显示研究区存在NE向明显地壳增厚(26~30 km),结合华南地区应力场分析,地壳增厚可能与中下地壳局部的NE向挤压作用相关。

摘要： FCZ地区阜三段(E_(1)f_(3))为三角洲前缘—滨浅湖沉积体系,主要发育滩坝型砂体,单层薄且横向变化快,导致滩坝砂体预测难度大,制约了滩坝型隐蔽油气藏勘探进度。通过细分砂层组开展沉积特征研究;利用近似古地貌恢复技术,明确特殊地质体有利发育区;综合运用频谱延拓与导向滤波处理技术,提升地震资料分辨率。基于地震属性深度学习与质控,优化砂体定性描述技术,融合敏感属性表征两套反射砂层组的平面展布特征。综合运用波形结构差异性与曲线重构声波阻抗开展联合反演,预测两套反射砂层组厚度。综合构造背景、油源输导、侧向封挡等条件,优选圈闭并提交有利目标,指导实际生产工作。

摘要： 大地电磁测深法因其频带范围宽,不受高阻屏蔽,探测深度可达上百公里,一直作为深部矿产资源探测、岩石圈动力环境等领域的重要研究手段之一.面波的频散特征是研究壳幔结构及横向不均匀性的一种强有力手段.基于NSGA-SBX联合反演大地电磁视电阻率和地震面波频散数据，该方法避免了这两种方法权重因子的选取和因二进制编码对反演精度的影响。本文设计速度与电阻模型含有相同和不相同间断的面模型算例，通过联合反演的结果可以发现联合方法仅在这两种物性参数有结构耦合关系的情况下，加强这种耦合关系，从而提高反演精度，而在结构不相似的地方，反演结果保留单方法的优势，不影响单独反演的结果。

摘要： 联合反演作为综合地球物理重要的定量解释工具,是当前地球物理研究的热点之一.在以往正则化联合反演研究中,不同地球物理方法间的模型约束采用相同的正则化因子,尽管此方案能够在一定程度上提高了解的稳定性和计算效率,但如果不同方法之间的模型约束泛函值存在数量级的差异,正则化因子就在更大程度上取决于模型约束泛函值较大的地球物理方法,不利于全面发挥各种方法的模型约束作用.鉴于此,本文以大地电磁测深法(MT)和重力正则化反演为例,提出了双正则化因子方案,更好地发挥了不同方法间模型稳定泛函的约束作用,通过模型试验验证了双正则化因子的有效性.模型实验表明:将双正则化因子引入到MT与重力联合反演后,相比于传统单正则化因子的联合反演,其更能发挥各种方法中模型约束项的约束效果,使得反演结果更趋合理.

摘要： 大地电阻率分布信息是影响磁共振地下水探测反演结果准确性的重要因素.在众多电磁法勘探技术中,瞬变电磁法具有高分辨率、高效率和大探测深度等优势,能准确探测地下几百米范围内的电阻率分布信息.因此磁共振与瞬变电磁联合解释方法具有重要意义.然而,利用单一测点拼接的磁共振与瞬变电磁联合解释方法进行模拟二维反演时存在解释结果不唯一,容易出现错误异常体等问题,尤其在复杂地质情况下,同一测线上相邻测点探测结果连续性差,解释结果偏离实际.基于此,本文提出磁共振与瞬变电磁横向约束联合反演方法(Laterally Constrained Inversion,简称LCI):引入外推积分法(quadrature with extrapolation,简称QWE),解决了传统正演过程中基于直接数值积分方法导致的求解效率低的难题,以保证联合反演方法顺利实施,进而以相邻测点地下结构应具备连续性为依据,引入横向约束反演思想.通过在联合反演目标函数中加入相邻测点间各模型参数约束矩阵,提高磁共振解释结果准确性,加强探测剖面地质结构和含水模型连续性.本文经过理论模型证实,LCI方法能有效提高传统一维反演解释结果的稳定性和唯一性.最后,针对安徽黄山野外实际探测数据进行横向约束联合反演,反演结果与地质资料基本一致,充分验证了磁共振与瞬变电磁LCI联合反演方法的实用性.本文的研究成果将为磁共振与瞬变电磁空间约束联合反演奠定基础.

摘要： 地球物理联合反演应用多种观测数据,通过反演地质体的物性和几何参数来求得同一个地质、地球物理模型。由于重力方法有较高的横向分辨率，重力场能直接反映地球介质密度变化和各种环境下动力学特征，地震方法则有较好的垂向分辨率，在划分地球圈层结构具有优势。因此，在所有不同种类地球物理方法中，地震与重力存在优势互补，最适合进行方法间的联合反演。利用重震联合反演方法对白云鄂博矿区一条北西向延伸约40km的重力与地震剖面进行深部构造解译。此次的重力-地震联合反演精细化了矿区深部构造分层，给出了地层起伏形态，重新解译了断裂展布位置，其中深大断裂与地质研究成果能够良好吻合，并更加精确。本区存在多期活动且能够沟通的深大断裂，可能为稀土碳酸岩岩浆交代白云岩（灰岩）提供了流体通道，对“先沉积铁、后交代稀土、多期多元”的成矿模式具有指导意义。

摘要： 青藏高原东北缘位于青藏块体、鄂尔多斯地块、阿拉善地块的交界区域.复杂的构造环境和强烈的地震活动性使青藏高原东北缘成为一个地震学研究的热点地区.近年来出现的背景噪声互相关方法，能够提取到面波经验格林函数，利用提取到的经验格林函数进行层析成像得到地壳上地幔高分辨率速度结构。本次研究数据主要来自于2012年1月-2013年12月布设在青藏高原东北缘地区的135个数字化测震台站，台站主要分布在甘肃、陕西、宁夏地区，利用三分量连续波形记录，通过互相关方法提取得到Rayleigh波和Love面波经验格林函数。然后采用时频分析方法和相位匹配滤波方法得到面波群速度和相速度频散曲线，利用Debayle和Sambridge的面波层析成像方法得到了8-50s共22个周期的Rayleigh和Love面波群速度和相速度分布图像。最后利用Herman的反演程序，通过Rayleigh和Love波联合反演得到了不同深度的剪切波速度结构。研究结果表明，短周期面波速度受沉积层约束，在银川盆地、海原盆地、鄂尔多斯地块和阿拉善盆地东部表现为低速异常，青藏高原东部波速呈现高速分布，海原断裂带为波速的过渡带；而周期大于30s的群速庋分布表明青藏高原地区呈现低速异常，以祁连山北缘断裂为界，北部表现为高速，南部为低速异常；青藏高原中下地壳的低速异常可能与中下地壳介质的软化有关。

摘要： 岩石圈是地球浅部的刚性圈层,由地壳和岩石圈地幔组成.研究表明,岩石圈构造变形、地震活动、油气生产、岩浆运移、矿物沉积等诸多热点问题都与岩石圈结构和性质的不均匀性(如温度、密度和化学组分结构的跳变等)相关.而其热结构和化学组分结构是这种不均匀性的重要内容,是深入认识岩石圈物理化学性质的基础.本研究利用重力异常、大地水准面高、地形、地表热流、地震面波频散曲线等地球物理观测资料，通过地球物理和岩石学的热动力学模拟，反演了华南地区岩石圈以及上地幔的热（温度）-化学组分结构，并以此为基础分析该地区岩石圈热-化学组分结构特征和岩石圈动力学过程作用。研究表明，华南缝合带（华夏块体）的岩石圈较薄，约为85~150千米；华南块体岩石圈厚度由东至西逐渐增厚，在四川盆地以下达到最大厚度250千米。反演得到的华南缝合带以及江南造山带的地幔岩石圈化学组成以饱满型橄榄岩为主（镁值Mg#约为88~90），该化学组分与中/新生代地幔橄榄岩包体的地球化学、岩石学分析结果一致，表明地球物理观测联合反演能够有效的评估岩石圈的化学组成，为缺乏地球化学观测的地区提供重要的参考信息。考虑到华南东部岩石圈地幔化学组成以饱满型橄榄岩为主、岩石圈较薄，笔者认为该地区古老的岩石圈地幔经历了广泛的交代变质或重富集的过程，该过程可能由古太平洋的平俯冲作用所导致。平俯冲作用导致了古老的大陆岩石圈下部的拆沉，使得软流圈地幔熔融上涌，同时导致剩余的岩石圈地幔发生交代变质和重新富集作用。反演得到的四川盆地以下的岩石圈地幔化学组分同样以饱满型橄榄岩为主，该结果与传统认为的克拉通地区亏损型的地幔岩石圈不同。导致这一差异的原因可能是由于四川盆地以下存在小尺度地幔对流或者拆沉过程。此外，青藏高原的东向挤压可能导致了四川盆地地幔岩石圈物质组成的改变。

摘要： 浅地表地球物理勘探，一般勘探深度从几十厘米到几十米的范围，直接为工程设计或施工服务，要求更高的分辨率及更好的定性定量解释。本文主要研究基于交叉梯度的直流电阻率与地震走时的联合反演，从结构上对两种介质物性予以约束。从测试结果得出，电阻率法与地震走时是一种较为有效互补的勘探方法。参数选取是联合反演中最重要的问题，包括单独反演的阻尼及平滑参数，联合反演中交叉梯度权重以及两种参与反演的模型之间的权重。针对模型之间权重问题，地震走时具有较高分辨率及容易产生次生干扰异常，电阻率反演具有较稳定平滑结果及较低分辨率，开始迭代时以地震模型结构约束为主，以此提高电阻率反演的分辨率；随着迭代的次数增加逐渐增大电阻率模型结构相似性约束权重，提高地震走时反演稳定性。提高速度与电阻率模型更高的相似性，减小反演的多解性。

摘要： 目前直流电法与大地电磁法已成为电法勘探领域中两种比较经典的方法,其数据采集、数据处理、数值模拟与反演研究都较为成熟.另外,由于观测数据量不足等原因,使得地球物理的反演结果不可避免的存在多解性,因此联合反演已成为近些年地球物理领域的热门话题.联合反演可分为基于相同物性的联合反演和基于不同物性的联合反演,其中前者是基于相同的岩石物性差异,具有相同的物理基础,其观测数据必然存在一定的相关性,因此此类联合反演具备天然的合理性.本文选择直流电法与大地电磁方法联合反演.直流电法对于地下浅部信息分辨率较高,但探测深度较浅;而大地电磁法勘探深度大,因此结合两种方法各自的优点,使两者之间相互约束,可获得相较于单一反演方法更好的反演结果。

摘要： 本文对1998年1月10日张北地震震源参数进行研究，透过将CAP(Cut And Paste)方法中的格点搜索替换为邻域反演算法(Neighbourhood Algorithm,NA)，实现了InSAR数据与近远震波形数据的联合反演。对反演结果的分析显示，InSAR数据可帮助反演得准确的水平位置和走向等震源参数，而地震波形数据约束了震源机制解中的倾角与滑移角，有效压制了InSAR单独反演中的折中(trade off)效应。由此可见，InSAR与地震学数据的联合反演可提供更可靠的震源参数结果。

摘要： 激发极化衰减信息与水量的正相关关系为水量的定量估算提供了理论基础，为不良地质体中充填水量的足量估算提供可行途径，前向三维激发极化适用于TBM环境，为TBM施工提供了可用的超前地质预报方法，联合反演被在降低多解性方面有较好效果。

摘要： 本发明公开了联合叠前弹性反演参数及深度网络目标反演技术，涉及目标反演技术领域，该目标反演技术具体步骤如下：获得叠前弹性阻抗反演参数；获得划分角度的叠前参数；将井旁道数据作为输入，目标数据作为输出约束；建立训练网络；通过训练获得映射网络，直接对全区数据进行映射计算；输入全区数据参数，并通过映射网络直接获得目标数据体；本发明能够更好的针对地震目标参数进行预测，提高了预测的准确性，高效性，深度网络深度挖掘数据内部信息，提高信息的表征度，解决了通过拟合求解目标数据过程中存在的累积误差和近似误差，直接通过叠前反演获得的参数和不同角度数据增强了算法的可信度和真实性。

摘要： 本发明公开了一种基于激光雷达的联合观测与消光系数标定反演方法，在待测区域水平和垂直点分别布设激光雷达，利用两台激光雷达联合观测，一台用于垂直定点观测，获取大气边界层结构和污染传输时空分析，一台用于水平扫描，获取本区域内污染分布特征，并就水平扫描激光雷达基于Collis方法获取消光系数，将该消光系数作为初始值作为垂直激光雷达消光系数的标定，实现探测方式的优势互补，最终实现对区域污染分布进行时空立体监测，得到立体空间气溶胶消光系数分布特征，并能够有效地揭示城市地区较大区域内的颗粒物分布和传输情况，实现以三维数据分析颗粒物的传输、分布和浓度变化情况，具有覆盖范围大、探测效率高的优点。

摘要： 本发明提供一种基于分频联合反演的叠前高分辨率流体因子反演方法，该方法包括：基于孔弹性理论对传统形式的精确Zoeppritz方程进行改写得到包含流体因子项的新形式精确Zoeppritz方程；利用工区实际地震数据的分频处理结果、测井数据统计得到的多个储层参数的先验分布函数及推导的新形式Zoeppritz方程构建贝叶斯框架下的多波数据分频联合反演目标函数；利用泰勒级数展开对目标函数进行求解得到所述储层参数扰动量的求解表达式；对求解表达式进行迭代求解，利用求解结果对初始储层参数值进行迭代更新，直至达到迭代终止条件。本发明能够提高储层流体因子的预测精度。

摘要： 本发明公开了一种提高高陡构造速度反演效率的全波形联合反演方法，属于石油地球物理勘探领域，本发明通过输入初始速度场、炮记录及震源子波，建立观测系统；使用传统逆时偏移获得成像结果并作为反射系数模型，应用线性正演模拟和线性波场反传；计算两种棱柱波的梯度方向并相加求得棱柱波形反演的梯度方向，用线性搜索方法或者抛物拟合方法求取更新步长，用求取的梯度方向和步长更新速度；用常规全波形反演方法再次更新速度，判断是否满足误差条件，如果不满足误差条件，利用更新的速度作为输入速度再次更新速度，如果满足误差条件，输出反演结果。本发明降低了对初始模型中高陡构造信息是否缺失的依赖性，提高了对高陡构造的反演精度和效率。

摘要： 本发明首先利用测井数据正演分析确定频段划分原则，然后基于划分原则对实际数据进行分频处理，将其变换为不同频段对应的角道集数据，然后将这些角道集数据作为输入进行分频联合反演，并合理提升高频地震信息的权重，将最终的反演结果与基于原始角道集数据的反演结果对比可以发现，本发明提出的方法能够进一步有效的提升反演结果的分辨率，提升反演结果对薄层的刻画能力。

摘要： 本发明提供一种基于分频联合反演的叠前高分辨率流体因子反演方法，该方法包括：基于孔弹性理论对传统形式的精确Zoeppritz方程进行改写得到包含流体因子项的新形式精确Zoeppritz方程；利用工区实际地震数据的分频处理结果、测井数据统计得到的多个储层参数的先验分布函数及推导的新形式Zoeppritz方程构建贝叶斯框架下的多波数据分频联合反演目标函数；利用泰勒级数展开对目标函数进行求解得到所述储层参数扰动量的求解表达式；对求解表达式进行迭代求解，利用求解结果对初始储层参数值进行迭代更新，直至达到迭代终止条件。本发明能够提高储层流体因子的预测精度。

摘要： 本发明公开了一种提高高陡构造速度反演效率的全波形联合反演方法，属于石油地球物理勘探领域，本发明通过输入初始速度场、炮记录及震源子波，建立观测系统；使用传统逆时偏移获得成像结果并作为反射系数模型，应用线性正演模拟和线性波场反传；计算两种棱柱波的梯度方向并相加求得棱柱波形反演的梯度方向，用线性搜索方法或者抛物拟合方法求取更新步长，用求取的梯度方向和步长更新速度；用常规全波形反演方法再次更新速度，判断是否满足误差条件，如果不满足误差条件，利用更新的速度作为输入速度再次更新速度，如果满足误差条件，输出反演结果。本发明降低了对初始模型中高陡构造信息是否缺失的依赖性，提高了对高陡构造的反演精度和效率。

摘要： 本发明提供了一种基于深度学习的隧道电阻率极化率联合反演梯度优化方法，获取隧道三维地质数据；根据隧道三维地质数据以及预设的JointGradOptNet网络，得到隧道电阻率和极化率的联合反演梯度的优化预测结果；其中，JointGradOptNet网络的训练中，针对基于隧道前方常见的含水体构建隧道三维地质模型数据库，通过一次迭代的反演结果与初始均一地电模型的差得到网络的电阻率输入梯度与极化率输入梯度，通过地电模型与初始均一地电模型的差得到网络的电阻率目标梯度与极化率目标梯度；本发明实现了对隧道掌子面前方异常体较为精确的成像。

摘要： 本发明提供一种联合GNSS三维形变的区域陆地水储量变化反演方法及系统，包括根据质量负荷的格林函数理论，建立GNSS北向、东向和垂直形变与区域陆地水储量变化之间的关系；根据质量负荷格林函数理论建立区域陆地水储量变化与GNSS北向、东向和垂直形变之间的观测方程，并得到相应法方程；根据最小二乘联合平差反演区域陆地水储量变化，包括根据法方程并结合先验信息方程形成GNSS垂直和水平形变的联合反演模型，然后给定观测值噪声方差以及正则化参数初始值，利用方差分量估计方法迭代计算各类数据的最优权比，得到最优联合反演的区域陆地水储量变化。本发明提高了区域陆地水储量变化反演结果的精度和可靠性。

摘要： 一种基于余弦点积梯度约束的三维重磁联合反演方法，包括基于几何格架法的重力、重力梯度和磁法三维正演快速算法；依据余弦定理推导出三维余弦点积梯度函数，基于截断高斯牛顿法的三维联合反演算法，该方法既能避免存储雅克比矩阵，又能减少共轭梯度法求解线性方程时正演计算的次数。在传统重力和磁法数据中又加入了全张量重力梯度数据，能够进一步改善重磁联合反演的分辨能力；将上述正演、结构耦合方式以及联合反演算法相结合应用于重力、重力梯度和磁法数据中，充分发挥三方面的优势，有效地改善三维重磁联合反演的计算效率、内存消耗和分辨能力。