共查询到19条相似文献,搜索用时 468 毫秒
1.
三维高密度电阻率法可以有效地解决二维高密度电阻率法在探测空间形态上的不足。将高密度电阻率法二维勘探数据做三维反演,不仅实现了地下空间结构的三维电阻率成像,还克服了传统三维高密度电阻率法在经济成本上的限制。将此方法应用到岩溶勘探中,利用DUK-2A高密度电法仪系统采集数据;然后,把二维高密度数据拼接转化成三维可反演数据,再将反演得到的结果通过三维成像软件实现任意三维地电切片,从而获得丰富的地电结构信息;最后,结合三维切片与钻孔资料对溶洞进行解释,为工程的顺利施工提供科学依据,从而证明该方法的可行性和有效性。 相似文献
2.
Zohdy方法在三维电阻率测深数据反演中的应用 总被引:2,自引:0,他引:2
三维电阻率测深数据采用Zohdy方法近似反演,正演计算采用有限单元法,反演初始模型由测量视电阻率数据给出,通过比较实测视电阻率值和预测模型计算的视电阻率值对数差来修改模型网格电阻率.为了解决任意电极距测深数据的反演,采用大、小双网格剖分,大网格反映地下电性分布情况,小网格用于实际有限元正演计算.文章对加5%随机噪声的模型理论视电阻率测深断面数据和实测视电阻率测深断面数据进行反演.结果表明该反演方法需要的内存小、计算量少,并能够取得很好的反演效果. 相似文献
3.
4.
有限单元法三维电阻率最小二乘反演中存在问题的研究 总被引:6,自引:1,他引:5
针对有限单元法在三维电阻率最小二乘反演中关于深部单元的反演能力、反演稳定性和速度问题进行了讨论与改进.首先构造了全新的"体积因子"作为先验信息加入到电阻率反演中,使地质体赋存区域的电阻率变化更符合实际情况,提高了深度较大的地质体的分辨能力.反演计算中对实测电阻率数据和模型正演数据分别取对数,并改进迭代修正量的计算方法,将以前的加减修正量改为与修正量指数之积,这样,不但可以改进电阻率模型出现负值的缺陷,提高反演计算的稳定性,而且还可减少迭代次数;针对反演时间过长的缺陷,将三维全域反演改为三维局域反演,使每次迭代的时间由4 h减少到目前的15 min;对模型的反演计算表明,以上改进是非常有效的. 相似文献
5.
从高密度电阻率成像法到三维空间上的包气带水文学 总被引:5,自引:0,他引:5
周启友 《水文地质工程地质》2003,30(6):97-104
高密度电阻率成像法是新近出现的通过测定介质电传导特性的地球物理勘探方法。它方便、快速和多尺度,可对探测对象进行无损的、三维的和动态的监测。由于包气带中发生的各种物理、化学、生物和水文过程不可避免的都会引起介质电阻率的变化,通过高密度电阻率成像法可以对这些过程进行三维的动态监测,从而为我们在三维空间上认识各种包气带过程创造了条件。本文首先介绍了高密度电阻率成像法在半无限的介质表面上应用时的电极布置、测定和反演计算方法,然后介绍了影响介质电阻率的因素。最后,从电阻率空间分布的差异性和时间上的变化性方面,介绍了目前高密度电阻率成像法的应用领域和研究方向,以期能促进高密度电阻率成像法的广泛应用和推动三维空间上的包气带水文学的研究。 相似文献
6.
7.
基于积分方程法的大地电磁三维反演 总被引:1,自引:0,他引:1
采用积分方程法计算均匀半空间中三维异常体的大地电磁响应,然后利用阻尼最小二乘法反演得到地下三维异常体电阻率分布。在计算灵敏度矩阵时,将灵敏度矩阵分为线性项和体现反演过程非线性的偏微分项进行求解,提高了计算精度。理论模型反演结果表明,反演较准确,对初值依赖性小,且只针对异常体剖分单元进行反演,反演效率较高。 相似文献
8.
9.
为了减小地球物理反演的多解性,通常采用累加型目标函数,即数据拟合差项加上正则化项。基于累加型目标函数的反演存在如何优选正则化因子的问题,这个过程通常需要做大量的反演计算。本研究系统分析和实现了基于乘积型目标函数(数据拟合差项和正则化项相乘)的电阻率法三维反演,乘积型目标函数的反演不存在优选正则化因子的问题。三维正演采用非结构化有限单元法,三维反演采用有限内存拟牛顿方法。使用理论模型合成数据进行了三维反演,检验了基于乘积型目标函数的电阻率法三维反演的可行性与有效性。反演算例结果表明:基于乘积型目标函数的反演方法能够可靠恢复异常体的电阻率值、形态和位置。 相似文献
10.
11.
在三维电阻率的正演计算中往往涉及到快速、准确求解大型线性方程纽Ax=b的问题。通过采用有限差分法来构造出求解点电源三维地电场的大型稀疏对称线性方程组。并引入Lanczos迭代技术,构造出三对角阵方程组,然后采用正交分解法进行求解,它是Krylov子空间方法中的一种。与传统迭代算法相比,它占用内存少,收敛速度快且稳定。针对大型稀疏矩阵及MATLAB语言的特点,采用简单记录矩阵的非零元素值及其所在行、列值的方法存储大型稀疏矩阵,可大大节省机器内存,提高运算速度。理论分析和计算实例显示,此算法是地电三维正演计算的有效方法,为下一步的反演计算打好基础。 相似文献
12.
直流电阻率法因效率高、成本低等特点广泛应用于有色金属、煤田等矿产资源勘探和地质调查行业。但在井–地、地–井电阻率法勘探过程中,钻孔因素(井液电阻率和钻孔孔径)对不同探测模式下视电阻率响应的影响尚未明确,是否影响视电阻率资料的解释是值得探讨的一个问题。因此,提出采用混合网格有限元法实现直流电阻率三维正演。给出异常电位法满足的边值问题及有限元变分问题,应用三棱柱和四面体混合网格实现对计算区域的快速离散,并建立2种网格的线性插值基函数和单元系数矩阵;采用SSOR–PCG迭代算法求解异常电位满足的大型线性方程组,得到各观测点的响应。在保证计算精度的前提下,应用混合网格有效地离散钻孔地电模型,探讨钻孔因素对井–地、地–井观测方式视电阻率数据的影响特征。对于井–地观测方式:钻孔附近视电阻率受钻孔因素影响最大,严重影响视电阻率资料的合理解释,但随着测点收发距的增加,视电阻率响应逐渐趋于围岩电阻率;而且随着发射源深度的增加,钻孔对地表视电阻率响应的影响也逐渐减小。对于地–井观测方式:钻孔因素对浅部数据影响较大,对深部数据影响小,相比于井液电阻率,视电阻率响应更易受钻孔孔径的影响。直流电阻率法混合网格有限元三维正演算法的提出,对实际井–地、地–井电阻率法勘探具有理论指导意义,可结合钻孔信息进行正演模拟分析,选择合适的收发距,有效压制钻孔因素对实测资料的影响。 相似文献
13.
Randomized maximum likelihood is known in the petroleum reservoir community as a Bayesian history matching technique by means of minimizing a stochastic quadratic objective function. The algorithm is well established and has shown promising results in several applications. For linear models with linear observation operator, the algorithm samples the posterior density accurately. To improve the sampling for nonlinear models, we introduce a generalized version in its simplest form by re-weighting the prior. The weight term is motivated by a sufficiency condition on the expected gradient of the objective function. Recently, an ensemble version of the algorithm was proposed which can be implemented with any simulator. Unfortunately, the method has some practical implementation issues due to computation of low rank pseudo inverse matrices and in practice only the data mismatch part of the objective function is maintained. Here, we take advantage of the fact that the measurement space is often much smaller than the parameter space and project the prior uncertainty from the parameter space to the measurement space to avoid over fitting of data. The proposed algorithms show good performance on synthetic test cases including a 2D reservoir model. 相似文献
14.
15.
利用ICCG迭代技术加快电阻率三维正演计算 总被引:2,自引:2,他引:0
一般而言,有限差分法求解点源三维地电场正问题所形成的大型稀疏线性方程组Ax=b,直接解法的计算效率极低。本文从系数矩阵A的不完全Cholesky分解及矩阵特征值的特点等角度,说明了不完全Cholesky共轭梯度(ICCG)迭代技术可大大提高电阻率三维正演速度的内在原因。结合矩阵A的稀疏存储模式,使得内存需求也大大减少。 相似文献
16.
电阻率三维复杂结构的快速反演 总被引:3,自引:1,他引:2
运用共轭梯度迭代算法解三维反问题中的线性方程组,并结合求Jacobi矩阵 G 的Rodi算法,则每次反演迭代仅需一次正演计算,大大加快了计算速度,实现了直流电阻率三维快速反演。另外,由于避免了存储 G 和 G T G 所需的庞大存储量,以及在三维反演中加入光滑约束,有利于精细网格下的复杂模型反演。 相似文献
17.
从目标函数的线性化出发,采用奇异值分解法求广义逆矩阵并适当加阻尼,用迭代法直接解超定方程实现改进广义逆矩阵反演。进一步从理论与实际结合上阐明改进广义逆矩阵反演法的灵活、稳定,它涵盖了多种反演方法的优点,又具有自身特点,在实际应中适应性强,并能提供一些辅助信息,可更好地评价解释结果。 相似文献
18.
3D behaviour of contamination in landfill sites using 2D resistivity/IP imaging: case studies in Portugal 总被引:2,自引:0,他引:2
Aiming at defining a valid spatial contamination model, resistivity and induced polarization (IP) measurements were used to
investigate contamination plumes in the vicinity of two municipal landfills (Ovar and ílhavo). Previous geophysical surveys
and underground water samples confirmed the contamination. However 2D resistivity/IP surveys enabled in obtaining a more accurate
spatial model. The Ovar survey consisted of two profiles with nine Wenner soundings each; the ílhavo survey was carried out
along two individual lines using a Wenner standard pseudo-section. In both situations, negative IP values were found associated
with positive IP values, which can be explained mainly by 2D or 3D geometric effects caused by the presence of the conductive
plumes. The data were modelled using a 2D inversion program (RES2DINV) and the resulting resistivity and chargeability distributions
were displayed as pseudo-sections. The resistivity and chargeability pseudo-sections define the contamination plumes and the
sedimentary structure. These case studies illustrate the advantages of 2D resistivity/IP surveys for the mapping of shape
and dimension of contamination associated with landfills. 相似文献
19.
高精度、高效率的深海勘探是当前国家海洋资源勘探开发的研究热点之一,设计观测效率高、勘探成本低的探测方式,开发可精确计算海底复杂电性环境的正演模拟算法有助于推进该研究的进展。本文基于前人的研究经验,引入了全拖曳式深海直流电阻率观测方式,同时考虑到海底沉积环境的电各向异性,开发了基于非结构有限元方法的三维任意各向异性深海直流电阻率正演算法,实现了对海底任意电各向异性情况的仿真模拟。本文算法与层状一维模型半解析解的对比验证了算法的精度。对典型各向异性电性模型的仿真模拟和分析,证明了全拖曳式深海探测方式对海底沉积层和矿产资源电各向异性的高分辨率。在此基础上,总结了电导率各向异性主轴沿笛卡尔坐标旋转时视电阻率分布也发生同方向旋转的视电阻率分布特征,并根据海底地形对各向异性探测的影响进行了简单探讨。 相似文献