共查询到19条相似文献，搜索用时 93 毫秒
1.
2.
3.
为使磁测数据的解释结果能够刻画复杂的磁性体,探讨了应用三维磁化率成像技术来反演磁测数据.本文应用基于Occam反演的模型粗糙度约束的光滑成像算法来求解纯欠定的磁化率成像问题.在求解大规模线性方程组时,本文使用预优共轭梯度法并将深度加权函数作为预优因子,既提高了反演的收敛速度,又使磁化率分布在合理的深度上.采用绝对约束的方法对反演的模型值的取值范围进行了约束,取得了较好的约束效果.二维和三维典型模型的数值实验表明,本文给出的磁化率成像方法可以恢复真实模型的主要特征信息.将本文的磁化率成像反演技术应用于大冶铁矿狮子山矿段的实测数据的反演解释,获得了与地质信息具有很好对应关系的成像结果. 相似文献
4.
井中激电是二次找矿重要的地球物理勘探手段,快速而稳定的正反演算法有助于方法的推广和应用.本文在正演模拟中,给出了考虑井眼影响下的网格剖分方式;用右端项校正技术减小边界效应和源点奇异性引起的模拟误差;并采用循环Krylov子空间算法提高多线性方程组的求解效率.反演用Gauss-Newton法结合Jacobian-free Krylov迭代求解技术,给出了Jacobian矩阵向量积的简化计算方法;用不精确预处理共轭梯度法对模型修正量方程近似求解以减少计算量;采用不同于正演的反演网格剖分降低不适定性.数值算例验证了相关算法的有效性和可靠性. 相似文献
5.
正演计算是反演研究的基础,为了实现基于三维弹性波方程的全波形反演成像,发展准确、高效、低数值频散的三维正演模拟方法至关重要.为此,本文将修正保辛分部龙格-库塔格式与优化有限差分算子结合,发展了用于数值求解三维弹性波方程的修正时空优化保辛方法(MTSOS).新方法使用二级龙格-库塔格式达到了三阶时间精度,且更适用于求解非均匀介质情况下的弹性波方程,数值频散误差小于同精度保辛分部龙格-库塔(SPRK)方法的误差,提高了计算精度.波场模拟结果表明,三维MTSOS方法可以精确给出数值模拟结果,能够清晰模拟地震波传播过程中产生的各种震相、有效压制数值频散. 相似文献
6.
测井数据和地震数据是地震勘探中两种最重要的资料. 测井约束地震波形反演是在非线性波形反演的基础上,利用已知测井资料详细的垂直分辨能力和地震资料均匀密集的水平采样特点, 通过迭代反演来求取一个具有较高分辨率的速度参数.本文建立了测井约束反演模型,研究了测井约束下地震波形反演的同伦摄动求解方法.同伦摄动法作为一种新的、求解数学物理中各种非线性问题的有效方法,具有计算速度快、计算精度高的优点.这对于提高反演的精度和效率是十分有益的. 为了表征该方法的有效性和稳定性,分别对水平层状介质模型和逆冲断层带模型进行了数值模拟,并与Landweber迭代法相对比,结果表明该算法具有更好的收敛性,能够取得更为满意的反演效果. 相似文献
7.
《应用地球物理》2015,(3)
频率域全波形反演充分利用全波场的振幅、相位以及频率信息,采用较少的频率便能反演得到精度很高的速度模型。本文以有限单元法为基础,对起伏地形条件下二维声波频率域全波形反演进行了研究。在正演算法中,针对截断边界问题,并考虑多频率联合反演中计算区域采用同一套剖分网格的需求,提出了一种适用于起伏地形的衰减边界条件算法。该算法的核心思想是在控制方程波数项中引入衰减因子,通过一定方式调节衰减因子使得声波在衰减层中充分衰减,达到压制截断边界影响的目的。根据指数衰减规律,文中推导出了一种新的衰减因子计算公式,并给出了不同频率条件下衰减层厚度计算公式;在反演算法中,采用共轭梯度法求解高斯牛顿反演迭代方程组,避免直接求解雅克比矩阵和HeSSian矩阵带来的巨额计算量,并采用相同的反演模型,对比分析了不同初始模型和频率组合对全波形反演结果的影响。起伏地形模型数值模拟和全波形反演数值试验表明,本文提出的指数衰减边界条件算法和基于该算法的全波形反演算法具有很好的应用效果。 相似文献
8.
同时反演地层电阻率与侵入深度的数值方法 总被引：4，自引：2，他引：4
作者曾从电法测井的数学模型出发,在假定侵入深度为已知的条件下,给出了反演地层电阻率的一个迭代方法。本文在文献[1]的基础上,将侵入深度也作为反演参数之一,通过适当的变量代换,运用微分方程反问题的数值求解方法,得到了一个可以同时反演地层电阻率与侵入深度的迭代方法。 相似文献
9.
利用天然地震记录直达P波后续20s内的波形信息,研究了地壳、土地慢速度结构反演的方法、唯一性、精度及应注意的问题,并将求解有条件极值的惩罚数法和求解无条件极值的单纯形法引入到波形反演.数值计算证明,波形反演方法对记录的误差有压制作用,在深源、远震、各向同性水平层状近似合理的前提下,可用一个地震事件在射线平面内的二分量记录反演台站区地壳、上地幔P、S波速度结构.在新疆天山地区转换波流动台站观测中,选择出较好的记录进行波形反演,给出了该转换被测深剖面上6个台站下方的地壳、上地慢P、S波速度结构.通过塔里木盆地内一测点上波形反演和地震勘探结果的对比证明,波形反演方法具有满意的精度. 相似文献
10.
基于交叉梯度结构约束的不同数据的联合反演可以提高地球物理成像的可靠度,但是由于不同观测数据对地下介质的灵敏度不同以及不同反演系统灵敏度矩阵元素的数值可能存在较大的差异,导致构建和求解联合反演系统存在很大的挑战.针对以上问题,本文提出一种新的基于单独反演模型更新量与交叉梯度结构约束相结合的联合反演策略.该策略利用单独反演系统分别确定出各个模型的更新量,然后利用它们约束交叉梯度系统的反演,得出新的模型更新量.通过这样的方式,有效实现了数据拟合与结构约束的平衡,实现了基于交叉梯度联合反演的目标.新的联合反演策略不需要对原来单独反演程序修改即可实现联合反演,减小了联合反演实现的难度,极大地提高了联合反演的易实现性,而且避免了联合反演矩阵存储及结构过于复杂难以求解的问题.基于新的策略,本文首次实现了基于交叉梯度结构约束的三维地震走时与直流电阻率联合反演.合成模型测试表明,与单一成像相比,联合成像减少了地震走时反演中出现的干扰异常并提高了电阻率反演的分辨率. 相似文献
11.
12.
本文以常用的几种粘弹性体模拟地球介质的粘弹性,建立了该种介质中的波动方程组;在运用合理的近似处理方法基础上,得到了基本解;提出了校正其它衰减因素的功率谱切比雪夫拟合方法,由此给出了粘性Q值的概念和算法,最后讨论了该方法的可行性。 相似文献
13.
14.
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. 相似文献
15.
In mathematical statistical filtering the deconvolution problem can be solved by two different methods:
- 1 by inverse filtering
- 2 by calculating the prediction error.
16.
17.
根据积分变换方法得出了半空间内部作用简谐水平力时的Gree函数．首先，利用Hankel积分变换方法，直接对频域内的Biot波动方程进行求解，得出Biot波动方程的通解；利用通解和半空间内部作用水平力时边界上的边界条件，以及力作用面上的连续性条件，可以得出上述边值问题的解；对于边值问题在变换域内的解进行相应的逆变换，就可以得出频域内的Gree函数．本文得到的线弹性退化解与文献中的结果吻合．最后，文中给出了两个算例． 相似文献
18.