共查询到20条相似文献，搜索用时 109 毫秒
1.
可控源音频大地电磁数据的反演方法 总被引：18，自引：7，他引：18
从反演方程、构造目标函数和求解三方面对用于可控源音频大地电磁法(CSAMT)的实用反演方法中的四种进行了描述．水平层状地层CSAMT法资料的直接反演法首次尝试了一维空间的全资料CSAMT反演，效果较好，但该方法尚难应用于2D、3D复杂介质中；奥克姆反演方法既考虑了横向的光滑函数，又考虑了纵向的光滑函数，得到比较光滑的横向、纵向变化的背景电性结果，但有可能把一些小构造光滑掉．快速松驰反演算法和共轭梯度算法由于计算速度快，占内存少而被用于三维反演中，二者相比，快速松驰算法在求解雅可比矩阵时只做一次正演计算，在更新模型时解小型方程组，所以在速度上更胜一筹．在后三种算法中，由于复杂电性结构无解析解，正演计算都采用数值计算．数值计算的可靠性、速度影响着反演算法的有效性，这方面的研究也将是2D、3D复杂电性结构反演的研究方向之一. 相似文献
2.
地壳厚度和波速比是研究地壳结构和组分的两个重要参数,可为区域构造研究提供重要约束。接收函数被广泛地用于确定地壳厚度和波速比,例如H-κ方法或H-κ-c方法,但是该方法只能确定台站下方的地壳厚度和速度比,当地震台站分布稀疏时,很难约束台站间的横向不均匀性。另一方面,重力数据也可用于莫霍面的起伏变化研究,它在横向上覆盖很好,有较高的分辨率,但在纵向上分辨率相对较低。为此,本研究提出了一种联合反演算法求解莫霍面深度和地壳波速比参数。联合反演算法综合考虑了接收函数在纵向上的较高分辨率和重力数据在横向上的较高分辨率,同时拟合区域内所有台站上的接收函数和区域重力数据。模型测试表明联合反演算法较单一的接收函数反演更精确,特别是对于地壳厚度的确定。 相似文献
3.
使用Frankel提出的基于空间光滑地震活动性模型的地震危险性分析方法,选择华南、华北、川滇3个地区的地震记录,比较分析了高斯、幂律和地震分形分布光滑函数3种光滑函数在不同地区的适用性．结果表明,使用交叉验证法可以为高斯光滑函数选取合适的相关距离c值,光滑得到的地震活动性模型能够真实反映研究区域的地震活动特征,根据活动性模型计算得出的峰值加速度(PGA)分布也符合人们对研究区域地震危险性的认识．幂律光滑函数适用于地震活动性较强的地区,且具有容易求取光滑参数的优点．光滑程度较低的幂律光滑函数不适用于地震活动性弱的地区,在该类地区应选择光滑程度较高的高斯光滑函数．地震分形分布光滑函数不适用于地震活动较强且地震活动强度差异较大的地区,其容易过分高估高震级地震对地震危险性的影响,而忽略了低震级地震对地震危险性的贡献．但对于地震活动较弱且地震活动强度差异较小的地区,可使用地震分形分布光滑函数,且同样具有容易求取光滑参数的优点． 相似文献
4.
已有的反演算法研究主要针对地质勘探,对于污染场地的探测数据反演则未有涉及.针对污染场地的特点,建立两种场地模型,通过数值模拟比较两种经典反演算法—阻尼最小二乘法及Occam反演算法.结果表明,两种反演方法的反演结果均能正确反映污染区域的位置,但是对污染区域的范围光滑模型的反演结果优于阻尼最小二乘法的反演结果,并且所需迭代次数小于阻尼最小二乘法.为了进一步比较两种算法在实际污染场地的应用效果,对重庆一污染场地的电阻率法探测数据进行反演.反演结果表明在对实际场地的反演中两种算法可能会存在虚假异常. 相似文献
5.
2.5维复电阻率反演及其应用试验 总被引：1，自引：0，他引：1
自复电阻率法被提出以来,其区分矿与非矿方面的潜力就引起了人们的广泛关注,但反演问题一直没有得到较好的解决,已经严重影响和制约了该方法的应用和发展.对此,本文提出并实现了利用多个排列视电阻率和视相位数据的2.5维SIP联合反演方法.算法利用最小二乘原理构建了反演目标函数,为提高反演的稳定性,在反演方程中加入了Occam法的光滑模型约束.通过借助电场的偏导数形式,推导出了灵敏度矩阵的解析表达式,并应用互换定理对其进行精确求取.反演算法充分利用了电磁感应和激电效应中的异常信息,能够同时反演出二维地质断面上所有单元的四种复电阻率参数.最后,利用该反演程序对安徽某地区的SIP实测数据进行了反演成像,通过与已知钻井资料及CSAMT反演结果的对比分析表明,该反演算法具有良好的应用效果. 相似文献
6.
《应用地球物理》2019,(3)
通常基于不同稳定泛函的反演对目标体边界的刻画存在着差异,尤其对复杂地质构造。为此,本文提出了面向电性边界识别的基于不同稳定泛函的极值边界反演算法。首先选择最平滑、最小支撑两种稳定泛函研究自适应正则化反演算法对不同稳定泛函的适用性,然后通过引入极值边界反演算法,研发基于不同稳定泛函的电性边界识别技术。研究结果表明:自适应正则化反演方法适用于不同的稳定泛函,而且对初始模型的依赖性较低;极值边界反演算法得到的光滑和聚焦上下界并集的比值可以作为划分电性界面的有效依据。在此基础上,本文以内蒙古塔木素地区高放废物粘土岩处置库场址预选为例,将本文提出的反演算法应用于该区的大地电磁测深(MT)数据反演,利用光滑反演了解区域构造,利用聚焦反演了解局部细节,并结合地表地质和钻孔资料,基于地下地质体的电性边界的识别,圈定了目标层的厚度,得到的目标层的厚度满足高放废物处置库选址的要求。 相似文献
7.
8.
井间电阻率层析成象的某些进展 总被引：7，自引：2，他引：7
用有限元法作电阻率层析成象的核心问题是Ｊａｃｏｂｉ矩阵的求取和反演成象算法的研究。本文主要论述了这两方面的最新进展，即求取Ｊａｃｏｂｉ矩阵的格林函数法，用于反演成象计算的平滑度约束反演。文中通过理论模型的成象计算，比较了两种反演算法的成象效果。 相似文献
9.
邱忠坊 《CT理论与应用研究》1994,3(3):27-32
本文是文献［８］的继续，文中提出了地震勘探的带阻尼波动方程反演的新方法，应用此反演算法可以通过微机识别地层的波阻抗，这里把地层分成１００层，地震数据含有１０％的误差，且反演计算时间仅仅是正演计算时间的三倍．文献［１］－［４］的计算方法是基于波动方程的混合问题，而文献［８］所提出的反演算法是基于波动方程的哥西问题，许多例子充分说明后者更有效． 相似文献
10.
对地震层析成像非线性问题线性化处理之后,各种反演算法归纳成为对不适定方 程的求解．地震层析成像反演算法的解的物理意义是给出地质结构,因此对于解的可靠性及 分辨率研究非常重要．然而许多反演算法不能给出解的评价方法,因而对解的可信度产生怀 疑．本研究根据解估计的分辨率矩阵的原理,提出ＬＳＱＲ（Ｌｅａｓｔ Ｓｑｕａｒｅ ＱＲ）算法解协方差矩 阵的评价算法,用相关分析可以为那些在求解过程中得不到分辨率矩阵的反演方法提供解的 定量评价．并用本文提出的解的定量评价方法试评了一个实际地壳模型的地震层析成像的 速度重建结果． 相似文献
11.
12.
高分辨率Radon变换方法及其在地震信号处理中的应用 总被引：32，自引：19，他引：13
Radon变换方法在地震资料处理中广泛采用,在地震同相轴识别和估计方面具有良好效果．无论是倾斜叠加,还是广义Radon变换方法,一般采用最小二乘反演方法实现．目前,在提高反演算法的效率和分辨率方面仍值得研究．本文从倾斜叠加的定义出发,阐明Radon变换分辨率问题的来源和解决办法．采用最小二乘反演方法研究高分辨率抛物线Radon变换和双曲Radon变换时,给出稀疏约束预条件共轭梯度法求解的高分辨率Radon变换的实现方法,同阻尼最小二乘方法相比,分辨率和精度明显提高,文中给出了模型算例．根据有效波和多次波NMO后剩余时差不同,采用高分辨率抛物线和双曲Radon变换可以压制多次波,分别给出了方法原理,最后给出应用实例．研究表明,稀疏约束预条件共轭梯度法可以有效实现高分辨率Radon变换;数值算例表明,算法计算效率和精度较高,可以更好地实现多次波压制． 相似文献
13.
CT机和核磁共振成像(MRI)机之间的比较研究 总被引：5，自引：1，他引：4
ＣＴ和ＭＲＩ是当今应用最普遍的两大类临床影像诊断设备。本文应邀就普遍关心的问题对ＣＴ和ＭＲＩ两大设备在原理、性能、指标、功能、参数、敏感性、安全性、信号性质、图像标准、图像重建的方法学理论等诸多方面进行了有益和有意义的比较。并沿两类机器发展的轨迹追溯了两种影像在历史上的交汇点（源宗）。阐明了两套数学方法（Ｒａｄｏｎ变换和Ｆｏｕｒｉｅｒ变换）的平行发展根源于其信号性质的重大区别。ＭＲ是时域信号,而Ｘ 相似文献
14.
Radon transform is a powerful tool with many applications in different stages of seismic data processing, because of its capability to focus seismic events in the transform domain. Three-parameter Radon transform can optimally focus and separate different seismic events, if its basis functions accurately match the events. In anisotropic media, the conventional hyperbolic or shifted hyperbolic basis functions lose their accuracy and cannot preserve data fidelity, especially at large offsets. To address this issue, we propose an accurate traveltime approximation for transversely isotropic media with vertical symmetry axis, and derive two versions of Radon basis functions, time-variant and time-invariant. A time-variant basis function can be used in time domain Radon transform algorithms while a time-invariant version can be used in, generally more efficient, frequency domain algorithms. Comparing the time-variant and time-invariant Radon transform by the proposed basis functions, the time-invariant version can better focus different seismic events; it is also more accurate, especially in presence of vertical heterogeneity. However, the proposed time-invariant basis functions are suitable for a specific type of layered anisotropic media, known as factorized media. We test the proposed methods and illustrate successful applications of them for trace interpolation and coherent noise attenuation. 相似文献
15.
一类广义Radon变换的反演:迭代重建再投影的外插算法 总被引：1，自引：0，他引：1
本文对不完全数据下的一类广义Radon变换的反演问题,建议用迭代重建再投影的外插算法重建图像。数值模拟结果表明,迭代重建再投影的外插算法对重建图像是有效的。 相似文献
16.
基于部分动校正（NMO）后反射同相轴在CMP道集上的抛物线走时近似，给出了加权抛物Radon变换叠前地震数据重建方法（WPRT）. WPRT通过在迭代过程中引入变化着的权系数，拓展和改进了传统抛物Radon变换方法，使其可同时完成不规则采样的规则化和空道及近偏移距道重建，且有更高的计算效率. 文中给出了应用WPRT进行近偏移距和中偏移距的空地震道重建及数据规则化的算法实现. 理论模型和实际地震资料的地震数据重建结果显示了本文算法的优点. 相似文献
17.
18.
高分辨率Radon变换存在计算效率和分辨率不能兼得的困境.时间域算法可以获得很高的分辨率,但计算效率非常低;频率域算法具有良好计算效率,但分辨率不理想.为此发展了混合域高分辨率抛物Radon变换,即对频率域抛物Radon变换引入时变的稀疏权.本文给出了一种新的混合域高分辨率抛物Radon变换实现方法,并将该算法应用于叠前数据衰减多次波.文中给出了Radon变换和衰减多次波的流程.理论和实际数据算例表明本文方法既有较高的分辨率又有很高的计算效率. 相似文献
19.
The parabolic Radon transform has been widely used in multiple attenuation. To further improve the accuracy and efficiency of the Radon transform, we developed the 2- fdomain high-resolution Radon transform based on the fast and modified parabolic Radon transform presented by Abbad. The introduction of a new variable 2 makes the transform operator frequency-independent. Thus, we need to calculate the transform operator and its inverse operator only once, which greatly improves the computational efficiency. Besides, because the primaries and multiples are distributed on straight lines with different slopes in the 2-fdomain, we can easily choose the filtering operator to suppress the multiples. At the same time, the proposed method offers the advantage of high-resolution Radon transform, which can greatly improve the precision of attenuating the multiples. Numerical experiments suggest that the multiples are well suppressed and the amplitude versus offset characteristics of the primaries are well maintained. Real data processing results further verify the effectiveness and feasibility of the method. 相似文献
20.
Radon变换在实际应用中必须考虑抽样取值的问题,通信理论中常用的Shannon取样定理在地球物理学等领域中并不适用,因为它要求函数是频谱有限的.本义介绍利用小波理论建立新的取样定理的方法,这一方法的优点在于可以根据具体问题的要求选取适当的基本空间V_0,从而获得较佳的抽样插值效果. 相似文献