Systems and methods for quantitative seismic integrated modelling (QSIM) are disclosed for integrating the one, two and three-dimensional (1D, 2D, 3D) data from different geoscience domains within a framework in order to produce hi-resolution geocellular models that simulate realistic sub-surface reservoir properties. The QSIM systems and methods accurately leverage the seismically derived reservoir rock properties, integrating the geophysical, geological and engineering information through an optimum rock physics models and takes in consideration all the empirically constrained templates to correct, validate and quality check all the input data.

The present disclosure provides a method and a device for imaging diffracted waves based on azimuth-dip angle gathers and a storage medium, which relates to the technical field of seismic exploration, comprising firstly acquiring seismic data and generating target azimuth-dip angle gathers based on the seismic data, wherein the target azimuth-dip angle gathers are a set of all azimuth-dip angle gathers in which the Fresnel zones have been muted, and each of the azimuth-dip angle gathers represents a dip-angle gather corresponding to each azimuth angle; then detecting diffracted waves based on the target azimuth-dip angle gathers, and determining the type of the diffracted waves; and finally, imaging the diffracted waves based on the type of the diffracted waves to obtain a diffracted wave imaging result.

The present disclosure provides a method and a device for imaging diffracted waves based on azimuth-dip angle gathers and a storage medium, which relates to the technical field of seismic exploration, comprising firstly acquiring seismic data and generating target azimuth-dip angle gathers based on the seismic data, wherein the target azimuth-dip angle gathers are a set of all azimuth-dip angle gathers in which the Fresnel zones have been muted, and each of the azimuth-dip angle gathers represents a dip-angle gather corresponding to each azimuth angle; then detecting diffracted waves based on the target azimuth-dip angle gathers, and determining the type of the diffracted waves; and finally, imaging the diffracted waves based on the type of the diffracted waves to obtain a diffracted wave imaging result.

Methods, systems, and computer readable media for gridding of subsurface salt structures include determining a predetermined area lacks three dimensional seismic coverage, generating a two dimensional seismic top salt interpretation for the predetermined area, generating a bathymetry elevation of the predetermined area, determining that at least one two dimensional seismic line intersects a bathymetric feature of interest, and determining a correlation coefficient between the two dimensional seismic top salt interpretation and the bathymetry elevation. The method may further include determining the correlation coefficient is greater than a predetermined threshold value, and applying the bathymetry elevation as an additional control for gridding top of the subsurface salt structure. The step of gridding the top of the subsurface salt structure may further include applying at least one of kriging with external drift (KED), polygon-based approaches, regression-kriging, and other geostatistical methods.

A seismic time-frequency analysis method based on generalized Chirplet transform with time-synchronized extraction, which has higher level of energy aggregation in the time direction and can better describe and characterize the local characteristics of seismic signals, and is applicable to the time-frequency characteristic representation of both harmonic signals and pulse signals, comprising the steps of processing generalized Chirplet transform with time-synchronized extraction for each seismic signal to obtain a time spectrum by: carrying out generalized Chirplet transform, calculating group delay operator and carrying out time-synchronized extraction on seismic signals, thereby the boundary and heterogeneity structure of the rock slice are more accurately and clearly shown and subsequence seismic analysis and interpretation are facilitated.

A zero-offset wavefield synthesis workflow to calculate a synthesized zero-offset wavefield output without the commitment to an rms velocity field output to circumvent velocity uncertainty. Said zero-offset wavefield synthesis workflow comprises calculating a migration cube output. Rendering a demigration cube output from said migration cube output with a demigration cube calculation. Rendering said synthesized zero-offset wavefield output from said demigration cube output with a zero-offset wavefield synthesis procedure.

A zero-offset wavefield synthesis workflow to calculate a synthesized zero-offset wavefield output without the commitment to an rms velocity field output to circumvent velocity uncertainty. Said zero-offset wavefield synthesis workflow comprises calculating a migration cube output. Rendering a demigration cube output from said migration cube output with a demigration cube calculation. Rendering said synthesized zero-offset wavefield output from said demigration cube output with a zero-offset wavefield synthesis procedure.

Processes and systems for deblending blended seismic data with attenuated source signatures and source ghost are described. Processes and systems compute blended upgoing pressure wavefield based on blended pressure wavefield and blended vertical velocity wavefield recorded in a marine survey of a subterranean formation. Downgoing vertical velocity wavefield is computed based on near-field pressure wavefields generated by source elements of sources activated in the marine survey. Deblended wavefield is computed based on the blended upgoing pressure wavefield and the downgoing vertical velocity source wavefield. The deblended wavefield may be used to generate an image of the subterranean formation with the source signatures and source ghosts contained in the blended pressure wavefield and blended vertical velocity wavefield.

Processes and systems for deblending blended seismic data with attenuated source signatures and source ghost are described. Processes and systems compute blended upgoing pressure wavefield based on blended pressure wavefield and blended vertical velocity wavefield recorded in a marine survey of a subterranean formation. Downgoing vertical velocity wavefield is computed based on near-field pressure wavefields generated by source elements of sources activated in the marine survey. Deblended wavefield is computed based on the blended upgoing pressure wavefield and the downgoing vertical velocity source wavefield. The deblended wavefield may be used to generate an image of the subterranean formation with the source signatures and source ghosts contained in the blended pressure wavefield and blended vertical velocity wavefield.

Systems and methods include a computer-implemented method for presenting interpretation results of synchronized seismic data and fracture treatment times. A standard format seismic dataset of sensor readings obtained from a three-component sensor is generated. Coordinates and recording times corresponding to the sensor readings are added to the standard format seismic dataset. Synchronized seismic data is generated from the standard format seismic dataset by synchronizing seismic recording times with fracture treatment times. Quality-controlled synchronized seismic data is generated by removing dead traces and abnormal data samples from the synchronized seismic data. A time-frequency analysis is performed on the quality-controlled synchronized seismic data, including performing short-time Fourier transforms to analyze variations of Fourier spectra over time. Based on the time-frequency analysis, resonance frequencies are extracted from each frequency spectrum at different time samples. Interpretation results based are presented to a user.