共查询到19条相似文献，搜索用时 125 毫秒
1.
地表起伏地形在野外矿产资源勘察中不可避免,其对直流电阻率法勘探影响巨大.近年来,电阻率三维正演取得诸多进展,特别是应用非结构网格我们能够进行任意复杂地形和几何模型的电阻率三维数值模拟,但面向实际应用的起伏地形下电阻率三维反演依然困难.本文基于非结构化四面体网格,并考虑到应用GPS/GNSS时,区域地球物理调查中可非规则布设测网的实际特点,实现了任意地形(平坦或起伏)条件下、任意布设的偶极-偶极视电阻率数据的不完全Gauss-Newton三维反演.合成数据的反演结果表明了方法的有效性,可应用于复杂野外环境下的三维电法勘探. 相似文献
2.
为了有效解决目前大地电磁和地震走时资料单方法反演结果一致性不好的问题,同时克服基于岩石不同物性参数间关系耦合约束联合反演的局限性,本文研究了基于交叉梯度耦合约束的大地电磁与地震走时资料的三维联合反演算法.以较为成熟的天然地震走时资料三维正反演和大地电磁三维正反演算法为基础,实现了具有共同的反演网格,以交叉梯度结构耦合约束,并能同时获得电阻率和速度模型的三维联合反演算法.分别利用单棱柱体模型和双棱柱体模型合成数据进行了联合反演试算.结果表明:无论是单棱柱体模型还是双棱柱体模型,联合反演结果比单独反演对异常体的空间形态都有更好的恢复,其中单棱柱体模型反演的异常体电阻率更接近于真实电阻率,双棱柱体模型的联合反演结果不仅消除了围岩的部分电阻率假异常,而且增强了对异常体深部速度结构特征的恢复程度.联合反演还能同时改善电阻率和速度反向变化异常体的单独反演结果,进一步证明交叉梯度耦合不依赖于岩石物性关系,而强调地下结构的相似性,具有更普遍的适用性. 相似文献
3.
为推进大地电磁三维反演的实用化,本文实现了基于L-BFGS算法的带地形大地电磁三维反演.首先推导了大地电磁法三维反演的Tikhonov正则化目标函数以及Hessian矩阵逆矩阵近似表达式和计算方法,然后设计了一种既能保证空气电阻率固定不变又能保证模型平滑约束的协方差矩阵统一表达式,解决带地形反演问题.在反演算法中采用正则化因子冷却法以及基于Wolf条件的步长搜索策略,提升了反演的稳定性.利用开发的算法对多个带地形地电模型（山峰地形下的单个异常模型、峰-谷地形下的棋盘模型）的合成数据进行了三维反演,并与已有大地电磁三维反演程序（ModEM）进行对比,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了华南某山区大地电磁实测数据,得到该区三维电性结构,揭示了研究区以高阻介质为基底,中间以低阻不整合面和相对低阻介质连续分布,浅部覆盖高阻介质的电性结构特征,进一步验证了本文算法的实用性. 相似文献
4.
系统地论述了用有限单元法研究复杂地形条件下三维直流电阻率的正演计算技术.首先给出了三维构造中点源电场的边值问题以及相应的变分问题;然后利用有限单元法求解变分问题,采用四面体单元对研究区域进行剖分,在单元中进行三线性函数插值,将变分方程化为线性代数方程组;最后,考虑到节约计算时间,利用对称超松弛顸条件共轭梯度迭代算法求解大型线性方程组,得到了各节点的电位值,进而计算出地表的视电阻率.通过理论模型的计算检验了算法的可行性之后,给出了几种常见纯地形异常的数值模拟结果和一个组合模型的计算结果,其研究工作为研究三维直流电阻率反演奠定了基础. 相似文献
5.
在坑道电阻率超前探测工作中，视电阻率异常极值位置与异常体位置、厚度、电性以及产状相关.设计有效的反演方法对异常体的各种参数综合反演是非常必要的.首先利用有限元数值模拟技术分析并消除坑道空腔对视电阻率的影响，并将坑道超前探测视电阻率曲线转化为全空间视电阻率曲线.然后采用最小二乘法对超前探测视电阻率曲线进行快速反演.反演结果可准确反映围岩电阻率、异常体的电性、位置和产状；而由于同层等值现象，反演结果只能确定异常体的厚度和电阻率的组合值. 相似文献
6.
7.
8.
为了提高对大地电磁三维反演结果的分析和解释能力,完成构造界面识别、异常构造刻画等地质和地球物理解释,本文提出使用非监督的聚类方法分析大地电磁反演结果.根据三维反演模型电阻率值的分布和各种聚类分析方法的特点,选择使用K均值聚类方法对反演模型电阻率值进行聚类分析.在K均值聚类分析过程中,本文采用了RS指数指导选择聚类数目, Kaufman法进行中心初始化.通过对东北地区大地电磁数据三维反演结果使用K均值聚类分析方法,得到了东北地区电性岩石圈的厚度估计,结果表明东北地区岩石圈底部聚类电阻率大约为339Ω·m,其中松辽盆地岩石圈最薄,约为60 km;大兴安岭地区最厚,约为150 km;佳木斯地体厚度约为100 km;而长白山地区岩石圈厚度不易确定,可能受新生代构造活动影响,电阻率明显减小.聚类方法能够有效地帮助对大地电磁三维反演结果中的地质构造进行识别和归类. 相似文献
9.
CSAMT在山区金属矿勘查中,采用各种滤波和相位积分之类的处理方法,校正因地形起伏和局部电性不均匀引起的静态效应,往往难保奏效,开发消除静态效应的新方法是提高CSAMT资料处理与解释水平和方法应用效果的重要研究课题. 本文以如何消除地形影响为重点,对起伏地形下CSAMT二维大地三维源地电模型,采用加权余弦数值积分法,进行波数域电磁场二维有限单元法正演. 为模拟复杂地形地电模型,选取交叉对称网格三角形剖分法,实现了在国内常用赤道电偶极装置的CSAMT二维正演计算;在二维正演的基础上,开发了基于奥克姆反演法的CSAMT二维反演技术,研制出一套起伏地形下CSAMT二维正反演处理与解释方法技术系统. 通过理论模型试算和实测数据处理证实,本系统能有效地削减起伏地形影响. 在找矿应用中,该系统反演的电阻率断面,极大地消除了起伏地形影响和静态效应,突显出清晰的控矿构造和矿体的异常,取得了重要成效. 相似文献
10.
11.
以传统地震环境噪声面波成像方法研究地壳速度结构时,在一些极端的地形条件下,结果与真实结构会存在较大偏差.我们以地震波场三维正演模拟为基础,提出了一种地形校正方法.我们保留了传统噪声面波成像简单的两步反演法,在面波层析成像和一维速度结构反演的基础上,通过地震波场三维模拟近似估计地形和散射波场的影响,并据此校正瑞利波频散曲线,最终反演得到校正地形影响的S波速度结构.理论测试与在实际观测数据上的应用都证明了校正方法的有效性,同时也显示了地形校正的必要性. 相似文献
12.
Vikas C. Baranwal Antje Franke Ralph-Uwe Börner Klaus Spitzer 《Journal of Applied Geophysics》2011,75(2):363-372
We present a 2-D inversion code incorporating a damped least-squares and a minimum-model approach for plane wave electromagnetic (EM) methods using an adaptive unstructured grid finite element forward operator. Unstructured triangular grids permit efficient discretization of arbitrary 2-D model geometries and, hence, allow for modeling arbitrary topography. The inversion model is parameterized on a coarse parameter grid which constitutes a subset of the forward modeling grid. The mapping from parameter to forward modeling grid is obtained by adaptive mesh refinement. Sensitivities are determined by solving a modified sensitivity equation system arising from the derivative of the finite element equations with respect to the model parameters. Firstly, we demonstrate that surface topography may induce significant effects on the EM response and in the inversion result, and that it cannot be ignored when the scale length of topographic variations is in the order of magnitude of the skin depth. Secondly, the dependency of the inversion on the starting model is discussed for VLF and VLF-R data. Thirdly, we demonstrate the inversion of a synthetic data set obtained from a model with topography. Finally, the inversion approach is applied to field data collected in a region with undulating topography. 相似文献
13.
We describe an algorithm for inverting magnetotelluric (MT) data in the presence of strong bathymetry or topography. Instead of correcting distortions due to bathymetry or topography we incorporate them directly into the inversion. To achieve a high accuracy in computing MT responses we use finite‐difference approximations that permit sloping discrete boundaries inside model elements. The same approach is applicable to any seafloor electromagnetic calculation and can also be used to incorporate steep topography on land. We test our approach on various topographic features and compare our results to that of a finite element approach. Finally, we present inversion examples that illustrate the effectiveness of our inversion algorithm in recovering true subsurface structures in the presence of strong bathymetry and topography. 相似文献
14.
With the prospect of studying the relevance of the topographic core-mantle coupling to the variations of the Earth’s rotation and also its applicability to constraining the core surface flow, we investigate the variability of the topographic torque estimated by using core surface flow models accompanied by (a) uncertainty due to the non-uniqueness problem in the flow inversion, and (b) variance originating in that of geomagnetic secular variation models employed in the inversion. Various flow models and their variances are estimated by inverting prescribed geomagnetic models at the epoch 1980. The subsequent topographic torque is then calculated by using a core-mantle boundary topography model obtained by seismic tomography. The calculated axial and equatorial torques are found subject to the variability of order 1019 and 1020 Nm, respectively, on which (b) is more effective than (a). The variability of the torque is attributed even to (a) and (b) of the large-scale flows (degrees 2 and 3). Yet, it still seems unlikely for the decadal polar motion with the observed amplitude to be excited exclusively by the equatorial topographic torque associated with any of reasonable core surface flow models. It is also confirmed that, with the topography model adopted here, the axial topographic torque on a rigid annulus in the core (coaxial with the Earth’s rotation axis) associated with any of reasonable flow models is larger by two orders of magnitude than the plausible inertial torque on such cylinders. This implies that any core surface flow model consistent with the topographic coupling does not exist, unless the topography model is appropriately modified. Nevertheless, the topographic coupling might provide not only a weak constraint for explaining the decadal LOD variations, but also the possibility to probe the core surface flow and the core dynamics. 相似文献
15.
时间域航空电磁作为一种高效地球物理勘探技术特别适合我国地形复杂地区(沙漠、高山、湖泊、沼泽等)资源勘查.然而,这些地区地形起伏较大,对航空电磁响应有严重影响,忽略地形影响会给航空电磁数据解释造成很大误差.到目前为止人们对航空电磁地形效应特征研究十分有限.本文提出了基于非结构化网格的有限元法模拟带地形时间域航空电磁系统响应.该方法与基于结构化网格的有限差分相比能更好地模拟地形.首先通过傅里叶变换将2.5维问题转化成二维问题,利用伽辽金方法对二维问题进行离散.通过使用MUMPS求解器,得到波数域电磁响应.利用反傅里叶变换将波数域电磁响应变换到空间域,并利用正弦变换将其变换到时间域,得到2.5维时间域航空电磁响应.通过将本文的计算结果与半空间模型解析解及其他已发表的结果进行对比,检验了本文算法的精度.最后,我们系统分析了山峰和山谷地形对航空响应的影响特征.本文研究结果对航空电磁地形效应的识别和校正具有指导意义. 相似文献
16.
Geoelectrical profiling with multi-electrode systems has become an important tool for monitoring dike embankments bordering rivers. Profiles running perpendicular to the dike axis are affected by the dike topography, with the amplitude of this effect dependent on the surface geometry and the choice of the electrode configuration. Investigations using seven different electrode configurations have shown that some configurations are less sensitive to the topography than others.The topography correction method (TCM) is an important tool for processing data from measurements at river dikes. This method is generally recommended for flank angles steeper than 10°. The topography effect is calculated by two-dimensional finite element modelling. The resulting synthetic data of a homogeneous dike body are used to apply a topographic correction for each measurement.The topographic effect and correction procedure is demonstrated for synthetic dike data and for a data set from a river dike in Thai Binh province (Vietnam). The topography can be ignored for flank angles less than 25° if an averaged Half-Wenner electrode configuration is used. This configuration has proved to be less affected by undulated topography and the focusing effect of averaging the two data sets provides reliable structural information without the need for time-consuming data inversion. 相似文献
17.
Global ellipsoid-referenced topographic,bathymetric and stripping corrections to gravity disturbance
This paper revisits several aspects of defining and computing the anomalous gravity data for purposes of gravimetric inversion/interpretation.
Attention is paid to evaluation of a refined global topographic correction to the gravity disturbance based on the reference
ellipsoid (RE) and constant reference density for solid topography onshore and sea water density for liquid topography offshore.
The global bathymetric correction is discussed. Two issues associated with compilation and inversion of bathymetrically and
topographically corrected gravity disturbances in regions of negative ellipsoidal (geodetic) heights are pointed out: the
evaluation of normal gravity and the harmonic continuation of the gravity data. Stripping, the removal of an effect of a known
density contrast, is considered also for additional geological elements such as lakes, glaciers, sedimentary basins, isostatic
mountain roots, etc. The stripping corrections are discussed in the context of the gravimetric inverse problem. 相似文献
18.
Myeong-Jong Yi Jung-Ho Kim Yoonho Song Seong-Jun Cho Seung-Hwan Chung Jung-Hee Suh 《Geophysical Prospecting》2001,49(4):483-497
We have developed a three-dimensional inverse scheme for carrying out DC resistivity surveys, incorporating complicated topography as well as arbitrary electrode arrays. The algorithm is based on the finite-element approximation to the forward problem, so that the effect of topographic variation on the resistivity data is effectively evaluated and incorporated in the inversion. Furthermore, we have enhanced the resolving power of the inversion using the active constraint balancing method. Numerical verifications show that a correct earth image can be derived even when complicated topographic variation exists. By inverting the real field data acquired at a site for an underground sewage disposal plant, we obtained a reasonable image of the subsurface structures, which correlates well with the surface geology and drill log data. 相似文献