CN112444869B - Seismic data processing method and storage medium for suppressing exogenous interference waves - Google Patents
Seismic data processing method and storage medium for suppressing exogenous interference waves Download PDFInfo
- Publication number
- CN112444869B CN112444869B CN201910817071.3A CN201910817071A CN112444869B CN 112444869 B CN112444869 B CN 112444869B CN 201910817071 A CN201910817071 A CN 201910817071A CN 112444869 B CN112444869 B CN 112444869B
- Authority
- CN
- China
- Prior art keywords
- seismic data
- interference wave
- wave
- exogenous
- coordinates
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 24
- 238000007781 pre-processing Methods 0.000 claims abstract description 11
- 238000000034 method Methods 0.000 claims description 43
- 238000001914 filtration Methods 0.000 claims description 14
- 238000001514 detection method Methods 0.000 claims description 11
- 238000005516 engineering process Methods 0.000 claims description 9
- 230000003068 static effect Effects 0.000 claims description 5
- 230000002452 interceptive effect Effects 0.000 claims 4
- 230000001629 suppression Effects 0.000 abstract description 3
- 230000005284 excitation Effects 0.000 description 10
- 230000000694 effects Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a seismic data processing method and a storage medium for suppressing exogenous interference waves. A seismic data processing method for suppressing an exogenous interference wave, comprising: inputting initial seismic data; preprocessing the initial seismic data to obtain the seismic data to be processed; acquiring coordinates of at least four points of a homophase axis of the exogenous interference wave and corresponding travel time; constructing an equation set by utilizing coordinates of at least four points and corresponding travel time based on a propagation equation of the exogenous interference wave; obtaining coordinates of the external interference wave and propagation speed of the external interference wave by solving the equation set formed by at least four equations of at least four points; and processing the seismic data to be processed according to the coordinates of the exogenous interference wave and the propagation speed of the exogenous interference wave so as to suppress the exogenous interference wave. According to the seismic data processing method and the storage medium for suppressing the exogenous interference wave, which are disclosed by the embodiment of the invention, the source point position of the exogenous interference wave can be automatically obtained, and the suppression effect of the exogenous interference wave can be better improved.
Description
Technical Field
The invention belongs to the technical field of geophysical exploration seismic data processing, and particularly relates to a seismic data processing method and a storage medium for suppressing exogenous interference waves.
Background
During the seismic exploration data acquisition construction, other activities which cannot stop the operation, such as an oil well being mined, a drilling operation being drilled, a windmill running normally and the like, exist in a construction area, and vibration caused by the activities is recorded in the seismic exploration data and becomes an exogenous interference wave in the seismic exploration data, so that the seismic exploration data can be better analyzed by pressing in the seismic exploration data processing. Accurate source point positions of the exogenous interference waves are needed for suppressing the exogenous interference waves, so that calculating the accurate source point positions of the exogenous interference waves is a key problem for suppressing the exogenous interference waves in seismic exploration data processing.
At present, the method for determining the source point position of the external source interference wave is an observation positioning method, namely, when the arrival of the external source interference wave is observed through the browsing of seismic data, the seismic channel closest to the source point position of the external source interference wave is found, the coordinates of the seismic channel are read, and then the coordinates of the source point of the external source interference wave are estimated approximately. The method requires a data processor to have stronger professional technical capability, the process is time-consuming and labor-consuming, the source point position of the exogenous interference wave can only be positioned between two adjacent receiving channels and two adjacent receiving lines, and for two-dimensional seismic data or three-dimensional seismic data with larger line distance, the positioning error of the method in the vertical line direction is larger, so that the accuracy of the source point position of the exogenous interference wave is influenced, and the suppression effect of the exogenous interference wave is further influenced.
Therefore, there is a need to provide a simple and accurate seismic data processing method and storage medium for suppressing an external disturbance wave.
Disclosure of Invention
One of the technical problems to be solved by the invention is how to provide a simple and accurate seismic data processing method and a storage medium for suppressing exogenous interference waves.
In order to solve the above technical problems, an embodiment of the present application firstly provides a seismic data processing method for suppressing an external disturbance wave, including:
inputting initial seismic data;
Loading coordinate information on the initial seismic data and preprocessing the initial seismic data to obtain seismic data to be processed;
Aiming at the seismic data to be processed, acquiring each point coordinate (R x,Ry) of at least four points of an in-phase axis of an external interference wave and travel time T corresponding to each point coordinate;
Propagation equation based on exogenous interference wave Constructing an equation set by utilizing each point coordinate (R x,Ry) of at least four points and travel time T corresponding to each point coordinate, wherein v is the propagation speed of the exogenous interference wave, (P x,Py) is the coordinate of the exogenous interference wave, and T 0 is the vibration excitation time of the exogenous interference wave;
Solving the equation set to obtain coordinates of the exogenous interference wave and propagation speed of the exogenous interference wave;
and processing the seismic data to be processed according to the coordinates of the exogenous interference wave and the propagation speed of the exogenous interference wave so as to suppress the exogenous interference wave.
Preferably, loading the initial seismic data with coordinate information includes:
and loading the coordinate information of the wave detection points on the initial seismic data.
Preferably, preprocessing the initial seismic data includes:
And carrying out static correction, filtering and gain on the initial seismic data.
Preferably, the at least four points are located on the same phase axis.
Preferably, the at least four points are located on the same phase axis of the first arrival time.
Preferably, solving the equation set to obtain coordinates of the external interference wave and propagation speed of the external interference wave includes:
And solving the equation set by using a least square method to obtain the coordinates of the external disturbance wave and the propagation speed of the external disturbance wave.
Preferably, processing the seismic data to be processed according to the coordinates of the external disturbance wave and the propagation velocity of the external disturbance wave to suppress the external disturbance wave includes:
And processing the seismic data to be processed by utilizing an f-kx-ky cone filtering technology according to the coordinates of the exogenous interference wave and the propagation speed of the exogenous interference wave so as to suppress the exogenous interference wave.
Preferably, the seismic data to be processed are processed in a cross arrangement domain by utilizing an f-kx-ky cone filtering technology according to the coordinates of the exogenous interference wave and the propagation speed of the exogenous interference wave so as to suppress the exogenous interference wave.
Preferably, the method further comprises: and outputting target seismic data after suppressing the exogenous interference.
The embodiment of the invention also discloses a storage medium which stores a program for executing the method of any one of the above.
One or more embodiments of the above-described solution may have the following advantages or benefits compared to the prior art:
A seismic data processing method for suppressing an exogenous interference wave, comprising: inputting initial seismic data; the coordinate information is loaded on the initial seismic data, and preprocessing is carried out to obtain the seismic data to be processed, wherein the signal-to-noise ratio of the seismic data is improved relative to that of the initial seismic data, and the exogenous interference waves can be identified more easily; acquiring coordinates (Rx, ry) of at least four points of an in-phase axis of an external interference wave and travel time T corresponding to the coordinates of the at least four points aiming at the seismic data to be processed; propagation equation based on exogenous interference wave Constructing an equation set by utilizing the coordinates of each point of at least four points and travel time corresponding to the coordinates of each point, wherein the coordinates of an external interference wave, the propagation speed of the external interference wave and the vibration excitation time of the external interference wave in the constructed equation set are unknown, v is the propagation speed of the external interference wave, (P x,Py) is the coordinates of the external interference wave, and T 0 is the vibration excitation time of the external interference wave; obtaining coordinates of the external interference wave and propagation speed of the external interference wave by solving the equation set formed by at least four equations of at least four points; and processing the seismic data to be processed according to the coordinates of the exogenous interference wave and the propagation speed of the exogenous interference wave so as to suppress the exogenous interference wave. According to the seismic data processing method for suppressing the exogenous interference wave, disclosed by the application, a data processor is not required to estimate the source point position of the exogenous interference wave through professional technical capability, the source point position of the exogenous interference wave can be automatically obtained, the accuracy of the source point coordinates of the exogenous interference wave is high, the method is suitable for various types of exogenous interference waves, the propagation speed of the exogenous interference wave can be obtained while the source point coordinates of the exogenous interference wave are obtained, the suppressing effect of the exogenous interference wave can be better improved, and the quality of seismic data is improved.
Additional advantages, objects, and features of the invention will be set forth in part in the description which follows and in part will become apparent to those having ordinary skill in the art upon examination of the following or may be learned from practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
The accompanying drawings are included to provide a further understanding of the technical aspects or prior art of the present application, and are incorporated in and constitute a part of this specification. The drawings, which are used to illustrate the technical scheme of the present application, are not limited to the technical scheme of the present application.
FIG. 1 is a flow chart of a seismic data processing method for suppressing an external disturbance wave according to a first embodiment of the invention;
FIG. 2 is a schematic diagram of the propagation of an external disturbance wave in a plane;
FIG. 3 is a seismic monopulse recording containing an exogenous interference wave.
Detailed Description
The following will describe embodiments of the present application in detail with reference to the drawings and examples, thereby solving the technical problems by applying technical means to the present application, and realizing the corresponding technical effects can be fully understood and implemented accordingly. The embodiment of the application and the characteristics in the embodiment can be mutually combined on the premise of no conflict, and the formed technical scheme is within the protection scope of the application.
Example 1
A seismic data processing method for suppressing an exogenous interference wave, comprising: inputting initial seismic data; the coordinate information is loaded on the initial seismic data, and preprocessing is carried out to obtain the seismic data to be processed, wherein the signal-to-noise ratio of the seismic data is improved relative to that of the initial seismic data, and the exogenous interference waves can be identified more easily; acquiring coordinates (Rx, ry) of at least four points of an in-phase axis of an external interference wave and travel time T corresponding to the coordinates of the at least four points aiming at the seismic data to be processed; propagation equation based on exogenous interference waveConstructing an equation set by using at least four points, wherein the coordinates of the external interference wave, the propagation speed of the external interference wave and the vibration excitation time of the external interference wave in the constructed equation set are unknown, v is the propagation speed of the external interference wave, (P x,Py) is the coordinates of the external interference wave, and T 0 is the vibration excitation time of the external interference wave; obtaining coordinates of the external interference wave and propagation speed of the external interference wave by solving the equation set formed by at least four equations of at least four points; and processing the seismic data to be processed according to the coordinates of the exogenous interference wave and the propagation speed of the exogenous interference wave so as to suppress the exogenous interference wave. According to the seismic data processing method for suppressing the exogenous interference wave, disclosed by the application, a data processor is not required to estimate the source point position of the exogenous interference wave through professional technical capability, the source point position of the exogenous interference wave can be automatically obtained, the accuracy of the source point coordinates of the exogenous interference wave is high, the method is suitable for various types of exogenous interference waves, the propagation speed of the exogenous interference wave can be obtained while the source point coordinates of the exogenous interference wave are obtained, the suppressing effect of the exogenous interference wave can be better improved, and the quality of seismic data is improved.
The invention is further illustrated by the following specific examples.
As shown in fig. 1, a seismic data processing method for suppressing an external disturbance wave according to an embodiment of the present invention includes:
S1, inputting initial seismic data.
The initial seismic data may be seismic recordings obtained by a plurality of receivers deployed on the ground. As shown in fig. 2, there is an external disturbance wave P, such as a rotating windmill, within the survey area. The windmill generates vibration when rotating, the seismic wave propagates outwards from the point P, and the wave front is concentric. Five detection lines L1-L5 are distributed on the ground, each detection line is composed of a plurality of detectors with the same distance, the positions of the detectors are detection points, and the detection points are shown as inverted triangles in fig. 2. On the L3 line, the positions of R 1、R2……R3、R6……Rn are marked as the positions of the detection points, and the positions can record seismic data. The geophone can record the seismic record of the external disturbance wave P in addition to the normal seismic record.
Fig. 3 shows a seismic monopulse recording containing an exogenous interference wave, with the horizontal axis representing the position of the detector point and the vertical axis representing time. Wherein A is refractive wave, B, C is plane wave interference, D, E, F is effective reflection wave, G is a group of external interference waves which are approximately parallel from top to bottom and are hyperbolic, and the external interference waves are required to be suppressed in reflection seismic data processing. The same phase axis G shown in fig. 3 can be formed by connecting the seismic recordings of the external disturbance wave P recorded by the plurality of detectors in fig. 2.
S2, loading coordinate information on the initial seismic data and preprocessing the initial seismic data to obtain the seismic data to be processed.
Loading the initial seismic data with coordinate information, including: and loading the coordinate information of the wave detection points on the initial seismic data. The initial seismic data generally does not have the wave-detecting point coordinate information, and the initial seismic data is matched with the wave-detecting point coordinate information and can be used for calculating the coordinates and the propagation speed of the external source interference waves.
Additionally, preprocessing the initial seismic data includes: static correction, filtering and gain are carried out on the initial seismic data.
The processing procedures of static correction, filtering, gain and the like can improve the signal-to-noise ratio of the seismic data, and can more easily identify the exogenous interference waves. In particular embodiments, other pretreatment methods may be used herein to increase the signal-to-noise ratio of the seismic data.
The coordinate information is loaded on the initial seismic data, which is generally the coordinate information of the wave detection point, and the seismic data loaded with the coordinate information is further preprocessed to obtain the seismic data to be processed.
S3, acquiring coordinates (R x,Ry) of at least four points of the same phase axis of the external interference wave and travel time T corresponding to the coordinates of the at least four points according to the seismic data to be processed.
For the seismic data to be processed, the seismic data with exogenous interference waves can be found. Such as the same phase axis G shown in fig. 3. And acquiring coordinates and corresponding travel time of at least four points on the same phase axis, and constructing an equation set to solve the coordinates and the propagation speed of the exogenous interference wave. For example, coordinates of 5 points, 10 points, 20 points on the same phase axis and corresponding travel times are obtained.
In one embodiment, at least four points are located on the same phase axis. For example, coordinates of at least four points on the same hyperbola and corresponding travel times are obtained.
In one embodiment, at least four points are located on the same phase axis of the first arrival time. In seismic exploration, a seismic wave excited from a shot and arriving at a pickup first is called a first arrival wave, and this arrival time is called the first arrival time of the seismic wave. And finding out an in-phase shaft of the first arrival time of the seismic wave, and acquiring coordinates of at least four points and corresponding travel time on the in-phase shaft.
In one embodiment, at least four points may be taken on different hyperbolas. For example, a first point may be acquired on a first hyperbola, a second point may be acquired on a second hyperbola, a third point may be acquired on a third hyperbola, and a fourth point may be acquired on a fourth hyperbola, any two or any three of the first, second, third, and fourth hyperbolas may be the same hyperbola, or the first, second, third, and fourth hyperbolas may be four different hyperbolas.
S4, propagation equation based on exogenous interference waveAn equation set is constructed by using at least four points, wherein v is the propagation speed of the external interference wave, (P x,Py) is the coordinates of the external interference wave, and T 0 is the vibration excitation time of the external interference wave.
Because the unknown number is four, namely two coordinate values of the exogenous interference wave, the propagation speed of the exogenous interference wave and the vibration excitation time of the exogenous interference wave, an equation set of the propagation equations of the 4 exogenous interference waves can be constructed by utilizing four points. The coordinates of the external interference wave, the propagation speed of the external interference wave and the vibration excitation time of the external interference wave can be obtained by solving the equation set.
In order to obtain a more accurate, stable and reliable solution, an overdetermined system of propagation equations of more than 4 exogenous interference waves can be constructed by using coordinates of more than four points and corresponding travel times. For example, an overdetermined equation set consisting of 5 propagation equations, 10 propagation equations, and 20 propagation equations is constructed using coordinates of 5 points, 10 points, and 20 points, and corresponding travel times, respectively.
In one embodiment, the coordinate values of the first detector R1 (R x1,Ry1), the coordinate value of the second detector R2 (R x2,Ry2), the coordinate value of the third detector R3 (R x3,Ry3), the coordinate value of the fourth detector R4 (R x4,Ry4) and the coordinate value of the fifth detector R5 (R x5,Ry5) and their corresponding travel times T 1、T2、T3、T4、T5 are obtained on a same phase axis.
An overdetermined equation set of 5 propagation equations is constructed according to the coordinate values of the five detection points and the travel time, and the overdetermined equation set is as follows:
In the above overdetermined equation set, v is the propagation velocity of the external disturbance wave, (P x,Py) is the coordinate of the external disturbance wave, and T 0 is the vibration excitation time of the external disturbance wave. By solving the above overdetermined equation set, the propagation speed of the external interference wave and the coordinates of the external interference wave can be obtained, and the obtained solution is accurate, stable and reliable.
And S5, solving the equation set to obtain the coordinates of the external disturbance wave and the propagation speed of the external disturbance wave.
In one embodiment, solving the system of equations to obtain the coordinates of the external disturbance wave and the propagation velocity of the external disturbance wave includes:
And solving an equation set by using a least square method to obtain the coordinates of the external disturbance wave and the propagation speed of the external disturbance wave.
The least square method is a common algorithm for solving the overdetermined equation set, has wide application in engineering and calculation mathematics, and can realize the calculation of the application by utilizing various open source programs.
Of course, in other embodiments, other methods may be utilized to solve the system of overdetermined equations. The scheme of the application is not limited to the method of solving the system of overdetermined equations.
S6, processing the seismic data to be processed according to the coordinates of the exogenous interference wave and the propagation speed of the exogenous interference wave so as to suppress the exogenous interference wave.
In one embodiment, processing the seismic data to be processed to suppress the exogenous interference wave according to the coordinates of the exogenous interference wave and the propagation velocity of the exogenous interference wave includes:
And processing the seismic data to be processed by using an f-kx-ky cone filtering technology according to the coordinates of the exogenous interference wave and the propagation speed of the exogenous interference wave so as to suppress the exogenous interference wave.
Of course, in other embodiments, other techniques for suppressing the external interference wave may be used to process the seismic data to be processed to suppress the noise of the external interference wave. By suppressing the external interference waves to the seismic data to be processed, the quality of the seismic data can be improved.
Further, when the f-kx-ky cone filtering technology is used for processing the to-be-processed seismic data, the f-kx-ky cone filtering technology can be used for processing the to-be-processed seismic data in a cross arrangement domain according to the coordinates of the exogenous interference wave and the propagation speed of the exogenous interference wave so as to suppress the exogenous interference wave.
In one embodiment, the method for processing seismic data for suppressing an external disturbance wave processes the seismic data to be processed according to coordinates of the external disturbance wave and a propagation speed of the external disturbance wave to suppress the external disturbance wave, further includes: and outputting target seismic data after suppressing the exogenous interference.
The method has the advantages that the external source interference wave is suppressed by processing the to-be-processed seismic data, and normal seismic waves are easier to analyze by utilizing the target seismic data after suppressing the external source interference wave, so that the method is used for data processing of seismic exploration.
Example two
By way of a specific example, the results of applying the prior art method and the seismic data processing method for suppressing an external disturbance wave of the present application are compared.
Referring to fig. 2, the coordinates of the source point position P point of the external interference wave are 456, 1238, the propagation speed of the external interference wave is 2000m/s, the receiving distance of the detectors used for seismic observation is 300m, and the track distance between two adjacent detectors is 50 m. The seismic record may refer to fig. 3.
The observation positioning method in the prior art and the seismic data processing method for suppressing the exogenous interference wave are adopted to calculate the propagation speed of the exogenous interference wave and position the source point of the exogenous interference wave, and the realization process and the result are as follows:
(1) Observing and positioning method:
① A data processor browses the seismic records to find seismic data containing exogenous interference waves; ② Determining the vertex position of the hyperbola G in FIG. 3; ③ Reading the abscissa of the seismic traces R5 and R6 at the vertex positions, which are 450 and 500 respectively; ④ Estimating the abscissa of a source point of the exogenous interference wave to be 480; ⑤ Sorting and extracting longitudinal seismic records; ⑥ Repeating ①-④ processes, determining that the source point position of the external interference wave is between L1 and L2, and estimating the ordinate of the source point of the external interference wave to be 1345; ⑦ And measuring the propagation speed of the exogenous interference wave on the seismic record to be 1800m/s.
In the method, the propagation speed of the measured external interference wave is 1800m/s, and the estimated source point coordinates of the external interference wave are (480, 1345). Propagation speed error 11%, positioning error (24, 107) meters respectively.
(2) The application relates to a seismic data processing method (travel time inversion method) for suppressing exogenous interference waves
① A data processor browses the seismic records to find seismic data containing exogenous interference waves; ② Picking up the time of one of the reflection phase axes along a hyperbola G in fig. 3; ③ According to the method of the present application, the coordinates of the point P are calculated as (455, 1235), and the propagation velocity of the external disturbance wave is calculated as 2003m/s.
In the method, the propagation speed of the calculated external disturbance wave is 2003m/s, and the coordinates of the source point of the calculated external disturbance wave are 455, 1235. Propagation speed error 0.15%, positioning error (1, 3) meters respectively.
According to the seismic data processing method for suppressing the exogenous interference wave, disclosed by the application, a data processor is not required to estimate the source point position of the exogenous interference wave through professional technical capability, the source point position of the exogenous interference wave can be automatically obtained, the accuracy of the source point coordinates of the exogenous interference wave is high, the method is suitable for various types of exogenous interference waves, the propagation speed of the exogenous interference wave can be obtained while the source point coordinates of the exogenous interference wave are obtained, the suppressing effect of the exogenous interference wave can be better improved, and the quality of seismic data is improved.
Example III
An embodiment of the present invention discloses a storage medium storing a program for executing the method of any one of the above. Wherein the program comprises: inputting initial seismic data; loading coordinate information on the initial seismic data and preprocessing the initial seismic data to obtain seismic data to be processed; acquiring coordinates (Rx, ry) of at least four points of an in-phase axis of an external interference wave and travel time T corresponding to the coordinates of the at least four points aiming at the seismic data to be processed; propagation equation based on exogenous interference waveConstructing an equation set by using at least four points, wherein v is the propagation speed of the exogenous interference wave, (Px, py) is the coordinates of the exogenous interference wave, and T0 is the vibration excitation time of the exogenous interference wave; solving the equation set to obtain coordinates of the exogenous interference wave and propagation speed of the exogenous interference wave; the step of processing the seismic data to be processed according to the coordinates of the external interference wave and the propagation speed of the external interference wave to suppress the external interference wave is the same as the step of the above seismic data processing method for suppressing the external interference wave, and will not be described herein.
Furthermore, in one embodiment, loading the initial seismic data with coordinate information includes:
and loading the coordinate information of the wave detection points on the initial seismic data.
Further, in one embodiment, preprocessing the initial seismic data includes:
And carrying out static correction, filtering and gain on the initial seismic data.
Furthermore, in one embodiment, the at least four points are located on the same phase axis.
Furthermore, in one embodiment, the at least four points are located on the same phase axis of the first arrival time.
Further, in one embodiment, solving the system of equations to obtain the coordinates of the external disturbance wave and the propagation velocity of the external disturbance wave includes:
And solving the equation set by using a least square method to obtain the coordinates of the external disturbance wave and the propagation speed of the external disturbance wave.
Further, in one embodiment, processing the seismic data to be processed to suppress the exogenous interference wave according to the coordinates of the exogenous interference wave and the propagation velocity of the exogenous interference wave includes:
and processing the seismic data to be processed by utilizing an f-kx-ky cone filtering technology according to the coordinates of the exogenous interference wave and the propagation speed of the exogenous interference wave so as to suppress the exogenous interference wave.
Furthermore, in one embodiment, the seismic data to be processed is processed in a cross-shaped permutation zone to suppress the exogenous interference wave using an f-kx-ky pyramid filtering technique according to the coordinates of the exogenous interference wave and the propagation velocity of the exogenous interference wave.
Furthermore, in one embodiment, the method further comprises:
And outputting target seismic data after suppressing the exogenous interference.
The storage medium stores a program of the seismic data processing method for suppressing the exogenous interference wave, a data processor is not required to estimate the source point position of the exogenous interference wave through professional technology capability, the source point position of the exogenous interference wave can be automatically obtained, the accuracy of the source point coordinates of the exogenous interference wave is high, the storage medium is suitable for various types of exogenous interference waves, the propagation speed of the exogenous interference wave can be obtained while the source point coordinates of the exogenous interference wave are obtained, the suppression effect of the exogenous interference wave can be better improved, and the quality of seismic data is improved.
The present invention is not limited to the above-mentioned embodiments, and any changes or substitutions that can be easily understood by those skilled in the art within the scope of the present invention are intended to be included in the scope of the present invention. Therefore, the protection scope of the present invention should be subject to the protection scope of the claims.
Claims (9)
1. A method of seismic data processing for suppressing an external disturbance wave, comprising:
inputting initial seismic data;
Loading coordinate information on the initial seismic data and preprocessing the initial seismic data to obtain seismic data to be processed;
Aiming at the seismic data to be processed, acquiring each point coordinate of at least four points of an in-phase axis of an external interference wave and travel time corresponding to each point coordinate, wherein the at least four points are positioned on the same in-phase axis;
based on a propagation equation of the exogenous interference wave, constructing an equation set based on each point coordinate of the at least four points and travel time corresponding to each point coordinate;
Solving the equation set to obtain coordinates of the exogenous interference wave and propagation speed of the exogenous interference wave;
and processing the seismic data to be processed according to the coordinates of the exogenous interference wave and the propagation speed of the exogenous interference wave so as to suppress the exogenous interference wave.
2. The method for processing seismic data suppressing an external disturbance wave according to claim 1, wherein loading the initial seismic data with coordinate information includes:
and loading the coordinate information of the wave detection points on the initial seismic data.
3. The method for processing seismic data suppressing an external disturbance wave according to claim 1, wherein preprocessing the initial seismic data includes:
And carrying out static correction, filtering and gain on the initial seismic data.
4. The method for processing seismic data suppressing an external disturbance wave according to claim 1, wherein the at least four points are located on an in-phase axis of a first arrival time.
5. The seismic data processing method for suppressing an external interference wave according to claim 1, wherein solving the system of equations to obtain coordinates of the external interference wave and propagation velocity of the external interference wave comprises:
And solving the equation set by using a least square method to obtain the coordinates of the external disturbance wave and the propagation speed of the external disturbance wave.
6. The seismic data processing method for suppressing an external interfering wave according to claim 1, wherein processing the seismic data to be processed according to coordinates of the external interfering wave and a propagation velocity of the external interfering wave to suppress the external interfering wave comprises:
And processing the seismic data to be processed by utilizing an f-kx-ky cone filtering technology according to the coordinates of the exogenous interference wave and the propagation speed of the exogenous interference wave so as to suppress the exogenous interference wave.
7. The method for processing seismic data for suppressing an external disturbance wave according to claim 6,
And processing the seismic data to be processed in a cross arrangement domain by utilizing an f-kx-ky cone filtering technology according to the coordinates of the exogenous interference wave and the propagation speed of the exogenous interference wave so as to suppress the exogenous interference wave.
8. The seismic data processing method for suppressing an external disturbance wave according to claim 1, wherein the method further comprises:
And outputting target seismic data after suppressing the exogenous interference.
9. A storage medium storing a program for executing the method of any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910817071.3A CN112444869B (en) | 2019-08-30 | 2019-08-30 | Seismic data processing method and storage medium for suppressing exogenous interference waves |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910817071.3A CN112444869B (en) | 2019-08-30 | 2019-08-30 | Seismic data processing method and storage medium for suppressing exogenous interference waves |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112444869A CN112444869A (en) | 2021-03-05 |
| CN112444869B true CN112444869B (en) | 2024-06-18 |
Family
ID=74734939
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910817071.3A Active CN112444869B (en) | 2019-08-30 | 2019-08-30 | Seismic data processing method and storage medium for suppressing exogenous interference waves |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112444869B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101551463A (en) * | 2008-04-03 | 2009-10-07 | 中国石油天然气集团公司 | Noise suppression evaluation method for three-dimensional observation system |
| WO2017108669A1 (en) * | 2015-12-22 | 2017-06-29 | Shell Internationale Research Maatschappij B.V. | Method and system for generating a seismic gather |
Family Cites Families (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100349011C (en) * | 2005-06-03 | 2007-11-14 | 中国石油集团东方地球物理勘探有限责任公司 | Method for suppressing background noise not related to excitation source in seismic data enhancement |
| CN101598809A (en) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | A kind of self-adaptation is eliminated the method for linear programming noise and multiple reflection interference |
| CN102478671A (en) * | 2010-11-23 | 2012-05-30 | 中国石油天然气集团公司 | Method for suppressing controllable seismic-source harmonic-wave interference |
| CN102262242B (en) * | 2011-07-15 | 2013-10-16 | 中国石油天然气集团公司 | Method for accurately positioning position of interference source in seismic exploration data acquisition |
| RU2451307C1 (en) * | 2011-07-18 | 2012-05-20 | Закрытое акционерное общество "Научно-инженерный центр "СИНАПС" | Method of measuring coordinates microseismic sources |
| US9671511B2 (en) * | 2012-08-31 | 2017-06-06 | Cgg Services Sas | Horizontal streamer broadband marine seismic acquisition configuration and processing |
| CN102890291B (en) * | 2012-10-09 | 2015-09-02 | 中国海洋石油总公司 | The self-adaptation drawing method of external source interference in a kind of seismic data |
| CN103852790A (en) * | 2014-03-31 | 2014-06-11 | 中国石油大学(北京) | Morphological filtering method and system of offshore earthquake data |
| CN105093355B (en) * | 2015-06-24 | 2018-07-13 | 中国石油天然气集团公司 | A kind of seismic observation system evaluation method |
| US10578755B2 (en) * | 2015-12-11 | 2020-03-03 | Ion Geophysical Corporation | System and method for reconstructed wavefield inversion |
| CN105607121B (en) * | 2016-02-02 | 2016-12-21 | 中国矿业大学(北京) | A kind of coal karst collapse col umn recognition methods and device |
| CN105911585B (en) * | 2016-07-05 | 2018-05-15 | 中国石油集团东方地球物理勘探有限责任公司 | A kind of extracting method and device of earthquake record regular interference |
| CN109143368B (en) * | 2017-06-28 | 2020-04-07 | 中国石油化工股份有限公司 | Pre-mountain area earthquake denoising method |
-
2019
- 2019-08-30 CN CN201910817071.3A patent/CN112444869B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101551463A (en) * | 2008-04-03 | 2009-10-07 | 中国石油天然气集团公司 | Noise suppression evaluation method for three-dimensional observation system |
| WO2017108669A1 (en) * | 2015-12-22 | 2017-06-29 | Shell Internationale Research Maatschappij B.V. | Method and system for generating a seismic gather |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112444869A (en) | 2021-03-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110133715B (en) | Microseism seismic source positioning method based on first-arrival time difference and waveform superposition | |
| Verdon et al. | Microseismic monitoring using a fiber-optic distributed acoustic sensor array | |
| CN112213768B (en) | Ground micro-seismic positioning method and system combining seismic source mechanism inversion | |
| AU2007272702B2 (en) | Time lapse marine seismic surveying employing interpolated multicomponent streamer pressure data | |
| CN106443788A (en) | Using a rotation sensor measurement to attenuate noise acquired by a streamer-disposed sensor | |
| CN101881836A (en) | Be used for calculating the method for seismic properties according to seismic signal | |
| RU2339056C2 (en) | Generic three-dimension prediction of multiple events from surface | |
| CN110389377B (en) | Microseism offset imaging positioning method based on waveform cross-correlation coefficient multiplication | |
| WO2024119967A1 (en) | Multi-excitation-point active surface wave exploration method and device, and storage medium | |
| CN108089229B (en) | Detection device position determining method and device and computer storage medium | |
| Zhang et al. | Automated microseismic event location by amplitude stacking and semblance | |
| CN108693560B (en) | Scattered wave imaging method and system based on cross-correlation channel | |
| CN111352153B (en) | Microseism interference positioning method based on instantaneous phase cross-correlation weighting | |
| GB2450163A (en) | Detecting the location of seismic events without picking events in received seismic wave data | |
| CN112444869B (en) | Seismic data processing method and storage medium for suppressing exogenous interference waves | |
| CN110967751B (en) | Positioning method of micro-seismic event based on ground shallow well monitoring and storage medium | |
| CN110579799B (en) | Seismic acquisition observation method and system with equal travel time intervals | |
| CN107664771B (en) | A kind of microseism Full wave shape localization method based on likeness coefficient | |
| CN110780346A (en) | Separation method for advanced detection of complex seismic wave field in tunnel | |
| Yue et al. | The theory and application of the random dislocation directional seismic wave technique in a tunnel environment | |
| CN111965705B (en) | Method, device, equipment and medium for calibrating geological horizon recorded by earthquake single shot | |
| CN114114407A (en) | Surface wave and direct transverse wave pressing processing method for seismic wave detection | |
| CN101118287A (en) | Spectrum equalizing and non time difference combined method in seismic prospecting single cymoscope recording channel chamber | |
| Mao et al. | Microseismic event location using an improved global grid search and its extended method in a downhole monitoring system | |
| CN113075725B (en) | Method for determining minimum offset distance and system delay of multi-channel earthquake in deep water area |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |