CN109143344A - Diffracted wave separation method, device and electronic equipment - Google Patents
Diffracted wave separation method, device and electronic equipment Download PDFInfo
- Publication number
- CN109143344A CN109143344A CN201811358076.6A CN201811358076A CN109143344A CN 109143344 A CN109143344 A CN 109143344A CN 201811358076 A CN201811358076 A CN 201811358076A CN 109143344 A CN109143344 A CN 109143344A
- Authority
- CN
- China
- Prior art keywords
- sampled point
- point
- trace gather
- zero
- diffracted wave
- 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
- 238000000926 separation method Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 53
- 238000000605 extraction Methods 0.000 claims abstract description 23
- 229910052704 radon Inorganic materials 0.000 claims abstract description 14
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 14
- 230000002123 temporal effect Effects 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 8
- 238000001914 filtration Methods 0.000 claims description 8
- 239000000284 extract Substances 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 238000010030 laminating Methods 0.000 claims description 2
- 230000002547 anomalous effect Effects 0.000 description 6
- 230000006854 communication Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000003325 tomography 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)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides a kind of diffracted wave separation method, device and electronic equipments, wherein this method comprises: obtaining zero-offset trace gather；Zero-offset trace gather is handled, the instantaneous phase of seismic wave is obtained；Instantaneous phase includes multiple sampled points；Following steps are executed to each sampled point: sampled point being converted, the corresponding multiple energy values of sampled point are obtained；According to multiple energy values, the maximum energy value of sampled point is determined；The location information of Diffraction Point is determined according to maximum energy value；It is extracted according to diffracted wave of the location information to Diffraction Point, obtains extracting result；The extraction result of multiple sampled points is overlapped, the diffracted wave after being separated.The present invention passes through the hyperbolic feature using diffracted wave travel time curve, using improved mobile vertex Radon transform method, accurately obtain the location information of Diffraction Point, and diffracted wave is effectively extracted, effectively compacting noise and back wave can preferably realize guarantor's polarity separation of diffracted wave without velocity information.
Description
Technical field
The present invention relates to seismic exploration technique fields, set more particularly, to a kind of diffracted wave separation method, device and electronics
It is standby.
Background technique
Seismic wave is in communication process, can be with when encountering the geological anomalous bodies such as breakpoint, depositional termination line, crack, karst collapse col umn
The form continuation of spherical wave is propagated forward, forms diffracted wave, diffracted wave information can be used to realize the high-resolution of geological anomalous body
Detection has important research significance in the exploitation of the mineral resources such as coal, petroleum.
But diffraction wave energy is weaker, is often suppressed by back wave, it is difficult to ideal imaging results is obtained, using plane wavelength-division
Solution technology, processing accuracy are difficult to accurately hold, and there are noise jamming, it is difficult to be accurately positioned to Diffraction Point, utilize diffracted wave
Kinematics and dynamic characteristic handled, it is difficult to distinguish edge diffraction and point diffraction.
Summary of the invention
In view of this, the purpose of the present invention is to provide a kind of diffracted wave separation method, device and electronic equipments, with accurate
Obtain the location information of Diffraction Point, and diffracted wave effectively extracted, effectively suppress noise and back wave, it is preferable realize around
Guarantor's polarity of ejected wave separates.
In a first aspect, the embodiment of the invention provides a kind of diffracted wave separation methods, wherein include: acquisition zero-offset
Trace gather；Zero-offset trace gather is handled, the instantaneous phase of seismic wave is obtained；Instantaneous phase includes multiple sampled points；To every
A sampled point executes following steps: converting to sampled point, obtains the corresponding multiple energy values of sampled point；According to multiple energy
Value, determines the maximum energy value of sampled point；The location information of Diffraction Point is determined according to maximum energy value；According to location information to around
The diffracted wave of exit point extracts, and obtains extracting result；The extraction result of multiple sampled points is overlapped, after being separated
Diffracted wave.
With reference to first aspect, the embodiment of the invention provides the first possible embodiments of first aspect, wherein right
The zero-offset trace gather is handled, and the instantaneous phase of seismic wave is obtained, comprising: using median filtering method to zero-offset road
Collection is handled, the zero-offset trace gather that obtains that treated；To treated, zero-offset trace gather carries out Hilbert transform, obtains
To the orthogonal trace gather of treated zero-offset trace gather；According to treated zero-offset trace gather and orthogonal trace gather, earthquake is obtained
The instantaneous phase of wave.
With reference to first aspect, the embodiment of the invention provides second of possible embodiments of first aspect, wherein right
Sampled point is converted, and obtains the corresponding multiple energy values of sampled point, comprising: obtain multiple scanning curvature；According to multiple scannings
Curvature is utilized respectively improved mobile vertex Radon transform method and converts to sampled point, it is corresponding to obtain each scanning curvature
Energy value.
The possible embodiment of second with reference to first aspect, the embodiment of the invention provides the third of first aspect
Possible embodiment, wherein improved mobile vertex Radon transform method are as follows:
Wherein, x0Indicate the coordinate of sampled point；The time of t expression sampled point；Q indicates scanning curvature；After τ indicates conversion
The time of sampled point；xiIndicate the coordinate of the remaining sampled point in addition to sampled point；M indicates the road number before sample point coordinate；N
Indicate the road number after sample point coordinate.
With reference to first aspect, the embodiment of the invention provides the 4th kind of possible embodiments of first aspect, wherein position
Confidence breath includes the Taoist monastic name information, temporal information and scanning curvature information of Diffraction Point.
The 4th kind of possible embodiment with reference to first aspect, the embodiment of the invention provides the 5th kind of first aspect
Possible embodiment, wherein extracted according to diffracted wave of the location information to Diffraction Point, obtain extracting result, comprising: root
According to Taoist monastic name information, double branch Radon transforms are carried out using the first calculation formula；First calculation formula are as follows:
Wherein, x0Indicate the coordinate of the corresponding sampled point of location information；T indicate the corresponding sampled point of location information when
Between；xiIndicate the coordinate of the remaining sampled point in addition to sampled point；M indicates the road number before sample point coordinate；N indicates sampled point
Road number after coordinate；
According to temporal information and scanning curvature information, the region at (T, Q) is carried out using the second calculation formula to draw east anti-
Transformation obtains extracting result；Second calculation formula are as follows:
Wherein, τ indicates temporal information T；Q indicates scanning curvature information Q；x0Indicate the seat of the corresponding sampled point of location information
Mark；M indicates the road number before sample point coordinate；N indicates the road number after sample point coordinate.
Second aspect, the embodiment of the present invention also provide a kind of diffracted wave separator, wherein include: acquisition module, be used for
Obtain zero-offset trace gather；Processing module obtains the instantaneous phase of seismic wave for handling zero-offset trace gather；Wink
Shi Xiangwei includes multiple sampled points；Extraction module, for executing following steps to each sampled point: sampled point is converted,
Obtain the corresponding multiple energy values of sampled point；According to multiple energy values, the maximum energy value of sampled point is determined；According to ceiling capacity
It is worth the location information for determining Diffraction Point；It is extracted according to diffracted wave of the location information to Diffraction Point, obtains extracting result；Superposition
Module, for the extraction result of multiple sampled points to be overlapped, the diffracted wave after being separated.
In conjunction with second aspect, the embodiment of the invention provides the first possible embodiments of second aspect, wherein place
Manage module, comprising: filter module obtains that treated zero for handling zero-offset trace gather using median filtering method
Offset gather；Conversion module obtains that treated zero for carrying out Hilbert transform to treated zero-offset trace gather
The orthogonal trace gather of offset gather；Computing module, for obtaining earthquake according to treated zero-offset trace gather and orthogonal trace gather
The instantaneous phase of wave.
The third aspect, the embodiment of the present invention also provide a kind of electronic equipment, including memory, processor, deposit in memory
Contain the computer program that can be run on a processor, wherein processor realizes above-mentioned first aspect when executing computer program
The step of described method.
Fourth aspect, the embodiment of the present invention also provide a kind of meter of non-volatile program code that can be performed with processor
Calculation machine readable medium, wherein program code makes processor execute above-mentioned first aspect the method.
The embodiment of the present invention bring it is following the utility model has the advantages that
The present invention provides a kind of diffracted wave separation method, device and electronic equipments, wherein this method comprises: obtaining zero
Offset gather；Zero-offset trace gather is handled, the instantaneous phase of seismic wave is obtained；Instantaneous phase includes multiple samplings
Point；Following steps are executed to each sampled point: sampled point being converted, the corresponding multiple energy values of sampled point are obtained；According to
Multiple energy values determine the maximum energy value of sampled point；The location information of Diffraction Point is determined according to maximum energy value；According to position
Information extracts the diffracted wave of Diffraction Point, obtains extracting result；The extraction result of multiple sampled points is overlapped, is obtained
Diffracted wave after separation.The present invention is by the hyperbolic feature using diffracted wave travel time curve, using improved mobile vertex La Dong
Converter technique accurately obtains the location information of Diffraction Point, and is effectively extracted to diffracted wave, effectively suppresses noise and back wave,
Guarantor's polarity separation of diffracted wave can be preferably realized without velocity information.
Other features and advantages of the present invention will illustrate in the following description, alternatively, Partial Feature and advantage can be with
Deduce from specification or unambiguously determine, or by implementing above-mentioned technology of the invention it can be learnt that.
To enable the above objects, features and advantages of the present invention to be clearer and more comprehensible, better embodiment is cited below particularly, and match
Appended attached drawing is closed, is described in detail below.
Detailed description of the invention
It, below will be to specific in order to illustrate more clearly of the specific embodiment of the invention or technical solution in the prior art
Embodiment or attached drawing needed to be used in the description of the prior art be briefly described, it should be apparent that, it is described below
Attached drawing is some embodiments of the present invention, for those of ordinary skill in the art, before not making the creative labor
It puts, is also possible to obtain other drawings based on these drawings.
Fig. 1 is a kind of flow chart of diffracted wave separation method provided in an embodiment of the present invention；
Fig. 2 is a kind of flow chart handled each sampled point provided in an embodiment of the present invention；
Fig. 3 is the flow chart of another diffracted wave separation method provided in an embodiment of the present invention；
Fig. 4 is a kind of structural schematic diagram of zero-offset trace gather provided in an embodiment of the present invention；
Fig. 5 is a kind of structure of the corresponding maximum energy value of multiple sampled points of instantaneous phase provided in an embodiment of the present invention
Schematic diagram；
Fig. 6 is the structural schematic diagram of the diffracted wave after a kind of separation provided in an embodiment of the present invention；
Fig. 7 is a kind of structural schematic diagram of diffracted wave separator provided in an embodiment of the present invention；
Fig. 8 is the structural schematic diagram of a kind of electronic equipment provided in an embodiment of the present invention.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with attached drawing to the present invention
Technical solution be clearly and completely described, it is clear that described embodiments are some of the embodiments of the present invention, rather than
Whole embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art are not making creative work premise
Under every other embodiment obtained, shall fall within the protection scope of the present invention.
Currently, there are noise jamming for isolated diffracted wave using decomposition of plane wave technology, and do not point out accurately
The position of Diffraction Point utilizes the kinematics and dynamic characteristic of diffracted wave, it is difficult to realize that the diffracted wave under polarity Reversion mentions
It takes, and edge diffraction and point diffraction cannot be distinguished, be based on this, a kind of diffracted wave separation method provided in an embodiment of the present invention, dress
It sets and electronic equipment, the diffracted wave that can be applied under polarity Reversion extracts.
For convenient for understanding the present embodiment, first to a kind of diffracted wave separation method disclosed in the embodiment of the present invention
It describes in detail.
The flow chart of a kind of diffracted wave separation method shown in Figure 1, wherein this method specific steps include:
Step S102 obtains zero-offset trace gather；
Seismic wave is in communication process, can be with spherical wave when encountering the geological anomalous bodies such as breakpoint, depositional termination line, crack
Form continuation propagate forward, form diffracted wave, diffracted wave carries the effective information of a large amount of geological anomalous bodies, but diffracted wave
Energy is weaker, is usually suppressed by back wave, in zero-offset trace gather, the travel time curve of back wave and underground continuous reflection interface
Shape it is similar, be mostly local linear, and the travel time curve of diffracted wave is then similar to hyperbola, it is right using this feature
Diffracted wave in zero-offset trace gather is separated.
Step S104 handles zero-offset trace gather, obtains the instantaneous phase of seismic wave；Instantaneous phase includes multiple
Sampled point；
The step of handling zero-offset trace gather, obtaining the instantaneous phase of seismic wave include:
Zero-offset trace gather is filtered using median filtering method, the zero-offset trace gather that obtains that treated；
Hilbert transform is carried out to treated zero-offset trace gather, the zero-offset trace gather that obtains that treated it is orthogonal
Trace gather；
According to treated zero-offset trace gather and orthogonal trace gather, the instantaneous phase of seismic wave is obtained.
Noise Elimination can be carried out to zero-offset trace gather using median filtering method, effectively reduce noise to subsequent processing
The influence of process, using Hilbert transform pairs, treated that zero-offset trace gather is converted, so as to the instantaneous of seismic wave
The definition of parameter and it is calculated as possibility, realizes the extraction of the instantaneous phase of seismic wave truly, the ground being calculated
The instantaneous phase of seismic wave includes multiple sampled points.
Step S106 handles each sampled point, obtains the extraction result of each sampled point；
Referring to fig. 2, it is described in detail to the step of in step S106, processing each sampled point, specific steps
It is as follows:
Step S202, converts sampled point, obtains the corresponding multiple energy values of sampled point；
Step S204 determines the maximum energy value of sampled point according to multiple energy values；
Step S206 determines the location information of Diffraction Point according to maximum energy value；
Step S208 is extracted according to diffracted wave of the location information to Diffraction Point, obtains extracting result.
Above-mentioned steps S202 describes to handle each sampled point to step S208, and it is corresponding to obtain each sampled point
Extraction result treatment process, continue with description diffracted wave separation method the step of.
The extraction result of multiple sampled points is overlapped, the diffracted wave after being separated by step S108.
The present invention provides a kind of diffracted wave separation methods, and this method is by obtaining zero-offset trace gather；To zero-offset
Trace gather is handled, and the instantaneous phase of seismic wave is obtained；Instantaneous phase includes multiple sampled points；Each sampled point is executed following
Step: converting sampled point, obtains the corresponding multiple energy values of sampled point；According to multiple energy values, sampled point is determined
Maximum energy value；The location information of Diffraction Point is determined according to maximum energy value；According to location information to the diffracted wave of Diffraction Point into
Row extracts, and obtains the corresponding extraction result of sampled point；The extraction result of multiple sampled points is overlapped, after being separated around
Ejected wave.The present invention is by the hyperbolic feature using diffracted wave travel time curve, using improved mobile vertex Radon transform method, accurately
The location information of Diffraction Point is obtained, and diffracted wave is effectively extracted, effectively suppresses noise and back wave, is not necessarily to velocity information
It can preferably realize guarantor's polarity separation of diffracted wave.
Corresponding to above method embodiment, the embodiment of the present invention processes diffracted wave separation by taking a sampled point as an example
Journey is described in detail, as shown in figure 3, specific steps include:
Step S302 obtains zero-offset trace gather；
Step S304 handles zero-offset trace gather, obtains the instantaneous phase of seismic wave, and instantaneous phase includes multiple
Sampled point；
Noise treatment is carried out to zero-offset trace gather shown in Fig. 4 using median filtering method, the zero offset that obtains that treated
Away from trace gather d (x, t), Hilbert transform is carried out to d (x, t), the orthogonal trace gather d* of the zero-offset trace gather that obtains that treated (x,
T), d (x, t) and d* (x, t) are handled using following formula, obtains the instantaneous phase of seismic wave, instantaneous phase includes more
A sampled point, the formula are as follows:
θ (x, t)=arctan (d (x, t)/d*(x,t))；
Wherein, x indicates the corresponding coordinate of sampled point；T indicates the sampled point corresponding time；θ (x, t) is the instantaneous of seismic wave
Phase.
Step S306 converts one of sampled point, obtains the corresponding multiple energy values of the sampled point；
With sampled point x0For be illustrated, which is considered as diffraction vertex (i.e. hyperbola vertex), utilizes improvement
Mobile vertex Radon transform method to sampled point x0Conversion process is carried out, specific conversion regime is as follows:
Instantaneous phase θ (x, t) is two-dimensional array, when being collected to form zero-offset trace gather to data, due to can not
Data collection is carried out to each position, one data collecting point generally is set every a distance, to be arranged every ten meters
For one data collecting point, the position x of sampled point can be obtained into the corresponding Taoist monastic name of the sampled point divided by 10.
When being extracted to sampled point, by sampled point x to be extracted0The Taoist monastic name in place road is set as 0, be located at the Taoist monastic name it
The Taoist monastic name in preceding road is then set as -1 accordingly, and -2, the Taoist monastic name of -3 ...,-M, the road after the Taoist monastic name are then set as 1 accordingly,
2 ..., N.
To sampled point x0When carrying out conversion process, multiple scanning curvature are first obtained, according to multiple scanning curvature, respectively
Using improved mobile vertex Radon transform method to sampled point x0It is converted, obtains the corresponding energy of each scanning curvature
Value, improved mobile vertex Radon transform method are as follows:
Wherein, x0Indicate the coordinate of sampled point；The time of t expression sampled point；Q indicates scanning curvature；After τ indicates conversion
The time of sampled point；xiIndicate the coordinate of the remaining sampled point in addition to sampled point；M indicates the road number before sample point coordinate；N
Indicate the road number after sample point coordinate.
Step S308 determines the corresponding maximum energy value of the sampled point according to multiple energy values；
According to multiple energy values, the corresponding maximum energy value E (τ, x0) of the sampled point is obtained；
E(τ,x0)=maxq{m(τ,q,x0)}。
Step S310 determines the location information of the corresponding Diffraction Point of the sampled point according to maximum energy value；
Wherein, location information includes the corresponding Taoist monastic name information TR of the sampled point, temporal information T and scanning curvature information Q；Road
The Taoist monastic name in the road where number information TR, that is, corresponding sampled point of maximum energy value；Temporal information T, that is, corresponding turn of maximum energy value
The time in the domain La Dong after changing.
Step S312 extracts according to diffracted wave of the location information of the Diffraction Point to the Diffraction Point, obtains the sampling
The corresponding extraction result of point；
It is extracted according to diffracted wave of the location information to Diffraction Point, obtains the step of extracting result, comprising:
According to Taoist monastic name information TR, double branch Radon transforms are carried out to the sampled point at TR using the first calculation formula；By time domain
It is converted into the domain La Dong, the first calculation formula are as follows:
Wherein, x0Indicate the coordinate of the corresponding sampled point of location information；T indicate the corresponding sampled point of location information when
Between；xiIndicate the coordinate of the remaining sampled point in addition to sampled point；M indicates the road number before sample point coordinate；N indicates sampled point
Road number after coordinate.
According to temporal information T and scanning curvature information Q, using the second calculation formula in the domain La Dong being converted to
Region at (T, Q) carries out drawing eastern inverse transformation, obtains extracting result；Second calculation formula are as follows:
Wherein, τ indicates temporal information T；Q indicates scanning curvature information Q；x0Indicate the seat of the corresponding sampled point of location information
Mark；M indicates the road number before sample point coordinate；N indicates the road number after sample point coordinate.
The corresponding extraction result of all sampled points of instantaneous phase is overlapped by step S314, after being separated around
Ejected wave.
Each sampled point of instantaneous phase is extracted in the manner described above, it is corresponding to obtain each sampled point
Maximum energy value, as shown in figure 5, carrying out double branch Radon transforms and La Dong inverse transformation to the corresponding Diffraction Point of these maximum energy values
To the corresponding extraction of each sampled point as a result, the corresponding extraction result of each sampled point is overlapped, can obtain
Diffracted wave after separation, as shown in Figure 6.
It, can also be further to tomography breakpoint, karst collapse col umn and other geology according to the characteristic information of the diffracted wave after separation
Anomalous body is identified and is detected.
The embodiment of the present invention by according to scanning curvature, using improved mobile vertex Radon transform method to instantaneous phase
Each sampled point converted, each corresponding energy value of scanning curvature of each sampled point is obtained, from sampled point
Maximum energy value is chosen in corresponding multiple energy values, using the corresponding location information of maximum energy value as the position of Diffraction Point
Time domain is converted to the domain La Dong at Diffraction Point by information, is chosen the region in the domain La Dong at Diffraction Point corresponding (T, Q) and is drawn
Eastern inverse transformation realizes the extraction to diffracted wave hyperbola left valve and right valve, using the separate equations pre-established, retains diffracted wave
Polarity, obtain the diffracted wave under polarity Reversion, and diffracted wave can preferably be separated without velocity information.
Corresponding to above method embodiment, the embodiment of the invention also provides a kind of diffracted wave separators, such as Fig. 7 institute
Show, wherein include:
Module 70 is obtained, for obtaining zero-offset trace gather；
Processing module 71 obtains the instantaneous phase of seismic wave for handling zero-offset trace gather；Instantaneous phase packet
Include multiple sampled points；
Extraction module 72, for executing following steps to each sampled point:
Sampled point is converted, the corresponding multiple energy values of sampled point are obtained；
According to multiple energy values, the maximum energy value of sampled point is determined；
The location information of Diffraction Point is determined according to maximum energy value；
It is extracted according to diffracted wave of the location information to Diffraction Point, obtains extracting result；
Laminating module 73, for the extraction result of multiple sampled points to be overlapped, the diffracted wave after being separated.
Above-mentioned processing module 71, further includes: filter module, for using median filtering method to zero-offset trace gather at
Reason, the zero-offset trace gather that obtains that treated；Conversion module, for zero-offset trace gather to carry out Hilbert change to treated
It changes, the orthogonal trace gather for the zero-offset trace gather that obtains that treated；Computing module, for according to treated zero-offset trace gather and
Orthogonal trace gather obtains the instantaneous phase of seismic wave.
The embodiment of the present invention draws east to become by the hyperbolic feature using diffracted wave travel time curve, using improved mobile vertex
Method is changed, accurately obtains the location information of Diffraction Point, and effectively extracted to diffracted wave, effectively suppresses noise and back wave, nothing
Need velocity information that can preferably realize that guarantor's polarity separation of diffracted wave can also be further according to the diffracted wave that extraction obtains
Edge diffraction in geological anomalous body and point diffraction are distinguished.
Diffracted wave separator provided in an embodiment of the present invention has with diffracted wave separation method provided by the above embodiment
Identical technical characteristic reaches identical technical effect so also can solve identical technical problem.
The embodiment of the invention also provides a kind of electronic equipment as shown in figure 8, electronic equipment 8 includes memory 81, processing
Device 82, the computer program that can be run on processor 82 is stored in memory 81, and processor executes real when computer program
The step of method that existing foregoing invention embodiment provides.
Referring to Fig. 8, electronic equipment further include: bus 83 and communication interface 84, processor 82, communication interface 84 and memory
81 are connected by bus 83；Processor 82 is for executing the executable module stored in memory 81, such as computer program.
Wherein, memory 81 may include high-speed random access memory (RAM, Random Access Memory),
It may further include nonvolatile memory (non-volatile memory), for example, at least a magnetic disk storage.By at least
One communication interface 84 (can be wired or wireless) realizes the communication between the system network element and at least one other network element
Connection, can be used internet, wide area network, local network, Metropolitan Area Network (MAN) etc..
Bus 83 can be isa bus, pci bus or eisa bus etc..It is total that bus can be divided into address bus, data
Line, control bus etc..Only to be indicated with a four-headed arrow in Fig. 8, it is not intended that an only bus or one convenient for indicating
The bus of seed type.
Wherein, memory 81 is for storing program, and processor 82 executes program after receiving and executing instruction, and aforementioned
Method performed by invention any embodiment can be applied in processor 82, or be realized by processor 82.
Processor 82 may be a kind of IC chip, the processing capacity with signal.During realization, above-mentioned side
Each step of method can be completed by the integrated logic circuit of the hardware in processor 82 or the instruction of software form.Above-mentioned
Processor 82 can be general processor, including central processing unit (Central Processing Unit, abbreviation CPU), network
Processor (Network Processor, abbreviation NP) etc.；It can also be digital signal processor (Digital Signal
Processing, abbreviation DSP), specific integrated circuit (Application Specific Integrated Circuit, referred to as
ASIC), ready-made programmable gate array (Field-Programmable Gate Array, abbreviation FPGA) or other are programmable
Logical device, discrete gate or transistor logic, discrete hardware components.It may be implemented or execute in the embodiment of the present invention
Disclosed each method, step and logic diagram.General processor can be microprocessor or the processor is also possible to appoint
What conventional processor etc..The step of method in conjunction with disclosed in the embodiment of the present invention, can be embodied directly in hardware decoding processing
Device executes completion, or in decoding processor hardware and software module combination execute completion.Software module can be located at
Machine memory, flash memory, read-only memory, programmable read only memory or electrically erasable programmable memory, register etc. are originally
In the storage medium of field maturation.The storage medium is located at memory 81, and processor 82 reads the information in memory 81, in conjunction with
Its hardware completes the step of above method.
The embodiment of the present invention also provide it is a kind of with processor can be performed non-volatile program code it is computer-readable
Medium, program code make processor execute the method as described in above-mentioned inventive embodiments.
The computer-readable medium of the non-volatile program code provided in an embodiment of the present invention that can be performed with processor,
It is reached with inventive embodiments provided by the above embodiment technical characteristic having the same so also can solve identical technical problem
To identical technical effect.
Computer program product provided by the embodiment of the present invention, including storing the executable non-volatile journey of processor
The computer readable storage medium of sequence code, the instruction that program code includes can be used for executing previous methods as described in the examples
Method, specific implementation can be found in embodiment of the method, and details are not described herein.
It, can be with if the function is realized in the form of SFU software functional unit and when sold or used as an independent product
It is stored in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be embodied in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be a
People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention.
And storage medium above-mentioned includes: that USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic or disk.
Finally, it should be noted that embodiment described above, only a specific embodiment of the invention, to illustrate the present invention
Technical solution, rather than its limitations, scope of protection of the present invention is not limited thereto, although with reference to the foregoing embodiments to this hair
It is bright to be described in detail, those skilled in the art should understand that: anyone skilled in the art
In the technical scope disclosed by the present invention, it can still modify to technical solution documented by previous embodiment or can be light
It is readily conceivable that variation or equivalent replacement of some of the technical features；And these modifications, variation or replacement, do not make
The essence of corresponding technical solution is detached from the spirit and scope of technical solution of the embodiment of the present invention, should all cover in protection of the invention
Within the scope of.Therefore, protection scope of the present invention should be based on the protection scope of the described claims.
Claims (10)
1. a kind of diffracted wave separation method characterized by comprising
Obtain zero-offset trace gather；
The zero-offset trace gather is handled, the instantaneous phase of seismic wave is obtained；The instantaneous phase includes multiple samplings
Point；
Following steps are executed to each sampled point:
The sampled point is converted, the corresponding multiple energy values of the sampled point are obtained；
According to the multiple energy value, the maximum energy value of the sampled point is determined；
The location information of Diffraction Point is determined according to the maximum energy value；
The diffracted wave of the Diffraction Point is extracted according to the positional information, obtains extracting result；
The extraction result of the multiple sampled point is overlapped, the diffracted wave after being separated.
2. being obtained the method according to claim 1, wherein described handle the zero-offset trace gather
The instantaneous phase of seismic wave, comprising:
The zero-offset trace gather is handled using median filtering method, the zero-offset trace gather that obtains that treated；
Hilbert transform is carried out to treated the zero-offset trace gather, obtains treated the zero-offset trace gather
Orthogonal trace gather；
According to treated zero-offset trace gather and the orthogonal trace gather, the instantaneous phase of the seismic wave is obtained.
3. obtaining described adopt the method according to claim 1, wherein described convert the sampled point
The corresponding multiple energy values of sampling point, comprising:
Obtain multiple scanning curvature；
According to multiple scanning curvature, it is utilized respectively improved mobile vertex Radon transform method and the sampled point is turned
It changes, obtains the corresponding energy value of each scanning curvature.
4. according to the method described in claim 3, it is characterized in that, the improved mobile vertex Radon transform method are as follows:
Wherein, x0Indicate the coordinate of the sampled point；T indicates the time of the sampled point；Q indicates the scanning curvature；τ is indicated
The time of the sampled point after conversion；xiIndicate the coordinate of the remaining sampled point in addition to the sampled point；It is adopted described in M expression
Road number before sampling point coordinate；N indicates the road number after the sample point coordinate.
5. the method according to claim 1, wherein the location information includes the Taoist monastic name letter of the Diffraction Point
Breath, temporal information and scanning curvature information.
6. according to the method described in claim 5, it is characterized in that, it is described according to the positional information to the Diffraction Point around
Ejected wave extracts, and obtains extracting result, comprising:
According to the Taoist monastic name information, double branch Radon transforms are carried out using the first calculation formula；First calculation formula are as follows:
Wherein, x0Indicate the coordinate of the corresponding sampled point of the location information；T indicates that the location information is corresponding described
The time of sampled point；xiIndicate the coordinate of the remaining sampled point in addition to the sampled point；M is indicated before the sample point coordinate
Road number；N indicates the road number after the sample point coordinate；
According to the temporal information and the scanning curvature information, the region at (T, Q) is drawn using the second calculation formula
Eastern inverse transformation obtains the extraction result；Second calculation formula are as follows:
Wherein, τ indicates the temporal information T；Q indicates the scanning curvature information Q；x0Indicate the corresponding institute of the location information
State the coordinate of sampled point；M indicates the road number before the sample point coordinate；N indicates the road number after the sample point coordinate.
7. a kind of diffracted wave separator characterized by comprising
Module is obtained, for obtaining zero-offset trace gather；
Processing module obtains the instantaneous phase of seismic wave for handling the zero-offset trace gather；The instantaneous phase
Including multiple sampled points；
Extraction module, for executing following steps to each sampled point:
The sampled point is converted, the corresponding multiple energy values of the sampled point are obtained；
According to the multiple energy value, the maximum energy value of the sampled point is determined；
The location information of Diffraction Point is determined according to the maximum energy value；
The diffracted wave of the Diffraction Point is extracted according to the positional information, obtains extracting result；
Laminating module, for the extraction result of the multiple sampled point to be overlapped, the diffracted wave after being separated.
8. device according to claim 7, which is characterized in that the processing module, comprising:
Filter module, for being handled the zero-offset trace gather using median filtering method, the zero offset that obtains that treated
Away from trace gather；
Conversion module obtains that described treated for carrying out Hilbert transform to treated the zero-offset trace gather
The orthogonal trace gather of zero-offset trace gather；
Computing module, for obtaining the seismic wave according to treated zero-offset trace gather and the orthogonal trace gather
Instantaneous phase.
9. a kind of electronic equipment, including memory, processor, be stored in the memory to run on the processor
Computer program, which is characterized in that the processor realizes that the claims 1 to 6 are any when executing the computer program
The step of method described in item.
10. a kind of computer-readable medium for the non-volatile program code that can be performed with processor, which is characterized in that described
Program code makes the processor execute described any the method for claim 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811358076.6A CN109143344A (en) | 2018-11-14 | 2018-11-14 | Diffracted wave separation method, device and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811358076.6A CN109143344A (en) | 2018-11-14 | 2018-11-14 | Diffracted wave separation method, device and electronic equipment |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109143344A true CN109143344A (en) | 2019-01-04 |
Family
ID=64806037
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811358076.6A Pending CN109143344A (en) | 2018-11-14 | 2018-11-14 | Diffracted wave separation method, device and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109143344A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113474817A (en) * | 2019-03-21 | 2021-10-01 | 徕卡显微系统股份有限公司 | System for off-axis imaging of a sample surface and related method and computer program product |
| CN113640872A (en) * | 2021-08-12 | 2021-11-12 | 中国矿业大学(北京) | Diffracted wave separation method and device and electronic equipment |
| CN113687416A (en) * | 2021-08-19 | 2021-11-23 | 吉林大学 | Imaging method for inclined metal ore body |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012107792A1 (en) * | 2011-02-08 | 2012-08-16 | Total Sa | Method of analyzing seismic data |
| US20170023686A1 (en) * | 2015-07-22 | 2017-01-26 | Pgs Geophysical As | Prediction and Subtraction of Multiple Diffractions |
| US20170248721A1 (en) * | 2014-09-10 | 2017-08-31 | Cgg Services | Wave-field reconstruction using a reflection from a variable sea surface |
| CN107861156A (en) * | 2017-10-30 | 2018-03-30 | 中国矿业大学(北京) | The extracting method and device of diffracted wave |
-
2018
- 2018-11-14 CN CN201811358076.6A patent/CN109143344A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012107792A1 (en) * | 2011-02-08 | 2012-08-16 | Total Sa | Method of analyzing seismic data |
| US20170248721A1 (en) * | 2014-09-10 | 2017-08-31 | Cgg Services | Wave-field reconstruction using a reflection from a variable sea surface |
| US20170023686A1 (en) * | 2015-07-22 | 2017-01-26 | Pgs Geophysical As | Prediction and Subtraction of Multiple Diffractions |
| CN107861156A (en) * | 2017-10-30 | 2018-03-30 | 中国矿业大学(北京) | The extracting method and device of diffracted wave |
Non-Patent Citations (1)
| Title |
|---|
| 李晓峰 等: "绕射波分离成像技术在断裂识别中的应用", 《新疆石油地质》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113474817A (en) * | 2019-03-21 | 2021-10-01 | 徕卡显微系统股份有限公司 | System for off-axis imaging of a sample surface and related method and computer program product |
| CN113640872A (en) * | 2021-08-12 | 2021-11-12 | 中国矿业大学(北京) | Diffracted wave separation method and device and electronic equipment |
| CN113640872B (en) * | 2021-08-12 | 2022-03-08 | 中国矿业大学(北京) | Diffracted wave separation method and device and electronic equipment |
| CN113687416A (en) * | 2021-08-19 | 2021-11-23 | 吉林大学 | Imaging method for inclined metal ore body |
| CN113687416B (en) * | 2021-08-19 | 2022-08-30 | 吉林大学 | Imaging method for inclined metal ore body |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Duncanson et al. | An efficient, multi-layered crown delineation algorithm for mapping individual tree structure across multiple ecosystems | |
| CN105607121B (en) | A kind of coal karst collapse col umn recognition methods and device | |
| CN104932010B (en) | A kind of diffracted wave separation method based on the sparse Radon transformation of shortcut fringing | |
| CN105785439B (en) | The Forecasting Methodology and device of small scale heterogeneous geologic body spatial distribution position | |
| CN109029381A (en) | A kind of detection method of tunnel slot, system and terminal device | |
| CN104793253B (en) | Aviation electromagnetic data de-noising method based on mathematical morphology | |
| CN109143344A (en) | Diffracted wave separation method, device and electronic equipment | |
| CN109164483A (en) | Multi-component earthquake data vector denoising method and multi-component earthquake data vector denoise device | |
| CN104200212A (en) | Building outer boundary line extraction method based on onboard LiDAR (Light Detection and Ranging) data | |
| CN103758511B (en) | Method and device for identifying hidden reservoir through underground reverse time migration imaging | |
| CN103675897A (en) | Seismic diffracted wave separating and imaging method | |
| CN104280784A (en) | Method for recognizing small fracture through gravity | |
| CN107861156B (en) | The extracting method and device of diffracted wave | |
| CN109143343B (en) | Diffraction wave imaging method, device and electronic equipment | |
| CN106772593B (en) | The imaging method and device of diffracted wave | |
| CN103675896B (en) | A kind of diffracted wave and echo method for separate imaging | |
| CN105068131A (en) | Aeromagnetic data leveling method | |
| CN108629760A (en) | A kind of remote sensing image Changing Area Detection method and device | |
| CN111736217B (en) | Seismic attribute fusion method and device | |
| CN102801459B (en) | Multi-mode satellite signal receiving device and method thereof | |
| Ni et al. | Edge guidance network for semantic segmentation of high resolution remote sensing images | |
| Zhao et al. | Automatic extraction of yardangs using Landsat 8 and UAV images: A case study in the Qaidam Basin, China | |
| CN101782964A (en) | Weight loss feather extraction method based on acceleration transducer | |
| CN107749174A (en) | Parking detection method, device, equipment and the storage medium of anti-subway interference | |
| CN103513288B (en) | A kind of compensation direction filtering method of two-dimensional grid data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190104 |
|
| RJ01 | Rejection of invention patent application after publication |