Provided are a device, method, and computer-readable recording medium for detecting an earthquake in a microelectromechanical system (MEMS)-based auxiliary seismic observation network. The method includes performing detrending of removing a moving average from original acceleration data received from single sensors of an MEMS-based auxiliary seismic observation network to preprocess the acceleration data, calculating a short-term average/long-term average (STA/LTA) value using a filter parameter value specified on the basis of the preprocessed acceleration data, generating an event occurrence message or event end message on the basis of the calculated STA/LTA value and transmitting the event occurrence message or event end message, when the event occurrence message is generated, calculating an earthquake probability through an earthquake detection deep learning model using the preprocessed acceleration data as an input, and analyzing noise by calculating a power spectral density (PSD) from the original acceleration data which is merged at certain intervals.

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

A method for weathered layer correction of seismic data includes identifying arrival times in the seismic data corresponding to a weathered layer velocity gradient. A velocity model of the weathered layer is generated using the arrival times. The seismic data are time adjusted using the velocity model.

Methods of and systems for forming an image of a subterranean region of interest are disclosed. The method includes obtaining 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 source and receiver geometry of the observed seismic dataset. The method also includes forming a plurality of time-windowed trace pairs from the simulated and the observed seismic datasets, and forming an objective function based on a penalty function and a cross-correlation between the members of each pair. The method further includes determining a seismic velocity increment based on the extremum of the objective function and forming an updated seismic velocity model by combining the seismic velocity increment and the seismic velocity model, and forming the image of the subterranean region of interest based on the updated seismic velocity model.

Geologic modeling methods and systems disclosed herein employ an improved simulation meshing technique. One or more illustrative geologic modeling methods may comprise: obtaining a geologic model representing a faulted subsurface region in physical space; providing a set of background cells that encompass one or more partial faults within the subsurface region; defining a pseudo-extension from each unterminated edge of said one or more partial faults to a boundary of a corresponding background cell in said set; using the pseudo-extensions and the background cell boundaries to partition the subsurface region into sub-regions; deriving a simulation mesh in each sub-region based on the horizons in each sub-region; and outputting the simulation mesh.

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.

Systems and methods for simulating subterranean regions having multi-scale fracture geometries. Non-intrusive embedded discrete fracture modeling formulations are applied in conjunction with commercial simulators to efficiently and accurately model subsurface transient flow characteristics in regions having complex hydraulic fractures, complex natural fractures, or a combination of both, and geometries including corner point grids.

Method for characterizing subterranean formation is described. One method involves simulating a poroelastic pressure response of known fracture geometry utilizing a geomechanical model to generate a simulated poroelastic pressure response. Compiling a database of simulated poroelastic pressure responses. Measuring a poroelastic pressure response of the subterranean formation during a hydraulic fracturing operation to generate a measured poroelastic pressure response. Identifying a closest simulated poroelastic pressure response in the library of simulated poroelastic pressure response. Estimating a geometrical parameter of a fracture or fractures in the subterranean formation based on the closest simulated poroelastic pressure response.

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.