首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 468 毫秒
1.
三维高密度电阻率法可以有效地解决二维高密度电阻率法在探测空间形态上的不足。将高密度电阻率法二维勘探数据做三维反演,不仅实现了地下空间结构的三维电阻率成像,还克服了传统三维高密度电阻率法在经济成本上的限制。将此方法应用到岩溶勘探中,利用DUK-2A高密度电法仪系统采集数据;然后,把二维高密度数据拼接转化成三维可反演数据,再将反演得到的结果通过三维成像软件实现任意三维地电切片,从而获得丰富的地电结构信息;最后,结合三维切片与钻孔资料对溶洞进行解释,为工程的顺利施工提供科学依据,从而证明该方法的可行性和有效性。  相似文献   

2.
Zohdy方法在三维电阻率测深数据反演中的应用   总被引:2,自引:0,他引:2  
三维电阻率测深数据采用Zohdy方法近似反演,正演计算采用有限单元法,反演初始模型由测量视电阻率数据给出,通过比较实测视电阻率值和预测模型计算的视电阻率值对数差来修改模型网格电阻率.为了解决任意电极距测深数据的反演,采用大、小双网格剖分,大网格反映地下电性分布情况,小网格用于实际有限元正演计算.文章对加5%随机噪声的模型理论视电阻率测深断面数据和实测视电阻率测深断面数据进行反演.结果表明该反演方法需要的内存小、计算量少,并能够取得很好的反演效果.  相似文献   

3.
固定翼频率域航电测量效率高,空间分辨率高。以往我国三频航电实测数据已获得良好的解释成果,但均以定性解释和一维反演为主,二维反演应用较少。本文将快速近似反演、一维反演、二维反演引入到二连—东乌旗重点航电异常实测数据的解释中,并将三种方法的应用效果进行了对比。快速近似反演以视电阻率计算为基础,一维反演采用Brent方法,二维反演采用数据空间Occam反演,采用Rodi法求解雅克比矩阵。结果显示,快速近似反演可快速有效地提供大数据量地下电性分布信息,一维反演可提供稳定的、具一定深度的反演结果,二维反演则具有明显的横向、纵向分辨率优势。  相似文献   

4.
有限单元法三维电阻率最小二乘反演中存在问题的研究   总被引:6,自引:1,他引:5  
针对有限单元法在三维电阻率最小二乘反演中关于深部单元的反演能力、反演稳定性和速度问题进行了讨论与改进.首先构造了全新的"体积因子"作为先验信息加入到电阻率反演中,使地质体赋存区域的电阻率变化更符合实际情况,提高了深度较大的地质体的分辨能力.反演计算中对实测电阻率数据和模型正演数据分别取对数,并改进迭代修正量的计算方法,将以前的加减修正量改为与修正量指数之积,这样,不但可以改进电阻率模型出现负值的缺陷,提高反演计算的稳定性,而且还可减少迭代次数;针对反演时间过长的缺陷,将三维全域反演改为三维局域反演,使每次迭代的时间由4 h减少到目前的15 min;对模型的反演计算表明,以上改进是非常有效的.  相似文献   

5.
从高密度电阻率成像法到三维空间上的包气带水文学   总被引:5,自引:0,他引:5  
高密度电阻率成像法是新近出现的通过测定介质电传导特性的地球物理勘探方法。它方便、快速和多尺度,可对探测对象进行无损的、三维的和动态的监测。由于包气带中发生的各种物理、化学、生物和水文过程不可避免的都会引起介质电阻率的变化,通过高密度电阻率成像法可以对这些过程进行三维的动态监测,从而为我们在三维空间上认识各种包气带过程创造了条件。本文首先介绍了高密度电阻率成像法在半无限的介质表面上应用时的电极布置、测定和反演计算方法,然后介绍了影响介质电阻率的因素。最后,从电阻率空间分布的差异性和时间上的变化性方面,介绍了目前高密度电阻率成像法的应用领域和研究方向,以期能促进高密度电阻率成像法的广泛应用和推动三维空间上的包气带水文学的研究。  相似文献   

6.
文章简单描述了时域激电法原理、技术和测区地质概况及岩(矿)石电性特征.单极-偶极装置测深电阻率和充电率二维反演模型能够有效地圈定电性体在纵向和横向的分布特点,依据二维测深数据进行三维反演可获得电阻率和充电率异常曲线图,推测了不同深度范围内断裂和主要岩性分界,极化体空间分布位置和走向,钻探验证见铜矿化于划分的激电异常带,激电二维测深及其数据三维反演取得了较好的地质效果.  相似文献   

7.
基于积分方程法的大地电磁三维反演   总被引:1,自引:0,他引:1  
采用积分方程法计算均匀半空间中三维异常体的大地电磁响应,然后利用阻尼最小二乘法反演得到地下三维异常体电阻率分布。在计算灵敏度矩阵时,将灵敏度矩阵分为线性项和体现反演过程非线性的偏微分项进行求解,提高了计算精度。理论模型反演结果表明,反演较准确,对初值依赖性小,且只针对异常体剖分单元进行反演,反演效率较高。  相似文献   

8.
二维电阻率图像重建的一种新方法   总被引:13,自引:0,他引:13  
本文提出了一种新的二维电阻率成像方法,即电阻率层析成像.介绍了电阻率成像的几种供收方式,借鉴地震波和电磁波的反投影技术,应用电阻率反投影并联合使用常规算法,即导电α中心法作正演计算、非线性最小二乘法作反演修改,从而进行二维电阻率的图像重建.用Fortran语言编制了实现上述过程的程序,进行了大量数字模拟实验,获得较好的成像效果.  相似文献   

9.
为了减小地球物理反演的多解性,通常采用累加型目标函数,即数据拟合差项加上正则化项。基于累加型目标函数的反演存在如何优选正则化因子的问题,这个过程通常需要做大量的反演计算。本研究系统分析和实现了基于乘积型目标函数(数据拟合差项和正则化项相乘)的电阻率法三维反演,乘积型目标函数的反演不存在优选正则化因子的问题。三维正演采用非结构化有限单元法,三维反演采用有限内存拟牛顿方法。使用理论模型合成数据进行了三维反演,检验了基于乘积型目标函数的电阻率法三维反演的可行性与有效性。反演算例结果表明:基于乘积型目标函数的反演方法能够可靠恢复异常体的电阻率值、形态和位置。  相似文献   

10.
油气资源和矿产资源勘探中,常需要较大的勘探深度和高的深部分辨能力,充分利用现有钻孔的三维井地电阻率成像技术能够适应于这种要求。通过设计一个高阻水平板状体模型,利用Ansys有限元软件对模型进行三维井地电阻率法有限元正演计算,并将模拟的响应视电阻率作为反演的输入,结合平滑约束最小二乘反演,分析了异常特征和分布规律,其结论和认识具有一定实用价值。  相似文献   

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.
目前,激电数据的二维反演主要针对高密度激电数据,而对于常规的激电测深则很少涉及.据此,针对常规激电测深二维反演所涉及到的一些技术难点进行探讨.在正演过程中,根据实测数据实现有限元网格自动剖分,并采用双网格系统,用相邻节点电导率双线性连续进行有限元模拟,并对模拟的误差进行校正.在反演过程中,利用互换原理和Broyden的方法计算偏导数矩阵;采用与奇异值分解算法精度相当的共轭梯度算法求解偏导数矩阵方程.最后通过对模拟数据和实测数据进行反演试算,验证本计算方法是可行的.  相似文献   

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.
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.
高精度、高效率的深海勘探是当前国家海洋资源勘探开发的研究热点之一,设计观测效率高、勘探成本低的探测方式,开发可精确计算海底复杂电性环境的正演模拟算法有助于推进该研究的进展。本文基于前人的研究经验,引入了全拖曳式深海直流电阻率观测方式,同时考虑到海底沉积环境的电各向异性,开发了基于非结构有限元方法的三维任意各向异性深海直流电阻率正演算法,实现了对海底任意电各向异性情况的仿真模拟。本文算法与层状一维模型半解析解的对比验证了算法的精度。对典型各向异性电性模型的仿真模拟和分析,证明了全拖曳式深海探测方式对海底沉积层和矿产资源电各向异性的高分辨率。在此基础上,总结了电导率各向异性主轴沿笛卡尔坐标旋转时视电阻率分布也发生同方向旋转的视电阻率分布特征,并根据海底地形对各向异性探测的影响进行了简单探讨。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号