A computer-implemented method, and system implementing the method, are disclosed for computing a final model of elastic properties, using nonlinear direct prestack seismic inversion for Poisson impedance. User inputs and earth-model data is obtained over points of incidence of a survey region, at various angles of incidence. Various models are then computed that serve for lithology identification and fluid discrimination and take part in preliminary seismic exploration and reservoir characterization. Therefore, further refinement of these models is required due to changes in burial depths, compaction and overburden pressure, as they provide limitations for reservoirs on porous media. The further refinement using nonlinear direct prestack seismic model is performed on a system computer, which produces a final model of elastic properties. This model can then be applied for lithology prediction and fluid detection to identify potential targets of oil and gas exploration and estimating spots in unconventional shale gas applications.

Inversion for marine seismic imaging of a full reflected wavefield can include generating an image of a subsurface formation by full wavefield migrating a recorded seismic wavefield and generating numerically modeled data using the image. A mismatch between the numerically modeled data and the seismic wavefield can be determined. Responsive to determining that the mismatch exceeds an inversion match threshold, the image can be updated using the mismatch between the numerically modeled data and the seismic wavefield.

Methods and apparatuses for seismic exploration of an underground structure obtain improved images by integrating partial match filtering in an FWI. Filtered (auxiliary) data replaces one of the observed data and the synthetic data in the FWI's objective function to avoid cycle skipping.

A method for correcting physical positions of seismic sensors and/or seismic sources for a seismic data acquisition system. The method includes estimating a respective energy generated by each source element, which belongs to a source array; calculating a respective energy recorded by each individual seismic sensor, which belongs to a composite receiver; summing, for each individual seismic sensor, all the generated energies from the all the source elements; estimating a model of direct arrival waves that propagate from the source elements to the individual seismic sensors; calculating positions of the individual seismic sensors based on the model of direct arrival waves; comparing calculated positions of the individual seismic sensors with observed positions of the individual seismic sensors; selecting a best calculated position for each of the individual seismic sensors based on an objective function; and correcting the observed positions of the individual seismic sensors with corresponding best calculated positions.

A method is described for image-domain full waveform inversion. The method may include receiving, at a computer processor, a seismic dataset representative of the subsurface volume of interest and an initial earth model; performing, via the computer processor, an image domain full waveform inversion to generate an updated earth model; and performing, via the computer processor, seismic imaging of the seismic dataset using the updated earth model to generate a seismic image. The method may be executed by a computer system.

Methods and apparatuses for seismic exploration of an underground structure obtain improved images by integrating partial match filtering in an FWI. Filtered (auxiliary) data replaces one of the observed data and the synthetic data in the FWI's objective function to avoid cycle skipping.

Systems and methods of detecting marine seismic survey parameters are provided. A data processing system can obtain seismic data from seismic data acquisition units disposed on a seabed responsive to an acoustic signal propagated from an acoustic source through a water column. The data processing system can determine from the seismic data, a direct arrival time for the acoustic signal at each of the plurality of seismic data acquisition units, and can obtain an estimated depth value of each of the plurality of seismic data acquisition units and an estimated water column transit velocity of the acoustic signal. The data processing system can apply a depth model and a water column transit velocity model to the estimated depth value and to the estimated water column transit velocity determine an updated depth value and an updated water column transit velocity for each of the plurality of seismic data acquisition units.

The disclosure herein discloses an imaging method of internal defects in longitudinal sections of trees, and belongs to the field of nondestructive testing of trees. The method includes the following steps: with the propagation time of stress waves in a tree as input data, dividing an imaging plane into a predetermined number of grid cells to establish initial velocity distribution in the imaging plane; then performing multiple iterations using a linear propagation model; following each iteration, adjusting the velocity distribution in the imaging plane using the SIRT algorithm; constraining the velocity of each grid cell using maximum and minimum velocity constraints and fuzzy constraints based on grid cell groups, and ending iteration until the final velocity distribution is in good fit with the measured data; by comparing the velocity value of the grid cell at this moment with the reference value of the tested healthy tree, determining an abnormal grid cell; and then performing secondary smoothing processing on the grid cell imaging to obtain the defect location inside the tree. The method can accurately detect the defective area of the tree, and has less false detection areas and good imaging effect.

A method, including: obtaining initial estimates of a plurality of subsurface parameters; obtaining a recorded wavefield decomposed into a plurality of discrete components; performing, with a computer, a cascaded inversion where the initial estimates of the subsurface parameters are individually updated, wherein each of the subsurface parameters are updated using a different discrete component of the recorded wavefield of the plurality of discrete components; and generating, with the computer, updated subsurface models from the cascaded inversion for each of the subsurface parameters.

Systems and methods of detecting marine seismic survey parameters are provided. A data processing system can obtain seismic data from seismic data acquisition units disposed on a seabed responsive to an acoustic signal propagated from an acoustic source through a water column. The data processing system can determine from the seismic data, a direct arrival time for the acoustic signal at each of the plurality of seismic data acquisition units, and can obtain an estimated depth value of each of the plurality of seismic data acquisition units and an estimated water column transit velocity of the acoustic signal. The data processing system can apply a depth model and a water column transit velocity model to the estimated depth value and to the estimated water column transit velocity determine an updated depth value and an updated water column transit velocity for each of the plurality of seismic data acquisition units.