共查询到20条相似文献，搜索用时 31 毫秒
1.
磁异常的反演是地球物理勘探的重要手段,三维磁化率反演是磁异常定量解释中的一种重要方法.由于剩磁的存在使得磁化方向与地磁场方向产生偏差,从而影响了磁异常反演与解释的精度.本文基于磁异常模量反演和磁化强度矢量反演方法得到了一种新的磁化强度矢量反演方法.与以往的磁化强度矢量反演方法相比,该方法以磁异常模量反演得到的磁化率模型为约束,采用Lp范数正则化方法求解,提高了磁化强度矢量反演的精度和效率.本文通过模拟试验的反演计算,验证了这种磁化强度矢量反演方法的有效性.最后,将本文方法应用于新疆东天山卡拉塔格地区航磁数据的解释,获得了地下空间不同磁性差异的磁性体的空间分布特征,为进一步分析研究区隐伏矿床提供了重要信息. 相似文献
2.
3.
强剩磁的存在通常导致了总磁化强度方向未知,进而影响了磁异常的反演和解释.磁异常模量是一种受磁化方向影响小的转换量,可以在强剩磁条件下通过反演三维磁化强度大小分布来推测场源分布状态.我们提出了一种数据空间磁异常模量反演算法来减少剩磁的影响.与标准的模型空间L2范数正则化反演方法相比,我们的方法有两个优点:一是无需搜索正则化参数(需要反复求解非线性反演问题),因而可以减少计算时间;二是反演结果更加聚焦,深度分辨率更高,我们对此进行了原因分析.通过模型和实测数据测试证明了该算法的有效性和更好的反演效果. 相似文献
4.
3D stochastic inversion of magnetic data 总被引：1，自引：0，他引：1
5.
三维反演是磁测数据定量解释的重要方法,在金属矿勘探中扮演着重要的角色.但是在实际矿区的应用中,传统的磁总场异常反演方法依然存在两个问题：一是地面磁异常反演的深度分辨率较低,深部场源体的成像效果差;二是金属矿中可能包含强剩磁,反演结果可能是完全错误的.尽管前人对上述两个问题分别进行了广泛的研究,但尚未尝试同时解决这两个问题.本文在前人研究的基础上,提出了一种井地磁异常模量联合反演方法,该方法需要的控制参数少,无需加入额外的地质信息,且可用于多场源复杂磁异常的反演,具有较强的适用性.本文方法首先将地面和井中磁异常转化为模量数据,然后利用基于核函数或距离的加权函数将井地模量数据结合起来,使得该方法适用于联合反演.我们利用井地多种异常参量进行反演的模型试验表明,在强剩磁存在时,本文方法的效果优于其他方法,在减少剩磁影响的同时,也改善了深部成像效果,具有良好的应用前景. 相似文献
6.
随着重力和重力梯度测量技术的日趋成熟,基于重力和重力梯度数据的反演技术得到了广泛关注.针对反演多解性严重、计算效率低和内存消耗大等难点问题,本文开展了三维重力和重力梯度数据的联合反演研究,该方法结合重力和重力梯度两种数据,将L0范数正则化项加入到目标函数中,并在数据空间下采用改进的共轭梯度算法求解反演最优化问题.同时,本文摒弃了依赖先验信息的深度加权函数,引入了自适应模型积分灵敏度矩阵,用来克服因重力和重力梯度数据核函数随深度增加而衰减引起的趋肤效应问题.为了提高反演计算效率,本文又推导出基于规则网格化的重力和重力梯度快速正演计算方法.模拟试算表明,改进的共轭梯度法可以降低反演的迭代次数,提高反演的收敛速度;自适应模型积分灵敏度矩阵,可以有效解决趋肤效应,提高反演纵向分辨能力;数据空间和改进的共轭梯度算法结合,可以更好地降低反演求解方程的维度,避免存储灵敏度矩阵,有效地降低反演计算时间和内存消耗量.野外实例表明,该算法可以在普通计算机下快速地获得地下密度分布模型,表现出较强的稳定性和适用性. 相似文献
7.
We develop an inversion procedure using the total variation (TV) regularization method as a stabilizing function to invert surface gravity data to retrieve 3-D density models of geologic structures with sharp boundaries. The developed inversion procedure combines several effective algorithms to solve the TV regularized problem. First, a matrix form of the gradient vector is designed using the Kronecker product to numerically approximate the 3-D TV function. The piecewise polynomial truncated singular value decomposition (PP-TSVD) algorithm is then used to solve the TV regularized inverse problem. To obtain a density model with depth resolution, we use a sensitivity-based depth weighting function. Finally, we apply the Genetic Algorithm (GA) to select the best combination of the PP-TSVD algorithm and the depth weighting function parameters. 3-D simulations conducted with synthetic data show that this approach produces sub-surface images in which the structures are well separated in terms of sharp boundaries, without the need of a priori detailed density model. The method applied to a real dataset from a micro-gravimetry survey of Gotvand Dam, southwestern Iran, clearly delineates subsurface cavities starting from a depth of 40 m within the area of the dam reservoir. 相似文献
8.
Interactive 3D forward modeling of total field surface and three-component borehole magnetic data for the Daye iron-ore deposit (Central China) 总被引：1，自引：0，他引：1
Yushan Yang Yuanyuan Li Tianyou Liu Yinglin Zhan Jie Feng 《Journal of Applied Geophysics》2011,75(2):254-263
The ground magnetic response of deep ore bodies in the Daye iron-ore deposit is relatively weak, and sometimes concealed by the strong magnetic background of shallower sources. Apart from the low-quality ground magnetic data, another critical problem for reconstructing the deep skarn-type ore bodies is developing a versatile inversion scheme that can simultaneously resolve 3D sources with arbitrary shapes. In this case, we resort to interactive 3D forward modeling solution with the joint use of two data sets-total field surface and three-component borehole magnetic data. Joint inversion of the two data sets is expected to help resolve the ambiguity associated with either data set and greatly reduces the nonuniqueness of the magnetic inversion. Such nonuniqueness is especially severe when a 3-D distribution of magnetic susceptibility, instead of a simple body, is sought from the inversion.In this paper, we calculate the magnetic field on the surface and in the borehole caused by 3D arbitrarily-shaped bodies with the triple integral method. The complex 3D magnetic sources having arbitrary shapes are constructed with cross-sections, termination points and facets in our visualization technology. We specify, interactively and in a user-friendly environment, the outline of the sources in terms of geometric elements and their magnetic parameters. The method automatically fits the observations within a prescribed precision. If dissatisfied, the user can redefine the model parameters and proceed to a new inversion. The method's ability to interpret a complicated 3D geologic environment is demonstrated on synthetic models and real data profiles in the Daye iron-ore deposit in central China. The interactive forward modeling results in all tests demonstrate a good correlation of estimated magnetic sources with corresponding known geologic features. 相似文献
9.
密度界面反演作为了解地球内部结构的一种重要方法,长期以来都是重力学研究的主要内容.本文结合抛物线密度模型及频率域算法的优点,将抛物线密度函数应用于Parker-Oldenburg算法,经过理论推导得到了抛物线密度模型的频率域公式,从而建立了基于抛物线密度模型的三维密度界面重力异常正反演的算法和流程.理论模型数据试验表明本方法快速、有效,适用于大多数浅部比深部增加更快的实际地壳密度.研究中还利用该方法对川滇地区重力异常进行了反演,获得了该区的莫霍面深度分布,并与接收函数研究结果进行对比分析,进一步验证了本文方法的正确性和有效性. 相似文献
10.
İlknur Kaftan 《Acta Geophysica》2017,65(4):627-634
A genetic algorithm (GA) is an artificial intelligence method used for optimization. We applied a GA to the inversion of magnetic anomalies over a thick dike. Inversion of nonlinear geophysical problems using a GA has advantages because it does not require model gradients or well-defined initial model parameters. The evolution process consists of selection, crossover, and mutation genetic operators that look for the best fit to the observed data and a solution consisting of plausible compact sources. The efficiency of a GA on both synthetic and real magnetic anomalies of dikes by estimating model parameters, such as depth to the top of the dike (H), the half-width of the dike (B), the distance from the origin to the reference point (D), the dip of the thick dike (δ), and the susceptibility contrast (k), has been shown. For the synthetic anomaly case, it has been considered for both noise-free and noisy magnetic data. In the real case, the vertical magnetic anomaly from the Pima copper mine in Arizona, USA, and the vertical magnetic anomaly in the Bayburt–Sar?han skarn zone in northeastern Turkey have been inverted and interpreted. We compared the estimated parameters with the results of conventional inversion methods used in previous studies. We can conclude that the GA method used in this study is a useful tool for evaluating magnetic anomalies for dike models. 相似文献
11.
通过研究重力数据三维反演解的病态性,利用基于拉格朗日插值方法的Extrapolation Tikhonov正则化方法来解决反演中解的不唯一性和不稳定性问题,该方法最大限度的减少了因正则化参数的引入而在反演结果中介入的误差,同时详细讨论了基于三种选择原则的正则化双参数的具体选择方法,模型试算结果表明,与原Tikhonov方法相比,该方法提高了反演的拟合精度.其次,为了消除核函数随深度增加而快速衰减对反演结果的影响,本文改进了前人的重力数据三维反演深度加权函数,改进后的加权函数与原函数相比能更好的识别异常体底部密度分布特征,对于埋深较深的异常体具有较好的识别效果,更好的解决了由近地面趋肤效应作用引起的密度分布不均的问题.同时,利用上下限约束函数限制每一个立方体的密度差范围,并应用于多组人工合成模型.结果表明:该反演方法能准确地获得正演模型的预设参数范围和位置. 相似文献
12.
以等效源及位场物性反演为基础,本文提出一种新的求取重磁场源深度的方法.该方法将一层等效源以一定的间隔从浅部向深部移动,并将等效源作为初始模型进行反演,当反演拟合差最小时,停止反演,此时的等效源底深即为所求场源的中心深度.由于仅需要反演一层等效源,比传统的物性反演计算时间大大减少,并且不需要进行深度加权约束.理论模型数据处理结果表明该方法能够获得较准确的场源深度:以长宽比为7.5的薄板模型为例,深度计算误差约为1个点距(25 m);以长宽比为0.5~1.5的厚板模型为例,深度计算误差小于1个点距(25m).将该方法应用于实测航磁梯度数据,计算的磁源中心深度在200~250m之间,钻井资料显示该异常由埋藏深度在200~300m的闪长岩引起,计算结果与钻井资料较吻合. 相似文献
13.
基于光滑约束的最小二乘法是三维电阻率反演的主要方法,但该方法在某些情况下存在着多解性较强的问题,且普遍耗时较长,严重制约了三维反演方法的推广与发展.为改善上述问题,将表征模型参数变化范围的不等式约束作为先验信息引入最小二乘线性反演方法中,有效地改善了反演结果的精度,降低了反演的多解性问题.为了解决耗时较长的问题,基于预条件共轭梯度(PCG)算法和Cholesky分解法的特点提出了一套优化三维电阻率反演计算效率的计算方案.在该方案中,Cholesky分解法被用来求解敏感度矩阵计算中的多个点源场的正演问题,Cholesky分解法只需对总体系数矩阵进行一次分解,然后对不同的右端向量进行回代即可.将预条件共轭梯度法引入到三维电阻率反演方程的求解中,将雅可比迭代中的对角阵作为预处理矩阵,其具有求逆方便、无需内存空间的特点,有效地加快了收敛速度.对合成数据以及实测数据的反演算例表明,借助不等式约束和反演效率优化方案,最小二乘反演方法可得到较为精确的反演结果,有效地提高了反演计算效率,具有良好的推广前景. 相似文献
14.
为了降低单一地球物理方法反演的多解性及受噪声的影响程度,本文围绕重力、磁法和大地电磁法开展了三维联合反演的研究.重、磁采用基于对数障碍法的正则化反演算法,大地电磁使用limited-memory BroydenFletcher-Goldfarb-Shanno(L-BFGS)反演算法,引入交叉梯度函数实现了三种物性结构的相互耦合,最终开发出一套重磁电三维联合反演算法,并实现MPI并行加速计算.通过理论模型算例验证了算法的准确性,结果表明:不论是单棱柱体模型还是组合棱柱体模型,联合反演结果相较单独反演对于异常体的空间形态刻画以及物性数值恢复具有较好的提升;单棱柱体模型算例使得异常体的物性参数(密度、磁化率和电阻率)更加接近于真实的物性参数;组合棱柱体模型的联合反演结果不仅仅消除了围岩物性参数的假异常,而且还增强了异常体边界结构的恢复程度. 相似文献
15.
研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法(BICG)计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次“拟正演”问题,进而将PARDISO的快速正演方案应用于“拟正演”问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性. 相似文献
16.
本文介绍了多震相的层析成像的思路和算法，通过穿透和反射走时可以同时作出２维和３维慢度（速度的倒数）重建。我们分析了在穿透和反射数据中确定速度和深度的不确定性，并认识到深度扰动对反射走时异常比慢度扰动更敏感。由不同波类型所提供的对速度和深度的约束，这个算法实际上减少了在一般反射层析成像在速度和深度之间的不确定性，并且也避免了在穿透层析成像中的不确定问题。线性化反演是通过从反射界面深度由分离速度参数迭代进行的。使用一个快速的２－Ｄ和３－Ｄ射线跟踪算法来计算穿透和反射走时和对幔度及反射界面深度的偏导数。深度和速度都用立方Ｂ样条函数来进行参数化。合成例子表明，当同时考虑穿透和反射时间，层析成像的结果得到改进。这个方法也应用到英国煤炭测量局（ＢｒｉｔｉｓｈＣｏａｌＭｅａｓｕｒｅｓ）沿跨线排列所记录的逆ＶＳＰ数据组。通过使用波形配合技术，用同时确定时间延迟和叠加权，可以自动拾取旅行时间。所观察到的逆ＶＳＰ层析成像可比周围介质具有较低速度的两个断层区域成像。断层的位置由附近的反射测线所确定。本文还讨论了在复杂２－Ｄ和３－Ｄ非均匀各向同性介质中地震射线跟踪方法。界面的几何形状和水平速度场都通过使用非均匀步长立方Ｂ－样条节点 相似文献
17.
L. B
RSTING PEDERSEN 《Geophysical Prospecting》1977,25(2):199-230
The interpretation of potential field data from two-dimensional structures with a single interface of density or susceptibility contrasts is solved in terms of generalized matrix inversion. The model equations are derived, and important features of generalized matrix inversion are treated. The method is subsequently used to solve two gravity problems, an artificial one where the solution is known, and a geophysical one related to the crust-mantle interface. The solution is shown to compare well with the FFT results of Oldenburg. The method is also used to solve two magnetic problems, an artificial one with the solution known, and a geophysical one from the continental shelf of Greenland. The advantages and limitations of the method are finally discussed. 相似文献
18.
19.
In this paper, we describe a non‐linear constrained inversion technique for 2D interpretation of high resolution magnetic field data along flight lines using a simple dike model. We first estimate the strike direction of a quasi 2D structure based on the eigenvector corresponding to the minimum eigenvalue of the pseudogravity gradient tensor derived from gridded, low‐pass filtered magnetic field anomalies, assuming that the magnetization direction is known. Then the measured magnetic field can be transformed into the strike coordinate system and all magnetic dike parameters – horizontal position, depth to the top, dip angle, width and susceptibility contrast – can be estimated by non‐linear least squares inversion of the high resolution magnetic field data along the flight lines. We use the Levenberg‐Marquardt algorithm together with the trust‐region‐reflective method enabling users to define inequality constraints on model parameters such that the estimated parameters are always in a trust region. Assuming that the maximum of the calculated gzz (vertical gradient of the pseudogravity field) is approximately located above the causative body, data points enclosed by a window, along the profile, centred at the maximum of gzz are used in the inversion scheme for estimating the dike parameters. The size of the window is increased until it exceeds a predefined limit. Then the solution corresponding to the minimum data fit error is chosen as the most reliable one. Using synthetic data we study the effect of random noise and interfering sources on the estimated models and we apply our method to a new aeromagnetic data set from the Särna area, west central Sweden including constraints from laboratory measurements on rock samples from the area. 相似文献
20.
计算密度分界面的起伏变化在区域地质构造研究和石油矿产资源勘探中具有重要意义.已有密度界面反演方法更多侧重约束变密度界面反演算法,而对约束信息的准确性、研究区横向密度变化往往考虑不足,影响了最终反演结果的可信度.本文在变密度界面正反演算法基础之上,结合实际需求,提出已知深度信息约束下的变密度界面反演策略.该策略主要包括变密度约束反演算法、已知深度约束信息校验和分区变密度模型三个方面.其中反演算法提供了带已知深度约束信息的密度界面迭代反演方法;约束信息校验用于评估约束信息精度,通过调差降低约束信息的系统误差;在反演过程中引入水平密度分区以应对不同构造背景密度界面模型,提高反演结果的可信度.最后将本策略应用于南海莫霍面深度反演计算中,结果显示借助已知约束信息,利用分区密度模型能够获得更为可信的深度反演结果,验证了该策略的正确性. 相似文献