首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 125 毫秒
1.
本文把波动方程反问题与广义Radon 变换的反演相联系。在假定弱散射条件下,把波动方程反问题转化成广义Radon 变换的反演问题,即如何从一系列关于目标函数在某类子流形上的积分值,去求出目标函数.这种转化提供了一种研究波动方程反问题的途径。  相似文献   

2.
一维波动方程的奇性反演与小波   总被引:4,自引:4,他引:4       下载免费PDF全文
利用线性化的技巧及G.Beylkin引进的奇性反演的概念,使波动方程非均匀背景场的反演问题有了实质性的进展,但在进行线性化简时,往往会将数值小但奇性高的项略去,因而使反演结果失真,本文利用小波变换这一工具,在化简时保留了奇性的主要部分,从而使反演所得的结果从奇性分析的观点看来更为。  相似文献   

3.
高分辨率Radon变换方法及其在地震信号处理中的应用   总被引:32,自引:19,他引:13  
Radon变换方法在地震资料处理中广泛采用,在地震同相轴识别和估计方面具有良好效果.无论是倾斜叠加,还是广义Radon变换方法,一般采用最小二乘反演方法实现.目前,在提高反演算法的效率和分辨率方面仍值得研究.本文从倾斜叠加的定义出发,阐明Radon变换分辨率问题的来源和解决办法.采用最小二乘反演方法研究高分辨率抛物线Radon变换和双曲Radon变换时,给出稀疏约束预条件共轭梯度法求解的高分辨率Radon变换的实现方法,同阻尼最小二乘方法相比,分辨率和精度明显提高,文中给出了模型算例.根据有效波和多次波NMO后剩余时差不同,采用高分辨率抛物线和双曲Radon变换可以压制多次波,分别给出了方法原理,最后给出应用实例.研究表明,稀疏约束预条件共轭梯度法可以有效实现高分辨率Radon变换;数值算例表明,算法计算效率和精度较高,可以更好地实现多次波压制.  相似文献   

4.
由于CT技术的出现,Radon变换受到人们的重视,已被应用于许多领域,但是,还有一些问题(如地球物理中的层析成像问题),已知的数据往往是未知函数在曲线上的积分值,而不是古典Radon变换中的直线。这就使得人们不得不考虑所谓广义的Radon变换:即已知点函数在曲线上的积分值,要求恢复此函数。使我们感到欣慰的是,这方面的工作已取得了不少的进展。  相似文献   

5.
一类广义Radon变换的反演:迭代重建再投影的外插算法   总被引:1,自引:0,他引:1  
本文对不完全数据下的一类广义Radon变换的反演问题,建议用迭代重建再投影的外插算法重建图像。数值模拟结果表明,迭代重建再投影的外插算法对重建图像是有效的。  相似文献   

6.
本文从畸变的Born近似的微扰技术出发,给出了利用广义Radon变换和Fourier积分算子的理论反演介质间断性的原理.将声学的广义Radon变换与经典Radon变换进行类比,近似地导出了声学广义Radon变换的反演公式. 本文对于反射地震学的情况,提出了一种拟线性化方法,考虑了成象点的一次散射场,从某种程度上减少了Born近似对弱散射的苛求. 利用同一模型的理论记录和物理实验记录的反演计算结果对提出的方法进行了验证,并讨论了进一步提高成象精度的方法.  相似文献   

7.
数据严重缺少的Radon变换的反演   总被引:6,自引:0,他引:6       下载免费PDF全文
在许多实际问题中,由于客观条件的限制经常会遇到数据不全的Radon变换的反演问题.在地球物理学的一些领域中更会遇到数据严重缺少的Radon变换的反演问题──有限个点源及有限个方向投影问题.本文讨论此类反演问题,并给出具体的数值计算方法.  相似文献   

8.
本文将线性层析与Herglotz-Wiechert反演公式进行了比较.二维和三维速度结构的层析反演使用了沿射线的线积分,而且可以以Radon交换的形式写出.对球对称介质,Radon变换可简化为Abel交换.因而对于直的射线路径,走时的Abel交换正是一维球对称介质的层析算法.Herglotz-Wiechert,公式应用地震走时资料反演一维地球结构,且通过Abel交换用精确的射线轨进推导.这是有历史意义的,因为它意味若特定的类似于层析的一种方法从本世纪初就应用在地反学中了(Her目心,1907,,八“h口t,1910).用数值例子比较了He叼otz一Wiechert算法和沿宜射线的线性层析.由Her川侧tZ刃犷沮忱h.时算法于特定条件下(无低速层)可应用于非直射线路径,它与层析的联系在评价推广到布曲射线中的层析结果的非唯一性方面非常有用.  相似文献   

9.
分区连续函数的Radon变换的高精度反演法   总被引:4,自引:0,他引:4       下载免费PDF全文
Radon变换及其反演是CT技术的数学基础,常用的反演算法虽然在函数比较光滑的情况下都能获得较满意的结果,但当函数有较强的奇性时精度就比较差.对函数在每一子区域中是光滑,但在不同于区域交界处有间断的情况,本文给出了有高精度的反演算法.  相似文献   

10.
3D地震数据不规则采样缺失重建是地震勘探数据处理流程中的重要问题.本文提出了一种基于具有保幅特性的非均匀高阶抛物Radon变换(NHOPRT)地震数据重建方法.在最小二乘反演方程中引入Delaunay三角网格剖分来计算空间不规则加权系数,从而获得最接近完整规则数据的高阶抛物Radon变换域系数.在用SVD求解反演方程过程中,利用高阶抛物Radon变换算子在频率域为指数函数,具有线性可分解特性,将二维空间的高阶抛物Radon变换算子分解为两个独立的一维空间变换算子,减小了变换算子的矩阵大小,从而很大程度地提高了计算效率.理论模型和实际地震数据重建测试证明了本文方法的有效性以及实用性.  相似文献   

11.
即使采用分辨率很高的双曲Radon变换,对速度各向异性发育介质及长偏移距情况下的地震数据,其Radon域内能量仍不收敛.为了克服此难题,我们在Radon变换的积分路径中考虑了非双曲走时的影响,通过引入非双曲时差公式中的各向异性非椭圆率η参数,可以准确描述出长偏移距条件下来自同一层位的时距曲线,并推导了由偏移距、慢度、非椭圆率三参数控制的积分曲线正反变换公式,我们称之为各向异性Radon变换.离散化求解时,各向异性Radon变换是时变的,频率域快速算法已不适用,本文采用了最优相似系数加权Gauss-Seidel迭代算法,保持其计算精度的同时也有较高的计算效率.将此方法应用在模型数据以及实际长偏移距海上地震数据的多次波压制处理中,收到了较好的处理效果.  相似文献   

12.
可控源电磁法具有分辨率高及抗干扰能力强等特点,是一种重要的地电磁勘探方法.目前,可控源电磁法的高精度正演计算一直是其核心研究问题之一.传统积分方程法一般采用近似积分公式、简单矩形网格和近似的奇异性体积分计算技术,制约了体积分方程法处理复杂地下异常体的能力,降低了计算精度.针对上述问题,本文基于完全积分公式、四面体非结构化网格和奇异体积分的精确解析解来高精度求解复杂可控源电磁模型的正演响应.首先,从电场积分公式出发,推导了可控源电磁问题满足的积分方程;其次,借助于非结构化四面体网格离散技术,实现了地下复杂异常体的有效模拟.最后,利用散度定理把强奇异值体积分转换为一系列弱奇异性的面积分公式,并通过推导获得了这些弱奇异性的面积分公式的解析解,从而最终实现三维可控源电磁问题的高精度积分求解.以块状低阻体地电模型为测试模型,采用本文提出的积分方程方法获得的数值解与其他公开数值算法解进行对比分析,其对比结果具有高度的吻合性,验证了算法的正确性;同时,设计了球状及复杂地电模型进行算法收敛性测试,进一步验证算法的正确性以及能够处理地下复杂模型的能力.  相似文献   

13.
In this paper, linearized tomography and the Herglotz-Wiechert inverse formulation are compared. Tomographic inversions for 2-D or 3-D velocity structure use line integrals along rays and can be written in terms of Radon transforms. For radially concentric structures, Radon transforms are shown to reduce to Abel transforms. Therefore, for straight ray paths, the Abel transform of travel-time is a tomographic algorithm specialized to a one-dimensional radially concentric medium. The Herglotz-Wiechert formulation uses seismic travel-time data to invert for one-dimensional earth structure and is derived using exact ray trajectories by applying an Abel transform. This is of historical interest since it would imply that a specialized tomographic-like algorithm has been used in seismology since the early part of the century (seeHerglotz, 1907;Wiechert, 1910). Numerical examples are performed comparing the Herglotz-Wiechert algorithm and linearized tomography along straight rays. Since the Herglotz-Wiechert algorithm is applicable under specific conditions, (the absence of low velocity zones) to non-straight ray paths, the association with tomography may prove to be useful in assessing the uniqueness of tomographic results generalized to curved ray geometries.  相似文献   

14.
3D Radon变换及其反变换是X-CT三维图像重建理论的核心,它在其他许多学科领域也有广泛应用。3D Radon变换的表达式是一个三重积分,按照定义直接计算相当费时。为此,研究一种新的快速的方法实现3D Radon变换,对X-CT图像重建理论及相关领域的发展有重要意义。本文以算法仿真常用椭球模型为基础,通过求解椭球模型与空间任意平面的面积,实现了用解析的方法快速得到模型的Radon变换,进一步比较了它与传统方法的优缺点,最后根据Radon反变换重建出原物体模型;计算机仿真结果验证了这种方法的正确。  相似文献   

15.
瞬变电磁场时域格林函数解   总被引:3,自引:3,他引:0       下载免费PDF全文
近源时间域电磁场具有信号强、探测深度大、精度高等优点,但传统勘探电磁场理论中偶极子近似在近源会引起较大误差,导致这一优势的发挥受到了制约.开展直接时间域电磁场解析式研究,是解决这一问题的途径之一.本文提出在点电荷微元假设下,引入时域格林函数,求取瞬变电磁场时间域解析解.采用积分运算法,把电磁场阻尼波动方程的求解问题转化为求其格林函数积分形式解的问题;建立辅助路径解决奇点问题,利用复分析中的约当引理、留数定理和广义函数等理论和方法,推导计算出时间域格林函数的时空四重广义积分.得到达朗贝尔方程的直接时域格林函数精确解析式,与传统方法“比拟”出的公式具有相同的形式,验证了本文推导的时域格林函数解析公式的正确性;推导出扩散方程的直接时间域解析解.通过与时变点电荷源时间域的电磁响应近似表达式进行对比,得出本文所推导的公式计算精度较高的结论;建立了全空间回线源瞬变电磁场问题的直接时间域求解公式.为解决全场区瞬变电磁场精细探测直接时域解析问题提供了基础理论.  相似文献   

16.
地球构造反演问题的新途径   总被引:3,自引:0,他引:3       下载免费PDF全文
本文讨论地球内部构造反演问题的某些新途径.其内容如下:地球构造反问题与固有值反问题;反散射问题中介质间断性的成象与因果广义Radon变换,包含地球构造反问题的新提法,反散射问题的线性化,古典Radon变换与广义Radon变换,线性化反问题的渐近解,渐近解和偏移格式.  相似文献   

17.
Modelling the theoretical response of several important geophysical systems involves the solution of Poisson's equation with homogeneous Neumann boundary conditions (i.e. a zero normal gradient) imposed over either open or closed surfaces. A simple integral equation solution to this problem is derived from first principles. It is applicable to both types of surface and in this respect represents an improvement on existing integral equation techniques. However, the present surface integral equation displays a strong singularity of order 1/R3 which requires an appropriate interpretation for its implementation. A comparison of some numerical results with analytical data taken from the literature demonstrates that the proposed integral equation technique is suitably robust, accurate and efficient for practical application in geophysical interpretation.  相似文献   

18.
在文献[1]中我们利用复变函数论的方法只解决了二维的曲面上的位场转换问题.在本文中,我们利用单层位势理论将它化为一个奇异积分方程求解的问题,并且提出了用奇点分离法来解此奇异积分方程,这对二维和三维情况都适用.通过理论模型的计算,结果的精度是令人满意的.  相似文献   

19.
A numerical method is proposed to accurately and efficiently compute a direct steady-state solution of the nonlinear Richards equation. In the proposed method, the Kirchhoff integral transformation and a complementary transformation are applied to the governing equation in order to separate the nonlinear hyperbolic characteristic from the linear parabolic part. The separation allows the transformed governing equation to be applied to partially- to fully-saturated systems with arbitrary constitutive relations between primary (pressure head) and secondary variables (relative permeability). The transformed governing equation is then discretized with control volume finite difference/finite element approximations, followed by inverse transformation. The approach is compared to analytical and other numerical approaches for variably-saturated flow in 1-D and 3-D domains. The results clearly demonstrate that the approach is not only more computationally efficient but also more accurate than traditional numerical solutions. The approach is also applied to an example flow problem involving a regional-scale variably-saturated heterogeneous system, where the vadose zone is up to 1 km thick. The performance, stability, and effectiveness of the transform approach is exemplified for this complex heterogeneous example, which is typical of many problems encountered in the field. It is shown that computational performance can be enhanced by several orders of magnitude with the described integral transformation approach.  相似文献   

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

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