CN111142155B - Complex field weak signal recovery method based on vector surface element and electronic equipment - Google Patents
Complex field weak signal recovery method based on vector surface element and electronic equipment Download PDFInfo
- Publication number
- CN111142155B CN111142155B CN201911406030.1A CN201911406030A CN111142155B CN 111142155 B CN111142155 B CN 111142155B CN 201911406030 A CN201911406030 A CN 201911406030A CN 111142155 B CN111142155 B CN 111142155B
- Authority
- CN
- China
- Prior art keywords
- gather
- vector
- target track
- track
- complex field
- 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
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000011084 recovery Methods 0.000 title claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 29
- 230000006870 function Effects 0.000 claims description 72
- 238000012216 screening Methods 0.000 claims description 6
- 238000004321 preservation Methods 0.000 abstract description 12
- 238000003672 processing method Methods 0.000 abstract description 5
- 108091006146 Channels Proteins 0.000 description 25
- 230000000694 effects Effects 0.000 description 13
- 238000001228 spectrum Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 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
- 230000008569 process Effects 0.000 description 1
- 238000005070 sampling Methods 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/288—Event detection in seismic signals, e.g. microseismics
-
- 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
- G01V1/362—Effecting static or dynamic corrections; Stacking
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
A complex field weak signal recovery method based on a vector bin and an electronic device are disclosed. The method can comprise the following steps: step 1: determining a target channel, and further constructing a vector surface element channel set corresponding to the target channel; step 2: performing horizontal superposition on the vector surface element gather to obtain a horizontal superposition channel; and step 3: performing Hilbert transform on the vector surface element gather, and calculating a complex field cosine phase function gather; and 4, step 4: calculating a weighting function according to the complex field cosine phase function gather; and 5: calculating a processing result of the target track according to the weighting function and the horizontal superposition track; step 6: and (5) replacing the target tracks, and repeating the steps 1-5 until all the target tracks are completed. The invention constructs the vector surface element channel set, restores the weak signal submerged by noise by adopting a complex field processing method, improves the signal-to-noise ratio of the prestack data on the basis of amplitude preservation, frequency preservation and phase preservation, meets the technical requirement of fidelity denoising, and improves the quality of the original data.
Description
Technical Field
The invention belongs to the technical field of seismic exploration digital signal processing, and particularly relates to a complex field weak signal recovery method based on a vector surface element and electronic equipment.
Background
The method has the disadvantages of severe change of terrain altitude difference, poor excitation receiving condition and the like in a complex seismic exploration block, the signal-to-noise ratio of the acquired seismic data is very low, and the existing pre-stack denoising method and technology are not completely applicable and have great limitation and inadaptability.
The prior pre-stack denoising technology mainly has three problems, namely, effective weak signals can be damaged when noise is removed; secondly, the fidelity is not enough, and side effects are easy to be brought; thirdly, effective fruits are gathered on the prestack tracks, but no obvious effect is generated on the superposed section, namely the prestack effect and the poststack effect are generated. Therefore, a prestack fidelity denoising technique is urgently needed in the prestack denoising technical field, and it is necessary to develop a complex field weak signal recovery method based on a vector bin and an electronic device.
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
In view of this, the embodiment of the present disclosure provides a vector bin-based complex field weak signal recovery method and an electronic device, where a vector bin gather is constructed, a complex field processing method is adopted to recover a weak signal submerged by noise, so as to improve the signal-to-noise ratio of prestack data on the basis of amplitude preservation, frequency preservation and phase preservation, meet the requirements of fidelity denoising technology, improve the quality of original data, achieve an obvious effect on both the prestack gather and the stacking section, and substantially solve the technical bottleneck problem of conventional prestack denoising, namely "prestack effective and poststack ineffective".
The embodiment of the disclosure provides a complex field weak signal recovery method based on a vector bin, which includes: step 1: determining a target track xm(t), further constructing a vector surface element gather corresponding to the target road; step 2: performing horizontal superposition on the vector surface element gather to obtain a horizontal superposition channel; and step 3: performing Hilbert transform on the vector surface element gather, and calculating a complex field cosine phase function gather; and 4, step 4: calculating a weighting function according to the complex field cosine phase function gather; and 5: calculating a processing result of the target track according to the weighting function and the horizontal superposition track; step 6: and (5) replacing the target tracks, and repeating the steps 1-5 until all the target tracks are completed.
Preferably, the step 1 comprises: with the target track xm(t) the CMP bin in which the target track x is located is taken as a central CMP binmAnd (t) taking the offset distance and the azimuth angle as standards, screening seismic channels with the offset distance and the azimuth angle consistent with the standards from the central CMP surface element and the adjacent CMP surface elements, and constructing the vector surface element channel set.
Preferably, the horizontal superposition lane is:
wherein k ism(t) is the mth target track xm(t) corresponding horizontal superimposed lanes, xi(t) is the vector bin gather, i is the track sequence number, t is the time series, n is the number of tracks or the number of coverage.
Preferably, the set of complex-field cosine-phase functions is calculated by equation (2):
wherein cos θi(t) is a complex field cosine phase function gather, ai(t) is the instantaneous envelope and,xiand (t) is a vector bin gather.
Preferably, the weighting function is calculated by equation (3):
wherein, gm(t) is the weighting function, cos θ, corresponding to the mth target tracki(t) is the complex field cosine phase function gather, and n is the number of tracks or the number of covering times.
Preferably, the processing result of the target track is calculated by formula (4):
ym(t)=gm(t)km(t) (4)
wherein, ym(t) is the processing result corresponding to the mth target track, gm(t) is the weighting function corresponding to the mth target track, kmAnd (t) is a horizontal superposition track corresponding to the mth target track.
As a specific implementation manner of the embodiment of the present disclosure, an embodiment of the present disclosure further provides an electronic device, including:
a memory storing executable instructions;
a processor executing the executable instructions in the memory to perform the steps of: step 1: determining a target track xm(t), further constructing a vector surface element gather corresponding to the target road; step 2: performing horizontal superposition on the vector surface element gather to obtain a horizontal superposition channel; and step 3: performing Hilbert transform on the vector surface element gather, and calculating a complex field cosine phase function gather; and 4, step 4: calculating a weighting function according to the complex field cosine phase function gather; and 5: calculating a processing result of the target track according to the weighting function and the horizontal superposition track; step 6: and (5) replacing the target tracks, and repeating the steps 1-5 until all the target tracks are completed.
Preferably, the step 1 comprises: with the target track xm(t) the CMP bin in which the target track x is located is taken as a central CMP binmAnd (t) taking the offset distance and the azimuth angle as standards, screening seismic channels with the offset distance and the azimuth angle consistent with the standards from the central CMP surface element and the adjacent CMP surface elements, and constructing the vector surface element channel set.
Preferably, the horizontal superposition lane is:
wherein k ism(t) is the mth target track xm(t) corresponding horizontal superimposed lanes, xi(t) is the vector bin gather, i is the track sequence number, t is the time series, n is the number of tracks or the number of coverage.
Preferably, the set of complex-field cosine-phase functions is calculated by equation (2):
wherein cos θi(t) is a complex field cosine phase function gather, ai(t) is the instantaneous envelope and,xiand (t) is a vector bin gather.
Preferably, the weighting function is calculated by equation (3):
wherein, gm(t) is the weighting function, cos θ, corresponding to the mth target tracki(t) is the complex field cosine phase function gather, and n is the number of tracks or the number of covering times.
Preferably, the processing result of the target track is calculated by formula (4):
ym(t)=gm(t)km(t) (4)
wherein, ym(t) is the processing result corresponding to the mth target track, gm(t) is the weighting function corresponding to the mth target track, kmAnd (t) is a horizontal superposition track corresponding to the mth target track.
Additional features and advantages of the invention will be set forth in the detailed description which follows.
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 throughout.
Fig. 1 shows a flow chart of the steps of a vector bin based complex field weak signal recovery method according to the present invention.
Fig. 2 shows a schematic diagram of vector binning construction according to one embodiment of the present invention.
Fig. 3a, 3b, and 3c respectively show a target trace, a set of vector bin traces, and a set of complex-domain cosine-phase function traces, according to an embodiment of the invention.
FIGS. 4a, 4b show a schematic of an original velocity spectrum and an original CMP gather, respectively, in accordance with an embodiment of the present invention.
FIGS. 5a, 5b show a velocity spectrum and a CMP gather after processing, respectively, in accordance with an embodiment of the present invention.
FIGS. 6a, 6b show schematic diagrams of a residual velocity spectrum and residual gathers before and after processing, respectively, according to an embodiment of the invention.
Fig. 7a, 7b show schematic diagrams of an original horizontal overlay cross section, a processed horizontal overlay cross section, respectively, according to an embodiment of the invention.
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.
Fig. 1 shows a flow chart of the steps of a vector bin based complex field weak signal recovery method according to the present invention.
In this embodiment, the vector bin-based complex field weak signal recovery method according to the present invention may include: step 1: determining a target track xm(t), further constructing a vector surface element gather corresponding to the target road; step 2: performing horizontal superposition on the vector surface element gather to obtain a horizontal superposition channel; and step 3: performing Hilbert transform on the vector surface element gather, and calculating a complex field cosine phase function gather; and 4, step 4: calculating a weighting function according to the complex field cosine phase function gather; and 5: calculating a processing result of the target track according to the weighting function and the horizontal superposition track; step 6: and (5) replacing the target tracks, and repeating the steps 1-5 until all the target tracks are completed.
In one example, step 1 comprises: with the target track xm(t) the CMP bin is the center CMP bin, and the target track x is usedmAnd (t) taking the offset distance and the azimuth angle as standards, screening seismic channels with the offset distance and the azimuth angle consistent with the standards from the central CMP surface element and the adjacent CMP surface elements, and constructing a vector surface element channel set.
In one example, the horizontal overlay lanes are:
wherein k ism(t) is the mth target track xm(t) the corresponding horizontal overlay channel, xi (t) is the vector bin gather, i is the channel sequence number, t is the time series, and n is the channel number or the coverage number.
In one example, a complex-domain cosine phase function gather is computed by equation (2):
wherein cos θi(t) is a complex field cosine phase function gather, ai(t) is the instantaneous envelope and,xiand (t) is a vector bin gather.
In one example, the weighting function is calculated by equation (3):
wherein, gm(t) is the weighting function, cos θ, corresponding to the mth target tracki(t) is the complex field cosine phase function gather, and n is the number of tracks or the number of covering times.
In one example, the processing result of the target track is calculated by equation (4):
ym(t)=gm(t)km(t) (4)
wherein, ym(t) is the processing result corresponding to the mth target track, gm(t) is the weighting function corresponding to the mth target track, kmAnd (t) is a horizontal superposition track corresponding to the mth target track.
Fig. 2 shows a schematic diagram of vector binning construction according to one embodiment of the present invention.
Specifically, based on the vector bin channel set, a complex number field processing method is used for protecting strong signals, recovering weak signals, improving the signal-to-noise ratio of a target channel, and meeting the requirements of a fidelity denoising technology.
The complex field weak signal recovery method based on the vector bin can comprise the following steps:
step 1: determining a target track xm(t) with target track xm(t) the CMP bin is the center CMP bin, and the target track x is usedm(t) the offset distance and the azimuth angle are standard, seismic channels with the offset distance and the azimuth angle consistent with the standard are screened out from the central CMP surface element and the adjacent CMP surface elements, and a vector surface element channel set is constructed, wherein the surface elements are distributed in a grid mode as shown in FIG. 2, and the adjacent surface elements are 8 surface elements around the central surface element; in the construction process, certain errors are allowed to exist between the offset distance and the azimuth angle and the standard, an error range is set, the offset distance error and the azimuth angle error are controlled by the error range, the larger the error is, the more the number of channels obtained by the vector surface element is, the coverage times are increased, the signal-to-noise ratio is favorably improved, but the fidelity is favorably improved, and similarly, the number of adjacent CMP surface elements, namely the size of the vector surface element is also controlled by parameters;
step 2: performing horizontal superposition on the vector bin gather to obtain a horizontal superposition channel as a formula (1);
and step 3: performing Hilbert transform on the vector surface element gather, and calculating a complex field cosine phase function gather through a formula (2); the complex field cosine phase function gather is independent of the amplitude, the amplitude of the weak signal is equal to that of the strong signal, and the weak signal is not divided into the strong signal and the weak signal, namely the weak signal is effectively strengthened;
and 4, step 4: calculating a weighting function according to a complex field cosine phase function gather through a formula (3);
and 5: calculating a processing result of the target track according to the weighting function and the horizontal superposition track through a formula (4);
step 6: and (5) replacing the target tracks, and repeating the steps 1-5 until all the target tracks are completed.
According to the method, a vector surface element gather is constructed, a complex number field processing method is adopted to restore a weak signal submerged by noise, the signal-to-noise ratio of prestack data is improved on the basis of amplitude preservation, frequency preservation and phase preservation, the technical requirements of fidelity denoising are met, the quality of original data is improved, obvious effects are achieved on the prestack gather and a stacking section, and the technical bottleneck problems of 'prestack effect and no poststack effect' in conventional prestack denoising are basically solved.
Application example
To facilitate understanding of the solution of the embodiments of the present invention and the effects thereof, a specific application example is 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.
Fig. 3a, 3b, and 3c respectively show a target trace, a set of vector bin traces, and a set of complex-domain cosine-phase function traces, according to an embodiment of the invention.
The complex field weak signal recovery method based on the vector bin comprises the following steps:
step 1: determining a target track x as shown in FIG. 3am(t) with target track xm(t) the CMP bin is the center CMP bin, and the target track x is usedm(t) the offset distance and azimuth angle are standard, seismic channels with the offset distance and azimuth angle consistent with the standard are screened out from the central CMP surface element and the adjacent CMP surface elements, and a vector surface element channel set is constructed, as shown in FIG. 3 b;
step 2: performing horizontal superposition on the vector bin gather to obtain a horizontal superposition channel as a formula (1);
and step 3: performing hilbert transform on the vector bin gather, and calculating a complex field cosine phase function gather by using a formula (2), as shown in fig. 3 c;
and 4, step 4: calculating a weighting function according to a complex field cosine phase function gather through a formula (3);
and 5: calculating a processing result of the target track according to the weighting function and the horizontal superposition track through a formula (4);
step 6: and (5) replacing the target tracks, and repeating the steps 1-5 until all the target tracks are completed.
Take a certain three-dimensional work area in the western China as an example.
The name of the work area: AB 3D; full coverage area: 203 square kilometers; sampling interval: 2 ms; CMP bin size: 25m × 25 m;
processing parameters of this embodiment:
1. vector bin size: 3 x 3;
the CMP surface element where the target channel is located is taken as the center, 3 CMP surface elements are arranged In the Inline direction, 3 CMP surface elements are arranged In the Cross line direction, and 9 CMP surface elements are counted.
2. Offset error: 25 m;
the offsets are considered consistent within an error of-25 m to +25m, based on the offset of the target track.
3. Azimuth error: 30 degrees;
the azimuth is considered to be uniform within an error of-30 degrees to +30 degrees, based on the azimuth of the target track.
Fig. 4a and 4b are schematic diagrams of an original velocity spectrum and an original CMP gather, respectively, according to an embodiment of the present invention, and the signal-to-noise ratio of the acquired seismic signals is very low due to poor excitation reception conditions in a desert region.
Fig. 5a, 5b show a schematic representation of a processed velocity spectrum and CMP gather, respectively, showing a significant improvement in signal-to-noise ratio after weak signal recovery and a clear view of the velocity spectrum effective wave energy mass, in accordance with an embodiment of the present invention.
Fig. 6a, 6b show schematic diagrams of a residual velocity spectrum before and after processing and a residual gather, respectively, without significant wave event axes within the residual gather and without any significant wave energy boluses in the velocity spectrum, in accordance with an embodiment of the present invention.
Fig. 7a, 7b show schematic diagrams of an original horizontal overlay cross section, a processed horizontal overlay cross section, respectively, according to an embodiment of the invention.
The comparison result shows that the method has obvious effect no matter before or after the stack and shows extremely high fidelity.
In summary, the invention constructs the vector bin gather, reduces the weak signal submerged by noise by adopting a complex field processing method, improves the signal-to-noise ratio of prestack data on the basis of amplitude preservation, frequency preservation and phase preservation, meets the technical requirement of fidelity denoising, improves the quality of original data, has obvious effects on the prestack gather and the stacking section, and basically solves the technical bottleneck problems of 'prestack effect and no poststack effect' of the conventional prestack denoising.
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.
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 to perform the steps of: step 1: determining a target track xm(t), further constructing a vector surface element gather corresponding to the target road; step 2: performing horizontal superposition on the vector surface element gather to obtain a horizontal superposition channel; and step 3: performing Hilbert transform on the vector surface element gather, and calculating a complex field cosine phase function gather; and 4, step 4: calculating a weighting function according to the complex field cosine phase function gather; and 5: calculating a processing result of the target track according to the weighting function and the horizontal superposition track; step 6: and (5) replacing the target tracks, and repeating the steps 1-5 until all the target tracks are completed.
In one example, the steps1 comprises: with the target track xm(t) the CMP bin is the center CMP bin, and the target track x is usedmAnd (t) taking the offset distance and the azimuth angle as standards, screening seismic channels with the offset distance and the azimuth angle consistent with the standards from the central CMP surface element and the adjacent CMP surface elements, and constructing a vector surface element channel set.
In one example, the horizontal overlay lanes are:
wherein k ism(t) is the mth target track xm(t) corresponding horizontal superimposed lanes, xi(t) is the vector bin gather, i is the track sequence number, t is the time series, n is the number of tracks or the number of coverage.
In one example, a complex-domain cosine phase function gather is computed by equation (2):
wherein cos θi(t) is a complex field cosine phase function gather, ai(t) is the instantaneous envelope and,xiand (t) is a vector bin gather.
In one example, the weighting function is calculated by equation (3):
wherein, gm(t) is the weighting function, cos θ, corresponding to the mth target tracki(t) is the complex field cosine phase function gather, and n is the number of tracks or the number of covering times.
In one example, the processing result of the target track is calculated by equation (4):
ym(t)=gm(t)km(t) (4)
wherein, ym(t)For the processing result corresponding to the mth target track, gm(t) is the weighting function corresponding to the mth target track, kmAnd (t) is a horizontal superposition track corresponding to the mth target track.
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.
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 (8)
1. A complex field weak signal recovery method based on vector bin is characterized by comprising the following steps:
step 1: determining a target track xm(t), further constructing a vector surface element gather corresponding to the target road;
step 2: performing horizontal superposition on the vector surface element gather to obtain a horizontal superposition channel;
and step 3: performing Hilbert transform on the vector surface element gather, and calculating a complex field cosine phase function gather;
and 4, step 4: calculating a weighting function according to the complex field cosine phase function gather;
and 5: calculating a processing result of the target track according to the weighting function and the horizontal superposition track;
step 6: replacing the target tracks, and repeating the steps 1-5 until all the target tracks are completed;
wherein the weighting function is calculated by equation (3):
wherein, gm(t) is the weighting function, cos θ, corresponding to the mth target tracki(t) is a complex field cosine phase function gather, n is the number of tracks or the number of covering times, i is a track serial number, and t is a time sequence;
wherein the processing result of the target track is calculated by formula (4):
ym(t)=gm(t)km(t) (4)
wherein, ym(t) is the processing result corresponding to the mth target track, gm(t) is the weighting function corresponding to the mth target track, kmAnd (t) is a horizontal superposition track corresponding to the mth target track.
2. The vector bin-based complex field weak signal recovery method according to claim 1, wherein the step 1 comprises:
with the target track xm(t) the CMP bin in which the target track x is located is taken as a central CMP binmAnd (t) taking the offset distance and the azimuth angle as standards, screening seismic channels with the offset distance and the azimuth angle consistent with the standards from the central CMP surface element and the adjacent CMP surface elements, and constructing the vector surface element channel set.
3. The vector bin-based complex field weak signal recovery method according to claim 1, wherein the horizontal superposition channel is:
wherein k ism(t) is the mth target track xm(t) corresponding horizontal superimposed lanes, xi(t) is the vector bin gather, i is the track sequence number, t is the time series, n is the number of tracks or the number of coverage.
4. The vector bin-based complex field weak signal recovery method according to claim 1, wherein the complex field cosine phase function gather is calculated by formula (2):
5. 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 perform the steps of:
step 1: determining a target track xm(t), further constructing a vector surface element gather corresponding to the target road;
step 2: performing horizontal superposition on the vector surface element gather to obtain a horizontal superposition channel;
and step 3: performing Hilbert transform on the vector surface element gather, and calculating a complex field cosine phase function gather;
and 4, step 4: calculating a weighting function according to the complex field cosine phase function gather;
and 5: calculating a processing result of the target track according to the weighting function and the horizontal superposition track;
step 6: replacing the target tracks, and repeating the steps 1-5 until all the target tracks are completed;
wherein the weighting function is calculated by equation (3):
wherein, gm(t) is the addition corresponding to the mth target trackWeight function, cos θi(t) is a complex field cosine phase function gather, n is the number of tracks or the number of covering times, i is a track serial number, and t is a time sequence;
wherein the processing result of the target track is calculated by formula (4):
ym(t)=gm(t)km(t) (4)
wherein, ym(t) is the processing result corresponding to the mth target track, gm(t) is the weighting function corresponding to the mth target track, kmAnd (t) is a horizontal superposition track corresponding to the mth target track.
6. The electronic device of claim 5, wherein the step 1 comprises:
with the target track xm(t) the CMP bin in which the target track x is located is taken as a central CMP binmAnd (t) taking the offset distance and the azimuth angle as standards, screening seismic channels with the offset distance and the azimuth angle consistent with the standards from the central CMP surface element and the adjacent CMP surface elements, and constructing the vector surface element channel set.
7. The electronic device of claim 5, wherein the horizontal overlay channel is:
wherein k ism(t) is the mth target track xm(t) corresponding horizontal superimposed lanes, xi(t) is the vector bin gather, i is the track sequence number, t is the time series, n is the number of tracks or the number of coverage.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911406030.1A CN111142155B (en) | 2019-12-31 | 2019-12-31 | Complex field weak signal recovery method based on vector surface element and electronic equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911406030.1A CN111142155B (en) | 2019-12-31 | 2019-12-31 | Complex field weak signal recovery method based on vector surface element and electronic equipment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111142155A CN111142155A (en) | 2020-05-12 |
CN111142155B true CN111142155B (en) | 2020-09-29 |
Family
ID=70522688
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911406030.1A Active CN111142155B (en) | 2019-12-31 | 2019-12-31 | Complex field weak signal recovery method based on vector surface element and electronic equipment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111142155B (en) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102053265A (en) * | 2009-10-30 | 2011-05-11 | 中国石油化工股份有限公司 | Filtering method for eliminating wild value of seismic data |
CN102053276A (en) * | 2009-10-30 | 2011-05-11 | 中国石油化工股份有限公司 | Two-dimensional filtering method for a plurality of gathers of digital seismic signal |
WO2015053811A1 (en) * | 2013-10-11 | 2015-04-16 | Chevron U.S.A. Inc. | System and method for regularizing seismic data |
CN104849756A (en) * | 2015-03-31 | 2015-08-19 | 中国地质大学(北京) | Method for improving resolution ratio of seismic data and enhancing energy of valid weak signals |
CN105093305A (en) * | 2014-05-15 | 2015-11-25 | 中国石油化工股份有限公司 | Method for determining position of reflection interface of seismic data |
CN105277985A (en) * | 2014-07-24 | 2016-01-27 | 中国石油化工股份有限公司 | OVT-domain seismic data regularization method based on image processing |
CN105510975A (en) * | 2015-12-31 | 2016-04-20 | 中国石油天然气集团公司 | Method and device for improving signal-to-noise ratio of seismic data |
CN105629306A (en) * | 2014-10-27 | 2016-06-01 | 中国石油化工股份有限公司 | Method for establishing signal-to-noise ratio(SNR) model |
CN106094029A (en) * | 2016-08-24 | 2016-11-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | The method utilizing offset distance vector sheet geological data Predicating Reservoir Fractures |
CN106199701A (en) * | 2016-07-15 | 2016-12-07 | 中国石油天然气集团公司 | The method for reconstructing of irregular geological data and device |
-
2019
- 2019-12-31 CN CN201911406030.1A patent/CN111142155B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102053265A (en) * | 2009-10-30 | 2011-05-11 | 中国石油化工股份有限公司 | Filtering method for eliminating wild value of seismic data |
CN102053276A (en) * | 2009-10-30 | 2011-05-11 | 中国石油化工股份有限公司 | Two-dimensional filtering method for a plurality of gathers of digital seismic signal |
WO2015053811A1 (en) * | 2013-10-11 | 2015-04-16 | Chevron U.S.A. Inc. | System and method for regularizing seismic data |
CN105093305A (en) * | 2014-05-15 | 2015-11-25 | 中国石油化工股份有限公司 | Method for determining position of reflection interface of seismic data |
CN105277985A (en) * | 2014-07-24 | 2016-01-27 | 中国石油化工股份有限公司 | OVT-domain seismic data regularization method based on image processing |
CN105629306A (en) * | 2014-10-27 | 2016-06-01 | 中国石油化工股份有限公司 | Method for establishing signal-to-noise ratio(SNR) model |
CN104849756A (en) * | 2015-03-31 | 2015-08-19 | 中国地质大学(北京) | Method for improving resolution ratio of seismic data and enhancing energy of valid weak signals |
CN105510975A (en) * | 2015-12-31 | 2016-04-20 | 中国石油天然气集团公司 | Method and device for improving signal-to-noise ratio of seismic data |
CN106199701A (en) * | 2016-07-15 | 2016-12-07 | 中国石油天然气集团公司 | The method for reconstructing of irregular geological data and device |
CN106094029A (en) * | 2016-08-24 | 2016-11-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | The method utilizing offset distance vector sheet geological data Predicating Reservoir Fractures |
Non-Patent Citations (3)
Title |
---|
Offset vector tile domain wide-azimuth prestack seismic interpretation methods and applications;Zhan Shifan 等;《INTERPRETATION》;20180501;第6卷(第2期);第T337-T347页 * |
OVT域地震数据规则化技术及应用;李博;《石油物探》;20190131;第58卷(第1期);第53-62页 * |
基于曲波域稀疏约束的OVT域地震数据去噪方法研究;孙苗苗 等;《石油物探》;20190331;第58卷(第2期);第208-218页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111142155A (en) | 2020-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1257850B1 (en) | Processing seismic data | |
CN109001813B (en) | Method, device and system for suppressing multiple waves | |
CN109164483A (en) | Multi-component earthquake data vector denoising method and multi-component earthquake data vector denoise device | |
CN104932010B (en) | A kind of diffracted wave separation method based on the sparse Radon transformation of shortcut fringing | |
CN111880219B (en) | Diffracted wave imaging method and device based on azimuth-dip angle gather | |
CN102854533A (en) | Wave field separation principle based denoising method for increasing signal to noise ratio of seismic data | |
CN107942327A (en) | Single channel HRWS SAR imaging methods based on impulse phase coding | |
CN106094046A (en) | Time domain aviation electromagnetic data denoising method based on singular value decomposition and wavelet analysis | |
CN107942389A (en) | For suppressing method, system and the computer-readable medium of adjacent big gun interference | |
CN106019374B (en) | The tomograph imaging method of similarity is dissipated based on reflection groove wave frequency | |
CN107861156B (en) | The extracting method and device of diffracted wave | |
CN105301648A (en) | Method of acquiring common reflection surface stacking parameters | |
CN105510975A (en) | Method and device for improving signal-to-noise ratio of seismic data | |
CN103412325A (en) | Method for separating and denoising wave fields of p-waves and s-waves of earthquake | |
CN111142155B (en) | Complex field weak signal recovery method based on vector surface element and electronic equipment | |
WO2006111543A1 (en) | Method of processing seismic data for avo or avoa characterisation | |
CN107450103A (en) | A kind of terrible ripple drawing method based on boundary integral inverse operator | |
CN110618449B (en) | Method and system for processing seismic data | |
CN103217677A (en) | Single-channel SAR (synthetic aperture radar) moving target detecting method based on joint detection amount | |
CN101840001B (en) | Acquiring method and device of geological structure three-dimensional imaging data | |
CN104155688A (en) | High precision weighted stack method | |
CN102338884B (en) | Elliptic window direction band-pass amplitude-preserving filtering data processing method in geophysical prospecting | |
CN104614725B (en) | Scanning synthetic aperture radar image quality improving method and device | |
CN104765069A (en) | Method of suppressing adjacent shot interference in case of synchronous excitation acquisition | |
CN103513288A (en) | Compensation direction filtering method of two-dimensional grid 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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |