A method of fluid substitution, wherein an initial data set is provided, wherein a substituted data set is provided, wherein a rock physics model is provided, wherein the initial data set includes initial data of a geophysical parameter and initial fluid data, and wherein the substituted data set includes substituted fluid data. The method includes using the model and the initial data set to calculate first calculated data of the geophysical parameter, using the model and the substituted data set to calculate second calculated data of the geophysical parameter, calculating the difference between the first calculated data of the geophysical parameter and the second calculated data of the geophysical parameter, and applying the difference to the initial data of the geophysical parameter to produce substituted data of the geophysical parameter.

According to one embodiment, a near-field data is used to determine a taper length that can isolate the source signature at the top of near-field data with minimum interaction with the Green's function. In some embodiments, a range of taper lengths is selected and for each length after tapering the near-filed data, converting each filtered near-field data to its minimum-phase equivalents. Summing pairwise cross-correlation of all of the minimum-phase equivalent wavelets at the zero-lag provides an attribute that shows how much the tapered portions of the near-field data look alike. An acceptable taper size will be the one that has the highest summation value. Finally, the average of the minimum-phase equivalents of tapered near-field data with the selected taper size is the estimated source signature.

Disclosed are methods, systems, and computer-readable medium to perform operations including: generating, using a source wavelet and a current velocity model, modeled seismic data of the subterranean formation; applying a pre-condition to a seismic data residual calculated using the modeled seismic data and acquired seismic data from the subterranean formation; generating a velocity update using the source wavelet and the pre-conditioned seismic data residual; updating, using the velocity update, the current velocity model to generate an updated velocity model; determining that the current velocity model satisfies a predetermined condition; and responsively determining that the updated velocity model is the velocity model of the subterranean formation.

A method and apparatus for generating an image of a subsurface region including obtaining geophysical data/properties for the subsurface region; resampling the geophysical data/properties to generate a resampled data set; iteratively (a) inverting the resampled data set with an initial prior model to generate a new model; and (b) updating the new model based on learned information to generate an updated prior model; substituting the initial prior model in each iteration with the updated prior model from an immediately-preceding iteration; and determining an end point for the iteration. A final updated model may thereby be obtained, which may be used in managing hydrocarbons. Inversion may be based upon linear physics for the first one or more iterations, while subsequent iterations may be based upon non-linear physics.

Embodiments of the present invention help in the processing and interpretation of seismic survey data, by correlating or otherwise comparing or associating seismic data obtained from a seismic survey with flow information obtained from a well or borehole in the surveyed area. In particular, embodiments of the present invention allow for flow data representing a flow profile along a well that is being monitored by a distributed acoustic sensor to be determined, such that regions of higher flow in the well can be determined. For example, in the production zone the well will be perforated to allow oil to enter the well, but it has not previously been possible to determine accurately where in the production zone the oil is entering the well. However, by determining a flow rate profile along the well using the DAS then this provides information as to where in the perforated production zone oil is entering the well, and hence the location of oil bearing sands. This location can then be combined or otherwise correlated, used, or associated with petroleum reservoir location information obtained from the seismic survey, to improve the confidence and/or accuracy in the determined petroleum reservoir location.

A method for estimating sonic slowness comprising: obtaining (700) a plurality of sonic waveforms are received by a plurality of receivers of a logging tool after emission of a source sonic wave by a transmitter, obtaining (710) slowness models of the subterranean formation, a slowness model being defined by a at least one cell of constant slowness for at least one wave energy mode, computing (720), for each slowness model, a set of candidate travel times, a candidate travel time of a set of candidate travel times being computed for a wave energy mode and a position of a receiver of the plurality of receivers, computing (730) a relevance indicator for each set of candidate travel times based on the recorded sonic waveforms; searching (740) a match between the sets of candidate travel times and the recorded sonic waveforms by searching a relevance indicator which is optimum, computing (750) a sonic slowness estimate for the subterranean formation from a set of candidate travel times for which the relevance indicator is optimum.

In order to monitor the subsoil of the earth under a target zone, seismic waves coming from an identified mobile noise source are recorded by means of at least one pair of sensors disposed on either side of the target zone, time periods are selected corresponding to the alignments of the pairs of sensors with the noise source, a seismogram of the target zone is reconstructed by interferometry based on the recorded seismic waves and on the selected time periods and an image of the subsoil of the target zone is generated using the seismogram.

Methods are provided for estimating a quality of cement in the annuli of a multi-string wellbore. Wideband acoustic energy signals are generated and detected in the wellbore and are processed to obtain indications of wideband casing-formation phase slowness dispersions in the wellbore. The indications are compared to reference wideband model casing-formation phase slowness dispersions in order to estimate status of cement or lack of cement in the annuli at that location based on the results of the comparison.

A computer-implemented method is described for seismic facies identification including receiving a seismic dataset representative of a subsurface volume of interest; applying a model conditioned by petrophysical classifications to the seismic dataset to identify seismic facies and generate a classified seismic image; and identifying geologic features based on the classified seismic image. The method generates seismic facies probability volumes.

A method for determining a favorable time window of an infill well of an unconventional oil and gas reservoir, which comprises the following steps: S1, establishing a three-dimensional geological model with physical properties and geomechanical parameters; S2, establishing a natural fracture network model in combination with indoor core-logging-seismic monitoring; S3, calculating complex fractures in hydraulic fracturing of parent wells; S4, establishing an unconventional oil and gas reservoir model and calculating a current pore pressure field; S5, establishing a dynamic geomechanical model and calculating a dynamic geostress field; S6, calculating complex fractures in horizontal fractures of the infill well in different production times of the parent wells based on pre-stage complex fractures and the current geostress field; S7, analyzing a microseismic event barrier region and its dynamic changes in infill well fracturing; and S8, analyzing the productivity in different infill times, and determining an infill time window.