刘明忱
- 作品数：6 被引量：32H指数：4
- 供职机构：吉林大学地球探测科学与技术学院更多>>
- 发文基金：国家自然科学基金国家科技重大专项更多>>
- 相关领域：天文地球更多>>
- 基于delta波包叠加的时间域深度偏移被引量：4
- 2016年
- delta波包可由高斯波束经傅里叶逆变换得到,是高斯波束在时空域的表达.它最早出现在合成理论地震图的研究中,本文将其应用于偏移领域.通过delta波包叠加表达时间域格林函数,可将高斯波束偏移由频率域转换到时间域,再结合Rayleigh积分和激励时间成像条件,本文给出了基于delta波包叠加的深度偏移算法.该偏移算法可在时间域直接计算,但因包含褶积运算,成像时将耗费大量的计算时间.针对这一问题,本文提出了将褶积简化为乘积的近似公式.近似后的偏移算法,不仅保留了高斯波束偏移的优点,而且计算效率得到显著提升.文中通过两个数值算例验证了上述结论.
- 石秀林孙建国孙辉刘明忱刘志强
- 关键词：高斯波束
- 2D共炮时间域高斯波束偏移被引量：12
- 2015年
- 针对传统射线方法在奇异区成像精度不高,而2D频率域高斯波束叠前深度偏移需要计算成像点处每个频率的格林函数,影响计算效率的问题,本文通过使用复走时代替实走时,改变频率域下成像公式的积分顺序,给出了在时间域下进行高斯波束偏移的方法和计算公式.本文使用复杂数值模型验证了2D时间域高斯波束叠前偏移方法的正确性,并同传统射线偏移成像结果做了对比.对比结果表明时间域高斯波束偏移在成像精度上优于传统射线偏移.
- 高成孙建国齐鹏孙辉刘志强刘明忱
- 关键词：时间域叠前深度偏移
- 基于快速推进与波前构建的联合3D走时计算（英文）被引量：3
- 2017年
- 3D地震波走时计算是偏移、反偏移、层析等诸多地震勘探技术中的重要中间步骤。快速推进法计算3D地震波走时具有高效率、稳定性及适应能力强的特点,但快速推进法在震源附近区域的计算精度不高,降低了整个走时算法的计算精度。本文提出了一种联合3D走时计算方法来解决这一问题。该方法在震源附近小范围内使用计算精度较高的波前构建法计算走时,在剩余区域使用快速推进法计算走时,由于模型中绝大多数网格节点走时是通过快速推进法计算的,故新方法保留了快速推进法高效的特点,同时由于震源附近网格节点走时精度的提高,整个新算法的计算精度相对于快速推进法而言有了较大的改善。文中通过数值分析对上述结论进行了验证并使用三维岩丘模型验证了新方法的稳定性和适应能力。
- 孙辉孙建国孙章庆韩复兴刘明忱刘志强高正辉石秀林
- 关键词：地震波传播走时计算
- 基于自适应网格的仿真型有限差分地震波数值模拟被引量：6
- 2016年
- 在复杂山地和复杂海底条件下,地表和海底的剧烈起伏对地震波数值模拟提出了更高的要求.常规有限差分法采用矩形网格对模型进行网格剖分,由于矩形网格自身的限制,起伏地表或起伏海底只能由一系列阶梯状折线代替,从而引起人为虚假绕射波.此外,在模拟液-固界面的反射波时,如果界面与网格线不一致,则需要更密的网格才能得到精确的结果.为了解决上述问题,本文将自适应网格生成技术引入到起伏海底速度模型的网格剖分中,采用高阶仿真型有限差分法(MFD)对曲线坐标下的声波方程波进行了数值模拟.利用自适应网格生成技术对速度模型进行网格剖分不仅可以准确地描述模型边界,而且可以有效消除虚假绕射波.高阶仿真型有限差分法可以有效压制频散提高计算精度.模型试算结果表明,本文方法对复杂海底模型具有很好的适应性.
- 刘志强孙建国孙辉刘明忱高正辉石秀林
- 关键词：数值频散自适应网格
- 基于波前构建法的时间域深度偏移——delta波包途径被引量：5
- 2016年
- 作为一种基于射线的局部瞬态场,delta波包源自高斯波束,是高斯波束在时空域的对偶表示,具有高斯波束的全部优点和缺陷。基于delta波包叠加的时间域深度偏移,当射线穿过高速岩体时,受折射效应影响,密度降低,进而导致delta波包的分布密度降低,使成像质量变差,甚至无法成像。为了弥补这个缺陷,本文采用波前构建法计算射线路径。波前构建法能够以插入射线的方式保证均匀的射线分布,从而保证delta波包以均匀的分布密度覆盖整个成像靶区,进而提高成像质量。在具体实现上,采用链表结构替代以往使用的数组结构。Sigsbee 2A模型的数值试算表明,利用波前构建法可以改善高速体下方区域的成像质量,而利用链表存储波前信息要比利用数组至少节省9%的CPU耗时。
- 石秀林孙建国孙辉刘明忱刘志强黄兴国
- 关键词：波前构建
- 基于C语言的SEG-Y数据格式分析与转化被引量：3
- 2016年
- 在地震数据资料存储与交换中,SEG-Y数据格式应用最为广泛.然而,由于CPU架构之间的差异,同一个SEG-Y文件并不能在具有不同CPU架构的计算机上直接读取.为了实现在不同CPU架构的计算机之间进行地震数据共享,基于标准的SEG-Y数据格式和C语言可移植性强的特点,对SEG-Y数据格式的结构构成进行了分析,分别给出了不同结构部分详细的转化算法.最常用的两种地震数据存储格式是32位IBM和IEEE浮点数格式.在研究了这两种浮点数存储结构、表示范围、精度分析和转化差异的基础上,分别给出了两种相互转化的算法并比较了不同算法的转化效率.对比文献中给出的其他编程语言编写的SEG-Y数据某部分转化程序,给出了完整的SEG-Y数据格式转化的C语言代码.实现了SEG-Y数据在不同CPU架构的计算机中的正确读取、改写与存储.良好的可移植性保证了地震数据的共享与应用.
- 刘明忱孙建国
- 关键词：地震数据C语言IBMIEEE
