地球物理学进展 ›› 2018, Vol. 33 ›› Issue (5): 1796-1806.doi: 10.6038/pg2018BB0444

• 固体地球物理及空间物理学(大气、行星、地球动力学、重磁电及地震学、地热学) • 上一篇    下一篇

井间各向异性多波走时联合层析成像:以波场模拟拾取走时为例

何雷宇1(),白超英1,2,*()   

  1. 1. 长安大学地质工程与测绘学院地球物理系, 西安 710054
    2. 长安大学计算地球物理研究所, 西安 710054
  • 收稿日期:2017-12-22 修回日期:2018-08-06 出版日期:2018-10-20 发布日期:2019-01-11
  • 通讯作者: 白超英 E-mail:heleiyu258@outlook.com;baicy@chd.edu.cn
  • 作者简介:何雷宇,男,1994年2月生,甘肃秦安人,长安大学地质资源与地质工程专业博士研究生,主要从事地震正、反演方法及应用研究.

Multi-phase traveltime tomography in crosshole anisotropic media:taking the picked traveltime from wavefileld as an example

HE Lei-yu1(),BAI Chao-ying1,2,*()   

  1. 1. Department of Geophysics, School of Geology Engineering and Geomatics, Chang’an University, Xi’an 710054, China;
    2. Institute of Computational Geophysics, Chang’an University, Xi’an 710054, China;
  • Received:2017-12-22 Revised:2018-08-06 Online:2018-10-20 Published:2019-01-11
  • Contact: Chao-ying BAI E-mail:heleiyu258@outlook.com;baicy@chd.edu.cn

摘要:

为了解决复杂地学模型中正、反演问题(如起伏地表、不规则界面和不规则异常体),本文实现了三角网格剖分下各向异性介质模型中分区多步最短路径射线追踪算法,结合共轭梯度法求解带约束的阻尼最小二乘问题,结合直达波和反射波走时实现了二维各向异性中同时反演四个弹性参数的算法.鉴于各向异性介质中相、群速度关于各弹性参数(a11、a13、a33、a44)偏导数随相角不同而有所不同(每个敏感核函数随相角的变化而变化),而且偏导数的值有正值也有负值.据此,在qP波反演时我们对敏感核函数(Jacob偏导矩阵元素)进行了归一化处理.另外,为了解决qSV波反演时出现的多解性问题,我们提出了分步反演的策略.同时为了验证方法的有效性和正确性,文中采用从地震波场数值模拟中拾取的走时作为观测走时,数值模拟结果表明:波场模拟观测走时反演结果与射线法模拟观测走时反演结果类似,验证了反演过程中正、反演方法的正确性;qP波反演时对敏感核函数归一化处理后的反演结果要明显优于未做归一化的反演结果;qSV波分步反演的结果要优于四个弹性参数同时反演的结果.

关键词: 三角网格模型, 各向异性介质, 弹性参数反演, 井间层析成像, 核函数归一化, 最短路径

Abstract:

The cell in anisotropic seismic tomography have always been regular cell for a long time. Although the regular cell model are simple to understand and easy to realize by computer, it also has some bad aspects. Such as, it is hard to describe complex geological models, and the model cannot be divided into different parts by using regular cell. In order to overcome the difficulty in forward and inversion problem in complex geological model, including undulated topography and irregular subsurface interface, we in this paper realize a multistage triangular shortest-path method for multi-phase seismic ray tracing, and combine a nonlinear inversion solver (constrained and damped least squares problem solved by a conjugated gradient method) to simultaneously invert four elastic parameters in anisotropic media by using combined direct and reflected arrival times. For trading off in different parameter updating, we normalize specific kind of sensitive function of qP-wave inversion, because the partial derivatives of phase velocity respective to different elastic parameters (a11, a13, a33, a44) varying with different phase slowness angles, and the partial derivatives have positive and negative values. In addition, in order to solve the multi-solution problem of qSV-wave inversion, a strategy of stepwise inversion is proposed. For testing the efficiency and correctness of proposed inversion method, we use the picked traveltimes in synthetic seismograms as the observed traveltimes in the inversion. The results show that both inverted images by taking the picked traveltimes in synthetic seismograms as the observed data and by taking predicted traveltimes using ray tracing method as observed data have a similar imaging pictures, which verify the efficiency and correctness of the proposed inversion method. Furthermore, the inverted results by normalized kernel functions is better than that of not normalized kernel functions and the result of stepwise inversion is better than four elastic parameters simultaneous inversion in qSV-wave inversion.

Key words: triangular cell parameterization, anisotropic TI media, elastic parameter inversion, crosshole tomography, kernel function normalization, shortest path method

中图分类号: