共查询到18条相似文献，搜索用时 188 毫秒
1.
2.
用FPS定量解释法实现三维近似成像的探讨 总被引：1，自引：0，他引：1
目的 在激发极化数据解释中，过去常采用原始曲线的特征值法，现介绍一种用FPS定量解释法确定异常体空间位置的三维近似成像方法。方法 用FPS定量解释法确定异常体空间位置的三维近似成像方法。结果 该方法可以较准确的确定异常体的上下左右边界以及异常体在垂直测线方向上的延伸情况。结论 算例表明，这种作图解释成像的方法效果较好，值得推广应用。 相似文献
3.
4.
5.
重力勘探中复杂条件下的三维正演计算量大存储要求高,使得这种条件下重力勘探高效、精细正反演变得困难.针对这一问题,提出一种空间-波数混合域数值模拟方法,该方法将空间域引力位积分进行水平方向二维傅里叶变换,将三维空间域卷积问题转换为多个不同波数之间相互独立的空间垂向一维积分问题,一维积分垂向可离散为多个单元积分之和,每个单元采用二次形函数表征密度变化,可得出单元积分的解析表达式.该方法计算量和存储需求少,算法高度并行;保留垂向为空间域,优势之一在于可根据实际情况合理调整单元疏密程度,准确模拟任意复杂地形和密度异常体的重力异常,兼顾计算精度与计算效率;优势之二在于用形函数拟合求得积分的解析解,计算精度和效率高;充分利用一维形函数积分的高效和高精度,不同波数之间一维积分高度并行性及快速傅里叶变换的高效性,实现重力异常场三维数值模拟.设计棱柱体模型,通过数值解和解析解对比验证了该方法的正确性、适用性和高效性.针对任意复杂地形条件下的重力场及其张量的模拟问题,提出一种快速算法,对其有效性进行了验证.探究标准FFT法的截断效应对计算精度的影响,对比分析Gauss-FFT法和标准FFT扩边法两种方法的计算精度和效率,总结了二者的选取策略,结果表明选用标准FFT扩边法计算效率更高.实际地形的数值模拟表明本文算法适用于任意复杂地形的高效计算. 相似文献
6.
7.
采用刚度矩阵方法结合Hankel积分变换,求解了层状黏弹性半空间中球面SH、P和SV波的自由波场.首先,在柱坐标系下建立层状黏弹性半空间的反轴对称(柱面SH波)和轴对称(柱面P-SV波)情况精确动力刚度矩阵.进而由Hankel变换将空间域内的球面波展开为波数域内柱面波的叠加,然后将球面波源所在层的上下端面固定,求得固定层内的动力响应和固定端面反力,将固端反力反向施加到层状黏弹性半空间上,采用直接刚度法求得固端反力的动力响应,叠加固定层内和固端反力动力响应,求得波数域内球面波源动力响应.最后由Hankel积分逆变换求得频率-空间域内球面波源自由场,时域结果由傅里叶逆变换求得.文中验证了方法的正确性,并以均匀半空间和基岩上单一土层中球面SH、P和SV波为例分别在频域和时域内进行了数值计算分析.研究表明基岩上单一土层中球面波自由场与均匀半空间情况有着本质差异;基岩上单一土层中球面波位移频谱峰值频率与场地固有频率相对应,基岩面的存在使得基岩上单一土层地表点的位移时程非常复杂,振动持续时间明显增长;阻尼的增大显著降低了动力响应的峰值,同时也显著减少了波在土层的往复次数. 相似文献
8.
9.
拟解析近似方法是一种解决电磁场散射问题的快速求解积分方程的近似方法,它绕开了传统数值方法中的求解大型代数方程组或大型矩阵问题,适用于强散射和大扰动问题.本文应用孙建国提出的求解异常电场积分方程的张量拟解析近似理论公式,研究用其求解直流电场积分方程.利用接近实际的地电模型对异常电场进行模拟,研究了均匀场中异常球体的张量拟解析近似解;并对均匀场中的立方体异常体进行了数值计算.效果良好并具有很高的计算速度.研究结果为三维直流电场快速正反演模拟打下基础. 相似文献
10.
本文研究了有限单元法非均匀剖分对直流电阻率法数值模拟结果的影响.通过非均匀剖分对半空间地质体的视电阻率异常与均匀剖分的视电阻率对比,分析这种模拟方法的可行性,研究计算结果差异产生的原因.通过相邻单元的不同边长比例剖分的对比,研究在误差允许范围内非均匀剖分的单元边长比例极限.计算结果表明,测线方向非均匀剖分是误差产生的主要原因,XYZ三方向非均匀剖分所带来的误差总和是各方向非均匀剖分产生误差之和;当测线方向剖分单元不均匀剖分比例小于0.49∶0.01时,由于非均匀剖分产生的误差在容许范围内. 相似文献
11.
E.H. ELORANTA 《Geophysical Prospecting》1986,34(6):856-872
An integral equation method is described for solving the potential problem of a stationary electric current in a medium that is linear, isotropic and piecewise homogeneous in terms of electrical conductivity. The integral equations are Fredholm's equations of the ‘second kind’ developed for the potential of the electric field. In this method the discontinuity-surfaces of electrical conductivity are divided into ‘sub-areas’ that are so small that the value of their potential can be regarded as constant. The equations are applied to 3-D galvanic modeling. In the numerical examples the convergence is examined. The results are also compared with solutions derived with other integral equations. Examples are given of anomalies of apparent resistivity and mise-a-la-masse methods, assuming finite conductivity contrast. We show that the numerical solutions converge more rapidly than compared to solutions published earlier for the electric field. This results from the fact that the potential (as a function of the location coordinate) behaves more regularly than the electric field. The equations are applicable to all cases where conductivity contrast is finite. 相似文献
12.
Tarkeshwar Lal 《Pure and Applied Geophysics》1970,82(1):259-269
Summary The theory of resistivity sounding on a three layer earth comprising an inhomogeneous interstratum has been presented. The electrical conductivity in the intermediate layer, embedded between two layers of uniform conductivity, is assumed to follow either generalised power law or exponential variation. A point source of current is placed on the surface. The electric potential functions on the ground surface have been obtained. The expressions for apparent resistivity over the above models for the Wenner and Schlumberger electrode configurations have been derived. 相似文献
13.
首先给出柱坐标系中电导率分段线性变化的水平层的点源电场的二维边值问题,然后用变分法将边值问题转变为变分问题。用有限单元法解变分问题,将区域剖分成矩形单元,在单元中进行双线性函数插值,将变分方程化为线性代数方程组。解方程组,得各节点的电位值,由此可计算地表的视电阻率。 算例表明,本方法计算结果与精确解十分符合。本文还举了一个定量分析视电阻率年变化的例子。 本方法占用计算机内存约100K数量级。在MV/6000超小型计算机上计算一条电测深曲线的时间为几十秒钟。 相似文献
14.
大地电磁法是广泛应用于深部地质结构探测、油气和矿产资源勘查等领域的一种地球物理方法.电性各向异性对电磁观测数据有很大影响,但介质各向异性参数对不同模式视电阻率的影响还较少有较为系统的研究.本文基于Maxwell方程组,推导了二维大地电磁场在任意各向异性介质中电场和磁场相互耦合的变分方程,结合有限单元法及并行计算编写了二维大地电磁任意各向异性正演程序,采用三角形网格剖分.验证程序正确性后,以倾斜板状体作为模型来研究三个主轴电阻率及三个旋转欧拉角和四种模式的视电阻率之间的关系.结果表明,主轴各向异性时,xy模式视电阻率几乎只受x方向电阻率影响,yx模式视电阻率主要受y方向电阻率影响,但同时也受z方向电阻率一定影响;三个欧拉角中只有倾角不为零时,yx模式视电阻率受倾角大小的影响较大,xy模式视电阻率几乎不受倾角的影响;只有走向角不为零时,四种模式的视电阻率同时受x、y两个主轴电阻率和走向角的大小的影响. 相似文献
15.
16.
针对任意各向异性地层,利用极向型和环向型标量位函数,导出相应的直流电视电阻率和磁电阻率的磁场响应关系.计算了各向异性地层的直流电视电阻率和磁电阻率响应特征,重点分析了电阻率测深方法与磁电阻率探测方法对地下各向异性介质的探测能力.文中采用状态矩阵的分析方法,首先采用极向型和环向型标量位构造了各向异性层状介质电场与磁场的通解,利用各层界面电场、磁场的连续性及地面激励源的耦合条件,推导了不同层之间电磁场的状态矩阵,建立了空间电场与磁场的递归计算关系.其次,针对递归计算中指数项数值计算的不稳定性,借用状态矩阵的性质,导出了将不稳定指数计算项转化为稳定的指数项的转换关系.针对横向各向同性(TI)介质中极向位与环向位解耦的特点,导出了电磁场的直接积分解.最后,采用解析解验证了算法的正确性,给出了多层各向异性地层模型的视电阻率和磁场响应曲线,分析了直流电法探测裂缝性地层、估计裂缝分布性状的可能性. 相似文献
17.
基于非结构网格有限元方法开展了三维复杂地电模型的线源井地电法的高效正演模拟研究,探讨了通过求取电场响应导数来刻画目标体边界范围、采用差异场地形校正技术来消除地形影响等措施对井地电法成像的效果和精度的影响。并通过对比与解析解,验证了本文数值解算法的有效性。模型计算结果表明：积水巷道的空间位置和走向均会引起视电阻率的显著变化,视电阻率变化率的极值准确且清晰地指示了巷道边界的位置;电位的归一化总水平导数极大地提高了井地电法对目标体复杂边界位置的识别能力;地形对井地电场分布的影响也很大,其视电阻率响应与地形形状近似呈对称关系,利用差异场技术能有效地削弱地形对井地电法高精度成像的影响。 相似文献
18.
二维层状介质点电源边界单元解 总被引：1，自引：0，他引：1
本文用付氏变换将点源二维地电断面的三维边值问题变成二维边值问题,然后讨论了边界元法在解决层状介质中存在二维柱体时点源场视电阻率计算问题中的应用。由于选取的基本解避免了对水平层状界面的积分计算,从而大大减少了边界剖分节点数。在计算中设法略去了对边界上位场法向导数的求解,使未知量大为减少,从而便利了求解二维边值问题,并由付氏反变换求得三维空间中的位。本方法所需计算机内存少,原始数据的输入简单,且有较高的计算精度。 相似文献