A computer-implemented method may include obtaining seismic data acquired in a time-domain for a subterranean region of interest. The method may further include obtaining a property model for the subterranean region of interest. The method may further include determining one or more time shifts using a segment dynamic image warping function based on the seismic data and the property model. The method may further include determining an adjoint source operator using the derived time shift and one-way wave equation. The method may further include updating the property model using a gradient solver in a data-domain reflection traveltime inversion. The method may further include outputting the updated property model for the subterranean region of interest. The method may further include generating a seismic image for the subterranean region of interest using the updated property model.

A method is presented wherein inversion for formation anisotropic constants is achieved using borehole sonic data.

In a first step, a defined scope value is selected for each of a plurality of hydrodynamic input parameters. A simulated topographical result is generated using the selected scope values and a forward model. A detailed seismic interpretation is generated to represent specific seismic features or observed topography. A calculated a misfit value representing a distance between the simulated topographical result and a detailed seismic interpretation is minimized. An estimated optimized sand ratio and optimized hydrodynamic input parameters are generated. In a second step, a genetic algorithm is used to determine a proportion of each grain size in the estimated optimized sand ratio. A misfit value is used that is calculated from thickness and porosity data extracted from well data and a simulation result generated by the forward model to generate optimized components of different grain sizes. Optimized hydrodynamic input parameters and optimized components of different grain sizes are generated.

A method for determining a thermal maturity image of a subterranean region and a non-transitory computer readable medium, storing instructions for executing the method, are disclosed. The method includes, obtaining a seismic dataset for the subterranean region of interest, obtaining a thermal maturity value for a plurality of core samples taken from different positions within the subterranean region, and obtaining a plurality of well log types from the core sampling location. The method further includes determining a calibrated rock physics model based on the plurality of well log types, determining a pore fluid type based on the calibrated rock physics model, and determining a thermal maturity model based on the plurality of core samples, on the pore fluid type, and on the plurality of well logs. The method still further includes determining the thermal maturity image of the subterranean region based on the seismic dataset and thermal maturity model.

A method is disclosed for radiaiiy profiling shear velocities of flexural wave modes in a formation. The method includes establishing sensitivity kernels with two non-dimensionalized parameters and using said sensitivity kernels to perform an inversion for radial shear wave velocity profiles. This method may be used for LWD, MWD, or wireline logging operations.

A method for identifying a material behind a pipe string. The method may comprise disposing an acoustic logging tool into a wellbore, insonifying a pipe string within the wellbore with the acoustic logging tool, recording sonic or ultrasonic data. The method may further comprise inputting the sonic or ultrasonic data into trained a machine learning model and identifying the material behind the pipe string using the trained machine learning model.

The present invention belongs to the field of treatment for data identification and recording carriers, and specifically relates to a method and system for evaluating the filling characteristics of a deep paleokarst reservoir through well-to-seismic integration, which aims to solve the problems that by adopting the existing petroleum exploration technology, the reservoir with fast lateral change cannot be predicted, and the development characteristics of a carbonate cave type reservoir in a large-scale complex basin cannot be identified. The method comprises: acquiring data of standardized logging curves; obtaining a high-precision 3D seismic amplitude data body by mixed-phase wavelet estimation and maximum posteriori deconvolution and enhancing diffusion filtering. According to the method and the system, the effect of identifying the development characteristics of the carbonate karst cave type reservoir in the large-scale complex basin can be achieved, and the characterization precision is improved.

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 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.