CN107229073A - Seismic data processing technique and device - Google Patents
Seismic data processing technique and device Download PDFInfo
- Publication number
- CN107229073A CN107229073A CN201610172759.7A CN201610172759A CN107229073A CN 107229073 A CN107229073 A CN 107229073A CN 201610172759 A CN201610172759 A CN 201610172759A CN 107229073 A CN107229073 A CN 107229073A
- Authority
- CN
- China
- Prior art keywords
- mrow
- channel set
- seismic
- denoising
- seismic 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.)
- Pending
Links
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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
- G01V2210/24—Multi-trace filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
Abstract
The present invention relates to seismic prospecting data process field, seismic data processing technique and device are specifically disclosed.This method includes：The seismic channel set d (t, x) in t x domains is obtained, wherein, t represents TWT, and x represents offset distance；τ p conversion is carried out based on the weighting τ p transformation for mula after optimization.Using the present invention, random noise present in low SNR data can be substantially reduced, the signal to noise ratio of data is significantly improved, false lineups and aliasing present in existing τ p conversion process is eliminated, follow-up migration imaging result can be obviously improved.
Description
Technical field
The present invention relates to seismic prospecting data process field, more particularly, to a kind of seismic data processing technique and one
Plant seismic data process device.
Background technology
There is earth's surface topography and geomorphology, near-surface rock character condition in southern china and western exploration acreage and intricately descend ground
Texture is made simultaneously and the work area deposited, topographical elevation difference change acutely, low velocity layer thickness change greatly, high rock stratum crop out suddenly,
Subinverse covers that nappe structure is more, and these complicated seismic geological coditions propose bigger to existing earthquake-capturing and treatment technology
It is required that.Such complicated geological condition, the seismic signal signal to noise ratio collected is very low, and conventional denoising method and technology are not
It is applicable completely, with significant limitation and inadaptability.
τ-p conversion filtering techniques are a kind of conventional seismic data process means, still, in complicated mountain front exploratory area, are had
Effect signal is often submerged in noise, can't see effective lineups, or is only capable of seeing only a few lineups, and signal to noise ratio is extremely low
(general signal to noise ratio is less than 0.5), the data obtained after filtering technique processing are converted using existing τ-p and still suffer from false lineups and mixed
Folded phenomenon, it is impossible to the need for meeting seismic prospecting.
The content of the invention
The present invention proposes a kind of weighting τ-p conversion filtering techniques, and it can significantly improve the signal to noise ratio of data.
According to an aspect of the invention, it is proposed that a kind of seismic data processing technique, this method includes：Obtain t-x domains
Seismic channel set d (t, x), wherein, t represents TWT, and x represents offset distance；Obtained based on following formula after weighted τ-p conversion
τ-p numeric field datas m (τ, p)：
Wherein,N
The road number in the seismic channel set is represented, τ represents the intercept time, and p represents ray parameter, xnRepresent n-th in the seismic channel set
The offset distance in road, Δ x represents the difference of the offset distance of neighboring track.
According to another aspect of the invention, it is proposed that a kind of seismic data process device, the device includes：Seismic channel set is obtained
Unit is taken, the seismic channel set d (t, x) for obtaining t-x domains, wherein, t represents TWT, and x represents offset distance；Weight τ-p
Converter unit, for obtained based on following formula the τ-p numeric field datas m after weighted τ-p conversion (τ, p)：
Wherein,N
The road number in the seismic channel set is represented, τ represents the intercept time, and p represents ray parameter, xnRepresent n-th in the seismic channel set
The offset distance in road, Δ x represents the difference of the offset distance of neighboring track.
τ-p conversion filtering technique after weighting employed in the present invention, which can be substantially reduced in low SNR data, to be present
Random noise, significantly improve the signal to noise ratio of data, eliminate false lineups and aliasing present in existing τ-p conversion process existing
As follow-up migration imaging result can be obviously improved.Present invention is particularly suitable for the low-down seismic data of signal to noise ratio, typically
Refer to signal to noise ratio and be less than 0.5, existing commercial processes module is almost helpless to this, and apply the present invention, can be from this kind of extremely low letter
Make an uproar than searching effective seismic signal in data.
Brief description of the drawings
By the way that exemplary embodiment of the invention is described in more detail with reference to accompanying drawing, it is of the invention above-mentioned and its
Its purpose, feature and advantage will be apparent, wherein, in exemplary embodiment of the invention, identical reference number
Typically represent same parts.
Fig. 1 shows the flow chart of seismic data processing technique according to an embodiment of the invention.
Fig. 2 (a) is to cut open the velocity profile of the filter result obtained when signal to noise ratio is 0.6 using existing τ-p conversion and time
Face；Fig. 2 (b) is the velocity profile and time section of the filter result obtained to the identical logging data application present invention.
Fig. 3 (a) is to cut open the velocity profile of the filter result obtained when signal to noise ratio is 0.2 using existing τ-p conversion and time
Face；Fig. 3 (b) is the velocity profile and time section of the filter result obtained to the identical logging data application present invention.
Fig. 4 (a) is the velocity profile and time section for the original seismic data that signal to noise ratio is 1；Fig. 4 (b) is to Fig. 4 (a)
The velocity profile and time section for the filter result that the shown logging data application present invention is obtained.
Fig. 5 (a) is the velocity profile and time section for the original seismic data that signal to noise ratio is 0.6；Fig. 5 (b) is to Fig. 5
(a) velocity profile and time section for the filter result that the logging data application present invention shown in is obtained.
Fig. 6 (a) is the velocity profile and time section for the original seismic data that signal to noise ratio is 0.4；Fig. 6 (b) is to Fig. 6
(a) velocity profile and time section for the filter result that the logging data application present invention shown in is obtained.
Fig. 7 (a) is the velocity profile and time section for the original seismic data that signal to noise ratio is 0.2；Fig. 7 (b) is to Fig. 7
(a) velocity profile and time section for the filter result that the logging data application present invention shown in is obtained.
Fig. 8 (a) is the velocity profile and time section for the original seismic data that signal to noise ratio is 0.1；Fig. 8 (b) is to Fig. 8
(a) velocity profile and time section for the filter result that the logging data application present invention shown in is obtained.
Fig. 9 (a), (b), (c), (d) are followed successively by original CMP (CMP) seismic channel set, the filter obtained using the present invention
Data, the velocity profile for the noise being filtered out after ripple, the time section for the noise being filtered out.
Figure 10 (a), (b), (c) are Fig. 9 (a), (b), (d) corresponding frequency spectrum.
Figure 11 (a) is the eastern data region of no relief CMP trace gathers in the preceding storehouse of denoising；Figure 11 (b) is the eastern data region of no relief CMP in storehouse after denoising
Trace gather.
Figure 12 (a) is denoising preceding storehouse east data Ya Ken areas CMP trace gathers；Figure 12 (b) is storehouse east data Ya Ken area CMP after denoising
Trace gather.
Figure 13 (a) is the eastern data Tianshan District CMP trace gathers in the preceding storehouse of denoising；Figure 13 (b) is the eastern data Tianshan District CMP in storehouse after denoising
Trace gather.
Figure 14 (a) is stacked section before the eastern data denoising in storehouse；Figure 14 (b) is stacked section after the eastern data denoising in storehouse.
Figure 15 (a) is the eastern data Tianshan Mountainous stacked section amplification in the preceding storehouse of denoising；Figure 15 (b) is the eastern data day in storehouse after denoising
The stacked section amplification of mountain mountain area.
Embodiment
The preferred embodiment of the present invention is more fully described below with reference to accompanying drawings.Although showing the present invention in accompanying drawing
Preferred embodiment, however, it is to be appreciated that may be realized in various forms the present invention without the embodiment party that should be illustrated here
Formula is limited.On the contrary, these embodiments are provided so that the present invention is more thorough and complete, and can be by the present invention's
Scope intactly conveys to those skilled in the art.
Embodiment 1
Fig. 1 shows the flow chart of seismic data processing technique according to an embodiment of the invention.In the present embodiment
In, this method includes：
Step 101, the seismic channel set d (t, x) in t-x domains is obtained, wherein, t represents TWT, and x represents offset distance；
Step 102, based on following formula obtain the τ-p numeric field datas m after weighted τ-p conversion (τ, p)：
Wherein,
N represents the road number in the seismic channel set, and τ represents the intercept time, and p represents ray parameter, xnRepresent the earthquake
The offset distance in the n-th road in trace gather, Δ x represents the difference of the offset distance of neighboring track.
In the present embodiment, filtering technique (such as shown in formula (1)) is converted using the τ-p after weighting, compared to existing τ-p
Filtering technique (such as shown in formula (2)) is converted, random noise present in low SNR data can be substantially reduced, significantly improved
The signal to noise ratio of data, eliminates false lineups and aliasing present in existing τ-p conversion process, can be obviously improved follow-up
Migration imaging result.
In the present embodiment, this method can also include：Can to the τ-p numeric field datas m (τ, p) carry out denoising, obtain
τ-p numeric field datas after to denoising.
Further, this method can also include：τ-p numeric field datas after the denoising can be carried out corresponding anti-
τ-p convert (the anti-τ-p corresponding with weighting τ-p conversion is converted) to reconstruct seismic channel set；And can the earthquake based on reconstruct
Trace gather carries out migration imaging.
In the present embodiment, the seismic channel set used can be CMP seismic channel sets.
Embodiment 2
The invention also discloses a kind of seismic data process device.The device includes：Seismic channel set acquiring unit, for
To the seismic channel set d (t, x) in t-x domains, wherein, t represents TWT, and x represents offset distance；τ-p converter units are weighted, are used for
Based on following formula obtain the τ-p numeric field datas m after weighted τ-p conversion (τ, p)：
Wherein,N
The road number in the seismic channel set is represented, τ represents the intercept time, and p represents ray parameter, xnRepresent n-th in the seismic channel set
The offset distance in road, Δ x represents the difference of the offset distance of neighboring track.
In a kind of possible embodiment, the device can also include：Denoising unit, can be used for the τ-p domains
(τ p) carries out denoising, obtains the τ-p numeric field datas after denoising data m.
In a kind of possible embodiment, the device can also include：Anti- τ-p converter unit units, can be used for pair
τ-p numeric field datas after the denoising carry out corresponding anti-τ-p conversion to reconstruct seismic channel set；Migration imaging unit, can be with
Migration imaging is carried out for the seismic channel set based on reconstruct.
In a kind of possible embodiment, the seismic channel set is CMP trace gathers
Using example
For ease of understanding the scheme and its effect of the embodiment of the present invention, a concrete application example given below.This area
It should be understood to the one skilled in the art that the example is only for the purposes of understanding the present invention, its any detail is not intended to be limited in any way
The system present invention.
Fig. 2 (a) is to cut open the velocity profile of the filter result obtained when signal to noise ratio is 0.6 using existing τ-p conversion and time
Face.From Fig. 2 (a) velocity profile, many many false lineups, it can be seen that false same phase into strips on its time section
Axle, filter effect is unsatisfactory.Fig. 2 (b) be the velocity profile of filter result that is obtained to the identical logging data application present invention and when
Between section.Energy group is concentrated very much after being filtered from Fig. 2 (b) velocity profile, and without false energy group, time section does not have yet
There are false lineups.
Fig. 3 (a) is to cut open the velocity profile of the filter result obtained when signal to noise ratio is 0.2 using existing τ-p conversion and time
Face.From Fig. 3 (a) velocity profile, (such as Fig. 2 (a)) adds more false lineups during compared with signal to noise ratio 0.6, in the time
Also there are obvious stronger false lineups on section.Fig. 3 (b) is the speed of the filter result obtained to the identical logging data application present invention
Spend section and time section.Energy group is concentrated very much after being filtered from Fig. 3 (b) velocity profile, the energy group without vacation, when
Between section also still without false lineups.
It can be seen that significantly improving filter effect according to the τ-p filtering techniques after the optimization of the present invention, significantly carrying
While high s/n ratio, without aliasing, filter result is true and reliable.
Following a series of experiments can be used for the performance of the further checking present invention.Fig. 4 (a) be signal to noise ratio be 1 it is original
The velocity profile and time section of seismic data.Fig. 4 (b) is the filtering knot obtained to the logging data application present invention shown in Fig. 4 (a)
The velocity profile and time section of fruit.As can be seen that filter effect is very good, filtered seismic signal lineups are continuous, make an uproar
Sound is almost pressed completely, and the energy group of useful signal is concentrated on velocity profile.
Reduce signal to noise ratio, shown in such as Fig. 5 (a), it illustrates signal to noise ratio for the velocity profile of 0.6 original seismic data and
Time section.Fig. 5 (b) is to cut open the velocity profile of the filter result obtained to the logging data application present invention shown in Fig. 5 (a) and time
Face, it can be seen that filter effect is equally fine.
Signal to noise ratio is further reduced, shown in such as Fig. 6 (a), it illustrates speed of the signal to noise ratio for 0.4 original seismic data
Section and time section.As can be seen that useful signal is hardly visible in Fig. 6 (a) on velocity profile.Fig. 6 (b) is to Fig. 6
(a) velocity profile and time section for the filter result that the logging data application present invention shown in is obtained.As can be seen that having in Fig. 6 (b)
Effect signal is highlighted, and signal to noise ratio is significantly improved.
Signal to noise ratio is further reduced, shown in such as Fig. 7 (a), it illustrates speed of the signal to noise ratio for 0.2 original seismic data
Section and time section.Now effective seismic signal is can't see completely on seismic traces collection.Fig. 7 (b) is to Fig. 7
(a) velocity profile and time section for the filter result that the logging data application present invention shown in is obtained.As can be seen that being filtered in Fig. 7 (b)
Trace gather effect after ripple is still good.
Signal to noise ratio is further reduced, shown in such as Fig. 8 (a), it illustrates speed of the signal to noise ratio for 0.1 original seismic data
Section and time section.Fig. 8 (b) be the velocity profile of filter result that is obtained to the logging data application present invention shown in Fig. 8 (a) and
Time section.As can be seen that filter effect now is insecure, because having in the velocity profile of filtered earthquake record
The energy group of effect seismic signal is offset.
The experimental result of summary difference signal to noise ratio model, it can be seen that at least at signal to noise ratio as little as 0.2, application
The present invention remains to filtering and obtains effective seismic signal.
Fig. 9 (a), (b), (c), (d) are followed successively by original CMP (CMP) seismic channel set, the filter obtained using the present invention
Data, the velocity profile for the noise being filtered out after ripple, the time section for the noise being filtered out.As can be seen that the noise filtered out
In be practically free of any effective wave component, further demonstrate beneficial effects of the present invention.
Figure 10 (a), (b), (c) are Fig. 9 (a), (b), (d) corresponding frequency spectrum, are gone it can be seen that the present invention is good
Make an uproar performance.
Figure 11 (a) is the eastern data region of no relief CMP trace gathers in the preceding storehouse of denoising；Figure 11 (b) is the eastern data region of no relief CMP in storehouse after denoising
Trace gather.
Figure 12 (a) is denoising preceding storehouse east data Ya Ken areas CMP trace gathers；Figure 12 (b) is storehouse east data Ya Ken area CMP after denoising
Trace gather.
Figure 13 (a) is the eastern data Tianshan District CMP trace gathers in the preceding storehouse of denoising；Figure 13 (b) is the eastern data Tianshan District CMP in storehouse after denoising
Trace gather.
Figure 14 (a) is stacked section before the eastern data denoising in storehouse；Figure 14 (b) is stacked section after the eastern data denoising in storehouse.
Figure 15 (a) is the eastern data Tianshan Mountainous stacked section amplification in the preceding storehouse of denoising；Figure 15 (b) is the eastern data day in storehouse after denoising
The stacked section amplification of mountain mountain area.
The present invention can be system, method and/or computer program product.Computer program product can include computer
Readable storage medium storing program for executing, containing for making processor realize the computer-readable program instructions of various aspects of the invention.
Computer-readable recording medium can keep and store to perform the tangible of the instruction that equipment is used by instruction
Equipment.Computer-readable recording medium for example can be-- but be not limited to-- storage device electric, magnetic storage apparatus, optical storage
Equipment, electromagnetism storage device, semiconductor memory apparatus or above-mentioned any appropriate combination.Computer-readable recording medium
More specifically example (non exhaustive list) includes：Portable computer diskette, hard disk, random access memory (RAM), read-only deposit
It is reservoir (ROM), erasable programmable read only memory (EPROM or flash memory), static RAM (SRAM), portable
Compact disk read-only storage (CD-ROM), digital versatile disc (DVD), memory stick, floppy disk, mechanical coding equipment, for example thereon
Be stored with instruction punch card or groove internal projection structure and above-mentioned any appropriate combination.It is used herein above to calculate
Machine readable storage medium storing program for executing is not construed as instantaneous signal in itself, the electromagnetic wave of such as radio wave or other Free propagations, logical
Cross the electromagnetic wave (for example, the light pulse for passing through fiber optic cables) of waveguide or the propagation of other transmission mediums or transmitted by electric wire
Electric signal.
Computer-readable program instructions as described herein can be downloaded to from computer-readable recording medium each calculate/
Processing equipment, or outer computer is downloaded to or outer by network, such as internet, LAN, wide area network and/or wireless network
Portion's storage device.Network can be transmitted, be wirelessly transferred including copper transmission cable, optical fiber, router, fire wall, interchanger, gateway
Computer and/or Edge Server.Adapter or network interface in each calculating/processing equipment are received from network to be counted
Calculation machine readable program instructions, and the computer-readable program instructions are forwarded, for the meter being stored in each calculating/processing equipment
In calculation machine readable storage medium storing program for executing.
For perform the computer program instructions that operate of the present invention can be assembly instruction, instruction set architecture (ISA) instruction,
Machine instruction, machine-dependent instructions, microcode, firmware instructions, condition setup data or with one or more programming languages
Source code or object code that any combination is write, programming language of the programming language including object-oriented-such as
Smalltalk, C++ etc., and conventional procedural programming languages-such as " C " language or similar programming language.Computer
Readable program instructions can perform fully on the user computer, partly perform on the user computer, as one solely
Vertical software kit is performed, part is performed or completely in remote computer on the remote computer on the user computer for part
Or performed on server.In the situation of remote computer is related to, remote computer can be by network-bag of any kind
LAN (LAN) or wide area network (WAN)-be connected to subscriber computer are included, or, it may be connected to outer computer is (such as sharp
With ISP come by Internet connection).In certain embodiments, by using computer-readable program instructions
Status information carry out personalized customization electronic circuit, such as PLD, field programmable gate array (FPGA) or can
Programmed logic array (PLA) (PLA), the electronic circuit can perform computer-readable program instructions, so as to realize each side of the present invention
Face.
Referring herein to method according to embodiments of the present invention, device (system) and computer program product flow chart and/
Or block diagram describes various aspects of the invention.It should be appreciated that each square frame and flow chart of flow chart and/or block diagram and/
Or in block diagram each square frame combination, can be realized by computer-readable program instructions.
These computer-readable program instructions can be supplied to all-purpose computer, special-purpose computer or other programmable datas
The processor of processing unit, so as to produce a kind of machine so that these instructions are passing through computer or other programmable datas
During the computing device of processing unit, work(specified in one or more of implementation process figure and/or block diagram square frame is generated
The device of energy/action.Can also be the storage of these computer-readable program instructions in a computer-readable storage medium, these refer to
Order causes computer, programmable data processing unit and/or other equipment to work in a specific way, so that, be stored with instruction
Computer-readable medium then includes a manufacture, and it is included in one or more of implementation process figure and/or block diagram square frame
The instruction of the various aspects of defined function/action.
Computer-readable program instructions can also be loaded into computer, other programmable data processing units or other
In equipment so that perform series of operation steps on computer, other programmable data processing units or miscellaneous equipment, to produce
Raw computer implemented process, so that performed on computer, other programmable data processing units or miscellaneous equipment
Instruct function/action specified in one or more of implementation process figure and/or block diagram square frame.
Flow chart and block diagram in accompanying drawing show system, method and the computer journey of multiple embodiments according to the present invention
Architectural framework in the cards, function and the operation of sequence product.At this point, each square frame in flow chart or block diagram can generation
One module of table, program segment or a part for instruction, the module, program segment or a part for instruction are used comprising one or more
In the executable instruction for realizing defined logic function.In some realizations as replacement, the function of being marked in square frame
Can be with different from the order marked in accompanying drawing generation.For example, two continuous square frames can essentially be held substantially in parallel
OK, they can also be performed in the opposite order sometimes, and this is depending on involved function.It is also noted that block diagram and/or
The combination of each square frame in flow chart and the square frame in block diagram and/or flow chart, can use function as defined in execution or dynamic
The special hardware based system made is realized, or can be realized with the combination of specialized hardware and computer instruction.
It is described above various embodiments of the present invention, described above is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.In the case of without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes will be apparent from for the those of ordinary skill in art field.For example, those of ordinary skill in the art
It is understood that disclosed in this invention technical scheme can be applied not only to seismic prospecting data process field, may be used also
With applied in various signal transactings and image procossing, such as signal of communication processing, medical signals processing and image procossing etc.,
These are only the simple deformations of the embodiment to being enumerated in the present invention, belong to same invention with the embodiment that is enumerated in the present invention
Design, also belongs to protection scope of the present invention.The selection of term used herein, it is intended to best explain the original of each embodiment
Reason, practical application or the improvement to the technology in market, or other those of ordinary skill of the art is understood that this
Each embodiment that text is disclosed.
Claims (8)
1. a kind of seismic data processing technique, this method includes：
The seismic channel set d (t, x) in t-x domains is obtained, wherein, t represents TWT, and x represents offset distance；
Based on following formula obtain the τ-p numeric field datas m after weighted τ-p conversion (τ, p)：
<mrow>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msubsup>
<mi>m</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>N</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
Wherein,N is represented
Road number in the seismic channel set, τ represents the intercept time, and p represents ray parameter, xnRepresent the n-th road in the seismic channel set
Offset distance, Δ x represents the difference of the offset distance of neighboring track.
2. seismic data processing technique according to claim 1, this method also includes：
To the τ-p numeric field datas m, (τ p) carries out denoising, obtains the τ-p numeric field datas after denoising.
3. seismic data processing technique according to claim 2, this method also includes：
Corresponding anti-τ-p conversion is carried out to the τ-p numeric field datas after the denoising to reconstruct seismic channel set；
Seismic channel set based on reconstruct carries out migration imaging.
4. the seismic data processing technique according to any one in claim 1-3, wherein, the seismic channel set is CMP
Trace gather.
5. a kind of seismic data process device, the device includes：
Seismic channel set acquiring unit, the seismic channel set d (t, x) for obtaining t-x domains, wherein, t represents TWT, and x is represented
Offset distance；
Weight τ-p converter units, for obtained based on following formula the τ-p numeric field datas m after weighted τ-p conversion (τ, p)：
<mrow>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msubsup>
<mi>m</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>N</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
Wherein,N is represented
Road number in the seismic channel set, τ represents the intercept time, and p represents ray parameter, xnRepresent the n-th road in the seismic channel set
Offset distance, Δ x represents the difference of the offset distance of neighboring track.
6. seismic data process device according to claim 5, the device also includes：
Denoising unit, for (τ p) to carry out denoising, obtains the τ-p numeric field datas after denoising to the τ-p numeric field datas m.
7. seismic data process device according to claim 6, the device also includes：
Anti- τ-p converter unit units, for carrying out corresponding anti-τ-p conversion to the τ-p numeric field datas after the denoising with weight
Structure seismic channel set；
Migration imaging unit, migration imaging is carried out for the seismic channel set based on reconstruct.
8. the seismic data process device according to any one in claim 5-7, wherein, the seismic channel set is CMP
Trace gather.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610172759.7A CN107229073A (en) | 2016-03-24 | 2016-03-24 | Seismic data processing technique and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610172759.7A CN107229073A (en) | 2016-03-24 | 2016-03-24 | Seismic data processing technique and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107229073A true CN107229073A (en) | 2017-10-03 |
Family
ID=59931874
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610172759.7A Pending CN107229073A (en) | 2016-03-24 | 2016-03-24 | Seismic data processing technique and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107229073A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110082823A (en) * | 2019-05-09 | 2019-08-02 | 中国石油大学(北京) | A kind of seismic data interpolation method and device |
| CN111929731A (en) * | 2020-07-21 | 2020-11-13 | 中国海洋大学 | Surface consistency and non-consistency combined static correction method |
| CN111983685A (en) * | 2020-07-21 | 2020-11-24 | 中国海洋大学 | Tau-p domain surface non-uniformity long wavelength static correction method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103869364A (en) * | 2014-03-25 | 2014-06-18 | 中国石油大学(华东) | Multiple wave suppression method based on dual parabolic Radon transformation |
| CN104502974A (en) * | 2014-12-31 | 2015-04-08 | 中国石油天然气集团公司 | Suppressing multiple reflection combining method and device |
-
2016
- 2016-03-24 CN CN201610172759.7A patent/CN107229073A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103869364A (en) * | 2014-03-25 | 2014-06-18 | 中国石油大学(华东) | Multiple wave suppression method based on dual parabolic Radon transformation |
| CN104502974A (en) * | 2014-12-31 | 2015-04-08 | 中国石油天然气集团公司 | Suppressing multiple reflection combining method and device |
Non-Patent Citations (8)
| Title |
|---|
| ORHAN YILMAZ ET AL.: "Discrete plane-wave decomposition by least-mean-square-error method", 《GEOPHYSICS》 * |
| PAUL L. STOFFA ET AL.: "Direct mapping of seismic data to the domain of intercept time and ray parameter -A plane-wave decomposition", 《GEOPHYSICS》 * |
| 巩向博等: "各向异性Radon变换及其在多次波压制中的应用", 《地球物理学报》 * |
| 殷八斤等: "利用局部τ-p变换增强叠后地震信号", 《石油地球物理勘探》 * |
| 罗小佳等: "最优化相似加权Radon变换压制多次波", 《勘探地球物理进展》 * |
| 谢金娥等: "一种优化预测滤波压制随机噪声的方法", 《地球物理学进展》 * |
| 陈伟: "微地震波场分离技术研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 鲁娥等: "混合Ｒadon变换地震噪声压制的应用", 《物探与化探》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110082823A (en) * | 2019-05-09 | 2019-08-02 | 中国石油大学(北京) | A kind of seismic data interpolation method and device |
| CN110082823B (en) * | 2019-05-09 | 2020-08-14 | 中国石油大学(北京) | Seismic data interpolation method and device |
| CN111929731A (en) * | 2020-07-21 | 2020-11-13 | 中国海洋大学 | Surface consistency and non-consistency combined static correction method |
| CN111983685A (en) * | 2020-07-21 | 2020-11-24 | 中国海洋大学 | Tau-p domain surface non-uniformity long wavelength static correction method |
| CN111929731B (en) * | 2020-07-21 | 2021-09-28 | 中国海洋大学 | Surface consistency and non-consistency combined static correction method |
| CN111983685B (en) * | 2020-07-21 | 2021-11-12 | 中国海洋大学 | Static correction method for tau-p domain surface non-uniformity |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hata et al. | A practical method to estimate strong ground motions after an earthquake, based on site amplification and phase characteristics | |
| CN108037531B (en) | A kind of seismic inversion method and system based on the full variational regularization of broad sense | |
| CN105510976B (en) | A kind of many subwaves combination adaptive attenuation method | |
| CN110178056A (en) | Virtual source denoising is carried out using small echo cross-correlation | |
| CN103926623B (en) | Method for suppressing reverse time migration low frequency noise | |
| CN107272062A (en) | A kind of underground medium Q methods of estimation of data-driven | |
| CN107229073A (en) | Seismic data processing technique and device | |
| CN107340536A (en) | The method and apparatus that seismic inversion peels off coal seam strong reflection | |
| CN105676291A (en) | Multiple wave matching attenuation method based on optimized phase axis tracking | |
| CN107179543A (en) | The method and apparatus that regularization is carried out to pre stack data | |
| DeShon et al. | Summary of the North Texas earthquake study seismic networks, 2013–2018 | |
| CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
| CN105629300B (en) | The method for improving complicated structure offset data signal-to-noise ratio | |
| CN107179550A (en) | A kind of seismic signal zero phase deconvolution method of data-driven | |
| CN106019376A (en) | Seismic wave compensation method constructed through frequency driving space-variant Q value model | |
| CN106873036A (en) | A kind of denoising method combined based on well shake | |
| CN105044769B (en) | The method for improving the resolution ratio of seismic signal | |
| CN106324669A (en) | Method of separating surface-related multiples of different orders in seismic exploration data | |
| CN105510975A (en) | Method and device for improving signal-to-noise ratio of seismic data | |
| CN105277986A (en) | Vibroseis harmonic wave suppressing method based on adaptive matching filter operator | |
| CN107255833B (en) | The determination method and apparatus of Seismic Stacked Section | |
| CN105093282A (en) | Energy replacement surface wave suppressing method based on frequency constraint | |
| Bhattacharyya | Design of spatial filters and their application to high-resolution aeromagnetic data | |
| Zhao et al. | Background noise suppression for DAS-VSP records using GC-AB-unet | |
| Yang et al. | Automatically Extracting Surface‐Wave Group and Phase Velocity Dispersion Curves from Dispersion Spectrograms Using a Convolutional Neural Network |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171003 |
|
| RJ01 | Rejection of invention patent application after publication |