共查询到20条相似文献，搜索用时 15 毫秒
1.
针对层析成像模型矩形网格剖分存在的一些问题,提出了复杂结构三角网最小走时射线追踪层析成像方法。以Delaunay三角剖分的优化准则,根据模型点、线、面的几何结构关系,进行三角网格剖分;采用三角网波行面扩展法计算声波初至走时,追踪出任意接收点至激发点处的射线路径,射线路径包含了从接收点至源点的坐标、所在三角单元等信息,形成矩阵方程。应用稳定性较好的联合迭代重构技术求解矩阵方程,得到模型的速度分布。数值模拟结果表明,三角网射线层析成像方法分辨率高,成像结果更接近实际结构形态,有效地解决了矩形网剖分对复杂区域网格参数化灵活性差,速度界面描述精度低等问题。 相似文献
2.
《地球物理学进展》2017,(1)
海底地层速度结构是识别海洋天然气水合物储层的直接依据,本文应用地震反射走时层析成像建立了海底地层速度模型.采用不规则网格对模型进行离散化,使速度单元与反射界面单元完全耦合;利用基于不规则单元波前扩展和走时插值的射线追踪方法,精确确定反射射线路径和反射波走时;在反演中同时使用先验约束、平滑约束、归一化和正则化技术,提高了层析反演的稳定性和结果的可靠性.对南海北部神狐海域SH2井附近的二维地震测线资料,利用多域人机交互法拾取了反射走时,用反射走时层析成像方法获得了SH2井附近含水合物目标区的速度结构.该反演结果与测井声波速度和钻探结果一致,其中的高速带对应水合物储层,表明反射走时层析成像能够有效地得到海洋天然气水合物储层的速度结构,为海洋天然气水合物储层识别提供依据. 相似文献
3.
为克服矩形网速度间断面及边界描述精度差的缺点,本文提出了一种基于波前最小走时单元的三角网射线追踪算法.该算法用三角单元剖分速度模型,把与波前面相邻、处于振动状态的三角单元称为波前单元,从源点所在单元开始构造由波前单元组成的波前单元域;1)选域中最小走时三角单元为当前单元;2)它先与相邻三角单元振动传递,再进行相邻单元各自内部节点的传递、并将处于振动状态的相邻单元扩展为波前单元;3)然后确认当前单元的振动状态、决定其在波前单元域中的去留,以此三步循环进行波前单元域的扩展及节点的最小走时和次级源搜索,直至波前单元域在剖分区域的消亡;最后利用各节点最小走时和次级源,通过最小走时搜索,拾取从接收点到源点的射线路径.数值模拟表明,该算法的模型剖分灵活性强,对速度间断面的描述准确,追踪结果精度高. 相似文献
4.
在复杂山地和复杂海底条件下,地表和海底的剧烈起伏对地震波数值模拟提出了更高的要求.常规有限差分法采用矩形网格对模型进行网格剖分,由于矩形网格自身的限制,起伏地表或起伏海底只能由一系列阶梯状折线代替,从而引起人为虚假绕射波.此外,在模拟液-固界面的反射波时,如果界面与网格线不一致,则需要更密的网格才能得到精确的结果.为了解决上述问题,本文将自适应网格生成技术引入到起伏海底速度模型的网格剖分中,采用高阶仿真型有限差分法(MFD)对曲线坐标下的声波方程波进行了数值模拟.利用自适应网格生成技术对速度模型进行网格剖分不仅可以准确地描述模型边界,而且可以有效消除虚假绕射波.高阶仿真型有限差分法可以有效压制频散提高计算精度.模型试算结果表明,本文方法对复杂海底模型具有很好的适应性. 相似文献
5.
地震层析成像技术是研究地层速度结构的一个重要手段，该技术是利用震源到接收点的地震波走时或波形的观测数据，结合建立的数学物理模型来重建地层速度结构。而层析成像正演模拟的精度将直接影响反演出的速度与地层真实速度的拟合度，因此找出一种精度较高的层析成像正演模拟算法是非常重要的。本文针对三维层状介质，通过网格化模型，采用最小走时射线追踪方法展开层析成像正演研究，分别采用三角形网格、矩形网格和六边形网格进行模型参数化，并对比分析各自正演结果的精度，总结出针对三维层状介质模型的最佳网格划分方式。 相似文献
6.
在许多地震反演和偏移成像方法中,都要涉及到射线路径和旅行时的计算.本文将波前面三角形网格剖分和三维波前重建法射线追踪技术结合使用,实现了射线路径和旅行时的准确快速计算.三维波前重建法射线追踪过程中可以保证稳定合理的射线密度,克服了常规射线追踪方法存在阴影区的问题.波前面三角形网格剖分在描述和拆分波前面时更加准确有效,而且不需太多的网格数目,从而提高了射线追踪的精度和效率.该方法在三维复杂构造成像方面有独特的优势,目前在实际的Kirchhoff 偏移中的已经有相关应用. 相似文献
7.
近地表地层与人类生产生活密切相关,利用地震层析成像方法准确重建浅部地壳速度结构有助于开展高精度地震勘探、探查浅部矿产资源、规避潜在自然灾害,并利于城市地下空间建设.中国大陆地表条件复杂,尤其中西部盆岭结合带地形起伏剧烈,对浅部地壳精确速度建模构成严重挑战.本文系统论述了地震层析成像领域基于高频近似理论的走时成像方法和有限频层析成像方法,阐明两类方法的基本原理、存在问题和发展方向等.依据正演走时有无显式射线追踪,基于高频近似理论的走时成像方法分为传统走时层析成像方法和无射线路径的走时层析成像方法.基于射线追踪的传统走时层析成像方法,在浅层速度强烈变化时,因存在阴影区或多路径现象引起成像失真,严重影响成像效率;而无射线路径的层析成像方法通过程函方程走时场的正传和逆传直接计算敏感核,并利用伴随状态法获得目标函数的梯度,具有快速、稳健的优点.以上两种基于地震射线高频近似理论的走时成像方法由于未考虑地震波频率的带限性,存在波散射、波前愈合及反演约束差等问题.有限频层析成像方法克服了射线理论"无限高频"假设所带来的弊端,已成为重要的研究方向之一.该类方法主要分为射线有限频层析成像方法和基于波动方程的有限频层析成像方法.射线有限频层析成像方法能够提高成像的分辨率,但在方法本质上仍依赖于射线理论,较难处理较复杂的波现象问题;基于波动方程的有限频层析成像方法能准确处理复杂地质问题,提高成像可靠性并能以图像形式直观展示地球内部地震波的速度结构分布,但是该方法在实际应用中强烈依赖于数据中的低频信息及较精确的初始速度模型,其推广应用仍需进一步探索. 相似文献
8.
跨孔雷达全波形反演是一种使用全波形信息反演两钻孔之间地下信息的层析成像技术.常规的层析成像反演大部分采用射线追踪方法,其中基于初至时的射线追踪方法可以反演出速度剖面（介电常数剖面）,基于最大振幅的层析成像可以反演出衰减剖面（电导率剖面）.常规射线追踪方法有许多不足,究其原因是该方法仅使用了小部分的信号信息.为了进一步提高成像分辨率,本文全面推导了全波形跨孔雷达层析成像反演方法,该方法利用雷达波全幅度相位信息能够反演出地下高分辨率的介电常数和电导率图像.本文通过基于局域网的分布式并行算法,有效地解决了巨量数据正演计算问题.文中首先建立了基于单轴各向异性介质完全匹配层的时间域有限差分二维正演算法,进而通过应用包括时间维度在内的全波场信息与残场逆向传播的全波场信息乘积来计算梯度方向,通过求取以步长为自变量的目标函数的极值确定步长公式,并提出以第一次介电常数反演作为同步反演的初始模型,能够有效提高收敛速度.本文对多组模型进行成像实验,取得了较好的反演效果. 相似文献
9.
变阻尼约束层析成像及其在VSP资料中的应用(英文) 总被引：1，自引：0，他引：1
初至波走时层析成像已经取得了广泛的应用,然而,由于观测系统的限制,射线在模型中分布不均匀,导致层析结果的分辨能力不足。变阻尼约束方法应用不均匀的先验信息来匹配不均匀的数据分布,可以减小速度模型校正量与射线覆盖程度的相关性。本文将变阻尼约束方法应用于初至波旅行时层析成像中,并将平滑约束方法加入正则化方程组中来避免单独使用变阻尼约束带来的不稳定性,利用阿尔法滤波器对反演中间迭代结果进行平滑和去噪,采用LSQR算法求解线性方程组来提高收敛速度和压制误差传递。本文应用上述层析成像算法对VSP观测系统进行速度反演,分别应用于检测板速度模型数据和实际VSP资料速度反演中,结果表咀变阻尼约束层析成像可以改善射线不均匀覆盖带来的影响,从而提高速度反演结果的质量;VSP资料检波点附近的速度反演结果可靠性高。 相似文献
10.
地震初至波中包含着丰富的近地表速度结构信息,如何分阶段、分尺度地利用这些信息进行近地表速度建模是地震勘探中的一个关键问题.在速度反演的不同阶段,综合利用初至波中的不同信息(如走时、包络和波形等)进行联合反演,可以有效地降低反演对初始模型的依赖程度,提高近地表速度模型的反演精度.为此,本文提出了一种统一基于波动方程正演引擎的初至波多信息联合反演方法,该方法同时匹配观测和模拟的初至波走时、包络和波形信息.在不同的反演阶段选择不同的权重因子调节不同信息的权重,这样不仅降低了反演对初始速度模型的依赖,而且自然地实现了多尺度反演.在每轮反演迭代中,一次正演模拟的波场同时应用于初至波走时、包络和波形匹配,无需额外的射线追踪.同时,联合反演方法在一定程度上缓解了串联反演中目标函数漂移问题,提高了近地表速度建模的精度. 相似文献
11.
地震走时层析成像是反演地层各向异性参数分布的有效方法,但是关于地震各向异性介质走时层析成像的研究并不多,其技术远远没有达到成熟的阶段.在野外数据采集时,地表反射波观测方式相对井间和垂直地震剖面观测方式的成本更低,利用qP反射波走时反演各向异性参数具有更加广泛的实用价值.本文实现的TI介质地震走时层析成像方法结合了TI介质反射波射线追踪算法、走时扰动方程和非线性共轭梯度算法,它可以对任意强度的TI介质模型进行反演,文中尝试利用qP反射波走时重建TI介质模型的参数图像.利用qP反射波对层状介质模型和块状异常体模型进行走时反演,由于qP波相速度对弹性模量参数和Thomsen参数的偏微分不同,所以可以分别反演弹性模量参数和Thomsen参数.数值模拟结果表明:利用qP反射波可以反演出TI介质模型的弹性模量参数与Thomsen参数,不同模型的走时迭代反演达到了较好的收敛效果,与各向同性介质走时反演结果相比较,各向异性介质走时反演结果具有较好的识别能力. 相似文献
12.
本文使用最小二乘线性迭代反演方法对跨孔雷达直达波初至时数据进行反演,每次迭代过程中,用有限差分法求解走时程函方程,并用高精度快速推进方法(HAFMM)进行波前扩展,通过追踪波前避免了进行射线追踪.为了验证该方案,我们对三组合成数据进行了测试,分析了单位矩阵算子、一阶差分算子和拉普拉斯算子等三种不同模型参数加权算子对模型的约束和平滑效果;讨论了FMM和HAFMM对反演精度的影响;测试了LSQR,GMRES和BICGSTAB等三种矩阵反演算法的反演效果.此外,我们还对一组野外实测数据进行了反演,对比了基于本方案以及基于平直射线追踪和弯曲射线追踪的走时层析成像反演效果.对比分析结果表明,使用拉普拉斯算子和HAFMM进行反演能较好地进行目标体重建,而三种矩阵反演方法对反演效果的影响差别不大;并且通过对波前等时线图的分析可以定性地判断异常体的性质和位置;而在对实测数据目标体的重建上,本方案能达到甚至优于弯曲射线算法的重建效果. 相似文献
13.
基于纵横波初至走时数据的层析成像方法越来越广泛地被应用于揭示不同构造域壳幔速度结构特征.我们从同一地质体的纵横波速度属性相关这一基本思想出发,提出一种相关反演成像的方法:纵横波速度反演交替进行,在迭代反演过程中每通过一次反演获得相应的纵波速度(或横波速度)结构后,更新相应的纵横波速度比模型以及相应的横波(或纵波)速度反演的初始模型,然后继续开展后续横波(或纵波)速度反演工作.在反演过程中依据纵横波速度的相关性信息和射线路径长度将走时残差以不同权重分配到射线路径经过的单元,依据网格节点周围平均的慢度扰动更新速度模型.正反演过程分别基于有限差分走时计算方法和反投影成像方法.两种典型模型试验表明,该技术应用于上地壳速度结构反演成像过程,可有效提高反演结果的可靠性,在很大程度上避免了常规单独反演纵波和横波速度过程容易带来的畸变和失真.该方法应用于重建青藏高原西部札达—泉水沟深地震测深(DSS)剖面下方的上地壳速度结构,揭示出与青藏高原西缘板块碰撞相关的上地壳速度结构特征. 相似文献
14.
对于微地震正演模拟,本文以射线追踪的原理为基础,对两点间的射线追踪问题进行了研究,应用二分算法、改进二分算法和微变网格算法对水平层状匀速模型、弯曲层状匀速模型和复杂地质模型进行射线追踪,使得计算效率和适用范围都得到了很大的改善.文中对每种算法误差范围和计算效率进行了对比验证,对于不同的地质模型,选用合适的算法才能在计算速度和精度上得到双重保证,最后正演模拟了多波三分量记录.在模型建立上引入了超薄层概念,并在前人模拟的直达波、透射波、反射波基础上拟了折射波,使正演模拟的多波场信息更丰富.文中的应用实例及模型结果表明:与二分法相比,改进二分法能够对弯曲界面进行射线追踪,并能保证结果的精度.弯曲层状模型中,改进二分法与微变网格法相比计算速度有显著提高,能够应用到资料的反演中. 相似文献
15.
基于华南及周边地区106个宽频带地震台站多年记录的MS≥5.0中浅源地震事件,开展瑞利面波层析成像和速度结构特征研究,获得了华南大陆及南海北部地区10~100s瑞利波群速度分布图像和典型剖面下方地壳上地幔速度结构,为理解该地区构造演化和深部过程提供约束.考虑到实际地震射线的覆盖情况以及华南地区主要构造的主体展布特征,本文同时采用传统的规则网格剖分和平行主要构造走向的非规则网格剖分方法,分别进行分格频散反演,开展了不同参数化方案对反演结果影响的对比分析研究.基于瑞利面波层析成像结果,进行了典型剖面横波速度结构反演,重建了华南地区由内陆至南海海域主要构造单元的壳幔横波速度结构.研究结果表明,扬子和华夏块体地壳上地幔结构特征差异显著,扬子块体地壳和岩石圈厚度均大于华夏地块,且扬子块体上地幔顶部速度较华夏块体低,岩石圈厚度在雪峰山造山带下方发生过渡和转换;南海北部陆缘和南海海盆上地幔速度较高且形态相对完整,表现为非火山型大陆边缘和已停止扩张海盆的壳幔结构特征. 相似文献
16.
《地球物理学进展》2020,(2)
网格层析成像采用离散网格点代替空间连续介质模型,通过更新网格点速度值对模型进行迭代优化,从而达到速度误差收敛的目的.常规网格层析是基于规则矩形网格进行网格剖分,横向及纵向网格采用等间距步长,该方法在网格剖分过程中没有考虑地层构造特点.本文提出一种基于层位约束的不规则网格层析成像方法.该方法通过截取时窗的策略对特殊构造部分采取加密网格的策略,从而实现根据地质构造指导网格剖分.对于简单的地质构造,本文提出方法与常规方法效果相当.对于复杂构造地层,尤其对于小尺度特殊地质构造体,本文提出方法更具针对性,对于特殊构造体速度建模精度更高.由实际数据应用证明,对于复杂地质构造,本文方法能够大幅提高速度场建模精度,为叠前深度偏移能够准确成像保驾护航. 相似文献
17.
18.
为了在复杂地表条件下实现地震波走时计算,提出了一种基于线性插值和窄带技术的走时计算新方法.其中,线性插值用于局部走时计算,窄带技术用于局部波前捕获和追踪.为了逼近起伏地表,采用三角网和矩形网相结合的方法对速度模型进行剖分.为了得到局部走时计算公式,利用费马（Fermat）原理和关于入射点位置的限定条件.有关编程实践和数值试验表明：新方法不仅可以有效、灵活地处理地表高程的剧烈变化,而且还具有很好的适应性和稳定性,得到的计算结果满足波前传播规律. 相似文献
19.
井间地震数据直达波走时层析成像 总被引：1，自引：0，他引：1
本文利用矢量射线追踪正演模拟技术计算地震波直达波传播的路径及走时,进而利用射线走时及路径的内插,发展了弯曲射线迭代反演技术.该方法可用来重建井间地层的速度图像.基于所发展的方法,我们对较为复杂的典型地质模型进行了井间速度重建.结果表明该方法是一项快速、高精度的走时层析成像技术. 相似文献
20.
采用新近提出的多次波射线追踪正演算法,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,分析讨论了利用多震相走时资料进行联合反演成像的方法及技术.考虑到不同震相走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数; 由于同时反演速度模型和反射界面起伏中不同模型参数变化对走时影响程度的不同, Jacobi偏导矩阵元素中引入了不同参数的归一化因子; 另外,为了克服射线密度过大(或过小)区域速度模型的过度(或欠)更新问题,反演算法中引入了等权射线密度的概念.几种数值模拟实例表明(含噪声敏感性试验): 多波走时的联合或同时反演成像技术是一种提高走时成像空间分辨率,进而降低重建模型失真度行之有效的方法. 相似文献