CN108732617B - Seismic data waste channel rejection method and device - Google Patents
Seismic data waste channel rejection method and device Download PDFInfo
- Publication number
- CN108732617B CN108732617B CN201710265486.5A CN201710265486A CN108732617B CN 108732617 B CN108732617 B CN 108732617B CN 201710265486 A CN201710265486 A CN 201710265486A CN 108732617 B CN108732617 B CN 108732617B
- Authority
- CN
- China
- Prior art keywords
- seismic
- channel
- peak
- trace
- trough
- 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
- 239000002699 waste material Substances 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000004590 computer program Methods 0.000 claims description 2
- 239000003245 coal Substances 0.000 abstract description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 abstract description 4
- 239000011707 mineral Substances 0.000 abstract description 4
- 239000003079 shale oil Substances 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
A method and apparatus for removing the waste channels from seismic channel are disclosed. The method comprises the following steps: 1) identifying the peak and valley points of the single-channel data; 2) counting the peak and trough number value K of each seismic channel in the seismic data and the average value Ave of the peak and trough number of the seismic channels; 3) comparing the peak-trough number value K of each trace with a set threshold F, and determining the seismic trace as a waste trace if K > is F; 4) and removing the determined waste channel. According to the method and the device for eliminating the seismic channel data waste channel, the waste channel can be effectively positioned, and exploration and development services of mineral products, conventional oil gas, shale oil gas and coal bed gas are provided.
Description
Technical Field
The invention relates to the field of seismic exploration and development of mineral products, oil gas, shale oil gas and coal bed gas, in particular to a seismic data waste channel removing method and device.
Background
In the seismic data processing, the data of each shot gather is subjected to waste channel elimination firstly, so that the later seismic data processing is more accurate and reliable.
The NEWS earthquake comprehensive interpretation system (hereinafter referred to as NEWS system) is earthquake comprehensive interpretation software independently developed and developed by Chinese petrifaction. The prestack inversion in the NEWS system uses common shot gather data or some other prestack data, which often has a part of waste traces, so we need to develop a reliable and effective method for removing waste traces to solve the above problems.
The information disclosed in this background section is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person skilled in the art.
Disclosure of Invention
The invention aims to solve the problems in the prior art and provides a method and a device for eliminating the seismic data waste channel based on extreme point analogy.
According to one aspect of the invention, a seismic channel data waste channel rejection method is provided. The method may comprise the steps of:
1) identifying the peak and valley points of the single-channel data;
2) counting the peak and trough number value K of each seismic channel in the seismic data and the average value Ave of the peak and trough number of the seismic channels;
3) comparing the peak-trough number value K of each trace with a set threshold F, and determining the seismic trace as a waste trace if K > is F;
4) and removing the determined waste channel.
Preferably, the peak identification formula is:
wherein, XnFor any number of samples, X, in a single seismic trace sequencen、Xn+1、Xn+2The number of samples is three in succession.
Preferably, the trough identification formula is:
wherein, XnFor any number of samples, X, in a single seismic trace sequencen、Xn+1、Xn+2The number of samples is three in succession.
Preferably, the average Ave of the number of peaks and troughs of the seismic trace in the step 2) is calculated according to the following formula:
k is the peak-trough number value of each channel, and n is the channel number of the seismic channel in the seismic data volume.
Preferably, the threshold F is set between 15% and 25% of the average value Ave of the number of peaks and troughs of the seismic trace.
According to another aspect of the invention, a seismic channel waste channel rejection device is provided. The apparatus comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program performs the steps of:
1) identifying the peak and valley points of the single-channel data;
2) counting the peak and trough number value K of each seismic channel in the seismic data and the average value Ave of the peak and trough number of the seismic channels;
3) comparing the peak-trough number value K of each trace with a set threshold F, and determining the seismic trace as a waste trace if K > is F;
4) and removing the determined waste channel.
Preferably, the peak identification formula is:
wherein, XnFor any number of samples, X, in a single seismic trace sequencen、Xn+1、Xn+2The number of samples is three in succession.
Preferably, the trough identification formula is:
wherein, XnFor any number of samples, X, in a single seismic trace sequencen、Xn+1、Xn+2The number of samples is three in succession.
Preferably, the average Ave of the number of peaks and troughs of the seismic trace in the step 2) is calculated according to the following formula:
k is the peak-trough number value of each channel, and n is the channel number of the seismic channel in the seismic data volume.
Preferably, the threshold F is set between 15% and 25% of the average value Ave of the number of peaks and troughs of the seismic trace.
According to the method and the device for eliminating the seismic channel data waste channel, the waste channel can be effectively positioned, and exploration and development services of mineral products, conventional oil gas, shale oil gas and coal bed gas are provided.
The method and apparatus of the present invention have other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the invention.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent by describing in more detail exemplary embodiments thereof with reference to the attached drawings, in which like reference numerals generally represent like parts.
FIG. 1 is a flow chart of a seismic trace data culling method according to an exemplary embodiment of the invention.
Detailed Description
The invention will be described in more detail below with reference to the accompanying drawings. While the preferred embodiments of the present invention are shown in the drawings, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
The seismic trace data culling method according to an exemplary embodiment of the invention is described in detail below with reference to fig. 1. The method mainly comprises the following steps:
step 1: and identifying the peak and valley points of the single-channel data.
A single seismic trace is a sequence of floating-point numbers, any point in the sequence being described as Xn. The wave crest is characterized in that the middle number of the three continuous sample points is the maximum value of the three samples, and the wave trough is characterized in that the three continuous sample points areThe median number is the minimum of the three. The values of the three numbers are the same and are not judged.
The peak identification formula is:
wherein, XnFor any number of samples, X, in a single seismic trace sequencen、Xn+1、Xn+2The number of samples is three in succession.
The trough identification formula is:
wherein, XnFor any number of samples, X, in a single seismic trace sequencen、Xn+1、Xn+2The number of samples is three in succession.
Identifying wave crests by using a formula (1), identifying wave troughs by using a formula (2), and recording the quantity K and the positions of the wave crests and the wave troughs of each channel. And determining the specific position of the data by using the line number, the shot number and the track number in the seismic data.
Step 2: and counting the peak and trough quantity value K of each seismic channel in the seismic data and the average value Ave of the peak and trough quantity of the seismic channels.
Calculating the average value Ave of the number of peaks and troughs of the seismic trace according to the following formula:
k is the peak-trough number value of each channel, and n is the channel number of the seismic channel in the seismic data volume.
And step 3: and comparing the peak-trough number value K of each trace with a set threshold value F, and determining the seismic trace as a waste trace if K > -F.
The wave crest and valley point of the waste channel in the seismic data is obviously more than that of the normal channel, so that the waste channel can be positioned by setting a threshold value F.
The threshold value F can be between 15% and 25% of the average value Ave of the number of wave crests and wave troughs of the seismic channel in the data body, then the value K of the number of wave crests and wave troughs of each channel is compared with F, if K is larger than or equal to F, the channel is a waste channel; if K is less than F, then this track is a normal track.
For example, F is set to 20% of the average Ave of the number of peaks and troughs of the seismic trace, i.e., F ═ Ave × 20%.
If K > is F, the channel is a waste channel; if K < F, this track is a normal track.
Since we have recorded the position of each track in the first step, the waste track can be successfully located according to the calculation result.
And 4, step 4: and removing the determined waste channel.
According to the seismic channel waste channel removing method and device, waste channels can be effectively positioned, and exploration and development services for mineral products, conventional oil gas, shale oil gas and coal bed gas are provided.
It will be appreciated by persons skilled in the art that the above description of embodiments of the invention is intended only to illustrate the benefits of embodiments of the invention and is not intended to limit embodiments of the invention to any examples given.
Having described embodiments of the present invention, the foregoing description is intended to be exemplary, not exhaustive, and not limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments. The terminology used herein is chosen in order to best explain the principles of the embodiments, the practical application, or improvements made to the technology in the marketplace, or to enable others of ordinary skill in the art to understand the embodiments disclosed herein.
Claims (8)
1. A seismic data waste channel rejection method is characterized by comprising the following steps:
1) identifying the peak and valley points of the single-channel data;
2) counting the peak and trough number value K of each seismic channel in the seismic data and the average value Ave of the peak and trough number of the seismic channels;
3) comparing the peak-trough number value K of each trace with a set threshold F, and determining the seismic trace as a waste trace if K > is F;
4) removing the determined waste channel;
wherein the set threshold F is between 15% and 25% of the average value Ave of the number of the wave crests and the wave troughs of the seismic channel.
2. The seismic data culling method of claim 1, wherein a peak identification formula is:
wherein, XnFor any number of samples, X, in a single seismic trace sequencen、Xn+1、Xn+2The number of samples is three in succession.
3. The seismic data culling method of claim 1, wherein a trough identification formula is:
wherein, XnFor any number of samples, X, in a single seismic trace sequencen、Xn+1、Xn+2The number of samples is three in succession.
4. The seismic data culling method of claim 1, wherein the average seismic trace peak-trough number Ave in step 2) is calculated according to the following formula:
k is the peak-trough number value of each channel, and n is the channel number of the seismic channel in the seismic data volume.
5. A seismic data culling apparatus, the apparatus comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program performs the steps of:
1) identifying the peak and valley points of the single-channel data;
2) counting the peak and trough number value K of each seismic channel in the seismic data and the average value Ave of the peak and trough number of the seismic channels;
3) comparing the peak-trough number value K of each trace with a threshold value F, and determining the seismic trace as a waste trace if K > - < F >;
4) removing the determined waste channel;
wherein the set threshold F is between 15% and 25% of the average value Ave of the number of the wave crests and the wave troughs of the seismic channel.
6. The seismic data culling apparatus of claim 5, wherein a peak identification formula is:
wherein, XnFor any number of samples, X, in a single seismic trace sequencen、Xn+1、Xn+2The number of samples is three in succession.
7. The seismic data culling apparatus of claim 5, wherein the trough identification formula is:
wherein, XnFor any number of samples, X, in a single seismic trace sequencen、Xn+1、Xn+2Is made continuousThree sample numbers.
8. The seismic data waste channel rejection device according to claim 5, wherein the average value Ave of the number of peaks and troughs of the seismic channel in step 2) is calculated according to the following formula:
k is the peak-trough number value of each channel, and n is the channel number of the seismic channel in the seismic data volume.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710265486.5A CN108732617B (en) | 2017-04-21 | 2017-04-21 | Seismic data waste channel rejection method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710265486.5A CN108732617B (en) | 2017-04-21 | 2017-04-21 | Seismic data waste channel rejection method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108732617A CN108732617A (en) | 2018-11-02 |
| CN108732617B true CN108732617B (en) | 2020-04-07 |
Family
ID=63933882
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710265486.5A Active CN108732617B (en) | 2017-04-21 | 2017-04-21 | Seismic data waste channel rejection method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108732617B (en) |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7149630B2 (en) * | 2005-01-13 | 2006-12-12 | Bp Corporation North America Inc. | Method of DMO calculation for use in seismic exploration |
| WO2012146894A2 (en) * | 2011-04-26 | 2012-11-01 | Kenneth Rayvenor Lusty Armitage | Oil & gas exploration and production |
| CN104316963A (en) * | 2014-10-21 | 2015-01-28 | 中国石油天然气股份有限公司 | Abnormal track identification method |
| CN104375188B (en) * | 2014-11-19 | 2017-02-08 | 中国石油天然气集团公司 | Seismic wave transmission attenuation compensation method and device |
| CN104459779B (en) * | 2014-12-09 | 2017-07-28 | 中国石油天然气集团公司 | A kind of method and apparatus of the automatic identification in Unusual Seismic road |
-
2017
- 2017-04-21 CN CN201710265486.5A patent/CN108732617B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108732617A (en) | 2018-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kennel et al. | Method to distinguish possible chaos from colored noise and to determine embedding parameters | |
| EP2405332A1 (en) | Method for determining a touch event and touch sensitive device | |
| EP3106970A1 (en) | Signal adaptive filtering for touch detection | |
| EP1881444A1 (en) | Pattern collation device and pattern collating method thereof, and pattern collating program | |
| ATE407409T1 (en) | AUTOMATIC GESTURE DETECTION | |
| CA2587947A1 (en) | Method for processing at least two sets of seismic data | |
| EP1887520A1 (en) | Image processing device and object region tracking program | |
| CN110969092B (en) | Pulse signal denoising method and device and terminal equipment | |
| CN108345586B (en) | Text duplicate removal method and system | |
| US20160077641A1 (en) | Adaptive Thresholding for Touch Screen Input | |
| CN105094708A (en) | Method and apparatus for predicting disk capacity | |
| US9119546B2 (en) | R-PEAK detection apparatus and control method thereof | |
| US5023918A (en) | Cross-product filter | |
| KR101384317B1 (en) | Apparatus and method for blocking the objectionable multimedia based on multimodal and multiscale features | |
| CN108107086A (en) | A kind of gas detection method and gas sensor based on array gas sensor | |
| CN108732617B (en) | Seismic data waste channel rejection method and device | |
| CN109557583B (en) | Seismic attribute extraction method and system | |
| US9739873B2 (en) | Range sidelobe suppression | |
| CN110796130A (en) | Method, device and computer storage medium for character recognition | |
| Matías et al. | Support vector machines and gradient boosting for graphical estimation of a slate deposit | |
| CN104781698A (en) | System and method for isolating signal in seismic data | |
| CN111222051B (en) | Training method and device for trend prediction model | |
| CN109977737A (en) | A kind of character recognition Robust Method based on Recognition with Recurrent Neural Network | |
| CN111951791B (en) | Voiceprint recognition model training method, electronic equipment and storage medium | |
| Kim et al. | Non-contact gesture recognition using the electric field disturbance for smart device application |
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 |