| 曲线坐标系程函方程的求解方法研究 |
| |
| 引用本文： | 刘一峰,兰海强.曲线坐标系程函方程的求解方法研究[J].地球物理学报,2012,55(6):2014-2026. |
| |
| 作者姓名： | 刘一峰 兰海强 |
| |
| 作者单位： | 1. 中国科学院地质与地球物理研究所, 岩石圈演化国家重点实验室, 北京 100029;
2. 中海石油(中国)有限公司北京研究中心, 北京 100027;
3. 中国科学院研究生院, 北京 100049 |
| |
| 摘 要： | 笛卡尔坐标系中经典的程函方程在静校正、叠前偏移、走时反演、地震定位、层析成像等许多地球物理工作都有应用,然而用其计算起伏地表的地震波走时时却比较困难.我们通过把曲线坐标系中的矩形网格映射到笛卡尔坐标系的贴体网格推导出了曲线坐标中的程函方程,此时,曲线坐标系的程函方程呈现为各向异性的程函方程(尽管在笛卡尔坐标系中介质是各向同同性的).然后尝试用求解各向同性程函方程的快速推进法和Lax-Friedrichs快速扫描算法来分别求解该方程.数值试验表明未加考虑各向异性程函方程与各向同性程函方程的差别而把求解各向同性程函方程的快速推进法直接拓展到曲线坐标中的程函方程的做法是错误的,而Lax-Friedrichs快速扫描算法总能稳定地求解曲线坐标系的程函方程,进而有效地处理了地表起伏的情况,得到稳定准确的计算结果.
|
| 关 键 词： | 程函方程 起伏地表 各向异性 快速推进法 Lax-Friedrichs快速扫描算法 |
| 收稿时间： | 2012-02-05 |
Study on the numerical solutions of the eikonal equationin curvilinear coordinate system |
| |
| LIU Yi-Feng
,
LAN Hai-Qiang.Study on the numerical solutions of the eikonal equationin curvilinear coordinate system[J].Chinese Journal of Geophysics,2012,55(6):2014-2026. |
| |
| Authors: | LIU Yi-Feng LAN Hai-Qiang |
| |
| Institution: | 1. State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China;
2. Beijing Research Center of CNOOC, Beijing 100027, China;
3. Graduate School of Chinese Academy of Sciences, Beijing 100049, China |
| |
| Abstract: | The classical eikonal equation is commonly used in Cartesian coordinate system for problems that involve static correction,prestack migration,earthquake location and seismic tomography,but is less effective for calculating travel times in an earth model that has an irregular surface.We transform the eikonal equation to a curvilinear coordinate system that conforms with the boundaries of the domain.The eikonal equation in the curvilinear coordinate system displays the mathematical form of an anisotropic eikonal equations(even though the medium is isotropic in the Cartesian coordinate system).Then,we try to solve the anisotropic eikonal equations with two numerical methods,a simple modification of the fast marching method and Lax-Friedrichs fast sweeping method,respectively.Through a number of different numerical experiments we find: one may compute wrong solutions by extending fast marching methods designed for isotropic eikonal equations to anisotropic eikonal equations without taking into account the differences between them;the Lax-Friedrichs fast sweeping method is stable and valid for the eikonal equation in the curvilinear coordinate system,which has a significant meaning in the direction of solution of the eikonal equation in the curvilinear coordinate system. |
| |
| Keywords: | Eikonal equation Topography Anisotropy Fast marching method Lax-Friedrichs fast sweeping method |
| 本文献已被 CNKI 万方数据 等数据库收录！ |
| 点击此处可从《地球物理学报》浏览原始摘要信息 |
| 点击此处可从《地球物理学报》下载免费的PDF全文 |
|
