To solve the problems in restoring sedimentary facies and predicting reservoirs in loose gas-bearing sediment,based on seismic sedimentologic analysis of the first 9-component S-wave 3D seismic dataset of China,a four...To solve the problems in restoring sedimentary facies and predicting reservoirs in loose gas-bearing sediment,based on seismic sedimentologic analysis of the first 9-component S-wave 3D seismic dataset of China,a fourth-order isochronous stratigraphic framework was set up and then sedimentary facies and reservoirs in the Pleistocene Qigequan Formation in Taidong area of Qaidam Basin were studied by seismic geomorphology and seismic lithology.The study method and thought are as following.Firstly,techniques of phase rotation,frequency decomposition and fusion,and stratal slicing were applied to the 9-component S-wave seismic data to restore sedimentary facies of major marker beds based on sedimentary models reflected by satellite images.Then,techniques of seismic attribute extraction,principal component analysis,and random fitting were applied to calculate the reservoir thickness and physical parameters of a key sandbody,and the results are satisfactory and confirmed by blind testing wells.Study results reveal that the dominant sedimentary facies in the Qigequan Formation within the study area are delta front and shallow lake.The RGB fused slices indicate that there are two cycles with three sets of underwater distributary channel systems in one period.Among them,sandstones in the distributary channels of middle-low Qigequan Formation are thick and broad with superior physical properties,which are favorable reservoirs.The reservoir permeability is also affected by diagenesis.Distributary channel sandstone reservoirs extend further to the west of Sebei-1 gas field,which provides a basis to expand exploration to the western peripheral area.展开更多
Accurate prediction of formation pore pressure is essential to predict fluid flow and manage hydrocarbon production in petroleum engineering.Recent deep learning technique has been receiving more interest due to the g...Accurate prediction of formation pore pressure is essential to predict fluid flow and manage hydrocarbon production in petroleum engineering.Recent deep learning technique has been receiving more interest due to the great potential to deal with pore pressure prediction.However,most of the traditional deep learning models are less efficient to address generalization problems.To fill this technical gap,in this work,we developed a new adaptive physics-informed deep learning model with high generalization capability to predict pore pressure values directly from seismic data.Specifically,the new model,named CGP-NN,consists of a novel parametric features extraction approach(1DCPP),a stacked multilayer gated recurrent model(multilayer GRU),and an adaptive physics-informed loss function.Through machine training,the developed model can automatically select the optimal physical model to constrain the results for each pore pressure prediction.The CGP-NN model has the best generalization when the physicsrelated metricλ=0.5.A hybrid approach combining Eaton and Bowers methods is also proposed to build machine-learnable labels for solving the problem of few labels.To validate the developed model and methodology,a case study on a complex reservoir in Tarim Basin was further performed to demonstrate the high accuracy on the pore pressure prediction of new wells along with the strong generalization ability.The adaptive physics-informed deep learning approach presented here has potential application in the prediction of pore pressures coupled with multiple genesis mechanisms using seismic data.展开更多
The South Yellow Sea basin is filled with Mesozoic-Cenozoic continental sediments overlying pre-Palaeozoic and Mesozoic-Palaeozoic marine sediments.Conventional multi-channel seismic data cannot describe the velocity ...The South Yellow Sea basin is filled with Mesozoic-Cenozoic continental sediments overlying pre-Palaeozoic and Mesozoic-Palaeozoic marine sediments.Conventional multi-channel seismic data cannot describe the velocity structure of the marine residual basin in detail,leading to the lack of a deeper understanding of the distribution and lithology owing to strong energy shielding on the top interface of marine sediments.In this study,we present seismic tomography data from ocean bottom seismographs that describe the NEE-trending velocity distributions of the basin.The results indicate that strong velocity variations occur at shallow crustal levels.Horizontal velocity bodies show good correlation with surface geological features,and multi-layer features exist in the vertical velocity framework(depth:0–10 km).The analyses of the velocity model,gravity data,magnetic data,multichannel seismic profiles,and drilling data showed that high-velocity anomalies(>6.5 km/s)of small(thickness:1–2 km)and large(thickness:>5 km)scales were caused by igneous complexes in the multi-layer structure,which were active during the Palaeogene.Possible locations of good Mesozoic and Palaeozoic marine strata are limited to the Central Uplift and the western part of the Northern Depression along the wide-angle ocean bottom seismograph array.Following the Indosinian movement,a strong compression existed in the Northern Depression during the extensional phase that caused the formation of folds in the middle of the survey line.This study is useful for reconstructing the regional tectonic evolution and delineating the distribution of the marine residual basin in the South Yellow Sea basin.展开更多
The Solomon Sea Basin is a Cenozoic back-arc spreading basin within the convergence system of the Pacific and Indo-Australian plates.Against the background of subduction polarity reversal,the current Solomon Sea Basin...The Solomon Sea Basin is a Cenozoic back-arc spreading basin within the convergence system of the Pacific and Indo-Australian plates.Against the background of subduction polarity reversal,the current Solomon Sea Basin gradually formed a rhombic morphology with the subduction of the basin along the New Britain Trench and the Trobriand Trough.By analyzing the vertical gravity gradient,natural earthquake and seismic reflection data,this study determines the structural characteristics of the Solomon Sea Basin.It was found that the tectonics of the basin are characterized by the original expansion structure within the central part in addition to the structure induced by the latest subduction along the basin margin.The original spreading structure of the basin presented an east–west linear graben and horst controlled by normal faults during the basin expansion period.As a result of the subduction and slab-pull of the Solomon Sea Basin,extensional structure belts parallel to the New Britain Trench formed along the basin margin.展开更多
At present,the acquisition of seismic data is developing toward high-precision and high-density methods.However,complex natural environments and cultural factors in many exploration areas cause difficulties in achievi...At present,the acquisition of seismic data is developing toward high-precision and high-density methods.However,complex natural environments and cultural factors in many exploration areas cause difficulties in achieving uniform and intensive acquisition,which makes complete seismic data collection impossible.Therefore,data reconstruction is required in the processing link to ensure imaging accuracy.Deep learning,as a new field in rapid development,presents clear advantages in feature extraction and modeling.In this study,the convolutional neural network deep learning algorithm is applied to seismic data reconstruction.Based on the convolutional neural network algorithm and combined with the characteristics of seismic data acquisition,two training strategies of supervised and unsupervised learning are designed to reconstruct sparse acquisition seismic records.First,a supervised learning strategy is proposed for labeled data,wherein the complete seismic data are segmented as the input of the training set and are randomly sampled before each training,thereby increasing the number of samples and the richness of features.Second,an unsupervised learning strategy based on large samples is proposed for unlabeled data,and the rolling segmentation method is used to update(pseudo)labels and training parameters in the training process.Through the reconstruction test of simulated and actual data,the deep learning algorithm based on a convolutional neural network shows better reconstruction quality and higher accuracy than compressed sensing based on Curvelet transform.展开更多
The Kuiyang-ST2000 deep-towed high-resolution multichannel seismic system was designed by the First Institute of Oceanography,Ministry of Natural Resources(FIO,MNR).The system is mainly composed of a plasma spark sour...The Kuiyang-ST2000 deep-towed high-resolution multichannel seismic system was designed by the First Institute of Oceanography,Ministry of Natural Resources(FIO,MNR).The system is mainly composed of a plasma spark source(source level:216 dB,main frequency:750 Hz,frequency bandwidth:150-1200 Hz)and a towed hydrophone streamer with 48 channels.Because the source and the towed hydrophone streamer are constantly moving according to the towing configuration,the accurate positioning of the towing hydrophone array and the moveout correction of deep-towed multichannel seismic data processing before imaging are challenging.Initially,according to the characteristics of the system and the towing streamer shape in deep water,travel-time positioning method was used to construct the hydrophone streamer shape,and the results were corrected by using the polynomial curve fitting method.Then,a new data-processing workflow for Kuiyang-ST2000 system data was introduced,mainly including float datum setting,residual static correction,phase-based moveout correction,which allows the imaging algorithms of conventional marine seismic data processing to extend to deep-towed seismic data.We successfully applied the Kuiyang-ST2000 system and methodology of data processing to a gas hydrate survey of the Qiongdongnan and Shenhu areas in the South China Sea,and the results show that the profile has very high vertical and lateral resolutions(0.5 m and 8 m,respectively),which can provide full and accurate details of gas hydrate-related and geohazard sedimentary and structural features in the South China Sea.展开更多
Paleostress plays a significant role in controlling the formation, accumulation, and distribution of reservoirs, and this could be an important factor in controlling the production of hydrocarbons from the unconventio...Paleostress plays a significant role in controlling the formation, accumulation, and distribution of reservoirs, and this could be an important factor in controlling the production of hydrocarbons from the unconventional reservoirs. In this study, we will use 3D seismic reflection data to perform the slip-tendency-based stress inversion to determine the stress field in the basement of the northern slope area in the Bongor Basin. The dataset for this technique is easily available in the oil and gas companies. The stress inversion results from the basement of the northern slope area of Bongor basin show that the maximum principal stress axis (σ1) is oriented vertically, the intermediate principal stress axis (σ2) is oriented N18° and the minimum principal stress axis (σ3) is oriented N105°, and σ2/σ1 = 0.60 and σ3/σ1 = 0.29. The findings of this paper provide significant information to understand the fault reactivation at the critical stage of hydrocarbon accumulation and the regional tectonic evolution.展开更多
The Pearl River Estuary(PRE) is located at the onshore-offshore transition zone between South China and South China Sea Basin, and it is of great significant value in discussing tectonic relationships between South Ch...The Pearl River Estuary(PRE) is located at the onshore-offshore transition zone between South China and South China Sea Basin, and it is of great significant value in discussing tectonic relationships between South China block and South China Sea block and seismic activities along the offshore active faults in PRE. However, the researches on geometric characteristics of offshore faults in this area are extremely lacking. To investigate the offshore fault distribution and their geometric features in the PRE in greater detail, we acquired thirteen seismic reflection profiles in 2015. Combining the analysis of the seismic reflection and free-air gravity anomaly data, this paper revealed the location, continuity, and geometry of the littoral fault zone and other offshore faults in PRE. The littoral fault zone is composed of the major Dangan Islands fault and several parallel, high-angle, normal faults, which mainly trend northeast to northeast-to-east and dip to the southeast with large displacements. The fault zone is divided into three different segments by the northwest-trending faults. Moreover, the basement depth around Dangan Islands is very shallow, while it suddenly increases along the islands westward and southward. These has resulted in the islands and neighboring areas becoming the places where the stress accumulates easily. The seismogenic pattern of this area is closely related to the comprehensive effect of intersecting faults together with the low velocity layer.展开更多
The field seismic data is disturbed by the interferential information, which has low signal to noise ratio (SNR). That is disadvantage for seismic data interpretation. So it is important to remove the noise of seismic...The field seismic data is disturbed by the interferential information, which has low signal to noise ratio (SNR). That is disadvantage for seismic data interpretation. So it is important to remove the noise of seismic data. Independent component analysis (ICA) can remove most of the noise interference. However, ICA has some defects in noise reduction, because it needs some conditions that seismic data is independent reciprocally for denoising. To solve these defects, this paper proposes an improved ICA algorithm to noise reduction. Through simulation experiments, it can be obtained that the best decomposition levels of the new algorithm is 3. At last, the proposed improved ICA is applied to deal with the actual seismic data. The results show that it can effectively eliminate most of seismic noise such as random noise, linear interference, surface waves, and so on. The improved ICA is not only easy to denoising, but also has excellent mathematical theoretical properties.展开更多
A new method is introduced to suppress the noise in seismic data processing. Based on the subtle difference in shape between the noise and the actual signal, we introduce morphologic filtering into seismic data proces...A new method is introduced to suppress the noise in seismic data processing. Based on the subtle difference in shape between the noise and the actual signal, we introduce morphologic filtering into seismic data processing. From the shape and the S/N we can see that the effect of morphologic filtering is superior to other methods like id-value filtering, neighbor average filtering, etc. The SNR of the signal after morphological filtering is comparatively great. In addition, the precision of the seismic data after morphological filtering is high. The characteristics of the actual signal, such as frequency and amplitude, are preserved. We give an example of the real seismic data processing using morphological filtering, in which the actual signal is retained, while the random high intensity noise was removed.展开更多
Comprehensive inversion of logging and seismic data presented in this paper is a method to improve seismic data resolution. It involves using ample high-frequency information and complete low-frequency information of ...Comprehensive inversion of logging and seismic data presented in this paper is a method to improve seismic data resolution. It involves using ample high-frequency information and complete low-frequency information of known logging to make up for the lack of limited bandwidth of practical seismic recording, obtaining an approximate reflection coefficient sequence (or wave impedance) of high resolution by iterative inversion and providing more reliable seismic evidence for further lithologic interpretation and lateral tracking, correlation and prediction of thin reservoir. The comprehensive inversion can be realized in the following steps: (1) to establish an initial model of higher resolution; (2) to obtain wavelets, and (3) to constrain iterative inversion. The key to this inversion lies in building an initial model. It is assumed from our experience that when the initial model is properly given, iterative inversion can be quickly converged to the ideal result.展开更多
A novel time-domain identification technique is developed for the seismic response analysis of soil-structure interaction.A two-degree-of-freedom (2DOF) model with eight lumped parameters is adopted to model the frequ...A novel time-domain identification technique is developed for the seismic response analysis of soil-structure interaction.A two-degree-of-freedom (2DOF) model with eight lumped parameters is adopted to model the frequency- dependent behavior of soils.For layered soil,the equivalent eight parameters of the 2DOF model are identified by the extended Kalman filter (EKF) method using recorded seismic data.The polynomial approximations for derivation of state estimators are applied in the EKF procedure.A realistic identification example is given for the layered-soil of a building site in Anchorage,Alaska in the United States.Results of the example demonstrate the feasibility and practicality of the proposed identification technique.The 2DOF soil model and the identification technique can be used for nonlinear response analysis of soil-structure interaction in the time-domain for layered or complex soil conditions.The identified parameters can be stored in a database tor use in other similar soil conditions,lfa universal database that covers information related to most soil conditions is developed in the thture,engineers could conveniently perform time history analyses of soil-structural interaction.展开更多
The Atlas region in northwest Africa is characterized by the Quaternary volcanism and elevated topography with past complex tectonic mitigation between the African and European plates.Geodynamics of this atypical regi...The Atlas region in northwest Africa is characterized by the Quaternary volcanism and elevated topography with past complex tectonic mitigation between the African and European plates.Geodynamics of this atypical region has left indubitably imprints in crustal architectonics,mainly regarding the crustal thickness as well as the crustal density structure.The knowledge of crustal thickness variations is of a significant interest,since it provides a crucial constraint to geodynamic and geophysical modelling of this region.In this study,we use gravity,topographic,bathymetric and sediment data together with results of seismic surveys to image the Moho topography beneath the Atlas region.The Bouguer gravity anomalies used for a gravimetric Moho recovery are obtained from the free-air gravity anomalies after subtracting the gravitational contributions of topography,bathymetry and sediments.The regional gravimetric Moho inversion constrained on seismic data is carried out by applying a regularized inversion technique based on Gauss-Newton’s formulation of improved Bott’s method,while adopting Earth’s spherical approximation.The numerical result reveals relatively significant Moho depth variations in the Moroccan Atlas,with minima of approximately 24 km along continental margins of the Mediterranean Sea and maxima exceeding 51 km beneath the Rif Cordillera.The Moho depth beneath the West African Craton varies from 32 km in its southern margin to 45 km beneath the Middle Atlas.The Tell Atlas is characterized by the shallow Moho depth of approximately 22 km and further deepening to 42 km towards the northern edge of the Aures Mountains.Our findings indicate a limited tectonic shortening of the High Atlas with the crustal thickness mostly within 36-42 km.Topographic discrepancies between the Rif Cordillera and the Atlas Mountains suggest that the hypothesis of isostatic compensation cannot be fully established.展开更多
On 21 May 2021(UTC),an MW 7.4 earthquake jolted the east Bayan Har block in the Tibetan Plateau.The earthquake received widespread attention as it is the largest event in the Tibetan Plateau and its surroundings since...On 21 May 2021(UTC),an MW 7.4 earthquake jolted the east Bayan Har block in the Tibetan Plateau.The earthquake received widespread attention as it is the largest event in the Tibetan Plateau and its surroundings since the 2008 Wenchuan earthquake,and especially in proximity to the seismic gaps on the east Kunlun fault.Here we use satellite interferometric synthetic aperture radar data and subpixel offset observations along the range directions to characterize the coseismic deformation of the earthquake.Range offset displacements depict clear surface ruptures with a total length of~170 km involving two possible activated fault segments in the earthquake.Coseismic modeling results indicate that the earthquake was dominated by left-lateral strike-slip motions of up to 7 m within the top 12 km of the crust.The well-resolved slip variations are characterized by five major slip patches along strike and 64%of shallow slip deficit,suggesting a young seismogenic structure.Spatial-temporal changes of the postseismic deformation are mapped from early 6-day and 24-day InSAR observations,and are well explained by time-dependent afterslip models.Analysis of Global Navigation Satellite System(GNSS)velocity profiles and strain rates suggests that the eastward extrusion of plateau is diffusely distributed across the east Bayan Har block,but exhibits significant lateral heterogeneities,as evidenced by magnetotelluric observations.The block-wide distributed deformation of the east Bayan Har block along with the significant co-and post-seismic stress loadings from the Madoi earthquake imply high seismic risks along regional faults,especially the Tuosuo Lake and Maqên-Maqu segments of the Kunlun fault that are known as seismic gaps.展开更多
A novel technique for automatic seismic data processing using both integral and local feature of seismograms was presented in this paper. Here, the term integral feature of seismograms refers to feature which may depi...A novel technique for automatic seismic data processing using both integral and local feature of seismograms was presented in this paper. Here, the term integral feature of seismograms refers to feature which may depict the shape of the whole seismograms. However, unlike some previous efforts which completely abandon the DIAL approach, i.e., signal detection, phase identifi- cation, association, and event localization, and seek to use envelope cross-correlation to detect seismic events directly, our technique keeps following the DIAL approach, but in addition to detect signals corresponding to individual seismic phases, it also detects continuous wave-trains and explores their feature for phase-type identification and signal association. More concrete ideas about how to define wave-trains and combine them with various detections, as well as how to measure and utilize their feature in the seismic data processing were expatiated in the paper. This approach has been applied to the routine data processing by us for years, and test results for a 16 days' period using data from the Xinjiang seismic station network were presented. The automatic processing results have fairly low false and missed event rate simultaneously, showing that the new technique has good application prospects for improvement of the automatic seismic data processing.展开更多
This paper is concerned with anisotropic effects on seismic data and signal analysis for transversely isotropic rock media with vertical anisotropy. It is understood that these effects are significant in many practica...This paper is concerned with anisotropic effects on seismic data and signal analysis for transversely isotropic rock media with vertical anisotropy. It is understood that these effects are significant in many practical applications, e.g. earthquake forecasting, materials exploration inside the Earth’s crust, as well as various practical works in oil industry. Under the framework of the most accepted anisotropic media model (i.e. VTI media, transverse isotropy with a vertical axis symmetry), with applications of a set of available anisotropic rock parameters for sandstone and shale, we have performed numerical calculations of the anisotropic effects. We show that for rocks with strong anisotropy, the induced relative depth error can be significantly large. Nevertheless, with an improved understanding of the seismic-signal propagation and proper data processing, the error can be reduced, which in turn may enhance the probability of forecasting accurately the various wave propagations inside the Earth’s crust, e.g. correctly forecasting the incoming earthquakes from the center of the Earth.展开更多
Oil and gas seismic exploration have to adopt irregular seismic acquisition due to the increasingly complex exploration conditions to adapt to complex geological conditions and environments.However,the irregular seism...Oil and gas seismic exploration have to adopt irregular seismic acquisition due to the increasingly complex exploration conditions to adapt to complex geological conditions and environments.However,the irregular seismic acquisition is accompanied by the lack of acquisition data,which requires high-precision regularization.The sparse signal feature in the transform domain in compressed sensing theory is used in this paper to recover the missing signal,involving sparse transform base optimization and threshold modeling.First,this paper analyzes and compares the effects of six sparse transformation bases on the reconstruction accuracy and efficiency of irregular seismic data and establishes the quantitative relationship between sparse transformation and reconstruction accuracy and efficiency.Second,an adaptive threshold modeling method based on sparse coefficient is provided to improve the reconstruction accuracy.Test results show that the method has good adaptability to different seismic data and sparse transform bases.The f-x domain reconstruction method of effective frequency samples is studied to address the problem of low computational efficiency.The parallel computing strategy of curvelet transform combined with OpenMP is further proposed,which substantially improves the computational efficiency under the premise of ensuring the reconstruction accuracy.Finally,the actual acquisition data are used to verify the proposed method.The results indicate that the proposed method strategy can solve the regularization problem of irregular seismic data in production and improve the imaging quality of the target layer economically and efficiently.展开更多
Sandy debris flow deposits are present in Unit I during Miocene of Gas Field A in the Baiyun Depression of the South China Sea. The paucity of well data and the great variability of the sedimentary microfacies make it...Sandy debris flow deposits are present in Unit I during Miocene of Gas Field A in the Baiyun Depression of the South China Sea. The paucity of well data and the great variability of the sedimentary microfacies make it difficult to identify and predict the distribution patterns of the main gas reservoir, and have seriously hindered further exploration and development of the gas field. Therefore, making full use of the available seismic data is extremely important for predicting the spatial distribution of sedimentary microfacies when constructing three-dimensional reservoir models. A suitable reservoir modeling strategy or workflow controlled by sedimentary microfacies and seismic data has been developed. Five types of seismic attributes were selected to correlate with the sand percentage, and the root mean square (RMS) amplitude performed the best. The relation between the RMS amplitude and the sand percentage was used to construct a reservoir sand distribution map. Three types of main sedimentary microfacies were identified: debris channels, fan lobes, and natural levees. Using constraints from the sedimentary microfacies boundaries, a sedimentary microfacies model was constructed using the sequential indicator and assigned value simulation methods. Finally, reservoir models of physical properties for sandy debris flow deposits controlled by sedimentary microfacies and seismic inversion data were established. Property cutoff values were adopted because the sedimentary microfacies and the reservoir properties from well-logging interpretation are intrinsically different. Selection of appropriate reservoir property cutoffs is a key step in reservoir modeling when using simulation methods based on sedimentary microfacies control. When the abnormal data are truncated and the reservoir properties probability distribution fits a normal distribution, microfacies-controlled reservoir property models are more reliable than those obtained from the sequence Gauss simulation method. The cutoffs for effective porosity of the debris channel, fan lobe, and natural levee facies were 0.2, 0.09, and 0.12, respectively; the corresponding average effective porosities were 0.24, 0.13, and 0.15. The proposed modeling method makes full use of seismic attributes and seismic inversion data, and also makes the property data of single-well depositional microfacies more conformable to a normal distribution with geological significance. Thus, the method allows use of more reliable input data when we construct a model of a sandy debris flow.展开更多
In this paper, a new concept called numerical structure of seismic data is introduced and the difference between numerical structure and numerical value of seismic data is explained. Our study shows that the numerical...In this paper, a new concept called numerical structure of seismic data is introduced and the difference between numerical structure and numerical value of seismic data is explained. Our study shows that the numerical seismic structure is closely related to oil and gas-bearing reservoir,so it is very useful for a geologist or a geophysicist to precisely interpret the oil-bearing layers from the seismic data. This technology can be applied to any exploration or production stage. The new method has been tested on a series of exploratory or development wells and proved to be reliable in China. Hydrocarbon-detection with this new method for 39 exploration wells on 25 structures indicates a success ratio of over 80 percent. The new method of hydrocarbon prediction can be applied for: (1) depositional environment of reservoirs with marine facies, delta, or non-marine facies (including fluvial facies, lacustrine facies); (2) sedimentary rocks of reservoirs that are non-marine clastic rocks and carbonate rock; and (3) burial depths range from 300 m to 7000 m, and the minimum thickness of these reservoirs is over 8 m (main frequency is about 50 Hz).展开更多
基金Supported by the CNPC Science and Technology Projects(2022-N/G-47808,2023-N/G-67014)RIPED International Cooperation Project(19HTY5000008).
文摘To solve the problems in restoring sedimentary facies and predicting reservoirs in loose gas-bearing sediment,based on seismic sedimentologic analysis of the first 9-component S-wave 3D seismic dataset of China,a fourth-order isochronous stratigraphic framework was set up and then sedimentary facies and reservoirs in the Pleistocene Qigequan Formation in Taidong area of Qaidam Basin were studied by seismic geomorphology and seismic lithology.The study method and thought are as following.Firstly,techniques of phase rotation,frequency decomposition and fusion,and stratal slicing were applied to the 9-component S-wave seismic data to restore sedimentary facies of major marker beds based on sedimentary models reflected by satellite images.Then,techniques of seismic attribute extraction,principal component analysis,and random fitting were applied to calculate the reservoir thickness and physical parameters of a key sandbody,and the results are satisfactory and confirmed by blind testing wells.Study results reveal that the dominant sedimentary facies in the Qigequan Formation within the study area are delta front and shallow lake.The RGB fused slices indicate that there are two cycles with three sets of underwater distributary channel systems in one period.Among them,sandstones in the distributary channels of middle-low Qigequan Formation are thick and broad with superior physical properties,which are favorable reservoirs.The reservoir permeability is also affected by diagenesis.Distributary channel sandstone reservoirs extend further to the west of Sebei-1 gas field,which provides a basis to expand exploration to the western peripheral area.
基金funded by the National Natural Science Foundation of China(General Program:No.52074314,No.U19B6003-05)National Key Research and Development Program of China(2019YFA0708303-05)。
文摘Accurate prediction of formation pore pressure is essential to predict fluid flow and manage hydrocarbon production in petroleum engineering.Recent deep learning technique has been receiving more interest due to the great potential to deal with pore pressure prediction.However,most of the traditional deep learning models are less efficient to address generalization problems.To fill this technical gap,in this work,we developed a new adaptive physics-informed deep learning model with high generalization capability to predict pore pressure values directly from seismic data.Specifically,the new model,named CGP-NN,consists of a novel parametric features extraction approach(1DCPP),a stacked multilayer gated recurrent model(multilayer GRU),and an adaptive physics-informed loss function.Through machine training,the developed model can automatically select the optimal physical model to constrain the results for each pore pressure prediction.The CGP-NN model has the best generalization when the physicsrelated metricλ=0.5.A hybrid approach combining Eaton and Bowers methods is also proposed to build machine-learnable labels for solving the problem of few labels.To validate the developed model and methodology,a case study on a complex reservoir in Tarim Basin was further performed to demonstrate the high accuracy on the pore pressure prediction of new wells along with the strong generalization ability.The adaptive physics-informed deep learning approach presented here has potential application in the prediction of pore pressures coupled with multiple genesis mechanisms using seismic data.
基金The National Natural Science Foundation of China under contract No.41806048the Open Fund of the Hubei Key Laboratory of Marine Geological Resources under contract No.MGR202009+2 种基金the Fund from the Key Laboratory of Deep-Earth Dynamics of Ministry of Natural Resource,Institute of Geology,Chinese Academy of Geological Sciences under contract No.J1901-16the Aoshan Science and Technology Innovation Project of Pilot National Laboratory for Marine Science and Technology(Qingdao)under contract No.2015ASKJ03-Seabed Resourcesthe Fund from the Korea Institute of Ocean Science and Technology(KIOST)under contract No.PE99741.
文摘The South Yellow Sea basin is filled with Mesozoic-Cenozoic continental sediments overlying pre-Palaeozoic and Mesozoic-Palaeozoic marine sediments.Conventional multi-channel seismic data cannot describe the velocity structure of the marine residual basin in detail,leading to the lack of a deeper understanding of the distribution and lithology owing to strong energy shielding on the top interface of marine sediments.In this study,we present seismic tomography data from ocean bottom seismographs that describe the NEE-trending velocity distributions of the basin.The results indicate that strong velocity variations occur at shallow crustal levels.Horizontal velocity bodies show good correlation with surface geological features,and multi-layer features exist in the vertical velocity framework(depth:0–10 km).The analyses of the velocity model,gravity data,magnetic data,multichannel seismic profiles,and drilling data showed that high-velocity anomalies(>6.5 km/s)of small(thickness:1–2 km)and large(thickness:>5 km)scales were caused by igneous complexes in the multi-layer structure,which were active during the Palaeogene.Possible locations of good Mesozoic and Palaeozoic marine strata are limited to the Central Uplift and the western part of the Northern Depression along the wide-angle ocean bottom seismograph array.Following the Indosinian movement,a strong compression existed in the Northern Depression during the extensional phase that caused the formation of folds in the middle of the survey line.This study is useful for reconstructing the regional tectonic evolution and delineating the distribution of the marine residual basin in the South Yellow Sea basin.
基金supported by the National Natural Science Foundation of China(Grant Nos.91858215 and 41906048)。
文摘The Solomon Sea Basin is a Cenozoic back-arc spreading basin within the convergence system of the Pacific and Indo-Australian plates.Against the background of subduction polarity reversal,the current Solomon Sea Basin gradually formed a rhombic morphology with the subduction of the basin along the New Britain Trench and the Trobriand Trough.By analyzing the vertical gravity gradient,natural earthquake and seismic reflection data,this study determines the structural characteristics of the Solomon Sea Basin.It was found that the tectonics of the basin are characterized by the original expansion structure within the central part in addition to the structure induced by the latest subduction along the basin margin.The original spreading structure of the basin presented an east–west linear graben and horst controlled by normal faults during the basin expansion period.As a result of the subduction and slab-pull of the Solomon Sea Basin,extensional structure belts parallel to the New Britain Trench formed along the basin margin.
基金This study was supported by the National Natural Science Foundation of China under the project‘Research on the Dynamic Location of Receiver Points and Wave Field Separation Technology Based on Deep Learning in OBN Seismic Exploration’(No.42074140).
文摘At present,the acquisition of seismic data is developing toward high-precision and high-density methods.However,complex natural environments and cultural factors in many exploration areas cause difficulties in achieving uniform and intensive acquisition,which makes complete seismic data collection impossible.Therefore,data reconstruction is required in the processing link to ensure imaging accuracy.Deep learning,as a new field in rapid development,presents clear advantages in feature extraction and modeling.In this study,the convolutional neural network deep learning algorithm is applied to seismic data reconstruction.Based on the convolutional neural network algorithm and combined with the characteristics of seismic data acquisition,two training strategies of supervised and unsupervised learning are designed to reconstruct sparse acquisition seismic records.First,a supervised learning strategy is proposed for labeled data,wherein the complete seismic data are segmented as the input of the training set and are randomly sampled before each training,thereby increasing the number of samples and the richness of features.Second,an unsupervised learning strategy based on large samples is proposed for unlabeled data,and the rolling segmentation method is used to update(pseudo)labels and training parameters in the training process.Through the reconstruction test of simulated and actual data,the deep learning algorithm based on a convolutional neural network shows better reconstruction quality and higher accuracy than compressed sensing based on Curvelet transform.
基金Supported by the National Key R&D Program of China(No.2016YFC0303900)the Laoshan Laboratory(Nos.MGQNLM-KF201807,LSKJ202203604)the National Natural Science Foundation of China(No.42106072)。
文摘The Kuiyang-ST2000 deep-towed high-resolution multichannel seismic system was designed by the First Institute of Oceanography,Ministry of Natural Resources(FIO,MNR).The system is mainly composed of a plasma spark source(source level:216 dB,main frequency:750 Hz,frequency bandwidth:150-1200 Hz)and a towed hydrophone streamer with 48 channels.Because the source and the towed hydrophone streamer are constantly moving according to the towing configuration,the accurate positioning of the towing hydrophone array and the moveout correction of deep-towed multichannel seismic data processing before imaging are challenging.Initially,according to the characteristics of the system and the towing streamer shape in deep water,travel-time positioning method was used to construct the hydrophone streamer shape,and the results were corrected by using the polynomial curve fitting method.Then,a new data-processing workflow for Kuiyang-ST2000 system data was introduced,mainly including float datum setting,residual static correction,phase-based moveout correction,which allows the imaging algorithms of conventional marine seismic data processing to extend to deep-towed seismic data.We successfully applied the Kuiyang-ST2000 system and methodology of data processing to a gas hydrate survey of the Qiongdongnan and Shenhu areas in the South China Sea,and the results show that the profile has very high vertical and lateral resolutions(0.5 m and 8 m,respectively),which can provide full and accurate details of gas hydrate-related and geohazard sedimentary and structural features in the South China Sea.
文摘Paleostress plays a significant role in controlling the formation, accumulation, and distribution of reservoirs, and this could be an important factor in controlling the production of hydrocarbons from the unconventional reservoirs. In this study, we will use 3D seismic reflection data to perform the slip-tendency-based stress inversion to determine the stress field in the basement of the northern slope area in the Bongor Basin. The dataset for this technique is easily available in the oil and gas companies. The stress inversion results from the basement of the northern slope area of Bongor basin show that the maximum principal stress axis (σ1) is oriented vertically, the intermediate principal stress axis (σ2) is oriented N18° and the minimum principal stress axis (σ3) is oriented N105°, and σ2/σ1 = 0.60 and σ3/σ1 = 0.29. The findings of this paper provide significant information to understand the fault reactivation at the critical stage of hydrocarbon accumulation and the regional tectonic evolution.
基金supported by the National Natural Science Foundation of China(Nos.41506046,41376060,41706054)the Opening Foundation of Key Laboratory of Ocean and Marginal Sea Geology,CAS(No.MSGL15-05)+1 种基金WPOS(No.XDA11030102-02)the Strategic Priority Research Program of the Chinese Academy of Sciences(No.XDA13010101)
文摘The Pearl River Estuary(PRE) is located at the onshore-offshore transition zone between South China and South China Sea Basin, and it is of great significant value in discussing tectonic relationships between South China block and South China Sea block and seismic activities along the offshore active faults in PRE. However, the researches on geometric characteristics of offshore faults in this area are extremely lacking. To investigate the offshore fault distribution and their geometric features in the PRE in greater detail, we acquired thirteen seismic reflection profiles in 2015. Combining the analysis of the seismic reflection and free-air gravity anomaly data, this paper revealed the location, continuity, and geometry of the littoral fault zone and other offshore faults in PRE. The littoral fault zone is composed of the major Dangan Islands fault and several parallel, high-angle, normal faults, which mainly trend northeast to northeast-to-east and dip to the southeast with large displacements. The fault zone is divided into three different segments by the northwest-trending faults. Moreover, the basement depth around Dangan Islands is very shallow, while it suddenly increases along the islands westward and southward. These has resulted in the islands and neighboring areas becoming the places where the stress accumulates easily. The seismogenic pattern of this area is closely related to the comprehensive effect of intersecting faults together with the low velocity layer.
基金Funded by the Project of China Geological Survey (No.1212010916040)the Sichuan Science and Technology Program (No.2017JY0051)the Sichuan Science and Technology Program (No.2018GZ0200)
文摘The field seismic data is disturbed by the interferential information, which has low signal to noise ratio (SNR). That is disadvantage for seismic data interpretation. So it is important to remove the noise of seismic data. Independent component analysis (ICA) can remove most of the noise interference. However, ICA has some defects in noise reduction, because it needs some conditions that seismic data is independent reciprocally for denoising. To solve these defects, this paper proposes an improved ICA algorithm to noise reduction. Through simulation experiments, it can be obtained that the best decomposition levels of the new algorithm is 3. At last, the proposed improved ICA is applied to deal with the actual seismic data. The results show that it can effectively eliminate most of seismic noise such as random noise, linear interference, surface waves, and so on. The improved ICA is not only easy to denoising, but also has excellent mathematical theoretical properties.
文摘A new method is introduced to suppress the noise in seismic data processing. Based on the subtle difference in shape between the noise and the actual signal, we introduce morphologic filtering into seismic data processing. From the shape and the S/N we can see that the effect of morphologic filtering is superior to other methods like id-value filtering, neighbor average filtering, etc. The SNR of the signal after morphological filtering is comparatively great. In addition, the precision of the seismic data after morphological filtering is high. The characteristics of the actual signal, such as frequency and amplitude, are preserved. We give an example of the real seismic data processing using morphological filtering, in which the actual signal is retained, while the random high intensity noise was removed.
文摘Comprehensive inversion of logging and seismic data presented in this paper is a method to improve seismic data resolution. It involves using ample high-frequency information and complete low-frequency information of known logging to make up for the lack of limited bandwidth of practical seismic recording, obtaining an approximate reflection coefficient sequence (or wave impedance) of high resolution by iterative inversion and providing more reliable seismic evidence for further lithologic interpretation and lateral tracking, correlation and prediction of thin reservoir. The comprehensive inversion can be realized in the following steps: (1) to establish an initial model of higher resolution; (2) to obtain wavelets, and (3) to constrain iterative inversion. The key to this inversion lies in building an initial model. It is assumed from our experience that when the initial model is properly given, iterative inversion can be quickly converged to the ideal result.
文摘A novel time-domain identification technique is developed for the seismic response analysis of soil-structure interaction.A two-degree-of-freedom (2DOF) model with eight lumped parameters is adopted to model the frequency- dependent behavior of soils.For layered soil,the equivalent eight parameters of the 2DOF model are identified by the extended Kalman filter (EKF) method using recorded seismic data.The polynomial approximations for derivation of state estimators are applied in the EKF procedure.A realistic identification example is given for the layered-soil of a building site in Anchorage,Alaska in the United States.Results of the example demonstrate the feasibility and practicality of the proposed identification technique.The 2DOF soil model and the identification technique can be used for nonlinear response analysis of soil-structure interaction in the time-domain for layered or complex soil conditions.The identified parameters can be stored in a database tor use in other similar soil conditions,lfa universal database that covers information related to most soil conditions is developed in the thture,engineers could conveniently perform time history analyses of soil-structural interaction.
基金conducted under the HK science project 1-ZE8F:Remote-sensing data for studding the Earth’s and planetary inner structure.
文摘The Atlas region in northwest Africa is characterized by the Quaternary volcanism and elevated topography with past complex tectonic mitigation between the African and European plates.Geodynamics of this atypical region has left indubitably imprints in crustal architectonics,mainly regarding the crustal thickness as well as the crustal density structure.The knowledge of crustal thickness variations is of a significant interest,since it provides a crucial constraint to geodynamic and geophysical modelling of this region.In this study,we use gravity,topographic,bathymetric and sediment data together with results of seismic surveys to image the Moho topography beneath the Atlas region.The Bouguer gravity anomalies used for a gravimetric Moho recovery are obtained from the free-air gravity anomalies after subtracting the gravitational contributions of topography,bathymetry and sediments.The regional gravimetric Moho inversion constrained on seismic data is carried out by applying a regularized inversion technique based on Gauss-Newton’s formulation of improved Bott’s method,while adopting Earth’s spherical approximation.The numerical result reveals relatively significant Moho depth variations in the Moroccan Atlas,with minima of approximately 24 km along continental margins of the Mediterranean Sea and maxima exceeding 51 km beneath the Rif Cordillera.The Moho depth beneath the West African Craton varies from 32 km in its southern margin to 45 km beneath the Middle Atlas.The Tell Atlas is characterized by the shallow Moho depth of approximately 22 km and further deepening to 42 km towards the northern edge of the Aures Mountains.Our findings indicate a limited tectonic shortening of the High Atlas with the crustal thickness mostly within 36-42 km.Topographic discrepancies between the Rif Cordillera and the Atlas Mountains suggest that the hypothesis of isostatic compensation cannot be fully established.
基金supported by the Natural Science Foundation of Jiangsu Province(Grant No.SBK2020043202)by Key Laboratory of Geospace Environment and Geodesy,Ministry of Education,Wuhan University(No.19-01-08).
文摘On 21 May 2021(UTC),an MW 7.4 earthquake jolted the east Bayan Har block in the Tibetan Plateau.The earthquake received widespread attention as it is the largest event in the Tibetan Plateau and its surroundings since the 2008 Wenchuan earthquake,and especially in proximity to the seismic gaps on the east Kunlun fault.Here we use satellite interferometric synthetic aperture radar data and subpixel offset observations along the range directions to characterize the coseismic deformation of the earthquake.Range offset displacements depict clear surface ruptures with a total length of~170 km involving two possible activated fault segments in the earthquake.Coseismic modeling results indicate that the earthquake was dominated by left-lateral strike-slip motions of up to 7 m within the top 12 km of the crust.The well-resolved slip variations are characterized by five major slip patches along strike and 64%of shallow slip deficit,suggesting a young seismogenic structure.Spatial-temporal changes of the postseismic deformation are mapped from early 6-day and 24-day InSAR observations,and are well explained by time-dependent afterslip models.Analysis of Global Navigation Satellite System(GNSS)velocity profiles and strain rates suggests that the eastward extrusion of plateau is diffusely distributed across the east Bayan Har block,but exhibits significant lateral heterogeneities,as evidenced by magnetotelluric observations.The block-wide distributed deformation of the east Bayan Har block along with the significant co-and post-seismic stress loadings from the Madoi earthquake imply high seismic risks along regional faults,especially the Tuosuo Lake and Maqên-Maqu segments of the Kunlun fault that are known as seismic gaps.
文摘A novel technique for automatic seismic data processing using both integral and local feature of seismograms was presented in this paper. Here, the term integral feature of seismograms refers to feature which may depict the shape of the whole seismograms. However, unlike some previous efforts which completely abandon the DIAL approach, i.e., signal detection, phase identifi- cation, association, and event localization, and seek to use envelope cross-correlation to detect seismic events directly, our technique keeps following the DIAL approach, but in addition to detect signals corresponding to individual seismic phases, it also detects continuous wave-trains and explores their feature for phase-type identification and signal association. More concrete ideas about how to define wave-trains and combine them with various detections, as well as how to measure and utilize their feature in the seismic data processing were expatiated in the paper. This approach has been applied to the routine data processing by us for years, and test results for a 16 days' period using data from the Xinjiang seismic station network were presented. The automatic processing results have fairly low false and missed event rate simultaneously, showing that the new technique has good application prospects for improvement of the automatic seismic data processing.
文摘This paper is concerned with anisotropic effects on seismic data and signal analysis for transversely isotropic rock media with vertical anisotropy. It is understood that these effects are significant in many practical applications, e.g. earthquake forecasting, materials exploration inside the Earth’s crust, as well as various practical works in oil industry. Under the framework of the most accepted anisotropic media model (i.e. VTI media, transverse isotropy with a vertical axis symmetry), with applications of a set of available anisotropic rock parameters for sandstone and shale, we have performed numerical calculations of the anisotropic effects. We show that for rocks with strong anisotropy, the induced relative depth error can be significantly large. Nevertheless, with an improved understanding of the seismic-signal propagation and proper data processing, the error can be reduced, which in turn may enhance the probability of forecasting accurately the various wave propagations inside the Earth’s crust, e.g. correctly forecasting the incoming earthquakes from the center of the Earth.
基金supported by the National Science and Technology Major project(No.2016ZX05024001003)the Innovation Consortium Project of China Petroleum,and the Southwest Petroleum University(No.2020CX010201).
文摘Oil and gas seismic exploration have to adopt irregular seismic acquisition due to the increasingly complex exploration conditions to adapt to complex geological conditions and environments.However,the irregular seismic acquisition is accompanied by the lack of acquisition data,which requires high-precision regularization.The sparse signal feature in the transform domain in compressed sensing theory is used in this paper to recover the missing signal,involving sparse transform base optimization and threshold modeling.First,this paper analyzes and compares the effects of six sparse transformation bases on the reconstruction accuracy and efficiency of irregular seismic data and establishes the quantitative relationship between sparse transformation and reconstruction accuracy and efficiency.Second,an adaptive threshold modeling method based on sparse coefficient is provided to improve the reconstruction accuracy.Test results show that the method has good adaptability to different seismic data and sparse transform bases.The f-x domain reconstruction method of effective frequency samples is studied to address the problem of low computational efficiency.The parallel computing strategy of curvelet transform combined with OpenMP is further proposed,which substantially improves the computational efficiency under the premise of ensuring the reconstruction accuracy.Finally,the actual acquisition data are used to verify the proposed method.The results indicate that the proposed method strategy can solve the regularization problem of irregular seismic data in production and improve the imaging quality of the target layer economically and efficiently.
基金partly supported by the National Natural Science Foundation of China(grants no.41272132 and 41572080)the Fundamental Research Funds for central Universities(grant no.2-9-2013-97)the Major State Science and Technology Research Programs(grants no.2008ZX05056-002-02-01 and 2011ZX05010-001-009)
文摘Sandy debris flow deposits are present in Unit I during Miocene of Gas Field A in the Baiyun Depression of the South China Sea. The paucity of well data and the great variability of the sedimentary microfacies make it difficult to identify and predict the distribution patterns of the main gas reservoir, and have seriously hindered further exploration and development of the gas field. Therefore, making full use of the available seismic data is extremely important for predicting the spatial distribution of sedimentary microfacies when constructing three-dimensional reservoir models. A suitable reservoir modeling strategy or workflow controlled by sedimentary microfacies and seismic data has been developed. Five types of seismic attributes were selected to correlate with the sand percentage, and the root mean square (RMS) amplitude performed the best. The relation between the RMS amplitude and the sand percentage was used to construct a reservoir sand distribution map. Three types of main sedimentary microfacies were identified: debris channels, fan lobes, and natural levees. Using constraints from the sedimentary microfacies boundaries, a sedimentary microfacies model was constructed using the sequential indicator and assigned value simulation methods. Finally, reservoir models of physical properties for sandy debris flow deposits controlled by sedimentary microfacies and seismic inversion data were established. Property cutoff values were adopted because the sedimentary microfacies and the reservoir properties from well-logging interpretation are intrinsically different. Selection of appropriate reservoir property cutoffs is a key step in reservoir modeling when using simulation methods based on sedimentary microfacies control. When the abnormal data are truncated and the reservoir properties probability distribution fits a normal distribution, microfacies-controlled reservoir property models are more reliable than those obtained from the sequence Gauss simulation method. The cutoffs for effective porosity of the debris channel, fan lobe, and natural levee facies were 0.2, 0.09, and 0.12, respectively; the corresponding average effective porosities were 0.24, 0.13, and 0.15. The proposed modeling method makes full use of seismic attributes and seismic inversion data, and also makes the property data of single-well depositional microfacies more conformable to a normal distribution with geological significance. Thus, the method allows use of more reliable input data when we construct a model of a sandy debris flow.
基金Mainly presented at the 6-th international meeting of acoustics in Aug. 2003, and The 1999 SPE Asia Pacific Oil and GasConference and Exhibition held in Jakarta, Indonesia, 20-22 April 1999, SPE 54274.
文摘In this paper, a new concept called numerical structure of seismic data is introduced and the difference between numerical structure and numerical value of seismic data is explained. Our study shows that the numerical seismic structure is closely related to oil and gas-bearing reservoir,so it is very useful for a geologist or a geophysicist to precisely interpret the oil-bearing layers from the seismic data. This technology can be applied to any exploration or production stage. The new method has been tested on a series of exploratory or development wells and proved to be reliable in China. Hydrocarbon-detection with this new method for 39 exploration wells on 25 structures indicates a success ratio of over 80 percent. The new method of hydrocarbon prediction can be applied for: (1) depositional environment of reservoirs with marine facies, delta, or non-marine facies (including fluvial facies, lacustrine facies); (2) sedimentary rocks of reservoirs that are non-marine clastic rocks and carbonate rock; and (3) burial depths range from 300 m to 7000 m, and the minimum thickness of these reservoirs is over 8 m (main frequency is about 50 Hz).