CN104459779A - Method and device for automatically identifying abnormal seismic traces - Google Patents
Method and device for automatically identifying abnormal seismic traces Download PDFInfo
- Publication number
- CN104459779A CN104459779A CN201410746666.1A CN201410746666A CN104459779A CN 104459779 A CN104459779 A CN 104459779A CN 201410746666 A CN201410746666 A CN 201410746666A CN 104459779 A CN104459779 A CN 104459779A
- Authority
- CN
- China
- Prior art keywords
- value
- sampled point
- seismic
- road
- seismic trace
- 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the invention provides a method and device for automatically identifying abnormal seismic traces. The method includes the steps that true values of seismic trace sampling points are obtained from collected single-explosion seismic records, and a maximum value threshold parameter value is set according to the true values of the seismic trace sampling points; absolute values of the true values of the seismic trace sampling points are obtained and compared with the maximum value threshold parameter value; when the absolute values of the true values of the seismic trace sampling points are larger than the maximum value threshold parameter value, the seismic traces where the sampling points are located are the maximum-value-type abnormal seismic traces. By means of the technical scheme, the accuracy for identifying the abnormal seismic traces can be greatly improved.
Description
Technical field
The present invention relates to technical field of geophysical exploration, particularly relate to a kind of method and apparatus of automatic identification of Unusual Seismic road.
Background technology
Carrying out in the earthquake-capturing exploration process of field, the impact such as various external world or equipment faults itself to be subject to unavoidably, as being subject to impacts such as destroying unintentionally when the impact of extraneous natural weather, the impact of earthquake data acquisition stability of instrument, field seismometer wrong or embedding ropy impact and people, animal are walked about.It is unfavorable that these factors all can bring to the duty of collecting device, and then the data causing wave detector to receive show abnormal in seismologic record, reduce quality of seismic data.In order to ensure the quality of seismologic record, just need to carry out Timeliness coverage and rectification to various Unusual Seismic road.When constructing, the on-the-spot Quality Control personnel of main dependence check that papery list big gun seismologic record finds various types of Unusual Seismic road in real time in the past, then by its notice field line inspection personnel, rectify and improve in time it.
Along with the fast development of high density highly effective gathering technology, the single shot record number of channels of collection project is increased considerably, in field construction of today, single shot record receives number of channels often will reach thousands of road, road up to ten thousand or even tens thousand of road, and the on-the-spot Quality Control time leaving every big gun seismologic record for only have short a few second, and have every day thousands of big gun even several ten thousand big gun seismologic records need check.If continue to adopt the mode of artificial Quality Control the single shot record of collection in worksite to be carried out to the inspection in Unusual Seismic road, both can give the Quality Control personnel work load bringing tool large, the omission in a large amount of Quality Control problem can be caused again.Therefore, this mode can not meet artillery number big data quantity collection project to the requirement of the real-time Quality Control in Unusual Seismic road.
The existing noise road automatic identifying method based on Standard Average Amplitude and standard amplitude road at present, can identify the Unusual Seismic road of maximum value type, but the counting yield of the method is difficult to meet the Quality Control efficiency requirements of some highly effective gathering projects to image data; And when identifying maximum value type Unusual Seismic road, also some normal earthquake roads can be mistaken for abnormal road.Therefore, the method is mainly as an aid of indoor earthquake data process, and the accuracy of identification in its abnormal road does not reach the Quality Control accuracy requirement of field condition to image data.
Summary of the invention
The object of the application is to provide a kind of method and apparatus of automatic identification of Unusual Seismic road, to improve the recognition accuracy in Unusual Seismic road.
To achieve these goals, this application provides a kind of automatic knowledge method for distinguishing of Unusual Seismic road, the method comprises:
Obtain seismic trace sampled point true value from the single big gun seismologic record gathered, maximum value critical parameter value is set according to described seismic trace sampled point true value;
Obtain the absolute value of seismic trace sampled point true value, the absolute value of more described seismic trace sampled point true value and described maximum value critical parameter;
When the absolute value of described seismic trace sampled point true value is greater than described maximum value critical parameter, the seismic trace at described sampled point place is the Unusual Seismic road of maximum value type.
In a preferred embodiment, arrange maximum value critical parameter value according to described seismic trace sampled point true value to comprise:
The seismic trace sampled point true value that geophone offset is less than geophone offset threshold value is obtained from described seismic trace sampled point true value;
Calculate the absolute value of described sampled point true value, the absolute value obtaining described sampled point true value is less than or equal to the absolute value of the sampled point true value of the first threshold value;
Maximal value in the absolute value of described sampled point true value is multiplied by threshold value coefficient, as maximum value critical parameter.
In a preferred embodiment, the method also comprises: the Unusual Seismic road identification seismic trace except the Unusual Seismic road except described maximum value type being fallen to arrangement type.
In a preferred embodiment, the Unusual Seismic road identification falling arrangement type described in comprises:
The seismic trace sampled point true value that the Unusual Seismic road that is identified as maximum value type is corresponding is removed from the seismic trace sampled point true value that the single big gun seismologic record of described collection obtains;
Add up the sampled point number that the seismic trace sampled point true value of residue seismic trace is equal continuously, when described sampled point number is greater than to fall arrangement sampled point threshold value, then described sampled point place seismic trace is the Unusual Seismic road falling arrangement type.
In a preferred embodiment, the method also comprises: the Unusual Seismic road identification seismic trace except the Unusual Seismic road falling arrangement type except described maximum value being carried out to mono-tone interference type.
In a preferred embodiment, the Unusual Seismic road identification of described mono-tone interference type comprises:
Remove from the seismic trace sampled point true value that the single big gun seismologic record of described collection obtains and be identified as maximum value and fall seismic trace sampled point true value corresponding to the Unusual Seismic road of arrangement type;
Obtain the position of peak value in sampled point true value corresponding to each road seismic trace of residue;
Calculate the time difference between all adjacent peak positions, and calculate the mean value of the described time difference;
Calculate the difference of described each time difference and described time difference mean value, and calculate the absolute value of described difference, obtain the maximal value in described absolute value;
When described maximal value is less than mono-tone interference time difference threshold value, then described maximal value place seismic trace is the Unusual Seismic road of mono-tone interference type.
In a preferred embodiment, the method also comprises: carry out to the seismic trace except the Unusual Seismic road that described maximum value falls arrangement and mono-tone interference type the Unusual Seismic road identification being connected in series road type.
In a preferred embodiment, the Unusual Seismic road identification of described serial connection road type comprises:
Remove from the seismic trace sampled point true value that the single big gun seismologic record of described collection obtains and be identified as maximum value and fall seismic trace sampled point true value corresponding to the Unusual Seismic road of arrangement and mono-tone interference type;
Symbol together in the sampled point true value of position in the same time mutually before the seismic trace sampled point true value that more remaining each seismic trace is corresponding is adjacent, obtains the rate value that the identical sampled point of described symbol accounts for total sampled point;
When described rate value is greater than serial connection road sampled point threshold value, then described sampled point place seismic trace is the Unusual Seismic road of serial connection road type.
In a preferred embodiment, the method also comprises: do not jump the seismic trace except the Unusual Seismic road that described maximum value falls arrangement, mono-tone interference and serial connection road type, lazy Unusual Seismic road identification of jumping type.
In a preferred embodiment, describedly not jump, lazy Unusual Seismic road identification of jumping type comprise:
Window width when curve extends fixing to deep layer direction is cut, window when obtaining first arrival according to the top of seismic trace sampled point true value first arrival position pickup;
Remove from the seismic trace sampled point true value that the single big gun seismologic record of described collection obtains and be identified as seismic trace sampled point true value corresponding to Unusual Seismic road that maximum value falls arrangement, mono-tone interference and serial connection road type;
Sampled point true value during the seismic trace first arrival of statistics residue within the scope of window, calculates the absolute value of described sampled point true value, and calculates the mean value of the absolute value of all sampled point true value, obtain described mean value and the first number percent product;
Each seismic trace arranges a cluster counters, and described counter original state is zero, the absolute value of more described each sampled point true value and described mean value and the first number percent product;
When the absolute value of described sampled point true value is less than described mean value and the first number percent product, described counter adds 1, otherwise described counter is re-set as null value, when the value of described counter is greater than weak amplitude road sampled point threshold value, then described sampled point place seismic trace be do not jump, lazy Unusual Seismic road of jumping type.
The application also provides a kind of device of automatic identification of Unusual Seismic road on the other hand, and this device comprises:
Maximum value critical parameter value acquiring unit, for obtaining seismic trace sampled point true value from the single big gun seismologic record gathered, arranges maximum value critical parameter value according to described seismic trace sampled point true value;
Comparing unit, for obtaining the absolute value of seismic trace sampled point true value, the absolute value of more described seismic trace sampled point true value and described maximum value critical parameter;
Maximum value type identification unit, for when the absolute value of described seismic trace sampled point true value is greater than described maximum value critical parameter, the seismic trace at described sampled point place is the Unusual Seismic road of maximum value type.
In a preferred embodiment, this device also comprises:
Fall arrangement type recognition unit, for falling the Unusual Seismic road identification of arrangement type to the seismic trace except the Unusual Seismic road except described maximum value type.
In a preferred embodiment, this device also comprises:
Mono-tone interference type identification unit, for carrying out the Unusual Seismic road identification of mono-tone interference type to the seismic trace except the Unusual Seismic road falling arrangement type except described maximum value.
In a preferred embodiment, this device also comprises:
Serial connection road type identification unit, for carrying out to the seismic trace except the Unusual Seismic road that described maximum value falls arrangement and mono-tone interference type the Unusual Seismic road identification being connected in series road type.
In a preferred embodiment, this device also comprises:
Do not jump, lazy jump type identification unit, for not jumping the seismic trace except the Unusual Seismic road falling arrangement, mono-tone interference and serial connection road type except described maximum value, the Unusual Seismic road identification of lazy jumping type.
As can be seen here, the technical scheme of the embodiment of the present application utilizes seismic trace sampled point true value to pre-set maximum value critical parameter value, and the seismic trace sampled point true value of described maximum value critical parameter value and acquisition is compared and statistics, can identify that the Unusual Seismic road of maximum value type identifies automatically accurately.Compared with prior art, substantially increase the recognition accuracy in the Unusual Seismic road of maximum value type, and the erroneous judgement in the Unusual Seismic road to other types can be avoided.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present application or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, the accompanying drawing that the following describes is only some embodiments recorded in the application, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is the process flow diagram of the first embodiment of the automatic knowledge method for distinguishing in a kind of Unusual Seismic of the application road;
Fig. 2 is field seismic data transmission procedure schematic diagram;
Fig. 3 is the process flow diagram of the second embodiment of the automatic knowledge method for distinguishing in a kind of Unusual Seismic of the application road;
Fig. 4 is the process flow diagram of the 3rd embodiment of the automatic knowledge method for distinguishing in a kind of Unusual Seismic of the application road;
Fig. 5 is the process flow diagram of the 4th embodiment of the automatic knowledge method for distinguishing in a kind of Unusual Seismic of the application road;
Fig. 6 is the process flow diagram of the 5th embodiment of the automatic knowledge method for distinguishing in a kind of Unusual Seismic of the application road;
Fig. 7 is the process flow diagram of the 6th embodiment of the automatic knowledge method for distinguishing in a kind of Unusual Seismic of the application road;
Fig. 8 is the schematic diagram of the device of the automatic identification in a kind of Unusual Seismic road that the embodiment of the present application provides.
Embodiment
Technical scheme in the application is understood better in order to make those skilled in the art person, below in conjunction with the accompanying drawing in the embodiment of the present application, technical scheme in the embodiment of the present application is clearly and completely described, obviously, described embodiment is only some embodiments of the present application, instead of whole embodiments.Based on the embodiment in the application, those of ordinary skill in the art are not making the every other embodiment obtained under creative work prerequisite, all should belong to the scope of the application's protection.
The specific implementation of the embodiment of the present application is described below in detail with several concrete example.
Below the first embodiment of the automatic knowledge method for distinguishing in a kind of Unusual Seismic of the application road is first introduced.By reference to the accompanying drawings 1, this embodiment comprises:
S101: obtain seismic trace sampled point true value from the single big gun seismologic record gathered, maximum value critical parameter value is set according to described seismic trace sampled point true value.
In actual applications, need to utilize single wave detector or string combination to obtain single big gun seismologic record, utilize earthquake data acquisition instrument 201 by single big gun seismologic record copy by local network transport to data Quality Control main frame 202, data Quality Control main frame reads in time to it again, field seismic data transmission procedure as shown in Figure 2, wherein also comprises the instrument truck 203 placing earthquake data acquisition instrument 201;
Further, data Quality Control main frame 202 and earthquake data acquisition instrument 201 are carried out LAN (Local Area Network) online, in order to reach the object of the Unusual Seismic road in seismologic record being carried out to the real-time Quality Control of robotization, while making earthquake data acquisition instrument 201 often gather big gun data, can send a seismologic record copy to data Quality Control main frame 202 by LAN (Local Area Network), data Quality Control main frame 202 again reading is in time sent to the geological data copy of the machine and carries out follow-up robotization Quality Control.
Further, described maximum value critical parameter is according to single big gun seismologic record, seismic trace sampled point true value geophone offset in record being less than to geophone offset threshold value carries out pointwise screening, absolute value is greater than the first threshold value sampled point to reject, from all the other sampled points, obtain the maximal value of sampled point absolute value, this maximal value is multiplied by threshold value coefficient as maximum value critical parameter.
Further, the object arranging geophone offset threshold value reduces the calculated amount of search maximum value sampled point true value, improve search efficiency, according to energy transmission rule, in big gun record, the sampled point of the strongest amplitude only can occur within the scope of nearly geophone offset, the geophone offset threshold value of the embodiment of the present application is set to the geophone offset from shot point unidirectional spread the 20th road seismic trace from the near to the remote, but in actual applications, can set according to the actual conditions of the seismologic record gathered, the embodiment of the present application is not as limit.
Further, the first threshold value is here the maximum threshold value of acquisition instrument normal sample.The maximum threshold value of acquisition instrument normal sample can obtain by instrument configuration index, or takes an overall view of and examine nearly geophone offset seismic trace waveform, when waveform have cut top phenomenon time, cut top position sampling point absolute value and be the maximum threshold value Imax of normal sample.
Further, threshold value coefficient must be greater than 1, is usually set to 100, but in actual applications, can set according to the actual conditions of the seismologic record gathered, the embodiment of the present application is not as limit.By threshold coefficient being multiplied by the maximum threshold value of normal sample of acquisition, namely obtain maximum value critical parameter value.
S102: the absolute value obtaining seismic trace sampled point true value, the absolute value of more described seismic trace sampled point true value and described maximum value critical parameter.
Obtain the absolute value of each sampled point true value of the seismic trace of single big gun seismologic record of current reading, and compare with the maximum value critical parameter that obtains in step S101 according to this.
S103: when the absolute value of described seismic trace sampled point true value is greater than described maximum value critical parameter, the seismic trace at described sampled point place is the Unusual Seismic road of maximum value type.
A yard phenomenon is lost because instrumentation cables exists in the process of transmission data, thus cause the seismic trace sampled point true value distortion that collects, and this distortion to be usually expressed as absolute value more much larger than normal sample data absolute value, therefore by utilizing maximum value critical parameter, when judging that the absolute value of described seismic trace sampled point true value is greater than described maximum value critical parameter, the seismic trace at described sampled point place is the Unusual Seismic road of maximum value type.
As can be seen here, the technical scheme of the embodiment of the present application utilizes seismic trace sampled point true value to pre-set maximum value critical parameter value, and the seismic trace sampled point true value of described maximum value critical parameter value and acquisition is compared and statistics, can identify that the Unusual Seismic road of maximum value type identifies automatically accurately.Compared with prior art, substantially increase the recognition accuracy in the Unusual Seismic road of maximum value type, and the erroneous judgement in the Unusual Seismic road to other types can be avoided.
On the basis that above-mentioned first method is implemented, below provide the second embodiment of the automatic knowledge method for distinguishing in a kind of Unusual Seismic of the application road.By reference to the accompanying drawings 3, this embodiment comprises:
S301: remove the seismic trace sampled point true value that the Unusual Seismic road that is identified as maximum value type is corresponding from the seismic trace sampled point true value that the single big gun seismologic record of described collection obtains.
Further, to be identified as the Unusual Seismic selective trace muting of maximum value type, also namely remove all sampled point true value in maximum value type seismic trace, obtain seismic trace sampled point true value in remaining single big gun seismologic record, carry out the identification in the Unusual Seismic road of arrangement type.
S302: add up the sampled point number that the seismic trace sampled point true value of residue seismic trace is equal continuously, when described sampled point number is greater than to fall arrangement sampled point threshold value, then described sampled point place seismic trace is the Unusual Seismic road falling arrangement type.
Further, the sampled point threshold value that falls to arrange here is generally 40, but in actual applications, can set according to the actual conditions of the seismologic record gathered, the embodiment of the present application is not as limit.
In addition, what needs further illustrated is, seismic trace sampled point true value there will be the phenomenon that continuous a few road falls arrangement, and the usual performance characteristic of seismic trace falling arrangement is continuously identical, in this case, if first carry out the Unusual Seismic road identification being connected in series road type, then can be mistaken for serial connection seismic trace.Therefore, in order to ensure, in the wild in quality control process, accurately to identify the type in abnormal road, will prior to the Unusual Seismic road identification step to serial connection road type to the identification step in the Unusual Seismic road falling arrangement type.
As can be seen here, the technical scheme of the embodiment of the present application is after eliminating maximum value type Unusual Seismic road, what utilization pre-set arrange sampled point threshold value, and the statistics to the equal continuously sampled point number of numerical value in each seismic trace sampled point true value of the residue seismic trace behind the Unusual Seismic road removing maximum value type, when judging that described sampled point number is greater than to fall arrangement sampled point threshold value, can judge that described sampled point place seismic trace is the Unusual Seismic road falling arrangement type.Compared with prior art, substantially increase the recognition accuracy in the Unusual Seismic road of arrangement type, and the erroneous judgement in the Unusual Seismic road to other types can be avoided.
On the basis that first, second method above-mentioned is implemented, below provide the 3rd embodiment of the automatic knowledge method for distinguishing in a kind of Unusual Seismic of the application road.By reference to the accompanying drawings 4, this embodiment comprises:
S401: remove from the seismic trace sampled point true value that the single big gun seismologic record of described collection obtains and be identified as maximum value and fall seismic trace sampled point true value corresponding to the Unusual Seismic road of arrangement type.
Further, to be identified as maximum value and fall the Unusual Seismic selective trace muting of arrangement type, also namely removal maximum value falls all sampled point true value in arrangement type seismic trace, obtain seismic trace sampled point true value in remaining single big gun seismologic record, carry out the identification in the Unusual Seismic road of mono-tone interference type.
S402: the position obtaining peak value in sampled point true value corresponding to each road seismic trace of residue.
To remaining each seismic trace, need the position obtaining peak value in sampled point true value corresponding to each road seismic trace, i.e. P
i(i=1,2 ..., M), wherein, P
ibe i-th peak value position in this road; M is the peak value number in this road.
S403: calculate the time difference between all adjacent peak positions, and calculate the mean value of the described time difference.
Calculate between all adjacent peak positions time difference T
i, and calculate all T in this road
imean value T.
T
i＝P
i+1-P
i
Wherein, i=1,2 ..., M-1; P
i+1for the i-th+1 peak value position of this road seismic trace; P
ifor i-th peak value position of this road seismic trace.
S404: the difference calculating described each time difference and described time difference mean value, and calculate the absolute value of described difference, obtain the maximal value in described absolute value.
Further, each T is calculated
iwith the absolute value DT of the difference of mean value T
i, finally search out all absolute value DT again
imaximal value DTmax.
DT
i＝︳T-T
i︳
S405: when described maximal value is less than mono-tone interference time difference threshold value, then described maximal value place seismic trace is the Unusual Seismic road of mono-tone interference type.
Further, mono-tone interference time difference threshold value here can be set to 2ms, but in actual applications, can set according to the actual conditions of the seismologic record gathered, the embodiment of the present application is not as limit.
In addition, need to further illustrate, in producing in actual field, be subject to the impacts such as the interference of neighbouring generator, commercial power induction, geological data there will be the phenomenon that continuous a few road shows as mono-tone interference road, and the performance characteristic in mono-tone interference road is identical.In order to avoid being mistaken for serial connection seismic trace, will prior to the Unusual Seismic road identification step to serial connection road type to the identification step in the Unusual Seismic road of mono-tone interference type.
As can be seen here, the technical scheme of the embodiment of the present application is eliminating after maximum value falls arrangement type Unusual Seismic road, utilize the mono-tone interference time difference threshold value pre-set, and the statistics to peak value position in the sampled point true value of residue seismic trace, obtain the time difference between adjacent peak position, and calculate the mean value of the described time difference, calculate the difference of described each time difference and described time difference mean value further, calculate the absolute value of described difference again, obtain the maximal value in described absolute value; When described maximal value is less than mono-tone interference time difference threshold value, then described maximal value place seismic trace is the Unusual Seismic road of mono-tone interference type.Compared with prior art, substantially increase the recognition accuracy in the Unusual Seismic road of mono-tone interference type, and the erroneous judgement in the Unusual Seismic road to other types can be avoided.
On the basis that above-mentioned first, second, third method is implemented, below provide the 4th embodiment of the automatic knowledge method for distinguishing in a kind of Unusual Seismic of the application road.By reference to the accompanying drawings 5, this embodiment comprises:
S501: remove from the seismic trace sampled point true value that the single big gun seismologic record of described collection obtains and be identified as maximum value and fall seismic trace sampled point true value corresponding to the Unusual Seismic road of arrangement and mono-tone interference type.
Further, arrangement and the Unusual Seismic selective trace muting of mono-tone interference type is fallen by being identified as maximum value, also namely removal maximum value falls all sampled point true value in arrangement and mono-tone interference type seismic trace, obtain seismic trace sampled point true value in remaining single big gun seismologic record, carry out the identification in the Unusual Seismic road being connected in series road type.
S502: the symbol together in the sampled point true value of position in the same time mutually before the seismic trace sampled point true value that more remaining each seismic trace is corresponding is adjacent, obtains the rate value that the identical sampled point of described symbol accounts for total sampled point.
Further, the symbol of sampled point true value here, the i.e. positive number of described sampled point true value or negative, and obtain the rate value that the identical sampled point of symbol accounts for total sampled point.
S503: when described rate value is greater than serial connection road sampled point threshold value, then described sampled point place seismic trace is the Unusual Seismic road of serial connection road type.
Further serial connection road sampled point threshold value is here generally 40, but in actual applications, can set according to the actual conditions of the seismologic record gathered, the embodiment of the present application is not as limit.
Further, in the wild in gatherer process, often have the situation of wave detector wrong, twice or multiple tracks string are connected on together, and the waveform that seismologic record shows as serial connection road is almost identical.
As can be seen here, the technical scheme of the embodiment of the present application is eliminating after maximum value falls arrangement and mono-tone interference type Unusual Seismic road, utilize the serial connection road sampled point threshold value pre-set, and in the sampled point true value of residue seismic trace with the statistics of adjacent twice sampled point of same-sign in the sampled point true value of mutually position in the same time, obtain the rate value that the identical sampled point of described symbol accounts for total sampled point; When described rate value is greater than serial connection road sampled point threshold value, then described sampled point place seismic trace is the Unusual Seismic road of serial connection road type.Compared with prior art, substantially increase the recognition accuracy in the Unusual Seismic road of serial connection road type, and the erroneous judgement in the Unusual Seismic road to other types can be avoided.
On the basis that above-mentioned first, second, third and fourth method is implemented, below provide the 5th embodiment of the automatic knowledge method for distinguishing in a kind of Unusual Seismic of the application road.By reference to the accompanying drawings 6, this embodiment comprises:
S601: cut window width when curve extends fixing to deep layer direction according to the top of seismic trace sampled point true value first arrival position pickup, window when obtaining first arrival.
Further, during here fixing, window width is generally 200ms, but in actual applications, can set according to the actual conditions of the seismologic record gathered, the embodiment of the present application is not as limit.
S602: remove from the seismic trace sampled point true value that the single big gun seismologic record of described collection obtains and be identified as seismic trace sampled point true value corresponding to Unusual Seismic road that maximum value falls arrangement, mono-tone interference and serial connection road type.
Further, to be identified as maximum value and fall the Unusual Seismic selective trace muting of arrangement, mono-tone interference and serial connection road type, also namely removal maximum value falls all sampled point true value in arrangement, mono-tone interference and serial connection road type seismic trace, obtain seismic trace sampled point true value in remaining single big gun seismologic record, do not jump, lazy Unusual Seismic road of jumping type, also namely carry out the identification in the Unusual Seismic road of weak amplitude type.
S603: the sampled point true value of seismic trace sampled point true value when first arrival within the scope of window of statistics residue seismic trace, calculate the absolute value of described sampled point true value, and calculate the mean value of the absolute value of all sampled point true value, obtain described mean value and the first number percent product.
S604: each seismic trace arranges a cluster counters, and described counter original state is zero, the absolute value of more described each sampled point true value and described mean value and the first number percent product.
S605: when the absolute value of described sampled point true value is less than described mean value and the first number percent product, described counter adds 1, otherwise described counter is re-set as null value, when the value of described counter is greater than weak amplitude road sampled point threshold value, then described sampled point place seismic trace be do not jump, lazy Unusual Seismic road of jumping type.
Further, the first number percent is here generally 0.1%, but in actual applications, can set according to the actual conditions of the seismologic record gathered, the embodiment of the present application is not as limit.
Further, weak amplitude road sampled point threshold value is here generally 40, but in actual applications, can set according to the actual conditions of the seismologic record gathered, the embodiment of the present application is not as limit.
Further, in field acquisition process, when adopting the mode acquiring seismic data of receiver pattern, for in a string multiple geophone, if when only having minority wave detector normally work or all do not work, seismologic record will show as do not jump, lazy phenomenon of jumping, such seismic trace energy is very weak, therefore, by carrying out energy statistics to the seismologic record within the scope of window during first arrival but not within the scope of whole road, and increased substantially by the mode of point-by-point comparison do not jump, lazy accuracy of identification of jumping the Unusual Seismic road of type.
As can be seen here, the technical scheme of the embodiment of the present application is eliminating after maximum value falls arrangement, mono-tone interference and type Unusual Seismic road, serial connection road, the sampled point true value of window to remaining seismic trace when pre-setting first arrival is utilized to screen, calculate the absolute value of sampled point true value, and obtain the mean value of the absolute value of all sampled point true value, obtain described mean value and the first number percent product further; Again by arranging a cluster counters for each seismic trace, and described counter original state is zero, the absolute value of more described each sampled point true value and described mean value and the first number percent product; When the absolute value of described sampled point true value is less than described mean value and the first number percent product, described counter adds 1, otherwise described counter is null value, when the value of described counter is greater than weak amplitude road sampled point threshold value, then described sampled point place seismic trace be do not jump, lazy Unusual Seismic road of jumping type.Compared with prior art, substantially increase do not jump, lazy recognition accuracy of jumping the Unusual Seismic road of type, and the erroneous judgement in the Unusual Seismic road to other types can be avoided.
On the basis that above-mentioned first, second, third, fourth and fifth method is implemented, below provide the 6th embodiment of the automatic knowledge method for distinguishing in a kind of Unusual Seismic of the application road.By reference to the accompanying drawings 7, this embodiment comprises:
S701: the Unusual Seismic road of maximum value type is identified.
S702: the Unusual Seismic road falling arrangement type is identified.
S703: the Unusual Seismic road of mono-tone interference type is identified.
S704: the Unusual Seismic road of serial connection road type is identified.
S705: the Unusual Seismic road of not jumping, lazyness jumps type is identified.
As can be seen here, the technical scheme of the embodiment of the present application is by arranging rational recognition sequence and pointed recognition methods, effectively identify Unusual Seismic road, and determine affiliated type exactly, simultaneously because identifying is only carried out in time domain, and only added up by simple sample value and compare, not even with carrying out the computings such as interpolation, therefore, the monitoring efficiency identifying Unusual Seismic road can greatly be improved.Compared with prior art, substantially increase maximum value fall arrangement, mono-tone interference, serial connection road and do not jump, lazy recognition accuracy of jumping type Unusual Seismic road, and the erroneous judgement in the Unusual Seismic road to other types can be avoided.In seismic prospecting highly effective gathering project, there is very wide application prospect.
The embodiment of the present application is the corresponding device 800 providing a kind of automatic identification of Unusual Seismic road also, and as shown in Figure 8, this device comprises: maximum value critical parameter value acquiring unit 801, comparing unit 802 and maximum value type identification unit 803.
Wherein, maximum value critical parameter value acquiring unit 801, for obtaining seismic trace sampled point true value from the single big gun seismologic record gathered, arranges maximum value critical parameter value according to described seismic trace sampled point true value;
Comparing unit 802 for obtaining the absolute value of seismic trace sampled point true value, the absolute value of more described seismic trace sampled point true value and described maximum value critical parameter;
Maximum value type identification unit 803 is for when the absolute value of described seismic trace sampled point true value is greater than described maximum value critical parameter, and the seismic trace at described sampled point place is the Unusual Seismic road of maximum value type.
This device 800 also comprises in a preferred embodiment: fall arrangement type recognition unit, for falling the Unusual Seismic road identification of arrangement type to the seismic trace except the Unusual Seismic road except described maximum value type.
This device 800 also comprises in a preferred embodiment: mono-tone interference type identification unit, for carrying out the Unusual Seismic road identification of mono-tone interference type to the seismic trace except the Unusual Seismic road falling arrangement type except described maximum value.
This device 800 also comprises in a preferred embodiment: serial connection road type identification unit, for carrying out to the seismic trace except the Unusual Seismic road that described maximum value falls arrangement and mono-tone interference type the Unusual Seismic road identification being connected in series road type.
This device 800 also comprises in a preferred embodiment: do not jump, lazy jump type identification unit, for not jumping the seismic trace except the Unusual Seismic road falling arrangement, mono-tone interference and serial connection road type except described maximum value, the Unusual Seismic road identification of lazy jumping type.
As can be seen here, the technical scheme of the embodiment of the present application is by arranging rational recognition sequence and pointed recognition methods, effectively identify Unusual Seismic road, and determine affiliated type exactly, simultaneously because identifying is only carried out in time domain, and only added up by simple sample value and compare, not even with carrying out the computings such as interpolation, therefore, Unusual Seismic road can accurately be identified.Compared with prior art, substantially increase maximum value fall arrangement, mono-tone interference, serial connection road and do not jump, lazy recognition accuracy of jumping type Unusual Seismic road, and the erroneous judgement in the Unusual Seismic road to other types can be avoided.In seismic prospecting highly effective gathering project, there is very wide application prospect.
Each embodiment in this instructions all adopts the mode of going forward one by one to describe, and what each embodiment stressed is the difference with other embodiments, between each embodiment identical similar part mutually see.Especially, for system embodiment, because it is substantially similar to embodiment of the method, so description is fairly simple, relevant part illustrates see the part of embodiment of the method.
Although depict the application by embodiment, those of ordinary skill in the art know, the application has many distortion and change and do not depart from the spirit of the application, and the claim appended by wishing comprises these distortion and change and do not depart from the spirit of the application.
Claims (15)
1. the automatic knowledge method for distinguishing in Unusual Seismic road, it is characterized in that, the method comprises:
Obtain seismic trace sampled point true value from the single big gun seismologic record gathered, maximum value critical parameter value is set according to described seismic trace sampled point true value;
Obtain the absolute value of seismic trace sampled point true value, the absolute value of more described seismic trace sampled point true value and described maximum value critical parameter;
When the absolute value of described seismic trace sampled point true value is greater than described maximum value critical parameter, the seismic trace at described sampled point place is the Unusual Seismic road of maximum value type.
2. method according to claim 1, is characterized in that, describedly arranges maximum value critical parameter value according to described seismic trace sampled point true value and comprises:
The seismic trace sampled point true value that geophone offset is less than geophone offset threshold value is obtained from described seismic trace sampled point true value;
Calculate the absolute value of described sampled point true value, the absolute value obtaining described sampled point true value is less than or equal to the absolute value of the sampled point true value of the first threshold value;
Maximal value in the absolute value of described sampled point true value is multiplied by threshold value coefficient, as maximum value critical parameter.
3. method according to claim 1 and 2, is characterized in that, the method also comprises: the Unusual Seismic road identification seismic trace except the Unusual Seismic road except described maximum value type being fallen to arrangement type.
4. method according to claim 3, is characterized in that, described in fall arrangement type the identification of Unusual Seismic road comprise:
The seismic trace sampled point true value that the Unusual Seismic road that is identified as maximum value type is corresponding is removed from the seismic trace sampled point true value that the single big gun seismologic record of described collection obtains;
The sampled point number that in the per pass seismic trace of statistics residue seismic trace, sampled point true value is equal continuously, when described sampled point number is greater than to fall arrangement sampled point threshold value, then described sampled point place seismic trace is the Unusual Seismic road falling arrangement type.
5. method according to claim 4, is characterized in that, the method also comprises: the Unusual Seismic road identification seismic trace except the Unusual Seismic road falling arrangement type except described maximum value being carried out to mono-tone interference type.
6. method according to claim 5, is characterized in that, the Unusual Seismic road identification of described mono-tone interference type comprises:
Remove from the seismic trace sampled point true value that the single big gun seismologic record of described collection obtains and be identified as maximum value and fall seismic trace sampled point true value corresponding to the Unusual Seismic road of arrangement type;
Obtain the position of peak value in sampled point true value corresponding to each road seismic trace of residue;
Calculate the time difference between all adjacent peak positions, and calculate the mean value of the described time difference;
Calculate the difference of described each time difference and described time difference mean value, and calculate the absolute value of described difference, obtain the maximal value in described absolute value;
When described maximal value is less than mono-tone interference time difference threshold value, then described maximal value place seismic trace is the Unusual Seismic road of mono-tone interference type.
7. method according to claim 6, is characterized in that, the method also comprises: carry out to the seismic trace except the Unusual Seismic road that described maximum value falls arrangement and mono-tone interference type the Unusual Seismic road identification being connected in series road type.
8. method according to claim 7, is characterized in that, the Unusual Seismic road identification of described serial connection road type comprises:
Remove from the seismic trace sampled point true value that the single big gun seismologic record of described collection obtains and be identified as maximum value and fall seismic trace sampled point true value corresponding to the Unusual Seismic road of arrangement and mono-tone interference type;
Symbol together in the sampled point true value of position in the same time mutually before the seismic trace sampled point true value that more remaining each seismic trace is corresponding is adjacent, obtains the rate value that the identical sampled point of described symbol accounts for total sampled point;
When described rate value is greater than serial connection road sampled point threshold value, then described sampled point place seismic trace is the Unusual Seismic road of serial connection road type.
9. method according to claim 8, is characterized in that, the method also comprises: do not jump the seismic trace except the Unusual Seismic road that described maximum value falls arrangement, mono-tone interference and serial connection road type, lazy Unusual Seismic road identification of jumping type.
10. method according to claim 9, is characterized in that, describedly not to jump, lazy Unusual Seismic road identification of jumping type comprises:
Window width when curve extends fixing to deep layer direction is cut, window when obtaining first arrival according to the top of seismic trace sampled point true value first arrival position pickup;
Remove from the seismic trace sampled point true value that the single big gun seismologic record of described collection obtains and be identified as seismic trace sampled point true value corresponding to Unusual Seismic road that maximum value falls arrangement, mono-tone interference and serial connection road type;
The sampled point true value of seismic trace sampled point true value when first arrival within the scope of window of statistics residue seismic trace, calculates the absolute value of described sampled point true value, and calculates the mean value of the absolute value of all sampled point true value, obtain described mean value and the first number percent product;
Each seismic trace arranges a cluster counters, and described counter original state is zero, the more described absolute value of each sampled point true value and the product of described mean value and the first number percent;
When the absolute value of described sampled point true value is less than the product of described mean value and the first number percent, described counter adds 1, otherwise described counter is re-set as null value, when the value of described counter is greater than weak amplitude road sampled point threshold value, then described sampled point place seismic trace be do not jump, lazy Unusual Seismic road of jumping type.
The device of the automatic identification in 11. 1 kinds of Unusual Seismic roads, is characterized in that, this device comprises:
Maximum value critical parameter value acquiring unit, for obtaining seismic trace sampled point true value from the single big gun seismologic record gathered, arranges maximum value critical parameter value according to described seismic trace sampled point true value;
Comparing unit, for obtaining the absolute value of seismic trace sampled point true value, the absolute value of more described seismic trace sampled point true value and described maximum value critical parameter;
Maximum value type identification unit, for when the absolute value of described seismic trace sampled point true value is greater than described maximum value critical parameter, the seismic trace at described sampled point place is the Unusual Seismic road of maximum value type.
12. devices according to claim 11, is characterized in that, this device also comprises:
Fall arrangement type recognition unit, for falling the Unusual Seismic road identification of arrangement type to the seismic trace except the Unusual Seismic road except described maximum value type.
13. devices according to claim 12, is characterized in that, this device also comprises:
Mono-tone interference type identification unit, for carrying out the Unusual Seismic road identification of mono-tone interference type to the seismic trace except the Unusual Seismic road falling arrangement type except described maximum value.
14. devices according to claim 13, is characterized in that, this device also comprises:
Serial connection road type identification unit, for carrying out to the seismic trace except the Unusual Seismic road that described maximum value falls arrangement and mono-tone interference type the Unusual Seismic road identification being connected in series road type.
15. devices according to claim 14, is characterized in that, this device also comprises:
Do not jump, lazy jump type identification unit, for not jumping the seismic trace except the Unusual Seismic road falling arrangement, mono-tone interference and serial connection road type except described maximum value, the Unusual Seismic road identification of lazy 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 true CN104459779A (en) | 2015-03-25 |
| CN104459779B 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) |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105785437A (en) * | 2016-04-26 | 2016-07-20 | 中国石油天然气集团公司 | Abnormal seismic trace automatic discriminating method and device |
| CN106842302A (en) * | 2016-12-26 | 2017-06-13 | 中国石油天然气股份有限公司 | Method and device for editing first arrivals in batches |
| CN107102357A (en) * | 2016-02-23 | 2017-08-29 | 中国石油化工股份有限公司 | Eliminate the Processing Seismic Data and device of mono-tone interference |
| CN107272066A (en) * | 2017-06-22 | 2017-10-20 | 东华理工大学 | A kind of noisy seismic signal first-arrival traveltime pick-up method and device |
| CN107664778A (en) * | 2017-08-24 | 2018-02-06 | 中国石油天然气股份有限公司 | Method and device for determining seismic data quality |
| CN107884826A (en) * | 2017-12-13 | 2018-04-06 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | The method of dead trace recognition methods and optimization seismic data process in geological data |
| CN108181648A (en) * | 2017-11-23 | 2018-06-19 | 中国石油天然气集团公司 | The quality control method and device of continuous record aliased seismic gathered data |
| CN108732617A (en) * | 2017-04-21 | 2018-11-02 | 中国石油化工股份有限公司 | A kind of useless selective trace muting method and device of seismic data |
| CN109001798A (en) * | 2017-06-07 | 2018-12-14 | 中国石油化工股份有限公司 | A kind of method and system in automatic identification exception road |
| CN109100785A (en) * | 2018-06-15 | 2018-12-28 | 中国石油天然气集团有限公司 | The quality control method and device of continuously record data |
| CN109425894A (en) * | 2017-08-23 | 2019-03-05 | 中国石油化工股份有限公司 | A kind of seismic anomaly road detection method and device |
| CN109557583A (en) * | 2017-09-26 | 2019-04-02 | 中国石油化工股份有限公司 | A kind of seismic attribute extracting method and system |
| CN110308483A (en) * | 2019-05-23 | 2019-10-08 | 中国石油天然气股份有限公司 | Reflection coefficient solving method and device based on multitask Bayes compressed sensing |
| 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 |
| CN112305593A (en) * | 2019-07-23 | 2021-02-02 | 中国石油天然气集团有限公司 | Abnormal first arrival wave identification method and device |
Family Cites Families (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 |
| CN103439738B (en) * | 2013-08-27 | 2016-09-28 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Seismic prospecting single shot record exception road recognition methods |
| CN104122587B (en) * | 2014-06-30 | 2016-10-19 | 中国石油天然气集团公司 | A kind of abnormal first arrival recognition methods based on migration in offset domain and system |
-
2014
- 2014-12-09 CN CN201410746666.1A patent/CN104459779B/en active Active
Non-Patent Citations (3)
| Title |
|---|
| 孙哲，等: "VSP优化预测反褶积与VSP子波替换法反褶积", 《石油地球物理勘探》 * |
| 梁晓峰，等: "FTP实时质量监控技术", 《石油地球物理勘探》 * |
| 石双虎，等: "高效地震勘探数据采集智能化质控技术", 《石油地球物理勘探》 * |
Cited By (21)
| 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 |
| CN105785437A (en) * | 2016-04-26 | 2016-07-20 | 中国石油天然气集团公司 | Abnormal seismic trace automatic discriminating method and device |
| CN106842302A (en) * | 2016-12-26 | 2017-06-13 | 中国石油天然气股份有限公司 | Method and device for editing first arrivals in batches |
| CN106842302B (en) * | 2016-12-26 | 2019-06-11 | 中国石油天然气股份有限公司 | Method and device for editing first arrivals in batches |
| CN108732617A (en) * | 2017-04-21 | 2018-11-02 | 中国石油化工股份有限公司 | A kind of useless selective trace muting method and device of seismic data |
| CN109001798A (en) * | 2017-06-07 | 2018-12-14 | 中国石油化工股份有限公司 | A kind of method and system in automatic identification exception road |
| CN107272066A (en) * | 2017-06-22 | 2017-10-20 | 东华理工大学 | A kind of noisy seismic signal first-arrival traveltime pick-up method and device |
| 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 |
| CN107664778A (en) * | 2017-08-24 | 2018-02-06 | 中国石油天然气股份有限公司 | Method and device for determining seismic data quality |
| CN109557583A (en) * | 2017-09-26 | 2019-04-02 | 中国石油化工股份有限公司 | A kind of seismic attribute extracting method and system |
| CN108181648B (en) * | 2017-11-23 | 2019-07-09 | 中国石油天然气集团公司 | The quality control method and device of continuous record aliased seismic acquisition data |
| CN108181648A (en) * | 2017-11-23 | 2018-06-19 | 中国石油天然气集团公司 | The quality control method and device of continuous record aliased seismic gathered data |
| CN107884826A (en) * | 2017-12-13 | 2018-04-06 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | The method of dead trace recognition methods and optimization seismic data process in geological data |
| CN109100785A (en) * | 2018-06-15 | 2018-12-28 | 中国石油天然气集团有限公司 | The quality control method and device of continuously record 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 solving method and device based on multitask Bayes compressed sensing |
| CN112305593A (en) * | 2019-07-23 | 2021-02-02 | 中国石油天然气集团有限公司 | Abnormal first arrival wave identification method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104459779B (en) | 2017-07-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104459779A (en) | Method and device for automatically identifying abnormal seismic traces | |
| CN106022592B (en) | Electricity consumption behavior abnormity detection and public security risk early warning method and device | |
| CN104062677B (en) | Multifunctional comprehensive integrated high-precision intelligent micro-seismic monitoring system | |
| CN105931458B (en) | A kind of method of road traffic flow detection device reliability assessment | |
| CN112134615B (en) | Monitoring system, method, terminal and readable storage medium based on optical fiber link | |
| CN103596208A (en) | Method and system for judging fault of network element | |
| CN102735992A (en) | Surface potential pilot frequency comparison-based ground grid defect identification method and system | |
| CN109100782B (en) | Appointed shot gather record acquisition method, electronic equipment and computer storage medium | |
| CN105277968A (en) | Seismic excitation source automatic excitation control method and device | |
| CN116520433B (en) | Coal mine working face directional vibration pickup method for multichannel signal fusion analysis | |
| CN110580492A (en) | Track circuit fault precursor discovery method based on small fluctuation detection | |
| CN116986246A (en) | Intelligent inspection system and method for coal conveying belt | |
| CN104597486B (en) | Utilize the method for first arrival time quick detection shot gather data | |
| CN105425105A (en) | Method and system for comprehensive positioning of fault line of power grid | |
| 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 | |
| CN105785437A (en) | Abnormal seismic trace automatic discriminating method and device | |
| CN104459461B (en) | A kind of fault data auto-screening method based on power frequency Sine-Fitting coefficient correlation | |
| Romeu Petit et al. | Development of an earthquake early warning system based on earthworm: Application to Southwest Iberia | |
| CN104519511A (en) | Method and device for detecting scene breaks of communication network cells | |
| CN109444660A (en) | Method for identifying faults and interferences of power transmission line | |
| CN112379412B (en) | Quality monitoring method and device for collecting seismic data | |
| CN115408927A (en) | Data processing method and device for predicting rock mass parameters | |
| CN115144704A (en) | Fault detection method, device and equipment for cable production and storage medium | |
| CN104462820B (en) | A kind of method of power network shaft tower coordinate error detection |
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 |