Reducing migration distortions in migrated images of the Earth's subsurface. Recorded seismic data may be migrated, using a migration velocity model, to generate a migration image comprising ADCIGs with distortions. Synthetic seismic data may be generated, using the migration velocity model, for a grid of point scatterers. The synthetic seismic data may be migrated, using the migration velocity model, to generate impulse responses for the point scatterers. The impulse responses are used as point spread functions (PSFs) which approximate the blurring operator, e.g., the Hessian. An optimal reflectivity model may be selected using image-domain least-squares migration (LSM), based on the PSFs, with regularization of the difference between the migration image and a reflectivity model and a total variation (TV) regularization of the reflectivity model in the spatial and angular domains. An image of the optimal reflectivity model may be generated with reduced migration distortions compared to the original migration image.

A method is described for generating a subsurface model using stochastic full waveform inversion by receiving a seismic dataset representative of a subsurface volume of interest; performing stochastic full waveform inversion of the seismic dataset to generate a long wavelength subsurface model; and performing full waveform inversion of the seismic dataset using the long wavelength subsurface model as a starting model to generate an improved subsurface model. The method may further include performing seismic imaging of the seismic dataset using the improved subsurface model to generate a seismic image and identifying geologic features based on the seismic image. The method may be executed by a computer system.

One or more first sensors may be configured to sense seismic signals and one or more second sensors may be configured to sense electromagnetic crosstalk signals. The second sensors are not responsive to the seismic signals. The data from the first and second sensors may be recorded as first data and second data, respectively. The first data may be modified based on the second data to remove the electromagnetic crosstalk noise form the seismic data.

A method for seismic processing includes steps of seismic signal forward propagation and seismic data back propagation. The subsurface medium image is created after correlating and summarizing forward and backward propagation results. To address migration footprint and noise due to the incomplete data acquisition aperture and migration approximation in the migration operator, the iteration inversion strategy incorporates tensor flow calculated from seismic image. A regularization operator based on structure tensor of image is applied to seismic image inversion.

The present disclosure belongs to the technical field of seismic exploration imaging, and relates to an interpretive-guided velocity modeling seismic imaging method and system, a medium and a device. The method comprises the following steps: S1. performing first imaging on a given initial velocity model to obtain a first imaging result; S2. performing relative wave impedance inversion on the first imaging result to obtain a relative wave impedance profile; S3. performing Curvelet filtering on the relative wave impedance profile to obtain a first interpretation scheme; S4. superposing the first interpretation scheme and the initial velocity model to obtain a new migration velocity field; S5. performing second imaging on a new migration velocity field to obtain a second imaging result; and S6. repeating steps S2-S4 for the obtained second imaging result until a final seismic imaging result is obtained.

Aspects of the subject technology relate to systems and methods for improving cluster and surface efficiency in hydraulic fracturing by utilizing a stage optimization tool. Systems and methods are provided for receiving one or more perforation parameters of a wellbore, generating a perforation schema based on the one or more perforation parameters, training a stage optimization model based on the perforation schema to generate an optimized perforation schema, estimating a pressure of the wellbore based on the optimized perforation schema, and updating the optimized perforation schema until the estimated pressure is less than a predetermined pressure limit.

Methods and systems for processing seismic data are presented. Primary wave (P) seismic data (PP data) and shear wave (P) seismic data (PS data) are jointly inverted as part of a nonlinear tomography process which adheres to one or more co-depthing constraints.

This disclosure relates to geophysical exploration and seismic imaging, and more particularly to a seismic imaging method, system, and device based on pre-stack high-angle fast Fourier transform (FFT). The method includes: acquiring seismic data acquired during seismic exploration; extracting a common shot point gather from the seismic data followed by conversion into a frequency wavenumber domain common offset gather; calculating wave propagation angles; dividing an imaging region into a first region and a second region; solving constant coefficients of the first region and the second region; performing frequency-division layer-by-layer wavefield continuation on a frequency-wave number domain common offset gather to obtain imaging results at different depths and frequencies; subjecting the imaging results to integration followed by transformation to a spatial domain to obtain common offset imaging profiles; and subjecting the common offset imaging profiles to superposition obtain final imaging results.

Systems and methods include a method for generating a high-resolution advanced three-dimensional (3D) transient model that models multiple wells by integrating pressure transient data into a static geological model. A crude 3D model is generated from a full-field geological model that models production for multiple wells in an area. A functional 3D model is generated from the crude 3D model. An intermediate 3D model is generated by calibrating the functional 3D model with single-well data. An advanced 3D transient model is generated by calibrating multi-well data in the functional 3D model.

Methods and systems are disclosed for forming an image of a subterranean region of interest. The method includes receiving an observed seismic dataset and a seismic velocity model for the subterranean region of interest, and generating a simulated seismic dataset based on the seismic velocity model and the geometry of the observed seismic dataset. The method further includes determining a transformed observed seismic dataset by applying a nonlinear amplitude transform to the observed seismic dataset and determining a transformed simulated seismic by applying the same transform to the simulated seismic dataset. The method still further includes forming an objective function based on the transformed observed seismic and the transformed simulated seismic dataset, and determining an updated seismic velocity model based upon finding an extremum of the objective function.