34results about "Seismic signal processing" patented technology

The invention provides a theoretical model and a test model for compressional wave velocity and shear wave velocity convention and a test model for wave velocity and density conversion based on elastic wave propagation theory and basic parameters including the compressional wave velocity, the shear wave velocity and density through theoretical analysis and experimental study to quantitatively describe basic physical parameters and models of coal-measure sedimentary rocks. The invention constructs correlation between dynamic elasticity mechanics parameters and static elasticity mechanics parameters of coal-measure sedimentary rocks and between physical mechanics parameters and acoustic velocity of coal-measure rocks and models thereof based on acoustic velocity experiment and three parameters of rocks including the compressional wave velocity, the shear wave velocity and the density, by measurement of static mechanics parameters of coal-measure rocks and statistics and regression analysis of test data. The invention further provides a method for obtaining seismic wave interval velocity parameter by interval velocity inversion of wave impedance, develops rock physical mechanics parameter analytical software based on seismic information and realizes in situ physical mechanics parameter calculation of coal and rocks.

A system and method for processing data on a peripheral device that is operatively coupled to a host computing system via a peripheral bus. The compression of input data transmitted to the peripheral device and/or the size of the storage provided on the peripheral device may enhance the efficiency of the processing of the data on the peripheral device by obviating a bottleneck caused by the relatively slow transfer of data between the host computing system and the peripheral device.

The invention provides a method for optimizing a design of a vertical seismic profile observation system for seismic exploration, which comprises the following steps: establishing a seismic geological model by known data; distributing more than two different observation systems; forwarding a response of a seismic wave field; calculating a position of a reflective spot, incident and reflection angles and length or ray impact frequency, and counting characteristic parameters again to construct an analysis map; and determining the optimal scheme for actual exploration according to a geological target. The method adopts multi-parameter analysis based on a spacial model, has strong rationality of results, adapts to details, can improve data acquisition quality, is suitable for homogeneous spaces and one-dimensional, two-dimensional or three-dimensional models, can be used for non-zero spacing, Walkaway and 3D-3C VSP exploration designs, and also can be used for analyzing VSP exploration designs of longitudinal waves, transverse waves and converted waves.

Method of imaging seismic data recorded using a set of seismic receivers positioned to receive seismic responses from a subsurface formation upon activating a seismic source. The seismic responses are direction-sensitive responses. A virtual signal received by receiver k from the virtual source at the position of receiver m is determined by processing involving cross correlating at least part of the responses of a receiver m with at least part of the responses of a seismic receiver k. Part of direction-sensitive responses or part of the virtual signal that contains a component of the wave field at receiver m traveling in a first direction and part of the direction-sensitive responses or part of the virtual signal that contains a component of the wave field at the receiver k traveling in a second direction may be removed.

A sensor device, comprising two symmetrically disposed sonolucent wedges (5), and a connecting piece for fixedly connecting the two sonolucent wedges (5); the upper surfaces of the sonolucent wedges (5) are provided with inclined planes; installation holes are formed on the inclined planes; transducers (3) are installed in respective installation holes; one transducer (3) is used to generate ultrasonic waves, and the other transducer (3) is used to receive the ultrasonic waves generated by the previous transducer (3). The residual stress detection system comprises a sensor device, an ultrasonic transmission card, and a data acquisition card.

The invention discloses a method for acquiring earthquake early warning magnitude, and the objective of the invention is to provide a novel method combining two earthquake early warning magnitude acquisition means of prediction and actual measurement to overcome shortcomings of existing early warning magnitude calculation methods. The technical scheme of the invention is that the method comprises the following steps: 1 P waves and S waves of earthquake are acquired by using a plurality of earthquake wave acquisition stations which are close to the epicenter; 2 feature information of first n seconds of the earthquake waves acquired by each of the earthquake wave acquisition stations after the earthquake is extracted, and earthquake early warning magnitude of each of the stations is calculated by using an earthquake magnitude prediction formula; 3 earthquake magnitude of each of the stations is calculated by using an actual measurement formula; 4 if the earthquake magnitude of one station calculated by using the earthquake magnitude actual measurement formula is greater than the earthquake early warning magnitude calculated by using the earthquake magnitude prediction formula or the end of the earthquake waves is monitored by the station, then the earthquake magnitude calculated by using the earthquake magnitude actual measurement formula is selected as the earthquake early warning magnitude for the station; otherwise the predicted earthquake magnitude is selected as the early warning magnitude; 5 the average value of the earthquake early warning magnitude of each of the earthquake wave acquisition stations is selected as comprehensive earthquake early warning magnitude of the earthquake early warning.

The invention relates to a prediction method for fractured oil-gas reservoirs based on three-dimensional seismic data and a fractional-frequency ant body technology. The method comprises the following steps: transforming the seismic data to a frequency domain to obtain the scope of an effective seismic frequency band; starting from the minimum value of the effective seismic frequency band to increase at intervals to obtain a fractional-frequency data series; computing the phase spectrum of each sampling point to obtain a discrete fractional-frequency phase data volume; adopting an ant body algorithm to generate a series of fractional-frequency ant body data; conducting horizon picking on an objective reservoir to obtain the time structure data of each time point location of seismic reflected waves; respectively extracting numerical values corresponding to an objective interval from a series of fractional-frequency ant body data to obtain a group of fractured reservoir detection results of different frequency bands of the objective interval. The computing method is simple and convenient and can be used for detecting smaller scale of faults and fracture zones, the precision is higher, and the geological effect is obvious.

The invention provides an elastic wave vector imaging method and device, a storage medium and an apparatus. The method comprises the steps of building a conversion relation between group angles and polarizing angles by using the P wave speed of a set direction of an anisotropic medium of each point in an imaging area, the S wave speed of each set direction and anisotropic parameters; using an elastic wave equation to perform wave field depth prolongation of multi-component earthquake data in the imaging area to obtain a focus forward propagation wave field and a shot record back propagation wave field and corresponding Poynting vectors; according to set time intervals, decomposing the focus forward propagation wave field and the shot record back propagation wave field by using the Poyntingvectors and the conversion relation between group angles and polarizing angles; performing vector imaging of elastic waves in anisotropic media by using wave field decomposition results correspondingto all shots in the multi-component earthquake data. The imaging efficiency and accuracy can be improved.

The invention discloses a well-ground seismic data joint acquisition system based on distributed optical fiber acoustic wave sensing and a well drive data processing method. In order to solve the problems of inconsistent energy, inconsistent frequency spectrum and incompletely consistent coupling between a seismic source and the ground caused by repeated excitation of each seismic source point when seismic data in a well and seismic data on the ground are respectively acquired in the prior art, the invention discloses an in-well-ground seismic data combined three-dimensional acquisition systemformed by a distributed optical fiber sound wave sensing seismic data acquisition unit and a ground seismic data acquisition unit. In-well-ground combined three-dimensional exploration is carried out, ground and in-well seismic data are synchronously collected, the high-density, high-benefit, high-resolution and low-cost well-ground combined three-dimensional seismic exploration technology is achieved, and oil and gas resource exploration and comprehensive evaluation are carried out.

Methods and systems for acoustically determining reservoir parameters of subterranean formations. A tool comprising a plurality of acoustic sources and configured for acoustic measurements is deployed within a wellhole. Acquired acoustic data are processed and utilized for deriving key parameters for the formations.

The invention provides a carbonatite reservoir stratum inversion method and system with transverse sedimentary argillaceous influence elimination. The method comprises: according to processed base data, forward rock physical modeling is carried out to obtain a forward modeling result including physical elastic rock characteristics of sedimentary argillaceous limestone, limestone cave filling clay, pure limestone, and the reservoir stratum; on the basis of the physical elastic rock characteristics of the sedimentary argillaceous limestone, a correspondence relation between the carbonatite sedimentary argillaceous characteristics curve and elastic impedance is constructed; according to the correspondence relations, corresponding to wells in a research area, between the carbonatite sedimentary argillaceous characteristics curves and elastic impedances, a carbonatite sedimentary argillaceous background impedance body is generated; geostatistics inversion is carried out based on the processed base data to obtain a carbonatite reservoir stratum impedance body; and the carbonatite reservoir stratum impedance body and the carbonatite sedimentary argillaceous background impedance body are coupled to obtain a carbonatite reservoir stratum inversion result with transverse sedimentary argillaceous changing influence elimination. Therefore, the predicted reservoir stratum and the actual reservoir stratum match each other; and the precision is high.

The invention provides a common detection point data gather velocity analysis method. The method comprises the steps that 1 a common detection point data gather to be subjected to velocity analysis is transformed to acquire a pseudo-co-center point data gather; and 2 velocity analysis is carried out on the pseudo-co-center point data gather through a co-center point velocity analysis technology to acquire superposition velocity based on the common detection point data gather. According to the invention, the method is used for common detection point data gather velocity analysis in multi-wave seismic data processing to acquire accurate superposition velocity based on the common detection point data gather, and an accurate common detection point data gather superposition result is provided for acquiring accurate conversion wave remaining static correction amount.

Provided is a self-adaptive realization method of a discontinuous Galerkin finite element earthquake value simulation algorithm. The method includes the following steps: firstly, defining a mesh generation density and conducting non-structural mesh generation on a calculation area according to the mesh density; secondly, determining a polynomial expansion order p of each mesh unit; thirdly, decoupling two adjacent mesh units, and selecting an appropriate numerical flux format as a wave field exchanged mode of two adjacent mesh units; fourthly, performing numeralization on a decoupling equation, and obtaining a space discrete equation of each mesh unit; fifthly, solving an inner product among wave field functions on a common boundary of two adjacent units, and obtaining a boundary value mapping matrix; sixthly, selecting an appropriate time integral format; and seventhly, utilizing the boundary value mapping matrix to carry out mapping on adjacent unit wave fields different in order, calculating space discrete portions, applying hypocenter conditions, updating wave field components, and obtaining a simulation result.

A method for identifying seismic horizons in digitized seismic images includes deploying a swarm of agents at an initial position in a seismic image to be analyzed, where the swarm of agents includes picking agents that define a direction for the swarm and averaging agents that smooth the direction of the swarm, identifying a direction to follow in the seismic image for each picking agent, and advancing each picking agent in the identified direction, and averaging, by the averaging agents, the directions identified by the picking agents, wherein if an information concentration measured by a picking agent at a current time step is greater than a previous time step, the picking agent keeps a previous direction, otherwise the picking agent changes direction according to the average current state of the set of averaging agents within its neighborhood.

A method for acquiring and improving seismic data includes activating a first seismic source located below a geophysical surface at a first depth and a second seismic source located below the geophysical surface at a second depth, wherein the second depth is below the first depth, acquiring seismic data with a seismic receiver in conjunction with activating the first seismic source and the second seismic source, and aligning primary reflections within the seismic data to provide improved seismic data. The method may also include determining changes to the regions below the second depth by comparing improved seismic data corresponding to a first acquisition event with improved seismic data corresponding to a second acquisition event. A corresponding system is also disclosed herein.

The invention provides a delta facies sand body spatial form characterization method which comprises the following steps: S1, preprocessing observed seismic data through a variable coefficient guideddirection optimization spatial filtering technology to improve a signal-to-noise ratio; S2, under the guidance of the longitudinal superposition characteristics of the sand body gene units identifiedby the logging data, acquiring position characteristic parameters required for characterizing the spatial form of the sand body from the preprocessed seismic data; and S3, carrying out spatial direction sensitivity research on the sand body posture characteristic parameters in inclination angles, azimuth angles and time directions to obtain optimized parameter values of equalization distribution processing. The challenge faced by the conventional method for characterizing the gene unit in the delta facies sand body is solved, and the increasingly increased high-precision characterization requirement of the reservoir in the oil and gas field development is met.

The invention discloses an RMS-SVD multiple suppression method and device, electronic equipment and a medium. The method comprises the following steps: picking up the superposition velocity of multiples, and further obtaining the root-mean-square velocity of the multiples; performing dynamic correction processing on the target data gather through the multiple wave root mean square velocity to obtain a dynamically corrected target data gather; determining the time window range of multiple waves, and further determining the time window length of SVD filtering; sVD filtering is carried out on the target data gather after dynamic correction, and a filtered target data gather is obtained; and carrying out reverse dynamic correction processing on the filtered target data gather to obtain a multiple suppressed target data gather. According to the method, the multiple energy on the target data gather can be effectively suppressed and eliminated through the difference between the speed of the multiple and the speed of the primary reflected wave, and the signal-to-noise ratio of the target data gather is improved.

The invention discloses a method for highlighting small-scale fracture-cavity information under a strong seismic reflection interface. In the method of the present invention, for the seismic data containing strong reflection interface interference, under the condition of ensuring the stability of the wavelet, the seismic data volume highlighting the small-scale fracture-cavity information is generated by reconstructing the reflection coefficient. According to the method of the present invention, the masking of reservoir seismic information by the strong reflection interface can be effectively weakened, and the seismic response characteristics of the reservoir can be highlighted, thereby improving reservoir prediction accuracy and drilling success rate, and providing powerful data support for guiding drilling deployment.

The invention relates to the technical field of oil and gas seismic exploration, in particular to a method for realizing pre-stack three-dimensional reverse time migration by utilizing a checkpoint technology based on an FPGA. The method comprises the steps of setting checkpoint parameters on a CPU platform through software design according to 3D-RTM calculation tasks with different scales and the resource quantity of a selected FPGA chip; and enabling the wave field data calculation time and the required data transmission time to be dynamic so as to achieve the effect of eliminating the I/O bottleneck. The resource utilization rate of the FPGA is improved, the calculation time of the whole calculation task is greatly shortened, and then the acceleration efficiency of the whole heterogeneous calculation is improved. According to the invention, the optimal FPGA and CPU data interaction scheme is set according to the characteristics of data volumes with different scales, so that the whole computing resource utilization rate is optimal. Meanwhile, dynamic adjustment parameters are realized at a CPU platform end through a software architecture, so that the expansibility of the method is greatly enhanced, the utilization rate of computing resources of any FPGA selected for 3D-RTM computing tasks with different scales is improved, and higher compatibility is achieved.

The invention provides an OVT gather processing method and device for PS wave seismic trace data and electronic equipment. The method comprises the steps that PS wave seismic trace data are mapped on a surface element grid of the CCP, a projection point of an imaging point of the PS wave seismic trace data relative to the earth surface is used as a center, a virtual shot point of the surface element grid relative to the PS wave seismic trace data is determined, and the virtual shot point and an original shot point of the surface element grid relative to the PS wave seismic trace data have the same shot-geophone distance. And based on the travel time of the PS wave seismic trace data at the virtual shot point reflected by the surface element grid, correcting the time distance curve to hyperbolic time distance relation of the PS wave seismic trace data, so that the PS wave seismic trace data is positioned on the common center point of the corrected seismic trace data relative to the conversion point reflected by the surface element grid. PS wave seismic trace data is mapped on the offset line grid of the OVT, OVT gather extraction is carried out on the PS wave seismic trace data, an OVT gather of the PS wave seismic trace data is obtained, and extraction information comprises offset, azimuth angle, coordinates of CMP points and travel time.

The invention discloses a first arrival chromatography near-surface modeling effective model depth control method and system. The method comprises the following steps: obtaining ray maximum depth andoptimization speed models of the first iteration and the second iteration; setting a threshold value, and carrying out the Nth iteration: 1, carrying out the preprocessing of the maximum depth of theray of the (N-2) th iteration and the maximum depth of the ray of the (N-1) th iteration; 2, calculating a depth difference; 3, calculating the effective model depth of the Nth iteration; 4, obtainingan optimization speed model; 5, whether N is equal to a threshold value or not is judged, if yes, a near-surface velocity model is obtained, and if not, the step 6 is executed; and 6, optimizing whether the speed model reaches a preset target or not, if so, obtaining a near-surface speed model, and if not, carrying out (N + 1) th iteration. According to the method, the effective model depth of the iteration is determined through the maximum ray depth of the latest secondary iteration in the first arrival chromatography iteration process, and the speed units outside the effective model depth are excluded, so that the forward modeling efficiency is improved.

The invention relates to a method for updating the stratum seismic velocity. The method comprises steps that the seismic velocity of a drilled area is updated according to seismic time horizon interpretation result data and the real depth of each horizon, and an initial seismic velocity model of the drilled area is obtained; pre-stack depth migration is carried out by utilizing the initial seismicvelocity model of the target area to obtain a common imaging point gather of the target area; based on the common imaging point gather, reflection points with different depths on a channel corresponding to a well track of the undrilled area are selected; for each reflection point, imaging depths corresponding to different reflection angles and ray lengths of corresponding reflection rays in eachhorizon are respectively obtained, and the initial seismic velocity of the stratum of the corresponding undrilled area is updated to obtain the updated seismic velocity of the undrilled area. The method is advantaged in that the number of parameters of traditional seismic tomography can be greatly reduced, so speed modeling efficiency is improved, and a foundation is laid for correcting the geomechanical model of the stratum to be drilled in real time.