55 results about "Earthquake map" patented technology

There is provided herein a system and method of acquiring, processing, and imaging transient Controlled Source ElectroMagnetic (t-CSEM) data in ways that are similar to those used for seismic data. In particular, the instant invention exploits the time-distance characteristics of t-CSEM data to permit the design and execution of t-CSEM surveys for optimal subsequent processing and imaging. The instant invention illustrates how to correct t-CSEM data traces for attenuation and dispersion, so that their characteristics are more like those of seismic data and can be processed using algorithms familiar to the seismic processor. The resulting t-CSEM images, particularly if combined with corresponding seismic images, may be used to infer the location of hydrocarbon reservoirs.

The invention discloses an earthquake image structure guiding noise reduction method based on regularization mixed norm filtering. The earthquake image structure guiding noise reduction method includes the following steps that a gradient structure tensor is solved for an input three-dimensional earthquake image; regularization mixed norm filtering is conducted on the gradient structure tensor; a diffusion tensor is designed according to the eigenvalue and eigenvector of the filtered gradient structure tensor; continuity factors are calculated, the continuity factors at the position of a boundary fault feather edge and the like are close to zero, and the maintain performance of the structure is achieved; a sobel operator serves as a derivation operator so that divergence can be calculated. By means of the earthquake image structure guiding noise reduction method based on regularization mixed norm filtering, the textured edge information of earthquake-related data can be reserved, Gaussian noise, ultra Gaussian noise and sub Gaussian noise can be effectively suppressed, and therefore an efficient noise reduction method is achieved.

Reflected-transmitted non-converted and converted seismic waves as used for seismic imaging which is not limited by the angle of inclination of the seismic interfaces. These wave signals are recorded at the surface and are later used to obtain seismic imaging of geological discontinuities and the physical properties of the medium. Known duplex waves are used with transmitted waves for obtaining seismic images from much wider angles, ranging from horizontal to vertical. Two types of migration procedures are used: the first one being based on migration of transmitted compressional and converted duplex waves, and the second one being based on the joint migration of two interrelated fields, such as interferometric principles; for example, the first field of primary waves corresponding to the one time reflected monotypic compressional (non-converted) type of waves and a second wavefield corresponding to secondary waves of more complex origin, such as reflected waves and converted waves.

A method for computing amplitude independent gradient for seismic velocity inversion in a frequency domain includes receiving seismic data associated with a region. The region comprises one or more earth subsurface layers represented by a plurality of points, and each point is associated with a seismic velocity. Seismic velocities at the plurality of points are determined by iteratively updating the seismic velocities based on a plurality of gradient values, where each gradient value corresponds to a point and is determined by evaluating a gradient of an objective function at a location of thepoint. A seismic image of the one or more earth subsurface layers is displayed based on the determined seismic velocities.

Automatic propagation of real-world parent seismic images to efficiently generate a collection of realistic synthetic child training images to train a model for accurate automatic seismic interpretation. A 3D structural model in a present-day geological space (e.g., G_B) depicting subsurface locations of particles (e.g., in region B) may be transformed by a 3D coordinate space transformation (e.g., uvt_B) to a depositional space (e.g., G*_B) depicting past depositional locations of those particles (e.g., corresponding depositional region B). A real-world parent image depicting subsurface locations of particles (e.g., in region A) may be transformed, via a forward transformation (e.g., uvt_A), to a depositional seismic image in the depositional space of the three-dimensional structural model (e.g., G*_A=G*_B). A reverse transformation (e.g., uut_B⁻¹) may transform the depositional seismic image from the depositional space into synthetic child training images in the present-day geological space (e.g., G_B) for training the model.

The invention discloses a rigidity matrix calculation method for point source seismic response for seismic prediction, and the method comprises the steps: supposing that a point source dislocation layer is fixed, and solving a dynamic response and a fixed end face counterforce in the dislocation source layer; reversely applying counter-force of the fixed end face to the layered half space to obtain dynamic response of the fixed end face; superposing the dynamic response in the dislocation source layer and the dynamic response of the fixed end surface to obtain a total dynamic response; and combining the total dynamic response with a seismic source time function, carrying out Fourier-Hankel inverse transformation, and obtaining a seismic map of three components in a time-space domain. The method is a semi-analytical method based on combination of an accurate dynamic stiffness matrix and integral transformation, and has the advantages that the precision is not limited by the layer thickness, and error accumulation does not exist.

The invention discloses a seismic image completion method based on an information entropy norm. The method is applied to the field of geophysical seismic data processing. The problem that a low-rank matrix obtained through recovery based on a matrix rank reduction method is not a real low rank due to the limitation of an existing L1 norm is solved. According to the method, the objective function based on the information entropy norm replaces a nuclear norm adopted by original SVD, then the information entropy norm is used for regularizing a rank minimization item of the objective function, andfinally an iterative ADMM estimation method is used for approximate solution to obtain a low-rank matrix, so that steady recovery of a missing image is obtained.

The invention discloses a seismic image completion method of a fast double-interpolation POCS applied to the field of seismic data processing, and in order to solve the problem of low calculation efficiency in the prior art, a specific linear combination of first two iterative process results is used as an iterative contraction operator, and then double interpolation is carried out, so that complete seismic data is recovered; compared with a traditional POCS method, the method provided by the invention has a faster convergence speed, and the high efficiency and effectiveness of the fast dual-interpolation algorithm are verified on real seismic data.

The invention discloses a lithofacies threshold-based mixed-order anisotropic diffusion seismic image denoising method, and the method combines a low-order anisotropic diffusion equation and a high-order anisotropic diffusion equation, estimates a diffusion threshold based on gradient statistical information of stratum lithofacies, and controls weights of two anisotropic diffusion items through the lithofacies threshold. According to the method, the defect that the boundary is excessively smooth or a local block model is generated when the seismic image is denoised through a traditional anisotropic diffusion method is overcome, noise interference of the seismic image can be effectively removed through the method, and structural detail information can be recovered while the boundary of theimage is protected.

The invention provides a seismic imaging method, which belongs to the technical field of geophysical exploration. The seismic imaging method comprises the steps of acquiring time domain seismic post-stack data of a to-be-studied area and depth domain well data of a plurality of test wells in the to-be-studied area; establishing a structural frame model through time domain seismic post-stack data;preprocessing each piece of depth domain well data to convert the depth domain well data into time domain well data; correspondingly interpolating each piece of time domain well data into the construction framework model, and calculating to obtain an initial speed model; and optimizing and updating the initial velocity model, and performing pre-stack depth migration processing on the pre-stack seismic data to obtain a seismic image. By means of the seismic imaging method, the imaging precision can be higher.

A method is described for seismic facies identification including receiving a seismic dataset representative of a subsurface volume of interest; performing a machine learning algorithm on the seismicdataset to identify seismic facies and generate a classified seismic image; and identifying geologic features based on the classified seismic image. The method may be executed by a computer system.

The invention provides a method and a system for obtaining seismic parameter estimation, and the method comprises the steps: carrying out the compressed sensing of a sparse seismic signal in seismic map data, and constructing a seismic map database; performing space division on the complete signal of the seismic map database, and establishing a first index; performing compressed sensing on sparseseismic signals in the seismic wave travel time information, and constructing a seismic wave travel time table; performing space division on the complete signal of the seismic travel time table, and establishing a second index; inputting seismic event information, and performing parallel search in the seismic map database and the seismic travel time table through the first index and the second index to obtain a search result set; and determining matched seismic map data and seismic wave travel time information from the search result set, and obtaining seismic parameter estimation of the seismic event from the matched seismic map data and seismic wave travel time information. According to the method, the seismic map database and the seismic wave travel time table are constructed by adoptinga compressed sensing technology, the influence of model difference on a search result can be reduced, and the estimation of seismic parameters is obtained.

The invention provides an adversarial neural network high-resolution seismic fault detection method and system, and the method comprises the steps: training a target adversarial neural network based on a preset training sample set, and obtaining a trained target adversarial neural network, wherein the preset training sample set comprises seismic data and fault labels; the target adversarial neural network comprises a segmentation module, a feature fusion module and a discriminator module; the segmentation module is a module for obtaining fault features based on a preset training sample set; the feature fusion module is used for fusing fault features and seismic data into a global feature map, and performing seismic fault detection on the target seismic image based on the trained target adversarial neural network. The technical problem that in the prior art, seismic fault detection precision and resolution are low is solved.

The invention relates to an observation signal amplitude measurement method based on natural seismic source radiation intensity correction, and aims to solve the problem that the amplitude measurement is difficult to perform by using a recording signal of a large-magnitude natural earthquake on a seismic station due to the influence of the intensity directivity of a seismic source, and a seismic signal amplitude-magnitude relation model of a specific path from a seismic source to a seismic station is difficult to establish. The method is based on a known seismic source mechanism of a natural earthquake, a theoretical seismic map is generated by using a theoretical seismic map calculation method, a theoretical seismic map waveform containing a plurality of path arrival signals is obtained, and a plurality of path arrival signals are obtained by measuring an average value of the theoretical seismic maps of different azimuth angles and measuring a theoretical seismic map signal amplitude on a seismic station. And the average value is normalized to obtain an amplitude correction coefficient corresponding to the seismic station so as to correct the radiation intensity of the seismic source.

The invention discloses a seismic data interpolation method combining Gabor feature extraction and support vector regression. The method comprises the steps of pre-interpolation, Gabor feature extraction, feature vector extraction, regression function training and seismic image reconstruction which are executed in sequence. The method comprises the following steps: performing pre-interpolation complementation on a missing seismic image lacking a seismic trace to obtain a low-resolution seismic image, performing Gabor filtering on the low-resolution seismic image to obtain a feature image, obtaining a prediction feature vector by utilizing the low-resolution seismic image and feature image transformation, and training a regression function by training the feature vector and a label; and finally, inputting the prediction feature vector into a trained regression function, and obtaining a reconstructed seismic image through regression reconstruction and transformation. According to the invention, the seismic data interpolation algorithm combining Gabor filtering and SVR is designed, the feature image obtained by Gabor feature extraction and the regression reconstruction capability of SVR are fully utilized, and a clear and complete seismic image is obtained.

The invention discloses a seismic data orderliness assessment method and device. The method comprises the following steps: acquiring a three-dimensional seismic data volume to be assessed; performing point-by-point scanning on the three-dimensional seismic data according to the depth direction, the main survey line direction and the connection survey line direction, and calculating the image entropy of each sampling point to obtain an entropy data volume of the three-dimensional seismic data volume; and according to the entropy data volume, performing order evaluation on the three-dimensional seismic data volume. According to the method, the reliability of the seismic data generating the fracture attribute is evaluated based on the seismic image entropy, and then the seismic data is optimized by using the evaluation result, so that the purposes of adjusting the fracture attribute and balancing the fracture attribute are achieved, and the response form of the fracture attribute is optimized.

The embodiment of the invention provides a time-lapse seismic observation data analysis method and device, and the method comprises the steps: detecting the site conditions of an observation site, anddetermining the arrangement types of an excitation point and a receiving point according to the site conditions of the observation site; performing a site experiment to obtain basic parameters; whenthe arrangement type is a linear type, arranging the positions of the excitation point and the receiving point according to the basic parameters; when the arrangement type is a regular grid type, calculating observation parameters according to the basic parameters, and arranging the positions of the excitation point and the receiving point according to the observation parameters; selecting different time points to excite the excitation point, and obtaining seismic waves received by the receiving point; calculating seismic images of the observation site at different times according to the seismic waves received by the receiving point; and acquiring site data, establishing a corresponding relationship according to the site data and the seismic images, and performing data analysis through thecorresponding relationship. By the adoption of the method, the seismic data of an expansive soil bank slope and a dam landslide can be reflected in time and space, and then abnormity can be accurately recognized.

Disclosed is a method and associated computer program and apparatus for characterising changes within a subsurface volume between a first time and a second time. The method comprises obtaining first seismic data corresponding to the first time and processing this data to obtain a seismic image of the subsurface volume. This processing is reversed for relevant portions of the seismic image to obtain relevant portions of first seismic data. Changes within the subsurface volume between the first time and the second time are characterised by estimating the changes between second seismic data corresponding to the second time and the relevant portions of first seismic data.