共查询到20条相似文献，搜索用时 218 毫秒
1.
2.
本文前一部分叙述了点源二维电阻率法数值模拟正演计算的有限单元法,文中采用混合边界条件、用LL~T分解解线性方程组等优化措施,使二维有限单元法电阻率法正演计算的速度和精度都比目前国外流行的L.Rijo方法和程序有了进一步的提高。文中介绍了有关的方法并引入了若干计算结果。本文后一部分叙述了三维电阻率法数值模拟的积分方程近似解法,由于对K.Dieter等人提出的方法作了一些近似处理,并用迭代法求解积分方程表示式,实算结果说明在保证计算精度的条件下,提高了计算速度,从而提高了三维电场正演计算的有效性和实用性。所提出的方法容易推广到激发极化法的正演计算中。 相似文献
3.
4.
复电阻率法二维有限元数值模拟 总被引：9，自引：2，他引：9
伴随着复电阻率法的广泛应用，发展精确和快速的正演和反演算法成为复电阻率法研究的重点.本文采用基于三角单元剖分的有限单元法进行了复电阻率二维数值模拟研究.为了提高计算速度，对无穷远边界进行了近似处理.整个正演计算过程分为两步，首先采用有限单元法计算四个不同频率的视复电阻率数据，然后对前一步得到的视复电阻率数据采用递推算法计算视Cole-Cole参数.采用这种正演算法与一维正演的结果进行了对比，验证了本文方法的正确性.设计了两个二维极化模型，数值模拟结果表明视复电阻率和Cole-Cole视参数等值线断面图对于异常目标体都有比较明显的反映. 相似文献
5.
在阮百尧、徐世浙(1998,2001)工作基础上,给出四边形内三角网格单元剖分,电导率连续变化2.5维直流电阻率数值模拟方法.在二维起伏地形情况下,有限单元设计为四边形内剖分四个三角网格单元,单元内的电位值和电导率设计为线性变化,将研究区域远边界的第三类边界条件简化为第二类齐次边界条件,重新给出一组适用性较好的Fourier反变换系数,并采用分离式Cholesky分解法快速求解线性方程组.通过对3个地电模型的计算与对比分析,其结果与解析法的均方根误差小于0.2%,计算速度比较常规方法提高了5倍,能有效模拟起伏地形和复杂形态的地电体模型. 相似文献
6.
对于时间域航空电磁法二维和三维反演来说,最大的困难在于有效的算法和大的计算量需求.本文利用非线性共轭梯度法实现了时间域航空电磁法2.5维反演方法,着重解决了迭代反演过程中灵敏度矩阵计算、最佳迭代步长计算、初始模型选取等问题.在正演计算中,我们采用有限元法求解拉式傅氏域中的电磁场偏微分方程,再通过逆拉氏和逆傅氏变换高精度数值算法得到时间域电磁响应.在灵敏度矩阵计算中,采用了基于拉式傅氏双变换的伴随方程法,时间消耗只需计算两次正演,从而节约了大量计算时间.对于最佳步长计算,二次插值向后追踪法能够保证反演迭代的稳定性.设计两个理论模型,检验反演算法的有效性,并讨论了选择不同初始模型对反演结果的影响.模型算例表明:非线性共轭梯度方法应用于时间域航空电磁2.5维反演中稳定可靠,反演结果能够有效地反映地下真实电性结构.当选择的初始模型电阻率值与真实背景电阻率值接近时,能得到较好的反演结果,当初始模型电阻率远大于或远小于真实背景电阻率值时反演效果就会变差. 相似文献
7.
8.
9.
研讨了频率域电磁法中不同源装置的大地电磁测深、线源频率电磁测深和偶极源频率电磁测深阻抗视电阻率的源效应影响特征。在唯象分析的基础上，提出了几种电磁测深法阻抗视电阻率的相互换算法──源效应校正法（大地电磁测深二维ＴＥ极化视电阻率和其它两种电磁法的赤道装置二维阻抗视电阻车）。模型试验表明，利用这一源效应校正法可以由大地电磁二维视电阻率近似地计算出线源频率电磁二维阻抗视电阻率。这一方法被尝试应用于由线源频率电磁二维阻抗视电阻率估算偶极源频率电磁二维阻抗视电阻率。 相似文献
10.
电阻率法有限单元正演模拟中,采用第三类边界条件时为保证精度仍要求较大范围的计算域.无单元法为地球物理领域的新兴正演模拟方法,其计算效率低,但其中采用的移动最小二乘(MLS)形函数相比于有限单元法形函数具有良好的连续性,模拟精度高.本文将MLS形函数应用于电阻率法有限单元2.5维正演的第三类边界条件处理,提出电阻率法有限单元-移动最小二乘(FEM-MLS)耦合正演方法.通过不同正演方法的模型算例模拟结果对比,验证了本文算法的有效性,并讨论了各个参数选择对模拟结果的影响.本文数值模拟结果表明采用第三类边界条件时,在同等计算精度前提下,FEM-MLS耦合法相比于有限单元法可进一步缩小计算域并提高了计算效率,相比于采用较大计算域满足边界条件的有限单元法计算效率提高了约一倍,相比于采用相同小范围计算域的有限单元法平均精度提高了约一倍. 相似文献
11.
地表起伏地形在野外矿产资源勘察中不可避免,其对直流电阻率法勘探影响巨大.近年来,电阻率三维正演取得诸多进展,特别是应用非结构网格我们能够进行任意复杂地形和几何模型的电阻率三维数值模拟,但面向实际应用的起伏地形下电阻率三维反演依然困难.本文基于非结构化四面体网格,并考虑到应用GPS/GNSS时,区域地球物理调查中可非规则布设测网的实际特点,实现了任意地形(平坦或起伏)条件下、任意布设的偶极-偶极视电阻率数据的不完全Gauss-Newton三维反演.合成数据的反演结果表明了方法的有效性,可应用于复杂野外环境下的三维电法勘探. 相似文献
12.
2.5-D直流电阻率有限元数值模拟中,模型的剖分及加密主要通过手动实现.另外,采用的单元类型比较规则如矩形单元等,不易实现复杂模型的模拟.为解决上述问题,文中提出了一种自适应有限元算法.算法中采用稳健的后验误差估计来自动预测下一次网格的单元尺寸,直到设定的迭代条件满足为止.另外,采用非结构化三角形单元实现了任意复杂模型的灵活剖分.基于此,利用垂直接触面模型分析和对比了不同自适应策略的效率.通过对比发现,点源附近的单元得到了加密以消除源的奇异性.另外,对于任意一种策略,有限元结果均能最终收敛到精确解.最后,模拟了两个模型：2-D单个异常体模型和2-D地形模型. 相似文献
13.
To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner (AGMG-CG) to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling (3DDCXH). In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity. 相似文献
14.
Computations of secondary potential for 3D DC resistivity modelling using an incomplete Choleski conjugate-gradient method 总被引：4，自引：0，他引：4
An accurate and efficient 3D finite-difference (FD) forward algorithm for DC resistivity modelling is developed. In general, the most time-consuming part of FD calculation is to solve large sets of linear equations: Ax = b , where A is a large sparse band symmetric matrix. The direct method using complete Choleski decomposition is quite slow and requires much more computer storage. We have introduced a row-indexed sparse storage mode to store the coefficient matrix A and an incomplete Choleski conjugate-gradient (ICCG) method to solve the large linear systems. By taking advantage of the matrix symmetry and sparsity, the ICCG method converges much more quickly and requires much less computer storage. It takes approximately 15 s on a 533 MHz Pentium computer for a grid with 46 020 nodes, which is approximately 700 times faster than the direct method and 2.5 times faster than the symmetric successive over-relaxation (SSOR) conjugate-gradient method. Compared with 3D finite-element resistivity modelling with the improved ICCG solver, our algorithm is more efficient in terms of number of iterations and computer time. In addition, we solve for the secondary potential in 3D DC resistivity modelling by a simple manipulation of the FD equations. Two numerical examples of a two-layered model and a vertical contact show that the method can achieve much higher accuracy than solving for the total potential directly with the same grid nodes. In addition, a 3D cubic body is simulated, for which the dipole–dipole apparent resistivities agree well with the results obtained with the finite-element and integral-equation methods. In conclusion, the combination of several techniques provides a rapid and accurate 3D FD forward modelling method which is fundamental to 3D resistivity inversion. 相似文献
15.
直流电阻率法通过观测单极-双极装置视电阻率曲线最小值位置估计掌子面前方低阻异常体距离,在巷道超前探测及矿井水害预测与防治工作中有广泛应用,并形成多种经验预测模型.本文基于直流电法巷道超前探测的各向异性电阻率三维有限元数值模拟,给出超前探线性预测模型.然而,实际井下地质结构复杂,异常大小不定且任意分布,而且可能存在有电阻率各向异性,使得目前由简单模型实验或数值模拟获得的预测公式不确定性较大,而且可靠性难于评价.Monte Carlo方法用随机化思想解决非确定性问题,我们将该方法与并行算法结合,对大量不确定的随机模型进行三维数值模拟,以验证各预测模型的准确度及可靠性.结果表明,本文的超前探线性预测模型较其他预测模型更为准确及可靠.但介质各向异性的影响很大,尤其当围岩电阻率具有各向异性时,所有预测模型的准确度及可靠性较差,给安全生产带来很大隐患. 相似文献
16.
射频大地电磁法(RMT)是以潜艇天线发射的射线源等作为场源的一种地球物理勘探方法,目前被广泛应用于近地表工程和环境地球物理勘探.RMT数据解释常采用基于静态假设的大地电磁法(MT)程序,往往会反演出不真实的浅层目标体.为解决这一长期困扰RMT资料解释的问题,本文实现了考虑位移电流的RMT有限元数值模拟.为了处理任意复杂模型,如起伏地形,非结构的三角形网格被用于离散RMT模型.首先通过算例对比,验证了程序的正确性和可靠性.通过Dike模型讨论了空气层厚度对RMT数值解的影响,结果表明当空气层厚度大于1/4波长即可满足精度要求.以山脊模型为例计算了位移电流对RMT响应的影响,表明位移电流的影响会随着山脊高程的增加而增大.最后通过舒家店实际模型进一步验证了位移电流在RMT中的重要性. 相似文献
17.
Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell’s equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic (EM) fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling. 相似文献
18.
复电阻率法在矿产、油气勘探调查中发挥着重要作用.为了认识复杂构造的复电阻率法电磁场的变化规律,本文基于自适应有限元方法,采用非结构化网格,引入Cole-Cole模型,实现了电偶源2.5D复电阻率法电磁场正演,可以模拟复杂地形和地电结构,正演结果更符合野外实际地质情况.通过将本文的计算结果与半空间模型解析解、层状介质和起伏模型结果进行对比,验证了本文算法的正确性.最后,基于复杂地电模型,通过正演模拟,系统分析了地形、激电参数、复杂构造对复电阻率法电磁场的影响特征. 相似文献
19.
为解决传统有限元截断边界所引起的问题,本文提出了一种新的三维直流电阻率有限元-无限元耦合数值模拟方法.首先推导了无限元三维单元映射函数,然后提出了一种全新的最优的无限元形函数并与多种其他形函数进行了对比,随后将其与非结构化四面体有限元相结合,取代了传统的混合边界条件,使得电位在无限域内连续并在无限远处衰减为零,最终形成的左端矩阵稀疏对称并与场源位置无关.数值计算表明,该方法可以在近似测区大小的计算范围内得到与混合边界条件相当的计算精度,优于相同计算范围下齐次边界条件的解,有利于减少计算节点数;由于左端矩阵不随场源位置改变,有利于加速反演计算. 相似文献
20.
在海洋可控源电磁法勘探中,接收站常置于海底.在进行海洋电磁场模拟时,由于海水和海底介质存在显著电性差异,这给海底接收点处场值的求取带来困难.本文提出一种新的接收点插值算法,该算法考虑到海底电场法向分量不连续性问题,用法向电流分量进行插值以准确求取海底任意接收点处电磁场值.本文利用交错网格有限差分法实现了二维介质中频率域海洋可控源法(CSEM)正演.对构造走向做傅里叶变换,将三维电磁模拟问题转换为波数域2.5维问题,即三维场源激励下针对二维地电模型的电磁模拟问题.使用交错网格有限差分法,基于一次场/二次场分离方法导出波数域二次电场离散形式,并进一步求得波数域电磁场.采用本文提出的改进的插值算法可求得海底任意接收点处波数域电磁场,采用傅里叶逆变换对波数域电磁场进行积分可得到接收点处空间域电磁场.模型算例表明,与常规的线性插值和严格插值算法相比,本文提出的改进的插值算法具有更高的精度. 相似文献