CN106324671A - Method and device of examining shot point offset - Google Patents
Method and device of examining shot point offset Download PDFInfo
- Publication number
- CN106324671A CN106324671A CN201510379271.7A CN201510379271A CN106324671A CN 106324671 A CN106324671 A CN 106324671A CN 201510379271 A CN201510379271 A CN 201510379271A CN 106324671 A CN106324671 A CN 106324671A
- Authority
- CN
- China
- Prior art keywords
- arrival
- actual
- geological data
- seismic channel
- big gun
- 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
Abstract
The invention provides a method and device of examining shot point offset. The method comprises steps of: reading single shot earthquake data, wherein the single shot earthquake data comprises preset quantity of earthquake channel data; collecting the actual first arrival of each earthquake channel data in the single shot earthquake data; calculating the theoretical first arrival of each earthquake channel data in the single shot earthquake data; and comparing the actual first arrival and the theoretical first arrival corresponding to each earthquake channel data in the single shot earthquake data, and according to a comparison result, determining whether shot point offset exists in the single shot earthquake data. According to the method and the device provided by the embodiment of the invention, the shot point offset can be quickly and precisely examined according to actual earth surface conditions.
Description
Technical field
The application relates to seismic data acquisition technical field, particularly to a kind of method and device checking shot-point offset.
Background technology
In recent years, along with seismic prospecting ability step up with near surface geophysics technology engineering field study constantly deep
Enter so that seismic data acquisition is carried out under complex near surface conditions.During the field acquisition under complex near surface conditions,
While improving field data acquisition technical method and strengthening field construction quality management, in addition it is also necessary to a set of effective technology monitors
Means ensure the quality gathered, and wherein the monitoring of shot position is exactly an important step.
During seismic data acquisition, in-situ processing personnel must carry out effective monitoring to the geological data collected in time,
To find shot point and the drift condition of geophone station, thus reject blown-out shot.During complicated surface conditions often makes earthquake-capturing
Increasing the monitoring difficulty of shot-point offset, particularly in area with a varied topography, that near-surface velocity cross directional variations is big, this problem is outstanding
It is prominent.Therefore, the same day is blown out gained geological data in time carry out check and preliminary treatment, implement quality in site monitoring
Particularly important.Data in-situ processing is directly closed as effective link of earthquake data acquisition quality monitoring, the good and bad of its technological means
It is tied to the height of seismic acquisition quality.
Conventional linear NMO is to check for shot-point offset by sorting section geophone offset.The shortcoming of this method
Being: run into the area that surface infrastructure is complicated, top layer lateral variation in velocity is violent, actual application effect is difficult to meet needs.Separately
Outward, in mass seismic data checking process, each big gun is checked by Quality Control personnel, there will be careless omission unavoidably.Quality Control simultaneously
The geological data amplification that each arrangement can only be received by personnel determines side-play amount, and this is likely to result in the distortion of data.Therefore, root
Earth's surface, border situation factually, finds a kind of suitable, method of accurate quick and seems extremely important.
It should be noted that introduction to technical background above is intended merely to the convenient technical scheme to the application and carries out clear, complete
Illustrate, and facilitate the understanding of those skilled in the art to illustrate.Can not be merely because these schemes be in the background technology of the application
Part is set forth and thinks that technique scheme is known to those skilled in the art.
Summary of the invention
The purpose of the embodiment of the present application is to provide a kind of method and device checking shot-point offset, it is possible to according to condition of expressing one's feelings practically
Check shot-point offset rapidly and accurately.
A kind of method and device checking shot-point offset that the embodiment of the present application provides is achieved in that
A kind of method checking shot-point offset, described method includes:
Reading single big gun geological data, wherein, described single big gun geological data includes the geological data of multiple seismic channel；
Pick up the actual first arrival of each seismic channel in described single big gun geological data；
Calculate the theoretical first arrival of each seismic channel in described single big gun geological data；
The actual first arrival of each seismic channel in described single big gun geological data and theoretical first arrival are contrasted, and according to the result of contrast
Judge whether described single big gun geological data exists shot-point offset.
A kind of device checking shot-point offset, described device includes:
Geological data reads unit, is used for reading single big gun geological data, and wherein, described single big gun geological data includes multiple earthquake
The geological data in road；
Actual first break pickup unit, is used for picking up the actual first arrival of each seismic channel in described single big gun geological data；
Theoretical first arrival computing unit, is used for calculating the theoretical first arrival of each seismic channel in described single big gun geological data；
Shot-point offset identifying unit, is used for carrying out the actual first arrival of each seismic channel in described single big gun geological data and theoretical first arrival
Contrast, and judge whether described single big gun geological data exists shot-point offset according to the result of contrast.
A kind of method and device checking shot-point offset that the embodiment of the present application provides, the actual first arrival picked up by contrast and calculating
Theoretical first arrival judge whether single big gun geological data exists shot-point offset, and the actual first arrival of pickup can be modified place
Reason, can apply to the area that surface infrastructure is complicated, top layer lateral variation in velocity is violent, and all processes all can be automatically real
Existing, it is not necessary to Quality Control personnel check, it is possible to checks shot-point offset rapidly and accurately.
Accompanying drawing explanation
A kind of method flow diagram checking shot-point offset that Fig. 1 provides for the embodiment of the present application；
The functional block diagram of a kind of device checking shot-point offset that Fig. 2 provides for the embodiment of the present application.
Detailed description of the invention
For the technical scheme making those skilled in the art be more fully understood that in the application, below in conjunction with in the embodiment of the present application
Accompanying drawing, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described embodiment is only
It is some embodiments of the present application rather than whole embodiments.Based on the embodiment in the application, ordinary skill people
All other embodiments that member is obtained under not making creative work premise, all should belong to the scope of the application protection.
A kind of method flow diagram checking shot-point offset that Fig. 1 provides for the embodiment of the present application.Although be described below flow process include with
Multiple operations that particular order occurs, but it should be clearly understood that these processes can include more or less of operation, these
Operation can sequentially perform or executed in parallel (such as using parallel processor or multi-thread environment).As it is shown in figure 1, described method
May include that
S1: reading single big gun geological data, wherein, described single big gun geological data includes the geological data of multiple seismic channel.
On the oil-gas exploration Favorable Zones that geological work and other Geophysical Works primarily determine that, according to two dimension or stereo observing system
Arrange survey line and shot position, use dynamite source or controlled source earthquake-wave-exciting, and pass through cymoscope and seismic detector with the time
Discrete sampling mode records seismic wave field.Original geological data can be obtained by above-mentioned means.In upper original earthquake number
According to, on each geophone station, the seismic data of record is referred to as seismic channel.Seismic channel on one geophone station is a single track.Ground
Often excite one earthquake (also referred to as single big gun) during shake record, typically have thousands of cymoscopes receiving signal, the letter of reception
Number set be called single big gun collection record.The collection of the seismic channel of the repeatedly earthquake of cymoscope record is collectively referred to as many big guns collection record.Many big guns collection
Record is to be shown side by side on same geological data by multiple single big gun collection records.
The embodiment of the present application, when carrying out shot-point offset and checking, can check for the data of each big gun in many big guns collection record.
Specifically, the embodiment of the present application can read single big gun geological data, and described single big gun geological data includes the seismic channel of predetermined number
Data.Specifically, in described single big gun geological data, each seismic channel data can pass through xiT () represents, wherein i=1,2,3 ..., n,
Representing the number of channels of seismic channel in this list big gun geological data, t represents the time.
S2: pick up the actual first arrival of each seismic channel in described single big gun geological data.
The embodiment of the present application can judge described list by analyzing the first arrival of each seismic channel data in described single big gun geological data
Whether big gun geological data exists shot-point offset.At shot point after earthquake-wave-exciting, owing to the spread speed of various seismic waves is different,
The time traveling to observation station the most just has successively.When seismic wave arrives certain observation station, the cymoscope in observation station detects matter
The moment of some vibration is just properly termed as the first arrival time of ripple, is called for short first arrival.The application specifically can be picked up by following six steps
Take the actual first arrival of each seismic channel data in described single big gun geological data:
S21: choose the geological data of a seismic channel being not picked up actual first arrival in described single big gun geological data.
S22: the geological data chosen is carried out Hilbert transform.
Assume that the seismic channel data chosen can be designated as xiT (), wherein i represents the seismic channel that this seismic channel data chosen is corresponding
Number, t represents the time.The embodiment of the present application can carry out Hilbert transform according to the following formula to described seismic channel data:
Wherein, h (t) represents the seismic channel data after Hilbert-Huang Transform,Represent the shock response factor.
Seismic channel data after Hilbert transform, it keeps constant in the amplitude of each frequency component of frequency domain, but phase place will go out
Existing 90 ° of phase shifts.
S23: based on the geological data after the described geological data chosen and Hilbert transform, the earthquake chosen described in calculating
The instantaneous envelope that data are corresponding.
In certain embodiments, after described seismic channel data is carried out Hilbert transform, can be based on the described seismic channel chosen
Seismic channel data after data and Hilbert transform, the instantaneous envelope that the seismic channel data chosen described in calculating is corresponding.Specifically
Ground, the embodiment of the present application can by following formula calculate described in instantaneous envelope corresponding to the seismic channel data chosen:
Wherein, A (t) represent described in instantaneous envelope corresponding to the seismic channel data chosen, xiThe seismic channel chosen described in (t) representative
Data, h (t) represents the seismic channel data after Hilbert-Huang Transform.
The energy ratio of window time before and after S24: window division when described instantaneous envelope is carried out, and calculating.
In some embodiment, after the instantaneous envelope that the seismic channel data chosen described in calculating is corresponding, can be to described instantaneous bag
The energy ratio of window time before and after window division when network is carried out, and calculating.Specifically, it is assumed that time previous, the initial time of window is T1,
The termination time is T2, during later, window is when previous on the basis of window, and during holding, window window length is constant, when initial time and termination
Between all slide backward Δ t, then time previous, when window and later, the energy ratio of window can be represented by following formula:
Wherein, the energy ratio of window when R is before and after, A (t) be described in instantaneous envelope corresponding to the seismic channel data chosen, T1With
T2For the initial time of window time previous and the time of termination, the slip interval of window when Δ t is.
N number of there is long time window of studying in the same school mutually by being divided into by described instantaneous envelope, such that it is able to window time before and after calculating N-1
Energy ratio.Such as, window when described instantaneous envelope being divided into 10, then window time before and after just can calculating 9
Energy ratio.
S25: choose the time window that maximum in described energy ratio is corresponding, and by true for the maximum of described instantaneous envelope in window time this
The actual first arrival of the geological data chosen described in being set to.
In some embodiment, when calculating before and after after the energy ratio of window, maximum in described energy ratio can be chosen corresponding
Window time front and back.In 9 energy ratio as escribed above, the 2nd energy ratio maximum, then before and after this energy ratio correspondence
Time window window and window when the 3rd when being exactly the 2nd.Time in selecting described energy ratio before and after maximum correspondence after window, can
The actual first arrival of the seismic channel data to choose described in the maximum of described instantaneous envelope in window time before and after this is defined as.Need
Bright, in actual application scenarios, described before and after time window there may be overlap.Such as, time previous, window is from 10s to 30s,
During later, window is then probably from 12s to 32s, then time front and back, window possesses the overlap of 12s to 30s, then at described instantaneous bag
The maximum of network 10s to 32s time window in determine.
S26: repeat above-mentioned steps S21 to S25, until pickup obtains each seismic channel in described single big gun geological data
Actual first arrival.
After certain seismic channel data has picked up actual first arrival in described single big gun geological data, other can be not picked up actual first arrival
Seismic channel data perform step S21 to S25 step, until described single big gun geological data in all seismic channel data all pick up
Till complete actual first arrival.
S3: calculate the theoretical first arrival of each seismic channel in described single big gun geological data.
In some embodiment, in having picked up described single big gun geological data after the actual first arrival of each seismic channel data, can calculate
The theoretical first arrival of each seismic channel data in described single big gun geological data.Specifically, the embodiment of the present application can be by following two
Step calculates the theoretical first arrival of each seismic channel data in described single big gun geological data:
S31: according to actual first arrival and the geophone offset of each seismic channel of pickup, calculate speed parameter and time parameter.
The embodiment of the present application can obtain the first actual first arrival that in described single big gun geological data, the first seismic channel data is corresponding, first
Geophone offset and the second actual first arrival corresponding to the second seismic channel data, the second geophone offset, but speed can be calculated by following formula and join
Number and time parameter:
Wherein, t0Represent time parameter, v representation speed parameter, t1Represent the first actual first arrival, t2Represent the second actual first arrival, x1Generation
Table the first geophone offset, x2Represent the second geophone offset.
S32: according to preset rules, calculate each earthquake in described single big gun geological data by described speed parameter and time parameter
The theoretical first arrival in road.
After calculating speed parameter and time parameter, can be calculated by described speed parameter and time parameter according to preset rules
The theoretical first arrival of each seismic channel data in described single big gun geological data.Specifically, the embodiment of the present application can obtain described single big gun
The geophone offset that in geological data, each seismic channel data is corresponding, but can be by geophone offset corresponding for each seismic channel data generation respectively
Enter formulaCalculate the theoretical first arrival of each seismic channel data in described single big gun geological data, wherein, t0Represent
Described time parameter, v represents described speed parameter, tiRepresent the theoretical first arrival of i-th seismic channel data, xiRepresent i-th
The geophone offset that seismic channel data is corresponding.
S4: the actual first arrival of each seismic channel in described single big gun geological data and theoretical first arrival are contrasted, and according to contrast
Result judges whether described single big gun geological data exists shot-point offset.
The embodiment of the present application can be based on the above-mentioned actual first arrival determined and the theoretical first arrival of calculating, by contrast pass between the two
System, and judge whether described single big gun geological data exists shot-point offset according to the result of contrast.Specifically, the embodiment of the present application can
To preset first threshold and Second Threshold, described first threshold can be used to weigh described actual first arrival and described theoretical first arrival
Error relationship, described Second Threshold can be used to weigh the quantity of first arrival error occur.
Specifically, actual first arrival corresponding to each seismic channel data in described single big gun geological data and theoretical first arrival can be contrasted,
When error between described actual first arrival and theoretical first arrival is more than described first threshold, this seismic channel is defined as one and exists by mistake
The seismic channel of difference data.In actual application scenarios, described first threshold is often set to 0.2, i.e. when described actual first arrival and
When error between theoretical first arrival is more than 0.2, this seismic channel is defined as a seismic channel that there is error information.
After in described single big gun geological data, each seismic channel data carries out error judgement, when the earthquake of the existence error information determined
When the quantity in road is more than described Second Threshold, it is determined that described single big gun geological data exists shot-point offset.In actual application scenarios,
Described Second Threshold can also is that the rate of specific gravity that the quantity of the seismic channel that there is error information is shared in all seismic channel quantity.Tool
Body ground, this rate of specific gravity could be arranged to 60%, i.e. when the quantity of described error information is in all earthquakes of described single big gun geological data
When proportion value is more than 60% in the quantity of track data, then judge that described single big gun geological data exists shot-point offset.For existing
Single big gun seismic channel data of shot-point offset then needs to carry out Resurvey.
In actual application scenarios, due to by effect of noise, when the actual first arrival that each seismic channel data of pickup is corresponding,
The actual first arrival picked up there may be singular value.In the application one preferred embodiment, the actual first arrival after pickup can be done
Correcting process, to reduce the singular value in actual first arrival, thus optimizes final shot-point offset result of determination.Specifically, this Shen
By following four steps, the actual first arrival of pickup please can be modified by a preferred embodiment:
S27: calculate the average time difference of described actual first arrival.
According to the actual first arrival of described pickup, the embodiment of the present application can calculate the average time difference of described actual first arrival by following formula:
Wherein, tavRepresent the average time difference of described actual first arrival, tiRepresent each seismic channel data pair in described single big gun geological data
The actual first arrival answered, n represents the number of seismic channel in described single big gun geological data.
S28: according to the described average time difference, sets the 3rd threshold value.
In actual application scenarios, often the multiple that defines of the described average time difference being set to the 3rd threshold value, specifically, this defines
Multiple is 10 times.The most described 3rd threshold value could be arranged to 10tav。
S29: the 3rd corresponding for the 3rd seismic channel in described actual first arrival actual first arrival is contrasted with described 3rd threshold value,
When described 3rd actual first arrival is more than described three threshold value, from described actual first arrival, reject described 3rd actual first arrival.
After determining the 3rd threshold value, can be by actual first arrival corresponding for each seismic channel data in described single big gun geological data and institute
State the 3rd threshold value to contrast, and the actual first arrival that there is error bigger is weeded out.The embodiment of the present application is with the 3rd earthquake number of channels
As a example by according to, the 3rd corresponding for the 3rd seismic channel data in described actual first arrival actual first arrival can be carried out with described 3rd threshold value
Contrast, when described 3rd actual first arrival is more than described three threshold value, rejects described 3rd actual first arrival from described actual first arrival.
S210: interpolation calculation between the actual first arrival do not rejected is carried out, and using the result of calculating as described 3rd seismic channel pair
The actual first arrival answered.
After corresponding for described 3rd seismic channel data the 3rd actual first arrival being rejected, the embodiment of the present application can be by not rejecting
Actual first arrival carry out between interpolation calculation, and using the result that calculates as actual first arrival corresponding to described 3rd seismic channel data.This
Sample, less by other errors or there is not the actual first arrival of error and carry out interpolation calculation, just can obtain relatively reasonable reality
Border first arrival, thereby may be ensured that the smooth of actual first arrival in whole single big gun geological data, does not haves singular point, and then can protect
Demonstrate,prove the accuracy that final shot-point offset result judges.
The embodiment of the present application also provides for a kind of device checking shot-point offset.A kind of inspection big gun that Fig. 2 provides for the embodiment of the present application
The functional block diagram of the device of some skew.As in figure 2 it is shown, described device includes:
Geological data reads unit 100, is used for reading single big gun geological data, and wherein, described single big gun geological data includes multiple
The geological data of seismic channel；
Actual first break pickup unit 200, is used for picking up the actual first arrival of each seismic channel in described single big gun geological data；
Theoretical first arrival computing unit 300, is used for calculating the theoretical first arrival of each seismic channel in described single big gun geological data；
Shot-point offset identifying unit 400, is used for the actual first arrival of each seismic channel in described single big gun geological data and theoretical first arrival
Contrast, and judge whether described single big gun geological data exists shot-point offset according to the result of contrast.
In the application one preferred embodiment, described actual first break pickup unit 200 specifically may include that
Seismic channel data chooses module 201, is used for choosing an earthquake being not picked up actual first arrival in described single big gun geological data
The geological data in road；
Hilbert transform module 202, is used for the geological data chosen is carried out Hilbert transform；
Instantaneous envelope computing module 203, is used for based on the geological data after the described geological data chosen and Hilbert transform,
The instantaneous envelope that the seismic channel data chosen described in calculating is corresponding；
Energy ratio computing module 204, the energy ratio of window time when being used for described instantaneous envelope is carried out before and after window division, and calculating
Value；
Actual first arrival determines module 205, is used for choosing the time window that maximum in described energy ratio is corresponding, and by institute in window time this
State the maximum of instantaneous envelope be defined as described in the actual first arrival of seismic channel chosen；
Repeat module 206, be used for repeating described seismic channel data and choose module 201 and determine mould to described actual first arrival
The handling process of block 205, until pickup obtains the actual first arrival of each seismic channel in described single big gun geological data.
In the application one preferred embodiment, described shot-point offset identifying unit 400 specifically includes:
Threshold value setting module 401, is used for presetting first threshold and Second Threshold；
First determination module 402, is used for contrasting actual first arrival corresponding to each seismic channel in described single big gun geological data and theory
First arrival, when the error between actual first arrival and theoretical first arrival is more than described first threshold, is defined as an existence by this seismic channel
The seismic channel of error information；
Second determination module 403, is used for when the quantity of the seismic channel of the existence error information determined is more than described Second Threshold,
Judge that described single big gun geological data exists shot-point offset.
In order to be modified the actual first arrival of pickup, in another preferred embodiment of the application, described device can also include:
Average time difference computing unit 2051, is used for calculating the average time difference of described actual first arrival；
3rd threshold setting unit 2052, is used for, according to the described average time difference, setting the 3rd threshold value；
Actual first arrival culling unit 2053, is used for the 3rd corresponding for the 3rd seismic channel in described actual first arrival actual first arrival and institute
State the 3rd threshold value to contrast, when described 3rd actual first arrival is more than described three threshold value, from described actual first arrival, reject institute
State the 3rd actual first arrival；
Interpolation calculation unit 2054, interpolation calculation between being used for the actual first arrival do not rejected carried out, and using the result that calculates as
The actual first arrival that described 3rd seismic channel is corresponding.
Actual first arrival and the theoretical first arrival of calculating that the embodiment of the present application is picked up by contrast judge whether single big gun geological data is deposited
At shot-point offset, and can be modified processing to the actual first arrival of pickup, can apply to surface infrastructure complexity, top layer speed
Spend the area that cross directional variations is violent, and all processes all can be automatically obtained, it is not necessary to Quality Control personnel check, it is possible to quickly
Check shot-point offset exactly.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have many deformation and a change and
Without departing from spirit herein, it is desirable to appended claim includes that these deformation and change are without deviating from spirit herein.
Claims (10)
1. the method checking shot-point offset, it is characterised in that described method includes:
Reading single big gun geological data, wherein, described single big gun geological data includes the geological data of multiple seismic channel；
Pick up the actual first arrival of each seismic channel in described single big gun geological data；
Calculate the theoretical first arrival of each seismic channel in described single big gun geological data；
The actual first arrival of each seismic channel in described single big gun geological data and theoretical first arrival are contrasted, and according to the result of contrast
Judge whether described single big gun geological data exists shot-point offset.
2. the method checking shot-point offset as claimed in claim 1, it is characterised in that the described single big gun earthquake number of described pickup
According to, the actual first arrival of each seismic channel includes:
S21: choose the geological data of a seismic channel being not picked up actual first arrival in described single big gun geological data；
S22: the geological data chosen is carried out Hilbert transform；
S23: based on the geological data after the described geological data chosen and Hilbert transform, the earthquake chosen described in calculating
The instantaneous envelope that data are corresponding；
The energy ratio of window time before and after S24: window division when described instantaneous envelope is carried out, and calculating；
S25: choose the time window that maximum in described energy ratio is corresponding, and by true for the maximum of described instantaneous envelope in window time this
The actual first arrival of the geological data chosen described in being set to；
Repeat above-mentioned steps S21 to S25, until pickup obtains the reality of each seismic channel in described single big gun geological data
First arrival.
3. the method checking shot-point offset as claimed in claim 1, it is characterised in that the described single big gun earthquake number of described calculating
The theoretical first arrival of each seismic channel according to, including:
Actual first arrival according to pickup and the geophone offset of each seismic channel, calculate speed parameter and time parameter；
According to preset rules, calculate the reason of each seismic channel in described single big gun geological data by described speed parameter and time parameter
Opinion first arrival.
4. the method checking shot-point offset as claimed in claim 3, it is characterised in that the described actual first arrival according to pickup,
Calculate speed parameter and time parameter, including:
Obtain the first actual first arrival, the first geophone offset and the second seismic channel pair that in described single big gun geological data, the first seismic channel is corresponding
The second actual first arrival of answering, the second geophone offset；
Equation below is utilized to calculate speed parameter and time parameter:
Wherein, t0Represent time parameter, v representation speed parameter, t1Represent the first actual first arrival, t2Represent the second actual first arrival,
x1Represent the first geophone offset, x2Represent the second geophone offset.
5. the method checking shot-point offset as claimed in claim 3, it is characterised in that described according to preset rules, passes through
Described speed parameter and time parameter calculate the theoretical first arrival of each seismic channel in described single big gun geological data, including:
Obtain the geophone offset that in described single big gun geological data, each seismic channel is corresponding；
Geophone offset corresponding for each seismic channel is substituted into below equation respectively, is calculated each earthquake in described single big gun geological data
The theoretical first arrival in road:
Wherein, t0Representing described time parameter, v represents described speed parameter, tiRepresent the theoretical first arrival of i-th seismic channel, xi
Represent the geophone offset that i-th seismic channel is corresponding.
6. the method checking shot-point offset as claimed in claim 1, it is characterised in that the described single big gun earthquake number of described contrast
Actual first arrival that each seismic channel according to is corresponding and theoretical first arrival, and judge that described single big gun geological data is according to the result of contrast
No there is shot-point offset, including:
Preset first threshold and Second Threshold；
Contrast actual first arrival corresponding to each seismic channel in described single big gun geological data and theoretical first arrival, when actual first arrival and theory
When error between first arrival is more than described first threshold, this seismic channel is defined as a seismic channel that there is error information；
When the quantity of the seismic channel of the existence error information determined is more than described Second Threshold, it is determined that described single big gun geological data is deposited
At shot-point offset.
A kind of method checking shot-point offset the most as claimed in claim 1, it is characterised in that at the described single big gun of described contrast
Before actual first arrival that each seismic channel in geological data is corresponding and theoretical first arrival, described method also includes:
Calculate the average time difference of described actual first arrival；
According to the described average time difference, set the 3rd threshold value；
The 3rd corresponding for the 3rd seismic channel in described actual first arrival actual first arrival is contrasted with described 3rd threshold value, when described
When 3rd actual first arrival is more than described three threshold value, from described actual first arrival, reject described 3rd actual first arrival；
Interpolation calculation between the actual first arrival do not rejected is carried out, and using the result that calculates as reality corresponding to described 3rd seismic channel
Border first arrival.
8. the device checking shot-point offset, it is characterised in that described device includes:
Geological data reads unit, is used for reading single big gun geological data, and wherein, described single big gun geological data includes multiple earthquake
The geological data in road；
Actual first break pickup unit, is used for picking up the actual first arrival of each seismic channel in described single big gun geological data；
Theoretical first arrival computing unit, is used for calculating the theoretical first arrival of each seismic channel in described single big gun geological data；
Shot-point offset identifying unit, is used for carrying out the actual first arrival of each seismic channel in described single big gun geological data and theoretical first arrival
Contrast, and judge whether described single big gun geological data exists shot-point offset according to the result of contrast.
A kind of device checking shot-point offset the most as claimed in claim 8, it is characterised in that described actual first break pickup list
Unit specifically includes:
Seismic channel data chooses module, is used for choosing a seismic channel being not picked up actual first arrival in described single big gun geological data
Geological data；
Hilbert transform module, is used for the geological data chosen is carried out Hilbert transform；
Instantaneous envelope computing module, is used for based on the geological data after the described geological data chosen and Hilbert transform, meter
The instantaneous envelope that the geological data chosen described in calculation is corresponding；
Energy ratio computing module, the energy ratio of window time when being used for described instantaneous envelope is carried out before and after window division, and calculating；
Actual first arrival determines module, is used for choosing the time window that maximum in described energy ratio is corresponding, and by described wink in window time this
Time envelope maximum be defined as described in the actual first arrival of seismic channel data chosen；
Repeat module, be used for repeating described seismic channel data and choose module to described actual first arrival and determine the process of module
Flow process, until pickup obtains the actual first arrival of each seismic channel in described single big gun geological data.
A kind of device checking shot-point offset the most as claimed in claim 8, it is characterised in that described shot-point offset judges single
Unit specifically includes:
Threshold value setting module, is used for presetting first threshold and Second Threshold；
First determination module, is used for contrasting at the beginning of actual first arrival corresponding to each seismic channel in described single big gun geological data and theory
Extremely, when the error between actual first arrival and theoretical first arrival is more than described first threshold, this seismic channel is defined as one and exists by mistake
The seismic channel of difference data；
Second determination module, is used for when the quantity of the seismic channel of the existence error information determined is more than described Second Threshold, it is determined that
There is shot-point offset in described single big gun geological data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510379271.7A CN106324671A (en) | 2015-07-01 | 2015-07-01 | Method and device of examining shot point offset |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510379271.7A CN106324671A (en) | 2015-07-01 | 2015-07-01 | Method and device of examining shot point offset |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106324671A true CN106324671A (en) | 2017-01-11 |
Family
ID=57726755
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510379271.7A Pending CN106324671A (en) | 2015-07-01 | 2015-07-01 | Method and device of examining shot point offset |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106324671A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107576983A (en) * | 2017-11-02 | 2018-01-12 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Improve the method and data evaluation method of the single-shot earthquake record first break pickup degree of accuracy |
| CN110967758A (en) * | 2018-09-30 | 2020-04-07 | 中国石油化工股份有限公司 | Method and system for detecting spatial position of seismic acquisition excitation point |
| CN111352152A (en) * | 2018-12-21 | 2020-06-30 | 中国石油天然气集团有限公司 | Rapid quality control method and device for seismic data observation system |
| CN111965729A (en) * | 2019-05-20 | 2020-11-20 | 中国石油天然气集团有限公司 | Real-time monitoring method, system and device for vibroseis combination center |
| CN112241023A (en) * | 2019-07-16 | 2021-01-19 | 中国石油天然气集团有限公司 | Land three-dimensional shot point secondary positioning and quantitative deviation monitoring method and device |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101609160A (en) * | 2008-06-19 | 2009-12-23 | 中国石油天然气集团公司 | A kind of automatic seismic data first-arrival wave pickup method of strong background noise resistance |
| US20140219054A1 (en) * | 2013-02-04 | 2014-08-07 | King Fahd University Of Petroleum And Minerals | Method of first arrival picking of seismic refraction data |
| CN104597486A (en) * | 2013-10-31 | 2015-05-06 | 中国石油天然气集团公司 | Method for quickly detecting shot gather data through first arrival time |
-
2015
- 2015-07-01 CN CN201510379271.7A patent/CN106324671A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101609160A (en) * | 2008-06-19 | 2009-12-23 | 中国石油天然气集团公司 | A kind of automatic seismic data first-arrival wave pickup method of strong background noise resistance |
| US20140219054A1 (en) * | 2013-02-04 | 2014-08-07 | King Fahd University Of Petroleum And Minerals | Method of first arrival picking of seismic refraction data |
| CN104597486A (en) * | 2013-10-31 | 2015-05-06 | 中国石油天然气集团公司 | Method for quickly detecting shot gather data through first arrival time |
Non-Patent Citations (6)
| Title |
|---|
| JUAN I. SABBIONE ET AL.: "Automatic first-breaks picking: New strategies and algorithms", 《GEOPHYSICS》 * |
| 孙哲等: "地震数据高效采集实时质控技术", 《2015年物探技术研讨会》 * |
| 张伟等: "地震道瞬时强度比法拾取初至波", 《地球物理学进展》 * |
| 张艳: "黄土塬地区实验单炮初至拉平的实现", 《中国新技术新产品》 * |
| 徐钰等: "浅层折射波勘探中初至自动拾取新算法", 《石油地球物理勘探》 * |
| 陈良: "炮点空间位置检测技术研究", 《科技信息》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107576983A (en) * | 2017-11-02 | 2018-01-12 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Improve the method and data evaluation method of the single-shot earthquake record first break pickup degree of accuracy |
| CN107576983B (en) * | 2017-11-02 | 2019-06-04 | 中国石油集团东方地球物理勘探有限责任公司 | Improve the method and data evaluation method of the accuracy of single-shot earthquake record first break pickup |
| CN110967758A (en) * | 2018-09-30 | 2020-04-07 | 中国石油化工股份有限公司 | Method and system for detecting spatial position of seismic acquisition excitation point |
| CN110967758B (en) * | 2018-09-30 | 2022-03-08 | 中国石油化工股份有限公司 | Method and system for detecting spatial position of seismic acquisition excitation point |
| CN111352152A (en) * | 2018-12-21 | 2020-06-30 | 中国石油天然气集团有限公司 | Rapid quality control method and device for seismic data observation system |
| CN111965729A (en) * | 2019-05-20 | 2020-11-20 | 中国石油天然气集团有限公司 | Real-time monitoring method, system and device for vibroseis combination center |
| CN111965729B (en) * | 2019-05-20 | 2023-04-07 | 中国石油天然气集团有限公司 | Real-time monitoring method, system and device for vibroseis combination center |
| CN112241023A (en) * | 2019-07-16 | 2021-01-19 | 中国石油天然气集团有限公司 | Land three-dimensional shot point secondary positioning and quantitative deviation monitoring method and device |
| CN112241023B (en) * | 2019-07-16 | 2023-04-25 | 中国石油天然气集团有限公司 | Method and device for monitoring secondary positioning and quantitative offset of land three-dimensional shot point |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106324671A (en) | Method and device of examining shot point offset | |
| CN106199703B (en) | A kind of automatic positioning of microseism focus and Reliability Synthesis evaluation method | |
| CN106226818B (en) | Seismic data processing technique and device | |
| CN101609160B (en) | Automatic seismic data first-arrival wave pickup method with strong background noise resistance | |
| CN104155691B (en) | Converted wave anisotropic velocity analysis method and device | |
| CN107918156B (en) | Detect the method and device of subsea node acquisition earthquake data polarity | |
| CN104199103B (en) | Static correction method and static correction device for converted wave | |
| CN105204065B (en) | A kind of method and apparatus of first break picking ripple | |
| CN107728210B (en) | Multiple instruments acquire the determination method and apparatus that road is lacked in data | |
| CN104570076A (en) | Automatic seismic wave first-arrival picking method based on dichotomy | |
| CN103592683A (en) | Earth surface consistency vibration amplitude compensation method and device for converted wave seismic information | |
| CN103605157B (en) | The method of decay near surface scattering wave | |
| CN101598803A (en) | A kind of method that directly obtains stacked section of converted wave | |
| CN102540252A (en) | High-precision median stacking method on basis of cross-correlation | |
| CN105093319A (en) | Ground micro-seismic static correction method based on three-dimensional seismic data | |
| CN103592682A (en) | Seismic wave field energy balancing method and processing device | |
| CN105510975B (en) | Improve the method and device of geological data signal to noise ratio | |
| CN101852864B (en) | Method for processing mass seismic data by using surface consistent statistical spectrum analysis technology | |
| CN106842302B (en) | A kind of method and device of batch editor first arrival | |
| CN102830170A (en) | Control method and control device for obtaining coal sample transverse wave signal based on ultrasonic test | |
| US20230072038A1 (en) | Method and apparatus for removing tube wave interference from optical fiber acoustic wave sensing seismic data | |
| CN102692651A (en) | Preliminary wave residual static correction method with space-variant velocity | |
| CN106054252A (en) | Pre-stack time migration method and device | |
| CN107255832B (en) | A kind of inversion method of subsurface structure | |
| CN104865601B (en) | The method for determining the type of the earthquake record gathered in seismic prospecting |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170111 |
|
| WD01 | Invention patent application deemed withdrawn after publication |