CN104459779B - A kind of method and apparatus of the automatic identification in Unusual Seismic road - Google Patents
A kind of method and apparatus of the automatic identification in Unusual Seismic road Download PDFInfo
- Publication number
- CN104459779B CN104459779B CN201410746666.1A CN201410746666A CN104459779B CN 104459779 B CN104459779 B CN 104459779B CN 201410746666 A CN201410746666 A CN 201410746666A CN 104459779 B CN104459779 B CN 104459779B
- Authority
- CN
- China
- Prior art keywords
- value
- sampled point
- seismic
- road
- seismic channel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the present application provides a kind of method and apparatus of the automatic identification in Unusual Seismic road.This method includes：Seismic channel sampled point true value is obtained from the single-shot earthquake record of collection, maximum critical parameter value is set according to the seismic channel sampled point true value；Obtain the absolute value of seismic channel sampled point true value, the absolute value and the maximum critical parameter of relatively more described seismic channel sampled point true value；When the absolute value of the seismic channel sampled point true value is more than the maximum critical parameter, the seismic channel where the sampled point is the Unusual Seismic road of very big Value Types.The technical scheme provided using the embodiment of the present application can greatly improve the degree of accuracy of Unusual Seismic road identification.
Description
Technical field
The present invention relates to technical field of geophysical exploration, more particularly to a kind of method of the automatic identification in Unusual Seismic road
And device.
Background technology
, unavoidably will be by the shadow such as various external worlds or equipment faults itself in field earthquake-capturing exploration process is carried out
Ring, such as influenceed by extraneous natural weather, earthquake data acquisition stability of instrument influences, field seismometer wrong or embedding quality
The influence and people, animal of difference such as unintentionally destroy when walking about at the influence.These factors can be brought to the working condition of collecting device
It is unfavorable, and then the data for causing wave detector to receive show abnormal in earthquake record, reduce quality of seismic data.In order to
Ensure the quality of earthquake record, it is necessary to various Unusual Seismic roads are carried out to find and rectify and improve in time.In the past construction when, mainly according to
Papery single-shot earthquake record is checked by live Quality Control personnel to find various types of Unusual Seismic roads in real time, then it is wild by its notice
Outer line inspection personnel, are rectified and improved in time to it.
With the fast development of high density highly effective gathering technology so that the single shot record road number of collection project significantly increases
Plus, in field construction of today, single shot record reception channel number will often reach the even tens thousand of roads in thousands of roads, road up to ten thousand, and stay
The live Quality Control time to every big gun earthquake record only have it is short several seconds, and have daily thousands of big guns even tens of thousands of big gun earthquake records
Need to check.If continuing to carry out the inspection in Unusual Seismic road to the single shot record of collection in worksite by the way of artificial Quality Control,
Both it can be brought to Quality Control personnel and have big work load, the omission in a large amount of Quality Control problems can be caused again.Therefore, this mode has been
Through requirement of the artillery number big data quantity collection project to the real-time Quality Control in Unusual Seismic road can not be met.
The current existing noise road automatic identifying method based on Standard Average Amplitude and standard amplitude road, can be to very big
The Unusual Seismic road of Value Types is identified, but the computational efficiency of this method is difficult to meet some highly effective gathering projects to collection number
According to Quality Control efficiency requirements；And when very big Value Types Unusual Seismic road is identified, some normal earthquake roads can also be missed
It is judged to abnormal road.Therefore, the aid that this method is handled primarily as indoor earthquake data, the identification in its abnormal road
Precision does not reach Quality Control required precision of the field condition to gathered data.
The content of the invention
The purpose of the application is to provide a kind of method and apparatus of the automatic identification in Unusual Seismic road, to improve Unusual Seismic
The recognition accuracy in road.
To achieve these goals, this application provides a kind of method of the automatic identification in Unusual Seismic road, this method bag
Include：
Seismic channel sampled point true value is obtained from the single-shot earthquake record of collection, is set according to the seismic channel sampled point true value
Maximum critical parameter value；
Obtain the absolute value of seismic channel sampled point true value, the absolute value of relatively more described seismic channel sampled point true value and the pole
Big value critical parameter；
When the absolute value of the seismic channel sampled point true value is more than the maximum critical parameter, where the sampled point
Seismic channel be very big Value Types Unusual Seismic road.
In a preferred embodiment, maximum critical parameter value bag is set according to the seismic channel sampled point true value
Include：
The seismic channel sampled point true value that geophone offset is less than geophone offset threshold value is obtained from the seismic channel sampled point true value；
The absolute value of the sampled point true value is calculated, the absolute value for obtaining the sampled point true value is less than or equal to the first threshold
The absolute value of the sampled point true value of value；
Maximum in the absolute value of the sampled point true value is multiplied by threshold value coefficient, maximum critical parameter is used as.
In a preferred embodiment, this method also includes：To in addition to the Unusual Seismic road of the very big Value Types
Seismic channel carry out fall arrangement type Unusual Seismic road identification.
In a preferred embodiment, it is described fall arrangement type Unusual Seismic road identification include：
Removed in the seismic channel sampled point true value obtained from the single-shot earthquake record of the collection and be identified as maximum class
The corresponding seismic channel sampled point true value in Unusual Seismic road of type；
The continuous equal sampled point number of seismic channel sampled point true value of the remaining seismic channel of statistics, when the sampled point number
It is more than to fall arrangement sampled point threshold value, then seismic channel where the sampled point is the Unusual Seismic road for falling arrangement type.
In a preferred embodiment, this method also includes：Unusual Seismic to falling arrangement type except the maximum value
Seismic channel beyond road carries out the Unusual Seismic road identification of mono-tone interference type.
In a preferred embodiment, the Unusual Seismic road identification of the mono-tone interference type includes：
Removed in the seismic channel sampled point true value obtained from the single-shot earthquake record of the collection and be identified as maximum value and fall
The corresponding seismic channel sampled point true value in Unusual Seismic road of arrangement type；
Obtain the position of peak value in the corresponding sampled point true value of remaining each road seismic channel；
The time difference between all adjacent peak positions is calculated, and calculates the average value of the time difference；
The difference of each time difference and the time difference average value are calculated, and calculates the absolute value of the difference, institute is obtained
State the maximum in absolute value；
When the maximum is less than mono-tone interference time difference threshold value, then seismic channel where the maximum is mono-tone interference class
The Unusual Seismic road of type.
In a preferred embodiment, this method also includes：To falling arrangement and mono-tone interference type except the maximum value
Unusual Seismic road beyond seismic channel carry out concatenation road type Unusual Seismic road identification.
In a preferred embodiment, the Unusual Seismic road identification of concatenation road type includes：
Removed in the seismic channel sampled point true value obtained from the single-shot earthquake record of the collection and be identified as maximum value and fall
Arrangement seismic channel sampled point true value corresponding with the Unusual Seismic road of mono-tone interference type；
Compare the corresponding seismic channel sampled point true value of remaining each seismic channel it is adjacent thereto before together in mutually position in the same time
The symbol of the sampled point true value at place is put, the rate value that the symbol identical sampled point accounts for total sampled point is obtained；
When the rate value is more than concatenation road sampled point threshold value, then seismic channel where the sampled point is concatenation road type
Unusual Seismic road.
In a preferred embodiment, this method also includes：To falling arrangement, mono-tone interference and string except the maximum value
Connect the seismic channel beyond the Unusual Seismic road of type do not jumped, the lazy Unusual Seismic road identification for jumping type.
In a preferred embodiment, it is described do not jump, the lazy Unusual Seismic road identification for jumping type includes：
The top picked up according to seismic channel sampled point true value first arrival position cuts curve and extends window width when fixing to deep layer direction,
Obtain window during first arrival；
Removed in the seismic channel sampled point true value obtained from the single-shot earthquake record of the collection and be identified as maximum value and fall
Arrangement, mono-tone interference seismic channel sampled point true value corresponding with the Unusual Seismic road of concatenation road type；
Sampled point true value in the remaining seismic channel first arrival window scope of statistics, calculates the absolute value of the sampled point true value,
And the average value of the absolute value of all sampled point true value is calculated, obtain the average value and the first percentage product；
Each seismic channel sets one group of counter, and the counter original state is zero, relatively more described each sampled point
The absolute value of true value and the average value and the first percentage product；
When the absolute value of the sampled point true value is less than the average value and the first percentage product, the counter adds
1, otherwise the counter be re-set as null value, when the counter value be more than weak amplitude road sampled point threshold value when, then
Seismic channel where the sampled point is does not jump, the lazy Unusual Seismic road for jumping type.
On the other hand the application also provides a kind of device of the automatic identification in Unusual Seismic road, and the device includes：
Maximum critical parameter value acquiring unit is true for obtaining seismic channel sampled point from the single-shot earthquake record of collection
Value, maximum critical parameter value is set according to the seismic channel sampled point true value；
Comparing unit, the absolute value for obtaining seismic channel sampled point true value, relatively more described seismic channel sampled point true value
Absolute value and the maximum critical parameter；
Maximum type identification unit, for being more than maximum door when the absolute value of the seismic channel sampled point true value
During sill parameter, the seismic channel where the sampled point is the Unusual Seismic road of very big Value Types.
In a preferred embodiment, the device also includes：
Fall arrangement type recognition unit, for being carried out to the seismic channel in addition to the Unusual Seismic road of the very big Value Types
Fall the Unusual Seismic road identification of arrangement type.
In a preferred embodiment, the device also includes：
Mono-tone interference type identification unit, in addition to the Unusual Seismic road to falling arrangement type except the maximum value
Seismic channel carries out the Unusual Seismic road identification of mono-tone interference type.
In a preferred embodiment, the device also includes：
Road type identification unit is concatenated, for the Unusual Seismic to falling arrangement and mono-tone interference type except the maximum value
Seismic channel beyond road carries out the Unusual Seismic road identification of concatenation road type.
In a preferred embodiment, the device also includes：
Do not jump, lazyness jumps type identification unit, for falling arrangement, mono-tone interference and concatenation road type except the maximum value
Unusual Seismic road beyond seismic channel do not jumped, the lazy Unusual Seismic road identification for jumping type.
As can be seen here, the technical scheme of the embodiment of the present application pre-sets maximum threshold using seismic channel sampled point true value
Parameter value, and the maximum critical parameter value and the seismic channel sampled point true value of acquisition are compared and counted, Neng Gouzhun
The Unusual Seismic road of the true very big Value Types of identification carries out automatic identification.Compared with prior art, maximum class is substantially increased
The recognition accuracy in the Unusual Seismic road of type, and the erroneous judgement to other kinds of Unusual Seismic road can be avoided.
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 first embodiment of the method for the automatic identification in Unusual Seismic road of the application；
Fig. 2 is field seismic data transmission procedure schematic diagram；
Fig. 3 is a kind of flow chart of the second embodiment of the method for the automatic identification in Unusual Seismic road of the application；
Fig. 4 is a kind of flow chart of the 3rd embodiment of the method for the automatic identification in Unusual Seismic road of the application；
Fig. 5 is a kind of flow chart of the fourth embodiment of the method for the automatic identification in Unusual Seismic road of the application；
Fig. 6 is a kind of flow chart of 5th embodiment of the method for the automatic identification in Unusual Seismic road of the application；
Fig. 7 is a kind of flow chart of the sixth embodiment of the method for the automatic identification in Unusual Seismic road of the application；
Fig. 8 is a kind of schematic diagram of the device of the automatic identification in Unusual Seismic road 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 the method for the automatic identification in Unusual Seismic road of the application first below.With reference to accompanying drawing
1, the embodiment includes：
S101：Seismic channel sampled point true value is obtained from the single-shot earthquake record of collection, it is true according to the seismic channel sampled point
Value sets maximum critical parameter value.
In actual applications, it is necessary to obtain single-shot earthquake record using single wave detector or string combination, earthquake is utilized
Single-shot earthquake record copy is passed through local network transport to data Quality Control main frame 202, data Quality Control main frame by data acquisition instrument 201
It is read in time again, field seismic data transmission procedure is as shown in Fig. 2 wherein also include placing earthquake data acquisition instrument
The instrument cab 203 of device 201；
Further, that data Quality Control main frame 202 and earthquake data acquisition instrument 201 carried out into LAN is online, in order to up to
Unusual Seismic road in earthquake record automate the purpose of real-time Quality Control so that earthquake data acquisition instrument 201 is every
While gathering a big gun data, a earthquake record copy can be sent by LAN to data Quality Control main frame 202, data matter
Control main frame 202 reads and is sent to the geological data copy of the machine and carries out follow-up automation Quality Control in time again.
Further, described maximum critical parameter is, according to single-shot earthquake record, to be less than big gun to geophone offset in record
Examine the seismic channel sampled point true value away from threshold value and carry out pointwise screening, absolute value is rejected more than the first threshold value sampled point, from
The maximum of sampled point absolute value is obtained in remaining sampled point, the maximum is multiplied by threshold value coefficient and joined as maximum threshold
Number.
Further, the purpose for setting geophone offset threshold value is to reduce the calculating of search maximum value sampled point true value
Amount, improves search efficiency, according to energy transmission rule, and the sampled point of most strong amplitude only can be in the range of nearly geophone offset in big gun record
Occur, the geophone offset threshold value of the embodiment of the present application is set to the big gun of the 20th road seismic channel from shot point single ender from the near to the remote
Inspection away from, but in actual applications, can be set according to the actual conditions of the earthquake record of collection, the embodiment of the present application not with
This is limited.
Further, the first threshold value here is the maximum threshold value of acquisition instrument normal sample.Acquisition instrument is normally adopted
Sample maximum threshold value can be obtained in instrument configuration index, or is taken an overall view of and examined nearly geophone offset seismic channel waveform, when waveform has
When cutting top phenomenon, it is the maximum threshold value Imax of normal sample to cut top position sampled point absolute value.
Further, threshold value coefficient must be greater than 1, be usually arranged as 100, but in actual applications, can be according to adopting
The actual conditions of the earthquake record of collection are set, and the embodiment of the present application is not limited thereto.Obtained by the way that threshold coefficient is multiplied by
The maximum threshold value of normal sample obtained, that is, obtain maximum critical parameter value.
S102：Obtain the absolute value of seismic channel sampled point true value, relatively the absolute value of the seismic channel sampled point true value with
The maximum critical parameter.
Obtain the absolute value of each sampled point true value of the seismic channel of single-shot earthquake record currently read, and according to this with step
The maximum critical parameter obtained in rapid S101 is compared.
S103：When the absolute value of the seismic channel sampled point true value is more than the maximum critical parameter, the sampling
Seismic channel where point is the Unusual Seismic road of very big Value Types.
A yard phenomenon is lost because instrumentation cables exist during data are transmitted, so as to cause the seismic channel collected to be sampled
Point true value distortion, and to be usually expressed as absolute value more much larger than normal sample data absolute value for this distortion, therefore pass through
Using maximum critical parameter, when the absolute value for judging the seismic channel sampled point true value is more than the maximum critical parameter
When, the seismic channel where the sampled point is the Unusual Seismic road of very big Value Types.
As can be seen here, the technical scheme of the embodiment of the present application pre-sets maximum threshold using seismic channel sampled point true value
Parameter value, and the maximum critical parameter value and the seismic channel sampled point true value of acquisition are compared and counted, Neng Gouzhun
The Unusual Seismic road of the true very big Value Types of identification carries out automatic identification.Compared with prior art, maximum class is substantially increased
The recognition accuracy in the Unusual Seismic road of type, and the erroneous judgement to other kinds of Unusual Seismic road can be avoided.
On the basis of the implementation of above-mentioned first method, a kind of side of the automatic identification in Unusual Seismic road of the application presented below
The second embodiment of method.With reference to accompanying drawing 3, the embodiment includes：
S301：Remove and be identified as greatly in the seismic channel sampled point true value obtained from the single-shot earthquake record of the collection
The corresponding seismic channel sampled point true value in Unusual Seismic road of Value Types.
Further, the Unusual Seismic selective trace muting of maximum type will be identified as, namely with removing very big Value Types
All sampled point true value in road are shaken, seismic channel sampled point true value in remaining single-shot earthquake record is obtained, carry out class of permutation
The identification in the Unusual Seismic road of type.
S302：The continuous equal sampled point number of seismic channel sampled point true value of the remaining seismic channel of statistics, when the sampling
Point number is more than to fall arrangement sampled point threshold value, then seismic channel where the sampled point is the Unusual Seismic road for falling arrangement type.
Further, the sampled point threshold value that falls to arrange here is generally 40, but in actual applications, can be according to collection
The actual conditions of earthquake record are set, and the embodiment of the present application is not limited thereto.
Furthermore, it is necessary to further illustrate, seismic channel sampled point true value occurs that continuous several roads fall the phenomenon of arrangement, and
The usual performance characteristic of seismic channel for continuously falling arrangement is identical, in this case, if first carrying out the different of concatenation road type
Normal seismic channel identification, then can be mistaken for concatenating seismic channel.Therefore, in order to ensure in quality control process in the wild, to abnormal road
Type is accurately identified, will be prior to concatenating road type singularly to falling the identification step in Unusual Seismic road of arrangement type
Shake road identification step.
As can be seen here, the technical scheme of the embodiment of the present application is after very big Value Types Unusual Seismic road is eliminated, using pre-
What is first set arrange sampled point threshold value, and to the every of the remaining seismic channel behind the Unusual Seismic road of the very big Value Types of removal
The statistics of the continuous equal sampled point number of numerical value in individual seismic channel sampled point true value, when judging that the sampled point number is more than to fall
Arrange sampled point threshold value, it can be determined that seismic channel where the sampled point is the Unusual Seismic road for falling arrangement type.With it is existing
Technology is compared, and substantially increases the recognition accuracy in the Unusual Seismic road of arrangement type, and can be avoided to other types
Unusual Seismic road erroneous judgement.
On the basis of the implementation of above-mentioned first, second method, a kind of automatic knowledge in Unusual Seismic road of the application presented below
The 3rd embodiment of method for distinguishing.With reference to accompanying drawing 4, the embodiment includes：
S401：Remove and be identified as greatly in the seismic channel sampled point true value obtained from the single-shot earthquake record of the collection
It is worth, falls the corresponding seismic channel sampled point true value in Unusual Seismic road of arrangement type.
Further, maximum value will be identified as and has fallen the Unusual Seismic selective trace muting of arrangement type, namely removed very big
All sampled point true value for being worth, falling in arrangement type seismic channel, obtain seismic channel sampled point in remaining single-shot earthquake record true
Value, carries out the identification in the Unusual Seismic road of mono-tone interference type.
S402：Obtain the position of peak value in the corresponding sampled point true value of remaining each road seismic channel.
To remaining each seismic channel, it is necessary to obtain the position of peak value in the corresponding sampled point true value of each road seismic channel,
That is Pi(i=1,2 ... ..., M), wherein, PiAs i-th of peak value position in the road；M is the peak value number in the road.
S403：The time difference between all adjacent peak positions is calculated, and calculates the average value of the time difference.
Calculate the time difference T between all adjacent peak positionsi, and calculate all T in the roadiAverage value T.
Ti=Pi+1-Pi
Wherein, i=1,2 ... ..., M-1；Pi+1For the i+1 peak value position of the road seismic channel；PiIt is genuine for this
Shake i-th of peak value position in road.
S404：The difference of each time difference and the time difference average value are calculated, and calculates the absolute value of the difference, is obtained
Take the maximum in the absolute value.
Further, each T is calculatediWith the absolute value DT of average value T differencei, finally search again for out all absolute values
DTiMaximum DTmax.
DTi=︳ T-Ti︳
S405：When the maximum is less than mono-tone interference time difference threshold value, then seismic channel where the maximum is single-frequency
The Unusual Seismic road of interference type.
Further, mono-tone interference time difference threshold value here could be arranged to 2ms, but in actual applications, can basis
The actual conditions of the earthquake record of collection are set, and the embodiment of the present application is not limited thereto.
Furthermore, it is necessary to further illustrate, in being produced in actual field, by the interference of neighbouring generator, industry
Influenceed with electric induction etc., geological data occurs that continuous several roads show as the phenomenon in mono-tone interference road, and the table in mono-tone interference road
It is characterized in now identical.In order to avoid being mistaken for concatenating seismic channel, the identification to the Unusual Seismic road of mono-tone interference type is walked
Suddenly will be prior to the Unusual Seismic road identification step to concatenating road type.
As can be seen here, the technical scheme of the embodiment of the present application is being eliminated after maximum value falls arrangement type Unusual Seismic road,
Using the mono-tone interference time difference threshold value pre-set, and to peak value position in the sampled point true value of remaining seismic channel
Statistics, obtains the time difference between adjacent peak position, and calculates the average value of the time difference, further calculates each time difference
With the difference of the time difference average value, then the absolute value of the difference is calculated, obtain the maximum in the absolute value；When described
Maximum is less than mono-tone interference time difference threshold value, then seismic channel where the maximum is the Unusual Seismic of mono-tone interference type
Road.Compared with prior art, the recognition accuracy in the Unusual Seismic road of mono-tone interference type is substantially increased, and can be avoided
Erroneous judgement to other kinds of Unusual Seismic road.
Above-mentioned first, second, third method implementation on the basis of, a kind of Unusual Seismic road of the application presented below from
The dynamic fourth embodiment for knowing method for distinguishing.With reference to accompanying drawing 5, the embodiment includes：
S501：Remove and be identified as greatly in the seismic channel sampled point true value obtained from the single-shot earthquake record of the collection
It is worth, falls arrangement seismic channel sampled point true value corresponding with the Unusual Seismic road of mono-tone interference type.
Further, maximum value will be identified as and has fallen arrangement and the Unusual Seismic selective trace muting of mono-tone interference type, namely
Remove maximum value and fall arrangement and all sampled point true value in mono-tone interference type seismic channel, obtain remaining single-shot earthquake note
Seismic channel sampled point true value in record, carries out the identification in the Unusual Seismic road of concatenation road type.
S502：Compare the corresponding seismic channel sampled point true value of remaining each seismic channel it is adjacent thereto before together when identical
The symbol of the sampled point true value at position is carved, the rate value that the symbol identical sampled point accounts for total sampled point is obtained.
Further, the symbol of sampled point true value here, i.e., the positive number or negative of described sampled point true value, and obtain
Symbol identical sampled point accounts for the rate value of total sampled point.
S503：When the rate value is more than concatenation road sampled point threshold value, then seismic channel where the sampled point is concatenation
The Unusual Seismic road of road type.
Further concatenation road sampled point threshold value here is generally 40, but in actual applications, can be according to collection
The actual conditions of earthquake record are set, and the embodiment of the present application is not limited thereto.
Further, in the wild in gatherer process, often there are the situation of wave detector wrong, twice or multiple tracks wave detector concatenation
It is nearly identical in the waveform for, showing as concatenating together road in earthquake record.
As can be seen here, the technical scheme of the embodiment of the present application eliminate maximum value fall arrangement and mono-tone interference type it is different
After normal seismic channel, using the concatenation road sampled point threshold value pre-set, and in the sampled point true value of remaining seismic channel with
The statistics of adjacent twice sampled point of same-sign in the sampled point true value at phase in the same time position, obtains the symbol identical
Sampled point account for the rate value of total sampled point；Road sampled point threshold value is concatenated when the rate value is more than, then the sampled point institute
In Unusual Seismic road of the seismic channel for concatenation road type.Compared with prior art, concatenation road type is substantially increased singularly
The recognition accuracy in road is shaken, and the erroneous judgement to other kinds of Unusual Seismic road can be avoided.
On the basis of the implementation of above-mentioned first, second, third and fourth method, a kind of Unusual Seismic of the application presented below
5th embodiment of the method for the automatic identification in road.With reference to accompanying drawing 6, the embodiment includes：
S601：Cut curve and extend solid timing window to deep layer direction in the top picked up according to seismic channel sampled point true value first arrival position
Width, obtains window during first arrival.
Further, window width is generally 200ms during fixation here, but in actual applications, can be according to the ground of collection
The actual conditions of shake record are set, and the embodiment of the present application is not limited thereto.
S602：Remove and be identified as greatly in the seismic channel sampled point true value obtained from the single-shot earthquake record of the collection
It is worth, falls arrangement, mono-tone interference seismic channel sampled point true value corresponding with the Unusual Seismic road of concatenation road type.
Further, by be identified as maximum value fall arrangement, mono-tone interference and concatenate road type Unusual Seismic road pick
Remove, namely remove maximum value to fall arrangement, mono-tone interference and concatenate all sampled point true value in road type seismic channel, obtain surplus
Seismic channel sampled point true value in remaining single-shot earthquake record, is not jumped, lazyness jumps the Unusual Seismic road of type, namely carry out weak shake
The identification in the Unusual Seismic road of width type.
S603：Sampled point true value of the seismic channel sampled point true value of the remaining seismic channel of statistics in first arrival window scope, meter
The absolute value of the sampled point true value is calculated, and calculates the average value of the absolute value of all sampled point true value, the average value is obtained
With the first percentage product.
S604：Each seismic channel sets one group of counter, and the counter original state is zero, relatively more described each to adopt
The absolute value of sampling point true value and the average value and the first percentage product.
S605：When the absolute value of the sampled point true value is less than the average value and the first percentage product, the meter
Number device adds 1, and otherwise the counter is re-set as null value, when the value of the counter is more than weak amplitude road sampled point threshold value
When, then seismic channel where the sampled point be do not jump, the lazy Unusual Seismic road for jumping type.
Further, the first percentage here is generally 0.1%, but in actual applications, can be according to the earthquake of collection
The actual conditions of record are set, and the embodiment of the present application is not limited thereto.
Further, weak amplitude road sampled point threshold value here is generally 40, but in actual applications, can be according to collection
The actual conditions of earthquake record set, the embodiment of the present application is not limited thereto.
Further, during field acquisition, when gathering geological data by the way of detector array, for a string
, will table in earthquake record if there was only a small number of wave detector normal works or do not work all in multiple geophone
Now not jump, the lazy phenomenon jumped, such seismic channel energy is very weak, therefore, by first arrival window scope rather than whole
Earthquake record in the range of road carries out energy statisticses, and is greatly improved by way of point-by-point comparison and does not jump, lazy jump type
Unusual Seismic road accuracy of identification.
As can be seen here, the technical scheme of the embodiment of the present application eliminate maximum value fall arrangement, mono-tone interference and concatenation road
Behind type Unusual Seismic road, the sampled point true value of remaining seismic channel is screened using window when pre-setting first arrival, calculated
The absolute value of sampled point true value is obtained, and obtains the average value of the absolute value of all sampled point true value, further obtains described flat
Average and the first percentage product；Again by setting one group of counter for each seismic channel, and the counter original state is
Zero, the absolute value and the average value and the first percentage product of relatively more described each sampled point true value；When the sampled point is true
When the absolute value of value is less than the average value and the first percentage product, the counter adds 1, and otherwise the counter is zero
Value, when the value of the counter is more than weak amplitude road sampled point threshold value, then seismic channel where the sampled point not jump, it is lazy
Jump the Unusual Seismic road of type.Compared with prior art, substantially increase do not jump, the identification in the lazy Unusual Seismic road for jumping type it is accurate
True rate, and the erroneous judgement to other kinds of Unusual Seismic road can be avoided.
On the basis of the implementation of above-mentioned first, second, third, fourth and fifth method, the application one kind presented below is different
The sixth embodiment of the method for the automatic identification of normal seismic channel.With reference to accompanying drawing 7, the embodiment includes：
S701：The Unusual Seismic road of very big Value Types is identified.
S702：The Unusual Seismic road for falling arrangement type is identified.
S703：The Unusual Seismic road of mono-tone interference type is identified.
S704：The Unusual Seismic road for concatenating road type is identified.
S705：To not jumping, the lazy Unusual Seismic road for jumping type is identified.
As can be seen here, the technical scheme of the embodiment of the present application is by setting rational recognition sequence and with targetedly knowing
Other method, efficiently identifies out Unusual Seismic road, and determines affiliated type exactly, simultaneously because identification process is only in the time
Carry out, and only by simple sampling point Data-Statistics and compare on domain, even without the computings such as row interpolation are entered, accordingly, it is capable to significantly
Improve the monitoring efficiency in identification Unusual Seismic road.Compared with prior art, substantially increase maximum value fall arrangement, mono-tone interference,
Concatenation road and do not jump, the recognition accuracy in lazy jump type Unusual Seismic road, and can avoid to other kinds of Unusual Seismic
The erroneous judgement in road.In seismic prospecting highly effective gathering project, there is very wide application prospect.
The embodiment of the present application correspondingly provides a kind of device 800 of the automatic identification in Unusual Seismic road, as shown in figure 8, should
Device includes：Maximum critical parameter value acquiring unit 801, comparing unit 802 and maximum type identification unit 803.
Wherein, maximum critical parameter value acquiring unit 801 is used to adopt from the single-shot earthquake record acquisition seismic channel of collection
Sampling point true value, maximum critical parameter value is set according to the seismic channel sampled point true value；
Comparing unit 802 is used for the absolute value for obtaining seismic channel sampled point true value, relatively more described seismic channel sampled point true value
Absolute value and the maximum critical parameter；
Maximum type identification unit 803 is used for when the absolute value of the seismic channel sampled point true value is more than the maximum
During critical parameter, the seismic channel where the sampled point is the Unusual Seismic road of very big Value Types.
The device 800 also includes in a preferred embodiment：Fall arrangement type recognition unit, for except the pole
Seismic channel beyond the Unusual Seismic road of big Value Types carries out the Unusual Seismic road identification of arrangement type.
The device 800 also includes in a preferred embodiment：Mono-tone interference type identification unit, for except described
The seismic channel that maximum value falls beyond the Unusual Seismic road of arrangement type carries out the Unusual Seismic road identification of mono-tone interference type.
The device 800 also includes in a preferred embodiment：Road type identification unit is concatenated, for except the pole
Know in the Unusual Seismic road that the seismic channel be worth greatly, fallen beyond the Unusual Seismic road of arrangement and mono-tone interference type carries out concatenation road type
Not.
The device 800 also includes in a preferred embodiment：Do not jump, lazyness jumps type identification unit, for except institute
State the seismic channel that maximum value falls beyond arrangement, mono-tone interference and the Unusual Seismic road for concatenating road type do not jumped, lazy jump type
Unusual Seismic road identification.
As can be seen here, the technical scheme of the embodiment of the present application is by setting rational recognition sequence and with targetedly knowing
Other method, efficiently identifies out Unusual Seismic road, and determines affiliated type exactly, simultaneously because identification process is only in the time
Carry out, and only by simple sampling point Data-Statistics and compare on domain, even without the computings such as row interpolation are entered, accordingly, it is capable to accurately
Identification Unusual Seismic road.Compared with prior art, substantially increase maximum value and fall arrangement, mono-tone interference, concatenation road and not
Jump, lazyness jumps the recognition accuracy in type Unusual Seismic road, and the erroneous judgement to other kinds of Unusual Seismic road can be avoided.
In seismic prospecting highly effective gathering project, there is very wide application prospect.
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 (14)
1. a kind of method of the automatic identification in Unusual Seismic road, it is characterised in that this method includes：
Seismic channel sampled point true value is obtained from the single-shot earthquake record of collection, sets very big according to the seismic channel sampled point true value
It is worth critical parameter value；
Obtain the absolute value of seismic channel sampled point true value, the absolute value and the maximum of relatively more described seismic channel sampled point true value
Critical parameter；
When the absolute value of the seismic channel sampled point true value is more than the maximum critical parameter, the ground where the sampled point
Shake road is the Unusual Seismic road of very big Value Types；
Wherein, it is described to set maximum critical parameter value to include according to the seismic channel sampled point true value：Adopted from the seismic channel
The seismic channel sampled point true value that geophone offset is less than geophone offset threshold value is obtained in sampling point true value；
The absolute value of the sampled point true value is calculated, the absolute value for obtaining the sampled point true value is less than or equal to the first threshold value
The absolute value of sampled point true value；
Maximum in the absolute value of the sampled point true value is multiplied by threshold value coefficient, maximum critical parameter is used as.
2. according to the method described in claim 1, it is characterised in that this method also includes：To different except the very big Value Types
Seismic channel beyond normal seismic channel carries out the Unusual Seismic road identification of arrangement type.
3. method according to claim 2, it is characterised in that it is described fall arrangement type the identification of Unusual Seismic road include：
Removed in the seismic channel sampled point true value obtained from the single-shot earthquake record of the collection and be identified as very big Value Types
The corresponding seismic channel sampled point true value in Unusual Seismic road；
The continuous equal sampled point number of sampled point true value in the per pass seismic channel of the remaining seismic channel of statistics, when the sampled point
Number is more than to fall arrangement sampled point threshold value, then seismic channel where the sampled point is the Unusual Seismic road for falling arrangement type.
4. method according to claim 3, it is characterised in that this method also includes：To falling class of permutation except the maximum value
Seismic channel beyond the Unusual Seismic road of type carries out the Unusual Seismic road identification of mono-tone interference type.
5. method according to claim 4, it is characterised in that the Unusual Seismic road identification bag of the mono-tone interference type
Include：
Removed in the seismic channel sampled point true value obtained from the single-shot earthquake record of the collection and be identified as maximum value and fall arrangement
The corresponding seismic channel sampled point true value in Unusual Seismic road of type；
Obtain the position of peak value in the corresponding sampled point true value of remaining each road seismic channel；
The time difference between all adjacent peak positions is calculated, and calculates the average value of the time difference；
The difference of each time difference and the time difference average value are calculated, and calculates the absolute value of the difference, obtains described exhausted
To the maximum in value；
When the maximum is less than mono-tone interference time difference threshold value, then seismic channel where the maximum is mono-tone interference type
Unusual Seismic road.
6. method according to claim 5, it is characterised in that this method also includes：To except the maximum value fall arrangement and
Seismic channel beyond the Unusual Seismic road of mono-tone interference type carries out the Unusual Seismic road identification of concatenation road type.
7. method according to claim 6, it is characterised in that the Unusual Seismic road identification of concatenation road type includes：
Removed in the seismic channel sampled point true value obtained from the single-shot earthquake record of the collection and be identified as maximum value and fall arrangement
Seismic channel sampled point true value corresponding with the Unusual Seismic road of mono-tone interference type；
Compare the corresponding seismic channel sampled point true value of remaining each seismic channel first one adjacent thereto at phase in the same time position
Sampled point true value symbol, obtain the rate value that the symbol identical sampled point accounts for total sampled point；
When the rate value is more than concatenation road sampled point threshold value, then seismic channel where the sampled point is the different of concatenation road type
Normal seismic channel.
8. method according to claim 7, it is characterised in that this method also includes：To falling to arrange except the maximum value,
Seismic channel beyond mono-tone interference and the Unusual Seismic road of concatenation road type is not jumped, the lazy Unusual Seismic road for jumping type is known
Not.
9. method according to claim 8, it is characterised in that it is described do not jump, the lazy Unusual Seismic road identification bag for jumping type
Include：
The top picked up according to seismic channel sampled point true value first arrival position cuts curve and extends window width when fixing to deep layer direction, obtains
Window during first arrival；
Removed in the seismic channel sampled point true value obtained from the single-shot earthquake record of the collection and be identified as maximum value and fall row
Row, mono-tone interference seismic channel sampled point true value corresponding with the Unusual Seismic road of concatenation road type；
Sampled point true value of the seismic channel sampled point true value of the remaining seismic channel of statistics in first arrival window scope, calculates the sampling
The absolute value of point true value, and the average value of the absolute value of all sampled point true value is calculated, obtain the average value and the first percentage
Compare product；
Each seismic channel sets one group of counter, and the counter original state is zero, relatively more described each sampled point true value
Absolute value and the average value and the first percentage product；
When the absolute value of the sampled point true value is less than the product of the average value and the first percentage, the counter adds 1,
Otherwise the counter is re-set as null value, when the value of the counter is more than weak amplitude road sampled point threshold value, then institute
Seismic channel where stating sampled point is does not jump, the lazy Unusual Seismic road for jumping type.
10. a kind of device of the automatic identification in Unusual Seismic road, it is characterised in that the device includes：
Maximum critical parameter value acquiring unit, for obtaining seismic channel sampled point true value, root from the single-shot earthquake record of collection
According to the seismic channel sampled point true value, maximum critical parameter value is set；
Comparing unit, the absolute value for obtaining seismic channel sampled point true value compares the absolute of the seismic channel sampled point true value
Value and the maximum critical parameter；
Maximum type identification unit, joins for the absolute value when the seismic channel sampled point true value more than the maximum threshold
During number, the seismic channel where the sampled point is the Unusual Seismic road of very big Value Types；
Wherein, it is described to set maximum critical parameter value to include according to the seismic channel sampled point true value：Adopted from the seismic channel
The seismic channel sampled point true value that geophone offset is less than geophone offset threshold value is obtained in sampling point true value；
The absolute value of the sampled point true value is calculated, the absolute value for obtaining the sampled point true value is less than or equal to the first threshold value
The absolute value of sampled point true value；
Maximum in the absolute value of the sampled point true value is multiplied by threshold value coefficient, maximum critical parameter is used as.
11. device according to claim 10, it is characterised in that the device also includes：
Fall arrangement type recognition unit, for carrying out row to the seismic channel in addition to the Unusual Seismic road of the very big Value Types
The Unusual Seismic road identification of row type.
12. device according to claim 11, it is characterised in that the device also includes：
Mono-tone interference type identification unit, for the earthquake in addition to the Unusual Seismic road to falling arrangement type except the maximum value
Road carries out the Unusual Seismic road identification of mono-tone interference type.
13. device according to claim 12, it is characterised in that the device also includes：
Concatenate road type identification unit, for except the maximum value fall arrangement and mono-tone interference type Unusual Seismic road with
Outer seismic channel carries out the Unusual Seismic road identification of concatenation road type.
14. device according to claim 13, it is characterised in that the device also includes：
Do not jump, lazyness jumps type identification unit, for falling the different of arrangement, mono-tone interference and concatenation road type except the maximum value
Seismic channel beyond normal seismic channel is not jumped, the lazy Unusual Seismic road identification for jumping type.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410746666.1A CN104459779B (en) | 2014-12-09 | 2014-12-09 | A kind of method and apparatus of the automatic identification in Unusual Seismic road |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410746666.1A CN104459779B (en) | 2014-12-09 | 2014-12-09 | A kind of method and apparatus of the automatic identification in Unusual Seismic road |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104459779A CN104459779A (en) | 2015-03-25 |
| CN104459779B true CN104459779B (en) | 2017-07-28 |
Family
ID=52906124
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410746666.1A Active CN104459779B (en) | 2014-12-09 | 2014-12-09 | A kind of method and apparatus of the automatic identification in Unusual Seismic road |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104459779B (en) |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107102357A (en) * | 2016-02-23 | 2017-08-29 | 中国石油化工股份有限公司 | Eliminate the Processing Seismic Data and device of mono-tone interference |
| CN105785437B (en) * | 2016-04-26 | 2018-03-13 | 中国石油天然气集团公司 | A kind of method and device of the automatic discrimination in Unusual Seismic road |
| CN106842302B (en) * | 2016-12-26 | 2019-06-11 | 中国石油天然气股份有限公司 | A kind of method and device of batch editor first arrival |
| CN108732617B (en) * | 2017-04-21 | 2020-04-07 | 中国石油化工股份有限公司 | Seismic data waste channel rejection method and device |
| CN109001798B (en) * | 2017-06-07 | 2020-04-07 | 中国石油化工股份有限公司 | Method and system for automatically identifying abnormal lane |
| CN107272066B (en) * | 2017-06-22 | 2019-01-25 | 东华理工大学 | A kind of noisy seismic signal first-arrival traveltime pick-up method and device |
| CN109425894A (en) * | 2017-08-23 | 2019-03-05 | 中国石油化工股份有限公司 | A kind of seismic anomaly road detection method and device |
| CN107664778B (en) * | 2017-08-24 | 2019-11-08 | 中国石油天然气股份有限公司 | The determination method and apparatus of earthquake data quality |
| CN109557583B (en) * | 2017-09-26 | 2020-12-01 | 中国石油化工股份有限公司 | Seismic attribute extraction method and system |
| CN108181648B (en) * | 2017-11-23 | 2019-07-09 | 中国石油天然气集团公司 | The quality control method and device of continuous record aliased seismic acquisition data |
| CN107884826A (en) * | 2017-12-13 | 2018-04-06 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | The method of dead trace recognition methods and optimization seismic data process in geological data |
| CN109100785B (en) * | 2018-06-15 | 2020-09-08 | 中国石油天然气集团有限公司 | Quality control method and device for continuous recording data |
| CN110927782A (en) * | 2018-09-20 | 2020-03-27 | 中国石油化工股份有限公司 | Effective signal determination method and system for UNITE seismic acquisition instrument data |
| CN111381274A (en) * | 2018-12-29 | 2020-07-07 | 中国石油天然气股份有限公司 | Transmission fault channel identification method and device |
| CN110308483A (en) * | 2019-05-23 | 2019-10-08 | 中国石油天然气股份有限公司 | Reflection coefficient acquiring method and device based on multitask Bayes's compressed sensing |
| CN112305593B (en) * | 2019-07-23 | 2023-07-25 | 中国石油天然气集团有限公司 | Abnormal first arrival wave identification method and device |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4953142A (en) * | 1989-01-06 | 1990-08-28 | Marathon Oil Company | Model-based depth processing of seismic data |
| CN103439738A (en) * | 2013-08-27 | 2013-12-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Identification method for seismic exploration single-shot record abnormal channels |
| CN104122587A (en) * | 2014-06-30 | 2014-10-29 | 中国石油天然气集团公司 | Abnormity first appearance recognition method and system based on offset domain |
-
2014
- 2014-12-09 CN CN201410746666.1A patent/CN104459779B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4953142A (en) * | 1989-01-06 | 1990-08-28 | Marathon Oil Company | Model-based depth processing of seismic data |
| CN103439738A (en) * | 2013-08-27 | 2013-12-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Identification method for seismic exploration single-shot record abnormal channels |
| CN104122587A (en) * | 2014-06-30 | 2014-10-29 | 中国石油天然气集团公司 | Abnormity first appearance recognition method and system based on offset domain |
Non-Patent Citations (3)
| Title |
|---|
| FTP实时质量监控技术;梁晓峰，等;《石油地球物理勘探》;20081231;第43卷;101-103 * |
| VSP优化预测反褶积与VSP子波替换法反褶积;孙哲，等;《石油地球物理勘探》;20091031;第44卷(第5期);569-573 * |
| 高效地震勘探数据采集智能化质控技术;石双虎，等;《石油地球物理勘探》;20131231;第48卷;7-11，46 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104459779A (en) | 2015-03-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104459779B (en) | A kind of method and apparatus of the automatic identification in Unusual Seismic road | |
| FI107083B (en) | Method and arrangement for detecting the cause of a partial discharge | |
| CN104062677B (en) | Multifunctional comprehensive integrated high-precision intelligent micro-seismic monitoring system | |
| US6507181B1 (en) | Arrangement and method for finding out the number of sources of partial discharges | |
| CN104747163B (en) | A kind of method and device that Reservoir Fracture is recognized in tight sand | |
| CN105261136B (en) | The method and device of weather interference is shielded in a kind of fiber-optic monitoring warning system | |
| CN108919353A (en) | A kind of automatic classification of microseism waveform first arrival-time picks up and preferred method | |
| CN105527650A (en) | Automatic identification algorithm for microseismic signal and p wave first arrival at engineering scale | |
| CN112083263B (en) | Multi-source data fusion high-voltage cable line fault online positioning device and method based on traveling wave principle | |
| CN116520433B (en) | Coal mine working face directional vibration pickup method for multichannel signal fusion analysis | |
| CN105515888A (en) | Intelligent substation communication network anomaly detection method based on multi-dimensional entropy sequence classification | |
| CN105549083A (en) | Method of determining seismic observation system repeatability integration measurement and apparatus thereof | |
| CN104597486B (en) | Utilize the method for first arrival time quick detection shot gather data | |
| CN104977602B (en) | A kind of control method and device of earthquake data acquisition construction | |
| US20100271044A1 (en) | Remote detection of discharge on a power line network | |
| CN114089332A (en) | Geological disaster early warning system based on radar and remote sensing technology | |
| CN106019089A (en) | Method for carrying out partial discharging determination based on related relationship feature of alternate signals | |
| Romeu Petit et al. | Development of an earthquake early warning system based on earthworm: Application to Southwest Iberia | |
| CN109444660A (en) | Method for identifying faults and interferences of power transmission line | |
| CN104519511A (en) | Method and device for detecting scene breaks of communication network cells | |
| CN105866552B (en) | The measurement method of the impedance of Aircraft Cable Shield layer | |
| CN111123356B (en) | Abnormal track intelligent identification method based on first arrival information | |
| Kinali et al. | Detection of weak seismic signals in noisy environments from unfiltered, continuous passive seismic recordings | |
| Grebovic et al. | Analytics of lightning data collected by the direct measurement system and lightning location system | |
| CN112379412A (en) | Quality monitoring method and device for collecting seismic data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |