Kirchhoff向下延拓法在海洋合成多道地震走时反演中的应用*

doi:10.11978/2019087

热带海洋学报 ›› 2020, Vol. 39 ›› Issue (4): 80-90.doi: 10.11978/2019087CSTR: 32234.14.2019087

Kirchhoff向下延拓法在海洋合成多道地震走时反演中的应用*

张茂传^1,^2,³(), 徐敏^1,³(), 赵旭^1,^2,³, 张佳政^1,³, 查财财^1,^2,³, VithanaM.V.P.^1,², 狄会哲^1,^2,³, 曾信^1,³

1.中国科学院边缘海与大洋地质重点实验室, 南海海洋研究所, 南海生态环境工程创新研究院, 广东广州 510301
2.中国科学院大学, 北京 100049
3.南方海洋科学与工程广东省实验室(广州), 广东广州 511458

收稿日期:2019-09-16 修回日期:2020-01-17 出版日期:2020-07-20 发布日期:2020-07-27
通讯作者: 徐敏
作者简介:张茂传(1994—), 男, 江西赣州人, 硕士研究生, 主要从事海洋地球物理研究。E-mail: zhangmaochuan@scsio.ac.cn
基金资助:
国家自然科学基金项目(41676044);国家自然科学基金项目(91858207);国家自然科学基金项目(41890813);中国科学院战略性先导科技专项(XDA13010105);南方海洋科学与工程广东省实验室(广州)人才团队引进重大专项(GML2019ZD0205)

Application of Kirchhoff downward-continued method to travel-time inversion of synthetic multi-channel seismic data

Maochuan ZHANG^1,^2,³(), Min XU^1,³(), Xu ZHAO^1,^2,³, Jiazheng ZHANG^1,³, Caicai ZHA^1,^2,³, M. V. P. Vithana^1,², Huizhe DI^1,^2,³, Xin ZENG^1,³

1. Key Laboratory of Ocean and Marginal Sea Geology, South China Sea Institute of Oceanology, Innovation Academy of South China Sea Ecology and Environmental Engineering, Chinese Academy of Sciences, Guangzhou 510301, China
2. University of Chinese Academy of Sciences, Beijing 100049, China
3. Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou 511548, China

Received:2019-09-16 Revised:2020-01-17 Online:2020-07-20 Published:2020-07-27
Contact: Min XU
Supported by:
Foundation item: National Natural Science Foundation of China(41676044);Foundation item: National Natural Science Foundation of China(91858207);Foundation item: National Natural Science Foundation of China(41890813);The Strategic Priority Research Program of the Chinese Academy of Sciences(XDA13010105);Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory(Guangzhou)(GML2019ZD0205)

摘要/Abstract

摘要：

海洋多道地震数据建模和成像是获取洋壳速度和构造信息的重要手段。海水层的存在使得多道地震拖缆接收到的折射波走时信息仅仅存在于较远的炮检距, 近炮检距被强振幅海底反射波覆盖, 制约了走时数据拾取和反演效果。本文基于波动方程的Kirchhoff积分法, 成功实现了多道地震数据向下延拓, 获取到了更大炮检距区间的初至折射波走时拾取, 并将其应用于洋中脊新生洋壳2A/2B层的合成多道地震数据走时反演。比较向下延拓前后的走时拾取范围及走时反演结果表明, 向下延拓法能够保持地震波场的运动学和动力学特征不变, 在共炮集数据的更大炮检距范围内进行初至折射走时拾取, 从而增加反演的数据选择和浅层射线覆盖, 反演结果能更加准确地分辨出洋壳2A/2B层界面, 并得到更高的分辨率和更准确的速度结构剖面。

关键词: 海洋多道地震, Kirchhoff积分, 向下延拓, 走时反演, 洋壳2A/2B层

Abstract:

Velocity model building and seismic imaging from marine multi-channel seismic (MCS) data are essential to obtain the structural information of oceanic crust. Existence of the seawater layer makes first-arrival refraction phases recorded by streamer only existed in far offset, while most of them in near offset are masked by the strong amplitude seafloor reflection phase, affecting travel-time picking and inversion. Based on the Kirchhoff integral approach to wave equation, we successfully implement the downward continuation (DC) method to the MCS data, and obtain first-arrival travel-time picks within a wider offset range. The method is applied to the synthetic MCS data simulated based on layer 2A/2B structure of newborn crust along a mid-ocean ridge, and first-arrival picking range and inversion result are compared before and after the DC. The DC method can expand offset range for first-arrival refraction picking and increase the ray coverage in the shallow depth, without changing the kinematics and dynamic characteristics of seismic wave field. The inversion results can more precisely distinguish the layer 2A/2B interface, and obtain higher resolution and more accurate velocity profile.

Key words: marine multi-channel seismic, Kirchhoff integration, downward continuation, travel-time inversion, oceanic crust layer 2A/2B

中图分类号:

P738.4

张茂传, 徐敏, 赵旭, 张佳政, 查财财, VithanaM.V.P., 狄会哲, 曾信. Kirchhoff向下延拓法在海洋合成多道地震走时反演中的应用*[J]. 热带海洋学报, 2020, 39(4): 80-90.

Maochuan ZHANG, Min XU, Xu ZHAO, Jiazheng ZHANG, Caicai ZHA, M. V. P. Vithana, Huizhe DI, Xin ZENG. Application of Kirchhoff downward-continued method to travel-time inversion of synthetic multi-channel seismic data[J]. Journal of Tropical Oceanography, 2020, 39(4): 80-90.

图/表 6

图1

图2

图4

图5

图6

图7

参考文献 37

[1]	耿建华, 黄海贵, 马在田, 1996. Kirchhoff积分波场延拓基准面静校正方法研究[J]. 同济大学学报, 24(6):665-669.
	GENG JIANHUA, HUANG HAIGUI, MA ZAITIAN, 1996. Kirchhoff wavefield extrapolation datuming[J]. Journal of Tongji University, 24(6):665-669 (in Chinese with English abstract).
[2]	金丹, 阎贫, 唐群署, 等, 2011. Kirchhoff波场延拓在OBC记录海水层基准面校正中的应用[J]. 热带海洋学报, 30(6):84-89.
	JIN DAN, YAN PIN, TANG QUNSHU, et al, 2011. Application of Kirchhoff integral wave field extrapolation to water layer datuming for OBC record[J]. Journal of Tropical Oceanography, 30(6):84-89 (in Chinese with English abstract).
[3]	王祥春, 王延峰, 夏常亮, 等, 2012. Kirchhoff积分法OBS数据地震波场延拓[J]. 现代地质, 26(6):1231-1236.
	WANG XIANGCHUN, WANG YANFENG, XIA CHANGLIANG, et al, 2012. Extrapolating the OBS data using the Kirchhoff integral method[J]. Geoscience, 26(6):1231-1236 (in Chinese with English abstract).
[4]	ARNULF A F, HARDING A J, KENT G M, et al, 2014a. Anatomy of an active submarine volcano[J]. Geology, 42(8):655-658.
[5]	ARNULF A F, HARDING A J, KENT G M, et al, 2014b. Constraints on the shallow velocity structure of the Lucky Strike Volcano, Mid-Atlantic Ridge, from downward continued multichannel streamer data[J]. Journal of Geophysical Research: Solid Earth, 119(2):1119-1144.
[6]	ARNULF A F, HARDING A J, SINGH S C, et al, 2012. Fine-scale velocity structure of upper oceanic crust from full waveform inversion of downward continued seismic reflection data at the Lucky Strike Volcano, Mid-Atlantic ridge[J]. Geophysical Research Letters, 39(8):L08303.
[7]	ARNULF A F, HARDING A J, SINGH S C, et al, 2013. Nature of upper crust beneath the Lucky Strike volcano using elastic full waveform inversion of streamer data[J]. Geophysical Journal International, 196(3):1471-1491.
[8]	ARNULF A F, SINGH S C, HARDING A J, et al, 2011. Strong seismic heterogeneity in layer 2A near hydrothermal vents at the Mid-Atlantic ridge[J]. Geophysical Research Letters, 38(13):L13320.
[9]	BERRYHILL J R, 1979. Wave-equation datuming[J]. Geophysics, 44(8):1329-1344.
[10]	BERRYHILL J R, 1984. Wave-equation datuming before stack[J]. Geophysics, 49(11):2064-2066.
[11]	BEVC D, 1995. Imaging under rugged topography and complex velocity structure[M]. PhD book, Stanford University.
[12]	BEVC D, 1997. Flooding the topography: Wave-equation datuming of land data with rugged acquisition topography[J]. Geophysics, 62(5):1558-1569.
[13]	BOLTE J F B, 2001. Velocity independent CFP redatuming: A strategy for subsalt imaging[C]//2001 SEG Annual Meeting. San Antonio, Texas: Society of Exploration Geophysicists, 901-904.
[14]	CHRISTESON G L, PURDY G M, FRYER G J, 1994. Seismic constraints on shallow crustal emplacement processes at the fast spreading East Pacific Rise[J]. Journal of Geophysical Research: Solid Earth, 99(B9):17957-17973.
[15]	COLLINS J A, BLACKMAN D K, HARRIS A, et al, 2009. Seismic and drilling constraints on velocity structure and reflectivity near IODP Hole U1309D on the Central Dome of Atlantis Massif, Mid-Atlantic Ridge 30°N[J]. Geochemistry, Geophysics, Geosystems, 10(1):Q01010.
[16]	GHOSAL D, SINGH S C, MARTIN J, 2014. Shallow subsurface morphotectonics of the NW Sumatra subduction system using an integrated seismic imaging technique[J]. Geophysical Journal International, 198(3):1818-1831.
[17]	HARDING A J, ARNULF A F, BLACKMAN D K, 2016. Velocity structure near IODP Hole U1309D, Atlantis massif, from waveform inversion of streamer data and borehole measurements[J]. Geochemistry, Geophysics, Geosystems, 17(6):1990-2014.
[18]	HARDING A J, KENT G M, BLACKMAN D K, et al, 2007. A new method for MCS refraction data analysis of the uppermost section at a Mid-Atlantic Ridge core complex[C]// American Geophysical Union, fall meeting 2007: S12A-03.
[19]	HARDING A J, KENT G M, ORCUTT J A, 1993. A multichannel seismic investigation of upper crustal structure at 9°N on the East Pacific Rise: Implications for crustal accretion[J]. Journal of Geophysical Research: Solid Earth, 98(B8):13925-13944.
[20]	HENIG A S, BLACKMAN D K, HARDING A J, et al, 2012. Downward continued multichannel seismic refraction analysis of Atlantis Massif oceanic core complex, 30°N, Mid-Atlantic ridge[J]. Geochemistry, Geophysics, 13(5):Q0AG07.
[21]	KENT G M, SINGH S C, HARDING A J, et al, 2000. Evidence from three-dimensional seismic reflectivity images for enhanced melt supply beneath mid-ocean-ridge discontinuities[J]. Nature, 406(6796):614-618. doi: 10.1038/35020543 pmid: 10949299
[22]	LARKIN S P, LEVANDER A, 1996. Wave-equation datuming for improving deep crustal seismic images[J]. Tectonophysics, 264(1-4):371-379.
[23]	MARJANOVIĆ M, FUJI N, SINGH S C, et al, 2017. Seismic signatures of hydrothermal pathways along the east pacific rise between 9°16′ and 9°56′N[J]. Journal of Geophysical Research: Solid Earth, 122(12):10241-10262.
[24]	MCMECHAN G A, CHEN H W, 1990. Implicit static corrections in prestack migration of common-source data[J]. Geophysics, 55(6):757-760.
[25]	QIN YANFANG, SINGH S C, 2018. Insight into frontal seismogenic zone in the Mentawai locked region from seismic full waveform inversion of ultra-long offset streamer data[J]. Geochemistry, Geophysics, Geosystems, 19(11):4342-4365.
[26]	RESHEF M, 1991. Depth migration from irregular surfaces with depth extrapolation methods[J]. Geophysics, 56(1):119-122.
[27]	SHTIVELMAN V, CANNING A, 1988. Datum correction by wave-equation extrapolation[J]. Geophysics, 53(10):1311-1322.
[28]	SINGH S C, HARDING A J, KENT G M, et al, 2006. Seismic reflection images of the Moho underlying melt sills at the East Pacific Rise[J]. Nature, 442(7100):287-290. pmid: 16855587
[29]	VERA E E, DIEBOLD J B, 1994. Seismic imaging of oceanic layer 2A between 9°30′N and 10°N on the East Pacific Rise from two-ship wide-aperture profiles[J]. Journal of Geophysical Research: Solid Earth, 99(B2):3031-3041.
[30]	VITHANA M V P, XU MIN, ZHAN XU, et al, 2019. Geological and geophysical signatures of the East Pacific Rise 8°-10° N[J]. Solid Earth Sciences, 4(2):66-83. doi: 10.1016/j.sesci.2019.04.001
[31]	WANG HAIYANG, SINGH S C, CALANDRA H, 2014. Integrated inversion using combined wave-equation tomography and full waveform inversion[J]. Geophysical Journal International, 198(1):430-446. doi: 10.1093/gji/ggu138
[32]	WHITMORE N D, 1983. Iterative depth migration by backward time propagation[M]// SEG Technical Program Expanded Abstracts. Society of Exploration Geophysicists, 1983: 382-385.
[33]	WIGGINS J W, 1984. Kirchhoff integral extrapolation and migration of nonplanar data[J] Geophysics, 49(8):1239-1248. doi: 10.1190/1.1441752
[34]	XU MIN, 2012. Advanced geophysical studies of accretion of oceanic lithosphere in Mid-Ocean Ridges characterized by contrasting tectono-magmatic settings[D]. Cambridge, MA: Massachusetts Institute of Technology.
[35]	XU MIN, STEPHEN R A, CANALES J P, 2017. Waveform modeling of the seismic response of a mid-ocean ridge axial melt sill[J]. Marine Geophysical Researches, 38(4):373-391. doi: 10.1007/s11001-017-9303-x
[36]	YILMAZ O, LUCAS D, 1986. Prestack layer replacement[J]. Geophysics, 51(7):1355-1369. doi: 10.1190/1.1442186
[37]	ZELT C A, BARTON P J, 1998. Three-dimensional seismic refraction tomography: A comparison of two methods applied to data from the Faeroe Basin[J]. Journal of Geophysical Research: Solid Earth, 103(B4):7187-7210.

Kirchhoff向下延拓法在海洋合成多道地震走时反演中的应用*

Application of Kirchhoff downward-continued method to travel-time inversion of synthetic multi-channel seismic data

RichHTML

PDF (PC)

赞

可视化

摘要/Abstract

引用本文

使用本文

图/表 6

参考文献 37

相关文章 1

Metrics

本文评价

推荐阅读 0