CN114185094A - RMS-SVD multiple suppression method, device, electronic apparatus and medium - Google Patents
RMS-SVD multiple suppression method, device, electronic apparatus and medium Download PDFInfo
- Publication number
- CN114185094A CN114185094A CN202010961833.XA CN202010961833A CN114185094A CN 114185094 A CN114185094 A CN 114185094A CN 202010961833 A CN202010961833 A CN 202010961833A CN 114185094 A CN114185094 A CN 114185094A
- Authority
- CN
- China
- Prior art keywords
- target data
- svd
- data gather
- filtering
- multiples
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 230000001629 suppression Effects 0.000 title claims abstract description 29
- 238000001914 filtration Methods 0.000 claims abstract description 87
- 238000012937 correction Methods 0.000 claims abstract description 53
- 238000012545 processing Methods 0.000 claims abstract description 32
- 238000001228 spectrum Methods 0.000 claims description 28
- 230000008569 process Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000006835 compression Effects 0.000 claims description 6
- 238000007906 compression Methods 0.000 claims description 6
- 230000008030 elimination Effects 0.000 description 5
- 238000003379 elimination reaction Methods 0.000 description 5
- 238000003825 pressing Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 238000011410 subtraction method Methods 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- RLBIQVVOMOPOHC-UHFFFAOYSA-N parathion-methyl Chemical compound COP(=S)(OC)OC1=CC=C([N+]([O-])=O)C=C1 RLBIQVVOMOPOHC-UHFFFAOYSA-N 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/52—Move-out correction
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
Disclosed are an RMS-SVD multiple suppression method, an RMS-SVD multiple suppression device, an electronic apparatus and a medium. The method can comprise the following steps: picking up the superposition speed of the multiples to further obtain the root-mean-square speed of the multiples; performing dynamic correction processing on the target data gather through the multiple root mean square speed to obtain a dynamically corrected target data gather; determining the time window range of the multiple waves, and further determining the time window length of SVD filtering; SVD filtering is carried out on the target data gather after the dynamic correction, and a filtered target data gather is obtained; and performing reverse motion correction processing on the filtered target data gather to obtain a multiple-suppressed target data gather. The invention can effectively suppress and eliminate the multiple energy on the target data gather and improve the signal-to-noise ratio of the target data gather through the difference between the multiple and the primary reflected wave speed.
Description
Technical Field
The invention relates to the field of seismic data processing, in particular to an RMS-SVD multiple suppression method, an RMS-SVD multiple suppression device, electronic equipment and a medium.
Background
In seismic exploration, people usually know the condition of underground geological structures by using only primary reflection wave information in seismic data, but in seismic data acquired in the field, besides primary reflection waves, many other wave field information exist, so that people need to process the seismic data from the field in order to extract the primary reflection wave information of the seismic data in the field, and the purposes of suppressing and eliminating other interference wave field information and retaining and enhancing the primary reflection wave field information are achieved through series processing. Among the other wavefield information, many types of multiple reflections are often seen in seismic data, which is a daily multiple known as a seismic wave that undergoes multiple upward reflections during propagation through the subsurface medium. Depending on the position of the shallowest downward reflecting interface, multiples can be divided into free surface multiples and interbed multiples. Wherein the shallowest downward reflection of the free surface multiples occurs at the earth surface and the downward reflections of the interbed multiples all occur at the subsurface interface. To date, these common multiple reflections are not used as effective reflection information in most cases, but only as a source of interference which is difficult to handle. Meanwhile, as the multiple waves undergo multiple reflections in the underground, the multiple waves tend to have certain periodicity on a stacking section, and the inclination angle of the same phase axis of the multiple waves is doubled and the amplitude of the multiple waves is amplified along with the increase of the order, and the phenomenon of focusing or defocusing of local structures is accompanied. In general, free surface multiples are primarily present in marine seismic data and occupy a significant portion of the multiple wave energy, while interbed multiples are generally associated with strong wave impedance interfaces in the subsurface medium. Generally speaking, different types and orders of multiples in seismic data interfere with the primary wave, which seriously affects the seismic imaging quality, reduces the reliability of inversion results, and increases the complexity of seismic interpretation. Thus, multiples are typically suppressed as noise during the pre-stack processing stage of seismic data. How to eliminate or attenuate multiple waves generated by various types of multiple reflections in the process of seismic data processing has been a processing technique of great interest. It should be noted that there is no complete processing technique to adapt to the attenuation of various types of multiple waves, and only a processing method is selected through experimental analysis of actual data. At present, the multiple suppression methods generally applied in actual seismic processing are mainly classified into two types: one class suppresses multiples by using the difference between the characteristics and properties of multiples and primaries, and is called a filtering method; the other type firstly predicts a multi-wave model from the seismic data and then adaptively subtracts the multi-wave model from the seismic data, and the model is called wave equation prediction subtraction; in addition to the two multiples suppression methods, the DELPHI project group of the DELPHI theory of dalf theory of the netherlands developed in recent years a full waveform inversion method that can directly estimate the primaries and multiples, called sparse inversion method. It should be said that, the above three processing methods for suppressing multiples all have respective advantages and disadvantages, the filtering method mainly utilizes the periodicity of the multiples and the separability of the multiples and the primaries in different transform domains to attenuate the multiples, the method is easy to implement and has small calculation amount, but under the condition of complex exploration geological conditions, because the periodicity of the multiples is not obvious and the separability is low, the effect of suppressing the multiples is not ideal and effective signals can be damaged; the prediction subtraction method is established on the basis of a fluctuation theory describing the seismic wave propagation process, adopts a prior velocity model or original seismic data as a drive to predict multiples, and then performs adaptive subtraction from the seismic data to achieve the purpose of attenuating the multiples; compared with a filtering method, the wave equation prediction subtraction method has less assumptions on underground information, can be suitable for multiple suppression under complex geological conditions, but still has the main problems of near offset data loss, adaptive subtraction damage to effective signals and the like in practical application; the sparse inversion method is established on the same theoretical basis as a feedback iteration method, primary waves and multiple waves are estimated in a parameter inversion mode, the step of self-adaptive subtraction is avoided, effective signals can be protected to the maximum degree, in addition, the method can utilize the multiple waves to estimate missing near offset data, the limitation that the conventional prediction subtraction method cannot effectively suppress the shallow water free surface multiple waves is overcome, but the sparse inversion method also has the problems of large calculated amount, sensitivity to seismic amplitude and the like, and further development and improvement are needed.
Because the existing processing methods for suppressing and eliminating multiples have certain limitations, in actual processing work, the purpose of effectively eliminating the multiples is difficult to achieve by using some existing general multiple elimination methods, and people are difficult to achieve the effect of effectively eliminating the interbed multiples by using the conventional multiple suppression method.
Therefore, it is necessary to develop an RMS-SVD multiple suppression method, apparatus, electronic device and medium.
The information disclosed in this background section is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person skilled in the art.
Disclosure of Invention
The invention provides an RMS-SVD multiple suppression method, an RMS-SVD multiple suppression device, an electronic device and a medium, which can realize … … through … ….
In a first aspect, an embodiment of the present disclosure provides a method for RMS-SVD multiple suppression, including:
picking up the superposition speed of the multiples to further obtain the root-mean-square speed of the multiples;
performing dynamic correction processing on the target data gather through the multiple root mean square velocity to obtain a dynamically corrected target data gather;
determining the time window range of the multiple waves, and further determining the time window length of SVD filtering;
SVD filtering is carried out on the target data gather after the dynamic correction, and a filtered target data gather is obtained;
and performing reverse motion correction processing on the filtered target data gather to obtain a multiple-suppressed target data gather.
Preferably, the target data gather is a CRP data gather or a CMP data gather.
Preferably, the range of multiple time windows includes all flattened multiples.
Preferably, the SVD filter has a window length that is the difference between the maximum time and the minimum time of the range of the multiple window.
Preferably, performing SVD filtering on the motion corrected target data gather, and obtaining a filtered target data gather includes:
calculating a singular value spectrum of the seismic data in the time window, and determining the range of the multiple;
determining a filtering factor of SVD filtering;
and carrying out SVD high-pass filtering on the seismic data in the time window to obtain a filtered target data gather.
Preferably, the method further comprises the following steps:
the filter factor is adjusted until the multiple energy is removed.
Preferably, the filtered target data gather is subjected to inverse dynamic correction processing through the multiple root mean square velocity.
As a specific implementation of the embodiments of the present disclosure,
in a second aspect, embodiments of the present disclosure further provide an RMS-SVD multi-wave suppression device, including:
the speed acquisition module is used for picking up the superposition speed of the multiples so as to obtain the root-mean-square speed of the multiples;
the dynamic correction module is used for performing dynamic correction processing on the target data gather through the multiple root mean square speed to obtain a dynamically corrected target data gather;
the calculation module is used for determining the time window range of the multiple waves and further determining the time window length of SVD filtering;
the filtering module is used for carrying out SVD filtering on the target data gather after the dynamic correction to obtain a filtered target data gather;
and the reverse motion correction module is used for performing reverse motion correction processing on the filtered target data gather to obtain a multiple-suppressed target data gather.
Preferably, the target data gather is a CRP data gather or a CMP data gather.
Preferably, the range of multiple time windows includes all flattened multiples.
Preferably, the SVD filter has a window length that is the difference between the maximum time and the minimum time of the range of the multiple window.
Preferably, performing SVD filtering on the motion corrected target data gather, and obtaining a filtered target data gather includes:
calculating a singular value spectrum of the seismic data in the time window, and determining the range of the multiple;
determining a filtering factor of SVD filtering;
and carrying out SVD high-pass filtering on the seismic data in the time window to obtain a filtered target data gather.
Preferably, the method further comprises the following steps:
the filter factor is adjusted until the multiple energy is removed.
Preferably, the filtered target data gather is subjected to inverse dynamic correction processing through the multiple root mean square velocity.
In a third aspect, an embodiment of the present disclosure further provides an electronic device, where the electronic device includes:
a memory storing executable instructions;
a processor executing the executable instructions in the memory to implement the RMS-SVD multiple pressing method.
In a fourth aspect, the disclosed embodiments also provide a computer-readable storage medium, which stores a computer program, and when the computer program is executed by a processor, the computer program implements the RMS-SVD multiple suppression method.
The beneficial effects are that: .
The method and apparatus of the present invention have other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the invention.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent by describing in more detail exemplary embodiments thereof with reference to the attached drawings, in which like reference numerals generally represent like parts.
FIG. 1 shows a flow chart of the steps of a RMS-SVD multiple attenuation method according to the present invention.
FIGS. 2a, 2b, and 2c respectively illustrate a schematic of an original target data gather and its corresponding velocity spectrum, kinematically corrected target data gather, according to one embodiment of the invention.
FIGS. 3a, 3b show schematic diagrams of a multiple-suppressed target data gather and its corresponding velocity spectrum, respectively, according to one embodiment of the invention.
FIG. 4 illustrates a block diagram of an RMS-SVD multi-wave suppression device in accordance with an embodiment of the present invention.
Description of reference numerals:
201. a speed acquisition module; 202. a dynamic correction module; 203. a calculation module; 204. a filtering module; 205. and a reverse correction module.
Detailed Description
Preferred embodiments of the present invention will be described in more detail below. While the following describes preferred embodiments of the present invention, it should be understood that the present invention may be embodied in various forms and should not be limited by the embodiments set forth herein.
The RMS-SVD multiple pressing method is characterized in that the difference of multiples and effective waves on a velocity spectrum based on a CRP (or CMP) data gather is utilized, the root mean square velocity of the multiples is obtained by picking up the velocity spectrum of the multiples, the CRP (or CMP) data gather is subjected to dynamic correction by utilizing the root mean square velocity of the multiples, the multiples on the data gather are flattened, then the linear interference waves can be effectively eliminated by utilizing an SVD filtering method (also called a singular value filtering method), the SVD filtering is carried out on the flattened data gather to eliminate the multiples, finally, the root mean square velocity of the multiples is utilized to carry out inverse dynamic correction on the data gather from which the horizontal multiple energy is removed, and thus the CRP (or CMP) data gather from which the multiples are removed is obtained.
The invention provides an RMS-SVD multiple pressing method, which comprises the following steps:
and picking up the superposition speed of the multiples to further obtain the root-mean-square speed of the multiples.
Specifically, the velocity of superposition of multiples is picked up on the velocity spectrum of the target data gather, and then the picked-up multiple superposition velocity is converted into a multiple root-mean-square velocity by velocity conversion.
Performing dynamic correction processing on the target data gather through the multiple root mean square speed to obtain a dynamically corrected target data gather; in one example, the target data gather is a CRP data gather or a CMP data gather.
Specifically, the root-mean-square velocity of the multiples is used as the dynamic correction velocity of the data gather to dynamically correct the CRP (or CMP) data gather containing the multiples, so as to obtain a target data gather after dynamic correction, and thus, the purpose of leveling the multiples in the data gather is achieved.
Determining the time window range of the multiple waves, and further determining the time window length of SVD filtering; in one example, the range of multiple time windows contains all of the flattened multiples. In one example, the SVD filtering has a window length that is the difference between the maximum time and the minimum time of the range of the multiple time window.
Specifically, in the target data track set with the multiple flattened, a time window range of the multiple is determined, all the flattened multiple should be included in the multiple time window range, and the difference between the maximum time and the minimum time of the multiple time window range is used as the time window length of the SVD filtering.
SVD filtering is carried out on the target data gather after the dynamic correction, and a filtered target data gather is obtained; in one example, SVD filtering the motion corrected target data gather, obtaining a filtered target data gather includes: calculating a singular value spectrum of the seismic data in the time window, and determining the range of the multiple; determining a filtering factor of SVD filtering; and carrying out SVD high-pass filtering on the seismic data in the time window to obtain a filtered target data gather. In one example, further comprising: the filter factor is adjusted until the multiple energy is removed.
Specifically, SVD filtering is carried out on a target data gather with flattened multiples, firstly, a singular value spectrum of seismic data in a time window is calculated, the range of the multiples is determined on the singular value spectrum, then a filtering factor of the SVD filtering is determined, because the flattened multiples belong to strong energy, the multiples are represented as large singular values on the singular value spectrum, the elimination of the multiples is actually carrying out SVD high-pass filtering on the seismic data in the time window, if the multiples are not eliminated enough, the value of the filtering factor is increased, if the filtering is too excessive, the loss of effective wave information is reduced, the value of the filtering factor is reduced, and filtering parameters are continuously adjusted until the horizontal multiple wave energy can be effectively eliminated.
And performing reverse motion correction processing on the filtered target data gather to obtain a multiple-suppressed target data gather. In one example, a reverse motion correction process is performed on the filtered target data gather by a multiple root mean square velocity.
Specifically, the filtered target data gather is reverse corrected using the root mean square velocity of the multiples to obtain a multiple-suppressed target data gather.
The invention also provides an RMS-SVD multi-wave suppression device, comprising:
and the speed acquisition module is used for picking up the superposition speed of the multiples and further acquiring the root-mean-square speed of the multiples.
Specifically, the velocity of superposition of multiples is picked up on the velocity spectrum of the target data gather, and then the picked-up multiple superposition velocity is converted into a multiple root-mean-square velocity by velocity conversion.
The dynamic correction module is used for performing dynamic correction processing on the target data gather through the multiple root mean square speed to obtain a dynamically corrected target data gather; in one example, the target data gather is a CRP data gather or a CMP data gather.
Specifically, the root-mean-square velocity of the multiples is used as the dynamic correction velocity of the data gather to dynamically correct the CRP (or CMP) data gather containing the multiples, so as to obtain a target data gather after dynamic correction, and thus, the purpose of leveling the multiples in the data gather is achieved.
The calculation module is used for determining the time window range of the multiple waves and further determining the time window length of SVD filtering; in one example, the range of multiple time windows contains all of the flattened multiples. In one example, the SVD filtering has a window length that is the difference between the maximum time and the minimum time of the range of the multiple time window.
Specifically, in the target data track set with the multiple flattened, a time window range of the multiple is determined, all the flattened multiple should be included in the multiple time window range, and the difference between the maximum time and the minimum time of the multiple time window range is used as the time window length of the SVD filtering.
The filtering module is used for carrying out SVD filtering on the target data gather after the dynamic correction to obtain a filtered target data gather; in one example, SVD filtering the motion corrected target data gather, obtaining a filtered target data gather includes: calculating a singular value spectrum of the seismic data in the time window, and determining the range of the multiple; determining a filtering factor of SVD filtering; and carrying out SVD high-pass filtering on the seismic data in the time window to obtain a filtered target data gather. In one example, further comprising: the filter factor is adjusted until the multiple energy is removed.
Specifically, SVD filtering is carried out on a target data gather with flattened multiples, firstly, a singular value spectrum of seismic data in a time window is calculated, a range of the multiples is determined on the singular value spectrum, then a filtering factor of the SVD filtering is determined, the filtering factor is determined, firstly, a range of a singular value corresponding to a filtering object is determined on the singular value spectrum, then an upper bound value corresponding to the range is used as the filtering factor, if the energy of the corresponding range is filtered, high-pass filtering is carried out by taking the upper bound value as the filtering factor, if the energy of the corresponding range is strengthened, low-pass filtering is carried out by taking the upper bound value as the filtering factor, because the flattened multiples belong to strong energy, the multiples are represented as larger singular values on the singular value spectrum, the elimination of the multiples is actually carried out by SVD high-pass filtering on the seismic data in the time window, and if the multiples are not eliminated enough, the value of the filter factor is increased and decreased if the filtering is too high resulting in a loss of useful wave information, by continually adjusting the filter parameters until the horizontal multiple wave energy can be effectively eliminated.
And the reverse motion correction module is used for performing reverse motion correction processing on the filtered target data gather to obtain a multiple-suppressed target data gather. In one example, a reverse motion correction process is performed on the filtered target data gather by a multiple root mean square velocity.
Specifically, the filtered target data gather is reverse corrected using the root mean square velocity of the multiples to obtain a multiple-suppressed target data gather.
The present invention also provides an electronic device, comprising: a memory storing executable instructions; and a processor executing executable instructions in the memory to implement the RMS-SVD multiple compression method described above.
The present invention also provides a computer-readable storage medium storing a computer program which, when executed by a processor, implements the RMS-SVD multiple suppression method described above.
To facilitate understanding of the scheme of the embodiments of the present invention and the effects thereof, four specific application examples are given below. It will be understood by those skilled in the art that this example is merely for the purpose of facilitating an understanding of the present invention and that any specific details thereof are not intended to limit the invention in any way.
Example 1
FIG. 1 shows a flow chart of the steps of a RMS-SVD multiple attenuation method according to the present invention.
As shown in FIG. 1, the RMS-SVD multiple suppression method includes: step 101, picking up the superposition speed of multiples to further obtain the root-mean-square speed of the multiples; 102, performing dynamic correction processing on a target data gather through a multiple root mean square speed to obtain a dynamically corrected target data gather; step 103, determining the time window range of the multiple waves, and further determining the time window length of SVD filtering; 104, performing SVD filtering on the target data gather after the dynamic correction to obtain a filtered target data gather; and 105, performing reverse motion correction processing on the filtered target data gather to obtain a multiple-suppressed target data gather.
The velocity of superposition of multiples is picked up on the velocity spectrum of the target data gather, and then the picked-up multiple superposition velocity is converted into a multiple root-mean-square velocity by velocity conversion.
And (3) dynamically correcting the CRP (or CMP) data gather containing the multiples by taking the root-mean-square velocity of the multiples as the dynamic correction velocity of the data gather to obtain a target data gather after dynamic correction so as to achieve the aim of leveling the multiples in the data gather.
And determining the time window range of the multiple in the target data channel set with the multiple flattened, wherein all the flattened multiple should be contained in the multiple time window range, and the difference between the maximum time and the minimum time of the multiple time window range is used as the time window length of the SVD filtering.
SVD filtering is carried out on a target data gather with flattened multiples, firstly, a singular value spectrum of seismic data in a time window is calculated, the range of the multiples is determined on the singular value spectrum, then a filtering factor of the SVD filtering is determined, the flattened multiples belong to strong energy, therefore, the multiples are represented as larger singular values on the singular value spectrum, the elimination of the multiples is really carried out on the seismic data in the time window by SVD high-pass filtering, if the elimination of the multiples is insufficient, the value of the filtering factor is increased, if the filtering is too excessive, the loss of effective wave information is caused, the value of the filtering factor is reduced, and filtering parameters are continuously adjusted until the horizontal multiple wave energy can be effectively eliminated.
And performing inverse correction on the filtered target data gather by using the root-mean-square velocity of the multiples to obtain a multiple-suppressed target data gather.
On a CRP data gather in three-dimensional seismic data of a certain working area, the interference of multiple waves is serious, and in order to verify the effectiveness of the method, according to the technical scheme provided by the invention, the multiple wave pressing test treatment is carried out on Common Reflection Point (CRP) data gather data of a certain longitudinal measuring line of the working area by using the multiple wave removing method provided by the invention, and the result is as follows:
fig. 2a, 2b, and 2c show schematic views of a raw CRP data gather and its corresponding velocity spectrum, respectively, a kinematically corrected CRP data gather, in accordance with one embodiment of the present invention.
This is also confirmed by the strong multiples interference energy on the data gather, which is readily apparent from FIG. 2a, and the strong multiples energy bolus on the CRP velocity spectrum, which is shown in FIG. 2 b. By picking up the multiple velocity spectrum of the CRP, the multiple energy axes in the CRP gather are flattened out, as shown in FIG. 2 c.
FIGS. 3a, 3b show schematic diagrams of a multiple-suppressed target data gather and its corresponding velocity spectrum, respectively, according to one embodiment of the invention.
After the multiple removal processing is performed on the selected CRP data gather data, it can be readily seen that the multiple on the data gather of fig. 3a is removed and the multiple wave energy clusters on the velocity spectrum of fig. 3b are also missing, leaving only the primary reflected wave energy.
Comparing fig. 2a-c with fig. 3a-b, the following conclusions can be drawn: by utilizing the RMS-SVD multiple suppression method to process the CRP (or CMP) seismic data gather with serious multiple interference, the multiple interference influence in the data gather can be effectively eliminated, and the signal-to-noise ratio of the CRP (or CMP) seismic data gather can be improved to a certain extent.
Example 2
FIG. 4 illustrates a block diagram of an RMS-SVD multi-wave suppression device in accordance with an embodiment of the present invention.
As shown in fig. 4, the RMS-SVD multi-wave suppression device includes:
the speed acquisition module 201 is used for picking up the superposition speed of the multiples and further acquiring the root-mean-square speed of the multiples;
the dynamic correction module 202 is used for performing dynamic correction processing on the target data gather through the multiple root mean square speed to obtain a dynamically corrected target data gather;
the calculation module 203 determines the time window range of the multiples, and further determines the time window length of the SVD filtering;
the filtering module 204 is used for performing SVD filtering on the target data gather after the dynamic correction to obtain a filtered target data gather;
and the reverse motion correction module 205 performs reverse motion correction processing on the filtered target data gather to obtain a multiple-suppressed target data gather.
Alternatively, the target data gather is a CRP data gather or a CMP data gather.
Alternatively, the range of multiple time windows includes all of the flattened multiples.
Alternatively, the SVD filter has a window length that is the difference between the maximum time and the minimum time of the multiple window range.
As an alternative, performing SVD filtering on the target data gather after the dynamic correction to obtain a filtered target data gather includes:
calculating a singular value spectrum of the seismic data in the time window, and determining the range of the multiple;
determining a filtering factor of SVD filtering;
and carrying out SVD high-pass filtering on the seismic data in the time window to obtain a filtered target data gather.
As an alternative, the method further comprises the following steps:
the filter factor is adjusted until the multiple energy is removed.
Alternatively, the inverse motion correction process is performed on the filtered target data gather by the multiple root mean square velocity.
Example 3
The present disclosure provides an electronic device including: a memory storing executable instructions; and the processor runs the executable instructions in the memory to realize the RMS-SVD multiple pressing method.
An electronic device according to an embodiment of the present disclosure includes a memory and a processor.
The memory is to store non-transitory computer readable instructions. In particular, the memory may include one or more computer program products that may include various forms of computer-readable storage media, such as volatile memory and/or non-volatile memory. The volatile memory may include, for example, Random Access Memory (RAM), cache memory (cache), and/or the like. The non-volatile memory may include, for example, Read Only Memory (ROM), hard disk, flash memory, etc.
The processor may be a Central Processing Unit (CPU) or other form of processing unit having data processing capabilities and/or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of the disclosure, the processor is configured to execute the computer readable instructions stored in the memory.
Those skilled in the art should understand that, in order to solve the technical problem of how to obtain a good user experience, the present embodiment may also include well-known structures such as a communication bus, an interface, and the like, and these well-known structures should also be included in the protection scope of the present disclosure.
For the detailed description of the present embodiment, reference may be made to the corresponding descriptions in the foregoing embodiments, which are not repeated herein.
Example 4
The disclosed embodiments provide a computer-readable storage medium storing a computer program that, when executed by a processor, implements the RMS-SVD multiple suppression method.
A computer-readable storage medium according to an embodiment of the present disclosure has non-transitory computer-readable instructions stored thereon. The non-transitory computer readable instructions, when executed by a processor, perform all or a portion of the steps of the methods of the embodiments of the disclosure previously described.
The computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROMs and DVDs), magneto-optical storage media (e.g., MOs), magnetic storage media (e.g., magnetic tapes or removable disks), media with built-in rewritable non-volatile memory (e.g., memory cards), and media with built-in ROMs (e.g., ROM cartridges).
It will be appreciated by persons skilled in the art that the above description of embodiments of the invention is intended only to illustrate the benefits of embodiments of the invention and is not intended to limit embodiments of the invention to any examples given.
Having described embodiments of the present invention, the foregoing description is intended to be exemplary, not exhaustive, and not limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.
Claims (10)
1. An RMS-SVD multiple compression method, comprising:
picking up the superposition speed of the multiples to further obtain the root-mean-square speed of the multiples;
performing dynamic correction processing on the target data gather through the multiple root mean square velocity to obtain a dynamically corrected target data gather;
determining the time window range of the multiple waves, and further determining the time window length of SVD filtering;
SVD filtering is carried out on the target data gather after the dynamic correction, and a filtered target data gather is obtained;
and performing reverse motion correction processing on the filtered target data gather to obtain a multiple-suppressed target data gather.
2. The RMS-SVD multiple compression method as claimed in claim 1, wherein said target data gather is a CRP data gather or a CMP data gather.
3. The RMS-SVD multiple suppression method according to claim 1, wherein the multiple time window range includes all flattened multiples.
4. The RMS-SVD multiple compression method as claimed in claim 1, wherein the SVD filtering has a window length that is the difference between the maximum time and the minimum time of the multiple window range.
5. The RMS-SVD multiple suppression method of claim 1, wherein SVD filtering the motion corrected target data gather to obtain a filtered target data gather comprises:
calculating a singular value spectrum of the seismic data in the time window, and determining the range of the multiple;
determining a filtering factor of SVD filtering;
and carrying out SVD high-pass filtering on the seismic data in the time window to obtain a filtered target data gather.
6. The RMS-SVD multiple compression method as claimed in claim 5, further comprising:
the filter factor is adjusted until the multiple energy is removed.
7. The RMS-SVD multiple compression method as claimed in claim 1, wherein said inverse motion correction process is performed on said filtered target data gather by said multiple root mean square velocity.
8. An RMS-SVD multi-wave suppression device, comprising:
the speed acquisition module is used for picking up the superposition speed of the multiples so as to obtain the root-mean-square speed of the multiples;
the dynamic correction module is used for performing dynamic correction processing on the target data gather through the multiple root mean square speed to obtain a dynamically corrected target data gather;
the calculation module is used for determining the time window range of the multiple waves and further determining the time window length of SVD filtering;
the filtering module is used for carrying out SVD filtering on the target data gather after the dynamic correction to obtain a filtered target data gather;
and the reverse motion correction module is used for performing reverse motion correction processing on the filtered target data gather to obtain a multiple-suppressed target data gather.
9. An electronic device, characterized in that the electronic device comprises:
a memory storing executable instructions;
a processor executing the executable instructions in the memory to implement the RMS-SVD multiple suppression method of any one of claims 1-7.
10. A computer-readable storage medium, wherein a computer program is stored, which when executed by a processor implements the RMS-SVD multiple suppression method of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010961833.XA CN114185094A (en) | 2020-09-14 | 2020-09-14 | RMS-SVD multiple suppression method, device, electronic apparatus and medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010961833.XA CN114185094A (en) | 2020-09-14 | 2020-09-14 | RMS-SVD multiple suppression method, device, electronic apparatus and medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114185094A true CN114185094A (en) | 2022-03-15 |
Family
ID=80538995
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010961833.XA Pending CN114185094A (en) | 2020-09-14 | 2020-09-14 | RMS-SVD multiple suppression method, device, electronic apparatus and medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114185094A (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104977615A (en) * | 2015-07-01 | 2015-10-14 | 中国海洋石油总公司 | Model-statistics-pickup-based multiple suppression method of deep-sea OBC data |
| CN109001813A (en) * | 2018-07-16 | 2018-12-14 | 中国石油天然气股份有限公司 | Method, device and system for suppressing multiple waves |
| CN109116423A (en) * | 2018-07-11 | 2019-01-01 | 北京奥能恒业能源技术有限公司 | A kind of diffraction multiple wave drawing method and device |
| CN110568500A (en) * | 2019-08-21 | 2019-12-13 | 中国石油化工股份有限公司 | Method and device for removing multiples in seismic data |
-
2020
- 2020-09-14 CN CN202010961833.XA patent/CN114185094A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104977615A (en) * | 2015-07-01 | 2015-10-14 | 中国海洋石油总公司 | Model-statistics-pickup-based multiple suppression method of deep-sea OBC data |
| CN109116423A (en) * | 2018-07-11 | 2019-01-01 | 北京奥能恒业能源技术有限公司 | A kind of diffraction multiple wave drawing method and device |
| CN109001813A (en) * | 2018-07-16 | 2018-12-14 | 中国石油天然气股份有限公司 | Method, device and system for suppressing multiple waves |
| CN110568500A (en) * | 2019-08-21 | 2019-12-13 | 中国石油化工股份有限公司 | Method and device for removing multiples in seismic data |
Non-Patent Citations (4)
| Title |
|---|
| 刘保童: "波场去噪与分离方法研究", 31 October 2006, 陕西科学技术出版社, pages: 1 * |
| 张树林: "钻前压力预测的关键问题及解决方法", 海洋地质动态, no. 05 * |
| 朱四新: "琼东南盆地速度研究与PSDM应用", 31 August 2015, 地质出版社, pages: 1 * |
| 罗小明, 牛滨华, 于延庆, 黄新武: "海上某区二维地震多次波的识别和压制", 现代地质, no. 04 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2837953A1 (en) | Geologic quality factor inversion method | |
| CN112817047B (en) | Ocean earthquake self-adaptive ghost wave removing method and device, electronic equipment and medium | |
| CN112835103B (en) | Adaptive ghost wave removing and broadband quasi-zero phase deconvolution combined processing method and system | |
| CN105116442A (en) | Lithologic oil-gas reservoir weak-reflection seismic signal reconstruction method | |
| CN110967735A (en) | Self-adaptive ghost wave suppression method and system | |
| CN112817040B (en) | Broadband quasi-zero phase deconvolution processing method, device, electronic equipment and medium | |
| CN110618449B (en) | Method and system for processing seismic data | |
| Wang et al. | Robust singular value decomposition filtering for low signal-to-noise ratio seismic data | |
| CN102163322A (en) | Suppression method of co-channel interference in radar image based on Laplace operator | |
| CN112183407B (en) | Tunnel seismic wave data denoising method and system based on time-frequency domain spectral subtraction | |
| CN107255833B (en) | The determination method and apparatus of Seismic Stacked Section | |
| CN109212609B (en) | Near-surface noise suppression method based on wave equation continuation | |
| CN114185094A (en) | RMS-SVD multiple suppression method, device, electronic apparatus and medium | |
| CN107589453B (en) | A kind of comentropy filter and seismic data random noise attenuation method | |
| RU2126984C1 (en) | Method determining depth and speed parameters of medium and construction of its image by seismic data- prime system | |
| CN110837119B (en) | Seismic data processing method and system for enhancing inverse Q value compensation stability | |
| CN117890979B (en) | Method, system, equipment and medium for compensating self-adaptive weak reflection signals of seismic data | |
| CN111781646B (en) | Amplitude-preserving and fidelity-preserving linear coherent noise suppression method and processing terminal | |
| CN114740530B (en) | Medium-high frequency quasi-linear noise suppression method and device based on hyperbolic time window constraint | |
| CN117741750B (en) | Multi-channel pre-stack deconvolution method and system based on Radon transformation | |
| CN109459788A (en) | Ground interval quality factors calculation method and system | |
| CN117930350A (en) | Land short period interlayer multiple wave pressing method, device, electronic equipment and medium | |
| CN118210026A (en) | High-resolution multiple suppression method for seismic data, electronic equipment and medium | |
| CN116009092A (en) | Method, device, medium and equipment for suppressing remote source interference of seismic exploration | |
| CN114755723A (en) | Method, device, equipment and medium for determining polarity zero phasing of marine seismic data |
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 |