An adaptive physics-informed deep learning method for pore pressure prediction using seismic data

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.

Application of 9-component S-wave 3D seismic data to study sedimentary facies and reservoirs in a biogasbearing area:A case study on the Pleistocene Qigequan Formation in Taidong area,Sanhu Depression,Qaidam Basin,NW China

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.

Robust high frequency seismic bandwidth extension with a deep neural network trained using synthetic data

Geophysicists interpreting seismic reflection data aim for the highest resolution possible as this facilitates the interpretation and discrimination of subtle geological features.Various deterministic methods based on Wiener filtering exist to increase the temporal frequency bandwidth and compress the seismic wavelet in a process called spectral shaping.Auto-encoder neural networks with convolutional layers have been applied to this problem,with encouraging results,but the problem of generalization to unseen data remains.Most published works have used supervised learning with training data constructed from field seismic data or synthetic seismic data generated based on measured well logs or based on seismic wavefield modelling.This leads to satisfactory results on datasets similar to the training data but requires re-training of the networks for unseen data with different characteristics.In this work seek to improve the generalization,not by experimenting with network architecture(we use a conventional U-net with some small modifications),but by adopting a different approach to creating the training data for the supervised learning process.Although the network is important,at this stage of development we see more improvement in prediction results by altering the design of the training data than by architectural changes.The approach we take is to create synthetic training data consisting of simple geometric shapes convolved with a seismic wavelet.We created a very diverse training dataset consisting of 9000 seismic images with between 5 and 300 seismic events resembling seismic reflections that have geophysically motived perturbations in terms of shape and character.The 2D U-net we have trained can boost robustly and recursively the dominant frequency by 50%.We demonstrate this on unseen field data with different bandwidths and signal-to-noise ratios.Additionally,this 2D U-net can handle non-stationary wavelets and overlapping events of different bandwidth without creating excessive ringing.It is also robust in the presence of noise.The significance of this result is that it simplifies the effort of bandwidth extension and demonstrates the usefulness of auto-encoder neural network for geophysical data processing.

Irregular seismic data reconstruction based on exponential threshold model of POCS method

Irregular seismic data causes problems with multi-trace processing algorithms and degrades processing quality. We introduce the Projection onto Convex Sets （POCS） based image restoration method into the seismic data reconstruction field to interpolate irregularly missing traces. For entire dead traces, we transfer the POCS iteration reconstruction process from the time to frequency domain to save computational cost because forward and reverse Fourier time transforms are not needed. In each iteration, the selection threshold parameter is important for reconstruction efficiency. In this paper, we designed two types of threshold models to reconstruct irregularly missing seismic data. The experimental results show that an exponential threshold can greatly reduce iterations and improve reconstruction efficiency compared to a linear threshold for the same reconstruction result. We also analyze the anti- noise and anti-alias ability of the POCS reconstruction method. Finally, theoretical model tests and real data examples indicate that the proposed method is efficient and applicable.

Design of multi-channel data acquisition system based on LabVIEW and NI PXI-5105

A data acquisition system based on LabVIEW and NI PXI-5105 is presented for multi-channel data acquisition. It can realize the functions of parameter setting, data acquisition and storage, waveform display and data analysis using LabVIEW and NI-SCOPE device driver. The advantages of the system are that the setting is convenient, the operation is easy, the interface is friendly and the functions are practical. The experiment results show that the system has good stability and high reliability and is a powerful tool for multi-channel data acquisition.

Separation of P. and SV-wavefields from multi-componen seismic data

In multi-component seismic exploration, the horizontal and vertical components both contain P- and SV-waves. The P- and SV-wavefields in a seismic record can be separated by their horizontal and vertical displacements when upgoing P- and SV-waves arrive at the sea floor. If the sea floor P wave velocity, S wave velocity, and density are known, the separation can be achieved in ther-p domain. The separated wavefields are then transformed to the time domain. A method of separating P- and SV-wavefields is presented in this paper and used to effectively separate P- and SV-wavefields in synthetic and real data. The application to real data shows that this method is feasible and effective. It also can be used for free surface data.

Offshore Fault Geometrics in the Pearl River Estuary,Southeastern China:Evidence from Seismic Reflection Data

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 Analysis of Seismic Data Structure and Oil and Gas Prediction

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 indi- cates a success ratio of over 80 percent. The new method of hydrocarbon prediction can be applied for: (1) depositional environment of reservoirs with marine fades, delta, or non-marine fades (including fluvial facies, lacustrine fades); (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).

The improved ICA algorithm and its application in the seismic data denoising

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.

Method for Morphological Filtering in Seismic Data Processing

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.

Application of multi-scaled morphology in denoising seismic data

In this paper, multi-scaled morphology is introduced into the digital processing domain for land seismic data. First, we describe the basic theory of multi-scaled morphology image decomposition of exploration seismic waves; second, we illustrate how to use multi-scaled morphology for seismic data processing using two real examples. The first example demonstrates suppressing the surface waves in pre-stack seismic records using multi-scaled morphology decomposition and reconstitution and the other example demonstrates filtering different interference waves on the seismic record. Multi-scaled morphology filtering separates signal from noise by the detailed differences of the wave shapes. The successful applications suggest that multi-scaled morphology has a promising application in seismic data processing.

Predicting the distribution of reservoirs by applying the method of seismic data structure characteristics： Example from the eighth zone in Tahe Oilfield

Seismic data structure characteristics means the waveform character arranged in the time sequence at discrete data points in each 2-D or 3-D seismic trace. Hydrocarbon prediction using seismic data structure characteristics is a new reservoir prediction technique. When the main pay interval is in carbonate fracture and fissure-cavern type reservoirs with very strong inhomogeneity, there are some difficulties with hydrocarbon prediction. Because of the special geological conditions of the eighth zone in the Tahe oil field, we apply seismic data structure characteristics to hydrocarbon prediction for the Ordovician reservoir in this zone. We divide the area oil zone into favorable and unfavorable blocks. Eighteen well locations were proposed in the favorable oil block, drilled, and recovered higher output of oil and gas.

A Method of Improving Seismic Data Resolution：Comprehensive Inversion of Well logging and Seismic Data

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.

Time-domain identification of dynamic properties of layered soil by using extended Kalman filter and recorded seismic data

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.

Crustal thickness beneath Atlas region from gravity,topographic,sediment and seismic data

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.

Resolving co- and early post-seismic slip variations of the 2021 MW 7.4 Madoi earthquake in east Bayan Har block with a block-wide distributed deformation mode from satellite synthetic aperture radar data

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.

Microseismic source location using the Log-Cosh function and distant sensor-removed P-wave arrival data

Source location is the core foundation of microseismic monitoring.To date,commonly used location methods have usually been based on the ray-tracing travel-time technique,which generally adopts an L1 or L2 norm to establish the location objective function.However,the L1 norm usually achieves low location accuracy,whereas the L2 norm is easily affected by large P-wave arrival-time picking errors.In addition,traditional location methods may be affected by the initial iteration point used to find a local optimum location.Furthermore,the P-wave arrival-time data that have travelled long distances are usually poor in quality.To address these problems,this paper presents a microseismic source location method using the Log-Cosh function and distant sensor-removed P-wave arrival data.Its basic principles are as follows:First,the source location objective function is established using the Log-Cosh function.This function has the stability of the L1 norm and location accuracy of the L2 norm.Then,multiple initial points are generated randomly in the mining area,and the established Log-Cosh location objective function is used to obtain multiple corresponding location results.The average value of the 50 location points with the largest data field potential values is treated as the initial location result.Next,the P-wave travel times from the initial location result to triggered sensors are calculated,and then the P-wave arrival data with travel times exceeding 0.2 s are removed.Finally,the aforementioned location steps are repeated with the denoised P-wave arrival dataset to obtain a high-precision location result.Two synthetic events and eight blasting events from the Yongshaba mine,China,were used to test the proposed method.Regardless of whether the P-wave arrival data with long travel times were eliminated,the location error of the proposed method was smaller than that of the L1/L2 norm and trigger-time-based location method(TT1/TT2 method).Furthermore,after eliminating the Pwave arrival data with long travel distances,the location accuracy of these three location methods increased,indicating that the proposed location method has good application prospects.

A novel technique for automatic seismic data processing using both integral and local feature of seismograms

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.

Reconstruction method of irregular seismic data with adaptive thresholds based on different sparse transform bases

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.

Application of Catastrophe Theory in 3D Seismic Data Interpretation of Coal Mine

In order to detect fault exactly and quickly, cusp catastrophe theory is used to interpret 3D coal seismic data in this paper. By establishing a cusp model, seismic signal is transformed into standard form of cusp catastrophe and catastrophe parameters, including time-domain catastrophe potential, time-domain catastrophe time, frequency-domain catastrophe potential and frequency- domain degree, are calculated. Catastrophe theory is used in 3D seismic structural interpretation in coal mine. The results show that the position of abnormality of the catastrophe parameter profile or curve is related to the location of fault, and the cusp catastrophe theory is effective to automatically pick up geology information and improve the interpretation precision in 3D seismic data.