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

doi:10.6038/pg20160429

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

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

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

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-ping^1,2,4, LIU He-xiu¹, HUANG Jian-ping³, ZHAO Yan-lin¹, ZHENG Gui-juan⁵, DAI Shi-xin¹, YANG Tian-chun¹, WANG Qi-ren^1,4

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

摘要/Abstract

摘要：

使用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.

中图分类号:

P631

廖建平, 刘和秀, 黄建平, 赵延林, 郑桂娟, 戴世鑫, 杨天春, 王齐仁. 2016. 二维无限频率初至走时反演井间地震实际资料. 地球物理学进展, 31(4): 1628-1632. doi: 10.6038/pg20160429.

LIAO Jian-ping, LIU He-xiu, HUANG Jian-ping, ZHAO Yan-lin, ZHENG Gui-juan, DAI Shi-xin, YANG Tian-chun, WANG Qi-ren. 2016. 2D infinite frequency first-arrival travel time inversion on crosswell real field seismic data. Progress in Geophysics. 31(4): 1628-1632. doi: 10.6038/pg20160429.

参考文献

Abgrall R., Benamou J. D. 1999. Big raytracing and eikonal solver on unstructured grids: Application to the computation of a multivalued traveltime field in the Marmousi model[J]. Geophysics, 64(1): 230-239.
Bendar, J. B., Shin, C., et al. 2007. Comparison of waveform inversion, part 2: Phase approach[J]. Geophysical Prospecting, 55(4): 465-475.
Bernasconi. G., Drufuca. G.. 2001. 3-D traveltimes and amplitudes by gridded rays[J]. Geophysics, 66(1): 277-282.
Bunks C., Saleck F. M., Zaleski S., et al. 1995. Multiscale seismic waveform inversion[J].Geophysics, 60(5): 1457-1473.
Coultrip. R. L. 1993. High-accuracy wavefront tracing traveltime calculation[J].Geophysics, 58(2): 284-292.
Eppstein. M. J., Dougherty. D. E.. 1998. Optimal 3-D traveltime tomography[J].Geophysics, 63(3): 1053-1061.
Ettrich. N., Gajewski. D. 1998. Traveltime computation by perturation with FD-eikonal solvers in isotropic and weakly anisotropic media[J].Geophysics, 63(3): 1066-1078.
Fischer, R. and Lees, J. M. 1993. Shortest path ray tracing with sparse graphs[J].Geophysics, 58(7): 987-996.
Guust Nolet. 1991. SEISMIC TOMOGRAPHY[M](in Chinese). Feng R, Hao J Q Trans. Beijing: Geological Publishing House.
Langan, R.T., Lerche, et al. 1985.Tracing of rays through heterogeneous media: An accurate and efficient procedure[J].Geophysics, 50(9): 1456-1465.
Liao J P, Liu H X, Wang H Z, et al. 2011b. High resolution acoustic wave full waveform velocity inversion in frequency space domain-theoretical model[J]. Progress in Geophys. (in Chinese), 26(5): 1690-1695. DOI: 10.3969/j.issn.
Liu Yu-zhu, Ding Kong-yun, Dong Liang-guo. 2010.The Dependence of first arrival travel time tomography on initial model[J].OGP, 45(4): 502-511.
Liu Yu-zhu, Dong Liang-guo.2007.Analysis of influence factor of first-breaks tomography[J].OGP, 42(5): 544-553.
Luo. Y., Schuster. G. T. 1991.Wave-equation traveltime inversion[J].Geophysics, 56(5): 645-653.
Mao W. J., Stuart G. W. 1997. Rapid multi-wave-type ray tracing in complex 2-D and 3-D isotopic media[J].Geophysics, 62(1): 298-308.
Matthew A. Brzostowski and George A. McMechan. 1992. 3-D tomographic imaging of near-surface seismic velocity and attenuation[J].Geophysics, 57(3): 396-403.
Nemeth T, Normark E, Qin F. 1997. Dynamic smoothing in crosswell traveltime tomography[J].Geophysics, 62(1): 168~176.
Pica A., Diet J. P., Tarantola A. 1990. Nonlinear inversion of seismic reflection data in a laterally invariant medium[J].Geophysics, 55(3): 284-293.
Pratt, R. G., and Worthington, M. H. 1990. Inverse theory applied to multi-source cross-hole tomography. Part I: Acoustic wave-equation method[J].Geophysical Prospecting, 38(3): 287-310.
Pratt R. G and Goulty N R. 1991. Combining wave-equation imaging with traveltime tomography to form high-resolution images from crosshole data[J].Geophysics, 56(2): 208-224.
Pulliam, J. and Snieder, R.. 1998. Ray perturbation theory, dynamic ray tracing and the determination of Fresnel zones[J].Geophys. J. Int, 135(2): 463-469.
Qin. F.H., Luo. Y., Olsen. K. B., et al. 1992. Finite-difference solution of the eikonal equation along expanding wavefronts[J].Geophysics, 57(3): 478-487.
Sava. P., Fomel. S., et al. 2001. 3-D traveltime computation using Huygens wavefront tracing[J].Geophysics, 66(3): 883-889.
Scales J. A. 1987. Tomographic inversion via the conjugate gradient method[J].Geophysics, 52(2): 179-185.
Schuster G. T., and Quintus-Bosz A.. 1993. Wavepath eikonal traveltime inversion: Theory[J].Geophysics, 58(9): 1314-1323.
Sebudandi C.. Toint P. L.. 1993. Non-linear optimization for seismic traveltime tomography[J].Geophys. J. Int., 115(3): 929-940.
Sethian J. A., Popovic A. M.. 1999.3-D traveltime computation using the fast marching method[J].Geophysics, 64(2): 516-523.
Simmons. J. L., Backus. M. M.. 1992. Linearized tomographic inversion of first-arrival times[J].Geophysics, 57(11): 1482-1492.
Soupios. P. M., Papazachos. C. B. et al. 2001. Nonlinear 3-D traveltime inversion of crosshole data with an application in the area of the Middle Ural Mountains[J].Geophysics, 66(2): 627-636.
Tarantola. A., 1984, Linearized inversion of seismic reflection data[J].Geophys. Prosp., 32(6): 998-1015.
Tarantola. A.. 1984. Inversion of seismic reflection data in acoustic approximation[J].Geophysics, 49(8): 1259-1266.
Van Avendonk, H. J. A., Harding, A. J., et al. 2001. Hybrid shortest path and ray bending method for traveltime and raypath calculations[J].Geophysics, 66(2): 648-653.
Vasco D. W. and Peterson J. E., Majer E. L.. 1995. Beyond ray tomography: wavepaths and Fresnel volumes[J].Geophysics, 60(6): 1790-1804.
Williamson. P. R.. 1991. A guide to the limits of resolution imposed by scattering in ray tomography[J].Geophysis, 56(2): 202-207.
Woodward M. J.. 1992. Wave-equation tomography[J].Geophysics, 57(1): 15-26.
Xie Chun, Liu Yuzhu, Dong Liangguo, Yang Jizhong.First arrival tomography with finite frequency adjoint-state method based on acoustic wave equation[J].OGP, 2015, 50(2): 274-282
Yao. Z. S., Osypov. K. S., Robert R. G.. 1998. Traveltime tomography using regularized recursive least squares[J].Geophys. J. Int., 135(2): 545-553.
Yang Qin-yong, Hu Guang-hui, Wang Li-xin. 2014. Research status and development trend of full waveform inversion[J]. Geophysical Prospecting for Petroleum, 53(1): 77-83.
YANG Wu-yang, WANG Xi-wen, YONG Xue-shan et al. 2013.The review of seismic Full waveform inversion method[J].Progress in Geophysics, 28(2): 766-776, doi: 10.6038/pg20130225
Zelt, C.A., P.J. Barton. 1998. Three-dimensional seismic refraction tomography: A comparison of two methods applied to data from the Faeroe Basin[J]. J. Geophys. Res., 103(B4): 7187-7210.
Zelt, C. A., S. Haines, M. H. Powers, et al. 2013. Blind test of methods for obtaining 2D near-surface seismic velocity models from first-arrival traveltimes[J].Journal of Environmental and Engineering Geophysics, 18(3) : 183-194.
Zelt, C.A., Azaria, A., and Levander, A., 2006. 3D seismic refraction traveltime tomography at a groundwater contamination site[J].Geophysics, 71(5) : H67-H78.
Zelt, C.A., Sain, K., Naumenko, J.V., et al. 2003.Assessment of crustal velocity models using seismic refraction and reflection tomography[J].Geophysical Journal International, 153(3): 609-626.
Zhang, J., and Toksoz M, N.. 1998. Nonlinear refraction traveltime tomography[J].Geophysics, 63(5): 1726-1737.
Zhou. H. 1993. Traveltime tomography with a spatial-coherency filter[J].Geophysics, 58(5): 720-726.
廖建平, 刘和秀, 王华忠, 等. 2011b. 高分辨率的频率空间域声波全波形速度反演-理论模型[J]. 地球物理学进展, 26(5): 1690-1695.DOI: 10.3969/j.issn.
诺列. 1991. 地震层析技术[M]. 冯锐, 郝锦绮译.北京: 地质出版社.
刘玉柱, 董良国. 2007. 初至波层析影响因素分析[J].石油地球物理勘探, 42(5): 544-553.
刘玉柱, 丁孔芸, 董良国. 2010.初至波走时层析成像对初始模型的依赖性[J].石油地球物理勘探, 45(4): 502-511.
谢春, 刘玉柱, 董良国, 杨积忠.基于声波方程的有限频伴随状态法初至波走时层析[J].石油地球物理勘探, 2015, 50(2): 274-282.
杨勤勇, 胡光辉, 王立歆.2014. 全波形反演研究现状及发展趋势[J].石油物探, 53(1): 77-83.
杨午阳, 王西文, 雍学善, 等.2013.地震全波形反演方法研究综述[J].地球物理学进展, 28(2): 766-776, doi: 10.6038/pg20130225.

[1]	屈念念,李家斌,张西君,刘俊. 基于重磁资料研究贵州省深部构造特征[J]. 地球物理学进展, 2019, 34(5): 1785-1793.
[2]	安鹏,曹丹平,赵宝银,杨晓利,张明. 基于LSTM循环神经网络的储层物性参数预测方法研究[J]. 地球物理学进展, 2019, 34(5): 1849-1858.
[3]	陈稳,薛国强,陈卫营,周楠楠,卢云飞. SOTEM多分量激电响应特性分析[J]. 地球物理学进展, 2019, 34(5): 1859-1865.
[4]	豆辉,徐逸鹤. 黏性逆时偏移成像研究进展[J]. 地球物理学进展, 2019, 34(5): 1866-1877.
[5]	毛庆辉,王鹏,曾隽. 水力压裂微地震事件定位方法综述[J]. 地球物理学进展, 2019, 34(5): 1878-1886.
[6]	尚新民,王延光,崔庆辉,金昌昆,赵翠霞,滕厚华,赵胜天. 准南缘山前带精细近地表建模技术研究[J]. 地球物理学进展, 2019, 34(5): 1887-1892.
[7]	石战战,夏艳晴,王元君,周怀来,庞溯,池跃龙. 基于L₁-L₁范数稀疏表示的共偏移距道集地震数据重建方法[J]. 地球物理学进展, 2019, 34(5): 1893-1899.
[8]	王迪,陈敏,王健伟,盛志超,严曙梅,马美媛,李东哲,郑欣. 基于T₂和毛管压力分类转换的渗透率计算新方法[J]. 地球物理学进展, 2019, 34(5): 1900-1909.
[9]	王海波,张伟,张宏,肖虎,高云路,顾剑华,唐美斌,韩峰. 高精度可控震源在深反射地震采集中的应用[J]. 地球物理学进展, 2019, 34(5): 1910-1916.
[10]	王静, 张军华, 王延光, 杨勇, 冯德永, 武刚, 黄德峰. 特征值相干理论诠释及效果比较[J]. 地球物理学进展, 2019, 34(5): 1917-1923.
[11]	王维红,包培楠,陈国飞,林弘喆. 基于最小平方的横波速度拟合及应用[J]. 地球物理学进展, 2019, 34(5): 1924-1929.
[12]	王亚,王孟华,田建章,王成泉,彭玲丽,冯小英,周兴海,张浩. 一种适用于深潜山及内幕的有效储层地震预测方法——以廊固凹陷杨税务潜山为例[J]. 地球物理学进展, 2019, 34(5): 1930-1937.
[13]	吴曲波,潘自强,陈金勇,Mubarak,Hani. 利用瞬态瑞雷面波法探测浅层玄武岩三维分布——以沙特Sabkhah Ad Dumathah地区为例[J]. 地球物理学进展, 2019, 34(5): 1938-1944.
[14]	严小丽,康慧敏,王光杰,何国丽,李红领. AMT方法在鳌山卫花岗岩地区深部地热构造勘探中的应用[J]. 地球物理学进展, 2019, 34(5): 1945-1953.
[15]	尹奇峰,潘冬明,郭全仕,何情. 基于快速模拟退火算法的井中微地震事件定位反演[J]. 地球物理学进展, 2019, 34(5): 1954-1961.

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

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

PDF (PC)

可视化

被引次数

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

多维度评价

本文评价