CN114076985B - Method, device, equipment and storage medium for suppressing noise of controllable vibration source - Google Patents
Method, device, equipment and storage medium for suppressing noise of controllable vibration source Download PDFInfo
- Publication number
- CN114076985B CN114076985B CN202010834408.4A CN202010834408A CN114076985B CN 114076985 B CN114076985 B CN 114076985B CN 202010834408 A CN202010834408 A CN 202010834408A CN 114076985 B CN114076985 B CN 114076985B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- ovt
- noise
- gather
- channel
- 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 72
- 238000003860 storage Methods 0.000 title claims abstract description 13
- 230000001629 suppression Effects 0.000 claims abstract description 48
- 238000003825 pressing Methods 0.000 claims abstract description 46
- 230000000007 visual effect Effects 0.000 claims abstract description 32
- 238000001228 spectrum Methods 0.000 claims abstract description 29
- 238000006243 chemical reaction Methods 0.000 claims abstract description 14
- 230000006835 compression Effects 0.000 claims abstract description 12
- 238000007906 compression Methods 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims description 46
- 238000004590 computer program Methods 0.000 claims description 19
- 230000006870 function Effects 0.000 claims description 14
- 238000011068 loading method Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 8
- 238000004519 manufacturing process Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 2
- 210000001233 cdp Anatomy 0.000 description 108
- 238000004637 computerized dynamic posturography Methods 0.000 description 108
- 230000008569 process Effects 0.000 description 13
- 230000008901 benefit Effects 0.000 description 5
- 230000009467 reduction Effects 0.000 description 5
- 238000010183 spectrum analysis Methods 0.000 description 5
- 229910052704 radon Inorganic materials 0.000 description 4
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000002708 enhancing effect Effects 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 description 2
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 241000220317 Rosa Species 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011049 filling Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000000153 supplemental effect Effects 0.000 description 1
- 238000012360 testing method 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/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
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)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a method, a device, equipment and a storage medium for suppressing noise of a controllable vibration source, wherein the method comprises the following steps: establishing a five-dimensional offset vector chip OVT gather corresponding to a current pressed noise channel by taking the current pressed noise channel as a center; analyzing the five-dimensional spectrum of the five-dimensional OVT gather to obtain a visual inclination angle of the five-dimensional OVT gather; superposing the five-dimensional OVT gathers along the visual inclination angle to obtain a noise pressing model gather; and performing noise suppression on the initial domain seismic data corresponding to the noise compression model channel by utilizing conversion among different domains to obtain target seismic data corresponding to a controllable seismic source, so that effective signals are better recovered while effective suppression of noise in a black triangle area is realized, and meanwhile, effective signals outside the black triangle area are enhanced, and the signal-to-noise ratio of data is improved.
Description
Technical Field
The invention belongs to the technical field of exploration, and particularly relates to a method, a device, equipment and a storage medium for suppressing noise of a controllable seismic source.
Background
The controllable vibration source collection can not avoid noise such as vibration source mechanical interference, harmonic interference, scattering interference and the like. In particular to the mechanical interference of a seismic source, the scattering interference seriously reduces the signal-to-noise ratio of data, and the difficulty in realizing the exploration target is increased especially for ultra-deep exploration with thicker sand layers on the surface. For vibroseis noise, there are two main causes of its generation: firstly, interference caused by nonlinearity of a controllable vibration source mechanical device, a vibration device and a hydraulic servo system; and secondly, harmonic interference generated by the coupling effect of the vibration plate and the ground. The harmonic wave and the strong scattering noise in the desert area have strong energy, wide frequency band, wide time-space domain distribution range, poor regularity and difficult suppression.
In the 70 s of the 20 th century, the acquisition of controllable vibration sources enters the conventional production stage abroad, and a synchronous excitation pilot test of a plurality of vibration sources and different shots is carried out in 1979 by Sliverman et al, so that the method is mainly applied to a desert area.
The collection mode of the domestic controllable seismic source starts relatively late, but has been rapidly developed in recent years. In 2004, eastern geophysical companies in China have developed a small number of controllable seismic source acquisitions in Gobi, and in 2014, attempts have been made to apply the methods in the abdominal desert area. Because of the advantages of safe and environment-friendly controllable vibration source, high construction efficiency, low cost and the like, oriental companies gradually form a series of controllable vibration source efficient acquisition technologies such as alternate scanning, sliding scanning synchronous excitation, independent synchronous scanning and the like by conventional acquisition, and the controllable vibration source system is widely applied to regions such as deserts, gobi, mountain front zones and the like at home and abroad. Petrochemical industry in 2015 adopts controllable seismic source to collect in the Tarim desert area for the first time, and good effect is obtained.
Seriff et al first discussed the problem of higher harmonic distortion that occurs when the vibroseis coupled to the ground. Then Sallas analyzes the controllable source signal of the well logging receiving point and combines theoretical deduction to further study the higher harmonic distortion generated when the controllable source panel is coupled with the ground, and a method for inhibiting harmonic distortion generation by adopting panel acceleration information is provided, so that the reliability and stability of excitation signals are improved. However, the high energy scattering noise caused by the continuous vibration is not well pressed. The existing processing means adopt a filtering method, a correlation method, a model method, a statistical method and the like to press according to the characteristics of amplitude energy difference, time-frequency relation, correlation and the like, but effective signals cannot be recovered well when noise in a black triangle is pressed.
Therefore, how to achieve suppression of noise in a black triangle and better recovery of effective signals is a technical problem to be solved by those skilled in the art.
Disclosure of Invention
The invention mainly aims to provide a method, a device, equipment and a storage medium for suppressing noise of a controllable vibration source, so that effective signals can be better recovered while suppressing noise in a black triangle.
In view of the above problems, the present invention provides a method for suppressing noise of a vibroseis, including:
establishing a five-dimensional offset vector chip OVT gather corresponding to a current pressed noise channel by taking the current pressed noise channel as a center;
analyzing the five-dimensional spectrum of the five-dimensional OVT gather to obtain a visual inclination angle of the five-dimensional OVT gather;
superposing the five-dimensional OVT gathers along the visual inclination angle to obtain a noise pressing model gather;
and performing noise suppression on the initial domain seismic data corresponding to the noise compression model channels by utilizing conversion among different domains to obtain target seismic data corresponding to the controllable seismic source.
Further, in the method for suppressing noise of a vibroseis described above, the five-dimensional OVT gather includes dimensions including offset, azimuth, common depth point CDP, CDP line and time;
The establishing a five-dimensional OVT gather corresponding to the current pressed noise channel by taking the current pressed noise channel as a center comprises the following steps:
selecting a first set number of CDP lines and a second set number of CDP lines by taking the current pressed noise channel as a center to form an original three-dimensional CDP channel set with the time;
setting the offset and the azimuth;
the five-dimensional OVT gather is generated based on the offset, the azimuth, and the original three-dimensional CDP gather.
Further, in the method for suppressing noise of a vibroseis described above, generating the five-dimensional OVT gather based on the offset, the azimuth, and the original three-dimensional CDP gather includes:
forming an OVT bin based on the offset and the azimuth;
constructing an initial five-dimensional OVT gather based on the OVT surface elements;
performing the following operations on each lane in the original three-dimensional CDP lane set to produce the five-dimensional OVT lane set;
detecting whether a current track of the original three-dimensional CDP track set belongs to the initial five-dimensional OVT track set;
if the current channel of the original three-dimensional CDP channel set belongs to the initial five-dimensional OVT channel set, detecting whether the initial five-dimensional OVT grid corresponding to the current channel is in an idle state or not;
If the initial five-dimensional OVT grid corresponding to the current channel is in an idle state, loading the current channel into the initial five-dimensional OVT grid so as to update the initial five-dimensional OVT channel set;
and if the initial five-dimensional OVT grid corresponding to the current track is in an occupied state, selecting the occupied track or loading the current track to the initial five-dimensional OVT grid so as to update the initial five-dimensional OVT track set.
Further, in the method for suppressing noise of a vibroseis, before the analyzing the five-dimensional spectrum of the five-dimensional OVT gather to obtain the view angle of the five-dimensional OVT gather, the method further includes:
and manufacturing the five-dimensional spectrum by a dip angle scanning method.
Further, in the method for suppressing noise of a vibroseis, the step of superposing the five-dimensional OVT gather along the view inclination angle to obtain a noise-suppressing model gather includes:
and superposing the five-dimensional OVT gathers along the visual inclination angle by using a mean superposition calculation formula to obtain a noise pressing model gather.
Further, the method for suppressing noise of the controllable vibration source further includes:
if the five-dimensional OVT channel set does not have singular amplitude values, selecting a first calculation formula as the mean value superposition calculation formula;
If the five-dimensional OVT channel set has singular amplitude values, selecting a second calculation formula as the mean value superposition calculation formula;
the first calculation formula is
Wherein, N is the total number of g (x, y, t), g (x, y, t) is a two-dimensional amplitude function, x and y are the coordinate positions of g in CDP and line direction;
the second calculation formula is
Further, in the method for suppressing noise of a vibroseis described above, performing noise suppression on the initial domain seismic data corresponding to the noise suppression model channel by using conversion between different domains to obtain target seismic data corresponding to the vibroseis, including:
converting the initial domain seismic data corresponding to the noise pressing model channel into target domain seismic data of a target domain;
performing noise suppression processing on the target domain seismic data to obtain noise-reduced seismic data;
converting the noise-reduced seismic data to a time domain to obtain supplementary seismic data;
and combining the initial domain seismic data with the supplementary seismic data to obtain target seismic data corresponding to the controllable seismic source.
The invention also provides a noise suppression device of the controllable vibration source, which comprises:
the building module is used for building a five-dimensional offset vector chip OVT gather corresponding to the current pressed noise channel by taking the current pressed noise channel as a center;
The analysis module is used for analyzing the five-dimensional spectrum of the five-dimensional OVT gather to obtain the visual inclination angle of the five-dimensional OVT gather;
the superposition module is used for superposing the five-dimensional OVT gathers along the visual inclination angle to obtain a noise pressing model gather;
and the noise pressing module is used for performing noise pressing on the initial domain seismic data corresponding to the noise pressing model channels by utilizing conversion among different domains to obtain the target seismic data corresponding to the controllable seismic source.
The invention also provides a noise suppression device of the controllable vibration source, which comprises a memory and a controller;
the memory has stored thereon a computer program which, when executed by a controller, implements the steps of the method as described in any of the above.
The invention also provides a storage medium having stored thereon a computer program which, when executed by a controller, implements the steps of the method as described in any of the preceding claims.
One or more embodiments of the above-described solution may have the following advantages or benefits compared to the prior art:
according to the method, the device, the equipment and the storage medium for suppressing the noise of the controllable vibration source, the five-dimensional OVT gather corresponding to the current noise pressed channel is established by taking the current noise pressed channel as the center; analyzing the five-dimensional spectrum of the five-dimensional OVT gather to obtain a visual inclination angle of the five-dimensional OVT gather; and superposing the five-dimensional OVT trace sets along the visual inclination angle to obtain a noise pressing model trace, performing noise pressing on the initial domain seismic data corresponding to the noise pressing model trace by utilizing conversion among different domains, and recovering effective data to obtain target seismic data corresponding to a controllable seismic source, thereby realizing effective pressing on noise of a black triangle area, better recovering effective signals, enhancing effective signals outside the black triangle area, and improving the signal-to-noise ratio of data.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by 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 invention, and are incorporated in and constitute a part of this specification, illustrate the invention and together with the embodiments of the invention, serve to explain the invention, without limitation to the invention. In the drawings:
FIG. 1 is a flow chart of an embodiment of a method for suppressing noise of a vibroseis according to the present invention;
FIG. 2 is a flow chart of step 103 in FIG. 1;
FIG. 3 is a graph comparing CMP gather records before and after noise suppression in a region of northwest China;
FIG. 4 is an enlarged view showing the comparison of CMP gathers before and after noise suppression in a region of northwest;
FIG. 5 is a diagram showing the contrast of the internal and external spectrums of a black triangle after noise reduction in a region in the northwest;
FIG. 6 is a comparison of seismic superimposed sections before and after noise suppression in a region of northwest China;
FIG. 7 is an enlarged display contrast diagram of an earthquake superposition section before and after noise suppression in a certain region in the northwest;
fig. 8 is a schematic structural view of an embodiment of a noise suppressing apparatus for a vibroseis according to the present invention.
Detailed Description
The following will describe embodiments of the present invention in detail with reference to the drawings and examples, thereby solving the technical problems by applying technical means to the present invention, and realizing the technical effects can be fully understood and implemented accordingly. It should be noted that, as long as no conflict is formed, each embodiment of the present invention and each feature of each embodiment may be combined with each other, and the formed technical solutions are all within the protection scope of the present invention.
Example 1
In order to solve the technical problems in the prior art, the embodiment of the invention provides a noise suppression method for a controllable vibration source.
Fig. 1 is a flowchart of an embodiment of a method for suppressing noise of a vibroseis according to the present invention, as shown in fig. 1, the method for suppressing noise of a vibroseis according to the present embodiment may specifically include the following steps:
100. establishing a five-dimensional offset vector patch (Offset Vector Tile, OVT) gather corresponding to a current pressed noise channel by taking the current pressed noise channel as a center;
in this embodiment, the five-dimensional OVT gather includes dimensions including offset, azimuth, common depth point (Common Depth Point, CDP), CDP line, and time.
In one specific implementation, for a common shot gather, each seismic trace is placed in a CDP grid (CDP bin) according to the location of the offset coordinate center point, and a pre-stack three-dimensional CDP gather is arranged.
In each CDP trace set, if the near offset trace is compared to the far offset trace, the correlation is not necessarily good, although the reflection points are very close. If two adjacent offset tracks are compared, the correlation between them is much better. The correlation is better if the azimuth angles of the two tracks are also similar.
The rose is made according to the position distribution density of the center point of the offset coordinate, and then concentric circles (offset distance) with equal intervals and bisectors (azimuth angles) passing through the center of the circle are drawn, so that a plurality of sector-shaped surface elements are formed. One fan-shaped bin is an OVT bin.
OVT bins are defined by azimuth and offset, with two dimensions. An OVT bin, corresponding to a three-dimensional CDP gather (CDP line, CDP, time), is selected from the original three-dimensional CDP gather based on azimuth and offset. All OVT cells are combined to form a five-dimensional OVT body. By such decomposition, the pre-stack three-dimensional CDP gather can be changed into the five-dimensional OVT volume needed by us. It features that adjacent tracks have better correlation than original prestack three-dimensional CDP track set. This is critical to the success of five-dimensional noise suppression.
In the noise pressing process, each pressed noise channel needs to form a five-dimensional OVT channel set taking the pressed noise channel as the center. The size of the OVT volume depends on the user-specified noise suppression operator. The grid of the noise pressing operator is defined by the offset number No, the azimuth number Na, the number Nx of CDPs and the number Ny of lines.
Assuming that the current noise-suppressed trace is offset bin number ko, azimuth bin number ka, CDP number kx, line number ky, and centered on the trace, the five-dimensional OVT trace set (gOVT (ko, ka, kx, ky, t)) consists of:
CDP number: -Nx/2, -Nx/2-1, …, kx, …, nx/2-1, nx/2;
line number: -Ny/2, -Ny/2-1, …, ky, …, ny/2-1, ny/2;
offset bin number: -No/2, -No/2-1, …, ko, …, no/2-1, no/2;
azimuth bin number: na/2, -Na/2-1, …, ka, …, na/2-1, na/2.
In this embodiment, the current noise-pressed channel may be used as a center, and a first set number of CDPs and a second set number of CDP lines may be selected to form an original three-dimensional CDP gather with the time; and after setting the offset and the azimuth, generating the five-dimensional OVT gather based on the offset, the azimuth and the original three-dimensional CDP gather. In this way, the five-dimensional OVT gather is generated by the original three-dimensional CDP gather, so that the search space can be accurately limited, and the search efficiency is improved.
In practical applications, special attention is required to: 1) The number of tracks contained in the original three-dimensional CDP track set has a considerable part and does not belong to the five-dimensional OVT track set; 2) Some grid points in the five-dimensional OVT gather may have duplicate channels in the original three-dimensional CDP gather; 3) Some grid points in the five-dimensional OVT gather may not exist in the original three-dimensional CDP gather. In this embodiment, the process of performing trace gather extraction in the original three-dimensional CDP trace set cannot adopt conventional trace sorting, and the relatively economical and applicable method is a five-dimensional grid filling method.
Specifically, based on the offset, the azimuth, and the original three-dimensional CDP gather, the process of generating the five-dimensional OVT gather is as follows:
a. forming an OVT bin based on the offset and the azimuth;
b. constructing an initial five-dimensional OVT gather based on the OVT surface elements;
since one OVT bin corresponds to one three-dimensional CDP trace set, in this embodiment, all OVT bins may be combined to form an initial five-dimensional OVT trace set.
c. Detecting whether a current track of the original three-dimensional CDP track set belongs to the initial five-dimensional OVT track set; if yes, executing the step d, and if not, executing the step g;
specifically, whether the initial five-dimensional OVT gather includes the current track of the original three-dimensional CDP gather may be determined according to the identifier of each track in the initial five-dimensional OVT gather and the identifier of the current track of the original three-dimensional CDP gather, if the initial five-dimensional OVT gather includes the current track of the original three-dimensional CDP gather, it is indicated that the current track of the original three-dimensional CDP gather belongs to the initial five-dimensional OVT gather, otherwise, if the initial five-dimensional OVT gather does not include the current track of the original three-dimensional CDP gather, it is indicated that the current track of the original three-dimensional CDP gather does not belong to the initial five-dimensional OVT gather.
d. Detecting whether the initial five-dimensional OVT grid corresponding to the current channel is in an idle state or not; if yes, executing the step e, and if not, executing the step f;
if the current track of the original three-dimensional CDP track set belongs to the initial five-dimensional OVT track set, it is indicated that the current track of the original three-dimensional CDP track set is valid, but the corresponding content may already be added to the initial five-dimensional OVT grid, so in this embodiment, it may be detected whether the initial five-dimensional OVT grid corresponding to the current track of the original three-dimensional CDP track set is in an idle state.
e. Loading a current track of the original three-dimensional CDP track set into the initial five-dimensional OVT grid so as to update the initial five-dimensional OVT track set;
if the initial five-dimensional OVT grid corresponding to the current track of the original three-dimensional CDP gather is in an idle state, at this time, the current track of the original three-dimensional CDP gather may be loaded into the initial five-dimensional OVT grid to update the initial five-dimensional OVT gather.
f. Selecting an occupied channel or the current channel to be loaded to the initial five-dimensional OVT grid so as to update the initial five-dimensional OVT channel set;
and if the initial five-dimensional OVT grid corresponding to the current track of the original three-dimensional CDP track set is in an occupied state, selecting an occupied track or loading the current track to the initial five-dimensional OVT grid so as to update the initial five-dimensional OVT track set.
Specifically, the current trace of the original three-dimensional CDP trace set may be geometrically compared with the occupied trace in the initial five-dimensional OVT grid to determine which trace is good and which trace is selected again. For example, the amplitude may be used as a comparison target, and if the amplitude is large, this means that the track is relatively good, whereas if the amplitude is small, this means that the track is relatively bad.
g. Skipping the current track of the original three-dimensional CDP gather.
If the current channel of the original three-dimensional CDP channel set does not belong to the original five-dimensional OVT channel set, the current channel of the original three-dimensional CDP channel set can be directly skipped, next extraction is carried out, and the needed five-dimensional OVT channel set can be obtained after all channels in the original three-dimensional CDP channel set are traversed.
101. Analyzing the five-dimensional spectrum of the five-dimensional OVT gather to obtain a visual inclination angle of the five-dimensional OVT gather;
specifically, in the process of creating the noise-pressing model trace, it is necessary to obtain the apparent dip angles of the five-dimensional OVT trace set in four directions caused by dip angles, tendencies, bends, and the like of the stratum. This requires five-dimensional spectrum analysis to obtain the maximum correlation spectrum of each data on each sample point and the corresponding viewing inclination angles of the offset, azimuth angle, line and CDP directions. And the reliable five-dimensional noise pressing model channel for denoising can be established by superposition along the optimal visual inclination angle direction. In order to obtain the visual tilt angles, the program uses a tilt angle scanning method to manufacture a five-dimensional spectrum of a five-dimensional OVT gather, and obtains the visual tilt angles through spectrum analysis.
The scanning range of the inclination angle is defined by the time difference (millisecond) of the adjacent channels, and is converted into the scanning number ndipo in the offset direction, the scanning number ndipa in the azimuth direction, the scanning number nditxl in the CDP direction and the scanning number ndipsl in the CDP line direction in actual operation. For each sample point in the seismic trace, the total number of spectral values NS is: ns=ndipo ndipsla ndipsl
In this way, most industrial software calculates the similarity coefficient ECC of each spectrum value by using the correlation spectrum and the superposition spectrum, and selects (a, o, y, x, t) corresponding to the selected maximum similarity coefficient ECC value as the best view angle, in this embodiment, a normalized correlation spectrum is selected for calculation, and the calculation formula (1) corresponding to the normalized correlation spectrum is as follows:
ECC represents a similarity coefficient, and the optimal inclination angle direction composed of the optimal offset, the optimal azimuth angle, the optimal line and the optimal track direction is obtained by using a similarity coefficient discrimination criterion. In the calculation formula (1), the numerator is the square of the sum of the amplitudes, and the denominator is the sum of the squares of the number of times of coverage M and the amplitudes. fa, o, y, x, t represent amplitude values at azimuth angle a, offset o, line direction y, CDP direction x, time t after a certain apparent tilt correction.
102. Superposing the five-dimensional OVT gathers along the visual inclination angle to obtain a noise pressing model gather;
in this embodiment, a five-dimensional noise compression model trace is manufactured, and the five-dimensional OVT trace set is required to be superimposed in the five-dimensional OVT trace set along the optimal viewing angle direction by using a mean value superposition calculation method, so as to obtain the noise compression model trace. Because the superposition is done on a sample-by-sample basis, taking one sample from a five-dimensional OVT gather degrades into a four-dimensional volume. The offset and azimuth of the stacking process need not be further considered, and these two dimensions can be removed and projected onto a two-dimensional plane consisting of the CDP and line number to form a two-dimensional amplitude function g (x, y, t) for the current point in time t, where x and y are the coordinate positions of g in the CDP and line directions.
In actual data, there are often some singular amplitude values, which are characterized by their being much larger than normal amplitude values. Therefore, in order to aim at the embodiment of the singular amplitude value and the non-singular amplitude value, different calculation methods can be adopted to overlap the five-dimensional OVT gather along the view inclination angle, so as to obtain a five-dimensional press noise model gather.
If no singular amplitude value exists in the five-dimensional OVT channel set, the first calculation formula (2) can be selected as a mean superposition calculation formula:
Where N is the total number of g (x, y, t), g (x, y, t) is a two-dimensional amplitude function, and x and y are the coordinate positions of g in CDP and line directions.
If the singular amplitude values exist in the five-dimensional OVT channel set, the second calculation formula (3) can be selected as a mean superposition calculation formula:
wherein gmite (x, y, t) is the amplitude value of g (x, y, t) after the singular values are removed. d is the distance of the superimposed track to the current track, ipow is the distance weighting factor, and p is the amplitude normalization factor.
In addition, in the stacking process in this embodiment, the stacking of one CMP gather at a certain time point is to stack data of different offset distances and azimuth angles into one line, and the five-dimensional OVT gather in the embodiment also involves different lines and tracks, which is equivalent to stacking a plurality of CMP from different lines and tracks.
103. And performing noise suppression on the initial domain seismic data corresponding to the noise compression model channels by utilizing conversion among different domains to obtain target seismic data corresponding to the controllable seismic source.
Specifically, noise suppression can be performed on the initial domain seismic data corresponding to the noise suppression model channel according to the next process of fig. 2, so as to obtain the target seismic data corresponding to the controllable seismic source. Fig. 2 is a flowchart of step 103 in fig. 1.
1031. Converting the initial domain seismic data corresponding to the noise pressing model channel into target domain seismic data of a target domain;
in one embodiment, the same data is observed from different angles, with different results. For single-pass data, the place where the waveform is present may be a reflected signal, or may be a direct wave, a refracted wave, a plane wave, a multiple, etc., which are limited bandwidth signals. However, the signals have random noise interference therein, and the signals and the random noise mixed therein are indiscriminate in the time domain. If it is transformed into the frequency domain, something different is seen. The bandwidth of the signal is generally limited, while the bandwidth of random noise is almost infinite. By filtering, a lot of noise can be filtered out, even if the noise is mixed in the signal in the time domain, it is partly suppressed.
In the case of two-dimensional data, such as signals from the same cannon from detectors arranged on the same line, we see a result that is different from the time-space domain if transformed into the frequency-wavenumber domain, or Tau-P domain. Many methods of suppressing noise are implemented based on different domains.
Therefore, in this embodiment, the initial domain seismic data corresponding to the noise compression model trace may be converted into the target domain seismic data of the target domain. The target domain is preferably a frequency wave number domain or a Tau-P domain.
Let g (x, t) be a two-dimensional data volume, x is the spatial position and t is the time, it can change g (x, t) from the two-dimensional spatial time domain to the frequency-wavenumber domain by two-dimensional fourier transform according to the calculation formula (4):
where ω is the frequency, k is the wave number, G (ω, k) is the data volume of the frequency-wave number domain corresponding to G (x, t).
G (x, t) can also be transformed from the two-dimensional spatial time domain to the Tau-P domain according to equation (5):
where p is slowness, τ is intercept time, n is a distance weighting factor, if a linear Radon transform is chosen, f (t, x, τ, p) is defined as: f (t, x, τ, p) =t- τ -px;
the Delta function is defined as:
if hyperbolic Radon transform is chosen, f (t, x, τ, p) is defined as: f (t, x, τ, p) =t- τ -px 2 ；
The definition of the Delta function is：
1032. Performing noise suppression processing on the target domain seismic data to obtain noise-reduced seismic data;
after the initial domain seismic data corresponding to the noise compression model channel is converted into the target domain seismic data of the target domain, the noise reduction seismic data can be obtained by compressing or cutting off the target domain.
1033. Converting the noise-reduced seismic data to a time domain to obtain supplementary seismic data;
After the noise reduction seismic data are obtained, the noise reduction seismic data can be converted into a time domain, and effective signals are recovered to obtain supplementary seismic data.
Specifically, if the target domain is the frequency-wave number domain, the inverse transformation may be performed by using the calculation formula (6), and the noise-reduced seismic data may be converted into the time domain.
If the target domain is a Tau-P domain, inverse transformation can be performed by using the calculation formula (7), and the noise reduction seismic data can be converted into a time domain.
If a linear radon transform is chosen, b (t, x, τ, p) is defined as: b (t, x, τ, p) =τ -t+px;
the Delta function is defined as:
if hyperbolic Radon transform is chosen, b (t, x, τ, p) is defined as: b (t, x, τ, p) =τ -t+px 2 ；
The Delta function is defined as:
1034. and combining the initial domain seismic data with the supplementary seismic data to obtain target seismic data corresponding to the controllable seismic source.
In this embodiment, the initial domain seismic data and the supplemental seismic data are combined to obtain the target seismic data corresponding to the controllable seismic source.
According to the noise suppression method of the controllable vibration source, a five-dimensional OVT gather corresponding to a current pressed noise channel is established by taking the current pressed noise channel as a center; analyzing the five-dimensional spectrum of the five-dimensional OVT gather to obtain a visual inclination angle of the five-dimensional OVT gather; and superposing the five-dimensional OVT trace sets along the visual inclination angle to obtain a noise pressing model trace, performing noise pressing on the initial domain seismic data corresponding to the noise pressing model trace by utilizing conversion among different domains, and recovering effective data to obtain target seismic data corresponding to a controllable seismic source, thereby realizing effective pressing on noise of a black triangle area, better recovering effective signals, enhancing effective signals outside the black triangle area, and improving the signal-to-noise ratio of data.
It should be noted that, the method of the embodiment of the present invention may be performed by a single device, for example, a computer or a server. The method of the embodiment can also be applied to a distributed scene, and is completed by mutually matching a plurality of devices. In the case of such a distributed scenario, one of the devices may perform only one or more steps of the method of an embodiment of the present invention, and the devices interact with each other to complete the method.
Example two
The present embodiment describes the technical solution of the present invention with specific examples. According to the method and the device for acquiring the controllable seismic source in the northwest desert area in China, the black triangle noise can be effectively suppressed, and effective signals in the black triangle area can be better recovered. Fig. 3 is a comparison chart of CMP gather records before and after noise suppression in a northwest region, fig. 4 is an enlarged display comparison chart of CMP gathers before and after noise suppression in a northwest region, fig. 5 is a comparison chart of internal and external frequency spectrums of black triangles after noise suppression in a northwest region, and fig. 6 is a comparison chart of earthquake superposition sections before and after noise suppression in a northwest region; FIG. 7 is an enlarged view showing the cross section of the seismic superposition before and after noise suppression in a region of northwest China.
As shown in fig. 3-4, the effective signals in the black triangle area are consistent with the effective signals outside the black triangle area in phase, and the signal to noise ratio of the gather is greatly improved. Fig. 5 is a spectral analysis of the same time window inside and outside the suppressed black triangle, showing nearly uniform spectral morphology, which verifies that the effective signal is recovered in the black triangle after noise suppression. The cross-section comparison results are superimposed before and after denoising as shown in fig. 6. It can be seen that the overall signal to noise ratio of the profile is improved. In fig. 7, the middle and deep layers are displayed in an enlarged manner, so that the same phase axis is more continuous, the resolution is higher, and the internal horizon of the otto system can be tracked continuously.
Example III
In order to solve the technical problems in the prior art, the embodiment of the invention also provides a noise suppression device for the controllable vibration source.
Fig. 8 is a schematic structural diagram of an embodiment of a noise suppression device for a vibroseis according to the present invention, and as shown in fig. 8, the noise suppression device for a vibroseis according to the present embodiment includes a building module 80, a parsing module 81, a superposition module 82, and a noise suppression module 83.
The establishing module 80 is configured to establish a five-dimensional offset vector chip OVT gather corresponding to a current pressed noise channel with the current pressed noise channel as a center;
in this embodiment, the five-dimensional OVT gather includes dimensions including offset, azimuth, co-depth point CDP, CDP line and time;
The establishing module 80 is specifically configured to:
selecting a first set number of CDP lines and a second set number of CDP lines by taking the current pressed noise channel as a center to form an original three-dimensional CDP channel set with the time;
setting the offset and the azimuth;
the five-dimensional OVT gather is generated based on the offset, the azimuth, and the original three-dimensional CDP gather.
Specifically, based on the offset, the azimuth, and the original three-dimensional CDP gather, the process of generating the five-dimensional OVT gather is as follows:
a. forming an OVT bin based on the offset and the azimuth;
b. constructing an initial five-dimensional OVT gather based on the OVT surface elements;
since one OVT bin corresponds to one three-dimensional CDP trace set, in this embodiment, all OVT bins may be combined to form an initial five-dimensional OVT trace set.
c. Detecting whether a current track of the original three-dimensional CDP track set belongs to the initial five-dimensional OVT track set; if yes, executing the step d, and if not, executing the step g;
specifically, whether the initial five-dimensional OVT gather includes the current track of the original three-dimensional CDP gather may be determined according to the identifier of each track in the initial five-dimensional OVT gather and the identifier of the current track of the original three-dimensional CDP gather, if the initial five-dimensional OVT gather includes the current track of the original three-dimensional CDP gather, it is indicated that the current track of the original three-dimensional CDP gather belongs to the initial five-dimensional OVT gather, otherwise, if the initial five-dimensional OVT gather does not include the current track of the original three-dimensional CDP gather, it is indicated that the current track of the original three-dimensional CDP gather does not belong to the initial five-dimensional OVT gather.
d. Detecting whether the initial five-dimensional OVT grid corresponding to the current channel is in an idle state or not; if yes, executing the step e, and if not, executing the step f;
if the current track of the original three-dimensional CDP track set belongs to the initial five-dimensional OVT track set, it is indicated that the current track of the original three-dimensional CDP track set is valid, but the corresponding content may already be added to the initial five-dimensional OVT grid, so in this embodiment, it may be detected whether the initial five-dimensional OVT grid corresponding to the current track of the original three-dimensional CDP track set is in an idle state.
e. Loading a current track of the original three-dimensional CDP track set into the initial five-dimensional OVT grid so as to update the initial five-dimensional OVT track set;
if the initial five-dimensional OVT grid corresponding to the current track of the original three-dimensional CDP gather is in an idle state, at this time, the current track of the original three-dimensional CDP gather may be loaded into the initial five-dimensional OVT grid to update the initial five-dimensional OVT gather.
f. Selecting an occupied channel or the current channel to be loaded to the initial five-dimensional OVT grid so as to update the initial five-dimensional OVT channel set;
and if the initial five-dimensional OVT grid corresponding to the current track of the original three-dimensional CDP track set is in an occupied state, selecting an occupied track or loading the current track to the initial five-dimensional OVT grid so as to update the initial five-dimensional OVT track set.
Specifically, the current trace of the original three-dimensional CDP trace set may be geometrically compared with the occupied trace in the initial five-dimensional OVT grid to determine which trace is good and which trace is selected again. For example, the amplitude may be used as a comparison target, and if the amplitude is large, this means that the track is relatively good, whereas if the amplitude is small, this means that the track is relatively bad.
g. Skipping the current track of the original three-dimensional CDP gather.
If the current channel of the original three-dimensional CDP channel set does not belong to the original five-dimensional OVT channel set, the current channel of the original three-dimensional CDP channel set can be directly skipped, next extraction is carried out, and the needed five-dimensional OVT channel set can be obtained after all channels in the original three-dimensional CDP channel set are traversed.
The analysis module 81 is configured to analyze the five-dimensional spectrum of the five-dimensional OVT gather to obtain a viewing angle of the five-dimensional OVT gather;
specifically, in the process of creating the noise-pressing model trace, it is necessary to obtain the apparent dip angles of the five-dimensional OVT trace set in four directions caused by dip angles, tendencies, bends, and the like of the stratum. This requires five-dimensional spectrum analysis to obtain the maximum correlation spectrum of each data on each sample point and the corresponding viewing inclination angles of the offset, azimuth angle, line and CDP directions. And the reliable five-dimensional noise pressing model channel for denoising can be established by superposition along the optimal visual inclination angle direction. In order to obtain the visual tilt angles, the program uses a tilt angle scanning method to manufacture a five-dimensional spectrum of a five-dimensional OVT gather, and obtains the visual tilt angles through spectrum analysis.
The superposition module 82 is configured to superimpose the five-dimensional OVT gather along the view angle to obtain a noise pressing model gather;
in this embodiment, a five-dimensional noise compression model trace is manufactured, and the five-dimensional OVT trace set is required to be superimposed in the five-dimensional OVT trace set along the optimal viewing angle direction by using a mean value superposition calculation method, so as to obtain the noise compression model trace. Because the superposition is done on a sample-by-sample basis, taking one sample from a five-dimensional OVT gather degrades into a four-dimensional volume. The offset and azimuth of the stacking process need not be further considered, and these two dimensions can be removed and projected onto a two-dimensional plane consisting of the CDP and line number to form a two-dimensional amplitude function g (x, y, t) for the current point in time t, where x and y are the coordinate positions of g in the CDP and line directions.
In actual data, there are often some singular amplitude values, which are characterized by their being much larger than normal amplitude values. Therefore, in order to aim at the embodiment of the singular amplitude value and the non-singular amplitude value, different calculation methods can be adopted to overlap the five-dimensional OVT gather along the view inclination angle, so as to obtain a five-dimensional press noise model gather.
If no singular amplitude value exists in the five-dimensional OVT channel set, the first calculation formula (2) can be selected as a mean value superposition calculation formula; if the singular amplitude values exist in the five-dimensional OVT channel set, the second calculation formula (3) can be selected as a mean superposition calculation formula:
In addition, in the stacking process in this embodiment, the stacking of one CMP gather at a certain time point is to stack data of different offset distances and azimuth angles into one line, and the five-dimensional OVT gather in the embodiment also involves different lines and tracks, which is equivalent to stacking a plurality of CMP from different lines and tracks.
And the noise pressing module 83 is configured to perform noise pressing on the initial domain seismic data corresponding to the noise pressing model channel by using conversion between different domains, so as to obtain target seismic data corresponding to the controllable seismic source.
Specifically, the initial domain seismic data corresponding to the noise pressing model channel can be converted into target domain seismic data of a target domain; performing noise suppression processing on the target domain seismic data to obtain noise-reduced seismic data; converting the noise-reduced seismic data to a time domain to obtain supplementary seismic data; and combining the initial domain seismic data with the supplementary seismic data to obtain target seismic data corresponding to the controllable seismic source.
The noise suppression device of the controllable vibration source of the embodiment establishes a five-dimensional OVT gather corresponding to a current pressed noise channel by taking the current pressed noise channel as a center; analyzing the five-dimensional spectrum of the five-dimensional OVT gather to obtain a visual inclination angle of the five-dimensional OVT gather; and superposing the five-dimensional OVT trace sets along the visual inclination angle to obtain a noise pressing model trace, performing noise pressing on the initial domain seismic data corresponding to the noise pressing model trace by utilizing conversion among different domains, and recovering effective data to obtain target seismic data corresponding to a controllable seismic source, thereby realizing effective pressing on noise of a black triangle area, better recovering effective signals, enhancing effective signals outside the black triangle area, and improving the signal-to-noise ratio of data.
It should be noted that, the apparatus of the foregoing embodiment is configured to implement the corresponding method in the foregoing embodiment, and has the beneficial effects of the corresponding method embodiment, which is not described herein again.
Example IV
In order to solve the technical problems in the prior art, the embodiment of the invention also provides noise suppression equipment for the controllable vibration source.
The noise suppression device of the controllable vibration source of the embodiment comprises a memory and a controller;
the memory has stored thereon a computer program which, when executed by a controller, performs the steps of:
establishing a five-dimensional offset vector chip OVT gather corresponding to a current pressed noise channel by taking the current pressed noise channel as a center;
analyzing the five-dimensional spectrum of the five-dimensional OVT gather to obtain a visual inclination angle of the five-dimensional OVT gather;
superposing the five-dimensional OVT gathers along the visual inclination angle to obtain a noise pressing model gather;
and performing noise suppression on the initial domain seismic data corresponding to the noise compression model channels by utilizing conversion among different domains to obtain target seismic data corresponding to the controllable seismic source.
Further, in the above embodiment, the execution of the computer program when executed by the controller may further execute the following steps:
The establishing a five-dimensional OVT gather corresponding to the current pressed noise channel by taking the current pressed noise channel as a center comprises the following steps:
selecting a first set number of CDP lines and a second set number of CDP lines by taking the current pressed noise channel as a center to form an original three-dimensional CDP channel set with the time;
setting the offset and the azimuth;
the five-dimensional OVT gather is generated based on the offset, the azimuth, and the original three-dimensional CDP gather.
Further, in the above embodiment, the execution of the computer program when executed by the controller may further execute the following steps:
forming an OVT bin based on the offset and the azimuth;
constructing an initial five-dimensional OVT gather based on the OVT surface elements;
performing the following operations on each lane in the original three-dimensional CDP lane set to produce the five-dimensional OVT lane set;
detecting whether a current track of the original three-dimensional CDP track set belongs to the initial five-dimensional OVT track set;
if the current channel of the original three-dimensional CDP channel set belongs to the initial five-dimensional OVT channel set, detecting whether the initial five-dimensional OVT grid corresponding to the current channel is in an idle state or not;
if the initial five-dimensional OVT grid corresponding to the current channel is in an idle state, loading the current channel into the initial five-dimensional OVT grid so as to update the initial five-dimensional OVT channel set;
And if the initial five-dimensional OVT grid corresponding to the current track is in an occupied state, selecting the occupied track or loading the current track to the initial five-dimensional OVT grid so as to update the initial five-dimensional OVT track set.
Further, the computer program, when executed by the controller, may perform the following steps:
and manufacturing the five-dimensional spectrum by a dip angle scanning method.
Further, in the above embodiment, the execution of the computer program when executed by the controller may further execute the following steps:
and superposing the five-dimensional OVT gathers along the visual inclination angle by using a mean superposition calculation formula to obtain a noise pressing model gather.
Further, the computer program, when executed by the controller, may perform the following steps:
if the five-dimensional OVT channel set does not have singular amplitude values, selecting a first calculation formula as the mean value superposition calculation formula;
if the five-dimensional OVT channel set has singular amplitude values, selecting a second calculation formula as the mean value superposition calculation formula;
the first calculation formula is
Wherein, N is the total number of g (x, y, t), g (x, y, t) is a two-dimensional amplitude function, x and y are the coordinate positions of g in CDP and line direction;
The second calculation formula is
Further, in the above embodiment, the execution of the computer program when executed by the controller may further execute the following steps:
converting the initial domain seismic data corresponding to the noise pressing model channel into target domain seismic data of a target domain;
performing noise suppression processing on the target domain seismic data to obtain noise-reduced seismic data;
converting the noise-reduced seismic data to a time domain to obtain supplementary seismic data;
and combining the initial domain seismic data with the supplementary seismic data to obtain target seismic data corresponding to the controllable seismic source.
Example five
In order to solve the technical problems in the prior art, an embodiment of the present invention provides a storage medium.
The storage medium provided by the embodiment of the invention stores a computer program, and the computer program realizes the following steps when being executed by a processor:
establishing a five-dimensional offset vector chip OVT gather corresponding to a current pressed noise channel by taking the current pressed noise channel as a center;
analyzing the five-dimensional spectrum of the five-dimensional OVT gather to obtain a visual inclination angle of the five-dimensional OVT gather;
superposing the five-dimensional OVT gathers along the visual inclination angle to obtain a noise pressing model gather;
And performing noise suppression on the initial domain seismic data corresponding to the noise compression model channels by utilizing conversion among different domains to obtain target seismic data corresponding to the controllable seismic source.
Further, in the above embodiment, the computer program may further execute the following steps when executed by the processor:
the establishing a five-dimensional OVT gather corresponding to the current pressed noise channel by taking the current pressed noise channel as a center comprises the following steps:
selecting a first set number of CDP lines and a second set number of CDP lines by taking the current pressed noise channel as a center to form an original three-dimensional CDP channel set with the time;
setting the offset and the azimuth;
the five-dimensional OVT gather is generated based on the offset, the azimuth, and the original three-dimensional CDP gather.
Further, in the above embodiment, the computer program may further execute the following steps when executed by the processor:
forming an OVT bin based on the offset and the azimuth;
constructing an initial five-dimensional OVT gather based on the OVT surface elements;
performing the following operations on each lane in the original three-dimensional CDP lane set to produce the five-dimensional OVT lane set;
detecting whether a current track of the original three-dimensional CDP track set belongs to the initial five-dimensional OVT track set;
If the current channel of the original three-dimensional CDP channel set belongs to the initial five-dimensional OVT channel set, detecting whether the initial five-dimensional OVT grid corresponding to the current channel is in an idle state or not;
if the initial five-dimensional OVT grid corresponding to the current channel is in an idle state, loading the current channel into the initial five-dimensional OVT grid so as to update the initial five-dimensional OVT channel set;
and if the initial five-dimensional OVT grid corresponding to the current track is in an occupied state, selecting the occupied track or loading the current track to the initial five-dimensional OVT grid so as to update the initial five-dimensional OVT track set.
Further, in the above embodiment, the computer program may further execute the following steps when executed by the processor:
and manufacturing the five-dimensional spectrum by a dip angle scanning method.
Further, in the above embodiment, the computer program may further execute the following steps when executed by the processor:
and superposing the five-dimensional OVT gathers along the visual inclination angle by using a mean superposition calculation formula to obtain a noise pressing model gather.
Further, in the above embodiment, the computer program may further execute the following steps when executed by the processor:
if the five-dimensional OVT channel set does not have singular amplitude values, selecting a first calculation formula as the mean value superposition calculation formula;
If the five-dimensional OVT channel set has singular amplitude values, selecting a second calculation formula as the mean value superposition calculation formula;
the first calculation formula is
Wherein, N is the total number of g (x, y, t), g (x, y, t) is a two-dimensional amplitude function, x and y are the coordinate positions of g in CDP and line direction;
the second calculation formula is
Further, in the above embodiment, the computer program may further execute the following steps when executed by the processor:
converting the initial domain seismic data corresponding to the noise pressing model channel into target domain seismic data of a target domain;
performing noise suppression processing on the target domain seismic data to obtain noise-reduced seismic data;
converting the noise-reduced seismic data to a time domain to obtain supplementary seismic data;
and combining the initial domain seismic data with the supplementary seismic data to obtain target seismic data corresponding to the controllable seismic source.
It is to be understood that the same or similar parts in the above embodiments may be referred to each other, and that in some embodiments, the same or similar parts in other embodiments may be referred to.
It should be noted that in the description of the present invention, the terms "first," "second," and the like are used for descriptive purposes only and are not to be construed as indicating or implying relative importance. Furthermore, in the description of the present invention, unless otherwise indicated, the meaning of "plurality" means at least two.
Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing specific logical functions or steps of the process, and further implementations are included within the scope of the preferred embodiment of the present invention in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the present invention.
It is to be understood that portions of the present invention may be implemented in hardware, software, firmware, or a combination thereof. In the above-described embodiments, the various steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, may be implemented using any one or combination of the following techniques, as is well known in the art: discrete logic circuits having logic gates for implementing logic functions on data signals, application specific integrated circuits having suitable combinational logic gates, programmable Gate Arrays (PGAs), field Programmable Gate Arrays (FPGAs), and the like.
Those of ordinary skill in the art will appreciate that all or a portion of the steps carried out in the method of the above-described embodiments may be implemented by a program to instruct related hardware, where the program may be stored in a computer readable storage medium, and where the program, when executed, includes one or a combination of the steps of the method embodiments.
Furthermore, the functional units in the various embodiments of the present invention may be integrated into one processing module 32, or the units may exist alone physically, or two or more units may be integrated into one module. The integrated modules may be implemented in hardware or in software functional modules. The integrated modules may also be stored in a computer readable storage medium if implemented in the form of software functional modules and sold or used as a stand-alone product.
The above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, or the like.
In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiments or examples. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
Although the embodiments of the present invention are disclosed above, the embodiments are only used for the convenience of understanding the present invention, and are not intended to limit the present invention. Any person skilled in the art can make any modification and variation in form and detail without departing from the spirit and scope of the present disclosure, but the scope of the present disclosure is still subject to the scope of the present disclosure as defined by the appended claims.
Claims (10)
1. A method of noise suppression for a vibroseis, comprising:
establishing a five-dimensional offset vector chip OVT gather corresponding to a current pressed noise channel by taking the current pressed noise channel as a center;
analyzing the five-dimensional spectrum of the five-dimensional OVT gather to obtain a visual inclination angle of the five-dimensional OVT gather;
superposing the five-dimensional OVT gathers along the visual inclination angle to obtain a noise pressing model gather;
performing noise suppression on the initial domain seismic data corresponding to the noise compression model channels by utilizing conversion among different domains to obtain target seismic data corresponding to a controllable seismic source;
the method for establishing the five-dimensional offset vector sheet OVT gather corresponding to the current pressed noise channel by taking the current pressed noise channel as the center comprises the following steps:
Selecting a first set number of CDP lines and a second set number of CDP lines by taking the current pressed noise channel as a center to form an original three-dimensional CDP channel set with time;
setting an offset distance and an azimuth angle;
the five-dimensional OVT gather is generated based on the offset, the azimuth, and the original three-dimensional CDP gather.
2. The method of claim 1, wherein the five-dimensional OVT gather includes dimensions including offset, azimuth, common depth point CDP, CDP line and time;
the establishing a five-dimensional OVT gather corresponding to the current pressed noise channel by taking the current pressed noise channel as a center comprises the following steps:
selecting a first set number of CDP lines and a second set number of CDP lines by taking the current pressed noise channel as a center to form an original three-dimensional CDP channel set with the time;
setting the offset and the azimuth;
the five-dimensional OVT gather is generated based on the offset, the azimuth, and the original three-dimensional CDP gather.
3. The method of noise suppression of a vibroseis of claim 2, wherein generating the five-dimensional OVT gather based on the offset, the azimuth, and the original three-dimensional CDP gather comprises:
Forming an OVT bin based on the offset and the azimuth;
constructing an initial five-dimensional OVT gather based on the OVT surface elements;
performing the following operations on each lane in the original three-dimensional CDP lane set to produce the five-dimensional OVT lane set;
detecting whether a current track of the original three-dimensional CDP track set belongs to the initial five-dimensional OVT track set;
if the current channel of the original three-dimensional CDP channel set belongs to the initial five-dimensional OVT channel set, detecting whether the initial five-dimensional OVT grid corresponding to the current channel is in an idle state or not;
if the initial five-dimensional OVT grid corresponding to the current channel is in an idle state, loading the current channel into the initial five-dimensional OVT grid so as to update the initial five-dimensional OVT channel set;
and if the initial five-dimensional OVT grid corresponding to the current track is in an occupied state, selecting the occupied track or loading the current track to the initial five-dimensional OVT grid so as to update the initial five-dimensional OVT track set.
4. The method for suppressing noise of a vibroseis according to claim 2, wherein before the analyzing the five-dimensional spectrum of the five-dimensional OVT gather to obtain the dip angle of the five-dimensional OVT gather, further comprises:
and manufacturing the five-dimensional spectrum by a dip angle scanning method.
5. The method for suppressing noise of a controllable source according to claim 1, wherein the superimposing the five-dimensional OVT gather along the view angle to obtain a noise suppressing model gather includes:
and superposing the five-dimensional OVT gathers along the visual inclination angle by using a mean superposition calculation formula to obtain a noise pressing model gather.
6. The method of noise suppression of a vibroseis of claim 5, further comprising:
if the five-dimensional OVT channel set does not have singular amplitude values, selecting a first calculation formula as the mean value superposition calculation formula;
if the five-dimensional OVT channel set has singular amplitude values, selecting a second calculation formula as the mean value superposition calculation formula;
the first calculation formula is
Wherein, N is the total number of g (x, y, t), g (x, y, t) is a two-dimensional amplitude function, x and y are the coordinate positions of g in CDP and line direction;
the second calculation formula is
Wherein g mute (x, y, t) is the amplitude value of g (x, y, t) after the singular value is removed, d is the distance from the superimposed track to the current track, lpow is the distance weighting factor, p is the amplitude normalization factor, and M is the number of coverage times.
7. The method for suppressing noise of a vibroseis according to claim 1, wherein performing noise suppression on the initial domain seismic data corresponding to the noise suppression model channel by using conversion between different domains to obtain target seismic data corresponding to the vibroseis, comprises:
Converting the initial domain seismic data corresponding to the noise pressing model channel into target domain seismic data of a target domain;
performing noise suppression processing on the target domain seismic data to obtain noise-reduced seismic data;
converting the noise-reduced seismic data to a time domain to obtain supplementary seismic data;
and combining the initial domain seismic data with the supplementary seismic data to obtain target seismic data corresponding to the controllable seismic source.
8. A vibroseis noise suppression apparatus, comprising:
the building module is used for building a five-dimensional offset vector chip OVT gather corresponding to the current pressed noise channel by taking the current pressed noise channel as a center;
the analysis module is used for analyzing the five-dimensional spectrum of the five-dimensional OVT gather to obtain the visual inclination angle of the five-dimensional OVT gather;
the superposition module is used for superposing the five-dimensional OVT gathers along the visual inclination angle to obtain a noise pressing model gather;
the noise pressing module is used for performing noise pressing on the initial domain seismic data corresponding to the noise pressing model channels by utilizing conversion among different domains to obtain target seismic data corresponding to the controllable seismic source;
the establishing module is configured to establish a five-dimensional offset vector chip OVT gather corresponding to a current pressed noise channel with the current pressed noise channel as a center, and includes:
Selecting a first set number of CDP lines and a second set number of CDP lines by taking the current pressed noise channel as a center to form an original three-dimensional CDP channel set with time;
setting an offset distance and an azimuth angle;
the five-dimensional OVT gather is generated based on the offset, the azimuth, and the original three-dimensional CDP gather.
9. A vibroseis noise suppressing apparatus comprising a memory and a controller;
stored on the memory is a computer program which, when executed by a controller, implements the steps of the method according to any one of claims 1 to 7.
10. A storage medium having stored thereon a computer program, which when executed by a controller, implements the steps of the method according to any of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010834408.4A CN114076985B (en) | 2020-08-18 | 2020-08-18 | Method, device, equipment and storage medium for suppressing noise of controllable vibration source |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010834408.4A CN114076985B (en) | 2020-08-18 | 2020-08-18 | Method, device, equipment and storage medium for suppressing noise of controllable vibration source |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114076985A CN114076985A (en) | 2022-02-22 |
| CN114076985B true CN114076985B (en) | 2024-04-09 |
Family
ID=80281390
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010834408.4A Active CN114076985B (en) | 2020-08-18 | 2020-08-18 | Method, device, equipment and storage medium for suppressing noise of controllable vibration source |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114076985B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115184994B (en) * | 2022-06-08 | 2023-04-18 | 北京东方联创地球物理技术有限公司 | Multi-azimuth acquisition three-dimensional seismic data fusion processing method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000338256A (en) * | 1999-05-31 | 2000-12-08 | Railway Technical Res Inst | Surface wave phase speed detection system, and method for detecting phase speed of surface wave |
| CN105277985A (en) * | 2014-07-24 | 2016-01-27 | 中国石油化工股份有限公司 | OVT-domain seismic data regularization method based on image processing |
| CN110618457A (en) * | 2018-06-20 | 2019-12-27 | 中国石油化工股份有限公司 | Combined denoising method and device for vibroseis seismic data |
| CN110737022A (en) * | 2018-07-20 | 2020-01-31 | 中国石油化工股份有限公司 | Suppression method for vibroseis to excite noise of seismic data black triangle area |
| CN110837121A (en) * | 2018-08-16 | 2020-02-25 | 中国石油化工股份有限公司 | Method and system for eliminating random noise in seismic data |
| CN111352158A (en) * | 2018-12-24 | 2020-06-30 | 中国石油天然气股份有限公司 | Seismic signal enhancement method and device |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9435903B2 (en) * | 2013-10-11 | 2016-09-06 | Chevron U.S.A. Inc. | System and method for regularizing seismic data |
| US10191165B2 (en) * | 2015-01-13 | 2019-01-29 | Cgg Services Sas | Using an offset vector tile gather to image a subsurface |
-
2020
- 2020-08-18 CN CN202010834408.4A patent/CN114076985B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000338256A (en) * | 1999-05-31 | 2000-12-08 | Railway Technical Res Inst | Surface wave phase speed detection system, and method for detecting phase speed of surface wave |
| CN105277985A (en) * | 2014-07-24 | 2016-01-27 | 中国石油化工股份有限公司 | OVT-domain seismic data regularization method based on image processing |
| CN110618457A (en) * | 2018-06-20 | 2019-12-27 | 中国石油化工股份有限公司 | Combined denoising method and device for vibroseis seismic data |
| CN110737022A (en) * | 2018-07-20 | 2020-01-31 | 中国石油化工股份有限公司 | Suppression method for vibroseis to excite noise of seismic data black triangle area |
| CN110837121A (en) * | 2018-08-16 | 2020-02-25 | 中国石油化工股份有限公司 | Method and system for eliminating random noise in seismic data |
| CN111352158A (en) * | 2018-12-24 | 2020-06-30 | 中国石油天然气股份有限公司 | Seismic signal enhancement method and device |
Non-Patent Citations (2)
| Title |
|---|
| B.Delouis ; J.Charlety ; M.Vallée ; 许若冰 ; .用强震记录的近场谱快速测定中强地震矩震级M_w的方法(MWSYNTH).世界地震译丛.2009,(03),全文. * |
| 李飞 ; 段文胜 ; 赵锐锐 ; 罗莉莉 ; 党青宁 ; .应用OVT域体τ-p变换提高地震资料信噪比.石油地球物理勘探.2015,(03),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114076985A (en) | 2022-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Birtill et al. | The application of phased arrays to the analysis of seismic body waves | |
| US9766357B2 (en) | Seismic image dip decomposition estimation and recomposition | |
| CN104237940B (en) | Diffracted wave imaging method and device based on dynamic characteristics | |
| US8209126B2 (en) | Wavefront-defined Radon transform | |
| Rentsch et al. | Fast location of seismicity: A migration-type approach with application to hydraulic-fracturing data | |
| CN104570125A (en) | Method for improving imaging velocity model precision by utilizing well data | |
| CN104280777B (en) | Method for suppressing interference of seismic data multiples on land | |
| CN105277985A (en) | OVT-domain seismic data regularization method based on image processing | |
| CN102073064B (en) | Method for improving velocity spectrum resolution by using phase information | |
| Nowack et al. | Application of Gaussian-beam migration to multiscale imaging of the lithosphere beneath the Hi-CLIMB array in Tibet | |
| Amoroso et al. | S‐wave identification by polarization filtering and waveform coherence analyses | |
| CN101923175B (en) | Method for directly generating angle gathers by using wave-equation migration | |
| CN105301648A (en) | Method of acquiring common reflection surface stacking parameters | |
| CN114076985B (en) | Method, device, equipment and storage medium for suppressing noise of controllable vibration source | |
| CN102736108B (en) | True three-dimensional earthquake data noise suppressing method based on spline fitting | |
| CN104977615B (en) | A kind of multiple ripple drawing method of deep water OBC data based on modeling statistics pickup | |
| CN105510975A (en) | Method and device for improving signal-to-noise ratio of seismic data | |
| US20060253257A1 (en) | Method of processing seismic data for AVO or AVOA characterization | |
| Hu et al. | Wavefield reconstruction of teleseismic receiver function with the stretching‐and‐squeezing interpolation method | |
| AU2006237350B2 (en) | Seismic data processing method for RMO picking | |
| CN102830424A (en) | Method for calculating combination parameters of detector | |
| CN102338884A (en) | Elliptic window direction band-pass amplitude-preserving filtering data processing method in geophysical prospecting | |
| CN112630822B (en) | Method and system for processing OVT domain pre-stack seismic data | |
| Bakhtiari Rad et al. | Using common-reflection-surface stack for enhanced near-surface seismic reflection imaging: Examples from consolidated and unconsolidated environments | |
| CN109254322B (en) | Method for establishing residual static correction model channel based on heredity |
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 |