|
|
Fast first-arrival traveltime tomography of diving waves under rugged surface and its application to static correction |
YANG Hua-Chen1,2( ), GE Da-Ming1, WANG Zhong-Cheng1,2, WANG Lei1, YUAN Yong-Qi1 |
1. Geophysical Research Institute,Shengli Oilfield Company,SINOPEC,Dongying 257022,China 2. Working Station for Postdoctoral Scientific Research of Shengli Oilfield,Dongying 257000,China |
|
|
Abstract Ray tracing-based first-arrival traveltime tomography is widely used to construct near-surface velocity models to achieve the static correction of seismic data from complex near surface.However,this method necessitates the calculation of ray paths for first-arrival traveltimes and the iterative updating of initial velocity models.As a result,significant computational time is required when applying this method to measured 3D high-density seismic data.To address this issue,this study proposed a method for quickly building 3D near-surface velocity models utilizing diving wave traveltimes under rugged surface.Specifically,based on the ray and traveltime equations of diving waves corresponding to velocities subjected to lateral and vertical changes under rugged surface,the velocity distribution from the observation surface downward was determined using common offset gathers.The proposed method eliminates the need for ray tracing and iterative updates of initial velocity models,offering high modeling efficiency.Tests based on data from theoretical models verified the effectiveness of the proposed method.When applied to measured 3D seismic data,the proposed method yielded static correction results comparable to those obtained using the Fresnel-volume first-arrival traveltime tomography while significantly improving computational efficiency.
|
Received: 18 September 2024
Published: 22 April 2025
|
|
|
|
|
|
The ray path of the diving wave propagating in heterogeneous media with undulant topography
|
|
Velocity models a—the theoretical velocity model;c—models constructed by the first-arrival traveltime inversion based on horizontal observation surface;e—Fresnel first-arrival traveltime tomography model;g—model invesed by the proposed method;b、d、f、h—the velocity slices at y=1 km corresponding to the 3D velocity models shown at the left column
|
|
The slices of velocity relative errors at y=1 km a—the relative errors of the velocity model inverted by the first-arrival traveltime inversion based on horizontal observation surface;b—the relative errors of the velocity model inverted by the Fresnel first-arrival traveltime tomography;c—the relative errors of the velocity model inverted by the proposed method
|
|
The topography of the exploration area
|
|
Observation geometry
|
|
Inverted velocity models and slices at x=3 km a—the models constructed by the first-arrival traveltime inversion based on horizontal observation surface;c—Fresnel first-arrival traveltime tomography model;e—the model inverted by the proposed method;b、d、f—the slices at x=3 km corresponding to the models at the left column respectively
|
|
The 1 200th common shot gather a—the original seismic gather;b—the gather with elevation correction;c—the gather with long-wavelength static correction using the model inverted by the first-arrival traveltime inversion based on horizontal observation surface;d—the gather with long-wavelength static correction using the model inverted by the Fresnel first-arrival traveltime tomography;e—the gather with long-wavelength static correction using the model inverted by the proposed method
|
|
The stacked profiles of the 75th CMP line a—the stack profile without static correction;b—the stack profile with elevation correction;c—the stack profiles with long-wavelength static correction using the model inverted by the first-arrival traveltime inversion based on horizontal observation surface;d—the stack profiles with long-wavelength static correction using the model inverted by the Fresnel first-arrival traveltime tomography;e—the stack profiles with long-wavelength static correction using the model inverted by the proposed method
|
[1] |
张林, 杨勤勇, 张兵, 等. 复杂近地表初至波层析反演静校正技术研究[J]. 地球物理学进展, 2017, 32(2):816-821.
|
[1] |
Zhang L, Yang Q Y, Zhang B, et al. Tomography inversion by first breaks in areas with complex near surface[J]. Progress in Geophysics, 2017, 32(2):816-821.
|
[2] |
于豪. 折射波静校正与层析静校正技术适用性分析[J]. 地球物理学进展, 2012, 27(6):2577-2584.
|
[2] |
Yu H. Applicability analysis of refraction static correction and tomographic inversion static correction[J]. Progress in Geophysics, 2012, 27(6):2577-2584.
|
[3] |
Cox M. Static corrections for seismic reflection surveys[M]. Tulsa: Society of Exploration Geophysicists, 1999.
|
[4] |
周衍, 饶莹. 黄土塬覆盖区的层析反演静校正方法研究[J]. 地球物理学报, 2019, 62(11):4393-4400.
|
[4] |
Zhou Y, Rao Y. Tomographic static corrections in loess plateaus[J]. Chinese Journal of Geophysics, 2019, 62(11):4393-4400.
|
[5] |
Shi T K, Zhang J Z, Huang Z L, et al. A layer-stripping method for 3D near-surface velocity model building using seismic first-arrival times[J]. Journal of Earth Science, 2015, 26(4):502-507.
|
[6] |
林伯香, 孙晶梅, 徐颖, 等. 几种常用静校正方法的讨论[J]. 石油物探, 2006, 45(4):367-372,5.
|
[6] |
Lin B X, Sun J M, Xu Y, et al. Static correction approaches being frequently applied[J]. Geophysical Prospecting for Petroleum, 2006, 45(4):367-372,5.
|
[7] |
黄明忠, 冯泽元, 周大同. 炮检域迭代直接静校正方法及应用效果[J]. 勘探地球物理进展, 2008, 31(2):122-128,87.
|
[7] |
Huang M Z, Feng Z Y, Zhou D T. Direct statics iterated in shot and receiver domains and its application[J]. Progress in Exploration Geophysics, 2008, 31(2):122-128,87.
|
[8] |
肖永新, 杨海申, 崔士天, 等. 一项高效的折射静校正技术[J]. 石油地球物理勘探, 2019, 54(4):768-774,722.
|
[8] |
Xiao Y X, Yang H S, Cui S T, et al. An efficient refraction statics method for massive seismic data[J]. Oil Geophysical Prospecting, 2019, 54(4):768-774,722.
|
[9] |
王克斌, 赵灵芝, 张旭民. 折射静校正在苏里格气田三维处理中的应用[J]. 石油物探, 2003, 42(2):248-251.
|
[9] |
Wang K B, Zhao L Z, Zhang X M. Application of refraction statics in 3-D data processing in Sulige gas field[J]. Geophysical Prospecting for Petroleum, 2003, 42(2):248-251.
|
[10] |
方勇, 罗文山, 姜翠苹, 等. 库车山地地震资料层析静校正方法的应用[J]. 石油地球物理勘探, 2017, 52(S1):23-27.
|
[10] |
Fang Y, Luo W S, Jiang C P, et al. Application of tomographic static correction method for Kuqa mountain seismic data[J]. Oil Geophysical Prospecting, 2017, 52(S1):23-27.
|
[11] |
郝鹏亮. 基于深度学习的复杂地表初至拾取及层析静校正应用研究[D]. 北京: 中国石油大学(北京), 2023.
|
[11] |
Hao P L. Research on application of first break picking and tomographic static correction of complex surface based on deep learning[D]. Beijing: China University of Petroleum (Beijing), 2023.
|
[12] |
Palmer D. The generalized reciprocal method:An integrated approach to shallow refraction seismology[J]. Exploration Geophysics, 1990, 21(1/2):33-44.
|
[13] |
井西利, 杨长春, 李幼铭, 等. 地震静校正全局最优化问题的求解[J]. 地球物理学报, 2002, 45(5):707-713.
|
[13] |
Jing X L, Yang C C, Li Y M, et al. A global optimized algorithm for seismic residual statics corrections[J]. Chinese Journal of Geophysics, 2002, 45(5):707-713.
|
[14] |
Zhang J Z, Shi T K, Zhao Y S, et al. Static corrections in mountainous areas using Fresnel-wavepath tomography[J]. Journal of Applied Geophysics, 2014, 111:242-249.
|
[15] |
Gibson B S. Nonlinear least-squares inversion of traveltime data for a linear velocity-depth relationship[J]. Geophysics, 1979, 44(2):185.
|
[16] |
Yang H C, Zhang J Z, Ren K, et al. First-arrival traveltime inversion of seismic diving waves observed on undulant surface[J]. Geophysical Journal International, 2021, 225(2):1020-1031.
|
[17] |
Jin C K, Zhang J Z. Stereotomography of seismic data acquired on undulant topography[J]. Geophysics, 2018, 83(4):U35-U41.
|
[1] |
ZHANG Xing-Yan, ZENG Wei-Hui, LIU Jin-Peng, ZHANG Li-Xia, DU Guang-Hui, WANG Fa-Kun. Key technologies for processing marine seismic data from sparker sources and their applications[J]. Geophysical and Geochemical Exploration, 2025, 49(1): 166-176. |
[2] |
DUAN Ying, ZHANG Gao-Cheng, TAN Ya-Li. Seismic imaging of high-steep structural zone in Biyang Depression[J]. Geophysical and Geochemical Exploration, 2021, 45(4): 981-989. |
|
|
|
|