CN105204065B - A kind of method and apparatus of first break picking ripple - Google Patents
A kind of method and apparatus of first break picking ripple Download PDFInfo
- Publication number
- CN105204065B CN105204065B CN201510649394.8A CN201510649394A CN105204065B CN 105204065 B CN105204065 B CN 105204065B CN 201510649394 A CN201510649394 A CN 201510649394A CN 105204065 B CN105204065 B CN 105204065B
- Authority
- CN
- China
- Prior art keywords
- sampled point
- seismic channel
- represent
- point
- envelope
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the present application provides a kind of method and apparatus of first break picking ripple.This method includes carrying out convolution transform to the geological data of collection, obtains the geological data after convolution；The envelope for obtaining sampled point in seismic channel is calculated using the geological data after the geological data and the convolution of the collection；The instantaneous amplitude of sampled point in the seismic channel is determined according to the envelope of the reference sample point of sampled point predeterminated position in the seismic channel；The energy ratio of sampled point in the seismic channel is calculated using the instantaneous amplitude and preset energy ratio formula of sampled point in seismic channel；Calculate the difference of the energy ratio of sampled point and the energy ratio of latter sampled point of the sampled point in seismic channel；Compare the size of the absolute value of the difference；Using the maximum corresponding sampled point of absolute value of numerical value as the preliminary wave of seismic channel sampling point.The technical scheme provided using the embodiment of the present application can greatly improve the accuracy rate of first break picking ripple.
Description
Technical field
The present invention relates to geophysical exploration technical field of data processing, more particularly to a kind of method of first break picking ripple and
Device.
Background technology
Static corrections processing technology in seismic prospecting can change for solving seismic wave near surface medium hourage
The problem of influenceing back wave stacking image.In static corrections processing technology, first arrival is an important attribute of seismic data.On ground
During seismic data processing, the static correction of geological data can be generally solved the problems, such as using the attribute of first arrival.When after the earthquake,
Because the spread speed of various seismic waves is different, traveling to the time of observation station also just has successively.Seismological observation point is received at first
To ripple be just properly termed as preliminary wave.
The method of first break picking ripple mainly includes the method first break picking ripple that boundary energy is detected in the prior art.Specifically
, energy measuring can be carried out by the geological data to collection, regard the corresponding ripple of the high sampled point of energy in trace gather as institute
State the preliminary wave of trace gather.But recently as the extensive use of controlled source, the geological data of collection typically have big data quantity,
The characteristics of low signal-to-noise ratio.This results in the method that is detected in the prior art by boundary energy from big data quantity, low signal-to-noise ratio
In geological data during first break picking ripple, geological data has a large amount of energy interference information, such as later arrivals etc. in itself.So as to reduce
The accuracy rate of first break picking.
Therefore, a kind of method of effective first break picking ripple is needed badly in the prior art, so as to be later use preliminary wave
Attribute solves the problems, such as that the static correction of geological data provides effective foundation.
The content of the invention
The purpose of the application is to provide a kind of method and apparatus of first break picking ripple, to improve the accurate of first break picking ripple
Rate, solves the problems, such as that the static correction of geological data provides effective foundation for the attribute of later use preliminary wave.
To achieve these goals, this application provides a kind of method of first break picking ripple, methods described includes：
Convolution transform is carried out to the geological data of collection, the geological data after convolution is obtained；
Calculated using the geological data after the geological data and the convolution of the collection and obtain sampled point in seismic channel
Envelope；
Determined to sample in the seismic channel according to the envelope of the reference sample point of sampled point predeterminated position in the seismic channel
The instantaneous amplitude of point；
Sampled point in the seismic channel is calculated using the instantaneous amplitude and preset energy ratio formula of sampled point in seismic channel
Energy ratio；
Calculate the difference of the energy ratio of sampled point and the energy ratio of latter sampled point of the sampled point in seismic channel；
Compare the size of the absolute value of the difference；
Using the maximum corresponding sampled point of absolute value of numerical value as the preliminary wave of the seismic channel sampling point；
Wherein, the preset energy ratio formula is as follows：
In above formula, RgRepresent the energy ratio of g-th of sampled point in the seismic channel；1≤g≤n, g >=W, n represent described
The sampled point quantity of seismic channel；HrRepresent the instantaneous amplitude of r sampled points in the seismic channel；HgRepresent the wink of the first sampled point
When amplitude；Window sampled point quantity when W is represented, is positive integer.
In a preferred embodiment, the geological data after the geological data and the convolution using the collection
The calculation formula that calculating obtains the envelope of sampled point in seismic channel is as follows：
In above formula, ErRepresent the envelope of r-th of sampled point in the seismic channel；1≤r≤n, n represent the seismic channel
Sampled point quantity；xrRepresent the geological data of r-th of sampled point in the seismic channel；R-th is represented in the seismic channel to adopt
Geological data after the convolution of sampling point.
In a preferred embodiment, it is described according to the reference sample point of sampled point predeterminated position in the seismic channel
Envelope determines that the instantaneous amplitude of sampled point in the seismic channel includes：
The when window sampled point quantity of the seismic channel is obtained, window sampled point quantity is used as predetermined number when will be described；
Calculated using the envelope of the sampled point of the predetermined number before the first sampled point in the seismic channel and obtain described first
The energy scale variance of the average energy value of sampled point and first sampled point；
Calculate using the average energy value of first sampled point and the energy scale variance of first sampled point and obtain institute
State the instantaneous amplitude of sampled point in seismic channel.
In a preferred embodiment, the sampling of the predetermined number using before the first sampled point in the seismic channel
The envelope of point calculates the calculating for obtaining the average energy value of first sampled point and the energy scale variance of first sampled point
Formula is as follows：
In above formula, AgRepresent the average energy value of first sampled point；1≤g≤n, g >=W, g represent first sampling
The period of point, n represents the sampled point quantity of the seismic channel；ErRepresent the envelope of r-th of sampled point in the seismic channel；W generations
Window sampled point quantity during table, is positive integer；
In above formula, DgRepresent the energy scale variance of first sampled point；ErRepresent r-th of sampling in the seismic channel
The envelope of point；Window sampled point quantity when W is represented, is positive integer；AgRepresent the average energy value of first sampled point.
In a preferred embodiment, the average energy value and first sampled point using first sampled point
Energy scale variance calculate obtain the instantaneous amplitude of sampled point in the seismic channel calculation formula it is as follows：
In above formula, HgRepresent the instantaneous amplitude of first sampled point；AgRepresent the average energy value of first sampled point；
DgRepresent the energy scale variance of first sampled point.
A kind of device of first break picking ripple, described device includes：
Convolution transform module, carries out convolution transform for the geological data to collection, obtains the geological data after convolution；
First computing module, calculates for the geological data after the geological data using the collection and the convolution and obtains
The envelope of sampled point in seismic channel；
Determining module, described in being determined according to the envelope of the reference sample point of sampled point predeterminated position in the seismic channel
The instantaneous amplitude of sampled point in seismic channel；
Second computing module, for calculating institute using the instantaneous amplitude and preset energy ratio formula of sampled point in seismic channel
State the energy ratio of sampled point in seismic channel；
3rd computing module, for calculating the energy ratio of sampled point and latter sampled point of the sampled point in seismic channel
Energy ratio difference；
Comparison module, the size of the absolute value for comparing the difference；
Data processing module, for using the maximum corresponding sampled point of absolute value of numerical value as the seismic channel preliminary wave
Sampling point；
Wherein, the preset energy ratio formula is as follows：
In above formula, RgRepresent the energy ratio of g-th of sampled point in the seismic channel；1≤g≤n, g >=W, n represent described
The sampled point quantity of seismic channel；HrRepresent the instantaneous amplitude of r sampled points in the seismic channel；HgRepresent the wink of the first sampled point
When amplitude；Window sampled point quantity when W is represented, is positive integer.
In a preferred embodiment, the calculation formula of first computing module is as follows：
In above formula, ErRepresent the envelope of r-th of sampled point in the seismic channel；1≤r≤n, n represent the seismic channel
Sampled point quantity；xrRepresent the geological data of r-th of sampled point in the seismic channel；Represent r-th of sampling in the seismic channel
Geological data after the convolution of point.
In a preferred embodiment, the determining module includes：
Data capture unit, the when window sampled point quantity for obtaining the seismic channel, window sampled point quantity when will be described
It is used as predetermined number；
First computing unit, the envelope for the sampled point using the predetermined number before the first sampled point in the seismic channel
Calculating obtains the average energy value of first sampled point and the energy scale variance of first sampled point；
Second computing unit, for the average energy value and the energy mark of first sampled point using first sampled point
Quasi- variance calculates the instantaneous amplitude for obtaining sampled point in the seismic channel.
In a preferred embodiment, the calculation formula of first computing unit can be as follows：
In above formula, AgRepresent the average energy value of first sampled point；1≤g≤n, g >=W, g represent first sampling
The period of point, n represents the sampled point quantity of the seismic channel；ErRepresent the envelope of r-th of sampled point in the seismic channel；W generations
Window sampled point quantity during table, is positive integer；
In above formula, DgRepresent the energy scale variance of first sampled point；ErRepresent r-th of sampling in the seismic channel
The envelope of point；Window sampled point quantity when W is represented, is positive integer；AgRepresent the average energy value of first sampled point.
In a preferred embodiment, the calculation formula of second computing unit is as follows：
In above formula, HgRepresent the instantaneous amplitude of first sampled point；AgRepresent the average energy value of first sampled point；
DgRepresent the energy scale variance of first sampled point.
The technical scheme provided from above the embodiment of the present application, the embodiment of the present application using collection geological data and
Geological data after convolution can calculate the envelope for obtaining sampled point in seismic channel；Then, according to sampled point in the seismic channel
The envelope of the reference sample point of predeterminated position determines the instantaneous amplitude of sampled point in the seismic channel, and the instantaneous amplitude can disappear
Random disturbances before being reached except preliminary wave.And the instantaneous amplitude using sampled point in the seismic channel and later arrivals pair can be recognized
The preset energy ratio formula of the interference of preliminary wave calculates the energy ratio of sampled point in the seismic channel；Calculate the earthquake
The difference of the energy ratio of sampled point and the energy ratio of latter sampled point of the sampled point in road；Compare the exhausted of the difference
To the size of value；Using the maximum corresponding sampled point of absolute value of numerical value as the preliminary wave of the seismic channel sampling point.Utilize this
The technical scheme that application embodiment is provided can effectively recognize interference of the later arrivals to preliminary wave, accurate first break picking ripple.With showing
There is technology to compare, the accuracy rate of first break picking ripple can be greatly improved, be that the attribute of later use preliminary wave solves geological data
Static correction problem effective foundation is provided.
Brief description of the drawings
, below will be to embodiment or existing in order to illustrate more clearly of the embodiment of the present application or technical scheme of the prior art
There is the accompanying drawing used required in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments described in application, for those of ordinary skill in the art, are not paying the premise of creative labor
Under, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 is a kind of flow chart of the embodiment of the method for first break picking ripple of the application；
Fig. 2 is the flow chart of the instantaneous amplitude process of sampled point in the determination seismic channel that the embodiment of the present application is provided；
Fig. 3 is a kind of schematic diagram of the embodiment of the device of first break picking ripple of the application；
Fig. 4 is the schematic diagram for the determining module that the embodiment of the present application is provided.
Embodiment
In order that those skilled in the art more fully understand the technical scheme in the application, it is real below in conjunction with the application
The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described implementation
Example only some embodiments of the present application, rather than whole embodiments.Based on the embodiment in the application, this area is common
The every other embodiment that technical staff is obtained under the premise of creative work is not made, should all belong to the application protection
Scope.
Below implementing for the embodiment of the present application is described in detail with several specific examples.
Introduce a kind of first embodiment of method of the application first below.Fig. 1 is a kind of side of first break picking ripple of the application
The flow chart of the embodiment of method, although multiple operations that flow includes occurring with particular order are described below, but should understand
Understand, these processes can include more or less operations, these operations can sequentially be performed or performed parallel.With reference to accompanying drawing
1, the embodiment includes：
S110：Convolution transform is carried out to the geological data of collection, the geological data after convolution is obtained.
In actual applications, the geological data of collection can be rolled up after exciting and gathering geological data
Product transformation, obtains the geological data after convolution.Specifically, the convolution transform can include Hilbert transform (Hilbert
transform)。
In a specific embodiment, the geological data of the collection can include single-shot geological data.Specifically, institute
24 arrangements can be included by stating single-shot geological data, and trace gather is 7680, and the sampling interval is 2ms, and per pass sampling number is 2000
Geological data.
In addition, convolution transform described in the embodiment of the present application is not limited in above-mentioned Hilbert transform, actually should
In, other method can also be included, the embodiment of the present application is not limited thereto.
S120：Calculated using the geological data after the geological data and the convolution of the collection and obtain sampling in seismic channel
The envelope of point.
In certain embodiments, after step silo, it is possible to use after the geological data of the collection and the convolution
Geological data calculate and obtain the envelope of sampled point in seismic channel.Specifically, can be calculated using formula is calculated as below
The envelope of sampled point into the seismic channel：
In above formula, ErRepresent the envelope of r-th of sampled point in the seismic channel；1≤r≤n, n represent the seismic channel
Sampled point quantity；xrRepresent the geological data of r-th of sampled point in the seismic channel；Represent r-th of sampling in the seismic channel
Geological data after the convolution of point.
S130：Determined according to the envelope of the reference sample point of sampled point predeterminated position in the seismic channel in the seismic channel
The instantaneous amplitude of sampled point.
In certain embodiments, can be according to the ginseng of sampled point predeterminated position in the seismic channel after step S120
The envelope for examining sampled point determines the instantaneous amplitude of sampled point in the seismic channel.Fig. 2 is that the embodiment of the present application is provided definitely
Shake the flow chart of the instantaneous amplitude process of sampled point in road.Specifically, with reference to accompanying drawing 2, can include：
S131：The when window sampled point quantity of the seismic channel is obtained, window sampled point quantity is used as predetermined number when will be described.
In certain embodiments, the when window sampled point quantity of the seismic channel, window sampling number when will be described can be obtained
Amount is used as predetermined number.
S132：Calculated using the envelope of the sampled point of the predetermined number before the first sampled point in the seismic channel and obtain described
The energy scale variance of the average energy value of first sampled point and first sampled point.
In certain embodiments, after step S131, it is possible to use default before the first sampled point in the seismic channel
The envelope of the sampled point of quantity calculates the energy scale of the average energy value for obtaining first sampled point and first sampled point
Variance.And when the sampling number before first sampled point is less than the predetermined number, adopted using in the seismic channel first
The envelope of sampled point before sampling point calculates the energy mark of the average energy value for obtaining first sampled point and first sampled point
Quasi- variance.Specifically, first sampled point can include any one sampled point in the seismic channel.Specifically, can wrap
Include to utilize and the average energy value that formula calculating obtains first sampled point is calculated as below：
In above formula, AgRepresent the average energy value of first sampled point；1≤g≤n, g >=W, g represent first sampling
The period of point, n represents the sampled point quantity of the seismic channel；ErRepresent the envelope of r-th of sampled point in the seismic channel；W generations
Window sampled point quantity during table, is positive integer, in a specific embodiment, and window sampled point quantity can be 20 when described.
Specifically, the energy scale variance of first sampled point can be obtained using formula being calculated as below calculating：
In above formula, DgRepresent the energy scale variance of first sampled point；ErRepresent r-th of sampling in the seismic channel
The envelope of point；Window sampled point quantity when W is represented, is positive integer, in a specific embodiment, window sampled point quantity when described
Can be 20；AgRepresent the average energy value of first sampled point.
S133：Calculated using the average energy value of first sampled point and the energy scale variance of first sampled point
The instantaneous amplitude of sampled point into the seismic channel.
In certain embodiments, after step S132, it is possible to use the average energy value of first sampled point and described
The energy scale variance of first sampled point calculates the instantaneous amplitude for obtaining sampled point in the seismic channel.Specifically, can include
The instantaneous amplitude of first sampled point is obtained using formula being calculated as below calculating：
In above formula, HgRepresent the instantaneous amplitude of first sampled point；AgRepresent the average energy value of first sampled point；
DgRepresent the energy scale variance of first sampled point.
S140：Calculated and adopted in the seismic channel using the instantaneous amplitude and preset energy ratio formula of sampled point in seismic channel
The energy ratio of sampling point.
In certain embodiments, after step s 130, it is possible to use the instantaneous amplitude of sampled point and default in seismic channel
Energy ratio formula calculates the energy ratio of sampled point in the seismic channel.Specifically, the preset energy ratio formula is as follows：
In above formula, RgRepresent the energy ratio of g-th of sampled point in the seismic channel；1≤g≤n, g >=W, g represent described
The period of first sampled point, n represents the sampled point quantity of the seismic channel；HrRepresent the wink of r sampled points in the seismic channel
When amplitude；HgRepresent the instantaneous amplitude of the first sampled point；Window sampled point quantity when W is represented, is positive integer, specific real at one
Apply in example, window sampled point quantity can be 20 when described.
S150：Calculate the energy ratio of sampled point and the energy ratio of latter sampled point of the sampled point in seismic channel
Difference.
In certain embodiments, after step s 140, can calculate in seismic channel the energy ratio of sampled point with it is described
The difference of the energy ratio of latter sampled point of sampled point.Specifically, it is possible to use formula is calculated as below：
Ui=Ri-Ri+1
In above formula, UiRepresent the energy ratio of ith sample point and the ith sample point in the seismic channel latter individual
The difference of the energy ratio of sampled point (i+1 sampled point)；RiRepresent the energy ratio of the ith sample point；Ri+1Represent
The energy ratio of the i+1 sampled point.
S160：Compare the size of the absolute value of the difference.
In certain embodiments, after the step s 150, the size of the absolute value of the difference can be compared.
S170：Using the maximum corresponding sampled point of absolute value of numerical value as the preliminary wave of the seismic channel sampling point.
In certain embodiments, the corresponding sampled point of the maximum absolute value of numerical value in step S160 comparative results can be made
For the sampling point of the preliminary wave of the seismic channel.
As can be seen here, the technical scheme that a kind of embodiment of the method for first break picking ripple of the application is provided utilizes the earthquake gathered
Geological data after data and convolution can calculate the envelope for obtaining sampled point in seismic channel；Then, according in the seismic channel
The envelope of the reference sample point of sampled point predeterminated position determines the instantaneous amplitude of sampled point in the seismic channel, the instantaneous amplitude
The random disturbances before preliminary wave is reached can be eliminated.And the instantaneous amplitude using sampled point in the seismic channel and can recognize it is continuous
The energy ratio of in seismic channel sampled point is calculated the preset energy ratio formula of the interference of preliminary wave to ripple；Calculate institute
State the difference of the energy ratio of sampled point and the energy ratio of latter sampled point of the sampled point in seismic channel；Compare the difference
The size of the absolute value of value；Using the maximum corresponding sampled point of absolute value of numerical value as the preliminary wave of the seismic channel sampling point.
The technical scheme provided using the embodiment of the present application can effectively recognize interference of the later arrivals to preliminary wave, accurate first break picking
Ripple.Compared with prior art, the accuracy rate of first break picking ripple can be greatly improved, is that the attribute of later use preliminary wave solves ground
The static correction problem for shaking data provides effective foundation.
On the other hand the application also provides a kind of device of first break picking ripple, and Fig. 3 is a kind of first break picking ripple of the application
The schematic diagram of the embodiment of device, with reference to accompanying drawing 3, described device 300 can include：
Convolution transform module 310, can be used for carrying out convolution transform to the geological data of collection, obtains the earthquake after convolution
Data；
First computing module 320, the geological data that can be used for after the geological data and the convolution using the collection
Calculate the envelope for obtaining sampled point in seismic channel；
Determining module 330, can be used for the envelope according to the reference sample point of sampled point predeterminated position in the seismic channel
Determine the instantaneous amplitude of sampled point in the seismic channel；
Second computing module 340, can be used for utilizing the instantaneous amplitude and preset energy ratio of sampled point in seismic channel public
Formula calculates the energy ratio of sampled point in the seismic channel；
3rd computing module 350, the energy ratio and the sampled point that can be used for calculating sampled point in seismic channel is latter
The difference of the energy ratio of individual sampled point；
Comparison module 360, can be used for the size of the absolute value of difference described in comparison；
Data processing module 370, can be used for regarding the maximum corresponding sampled point of absolute value of numerical value as the seismic channel
Preliminary wave sampling point.
In a preferred embodiment, the calculation formula of first computing module can be as follows：
In above formula, ErRepresent the envelope of r-th of sampled point in the seismic channel；1≤r≤n, n represent the seismic channel
Sampled point quantity；xrRepresent the geological data of r-th of sampled point in the seismic channel；Represent r-th of sampling in the seismic channel
Geological data after the convolution of point.
Fig. 4 is the schematic diagram for the determining module that the embodiment of the present application is provided.In a preferred embodiment, with reference to accompanying drawing
4, the determining module 330 can include：
Data capture unit 331, can be used for the when window sampled point quantity for obtaining the seismic channel, and window is sampled when will be described
Point quantity is used as predetermined number；
First computing unit 332, can be used for the sampling for utilizing the predetermined number in the seismic channel before the first sampled point
The envelope of point calculates and obtains the average energy value of first sampled point and the energy scale variance of first sampled point；
Second computing unit 333, can be used for the average energy value and first sampled point using first sampled point
Energy scale variance calculate and obtain the instantaneous amplitude of sampled point in the seismic channel.
In a preferred embodiment, the calculation formula of first computing unit can be as follows：
In above formula, AgRepresent the average energy value of first sampled point；1≤g≤n, g >=W, g represent first sampling
The period of point, n represents the sampled point quantity of the seismic channel；ErRepresent the envelope of r-th of sampled point in the seismic channel；W generations
Window sampled point quantity during table, is positive integer, in a specific embodiment, and window sampled point quantity can be 20 when described；
In above formula, DgRepresent the energy scale variance of first sampled point；ErRepresent r-th of sampling in the seismic channel
The envelope of point；Window sampled point quantity when W is represented, is positive integer, in a specific embodiment, window sampled point quantity when described
Can be 20；AgRepresent the average energy value of first sampled point.
In a preferred embodiment, the calculation formula of second computing unit can be as follows：
In above formula, HgRepresent the instantaneous amplitude of first sampled point；AgRepresent the average energy value of first sampled point；
DgRepresent the energy scale variance of first sampled point.
In a preferred embodiment, the preset energy ratio formula can be as follows：
In above formula, RgRepresent the energy ratio of g-th of sampled point in the seismic channel；1≤g≤n, g >=W, g represent described
The period of first sampled point, n represents the sampled point quantity of the seismic channel；HrRepresent the wink of r sampled points in the seismic channel
When amplitude；HgRepresent the instantaneous amplitude of the first sampled point；Window sampled point quantity when W is represented, is positive integer, specific real at one
Apply in example, window sampled point quantity can be 20 when described.
As can be seen here, the technical scheme that a kind of method and apparatus embodiment of first break picking ripple of the application is provided utilizes collection
Geological data and convolution after geological data can calculate the envelope for obtaining sampled point in seismic channel；Then, according to described
The envelope of the reference sample point of sampled point predeterminated position determines the instantaneous amplitude of sampled point in the seismic channel, the wink in shake road
When amplitude can eliminate preliminary wave reach before random disturbances.And instantaneous amplitude using sampled point in the seismic channel and can be with
Recognize that later arrivals calculate the preset energy ratio formula of the interference of preliminary wave the energy ratio of sampled point in the seismic channel；Meter
Calculate the difference of the energy ratio of sampled point and the energy ratio of latter sampled point of the sampled point in the seismic channel；Compare
The size of the absolute value of the difference；It regard the maximum corresponding sampled point of absolute value of numerical value as the preliminary wave of the seismic channel
Sampling point.The technical scheme provided using the embodiment of the present application can effectively recognize interference of the later arrivals to preliminary wave, accurate pickup
Preliminary wave.Compared with prior art, the accuracy rate of first break picking ripple can be greatly improved, is the attribute solution of later use preliminary wave
Certainly the static correction problem of geological data provides effective foundation.
Each embodiment in this specification is described by the way of progressive, what each embodiment was stressed be with
Between the difference of other embodiment, each embodiment identical similar part mutually referring to.It is real especially for system
Apply for example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method
Part explanation.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and
Change is without departing from spirit herein, it is desirable to which appended claim includes these deformations and changed without departing from the application's
Spirit.
Claims (10)
1. a kind of method of first break picking ripple, it is characterised in that methods described includes：
Convolution transform is carried out to the geological data of collection, the geological data after convolution is obtained；
The envelope for obtaining sampled point in seismic channel is calculated using the geological data after the geological data and the convolution of the collection；
Sampled point in the seismic channel is determined according to the envelope of the reference sample point of sampled point predeterminated position in the seismic channel
Instantaneous amplitude；
The energy of sampled point in the seismic channel is calculated using the instantaneous amplitude and preset energy ratio formula of sampled point in seismic channel
Measure ratio；
Calculate the difference of the energy ratio of sampled point and the energy ratio of latter sampled point of the sampled point in seismic channel；
Compare the size of the absolute value of the difference；
Using the maximum corresponding sampled point of absolute value of numerical value as the preliminary wave of the seismic channel sampling point；
Wherein, the preset energy ratio formula is as follows：
In above formula, RgRepresent the energy ratio of g-th of sampled point in the seismic channel；1≤g≤n, g >=W, n represent the earthquake
The sampled point quantity in road；HrRepresent the instantaneous amplitude of r sampled points in the seismic channel；HgThe instantaneous of the first sampled point is represented to shake
Width；Window sampled point quantity when W is represented, is positive integer.
2. according to the method described in claim 1, it is characterised in that the geological data and the convolution using the collection
The calculation formula that geological data calculating afterwards obtains the envelope of sampled point in seismic channel is as follows：
In above formula, ErRepresent the envelope of r-th of sampled point in the seismic channel；1≤r≤n, n represent the sampled point of the seismic channel
Quantity；xrRepresent the geological data of r-th of sampled point in the seismic channel；Represent in the seismic channel r-th sampled point
Geological data after convolution.
3. according to the method described in claim 1, it is characterised in that described according to sampled point predeterminated position in the seismic channel
The envelope of reference sample point determines that the instantaneous amplitude of sampled point in the seismic channel includes：
The when window sampled point quantity of the seismic channel is obtained, window sampled point quantity is used as predetermined number when will be described；
Calculated using the envelope of the sampled point of the predetermined number before the first sampled point in the seismic channel and obtain first sampling
The average energy value of point and the energy scale variance of first sampled point；
Calculated and obtained describedly using the average energy value of first sampled point and the energy scale variance of first sampled point
Shake the instantaneous amplitude of sampled point in road.
4. method according to claim 3, it is characterised in that described using pre- before the first sampled point in the seismic channel
If the envelope of the sampled point of quantity calculates the energy mark of the average energy value for obtaining first sampled point and first sampled point
The calculation formula of quasi- variance is as follows：
In above formula, AgRepresent the average energy value of first sampled point；1≤g≤n, g >=W, g represent the point of first sampled point
Number, n represents the sampled point quantity of the seismic channel；ErRepresent the envelope of r-th of sampled point in the seismic channel；Window when W is represented
Sampled point quantity, is positive integer；
In above formula, DgRepresent the energy scale variance of first sampled point；ErRepresent in the seismic channel r-th sampled point
Envelope；Window sampled point quantity when W is represented, is positive integer；AgRepresent the average energy value of first sampled point.
5. method according to claim 3, it is characterised in that the average energy value and institute using first sampled point
State the first sampled point energy scale variance calculate obtain the instantaneous amplitude of sampled point in the seismic channel calculation formula it is as follows：
In above formula, HgRepresent the instantaneous amplitude of first sampled point；AgRepresent the average energy value of first sampled point；DgGeneration
The energy scale variance of first sampled point described in table.
6. a kind of device of first break picking ripple, it is characterised in that described device includes：
Convolution transform module, carries out convolution transform for the geological data to collection, obtains the geological data after convolution；
First computing module, calculates for the geological data after the geological data using the collection and the convolution and obtains earthquake
The envelope of sampled point in road；
Determining module, for determining the earthquake according to the envelope of the reference sample point of sampled point predeterminated position in the seismic channel
The instantaneous amplitude of sampled point in road；
Second computing module, for the instantaneous amplitude using sampled point in seismic channel and the calculating of preset energy ratio formula describedly
Shake the energy ratio of sampled point in road；
3rd computing module, for calculating the energy ratio of sampled point and the energy of latter sampled point of the sampled point in seismic channel
Measure the difference of ratio；
Comparison module, the size of the absolute value for comparing the difference；
Data processing module, for using the maximum corresponding sampled point of absolute value of numerical value as the preliminary wave of the seismic channel sample
Point；
Wherein, the preset energy ratio formula is as follows：
In above formula, RgRepresent the energy ratio of g-th of sampled point in the seismic channel；1≤g≤n, g >=W, n represent the earthquake
The sampled point quantity in road；HrRepresent the instantaneous amplitude of r sampled points in the seismic channel；HgThe instantaneous of the first sampled point is represented to shake
Width；Window sampled point quantity when W is represented, is positive integer.
7. device according to claim 6, it is characterised in that the calculation formula of first computing module is as follows：
In above formula, ErRepresent the envelope of r-th of sampled point in the seismic channel；1≤r≤n, n represent the sampled point of the seismic channel
Quantity；xrRepresent the geological data of r-th of sampled point in the seismic channel；Represent the volume of r-th of sampled point in the seismic channel
Geological data after product.
8. device according to claim 6, it is characterised in that the determining module includes：
Data capture unit, the when window sampled point quantity for obtaining the seismic channel, when will be described window sampled point quantity as
Predetermined number；
First computing unit, the envelope for the sampled point using the predetermined number before the first sampled point in the seismic channel is calculated
Obtain the average energy value of first sampled point and the energy scale variance of first sampled point；
Second computing unit, for the average energy value using first sampled point and the energy scale side of first sampled point
Difference calculates the instantaneous amplitude for obtaining sampled point in the seismic channel.
9. device according to claim 8, it is characterised in that the calculation formula of first computing unit is as follows：
In above formula, AgRepresent the average energy value of first sampled point；1≤g≤n, g >=W, g represent the point of first sampled point
Number, n represents the sampled point quantity of the seismic channel；ErRepresent the envelope of r-th of sampled point in the seismic channel；Window when W is represented
Sampled point quantity, is positive integer；
In above formula, DgRepresent the energy scale variance of first sampled point；ErRepresent in the seismic channel r-th sampled point
Envelope；Window sampled point quantity when W is represented, is positive integer；AgRepresent the average energy value of first sampled point.
10. device according to claim 8, it is characterised in that the calculation formula of second computing unit is as follows：
In above formula, HgRepresent the instantaneous amplitude of first sampled point；AgRepresent the average energy value of first sampled point；DgGeneration
The energy scale variance of first sampled point described in table.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510649394.8A CN105204065B (en) | 2015-10-09 | 2015-10-09 | A kind of method and apparatus of first break picking ripple |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510649394.8A CN105204065B (en) | 2015-10-09 | 2015-10-09 | A kind of method and apparatus of first break picking ripple |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105204065A CN105204065A (en) | 2015-12-30 |
| CN105204065B true CN105204065B (en) | 2017-10-17 |
Family
ID=54951840
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510649394.8A Active CN105204065B (en) | 2015-10-09 | 2015-10-09 | A kind of method and apparatus of first break picking ripple |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105204065B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108072896B (en) * | 2016-11-18 | 2019-08-27 | 中国石油化工股份有限公司 | A kind of fully automatically seismic wave first break pickup method and system |
| TWI647475B (en) * | 2017-03-13 | 2019-01-11 | 顏伯聰 | An earthquake prediction method and system |
| CN109254323A (en) * | 2017-07-14 | 2019-01-22 | 中国石油化工股份有限公司 | Seismic first breaks pick-up method and computer readable storage medium based on energy envelope |
| CN109757925B (en) * | 2018-12-27 | 2021-03-02 | 艾感科技(广东)有限公司 | Mattress based on strike signal early warning |
| CN112180433B (en) * | 2019-07-04 | 2023-08-22 | 中国石油天然气集团有限公司 | Method and device for picking up first arrival wave of earthquake |
| CN112305593B (en) * | 2019-07-23 | 2023-07-25 | 中国石油天然气集团有限公司 | Abnormal first arrival wave identification method and device |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101609160B (en) * | 2008-06-19 | 2011-08-03 | 中国石油天然气集团公司 | Automatic seismic data first-arrival wave pickup method with strong background noise resistance |
| CN101630015B (en) * | 2008-07-16 | 2012-01-11 | 中国石油天然气集团公司 | Method for improving pickup precision and efficiency of first-arrival wave |
| US9103935B2 (en) * | 2013-02-04 | 2015-08-11 | King Fahd University Of Petroleum And Minerals | Method of first arrival picking of seismic refraction data |
-
2015
- 2015-10-09 CN CN201510649394.8A patent/CN105204065B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105204065A (en) | 2015-12-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105204065B (en) | A kind of method and apparatus of first break picking ripple | |
| US10705238B2 (en) | Method and apparatus for processing seismic data | |
| CN105785439B (en) | The Forecasting Methodology and device of small scale heterogeneous geologic body spatial distribution position | |
| CN101609160B (en) | Automatic seismic data first-arrival wave pickup method with strong background noise resistance | |
| CN104570076B (en) | Automatic seismic wave first-arrival picking method based on dichotomy | |
| CN110907995B (en) | Reverse time migration method and device for VSP seismic data in well | |
| WO2009105964A1 (en) | Method of pre-stack two-dimension-like transformation of three-dimensional seismic record | |
| CN103210323B (en) | Process the system and method for geological data | |
| CN109765615A (en) | Stratum quality factor inversion method and device | |
| CN106680869A (en) | Microearthquake event detecting and positioning method and apparatus | |
| CN106646596A (en) | Undulatory sea surface ghost reflection suppression method and device | |
| CN105093318B (en) | A kind of adaptive wave equation wave field extrapolation static correcting method | |
| CN106324671A (en) | Method and device for checking shot point offset | |
| CN102338885B (en) | Three-component VSP data first arrival time automatic pick method | |
| CN105388512B (en) | A kind of calibration method and device of earthquake data acquisition exploration time | |
| CN104977602B (en) | A kind of control method and device of earthquake data acquisition construction | |
| CN106199705B (en) | The detection method of seismic acquisition data | |
| CN109782355A (en) | The detection method and device of OBS detection point drift | |
| CN106468784A (en) | A kind of denoising method based on poststack seismic data signal to noise ratio attribute and device | |
| CN107942391A (en) | A kind of seabed geophone station localization method and device | |
| CN111123361B (en) | Method and device for regularly reconstructing seismic data of vertical seismic profile and storage medium | |
| CN105445788A (en) | Automatic velocity spectrum interpretation method based on model and global optimization | |
| CN110824555B (en) | Seismic energy equalization method and device and computer readable storage medium | |
| CN105510986B (en) | Method and device for determining content of glauconite in glauconite sandstone | |
| CN110161560A (en) | A kind of detection method and device of microseismic event |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |