CN107402404A - The wave field separation method of VSP data - Google Patents
The wave field separation method of VSP data Download PDFInfo
- Publication number
- CN107402404A CN107402404A CN201610331007.0A CN201610331007A CN107402404A CN 107402404 A CN107402404 A CN 107402404A CN 201610331007 A CN201610331007 A CN 201610331007A CN 107402404 A CN107402404 A CN 107402404A
- Authority
- CN
- China
- Prior art keywords
- wave field
- lineups
- point
- vsp data
- field separation
- 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.)
- Granted
Links
- 238000000926 separation method Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 claims abstract description 31
- 238000001914 filtration Methods 0.000 claims description 25
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000009416 shuttering 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
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)
- Image Processing (AREA)
Abstract
A kind of wave field separation method of VSP data, comprises the following steps：Step 1：Enter traveling-wave field pickup to pending VSP data, obtain the information of multigroup wave field, record the position of the lineups of every group of wave field；Step 2：It is the dot-dash timing window according to the position of the point for the every bit on lineups for each lineups；Step 3：For each lineups, for the every bit on the lineups, according to the when window delimited for the point, one-dimensional median filter process is carried out to the point.This method eliminate existing method even up and reverse drawing process, and multigroup wave field information can be handled simultaneously, improve the efficiency of wave field separation.
Description
Technical field
The present invention relates to wave field separation technology, more particularly to a kind of wave field separation method of VSP data.
Background technology
Wave field separation is very important link step in the processing of VSP data, for VSP data, essence
Really efficient wave field separation is the basis of data processing, because containing many wave field information in VSP data,
And various wave fields are interweaved, it is necessary to this complex wave field progress efficiently accurate wave field separation, obtain
Relatively accurate single wave field is just advantageous to the next step processing of data and improves the utilization rate of VSP data.
For the accurate efficient of the wave field information such as upgoing wave, down going wave and the converted wave in VSP data trace gathers
Separation, is advantageous to obtain more preferable underground medium attribute information, is cutd open for the ensuing imaging offset
The precondition that face provides, also joint interpretation and inverting for well, and then improve the resolution of seismic data
Rate provides very big help.Existing medium filtering wave field separation technology needs to separate ripple to the need in data first
The lineups of field carry out evening up processing, reapply medium filtering principle and eliminate wave field, finally carry out again at reverse drawing
Reason reverts to the trace gather after separation.This method due to exist it is cumbersome even up and reverse drawing process, separative efficiency
Substantially reduce, and a kind of wave field information can only be removed every time.
The content of the invention
It is an object of the invention to provide a kind of wave field separation method of VSP data, existing method is that eliminated
Even up and reverse drawing process, and multilayer wave field information can be handled simultaneously, improve the efficiency of wave field separation.
The present invention uses solution below：
A kind of wave field separation method of VSP data, comprises the following steps：
Step 1：Enter traveling-wave field pickup to pending VSP data, obtain the information of multigroup wave field, record
The position of the lineups of every group of wave field；
Step 2：For each lineups, for the every bit on the lineups, according to the position of the point
For the dot-dash timing window；
Step 3：For each lineups, for the every bit on the lineups, delimited according to for the point
When window, one-dimensional median filter process is carried out to the point.
Preferably, one-dimensional median filter process is carried out according to below equation (1)：
G (x)=med f (x-k) ..., f (x) ..., f (x+k), k ∈ N } (1)
Wherein, g (x) is the output of one-dimensional medium filtering, and f (x) is to need to carry out the original signal of one-dimensional medium filtering, med
The intermediate value after sequence is represented, N is positive integer, the number put when representing described in window.
Preferably, N is equal to 5,7 or 9.
Preferably for each lineups, for the every bit on the lineups, with the position of the point
Centered on, the dot-dash timing window of selection same number above and below the point.
Compared with prior art, the beneficial effects of the present invention are the intermediate value improved in existing wave field separation method
Filter processing method, accurate efficient separation can be carried out simultaneously to multilayer complexity VSP wave fields.
Brief description of the drawings
Disclosure exemplary embodiment is described in more detail in conjunction with the accompanying drawings, the disclosure it is above-mentioned with
And other purposes, feature and advantage will be apparent.
Fig. 1 a to Fig. 1 c show the example of one-dimensional median filter method；
Fig. 2 shows the example of two dimension median filter method；
Fig. 3 a and Fig. 3 b respectively illustrate the circular shuttering and Crisscross template of two dimension median filter method；
Fig. 4 a and Fig. 4 b respectively illustrate the ripple that VSP data is carried out using traditional one-dimensional median filter method
Separation, the schematic diagram before being evened up to lineups and after evening up；
Fig. 5 shows the flow chart of the wave field separation method of the VSP data according to exemplary embodiment；
Fig. 6 shows the same phase picked up in the wave field separation method according to the VSP data of exemplary embodiment
Axle；
Fig. 7 shows same for one in the wave field separation method according to the VSP data of exemplary embodiment
The when window that phase axle is chosen；
Fig. 8 shows the wave field separation method using the VSP data according to exemplary embodiment, obtained ripple
Data after the separation of field.
Embodiment
Preferred embodiment of the present disclosure is more fully described below with reference to accompanying drawings.Although this is shown in accompanying drawing
Disclosed preferred embodiment, however, it is to be appreciated that may be realized in various forms the disclosure without should be by here
The embodiment of elaboration is limited.On the contrary, these embodiments are provided so that the disclosure is more thorough and complete,
And the scope of the present disclosure can be intactly communicated to those skilled in the art.
In order that those skilled in the art more fully understand the present invention, it is discussed in detail first conventional based on intermediate value
The wave field separation technology of filtering.
Medium filtering is that medium filtering can based on a kind of sequencing statistical theoretical nonlinear signal processing technology
Suppress and eliminate noise.The general principle of medium filtering is the value use any in digital picture or Serial No.
The Mesophyticum of each point value replaces in one neighborhood of the point, makes the numerical value of surrounding mutual relatively, so as to disappear
Except isolated noise spot, and then the purpose of smoothing processing.
Median filter method includes one-dimensional medium filtering and two dimension median filter.One-dimensional medium filtering usually exists
The point of same number is respectively taken before and after adjacent with pending point, then by this odd number numerical value of acquirement from it is small to
Sort successively greatly, the value of pending point is replaced with the intermediate value after sequence, is then carried out successively, as one
Tie up medium filtering.
Fig. 1 a-1c show the example of an one-dimensional medium filtering.In this example, it is each before and after numerical value 25
Two values (Fig. 1 a) are chosen, Fig. 1 b show the tables of data obtained after sequence, and now intermediate value becomes
For 13, numerical value 25 to be filtered then is replaced with 13, obtains the tables of data (Fig. 1 c) after medium filtering.
Shown in the output of one-dimensional medium filtering such as below equation (1)：
G (x)=med f (x-k) ..., f (x) ..., f (x+k), k ∈ N } (1)
Wherein, g (x) is the output of one-dimensional medium filtering, and f (x) is to need to carry out the original signal of one-dimensional medium filtering, med
The intermediate value after sequence is represented, N is positive integer, the number put when representing in window.
Two dimension median filter needs to select the two-dimentional sleiding form of certain structure, this module can be it is square, ten
Font (see Fig. 3 b) or circular (see Fig. 3 a) etc., numerical value in plate is ranked up according to the size of numerical value,
Generate the sequence of values of monotone increasing (or decline) (see Fig. 2).
Shown in the output of two dimension median filter such as below equation (2)：
G (x, y)=med { f (x-k, y-l), k, l ∈ W } (2)
Wherein, g (x, y) is the output of two dimension median filter, and f (x, y) is the original signal for needing to carry out two dimension median filter,
W is two dimension pattern plate, usually 5X 5, the regions of 3X 3, and it can be different shapes, such as wire,
Circular (see Fig. 3 a), cross (see Fig. 3 b), annular etc., med represent the intermediate value after sequence.
One-dimensional median filter method is used in the wave field separation of VSP data.Filtered using traditional intermediate value
Wave method needs the wave field lineups of pickup being pulled to that horizontal (Fig. 4 a and Fig. 4 b are shown respectively to enter lineups
Schematic diagram before row is evened up and after evening up), then the lineups after evening up are taken at medium filtering by road
Reason, eliminates the lineups chosen, and lineups reverse drawing finally is returned into reset condition again, so completes a step wave field
Separation.It can be seen that during from this, traditional median filter method is comparatively laborious, and is artificially evening up
With unnecessary mistake can be introduced in reverse drawing process.In addition, this method disappears again after can only eliminating a kind of wave field
Except another wave field, the abundant information received in wavelength complexity, data, wave field are interweaved serious
In the case of, it can be substantially reduced using this traditional median filter method separative efficiency, influence processing progress.
The defects of for being carried out using traditional median filter method present in the wave field separation of VSP data, this
Invention proposes a kind of wave field separation method of new VSP data, and it is carried out to traditional median filter method
Improve, handled for per pass seismic signal, multigroup wave field can be filtered simultaneously along lineups direction
Ripple.
Below with reference to Fig. 5 descriptions according to the wave field separation method of the VSP data of exemplary embodiment, it includes
Following steps：
Step 1：Traveling-wave field pickup is entered to the VSP data of pending wave field separation, obtains the letter of multigroup wave field
Breath, record the position of the lineups of every group of wave field.
Fig. 6, which is shown, enters the example that traveling-wave field picks up to the VSP data of pending wave field separation, black in Fig. 6
Color solid line represents the lineups of multigroup wave field of pickup, and every group of wave field includes a lineups, picked up altogether in Fig. 6
6 groups of wave fields are taken, correspondingly including 6 lineups.
Step 2：For each lineups, for the every bit on the lineups, being according to the position of the point should
Dot-dash timing window.
The wave field information picked up from VSP data is represented by discrete data point, for each lineups, pin
To the every bit on the lineups, (also referred to as filtered in some documents for the dot-dash timing window according to the position of the point
Ripple window).
When window determine the range of choice of point adjacent with pending point in above-mentioned one-dimensional medium filtering, pair when window in
The numerical value of point be ranked up, take the intermediate value after sequence to replace the numerical value of pending point.
Preferably, in this step, for any point on lineups, using the position of the point as
Center, the dot-dash timing window of the same number of square and lower section selection at that point, for one-dimensional in next step
Medium filtering.
Preferably, centered on the position of the point, side and lower section choose 5 altogether including the point at that point
Point, 7 points, 9: 11 points or 13 dot-dash timing windows.
Fig. 7 is shown for the lineups (lineups that arrow indicates in Fig. 7) shown in Fig. 6, same along this
The window when direction of phase axle is each point selection on lineups, the when window collection ultimately formed, such as arrow in Fig. 7
Shown in indicated black parallelogram wire frame (the when window of a lineups, other same phases are only shown herein
Axle is similar with its).
Step 3：For each lineups, for the every bit on the lineups, according to for the point delimit when
Window, one-dimensional median filter process is carried out to the point.
It is right according to the when window delimited for the point for the every bit on the lineups for each lineups
The point carries out one-dimensional median filter process.One-dimensional median filter process can be implemented using above formula (1),
I.e.：
G (x)=med f (x-k) ..., f (x) ..., f (x+k), k ∈ N } (1)
Wherein, g (x) is the output of one-dimensional medium filtering, and f (x) is to need to carry out the original signal of one-dimensional medium filtering, med
The intermediate value after sequence is represented, N is positive integer, represents the number for the interior point of when window delimited in step 2.
One-dimensional median filter process is carried out successively to multiple lineups of multigroup wave field, so as to reach while eliminate
Multiple wave fields, retain the purpose of single wave field, obtain the data after wave field separation.
Fig. 8 shows the wave field separation method using the VSP data according to exemplary embodiment, obtained ripple
Data after the separation of field.From figure 8, it is seen that what is obtained after wave field separation is single up reflected wave field,
The lineups of intersecting therewith representative remaining wave field knitted while eliminated.Illustrate that this method can be efficient
Wave field separation is carried out, effect meets the requirement of processing in next step.
, it is necessary to pick up the same phase of a wave field first in traditional wave field separation method based on medium filtering
Direction of principal axis, it is level then to draw lineups, then carries out median filter process, finally carries out reverse drawing processing again
Revert to the trace gather after separation.
Multiple wave fields of VSP data are distinguished according to the wave field separation method of the VSP data of exemplary embodiment
Carry out along directional medial filtering, and wave field separation can be carried out along multiple directions simultaneously in same trace gather,
Avoid evening up lineups and reverse drawing is handled, and then improve the precision and efficiency of wave field separation, be follow-up
The data basis provided is provided, is easy to subsequent treatment to obtain good effect.
Above-mentioned technical proposal is a kind of embodiment of the present invention, for those skilled in the art,
On the basis of principle disclosed by the invention, it is easy to make various types of improvement or deformation, and not only limit
In the description of the above-mentioned specific embodiment of the present invention, therefore description above is simply preferable, and and without limit
The meaning of property processed.
Claims (4)
1. a kind of wave field separation method of VSP data, comprises the following steps：
Step 1：Enter traveling-wave field pickup to pending VSP data, obtain the information of multigroup wave field, record
The position of the lineups of every group of wave field；
Step 2：For each lineups, for the every bit on the lineups, according to the position of the point
For the dot-dash timing window；
Step 3：For each lineups, for the every bit on the lineups, delimited according to for the point
When window, one-dimensional median filter process is carried out to the point.
2. the wave field separation method of VSP data according to claim 1, wherein according to below equation (1)
Carry out one-dimensional median filter process：
G (x)=med f (x-k) ..., f (x) ..., f (x+k), k ∈ N } (1)
Wherein, g (x) is the output of one-dimensional medium filtering, and f (x) is to need to carry out the original signal of one-dimensional medium filtering, med
The intermediate value after sequence is represented, N is positive integer, the number put when representing described in window.
3. the wave field separation method of VSP data according to claim 2, wherein N be equal to 5,7 or
9。
4. the wave field separation method of VSP data according to claim 1, wherein for each lineups,
For the every bit on the lineups, centered on the position of the point, above and below the point
Choose the dot-dash timing window of same number.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610331007.0A CN107402404B (en) | 2016-05-18 | 2016-05-18 | VSP data wave field separation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610331007.0A CN107402404B (en) | 2016-05-18 | 2016-05-18 | VSP data wave field separation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107402404A true CN107402404A (en) | 2017-11-28 |
| CN107402404B CN107402404B (en) | 2019-12-17 |
Family
ID=60394538
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610331007.0A Active CN107402404B (en) | 2016-05-18 | 2016-05-18 | VSP data wave field separation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107402404B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1404582A (en) * | 2000-01-21 | 2003-03-19 | 施鲁博格控股有限公司 | System and method for seismic wave field separation |
| CN101246218A (en) * | 2007-02-15 | 2008-08-20 | 中国石油化工股份有限公司 | Three-component VSP wave field separation method |
| CN102193107A (en) * | 2010-03-05 | 2011-09-21 | 西安石油大学 | Method for separating and denoising seismic wave field |
| US8077544B2 (en) * | 2008-03-28 | 2011-12-13 | Westerngeco L.L.C. | Dual-wavefield multiple attenuation |
| CN103149592A (en) * | 2013-03-07 | 2013-06-12 | 天津城市建设学院 | Method for separating variable offset vertical seismic profile (VSP) wave fields |
| CN104977615A (en) * | 2015-07-01 | 2015-10-14 | 中国海洋石油总公司 | Model-statistics-pickup-based multiple suppression method of deep-sea OBC data |
-
2016
- 2016-05-18 CN CN201610331007.0A patent/CN107402404B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1404582A (en) * | 2000-01-21 | 2003-03-19 | 施鲁博格控股有限公司 | System and method for seismic wave field separation |
| CN101246218A (en) * | 2007-02-15 | 2008-08-20 | 中国石油化工股份有限公司 | Three-component VSP wave field separation method |
| US8077544B2 (en) * | 2008-03-28 | 2011-12-13 | Westerngeco L.L.C. | Dual-wavefield multiple attenuation |
| CN102193107A (en) * | 2010-03-05 | 2011-09-21 | 西安石油大学 | Method for separating and denoising seismic wave field |
| CN103149592A (en) * | 2013-03-07 | 2013-06-12 | 天津城市建设学院 | Method for separating variable offset vertical seismic profile (VSP) wave fields |
| CN104977615A (en) * | 2015-07-01 | 2015-10-14 | 中国海洋石油总公司 | Model-statistics-pickup-based multiple suppression method of deep-sea OBC data |
Non-Patent Citations (2)
| Title |
|---|
| 李小红等: "中值在地震资料处理中的应用", 《石油物探》 * |
| 聂爱兰: "VSP波场分离的径向中值滤波方法", 《煤田地质与勘探》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107402404B (en) | 2019-12-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106846344B (en) | A kind of image segmentation optimal identification method based on the complete degree in edge | |
| CN102999886B (en) | Image Edge Detector and scale grating grid precision detection system | |
| CN104390981B (en) | The method of quick and precisely focusing scanning pathological section tissue based on image collecting device | |
| CN103886553B (en) | Method and system for non-local average value denoising of image | |
| CN102436650A (en) | image processing apparatus and image processing method | |
| CN103593850A (en) | SIFT parallelization system and method based on recursion Gaussian filtering on CUDA platform | |
| CN105513066B (en) | It is a kind of that the generic object detection method merged with super-pixel is chosen based on seed point | |
| CN103033134B (en) | For detection system and the method for damaged optical fiber geometric parameter | |
| WO2023070593A1 (en) | Line width measurement method and apparatus, computing processing device, computer program, and computer readable medium | |
| CN110726998B (en) | Method for measuring mining subsidence basin in mining area through laser radar scanning | |
| CN103743750B (en) | A kind of generation method of distribution diagram of surface damage of heavy calibre optical element | |
| CN103679672A (en) | Panorama image splicing method based on edge vertical distance matching | |
| CN108550166A (en) | A kind of spatial target images matching process | |
| CN103136752A (en) | Image magnification method based on edge extraction | |
| CN1703095B (en) | Apparatus and method for detecting letter box, and mpeg decoding device having the same | |
| CN102663733A (en) | Characteristic points matching method based on characteristic assembly | |
| CN112734822A (en) | Stereo matching algorithm based on infrared and visible light images | |
| CN106846354A (en) | A kind of Book Inventory method on frame converted based on image segmentation and random hough | |
| CN110599509B (en) | Edge detection method based on eight-direction fractional order differential operator | |
| JP2010016812A (en) | Image processing apparatus and method, and computer-readable medium | |
| CN110827355A (en) | Moving target rapid positioning method and system based on video image coordinates | |
| CN110706275A (en) | Local gravity anomaly extraction method based on satellite altimetry gravity data | |
| CN107402404A (en) | The wave field separation method of VSP data | |
| CN111368826B (en) | Open fire detection algorithm based on variable convolution kernel | |
| CN109215046B (en) | Laplacian operator edge detection method based on image interpolation operation |
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 |