首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 140 毫秒
1.
本文提出用于地下成像无源地震方法的分析,在该方法中应用环境地震噪声作为地下散射体的照明源。该成像算法能够以递归方式将新的数据融入到成像中,这使成像背景噪声随时间减少。在空间域不相干环境背景噪声的假设前提下推导出成像算法的点-扩展函数。点-扩展函数表征成像的分辨率,即接收排列长度和环境带宽的函数。  相似文献   

2.
像域层析速度建模方法利用偏移道集剩余曲率构建目标函数并迭代更新速度场.迭代更新速度场需反复进行叠前深度偏移,对于三维地震数据,偏移成像计算时间长,迭代更新过程成倍增加层析成像计算时间.本文由剩余时差拾取机理出发,提出隐式剩余时差概念,并建立隐式剩余时差与剩余速度之间的函数关系,利用该关系通过一次计算就可以对剩余速度进行精确求取并对速度模型进行精确更新,该方法避免了常规层析方法需要多次迭代的流程.由模型数据测试及实际数据应用证明,本文提出方法是有效可行的,该方法与常规层析方法进行比较,在保证提高层速度场建模精度的同时有效提高计算效率,为后续偏移成像节省大量时间.  相似文献   

3.
主动源与被动源地震数据插值及联合数据成像   总被引:1,自引:1,他引:0       下载免费PDF全文
本文提出了两种情况下主动源数据和被动源数据的插值方法,并研究了两种数据在偏移成像中的互补效果.基于互相关法被动源数据重构原理,本文提出了结合多域迭代去噪技术的重构方法.提出了两种时间域主动源和被动源数据的插值方法,分别是共炮点域能量匹配插值和共检波点域最小平方匹配插值.然后对获得的主动源和被动源联合地震数据进行叠前深度偏移成像.在被动源活跃度不是很高的地区进行被动源地震勘探时,少量的主动源地震数据可以有效控制和补充被动源数据的成像效果.在稀疏炮点的主动源勘探中,有效利用被动源的信息能够在成像中增加更多的细节信息,提高成像质量.  相似文献   

4.
最小二乘偏移能够消除非规则地震采集、带限子波等因素对偏移结果造成的不良影响,提高成像分辨率和振幅保真度,但高昂的计算成本严重制约了其实用性和应用前景.本文基于射线理论格林函数和地震波走时一阶近似,提出了一种快速的点扩散函数计算方法,不但可以大幅降低点扩散函数的计算成本,还能够适应任意的空间采样以保证点扩散函数的插值精度.以此为基础,本文进一步发展了成像域最小二乘偏移方法,能够高效、灵活地实现地下反射率的迭代更新,并最终获得分辨率更高、照明更加均衡的反演成像结果.文中给出的模型和实际数据的测试结果验证了上述方法的正确性和有效性.  相似文献   

5.
基于不完备投影数据重建的四种迭代算法比较研究   总被引:1,自引:0,他引:1  
在发射成像和穿透成像过程中,由于各种原因会造成投影数据不完备,若仍采用传统解析方法重建出的图像会产生伪影,而迭代算法则可以很好地改善图像质量。本文应用不同的迭代算法分别就真实标准线对测试卡的稀疏投影数据和有限角度投影数据进行重建,定量分析比较它们各自的优缺点,为几种迭代算法的工程或临床应用提供重要参考。  相似文献   

6.
给出一种稳健反演方法(也称迭代加权最小二乘IRLS),该方法对分离数据点极不敏感,因为它可以自动剔除异常的数据点,不需要人工仔细检查和控制数据质量。为了使用Tarantola-Valette的迭代广义离散反演方法进行稳健反演,将基于残差和信噪比的不同判据引入到协方差矩阵(作用相当于权重函数)中。在实际应用中,这个算法的作用相当于震源机制初步测定(PDFM),后者在海啸预警中用于快速估计强震震源参数。要反演的输入数据是长周期面波的频谱,输出的计算结果为地震矩张量,根据这些张量获得地震的震源几何形状和主应力轴。这种算法可以应用于任何非线性(迭代)反演计算中,解决污染数据集的分离数据点的问题。  相似文献   

7.
在工业CT的检测过程中,由于各种原因可能造成投影数据的不完备.在这种情况下,直接利用传统的解析重建方法(FBP)或者统计重建方法(EM)都不能给出理想的重建结果,从切向CT扫描方式获得的就是典型的不完备投影数据.本文利用迭代算法来改善重建图像的质量,首先通过已有的投影数据合理估计缺失的投影数据,并通过不断的迭代来改善这种估计,从而使得重建图像不断逼近真实物体.  相似文献   

8.
由投影重建图像的对称网格迭代算法   总被引:1,自引:0,他引:1  
本文对于工业CT检测中常用的代数迭代重建算法提出了改进,利用投影射线之间存在的几何对称结构,提出了图像重建的对称网格迭代算法(简写为SM-IRT).该算法简化了投影系数矩阵的计算,调整了迭代算法逐线校正的迭代顺序.对模拟数据和工业CT实测数据进行了重建图像的数值实验,结果表明:与常规算法比较,本文提出的新算法重建速度快,成像精度高.  相似文献   

9.
刘东甲 《地球物理学报》2017,60(8):3215-3228
本文提出了位场区域-剩余异常分离的空间域递减半径迭代法.由给定半径的圆周上八点位场值的算术平均导出一个新的八点圆周平均公式,它是一个由重磁数据计算区域异常的滤波器.该滤波器的传递函数有一主瓣和多个旁瓣,且半径越大,旁瓣越多,滤波器特性越差.由从大到小不同半径对应的传递函数的连乘积构造了递减半径迭代传递函数,它以大半径为参变量.递减半径迭代传递函数类似低通滤波器,其截止波数与大半径成反比.由递减半径迭代传递函数,给出空间域分离区域异常的递减半径线性迭代法,由重磁数据减去区域异常求得剩余异常.进一步,通过构造非线性修正系数,把空间域递减半径线性迭代法中线性迭代公式变成非线性迭代公式,得到空间域递减半径非线性迭代法.重磁理论模型及安徽泥河铁矿重磁资料的试验表明,空间域递减半径非线性迭代法可有效地压制假异常和高频干扰以及减少异常畸变,并能有效地分离出区域异常和剩余异常.  相似文献   

10.
在工业CT成像应用中,由于复杂结构件形状、结构、组分等特征的限制,使得射线能量与有效厚度不匹配,常规的扫描方式易使投影数据缺失,投影数据完备性较差。为此,本文通过改变轨迹来弥补投影数据缺失,推导任意轨迹的成像模式。首先,在任意轨迹的基础上进行投影几何关系建模,结合空间几何知识完成投影矩阵的刻画;其次,基于迭代重建算法对轨迹的无约束化,进行任意轨迹成像模式下的迭代重建;最后通过小锥角、大锥角以及厚度差异比较大的构件成像仿真实验,验证算法的正确性和可行性。对于厚度差异比较大的工件,该方法相对于传统的圆轨迹扫描,投影数据更完备,CT重建质量高。  相似文献   

11.
最小曲率位场分离方法研究   总被引:9,自引:7,他引:2       下载免费PDF全文
位场分离是位场数据处理和解释中的重点和难点之一.本文给出了单步长非原位和原位两种最小曲率位场分离差分迭代格式,并利用Fourier频谱分析理论研究了这两种迭代格式的收敛性.通过研究表明,单步长非原位迭代格式不收敛,只有单步长原位迭代格式收敛,但单步长原位迭代格式受迭代方向选择的影响,随着迭代次数的增大其影响逐渐消失.根据单步长非原位迭代格式的频谱特点,提出了叠加步长非原位和原位迭代格式,同样利用Fourier频谱分析理论研究了叠加步长非原位和原位迭代格式的收敛性.通过研究认为,一维叠加步长非原位迭代格式收敛,但二维叠加步长非原位迭代格式不收敛;不论是一维或二维,其原位迭代格式均收敛.进一步的理论研究表明,非原位迭代格式的频率响应是一个实偶函数,而原位迭代格式的频率响应是一个复函数;单步长迭代格式的频率响应具有一定的周期性,而叠加步长迭代格式的频率响应无周期性特征;叠加步长迭代格式比单步长迭代格式的收敛性好.  相似文献   

12.
Interferometric redatuming is a data‐driven method to transform seismic responses with sources at one level and receivers at a deeper level into virtual reflection data with both sources and receivers at the deeper level. Although this method has traditionally been applied by cross‐correlation, accurate redatuming through a heterogeneous overburden requires solving a multidimensional deconvolution problem. Input data can be obtained either by direct observation (for instance in a horizontal borehole), by modelling or by a novel iterative scheme that is currently being developed. The output of interferometric redatuming can be used for imaging below the redatuming level, resulting in a so‐called interferometric image. Internal multiples from above the redatuming level are eliminated during this process. In the past, we introduced point‐spread functions for interferometric redatuming by cross‐correlation. These point‐spread functions quantify distortions in the redatumed data, caused by internal multiple reflections in the overburden. In this paper, we define point‐spread functions for interferometric imaging to quantify these distortions in the image domain. These point‐spread functions are similar to conventional resolution functions for seismic migration but they contain additional information on the internal multiples in the overburden and they are partly data‐driven. We show how these point‐spread functions can be visualized to diagnose image defocusing and artefacts. Finally, we illustrate how point‐spread functions can also be defined for interferometric imaging with passive noise sources in the subsurface or with simultaneous‐source acquisition at the surface.  相似文献   

13.
积分法三维电阻率成像   总被引:37,自引:5,他引:37       下载免费PDF全文
二维或是三维电阻率反演成像研究,最关键的环节是在反演系数矩阵即敏感矩阵(或雅可比矩阵)的求取上.本文从微分方程的积分解出发,推导了表达式极为简单的三维雅可比系数矩阵,构造了成像方程.根据反演方程系数的稀疏特征,采用改进的降维高斯赛德尔迭代法来求解该反演方程,并通过内外迭代的结合,对大型稀疏欠定方程能很快收敛,得到可靠的解答.合成数据模型结果及实际资料的成像结果表明积分法不但实现起来极容易,成像结果的精度也相当高.  相似文献   

14.
反馈迭代法压制表面多次波效果分析   总被引:3,自引:0,他引:3  
地震数据处理通常假定反射数据仅由一次波组成,地震勘探资料中的表面多次波通常被视为相干噪音予以压制.基于波动方程预测的反馈迭代法为数据驱动的方法,可有效压制复杂地下介质的表面多次波,通常可由级数展开法和迭代法两种方法实现.文中分别从理论阐述以及应用效果分析方面对基于波动方程预测的反馈迭代法的两种实现途径进行研究,并采用GPU/CPU协同并行加速计算预测表面多次波,改善以往CPU计算效率低的状况.文中对含表面多次波的炮数据和复杂的SMAART模型进行了计算,并将利用级数展开法与迭代法压制表面多次波的效果进行对比分析,结果表明,迭代形式反馈迭代法压制表面多次波的方法效果更佳.由于实现迭代法和级数展开法需要完成的自适应匹配相减的次数不同,迭代法的计算成本略高于级数展开法.  相似文献   

15.
积分迭代法的正则性分析及其最优步长的选择   总被引:1,自引:1,他引:0       下载免费PDF全文
位场积分迭代法是一种解决位场大数据量、大深度向下延拓的有效方法.本文基于Kirsch正则化子理论,推导了积分迭代法对应的正则化滤子函数,并证明积分迭代法为一种求解位场向下延拓不适定反问题的正则化方法.针对积分迭代法迭代步长固定、迭代次数较多影响收敛速度的问题,提出该迭代法最优迭代步长的选择原理.理论模型和实测数据对比分...  相似文献   

16.
三维井间电磁场的正反演计算   总被引:21,自引:7,他引:21       下载免费PDF全文
提出了一种计算体积分方程的方法-改进型局域非线性迭代(MLNI),并用其对三维井间电磁场进行反演计算,该方法将井间大尺度散射体分为近场和远场区域两部分,它们的位置和尺寸均随场点位置的变化而改变,采用局域非线性近似计算近场区域的影响,将远场区域的影响作为外部激励源,采用迭代方法计算,该方法具有计算速度快,所需内存量少,收敛性好的优点,在反演中采用基于MLNI的非线性反演方法,利用这种方法可以将成像区域集中于一定范围内而不是整个三维空间,该反演方法由于考虑到了雅可比矩阵元素的非线性项,因而与传统的Born迭代反演方法相比更为精确和稳定,考虑到信息量和计算机内存的限制,第1次成像采用双重体元分割法进行较粗略的成像,然后缩小成像范围进行第2次较精确的成像,数值计算结果表明,MLNI是一种有效的计算井间大尺度异常体散射场的方法,将该方法用于反演过程能够得到较高分辨率的三维井间电导率图像。  相似文献   

17.
Zohdy方法近似反演三维电阻率测深数据。正演计算采用有限单元法。反演初始模型由测量视电阻率数据给出。通过比较实测视电阻率值和预测模型计算的视电阻率值对数差来修改模型网格电阻率.为了解决任意电极距测深数据的反演,采用大、小双网格剖分。大网格反映地下电性分布情况。小网格用于实际有限元正演计算.在电阻率调整公式中加入一个迭代系数,能够加快收敛速度.并对加5%随机噪声的模型理论视电阻率测深断面数据进行反演,得到的电阻率分布与模型电阻率基本一致.  相似文献   

18.
刘洋  张鹏  刘财  张雅晨 《地球物理学报》2018,61(4):1400-1412
人工地震方法由于受到野外观测系统和经济因素等的限制,采集的数据在空间方向总是不规则分布.但是,许多地震数据处理技术的应用(如:多次波衰减,偏移和时移地震)都基于空间规则分布条件下的地震数据体.因此,数据插值技术是地震数据处理流程中关键环节之一.失败的插值方法往往会引入虚假信息,给后续处理环节带来严重的影响.迭代插值方法是目前广泛应用的地震数据重建思路,但是常规的迭代插值方法往往很难保证插值精度,并且迭代收敛速度较慢,尤其存在随机噪声的情况下,插值地震道与原始地震道之间存在较大的信噪比差异.因此开发快速的、有效的迭代数据插值方法具有重要的工业价值.本文将地震数据插值归纳为数学基追踪问题,在压缩感知理论框架下,提出新的非线性Bregman整形迭代算法来求解约束最小化问题,同时在迭代过程中提出两种匹配的迭代控制准则,通过有效的稀疏变换对缺失数据进行重建.通过理论模型和实际数据测试本文方法,并且与常规迭代插值算法进行比较,结果表明Bregman整形迭代插值方法能够更加有效地恢复含有随机噪声的缺失地震信息.  相似文献   

19.
基于牛顿迭代法和遗传算法的CSAMT近场校正   总被引:1,自引:0,他引:1       下载免费PDF全文
栾晓东  底青云  雷达 《地球物理学报》2018,61(10):4148-4159
由于可控源音频大地电磁法(CSAMT)采用人工场源在大大增加信号强度的同时也带来了在近区产生非平面波效应的问题,其表现特征为在双对数坐标中卡尼亚电阻率曲线呈45°上升,即所谓的近场效应.本文首先提出了基于牛顿迭代法求解视电阻率方程的近场效应校正方法,通过对二层、三层理论模型试算验证了该方法的有效性;为进一步增强对噪声的压制能力,本文进一步引入了遗传反演的思想,将全频段误差最小作为目标函数,采用自适应正则化方法引入模型参数梯度最小作为稳定因子,由遗传算法求解得到校正后视电阻率,通过对理论模拟数据加随机噪声的校正结果,表明该方法在实现近场校正的同时能剔除噪声干扰,得到的结果与理论模拟曲线吻合度较高.通过对内蒙古曹四夭钼矿区实测CSAMT数据处理结果表明,牛顿法和遗传算法均能明显校正CSAMT因近场效应引起的假高阻异常,能更好地反映地电结构特征.  相似文献   

20.
介绍了利用地震计对阶跃电流标定或稳态正弦电流标定的响应波形获取地震计传递函数的方法。对阶跃标定电流的响应波形进行傅立叶变换(FFT)和扣除阶跃信号频谱的运算获得地震计的幅频特性;从稳态正弦标定电流的响应波形可直接提取幅频特性,然后用高斯牛顿迭代法拟合出地震计传递函数的各个系数。该方法的最大优点是只需要已知幅频特性或其有限离散点的数值就可以拟合出地震计的传递函数,较目前常用的需要同时已知幅频和相频特性的拟合方法更为方便和实用。  相似文献   

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

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