Seismic exploration of an underground formation uses seismic excitations to probe the formation's properties such as reflectivity that can be imaged using reverse time migration. Using an equal area spherical binning at reflection points improves and simplifies RTM imaging together with adaptability to the data acquisition geometry, while overcoming drawbacks of conventional cylindrical binning.

The present disclosure provides a method and system for stabilizing a Poynting vector of a seismic wavefield. The method includes: adjusting an amplitude of a time derivative of the seismic wavefield, and computing a Poynting vector of the adjusted time derivative of the seismic wavefield to obtain a first Poynting vector, where a difference between the amplitude of the first Poynting vector and the amplitude of a second Poynting vector is within a set range, and the second Poynting vector belongs to the seismic wavefield; and conducting operation on the second Poynting vector and the first Poynting vector to obtain a final Poynting vector of the seismic wavefield. The present disclosure addresses instability of Poynting vectors computation.

A method may include obtaining seismic data regarding a geological region of interest. The method may further include obtaining a property model regarding the geological region of interest. The method may further include determining an adjoint migration operator based on the property model. The method may further include updating the property model using the seismic data and a conjugate gradient solver in a least-squares reverse time migration to produce a first updated property model. The conjugate gradient solver is based on the adjoint migration operator. The method may further include updating the first updated property model using a threshold shrinkage function to produce a second updated property model. The threshold shrinkage function comprises a sign function and a maximum function that are applied to the first updated property model. The method may further include generating a seismic image of the geological region of interest using the second updated property model.

A method may include obtaining seismic data regarding a geological region of interest. The method may further include obtaining a property model regarding the geological region of interest. The method may further include determining an adjoint migration operator based on the property model. The method may further include updating the property model using the seismic data and a conjugate gradient solver in a least-squares reverse time migration to produce a first updated property model. The conjugate gradient solver is based on the adjoint migration operator. The method may further include updating the first updated property model using a threshold shrinkage function to produce a second updated property model. The threshold shrinkage function comprises a sign function and a maximum function that are applied to the first updated property model. The method may further include generating a seismic image of the geological region of interest using the second updated property model.

Methods and platforms for allowing efficient identification of 3D stratigraphic models that explain observed log data.

Methods and platforms for allowing efficient identification of 3D stratigraphic models that explain observed log data.

A method for identifying a reef-shoal reservoir in a faulted lacustrine basin based on a basement structure-paleogeomorphology-seismic facies progressive constraint, including: analyzing a basement structure of a work area; establishing a paleogeomorphology classification standard according to thickness, reflection structure and stratigraphic dip; based on well-seismic calibration and forward modeling, establishing a seismic facies classification standard for reef-shoal facies belts under different paleo-geomorphic conditions, and quantitatively predicting and describing a reservoir in the reef-shoal facies belts using seismic facies-controlled inversion; and according to analysis results of basement structure characteristic, paleogeomorphology classification and seismic facies, establishing a method for predicting a favorable reservoir.

A method for identifying a reef-shoal reservoir in a faulted lacustrine basin based on a basement structure-paleogeomorphology-seismic facies progressive constraint, including: analyzing a basement structure of a work area; establishing a paleogeomorphology classification standard according to thickness, reflection structure and stratigraphic dip; based on well-seismic calibration and forward modeling, establishing a seismic facies classification standard for reef-shoal facies belts under different paleo-geomorphic conditions, and quantitatively predicting and describing a reservoir in the reef-shoal facies belts using seismic facies-controlled inversion; and according to analysis results of basement structure characteristic, paleogeomorphology classification and seismic facies, establishing a method for predicting a favorable reservoir.

A method and apparatus for subsurface data processing includes determining a set of clusters based at least in part on measurement vectors associated with different depths or times in subsurface data, defining clusters in a subsurface data by classes associated with a state mode, reducing a quantity of the subsurface data based at least in part on the classes, and storing the reduced quantity of the subsurface data and classes with the state model in a training database for a machine learning process.

A method is claimed that includes obtaining a measured present-day value of at least one parameter for each member of a set of unvalidated geological layers arranged in order of increasing depth and iteratively selecting a member of the set as a current layer. For each current layer in turn, the method further determines an estimated archaic value of at least one parameter of the current layer based on its measured present-day value by applying an alternating cycle of decompaction followed by geomechnical modeling to predict a present-day value of the parameter of the current layer based on its estimated archaic value. The method still further determines a validated archaic value of at least one parameter of each current layer based on a difference between the predicted and the measured present-day values. A non-transitory computer readable medium storing instructions for validating the archaic value for each layer is claimed.