波前扩展有限差分地震波走时算法的C++语言描述

波前扩展有限差分地震波走时算法的C⁺⁺语言描述

杨昊^1,2，孙建国^1,2，韩复兴^1,2

1.吉林大学地球探测科学与技术学院，长春 130026;2.国土资源部应用地球物理综合解释理论开放实验室，长春 130026

收稿日期:2006-09-28 修回日期:1900-01-01 出版日期:2007-05-26 发布日期:2007-05-26
通讯作者: 杨昊

An C⁺⁺ Language Program Implement of Traveltime Calculation of Expanding Wavefronts Finite-Difference Method

YANG Hao^1,2, SUN Jian-guo^1,2, HAN Fu-xing^1,2

1.College of GeoExploration Science and Technology, Jilin University, Changchun, 130026,China;

Received:2006-09-28 Revised:1900-01-01 Online:2007-05-26 Published:2007-05-26
Contact: YANG Hao

摘要/Abstract

摘要： 回顾了波前扩展有限差分地震波走时算法及其主要优缺点，分析了该算法在程序实现过程中存在的一些问题。针对这些问题，给出了波前扩展有限差分地震波走时算法的C⁺⁺语言描述。目的是:（1）确保程序的稳定性。使用体波、首波、散射波相结合的走时算法，在保证走时计算精度的同时，较好地解决了前述算法在局部速度差异较大的情况下可能会出现负数开平方问题；（2）提高程序的执行效率。给出一种波前点的“列队”式存储方法，有效提高了波前最小走时点的查找效率；（3）为波前扩展有限差分地震波走时算法构建一个方便、可行的程序框架。对原始速度模型增加了边界，不仅方便了程序实现，而且省去了程序用于判断和处理模型边界所花费的时间。

关键词: 走时, 波前扩展, 程函方程, 列队, C⁺⁺

Abstract: We first review the finite-difference scheme proposed by Qin in 1988 which solves the eikonal equation along expanding wavefronts and analyze its advantages and disadvantages. We also show some problems about its programming implement. Then we show an implement method in the C⁺⁺ programming language, which aims at solving the problems. (1) To ensure the stability of the program, we also consider the refracted and diffracted waves to calculate the traveltimes. The new method improves the local algorithmic stabilities of Qin’s method without losing the precise. (2) To improve the efficiency of the program, we propose a queue storage structure of the wavefront points. (3) We build a convenient and doable program structure for the expanding wavefronts method. We append four boundaries to the original velocity model, which not only makes programming simple but also save the time dealing with the model boundaries.

Key words: traveltime, expanding wavefronts, eikonal equation, queue structure, C⁺⁺

中图分类号:

杨昊，孙建国，韩复兴. 波前扩展有限差分地震波走时算法的C⁺⁺语言描述[J]. J4, 2007, 37(3): 615-0619.

YANG Hao, SUN Jian-guo, HAN Fu-xing. An C⁺⁺ Language Program Implement of Traveltime Calculation of Expanding Wavefronts Finite-Difference Method[J]. J4, 2007, 37(3): 615-0619.

参考文献

相关文章 8

[1]	孙建国, 苗贺. 基于Chebyshev走时逼近的三维多次反射射线计算[J]. 吉林大学学报(地球科学版), 2018, 48(3): 890-899.
[2]	孙建国, 李懿龙, 孙章庆, 苗贺. 基于模型参数化的地震波走时与射线路径计算[J]. 吉林大学学报(地球科学版), 2018, 48(2): 343-349.
[3]	李桐林, 张镕哲, 朴英哲. 大地电磁测深与地震初至波走时交叉梯度反演[J]. 吉林大学学报(地球科学版), 2015, 45(3): 952-961.
[4]	曲昕馨，李桐林，王飞. 基于数字图像分割法的跨孔雷达走时层析成像[J]. 吉林大学学报(地球科学版), 2014, 44(4): 1340-1347.
[5]	张廉萍, 刘洪. 基于波前重建和李代数积分的地震波走时计算[J]. J4, 2010, 40(6): 1463-1470.
[6]	杨昊, 孙建国, 韩复兴, 马淑芳. 基于完全三叉树堆排序的波前扩展有限差分地震波走时快速算法[J]. J4, 2010, 40(1): 188-194.
[7]	孙建国,何洋. 基于波前构建的射线追踪：一种Java实现[J]. J4, 2007, 37(4): 814-0820.
[8]	滕吉文. 地球深部壳-幔边界的层束精细结构与物理属性研究[J]. J4, 2006, 36(01): 1-0023.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed