地球物理学进展 ›› 2016, Vol. 31 ›› Issue (4): 1628-1632.doi: 10.6038/pg20160429

• 应用地球物理学Ⅰ • 上一篇    下一篇

二维无限频率初至走时反演井间地震实际资料

廖建平1,2,4, 刘和秀1, 黄建平3, 赵延林1, 郑桂娟5, 戴世鑫1, 杨天春1, 王齐仁1,4   

  1. 1. 湖南科技大学地质系, 湘潭 411201;
    2. 煤炭资源与安全开采国家重点实验室(中国矿业大学), 北京 100083;
    3. 中国石油大学(华东)地学院, 青岛 266580;
    4. 湖南科技大学页岩气资源利用湖南省重点实验室, 湘潭 411201;
    5. 中国石油东方地球物理公司研究院, 涿州 072751
  • 收稿日期:2015-11-12 修回日期:2016-04-13 出版日期:2016-08-20 发布日期:2016-08-20
  • 通讯作者: 刘和秀,女,1977,博士,讲师,主要从事地震波正演及波形反演工作.(E-mail:lhx2001peking@sohu.com) E-mail:lhx2001peking@sohu.com
  • 作者简介:廖建平,男,1976,博士,副教授.研究方向为复杂介质中地震波传播及数值模拟;走时反演;二维、三维频率-空间域、时间-空间域声波、弹性波全波形反演方法研究.(E-mail:bobljp@163.com)
  • 基金资助:

    国家自然科学基金(41274126、51274097)、煤炭资源与安全开采国家重点实验室(中国矿业大学)开放基金(SKLCRSM11KFB01)、 湖南省自然科学基金项目(12JJ6035)、中国博士后科学基金(2012M541116)、湖南省教育厅科研项目(14C0462)以及湖南科技大学教研教改项目(G31385)联合资助.

2D infinite frequency first-arrival travel time inversion on crosswell real field seismic data

LIAO Jian-ping1,2,4, LIU He-xiu1, HUANG Jian-ping3, ZHAO Yan-lin1, ZHENG Gui-juan5, DAI Shi-xin1, YANG Tian-chun1, WANG Qi-ren1,4   

  1. 1. Department of Geology of Hunan University of Science and Technology, Xiangtan 411201, China;
    2. State Key Laboratory of Coal Resources and Safe Mining, China University of Mining & Technology, Beijing 100083, China;
    3. School of Geosciences of China University of petroleum, Qingdao 266580, China;
    4. Hunan Provincal Key Laboratory of Shale Gas Resource Utilization, Hunan University of Science and Technology, Xiangtan 411201, China;
    5. GRI.BGP.CNPC, ZhuoZhou 072751, China
  • Received:2015-11-12 Revised:2016-04-13 Online:2016-08-20 Published:2016-08-20

摘要:

使用Zelt和Barton的方法,通过一个计算效率高的有限差分求解eikonal方程,正演计算走时和射线路径.使用最小二乘QR分解法,求解稀疏线性系统方程组.使用正则化层析反演,结合用户给定的最小的、最平坦和最平滑的扰动限制,每一个加权因子随深度变化.结合数据残差和模型粗糙度的最小化,为数据残差提供一个最平滑的近似模型.该反演方法为非线性反演,需要一个初始模型,在每一次迭代时,需要计算新的射线路径.使用二维初至走时数据,对某油田二维井间地震实际资料进行无限频率初至走时层析反演.将反演所得到的速度与井的测井速度曲线相比较,二者吻合程度较高,表明该反演方法所得速度的分辨率比较高.证实了二维无限频率初至走时层析反演可以为全波形反演提供一个分辨率较高的长波长速度模型,从而为全波形反演井间地震实际资料提供了一个比较可靠的初始速度模型.

Abstract:

Using Zelt and Barton method, through an efficient implementation of an existing method for solving the eikonal equation by finite difference solution, the forward modeling calculate travel time and ray path. Using LSQR method, to solve the sparse linear system of equations. Using regularized tomographic inversion in which a combination with the user given minimal disturbance limits, the flatest and smoothest, each weighted factor changs with depth. Combined with the data errors and model roughness of minimizing, providing the smoothest model approximation for the data residuals. The inversion method is a nonlinear inversion, it needs an initial model, it needs to compute a new ray path in each iteration step. We use a two-dimensional first-arrival travel time data, on a certain oil field 2D crosswell real field seismic data to infinite frequency travel time inversion(IFTT). The match degree between the inversion of velocity and shot and receiver's well logging curve is high, which shows that the high resolution of the first-arrival travel time inversion method. It confirms that this 2D infinite frequency travel time inversion can provide a high resolution of long wavelength velocity model, which provides a reliable initial velocity model for the full waveform inversion of crosswell real field seismic data.

中图分类号: