共查询到20条相似文献，搜索用时 31 毫秒
1.
无线电波透视法是常用的工作面地质构造探测方法之一,目前普遍使用的SIRT方法层析分辨率不高。本文采用约束正则化方法,推导Tikhonov正则化和全变差正则化的最小化问题表达式,讨论影响层析结果的主要因素,对典型理论模型进行了层析成像实验。结果表明:正则化方法具有比SIRT方法更好的分辨率;射线条数越多、噪声水平越低,层析分辨率越高;Tikhonov正则化在正则参数增大时层析结果更光滑,减小时则更贴近异常,全变差正则化与其相反。最后对实际坑透数据进行层析,识别出的异常构造基本吻合已知疑似构造位置,从而说明正则化方法在无线电波透视应用中的可行性。 相似文献
2.
针对Mogi模型垂直位移与水平位移联合反演中的病态问题,改进火山形变总体最小二乘(Total Least Squares,TLS)联合反演的虚拟观测法,并使用方差分量估计(Variance Components Estimation,VCE)方法确定病态问题的正则化参数.将附有先验信息的参数作为观测方程,与垂直位移和水平位移的观测方程联合解算,推导了三类观测方程联合反演的求解公式及基于总体最小二乘方差分量估计确定正则化参数的表达式,给出了算法的迭代流程.通过算例实验,研究了总体最小二乘联合反演的虚拟观测法在火山Mogi模型形变反演中的应用;算例结果表明,三类数据的联合平差及方差分量估计方法可以确定权比因子并得到修正后的压力源参数,具有一定的实际参考价值. 相似文献
3.
4.
本文针对蒸汽辅助重力泄油(SAGD)生产中开发监测问题,发展了综合应用地震及重力数据反演储层密度的联合反演算法.通过测井数据建立纵波阻抗与密度的直接关系,并推导出这种关系下重力与纵波阻抗数据联合反演的计算方法,从而计算出蒸汽腔体密度分布规律.文中应用密度反演后验约束正则化方法,采用Tikhonov正则化模型,通过波阻抗数据作为约束进行联合反演,在算法上提高了稳定性,同时得到较高的反演精度.文中对SAGD生产中的理论模型进行了方法试算,并分析了算法的误差,最终应用于SAGD生产的实际数据中,通过最终反演结果分析,该方法取得了很好的应用效果. 相似文献
5.
针对探地雷达（GPR）双参数全波形反演中电导率反演精度差、双参数存在串扰现象、反演计算量大、易陷入局部极值等问题.作者将具有多参数调节功能的L-BFGS算法引入到GPR时间域全波形反演中,它避免了对Hessian矩阵的直接存储与精确求解,减小了存储量和计算量.结合参数调节因子的选取,有效减小了同步反演时介电常数与电导率的串扰影响,在不降低介电常数反演精度的前提下,提高电导率参数的反演精度.通过在反演目标函数中加载改进全变差正则化方法,提高了反演的稳定性,使目标体边缘轮廓更加清晰.首先以简单模型为例,对比了单尺度反演与多尺度串行反演策略的优劣,说明多尺度串行反演有利于逐步搜索全局最优解;而开展参数调节因子的选取实验,说明合适的参数调节因子可以有效改善介质电导率的反演精度;测试了不同正则化的反演效果,表明改进全变差正则化能提高反演稳定性,显著降低模型重构误差.最后,分别对含噪合成数据和实测数据进行了反演测试,说明本文提出的多尺度、双参数反演具有较强的鲁棒性,能提供更丰富的信息约束,重构图像界面清晰、反演效果好. 相似文献
6.
地球物理反演是获取地球信息的重要手段,其求解具有严重的不适定性.为获得稳定的反问题结果,通常需要在目标泛函中加入正则化约束项.正确地估计正则化参数一直是地球物理反问题中的难点.目前存在的选取方法需要根据大量的试验来确定正则化参数,工作量十分巨大,并且存在很大的经验性,很难得到最优的正则化参数.针对这个问题,本文提出了一种基于广义Stein无偏风险估计的正则化参数求取方法.该方法的具体思路是通过求解模型参数均方误差的广义Stein无偏风险估计函数,在反问题求解过程中自动求取正则化参数.本文模型测试结果表明,相比于目前常用的方法,通过该方法得到的正则化参数是最优的. 相似文献
7.
联合反演是综合地球物理研究的重要定量解释手段.本文在总结和分析重力与地震资料联合反演的研究现状基础上,利用改进的全局寻优的快速模拟退火算法,实现了重力和地震资料的约束同步联合反演.针对性地设计了密度和速度界面不完全一致的模型,理论模型的试验说明了方法的效果和适用性.结合最近完成的广东徐闻地区实际资料的处理和解释,表明该方法可准确确定复杂构造物性界面的密度和速度结构,在该地区的油气勘探中发挥了作用.在先验信息约束下,该联合反演方法要明显优于单独的重力反演. 相似文献
8.
复杂地表探区,尤其是盆山过渡区的油气勘探是我国也是世界上油气勘探的重点区域,但是此类区域油气地震勘探中满足精确地震成像的速度建模一直是个没有很好解决的问题.本文提出了一种综合性的数据域初至波走时与成像域反射波走时联合层析复杂地表浅中深层速度建模方法,并针对联合层析速度反演解的非唯一性问题,深入地分析了层析反演中正则化的本质意义,指出了建立构造特征正则化方法的具体技术路线,提出了联合层析的实现流程及策略.理论和实际数据试验表明,本文提出的数据域初至走时与成像域反射走时联合层析浅中深层速度建模技术避免了常规建模方法中浅层速度模型与中深层速度模型的融合问题,较好地解决了传统成像域反射层析对近地表模型的不可控更新问题,整体提升了深度域浅中深层速度模型的建模精度,进而提高了复杂地表、复杂构造区的地震成像质量. 相似文献
9.
基于交叉梯度结构约束的不同数据的联合反演可以提高地球物理成像的可靠度,但是由于不同观测数据对地下介质的灵敏度不同以及不同反演系统灵敏度矩阵元素的数值可能存在较大的差异,导致构建和求解联合反演系统存在很大的挑战.针对以上问题,本文提出一种新的基于单独反演模型更新量与交叉梯度结构约束相结合的联合反演策略.该策略利用单独反演系统分别确定出各个模型的更新量,然后利用它们约束交叉梯度系统的反演,得出新的模型更新量.通过这样的方式,有效实现了数据拟合与结构约束的平衡,实现了基于交叉梯度联合反演的目标.新的联合反演策略不需要对原来单独反演程序修改即可实现联合反演,减小了联合反演实现的难度,极大地提高了联合反演的易实现性,而且避免了联合反演矩阵存储及结构过于复杂难以求解的问题.基于新的策略,本文首次实现了基于交叉梯度结构约束的三维地震走时与直流电阻率联合反演.合成模型测试表明,与单一成像相比,联合成像减少了地震走时反演中出现的干扰异常并提高了电阻率反演的分辨率. 相似文献
10.
为解决地球物理反演中多解性的问题,综合多种地球物理信息的联合反演受到了广泛的关注.本文依据不同地球物理响应可能由相同异常体引起,而不同地球物理分布参数之间存在相关性等特点,提出了一种基于局部Pearson相关系数约束的联合反演方法.该方法假设每个局部区域模型参数的分布具有线性相关特性,在拟合不同类型观测数据时,对局部模型参数施加相关性约束,进行联合反演以减少多解性.本文采用交替迭代联合反演流程,改善了同一目标函数下联合反演收敛性和速度问题.基于新的联合反演方法和流程,我们测试了三维大地电磁和重力仿真数据的联合反演.结果表明,本文提出的基于局部相关性约束的联合反演方法,能充分利用大地电磁和重力观测数据信息,有效改善单一地球物理反演收敛性和多解性的问题,反演效果得到明显提升. 相似文献
11.
12.
反射地震走时层析本质上是一个病态问题,而正则化是改善问题病态程度的有效手段.反射地震走时层析最终可归结为线性方程组的求解,本文讨论了在线性方程组求解过程中正则化的作用和方式.正则化的作用有:(1)用超定分量约束欠定分量和零空间分量;(2)用先验信息约束欠定分量和零空间分量;(3)对射线的不均匀覆盖进行阻尼;(4)对数据的不准确性进行阻尼.正则化的加入方式有:(1)加法型(将正则化矩阵补在层析矩阵后面,包括导数型正则化和零阶正则化,一阶导数型正则化对应最平坦解,二阶导数型正则化对应最光滑解,零阶正则化对应紧约束解);(2)乘法型(将正则化矩阵与层析矩阵相乘,主要包括阻尼型正则化).并利用简单的模型对正则化的效果进行了试验,发现经各种正则化约束后,与未加任何正则化约束得到的速度模型比较,尽管恢复的异常体的幅度不如后者大,但得到的速度剖面要平滑得多,更利于后续的射线追踪正演和层析反演. 相似文献
13.
The technique of seismic amplitude-versus-angle inversion has been widely used to estimate lithology and fluid properties in seismic exploration. The amplitude-versus-angle inversion problem is intrinsically ill-posed and generally stabilized by the use of L2-norm regularization methods but with drawback of smoothing important boundaries between adjacent layers. In this study, we propose a sparse Bayesian linearized solution for amplitude-versus-angle inversion problem to preserve the sharp geological interfaces. In this regard, a priori constraint term with two regularization functions is presented: the sparse constraint regularization and the low-frequency model information. In addition, to obtain high-resolution reflectivity estimation, the model parameters decorrelation technique combined with dipole decomposition method is employed. We validate the applicability of the presented method by both synthetic and real seismic data from the Gulf of Mexico. The accuracy improvement of the presented method is also confirmed by comparing the results with the commonly used Bayesian linearized amplitude-versus-angle inversion. 相似文献
14.
电阻率法和背景噪声法是通过获得地下介质的电阻率和速度参数的分布来探究地球内部物质分布的非均匀性特征,联合反演可以有效地发挥两种方法的优势,减小反演多解性.本文采用有限内存拟牛顿（L-BFGS）算法实现了电阻率法和背景噪声法的单方法三维反演,然后基于电阻率和速度模型结构耦合约束,采用交替迭代方式实现了电阻率法和背景噪声法的三维联合反演.通过设计规则体组合模型和不规则体组合模型进行理论模型合成数据三维反演试算,结果表明：联合反演可以获得同时满足数据拟合差和结构更为相似的速度-电阻率模型;联合反演可以压制背景噪声单方法反演出现的假异常,获得更精细的速度结构;联合反演获得的电阻率模型对倾斜异常体、高阻覆层下方异常体、圈闭的高/低阻体等边界信息有明显的提升,有效克服电阻率法单方法反演的局限,提高深部电阻率的分辨率. 相似文献
15.
In this study, we apply Tikhonov’s regularization algorithm for a 3-D density inversion from the gravity-gradiometry data. To reduce the non-uniqueness of the inverse solution (carried out without additional information from geological evidence), we implement the depth-weighting empirical function. However, the application of an empirical function in the inversion equation brings the bias problem of the regularization factor when a traditional Tikhonov’s algorithm is applied. To solve the bias problem of regularization factor selection, we present a standardized solution that comprises two parts for solving a 3-D constrained inversion equation, specifically the linear matrix transformation and Tikhonov’s regularization algorithm. Since traditional regularization techniques become numerically inefficient when dealing with large number of data, we further apply methods which include the Simultaneous Iterative Reconstruction Technique (SIRT) and the wavelet compression combined with Least Squares QR-decomposition (LSQR). In our simulation study, we demonstrate that SIRT as well as the wavelet compression plus LSQR algorithm improve the computation efficiency, while provide results which closely agree with that obtained from applying Tikhonov’s regularization. In particular, the algorithm of wavelet compression plus LSQR shows the best computing efficiency, because it combines the advantages of coefficients compression of big matrix and fast solution of sparse matrix. Similar findings are confirmed from the vertical gravity gradient data inversion for detecting potential deposits at the Kauring (near Perth, Western Australia) testing site. 相似文献
16.
Different from the stacked seismic data, pre-stack data includes abundant information about shear wave and density. Through inversing the shear wave and density information from the pre-stack data, we can determine oil-bearing properties from different incident angles. The state-of-the-art inversion methods obtain either low vertical resolution or lateral discontinuities. However, the practical reservoir generally has sharp discontinuities between different layers in vertically direction and is horizontally smooth. Towards obtaining the practical model, we present an inversion method based on the regularized amplitude-versus-incidence angle (AVA) data to estimate the piecewise-smooth model from pre-stack seismic data. This method considers subsurface stratum as a combination of two parts: a piecewise smooth part and a constant part. To fix the ill-posedness in the inversion, we adopt four terms to define the AVA inversion misfit function: the data misfit itself, a total variation regularization term acting as a sparsing operator for the piecewise constant part, a Tikhonov regularization term acting as a smoothing operator for the smooth part, and the last term to smoothly incorporate a priori information for constraining the magnitude of the estimated model. The proposed method not only can incorporate structure information and a priori model constraint, but also is able to derive into a convex objective function that can be easily minimized using iterative approach. Compared with inversion results of TV and Tikhonov regularization methods, the inverted P-wave velocity, S-wave velocity and density of the proposed method can better delineate the piecewise-smooth characteristic of strata. 相似文献
17.
重磁电震等地球物理勘探成果的联合反演是解决复杂地质问题所用的综合地球物理解释方法,是目前可信的定量的综合地球物理解释技术,是联合应用多种地球物理信息,通过反演地质体的岩石物性和几何参数来求得同一个地下地质和地球物理模型的技术.本文综述了从上世纪80年代以来,地球物理场联合反演技术取得的长足发展,主要表现在两个方面:一是在地球物理场综合应用基础上,发展了两两地球物理参数之间互为约束的反演的技术,二是是反映了联合反演计算方法的进步,如线性与非线性反演比较、反演方法选择等.综述了联合反演的思路与对策,提出了联合反演的发展方向. 相似文献
18.
Inversion for seismic impedance is an ill-posed and band-limited problem so that inversion results are non-unique and unstable and low and high frequency components of inversion results are missed. Combining regularization with constraints of sonic log data and geological structure information can help to alleviate these problems. To achieve this, we developed an inversion method by constructing a new objective function which includes edge-preserving regularization and soft constraint based on Markov random field (MRF). The method directly introduces absolute constraints with prior impedance and sonic log data in the objective function and indirectly achieves relative constraints with geologic structures of layer interfaces and faults by adjusting the regularization parameter which is the scaling parameter δ. Moreover, we improved the inversion result using anisotropic diffusion smoothing method. Optimization approach utilized in inversion is a fast simulated annealing (FSA). We test the method on both synthetic and field data examples. Tests on 2-D synthetic data indicate that aspects of the discontinuity in the inversion results are significantly improved by adding δ values in faults and layer interfaces. We obtained better results by combining the first-order neighborhood and the third-order neighborhood of MRF. The inversion results of the field data provide more detailed information of the layers. The results of nearby faults were improved by introducing the geological structure constraints. 相似文献
19.
Techniques that incorporate regularization in space and time have been proposed to reduce inversion artifacts that may lead a misinterpretation of geophysical monitoring data. Applying this time regularization, however, may result in a model too smoothly carrying in the time domain. To alleviate this problem, we propose an algorithm for inverting time-lapse resistivity monitoring data.Here the time regularization is not considered to be constant between different time steps but is now allowed to vary depending on the degree of spatial resistivity changes occurring between different monitoring stages. Two methods are proposed to assign different time Lagrangian values, one based on a pre-estimation during execution time, and one using a-priori information. Both methods require a threshold to characterize the significance of the observed resistivity changes with time. We performed numerous numerical experiments using synthetic data to provide reasonable threshold values. Synthetic data tests illustrate that the new algorithm, named 4D Active Time Constrained (4D-ATC), produces in most cases improved time-lapse images when compared with existing techniques. Further the applicability of the new scheme is demonstrated with real data. Overall, the new algorithm is shown to be a useful tool for processing time-lapse resistivity data, which can be used with minor modifications to other types of time-lapse geophysical data. 相似文献
20.
常规三维大地电磁反演的正则项为L2范数,它以电阻率空间分布函数处处光滑为模型期望,弱化了算法对电性突变界面的分辨能力.本文实现了正则项为L1范数的三维大地电磁反演算法,让模型空间梯度向量更有机会取得稀疏解,在充分正则的迭代下能够有效突出模型真实电性界面.为避免L1范数零点不可导带来的求解困难,使用迭代重加权最小二乘法把原问题转换为一系列L2正则子问题迭代求解.每个子问题的极小方法使用改进型拟牛顿法,其下降方向既能保证正则项海塞矩阵的精确性,又能允许反演过程随迭代灵活更新正则因子.使用比值法或分段衰减法自适应更新正则因子以避免迭代早期陷入奇异解,从而提升反演收敛的稳定性并降低初始模型依赖度.合成的无噪数据反演表明L1正则算法的模型恢复效果优于L2正则;不同噪声水平的合成数据反演表明本文的算法具有稳健性;实测数据反演对比表明在合理的正则因子调整策略下,L1正则反演结果的模型分辨率优于L2正则.另外,不同初始模型的反演测试还表明,正则因子选取不合理时L1正则可能造成方块状假异常. 相似文献