The use of ocean bottom seismometers provides an effective means of studying the process and the dynamic of cold seeps by continuously recording micro-events produced by sub-seafloor fluid migration.We deployed a four...The use of ocean bottom seismometers provides an effective means of studying the process and the dynamic of cold seeps by continuously recording micro-events produced by sub-seafloor fluid migration.We deployed a four-component Ocean Bottom Seismometer(OBS)at an active site of the Haima cold seep from 6 November to 19 November in 2021.Here,we present the results of this short-term OBS monitoring.We first examine the OBS record manually to distinguish(by their distinctive seismographic signatures)four types of events:shipping noises,vibrations from our remotely operated vehicle(ROV)operations,local earthquakes,and short duration events(SDEs).Only the SDEs are further discussed in this work.Such SDEs are similar to those observed in other sea areas and are interpreted to be correlated with sub-seafloor fluid migration.In the OBS data collected during the 14-day monitoring period.We identify five SDEs.Compared to the SDE occurrence rate observed in other cold seep regions,five events is rather low,from which it could be inferred that fluid migration,and subsequent gas seepage,is not very active at the Haima site.This conclusion agrees with multi-beam and chemical observations at that site.Our observations thus provide further constraint on the seepage activity in this location.This is the first time that cold seep-related SDEs have been identified in the South China Sea,expanding the list of sea areas where SDEs are now linked to cold seep fluid migration.展开更多
The waveform inversion method is applied-- using synthetic ocean-bottom seismometer （OBS） data--to study oceanic crust structure. A niching genetic algorithm （NGA） is used to implement the inversion for the thickn...The waveform inversion method is applied-- using synthetic ocean-bottom seismometer （OBS） data--to study oceanic crust structure. A niching genetic algorithm （NGA） is used to implement the inversion for the thickness and P-wave velocity of each layer, and to update the model by minimizing the objective function, which consists of the misfit and cross-correlation of observed and synthetic waveforms. The influence of specific NGA method parameters is discussed, and suitable values are presented. The NGA method works well for various observation systems, such as those with irregular and sparse distribu- tion of receivers as well as single receiver systems. A strategy is proposed to accelerate the convergence rate by a factor of five with no increase in computational complex- ity; this is achieved using a first inversion with several generations to impose a restriction on the preset range of each parameter and then conducting a second inversion with the new range. Despite the successes of this method, its usage is limited. A shallow water layer is not favored because the direct wave in water will suppress the useful reflection signals from the crust. A more precise calculation of the air-gun source signal should be considered in order to better simulate waveforms generated in realistic situa- tions; further studies are required to investigate this issue.展开更多
Studies of converted S-wave data recorded on the ocean bottom seismometer(OBS)allow for the estimation of crustal S-wave velocity,from which is further derived the Vp/Vs ratio to constrain the crustal lithology and ge...Studies of converted S-wave data recorded on the ocean bottom seismometer(OBS)allow for the estimation of crustal S-wave velocity,from which is further derived the Vp/Vs ratio to constrain the crustal lithology and geophysical properties.Constructing a precise S-wave velocity model is important for deep structural research,and inversion of converted S-waves provides a potential solution.However,the inversion of the converted S-wave remains a weakness because of the complexity of the seismic ray path and the inconsistent conversion interface.In this study,we introduced two travel time correction methods for the S-wave velocity inversion and imaged different S-wave velocity structures in accordance with the corresponding corrected S-wave phases using seismic data of profile EW6 in the northeastern South China Sea(SCS).The two inversion models show a similar trend in velocities,and the velocity difference is<0.15 km/s(mostly in the range of 0–0.1 km/s),indicating the accuracy of the two travel time correction methods and the reliability of the inversion results.According to simulations of seismic ray tracing based on different models,the velocity of sediments is the primary influencing factor in ray tracing for S-wave phases.If the sedimentary layer has high velocities,the near offset crustal S-wave refractions cannot be traced.In contrast,the ray tracing of Moho S-wave reflections was not significantly impacted by the velocity of the sediments.The two travel time correction methods have their own advantages,and the application of different approaches is based on additional requirements.These works provide an important reference for future improvements in converted S-wave research.展开更多
The acquisition of seabed physical parameters is one of the focuses of marine acoustic researches.However,the activesource ocean bottom seismometer(OBS)detection method in the marine geophysical research is rarely use...The acquisition of seabed physical parameters is one of the focuses of marine acoustic researches.However,the activesource ocean bottom seismometer(OBS)detection method in the marine geophysical research is rarely used to acquire seabed physical parameters,and less work is performed in the Arctic.In this study,two active-source OBS data collected from the 9th and 11th Chinese National Arctic Research Expedition(CHINARE)are selected to obtain the physical parameters of seabed sediments.Two kinds of energy spark are used as the active sources,while the cost function inversion method is used based on the arrival time difference between the reflected and direct waves.The thickness and sound velocity of the sediment layers are obtained by inversion,and the empirical formula is used to calculate the physical parameters of the seabed sediment,which are compared with the measured results.The cost function inversion method based on the time difference of arrival of the reflected and direct waves is tested to be effective and feasible in the inversion of seabed parameters from active-source OBS data.The method is further applied to obtain the physical parameters of Chukchi seabed sediments,which provides the idea and reference for the application of marine geophysical activesource OBS detection technology in the inversion of polar seabed physical parameters.展开更多
Three-component Ocean Bottom Seismometers, portable land stations and marine air gun seismic sources were used to carry out an onshore-offshore deep seismic profile in northeastern South China Sea. This profile, orien...Three-component Ocean Bottom Seismometers, portable land stations and marine air gun seismic sources were used to carry out an onshore-offshore deep seismic profile in northeastern South China Sea. This profile, orientated in NNW-SSE, was as long as 500 km and perpendicular to the strike of regional tectonics. The offshore data were processed in Taiwan Ocean University using a number of available software and the onshore data were analyzed in South China Sea Institute of Oceanology by new-written programs and public software. Preliminary results show that the seismic data are in good quality and contain rich information of deep structure. Seismic phases, e.g. Pg, PmP and Pn, are identified in the offset range 5～220 kin, which will provide an important dataset for the deep crustal structure and oil-gas basin evolution studies of this region.展开更多
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.展开更多
The intense deformation zone in the central Indian Ocean, south of Indian continent is one of the most complex regions in terms of its structure and geodynamics. The deformation zone has been studied and debated in 19...The intense deformation zone in the central Indian Ocean, south of Indian continent is one of the most complex regions in terms of its structure and geodynamics. The deformation zone has been studied and debated in 1990s for its genesis. It was argued that deformation is mainly confined to sedimentary and oceanic crustal layers, while the large wave length geoidal anomalies, on which the deformation region lies, called for deeper sources. The inter connection between deeper and the shallower sources is found missing. The current study focuses on the complexities of this region by analyzing OBS （ocean bottom seismometer） data. The data acquired by five OBS systems along a 300 km long south-north profile in the CIOB （central Indian Ocean basin） have been modeled and the crustal and sub-crustal structure has been determined using 2-D tomographic inversion. Four subsurface layers are identified representing the sediment column, upper crustal layer, lower crustal layer and a sub-crustal layer （upper mantle layer）. A considerable variation in thickness as well as velocity at all interfaces from sedimentary column to upper mantle is observed which indicates that the tectonic forces have affected the entire crust and sub-crustal configuration. The sediments are characterized by higher velocities （2.1 kin/s） due to the increased confining pressure. Modeling results indicated that the velocity in upper crust is in the range of 5.7-6.2 km/s and the velocity of the lower crust varies from 7.0-7.6 km/s. The velocity of the sub-crustal layer is in the range of 7.8-8.4 km/s. This high-velocity layer is interpreted as magmatic under-plating with strong lateral variations. The base of the 7.0 km/s layer at 12-15 km depth is interpreted as the Moho.展开更多
自由表面多次波压制是海底地震仪(Ocean Bottom Seismometer,OBS)数据处理和成像中的难点,OBS数据多次波能量强,周期长,严重影响深层一次反射波的处理和成像.不同于常规拖缆观测系统,OBS数据站点一般相隔较远,仅仅利用检波点稀疏的波场...自由表面多次波压制是海底地震仪(Ocean Bottom Seismometer,OBS)数据处理和成像中的难点,OBS数据多次波能量强,周期长,严重影响深层一次反射波的处理和成像.不同于常规拖缆观测系统,OBS数据站点一般相隔较远,仅仅利用检波点稀疏的波场信息难以压制OBS数据中的自由表面多次波.本文采用拖缆数据与OBS数据联合,利用稀疏反演估计(Estimation of Primaries and Multiples by Sparse Inversion,EPSI)方法,研究了OBS数据自由表面多次波压制理论,分析了OBS多次波产生的机理,详细推导了拖缆数据与OBS数据联合预测OBS多次波的EPSI方法基本原理.通过利用拖缆数据的信息,实现了OBS检波点稀疏数据多次波的压制问题.EPSI方法通过稀疏反演直接估计一次反射波,避免了SRME(Surface Related Multiple Elimination)方法中自适应相减对有效信号的损害,保真了一次反射有效信号,理论模拟OBS数据验证了方法的有效性.展开更多
四分量海底地震仪(OBS)是一种新的接收仪器,其获得的地震数据信息量大,整体信噪比较高,频率成分丰富,具有较高的使用价值。OBS用于研究深部地壳结构已取得较好进展,但较少应用到油气勘探和盆地研究中。在南海西北部采用大容量气枪阵列(0...四分量海底地震仪(OBS)是一种新的接收仪器,其获得的地震数据信息量大,整体信噪比较高,频率成分丰富,具有较高的使用价值。OBS用于研究深部地壳结构已取得较好进展,但较少应用到油气勘探和盆地研究中。在南海西北部采用大容量气枪阵列(0.083m3)作为震源、利用OBS记录气枪震源的反射和折射信息,在油气勘探中是一种新的尝试。为正确使用OBS地震资料,必须采用不同于常规地震资料处理的方法,文中开拓性地使用了镜像叠前时间偏移处理方法。镜像叠前时间偏移处理方法是利用海面多次波进行成像的技术,包括OBS资料矢量保真处理、P分量(压力分量)和Z分量(垂直分量)的叠前处理和镜像叠前时间偏移处理等;其中矢量保真处理主要进行检波点位置二次重定位和检波点方向重定位及倾角校正,P、Z分量的叠前处理包括零相位化处理与地震道修改,P、Z矢量合并(P and Z summation)、剩余静校正、镜像叠加和多次波衰减等。采用镜像叠前时间偏移处理方法获得的剖面获得了较好的成像效果。展开更多
基金supported by the Key Research and Development Project of Guangdong Province(Grant:2020B1111510001)supported by the Project of Sanya Yazhou Bay Science and Technology City(Grant No:SCKJ-JYRC-2022-14)the National Natural Science Foundation of China(Grant No:92262304).
文摘The use of ocean bottom seismometers provides an effective means of studying the process and the dynamic of cold seeps by continuously recording micro-events produced by sub-seafloor fluid migration.We deployed a four-component Ocean Bottom Seismometer(OBS)at an active site of the Haima cold seep from 6 November to 19 November in 2021.Here,we present the results of this short-term OBS monitoring.We first examine the OBS record manually to distinguish(by their distinctive seismographic signatures)four types of events:shipping noises,vibrations from our remotely operated vehicle(ROV)operations,local earthquakes,and short duration events(SDEs).Only the SDEs are further discussed in this work.Such SDEs are similar to those observed in other sea areas and are interpreted to be correlated with sub-seafloor fluid migration.In the OBS data collected during the 14-day monitoring period.We identify five SDEs.Compared to the SDE occurrence rate observed in other cold seep regions,five events is rather low,from which it could be inferred that fluid migration,and subsequent gas seepage,is not very active at the Haima site.This conclusion agrees with multi-beam and chemical observations at that site.Our observations thus provide further constraint on the seepage activity in this location.This is the first time that cold seep-related SDEs have been identified in the South China Sea,expanding the list of sea areas where SDEs are now linked to cold seep fluid migration.
基金supported by the National Natural Science Foundation grant No.41174034the Major State Basic Research Development Program of China(973 Program)
文摘The waveform inversion method is applied-- using synthetic ocean-bottom seismometer （OBS） data--to study oceanic crust structure. A niching genetic algorithm （NGA） is used to implement the inversion for the thickness and P-wave velocity of each layer, and to update the model by minimizing the objective function, which consists of the misfit and cross-correlation of observed and synthetic waveforms. The influence of specific NGA method parameters is discussed, and suitable values are presented. The NGA method works well for various observation systems, such as those with irregular and sparse distribu- tion of receivers as well as single receiver systems. A strategy is proposed to accelerate the convergence rate by a factor of five with no increase in computational complex- ity; this is achieved using a first inversion with several generations to impose a restriction on the preset range of each parameter and then conducting a second inversion with the new range. Despite the successes of this method, its usage is limited. A shallow water layer is not favored because the direct wave in water will suppress the useful reflection signals from the crust. A more precise calculation of the air-gun source signal should be considered in order to better simulate waveforms generated in realistic situa- tions; further studies are required to investigate this issue.
基金The National Natural Science Foundation of China under contract Nos 42276062 and 42006071the Seismological Research Foundation for Youths of Guangdong Earthquake Agency under contract No.GDDZY202307+1 种基金the Strategic Priority Research Program of Chinese Academy of Sciences under contract No.XDA22020303the Science and Technology Planning Project of Guangdong Province-Guangdong Collaborative Innovation Center for Earthquake Prevention and Disaster Mitigation Technology under contract No.2018B020207011.
文摘Studies of converted S-wave data recorded on the ocean bottom seismometer(OBS)allow for the estimation of crustal S-wave velocity,from which is further derived the Vp/Vs ratio to constrain the crustal lithology and geophysical properties.Constructing a precise S-wave velocity model is important for deep structural research,and inversion of converted S-waves provides a potential solution.However,the inversion of the converted S-wave remains a weakness because of the complexity of the seismic ray path and the inconsistent conversion interface.In this study,we introduced two travel time correction methods for the S-wave velocity inversion and imaged different S-wave velocity structures in accordance with the corresponding corrected S-wave phases using seismic data of profile EW6 in the northeastern South China Sea(SCS).The two inversion models show a similar trend in velocities,and the velocity difference is<0.15 km/s(mostly in the range of 0–0.1 km/s),indicating the accuracy of the two travel time correction methods and the reliability of the inversion results.According to simulations of seismic ray tracing based on different models,the velocity of sediments is the primary influencing factor in ray tracing for S-wave phases.If the sedimentary layer has high velocities,the near offset crustal S-wave refractions cannot be traced.In contrast,the ray tracing of Moho S-wave reflections was not significantly impacted by the velocity of the sediments.The two travel time correction methods have their own advantages,and the application of different approaches is based on additional requirements.These works provide an important reference for future improvements in converted S-wave research.
基金supported by the National Key R&D Program of China(No.2021YFC2801200)the National Natural Science Foundation of China(No.42076224)the Fundamental Research Funds for the Central Universities(No.201964015)。
文摘The acquisition of seabed physical parameters is one of the focuses of marine acoustic researches.However,the activesource ocean bottom seismometer(OBS)detection method in the marine geophysical research is rarely used to acquire seabed physical parameters,and less work is performed in the Arctic.In this study,two active-source OBS data collected from the 9th and 11th Chinese National Arctic Research Expedition(CHINARE)are selected to obtain the physical parameters of seabed sediments.Two kinds of energy spark are used as the active sources,while the cost function inversion method is used based on the arrival time difference between the reflected and direct waves.The thickness and sound velocity of the sediment layers are obtained by inversion,and the empirical formula is used to calculate the physical parameters of the seabed sediment,which are compared with the measured results.The cost function inversion method based on the time difference of arrival of the reflected and direct waves is tested to be effective and feasible in the inversion of seabed parameters from active-source OBS data.The method is further applied to obtain the physical parameters of Chukchi seabed sediments,which provides the idea and reference for the application of marine geophysical activesource OBS detection technology in the inversion of polar seabed physical parameters.
基金Supported by SCSIO(LYQY200302)the Chinese Ministry of Science and Technology(G2000046701)+2 种基金the Guangdong Department of Science and Technology(2002C32604)the Guangdong Natural Science Foundation(021557)the National Natural Science Foundation of China(4000161958).
文摘Three-component Ocean Bottom Seismometers, portable land stations and marine air gun seismic sources were used to carry out an onshore-offshore deep seismic profile in northeastern South China Sea. This profile, orientated in NNW-SSE, was as long as 500 km and perpendicular to the strike of regional tectonics. The offshore data were processed in Taiwan Ocean University using a number of available software and the onshore data were analyzed in South China Sea Institute of Oceanology by new-written programs and public software. Preliminary results show that the seismic data are in good quality and contain rich information of deep structure. Seismic phases, e.g. Pg, PmP and Pn, are identified in the offset range 5～220 kin, which will provide an important dataset for the deep crustal structure and oil-gas basin evolution studies of this region.
基金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.
文摘The intense deformation zone in the central Indian Ocean, south of Indian continent is one of the most complex regions in terms of its structure and geodynamics. The deformation zone has been studied and debated in 1990s for its genesis. It was argued that deformation is mainly confined to sedimentary and oceanic crustal layers, while the large wave length geoidal anomalies, on which the deformation region lies, called for deeper sources. The inter connection between deeper and the shallower sources is found missing. The current study focuses on the complexities of this region by analyzing OBS （ocean bottom seismometer） data. The data acquired by five OBS systems along a 300 km long south-north profile in the CIOB （central Indian Ocean basin） have been modeled and the crustal and sub-crustal structure has been determined using 2-D tomographic inversion. Four subsurface layers are identified representing the sediment column, upper crustal layer, lower crustal layer and a sub-crustal layer （upper mantle layer）. A considerable variation in thickness as well as velocity at all interfaces from sedimentary column to upper mantle is observed which indicates that the tectonic forces have affected the entire crust and sub-crustal configuration. The sediments are characterized by higher velocities （2.1 kin/s） due to the increased confining pressure. Modeling results indicated that the velocity in upper crust is in the range of 5.7-6.2 km/s and the velocity of the lower crust varies from 7.0-7.6 km/s. The velocity of the sub-crustal layer is in the range of 7.8-8.4 km/s. This high-velocity layer is interpreted as magmatic under-plating with strong lateral variations. The base of the 7.0 km/s layer at 12-15 km depth is interpreted as the Moho.
文摘自由表面多次波压制是海底地震仪(Ocean Bottom Seismometer,OBS)数据处理和成像中的难点,OBS数据多次波能量强,周期长,严重影响深层一次反射波的处理和成像.不同于常规拖缆观测系统,OBS数据站点一般相隔较远,仅仅利用检波点稀疏的波场信息难以压制OBS数据中的自由表面多次波.本文采用拖缆数据与OBS数据联合,利用稀疏反演估计(Estimation of Primaries and Multiples by Sparse Inversion,EPSI)方法,研究了OBS数据自由表面多次波压制理论,分析了OBS多次波产生的机理,详细推导了拖缆数据与OBS数据联合预测OBS多次波的EPSI方法基本原理.通过利用拖缆数据的信息,实现了OBS检波点稀疏数据多次波的压制问题.EPSI方法通过稀疏反演直接估计一次反射波,避免了SRME(Surface Related Multiple Elimination)方法中自适应相减对有效信号的损害,保真了一次反射有效信号,理论模拟OBS数据验证了方法的有效性.
文摘四分量海底地震仪(OBS)是一种新的接收仪器,其获得的地震数据信息量大,整体信噪比较高,频率成分丰富,具有较高的使用价值。OBS用于研究深部地壳结构已取得较好进展,但较少应用到油气勘探和盆地研究中。在南海西北部采用大容量气枪阵列(0.083m3)作为震源、利用OBS记录气枪震源的反射和折射信息,在油气勘探中是一种新的尝试。为正确使用OBS地震资料,必须采用不同于常规地震资料处理的方法,文中开拓性地使用了镜像叠前时间偏移处理方法。镜像叠前时间偏移处理方法是利用海面多次波进行成像的技术,包括OBS资料矢量保真处理、P分量(压力分量)和Z分量(垂直分量)的叠前处理和镜像叠前时间偏移处理等;其中矢量保真处理主要进行检波点位置二次重定位和检波点方向重定位及倾角校正,P、Z分量的叠前处理包括零相位化处理与地震道修改,P、Z矢量合并(P and Z summation)、剩余静校正、镜像叠加和多次波衰减等。采用镜像叠前时间偏移处理方法获得的剖面获得了较好的成像效果。
