First-arrival seismic traveltime tomography(FAST)is a well-established technique to estimate subsurface velocity structures.Although several existing open-source packages are available for first-arrival traveltime tom...First-arrival seismic traveltime tomography(FAST)is a well-established technique to estimate subsurface velocity structures.Although several existing open-source packages are available for first-arrival traveltime tomography,most were written in compiled languages and lack sufficient extendibility for new algorithms and functionalities.In this work,we develop an open-source,selfcontained FAST package based on MATLAB,one of the most popular interpreted scientific programming languages,with a focus on ocean bottom seismometer refraction traveltime tomography.Our package contains a complete traveltime tomography workflow,including ray-tracing-based first-arrival traveltime computation,linearized inversion,quality control,and high-quality visualization.We design the package as a modular toolbox,making it convenient to integrate new algorithms and functionalities as needed.At the current stage,our package is most efficient for performing FAST for two-dimensional ocean bottom seismometer surveys.We demonstrate the efficacy and accuracy of our package by using a synthetic data example based on a modified Marmousi model.展开更多
Subsurface excavation results in the formation of a zone called excavation damaged zone(EDZ)around the tunnel wall.An EDZ is a major concern in the field of high-level radioactive waste disposal because it may act as ...Subsurface excavation results in the formation of a zone called excavation damaged zone(EDZ)around the tunnel wall.An EDZ is a major concern in the field of high-level radioactive waste disposal because it may act as a flow path after the closure of a repository.In this study,first-arrival traveltime tomography was repeatedly conducted on the EDZ at a depth of 350 m in the Horonobe Underground Research Laboratory.However,the acquired data was highly affected by the support structure on the drift wall.For proper visualization of the EDZ,information about the structure was incorporated into the inversion by modifying the model constraint.The synthetic study showed that the approach reproduced the EDZ in the model without the artifacts.The method was applied to field data,and the EDZ around the drift was detected.The inversion was extended to a time-lapse inversion to trace the changes in P-wave velocity in the EDZ.The synthetic study demonstrated that temporal changes in the P-wave velocity distribution could be detected.Data obtained from 12 surveys under open-drift conditions were analyzed by time-lapse inversion.The results indicated that the EDZ did not undergo sealing or evolution at the site for approximately seven years.展开更多
In order to improve the efficiency of 3D near-surface velocity model building, we develop a layer-stripping method using seismic first-arrival times. The velocity model within a Common Mid-Point (CMP) gather is assu...In order to improve the efficiency of 3D near-surface velocity model building, we develop a layer-stripping method using seismic first-arrival times. The velocity model within a Common Mid-Point (CMP) gather is assumed to be stratified into thin layers, and the velocity of each layer var- ies linearly with depth. The thickness and velocity of the top layer are estimated using minimum-offset first-arrival data in a CMP gather. Then the top layer is stripped and the second layer becomes a new top layer. After removing the effect of the top layer from the former first-arrival data, the new first-arrival data are obtained and then used to estimate the parameters of the second layer. In this manner, the velocity model, being regarded as that at a CMP location, is built layer-by-layer from the top to the bottom. A 3D near-surface velocity model is then formed using the velocity models at all CMP locations. The tests on synthetic and observed seismic data show that the layer-stripping method can be used to build good near-surface velocity models for static correction, and its computation speed is approximately hundred times faster than that of grid tomography.展开更多
Estimation of Thomsen's anisotropic parameters is very important for accuratetime-to-depth conversion and depth migration data processing. Compared with othermethods, it is much easier and more reliable to estimate a...Estimation of Thomsen's anisotropic parameters is very important for accuratetime-to-depth conversion and depth migration data processing. Compared with othermethods, it is much easier and more reliable to estimate anisotropic parameters that arerequired for surface seismic depth imaging from vertical seismic profile (VSP) data, becausethe first arrivals of VSP data can be picked with much higher accuracy. In this study, wedeveloped a method for estimating Thomsen's P-wave anisotropic parameters in VTImedia using the first arrivals from walkaway VSP data. Model first-arrival travel times arecalculated on the basis of the near-offset normal moveout correction velocity in VTI mediaand ray tracing using Thomsen's P-wave velocity approximation. Then, the anisotropicparameters 0 and e are determined by minimizing the difference between the calculatedand observed travel times for the near and far offsets. Numerical forward modeling, usingthe proposed method indicates that errors between the estimated and measured anisotropicparameters are small. Using field data from an eight-azimuth walkaway VSP in TarimBasin, we estimated the parameters 0 and e and built an anisotropic depth-velocity modelfor prestack depth migration processing of surface 3D seismic data. The results showimprovement in imaging the carbonate reservoirs and minimizing the depth errors of thegeological targets.展开更多
Transducers that are widely applied in cement bond evaluation tools, such as cement bond logs and variable density logs, cannot radiate acoustic energy directionally because of the characteristics of monopole sources....Transducers that are widely applied in cement bond evaluation tools, such as cement bond logs and variable density logs, cannot radiate acoustic energy directionally because of the characteristics of monopole sources. A phased arc array transmitter, which is a novel transducer that differs from monopole and dipole transducers, is presented in this study. To simulate the acoustic field generated by a phased arc array in a fluid-filled cased borehole with different channelings, a 3D finite-difference time-domain method is adopted. The acoustic field generated by a traditional monopole source is also simulated and compared with the field generated by the phased arc array transmitter. Numerical simulation results show that the phased arc array radiates energy directionally in a narrow angular range in the borehole, thereby compressing the acoustic energy into a narrow range in the casing pipe, the cement, and the formation. We present the analyses of first-arrival waveforms and the amplitudes of casing waves at different azimuthal angles for the two different sources. The results indicate that employing a directional source facilitates azimuthal identification and analysis of possible channeling behind the casing pipe.展开更多
To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, t...To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.展开更多
基金financially supported by the Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) (Grant No. GML2019 ZD0207)supported by the Guangzhou Municipal Science and Technology Bureau (Grant No. 202102021054)
文摘First-arrival seismic traveltime tomography(FAST)is a well-established technique to estimate subsurface velocity structures.Although several existing open-source packages are available for first-arrival traveltime tomography,most were written in compiled languages and lack sufficient extendibility for new algorithms and functionalities.In this work,we develop an open-source,selfcontained FAST package based on MATLAB,one of the most popular interpreted scientific programming languages,with a focus on ocean bottom seismometer refraction traveltime tomography.Our package contains a complete traveltime tomography workflow,including ray-tracing-based first-arrival traveltime computation,linearized inversion,quality control,and high-quality visualization.We design the package as a modular toolbox,making it convenient to integrate new algorithms and functionalities as needed.At the current stage,our package is most efficient for performing FAST for two-dimensional ocean bottom seismometer surveys.We demonstrate the efficacy and accuracy of our package by using a synthetic data example based on a modified Marmousi model.
基金This study was conducted as part of a commissioned project entitled“The project for validating sealing of the geological repository(2020 and 2021 FY)”funded by the Ministry of Economy,Trade,and Industry of Japan.
文摘Subsurface excavation results in the formation of a zone called excavation damaged zone(EDZ)around the tunnel wall.An EDZ is a major concern in the field of high-level radioactive waste disposal because it may act as a flow path after the closure of a repository.In this study,first-arrival traveltime tomography was repeatedly conducted on the EDZ at a depth of 350 m in the Horonobe Underground Research Laboratory.However,the acquired data was highly affected by the support structure on the drift wall.For proper visualization of the EDZ,information about the structure was incorporated into the inversion by modifying the model constraint.The synthetic study showed that the approach reproduced the EDZ in the model without the artifacts.The method was applied to field data,and the EDZ around the drift was detected.The inversion was extended to a time-lapse inversion to trace the changes in P-wave velocity in the EDZ.The synthetic study demonstrated that temporal changes in the P-wave velocity distribution could be detected.Data obtained from 12 surveys under open-drift conditions were analyzed by time-lapse inversion.The results indicated that the EDZ did not undergo sealing or evolution at the site for approximately seven years.
基金supported by the National Natural Science Foundation of China(Nos.41230318,41074077)the Specialized Research Fund for the Doctoral Program of Higher Education(No.20130132110023)the Fundamental Research Funds for the Central Universities of China(No.201413004)
文摘In order to improve the efficiency of 3D near-surface velocity model building, we develop a layer-stripping method using seismic first-arrival times. The velocity model within a Common Mid-Point (CMP) gather is assumed to be stratified into thin layers, and the velocity of each layer var- ies linearly with depth. The thickness and velocity of the top layer are estimated using minimum-offset first-arrival data in a CMP gather. Then the top layer is stripped and the second layer becomes a new top layer. After removing the effect of the top layer from the former first-arrival data, the new first-arrival data are obtained and then used to estimate the parameters of the second layer. In this manner, the velocity model, being regarded as that at a CMP location, is built layer-by-layer from the top to the bottom. A 3D near-surface velocity model is then formed using the velocity models at all CMP locations. The tests on synthetic and observed seismic data show that the layer-stripping method can be used to build good near-surface velocity models for static correction, and its computation speed is approximately hundred times faster than that of grid tomography.
基金supported by the National Science and Technology Major Project of China(No.2011ZX05046)
文摘Estimation of Thomsen's anisotropic parameters is very important for accuratetime-to-depth conversion and depth migration data processing. Compared with othermethods, it is much easier and more reliable to estimate anisotropic parameters that arerequired for surface seismic depth imaging from vertical seismic profile (VSP) data, becausethe first arrivals of VSP data can be picked with much higher accuracy. In this study, wedeveloped a method for estimating Thomsen's P-wave anisotropic parameters in VTImedia using the first arrivals from walkaway VSP data. Model first-arrival travel times arecalculated on the basis of the near-offset normal moveout correction velocity in VTI mediaand ray tracing using Thomsen's P-wave velocity approximation. Then, the anisotropicparameters 0 and e are determined by minimizing the difference between the calculatedand observed travel times for the near and far offsets. Numerical forward modeling, usingthe proposed method indicates that errors between the estimated and measured anisotropicparameters are small. Using field data from an eight-azimuth walkaway VSP in TarimBasin, we estimated the parameters 0 and e and built an anisotropic depth-velocity modelfor prestack depth migration processing of surface 3D seismic data. The results showimprovement in imaging the carbonate reservoirs and minimizing the depth errors of thegeological targets.
基金supported by the National Natural ScienceFoundation of China (Grant Nos. 11204380, 11374371, 11134011 and 61102102)National Science and Technology Major Project (Grant No. 2011ZX05020-009)PetroChina Innovation Foundation (2013D-5006-0304)
文摘Transducers that are widely applied in cement bond evaluation tools, such as cement bond logs and variable density logs, cannot radiate acoustic energy directionally because of the characteristics of monopole sources. A phased arc array transmitter, which is a novel transducer that differs from monopole and dipole transducers, is presented in this study. To simulate the acoustic field generated by a phased arc array in a fluid-filled cased borehole with different channelings, a 3D finite-difference time-domain method is adopted. The acoustic field generated by a traditional monopole source is also simulated and compared with the field generated by the phased arc array transmitter. Numerical simulation results show that the phased arc array radiates energy directionally in a narrow angular range in the borehole, thereby compressing the acoustic energy into a narrow range in the casing pipe, the cement, and the formation. We present the analyses of first-arrival waveforms and the amplitudes of casing waves at different azimuthal angles for the two different sources. The results indicate that employing a directional source facilitates azimuthal identification and analysis of possible channeling behind the casing pipe.
基金National Natural Science Foundation of China (49894190-024) and Geophysical Prospecting Key Laboratory Foundation of China National Petroleum Corporation.
文摘To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.