基于交叉梯度交替结构约束的二维地震走时与全通道直流电阻率联合反演

在利用不同的地球物理勘探方法对地下复杂介质成像时,因观测系统的非完备性及数据本身对某些岩石物性的不敏感性,单独成像的结果存在较大的不确定性和不一致性.对于地震体波走时成像与直流电阻率成像,均面临着成像阴影区问题.对于地震走时成像,地震射线对低速区域覆盖较差形成阴影区,造成低速区域分辨率降低.对于电阻率成像,电场线在高阻区域分布较少,造成高阻区域分辨率较低.为了提高地下介质成像的精度,Gallado和Meju(2003)提出了基于交叉梯度结构约束的联合地球物理成像方法.在要求不同的物性模型拟合各自对应的数据同时,模型之间的结构要求一致,即交叉梯度趋于零.为了更有效地实现基于交叉梯度的结构约束,我们提出了一种新的交替结构约束的联合反演流程,即交替反演不同的数据而且在反演一种数据时要求对应的模型与另一个模型结构一致.新的算法能够更容易地把单独的反演系统耦合在一起,而且也更容易建立结构约束和数据拟合之间的平衡.基于新的联合反演流程,我们测试了基于交叉梯度结构约束的二维跨孔地震走时和直流电阻率联合成像.合成数据测试表明,我们提出的交替结构约束流程能够很好地实现基于交叉梯度结构约束的联合成像.与单独成像结果相比,地震走时和全通道电阻率联合成像更可靠地确定了速度和电阻率异常.     

基于交叉梯度耦合的大地电磁与地震走时资料三维联合反演

为了有效解决目前大地电磁和地震走时资料单方法反演结果一致性不好的问题,同时克服基于岩石不同物性参数间关系耦合约束联合反演的局限性,本文研究了基于交叉梯度耦合约束的大地电磁与地震走时资料的三维联合反演算法.以较为成熟的天然地震走时资料三维正反演和大地电磁三维正反演算法为基础,实现了具有共同的反演网格,以交叉梯度结构耦合约束,并能同时获得电阻率和速度模型的三维联合反演算法.分别利用单棱柱体模型和双棱柱体模型合成数据进行了联合反演试算.结果表明:无论是单棱柱体模型还是双棱柱体模型,联合反演结果比单独反演对异常体的空间形态都有更好的恢复,其中单棱柱体模型反演的异常体电阻率更接近于真实电阻率,双棱柱体模型的联合反演结果不仅消除了围岩的部分电阻率假异常,而且增强了对异常体深部速度结构特征的恢复程度.联合反演还能同时改善电阻率和速度反向变化异常体的单独反演结果,进一步证明交叉梯度耦合不依赖于岩石物性关系,而强调地下结构的相似性,具有更普遍的适用性.     

基于交叉梯度约束的电阻率法和背景噪声法三维联合反演研究

电阻率法和背景噪声法是通过获得地下介质的电阻率和速度参数的分布来探究地球内部物质分布的非均匀性特征,联合反演可以有效地发挥两种方法的优势,减小反演多解性.本文采用有限内存拟牛顿（L-BFGS）算法实现了电阻率法和背景噪声法的单方法三维反演,然后基于电阻率和速度模型结构耦合约束,采用交替迭代方式实现了电阻率法和背景噪声法的三维联合反演.通过设计规则体组合模型和不规则体组合模型进行理论模型合成数据三维反演试算,结果表明：联合反演可以获得同时满足数据拟合差和结构更为相似的速度-电阻率模型;联合反演可以压制背景噪声单方法反演出现的假异常,获得更精细的速度结构;联合反演获得的电阻率模型对倾斜异常体、高阻覆层下方异常体、圈闭的高/低阻体等边界信息有明显的提升,有效克服电阻率法单方法反演的局限,提高深部电阻率的分辨率.     

地震走时和航空电磁联合反演三维速度和电阻率

地震走时数据可以反演近地表的速度,但不能反演出隐蔽层和低速层。航空电磁数据可以反演近地表的高电阻率和低电阻率,但是对垂直方向的分辨率低。联合三维地震走时数据和航空电磁数据反演了近地表的速度和电阻率结构,并采用棋盘模型测试和实际数据应用展示了联合反演出的速度模型优于单独的走时反演出的速度模型。结果表明:该算法可以应用到处理数据量大的勘探地震成像中,能提供优化的速度结构。     

采空区地震与电法探测数据的联合反演研究及应用

综合地球物理技术在采空区的探测中发挥了重要作用.目前通常采用单方法反演、仅对不同方法反演结果进行对比解释的综合勘探方式,单方法反演的多解性严重降低了其探测精度.如何提高采空区的探测精度,对采空区进行有效探测一直被认为是地球物理技术面临的首要难题.为了提高地震与电法技术的探测精度,基于交叉梯度联合反演理论,设计了地震初至折射走时数据和高密度电法数据的联合反演算法流程,对采空区理论模型和野外实际数据进行了联合反演处理.结果发现通过两者的联合反演,不仅可以提高采空区电阻率反演模型的成像效果,而且能够获得地震单方法反演难以成像的采空区低速异常体,从而提高了地震与电法技术对采空区的探测精度.表明地震与电法探测数据联合反演是一种提高采空区探测精度的有效方法.     

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

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

一种高效的基于交叉梯度结构约束的三维地震走时与直流电阻率联合反演策略

基于交叉梯度结构约束的不同数据的联合反演可以提高地球物理成像的可靠度,但是由于不同观测数据对地下介质的灵敏度不同以及不同反演系统灵敏度矩阵元素的数值可能存在较大的差异,导致构建和求解联合反演系统存在很大的挑战.针对以上问题,本文提出一种新的基于单独反演模型更新量与交叉梯度结构约束相结合的联合反演策略.该策略利用单独反演系统分别确定出各个模型的更新量,然后利用它们约束交叉梯度系统的反演,得出新的模型更新量.通过这样的方式,有效实现了数据拟合与结构约束的平衡,实现了基于交叉梯度联合反演的目标.新的联合反演策略不需要对原来单独反演程序修改即可实现联合反演,减小了联合反演实现的难度,极大地提高了联合反演的易实现性,而且避免了联合反演矩阵存储及结构过于复杂难以求解的问题.基于新的策略,本文首次实现了基于交叉梯度结构约束的三维地震走时与直流电阻率联合反演.合成模型测试表明,与单一成像相比,联合成像减少了地震走时反演中出现的干扰异常并提高了电阻率反演的分辨率.     

二维复杂层状介质中地震多波走时联合反演成像

采用新近提出的多次波射线追踪正演算法,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,分析讨论了利用多震相走时资料进行联合反演成像的方法及技术.考虑到不同震相走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数; 由于同时反演速度模型和反射界面起伏中不同模型参数变化对走时影响程度的不同, Jacobi偏导矩阵元素中引入了不同参数的归一化因子; 另外,为了克服射线密度过大(或过小)区域速度模型的过度(或欠)更新问题,反演算法中引入了等权射线密度的概念.几种数值模拟实例表明(含噪声敏感性试验): 多波走时的联合或同时反演成像技术是一种提高走时成像空间分辨率,进而降低重建模型失真度行之有效的方法.     

基于物性参数梯度累加约束的电阻率法和重力勘探数据二维联合反演

对不同地球物理方法的数据进行联合反演,通过模型参数的相互约束可以减小反演的多解性.本文对电阻率法和重力勘探数据开展了基于电阻率和剩余密度梯度累加约束的二维联合反演研究.电阻率法采用数据空间Occam反演,重力勘探采用基于对数障碍法的正则化反演.通过在电阻率法和重力勘探反演的目标函数中引入电阻率和剩余密度梯度累加约束项,利用交替迭代方式实现了电阻率法和重力勘探数据的二维联合反演.对理论模型合成数据进行了单方法反演和联合反演,对比了反演效果.结果表明：联合反演结果优于单方法反演结果,联合反演结果对异常体的物性值恢复和形态刻画效果更好.     

三角网格剖分下速度与反射界面的同时反演

为了解决复杂速度模型中的走时正、反演问题,例如:含不规则起伏地表、不规则地下波阻抗界面、以及不规则速度异常体的复杂地学模型,本文采用三角网格单元模型参数化下的分区多步改进型最短路径算法,实现了多震相地震射线的追踪计算,结合共轭梯度法求解带约束的阻尼最小二乘反演问题,实现了多震相走时联合同时反演成像的方法技术.当界面起伏较大时出现散射,从而造成散射点所在区域射线密度过密,导致该区域内速度和界面的过度更新.为了克服上述问题,我们在同时反演中引入了射线密度的概念,从而有效地解决了上述过度更新问题.数值模拟实验表明:采用三角网格单元进行模型参数化,可保证在复杂模型中的正演计算具有较高的计算精度;同时反演中可以准确地刻画不规则异常体和不规则反射界面.因此,本文提出的走时成像方法技术具有较广的实用价值.     

基于初至波走时层析成像的Tikhonov正则化与梯度优化算法

初至波走时层析成像是利用地震初至波走时和其传播的射线路径来反演地下介质速度的技术.该问题本质上是一个不适定问题,需要使用正则化方法并辅之以适当的最优化技巧.本文从数值优化的角度介绍了初至波走时层析成像的反演原理,建立了Tikhonov正则化层析成像反演模型并提出求解极小化问题的加权修正步长的梯度下降算法.该方法可以从速度模型的可行域中迭代找到一个最优解.数值试验表明,该方法是可行和有应用前景的.     

多震相走时联合三参数同时反演成像

采用新近研制的分区多步不规则最短路径多震相地震射线追踪正演技术,结合流行的子空间反演算法,提出了一种联合多震相走时资料进行地震三参数 （速度、反射界面和震源位置） 同时反演的方法技术.数值模拟反演实例、以及与双参数 （速度和反射界面或速度和震源位置） 同时反演的对比分析表明：三参数同时反演成像结果大体接近双参数同时反演成像的结果.另外,噪声敏感性试验表明：所提算法对到时数据中可容许的随机误差并不敏感,结果说明多震相走时的联合三参数同时反演成像方法技术不失为一种提高走时成像空间分辨率、进而降低重建模型参数失真度、行之有效的方法技术.     

波动方程初至波多信息联合反演方法

地震初至波中包含着丰富的近地表速度结构信息,如何分阶段、分尺度地利用这些信息进行近地表速度建模是地震勘探中的一个关键问题.在速度反演的不同阶段,综合利用初至波中的不同信息(如走时、包络和波形等)进行联合反演,可以有效地降低反演对初始模型的依赖程度,提高近地表速度模型的反演精度.为此,本文提出了一种统一基于波动方程正演引擎的初至波多信息联合反演方法,该方法同时匹配观测和模拟的初至波走时、包络和波形信息.在不同的反演阶段选择不同的权重因子调节不同信息的权重,这样不仅降低了反演对初始速度模型的依赖,而且自然地实现了多尺度反演.在每轮反演迭代中,一次正演模拟的波场同时应用于初至波走时、包络和波形匹配,无需额外的射线追踪.同时,联合反演方法在一定程度上缓解了串联反演中目标函数漂移问题,提高了近地表速度建模的精度.     

电阻率和速度随机分布的MT与地震联合反演

在已有研究成果的基础上,为了适应物性参数剧烈变化的复杂模型并满足联合反演的要求,开发了速度和电阻率随机分布共网格单元模型的建模技术.基于这种统一的物性随机分布的网格介质模型,利用有限元方法和改进的射线追踪法分别正演计算大地电磁场和地震走时,结合改进的模拟退火算法,研究实现了电阻率和速度随机分布条件下的大地电磁与地震资料的同步联合反演.对物性界面不完全一致和物性变化剧烈的带地形复杂模型的试验,表明了该方法在精细反演复杂电阻率和速度结构方面的效果,克服了以往研究局限于简单模型的不足.对地震资料品质差的地区开展的实际资料联合反演,表明了方法的适用性,先验信息约束下的联合反演提高了反演精度.     

求解二次方程的二维VTI介质qSV波和qSH波走时快速扫描算法

地震波走时计算在观测系统设计、偏移成像、速度模型走时反演和地震定位等方面起到重要作用.各向异性广泛存在于地球介质中,影响地震波传播的振幅和走时,忽略各向异性的影响将对成像、反演以及地震定位等造成一定的误差.因此对于高分辨率成像和反演,走时计算中考虑各向异性十分重要.快速扫描法不需要存储和追踪波前面信息,在各向异性初至波...     

多震相菲涅耳体射线走时同时反演成像

高频假设下的地震射线理论以及相应的地震成像理论表明,在射线稀疏条件下,不可能得到较高分辨率的构造成像;而有限频射线理论更符合实际地震的传播规律,即地震波的走时不仅与中心射线(传统的几何射线)上的速度分布有关,而且与中心射线附近一定范围(称其为第一菲涅耳体)内的速度异常分布有关.鉴于此,本文提出了计算多震相地震波菲涅耳体有限频射线的方法,并定义了走时敏感核函数,同时给出了利用多震相菲涅耳体有限频射线进行速度模型和反射界面同时反演成像的公式.利用多震相走时资料,使用传统射线层析成像方法与有限频射线层析成像方法进行了速度和界面的同时反演成像.结果表明,当射线密度较小时,无论是对速度模型的重建还是对反射界面几何形状的更新,有限频射线层析成像方法均优于传统射线层析成像方法, 而变频有限频射线层析成像则是实际地震层析成像的首选反演算法.      

直流电阻率法与大地电磁法的二维联合反演

本文对电性联合反演进行了深入研究,以减少地球物理反演的多解性.将直流电阻率(DC)与大地电磁(MT)数据加入到同一反演数据集中.引入Tikhonov正则化思想建立反演目标函数,使反演过程更加高效稳定.在解决正则化反演问题过程中,分别采用了二阶最大平滑稳定因子和改进的L-curve法,提高了反演结果的稳定性和正则化因子的求取精度;最后运用非线性共轭梯度法(NLCG)对反演目标函数实现最优化求解.经研究表明:联合反演方法与单一反演方法相比,能够更加有效的约束反演模型范围;反演算法快速稳定,提高了反演精度,减少了对地下地质结构认识的模糊性.     

基于最小支撑泛函模型约束的电阻率法三维聚焦反演

目前电阻率法三维反演方法大多是基于最小构造或最平缓模型约束的反演,这些反演算法能稳定收敛,但有时反演结果分辨率较差,不利于地质解释.本文在分析Zhdanov(2002)提出的基于最小支撑泛函聚焦反演方法的基础上,将最小支撑泛函引入到电阻率法三维反演的目标函数中,然后采用高斯牛顿法求解反演目标函数最优化问题,同时结合预条件共轭梯度法得到电阻率法三维聚焦反演结果.通过对几个典型模型的试算,并与传统光滑模型约束反演结果进行比较,表明本文反演方法结果与实际模型吻合的更好,分辨率更高,模型更聚焦,而且算法收敛速度较快.     

三维多层介质重力-地震同步联合反演

联合反演是地球物理勘探的重要解释手段,能够提高模型参数的反演精度.本文在归纳和分析重力与地震资料联合反演的研究和应用现状的基础上,利用三维多层介质模型的地震走时和重力正演公式,推导了地震走时和重力异常对界面深度的雅可比矩阵,实现了三维重力-地震同步联合反演界面成像.最后进行了数值理论模型模拟和实例计算,结果表明地震走时和重力同步联合反演很好的重建了三维多层介质界面.     

大地电磁、重力、磁法和地震初至波走时的交叉梯度二维联合反演研究

为了准确的探测和描绘地下复杂的地质结构,同时克服地球物理单一方法反演的多解性和单一参数反演模型的不一致性等问题,近年来基于交叉梯度联合反演的综合地球物理解释已经得到了广泛的关注和应用.本文首先研究了两种地球物理方法的交叉梯度联合反演算法,在此基础上,推导并实现了多种地球物理方法(大地电磁,重力,磁法,地震初至波走时)的多交叉梯度约束的二维联合反演算法;其次,我们设计了结构不一致模型和复杂模型,针对多物性联合反演算法的准确性和有效性进行了模拟试算,并对复杂模型的单独反演结果和联合反演结果进行了交叉梯度值和物性交会图的对比;最后,本文将成熟的卫星资料多光谱综合分析技术应用到联合反演中,将多物性参数反演模型结果图通过RGB(红-绿-蓝)模式进行合成,得到融合的RGB合成图.结果表明:通过对结构不一致模型和复杂模型的联合反演结果和单独反演结果的对比分析,可以得出联合反演得到的结果更接近真实模型,并从得到的交叉梯度值进一步证明了联合反演模型相似度高,也从物性交会图中得到联合反演的物性相关性更好的结论,反向证明了算法的正确性.最终从得到的RGB合成图像,我们可以更直观的分析反演结果,更有利于准确划分地下模型结构.