CN105445793A - Method and apparatus for determining bad trace data - Google Patents
Method and apparatus for determining bad trace data Download PDFInfo
- Publication number
- CN105445793A CN105445793A CN201510852243.2A CN201510852243A CN105445793A CN 105445793 A CN105445793 A CN 105445793A CN 201510852243 A CN201510852243 A CN 201510852243A CN 105445793 A CN105445793 A CN 105445793A
- Authority
- CN
- China
- Prior art keywords
- channel data
- seismic channel
- data
- seismic
- error correction
- 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
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000012937 correction Methods 0.000 claims abstract description 118
- 238000005070 sampling Methods 0.000 claims description 9
- 230000001143 conditioned effect Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 12
- 238000012216 screening Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000019771 cognition Effects 0.000 description 2
- 239000002360 explosive Substances 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method and apparatus for determining bad trace data. The method for determining bad trace data comprises: performing dynamic correction on the collected seismic trace data, and obtaining the seismic trace data after dynamic correction; acquiring first seismic trace data in a preset time window from the seismic trace data after dynamic correction; performing time difference correction on the seismic trace data after dynamic correction, and obtaining the seismic trace data after time difference correction; acquiring second seismic trace data in the preset time window from the seismic trace data after time difference correction; determining the correlation between the first seismic trace data and the second seismic trace data; and according to the correlation, determining whether the collected seismic trace data is bad trace data. The method and apparatus for determining bad trace data can avoid errors caused by artificial screening and guarantee no occurrence of bad trace data during the subsequent data processing process.
Description
Technical field
The application relates to seismic data processing technology field, particularly a kind of method and device determining bad track data.
Background technology
Often following three phases can be divided into the process that geological data processes:
1) first stage, in the earthquake data acquisition stage, the working contents of this one-phase is: in survey area, arranges two dimension or three-dimensional survey line; Use explosive source or vibroseis earthquake-wave-exciting, wherein, explosive source or vibroseis point can be called shot point; Equidistantly arrange that multiple wave detector is to receive seismic signal along survey line, the quantity of wave detector or geophone group can set as required, between each geophone group, arrangement mode can be different, as straddle spread, offend spread etc., also can set as required; After wave detector receives seismic signal, with constant duration discrete sampling geological data, and be recorded in digital form on tape.
2) second stage, seismic data processing stage, the working contents of this one-phase is: based on seismic wave propagation theory, utilize computing machine and Seismic data processing software, the geological data of field acquisition in the processing process first stage, becomes by geological data into the seismic cross-section that can reflect underground structure and can reflect the information sectional views such as seismic amplitude, frequency and velocity of propagation that subsurface lithologic changes.
3) three phases, the seismic data interpretation stage, the working contents of this one-phase is: the data such as the seismic cross-section that in analysis interpretation subordinate phase, seismic data process obtains, according to petroleum geology principle and seismic wave propagation theory, determine the structure of subterranean strata, thus find out the Favorable Zones of oily and propose drilling well position.
Can be found by the process of above-mentioned seismic prospecting, earthquake data acquisition is the basis of latter earthquake exploration steps, and the quality of earthquake data acquisition directly will affect the result of finally seismic exploration.In addition, earthquake data acquisition parameter is the key factor determining earthquake data quality.Therefore, how from the geological data gathered, to reject bad track data and will seem very important.
The mode of the bad track data in current identification geological data mainly contains two kinds, is respectively artificial cognition and machine pressing.But these two kinds of methods all exist certain drawback, and wherein, artificial cognition often introduces too much human error, thus make the identification of bad track data thorough not, finally can cause still having part bad track parameter to retain; Bad track data can not really be removed by machine pressing from geological data, but bad track data are suppressed in the seismic data, to reduce its inaccuracy in data processing, but along with the increase of subsequent processing steps and the changeable of data processing form, downtrodden bad track data often reappear, thus affect the accuracy of follow-up data process.
Above it should be noted that, just conveniently to the technical scheme of the application, clear, complete explanation is carried out to the introduction of technical background, and facilitate the understanding of those skilled in the art to set forth.Only can not think that technique scheme is conventionally known to one of skill in the art because these schemes have carried out setting forth in the background technology part of the application.
Summary of the invention
The object of the embodiment of the present application is to provide a kind of method and the device of determining bad track data, the error brought to avoid artificial screening and ensure no longer to occur bad track data in follow-up data processing procedure.
What the embodiment of the present application provided a kind ofly determines that the method for bad track data and device are achieved in that
Determine a method for bad track data, comprising:
Normal moveout correction process is carried out to the seismic channel data gathered, obtains the seismic channel data after normal moveout correction;
The first seismic channel data when presetting in window is obtained from the seismic channel data after described normal moveout correction;
Seismic channel data after described normal moveout correction is carried out TEC time error correction, obtains the seismic channel data after TEC time error correction;
Obtain from the seismic channel data after described TEC time error correction described default time window in the second seismic channel data;
Determine the correlativity of described first seismic channel data and described second seismic channel data;
According to the described correlativity determined, determine whether the seismic channel data of described collection is bad track data.
Determine a device for bad track data, comprising:
Normal moveout correction unit, for carrying out normal moveout correction process to the seismic channel data gathered, obtains the seismic channel data after normal moveout correction;
First seismic channel data acquiring unit, the first seismic channel data when presetting for obtaining from the seismic channel data after described normal moveout correction in window;
TEC time error correction unit, for the seismic channel data after described normal moveout correction is carried out TEC time error correction, obtains the seismic channel data after TEC time error correction;
Second seismic channel data acquiring unit, for obtain from the seismic channel data after described TEC time error correction described default time window in the second seismic channel data;
Correlation determination unit, for determining the correlativity of described first seismic channel data and described second seismic channel data;
Bad track data determination unit, for according to the described correlativity determined, determines whether the seismic channel data of described collection is bad track data.
A kind of method and device determining bad track data that the embodiment of the present application provides, by carrying out normal moveout correction and TEC time error correction to the seismic channel data gathered, thus can obtain the first seismic channel data and the second seismic channel data.Further, described first seismic channel data and described second seismic channel data are arranged in the stratum of different depth, by the correlativity between comparison first seismic channel data and the second seismic channel data, thus can determine whether the seismic channel data gathered is bad track data.A kind of method determining bad track data that the embodiment of the present application provides, can evade the error that artificial screening exists, and fundamentally can reject, bad track data to avoid again occurring bad track data in follow-up data processing procedure simultaneously.
With reference to explanation hereinafter and accompanying drawing, disclose in detail the particular implementation of the application, the principle specifying the application can adopted mode.Should be appreciated that, thus the embodiment of the application is not restricted in scope.In the spirit of claims and the scope of clause, the embodiment of the application comprises many changes, amendment and is equal to.
The feature described for a kind of embodiment and/or illustrate can use in one or more other embodiment in same or similar mode, combined with the feature in other embodiment, or substitutes the feature in other embodiment.
Should emphasize, term " comprises/comprises " existence referring to feature, one integral piece, step or assembly when using herein, but does not get rid of the existence or additional of one or more further feature, one integral piece, step or assembly.
Accompanying drawing explanation
Included accompanying drawing is used to provide the further understanding to the embodiment of the present application, which constitutes a part for instructions, for illustrating the embodiment of the application, and comes together to explain the principle of the application with text description.Apparently, the accompanying drawing in the following describes is only some embodiments of the application, for those of ordinary skill in the art, under the prerequisite not paying creative work, can also obtain other accompanying drawing according to these accompanying drawings.In the accompanying drawings:
A kind of method flow diagram determining bad track data that Fig. 1 provides for the embodiment of the present application;
The principle schematic of a kind of normal moveout correction that Fig. 2 provides for the embodiment of the present application;
A kind of functional block diagram determining the device of bad track data that Fig. 3 provides for the embodiment of the present application.
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 other embodiments all obtained under creative work prerequisite, all should belong to the scope of the application's protection.
A kind of method flow diagram determining bad track data that Fig. 1 provides for the embodiment of the present application.Although hereafter describe flow process to comprise the multiple operations occurred with particular order, but should have a clear understanding of, these processes can comprise more or less operation, and these operations can sequentially perform or executed in parallel (such as using parallel processor or multi-thread environment).As shown in Figure 1, described method comprises:
S1: normal moveout correction process is carried out to the seismic channel data gathered, obtains the seismic channel data after normal moveout correction.
In the embodiment of the present application, in order to correct the time difference brought because offset distance is different in the seismic channel data of collection, normal moveout correction process can be carried out to the seismic channel data gathered.Before carrying out normal moveout correction process to the seismic channel data gathered, the embodiment of the present application can adopt the static correction means of this area routine, first carries out static correction to the seismic channel data of described collection.After static correction, the shot point in seismic channel data and geophone station can be moved on reference field.For CRP gather, when formation velocity is identical, the reflection wave of same shot point is in the seismic trace of different offset distance, and its time arriving geophone station is different.Offset distance is larger, and the time of arrival then can be longer.
For the ease of subsequent treatment, the seismic channel data through static correction can be carried out normal moveout correction process in the embodiment of the present application, to be corrected to unified time the reflection interval of same shot point.Particularly, a kind of principle schematic of normal moveout correction that provides for the embodiment of the present application of Fig. 2.Fig. 2 shows exciting of CRP gather and detection process.As shown in Figure 2, S point is shot point, and P point is reflection spot, R point is geophone station, and M point is the vertical projection of P point on stratum, and the vertical range on reflection spot P distance ground is h, offset distance between shot point S and geophone station R is x, from P point to the firing time t of R point is so:
Wherein, v represents formation velocity.
In the embodiment of the present application, the object of normal moveout correction is to be that t reflection interval of the seismic trace of x is corrected to the reflection interval that offset distance is the seismic trace of 0 by offset distance, and namely M point is to the vertical reflection time of P point:
Wherein, T
0for M point is to the vertical reflection time of P point.
Like this, namely the dynamic correction value Δ t of this seismic trace can be expressed as:
Δt＝t-T
0
Eventually pass through abbreviation can obtain:
As can be seen here, the dynamic correction value that each seismic trace is corresponding is relevant to the offset distance x of this seismic trace, and offset distance is larger, and dynamic correction value is just larger.By carrying out the process of normal moveout correction to the seismic channel data gathered, thus can be corrected to unified time the reflection interval of same shot point, obtain the seismic channel data after normal moveout correction, so that subsequent treatment.
S2: obtain the first seismic channel data when presetting in window from the seismic channel data after described normal moveout correction.
In the embodiment of the present application, can adopt for the analysis of bad track data in seismic channel data the method analyzed piecemeal.Particularly, the embodiment of the present application can select one preset time window, described default time window duration can determine according to the duration of described seismic channel data.Such as, described default time window duration can be set to 10% of described seismic channel data duration.In the embodiment of the present application, the duration of described seismic channel data can be 2000ms, so described default time window can be just 0ms to 200ms, when this is default, the duration of window just can be 200ms.Determine described default time window after, just can obtain from the seismic channel data after described normal moveout correction described default time window in the first seismic channel data.Whether described first seismic channel data truly can reflect the attribute of the seismic channel data of described collection, be namely bad track data.
S3: the seismic channel data after described normal moveout correction is carried out TEC time error correction, obtains the seismic channel data after TEC time error correction.
In the embodiment of the present application, in order to determine whether the seismic channel data gathered is bad track data, after the seismic channel data of collection can being converted to shallow-layer data from deep layer data, again the shallow-layer data be converted to and original deep layer data are contrasted, thus whether the seismic channel data determining described collection is bad track data.After carrying out normal moveout correction process to the seismic channel data of described collection, the embodiment of the present application can carry out TEC time error correction to the seismic channel data after described normal moveout correction again, obtains the seismic channel data after TEC time error correction.Particularly, the embodiment of the present application can pre-determine the duration of TEC time error correction.This duration can be considered as the duration that the seismic channel data after described normal moveout correction needs to carry out movement.In the application one preferred embodiment, described TEC time error correction duration can be determined according to the following equation:
h≤L-2l-x
max/v
Wherein, h is described TEC time error correction duration, and L is the duration of the seismic channel data of described collection, l be described default time window duration, x
maxfor the maximum offset that the seismic channel data of described collection is corresponding, v is the speed of seismic wave propagation.
What above formula provided is a value range, according to the demand in practical application scene, the embodiment of the present application can determine a concrete value in this value range, such as 3000ms, thus the duration that this can be utilized to determine corrects duration, and the seismic channel data after described normal moveout correction is converted to second degree of depth from first degree of depth.In the embodiment of the present application, the TEC time error correction of carrying out the seismic channel data after described normal moveout correction can be upwards correct, thus the seismic channel data being positioned at deep layer can be converted to the seismic channel data being positioned at shallow-layer.Therefore, described first degree of depth is greater than described second degree of depth, like this, just the seismic channel data being positioned at described second depth can be defined as the seismic channel data after TEC time error correction.Seismic channel data after this TEC time error correction can, as with reference to data, utilize this reference data just can weigh the validity of the seismic channel data of collection, to determine whether the seismic channel data of described collection is bad track data.
S4: obtain from the seismic channel data after described TEC time error correction described default time window in the second seismic channel data.
After seismic channel data after obtaining described TEC time error correction, just can obtain the seismic channel data corresponding with described first seismic channel data from these data, thus can follow-up comparison process be carried out.Particularly, the embodiment of the present application can obtain from the seismic channel data after described TEC time error correction equally described default time window in seismic channel data, thus obtain the second seismic channel data for reference.Described default time window and step S2 in default time window should be consistent, such as choose in step s 2 default time window be 200ms to 400ms, so now also should obtain the seismic channel data of 200ms to 400ms in the seismic channel data after described TEC time error correction, like this, described first seismic channel data and described second seismic channel data just can one_to_one corresponding, thus facilitate follow-up comparison process.
S5: the correlativity determining described first seismic channel data and described second seismic channel data.
In the embodiment of the present application, after obtaining described first seismic channel data and described second seismic channel data, can be analyzed these two data.The object analyzed is to weigh the similarity degree between described two data.When these two data are relatively more consistent, then illustrate that the seismic channel data of described collection, in TEC time error correction process, significantly change does not occur, thus prove that this seismic channel data is normal seismic channel data; But, if described first seismic channel data differs comparatively large with described second seismic channel data, then illustrate in the process of TEC time error correction, the seismic channel data of described collection there occurs obvious change, then can illustrate that the seismic channel data of described collection exists instability, be bad track data.
In the embodiment of the present application, the correlation of described first seismic channel data and described second seismic channel data can be calculated.Particularly, a sampling parameter can be preset, carry out discrete to described first seismic channel data and described second seismic channel data.Such as, described sampling parameter can be set to 10ms, like this, just can obtain a sampled point from these two groups of data every 10ms, for the duration of 200ms, just can obtain 20 sampled points respectively from described first seismic channel data and described second seismic channel data.But can the data corresponding to these 20 sampled points calculate, to ask for the correlation between described first seismic channel data and described second seismic channel data.Particularly, the formula calculating correlation can be:
Wherein, σ is the correlation between described first seismic channel data and described second seismic channel data, and n is the number of sampled point, x
ifor the data that i-th sampled point in described first seismic channel data is corresponding, y
ifor the data that i-th sampled point in described second seismic channel data is corresponding,
for the mean value of each sampling number certificate in described first seismic channel data,
for the mean value of each sampling number certificate in described second seismic channel data.
In another preferred embodiment of the application, in order to simplify the contrast flow process of two groups of seismic channel data and improve the precision of contrast, the first root-mean-square amplitude of described first seismic channel data and the second RMS amplitude of described second seismic channel data can be calculated.Particularly, described first root-mean-square amplitude and described second root-mean-square amplitude can be calculated by following formula:
Wherein, X
rmsqfor described first RMS amplitude, X
rmssfor described second root-mean-square amplitude, n is the sampled point number of described first seismic channel data or the second seismic channel data, X
qkt () is a kth sampling number certificate in described first seismic channel data, X
skt () is a kth sampling number certificate in described second seismic channel data.
After calculating described first RMS amplitude and described second RMS amplitude, the ratio of described first root-mean-square amplitude and described second root-mean-square amplitude can be calculated, and the ratio of described first root-mean-square amplitude and described second root-mean-square amplitude is defined as the correlativity of described first seismic channel data and described second seismic channel data.
S6: according to the described correlativity determined, determines whether the seismic channel data of described collection is bad track data.
In the embodiment of the present application, after determining the correlativity between described first seismic channel data and described second seismic channel data, just based on this correlativity, can determine whether the seismic channel data of described collection is bad track data.Particularly, the embodiment of the present application can pre-set a condition, and this pre-conditioned critical condition that can be used as screening bad track, when described correlativity does not meet pre-conditioned, is defined as bad track data by the seismic channel data of described collection; When described correlativity meet described pre-conditioned time, the seismic channel data of described collection is defined as normal earthquake track data.
After determining bad track data, number record corresponding for this seismic trace can be got off, and the numbering of the seismic data belonging to this seismic trace can be inquired about, thus can know which seismic channel data in which seismic data is bad track data, thus from described seismic data, bad track data can be removed, to ensure follow-up accuracy of carrying out the process analyzed based on seismic data.
Can being found by above-described embodiment, a kind of method determining bad track data that the embodiment of the present application provides, by carrying out normal moveout correction and TEC time error correction to the seismic channel data gathered, thus the first seismic channel data and the second seismic channel data can be obtained.Further, described first seismic channel data and described second seismic channel data are arranged in the stratum of different depth, by the difference between comparison first seismic channel data and the second seismic channel data, thus can determine whether the seismic channel data gathered is bad track data.A kind of method determining bad track data that the embodiment of the present application provides, can evade the error that artificial screening exists, and fundamentally can reject, bad track data to avoid again occurring bad track data in follow-up data processing procedure simultaneously.
The embodiment of the present application also provides a kind of device determining bad track data.A kind of functional block diagram determining the device of bad track data that Fig. 3 provides for the embodiment of the present application.As shown in Figure 3, described device comprises:
Normal moveout correction unit 100, for carrying out normal moveout correction process to the seismic channel data gathered, obtains the seismic channel data after normal moveout correction;
First seismic channel data acquiring unit 200, the first seismic channel data when presetting for obtaining from the seismic channel data after described normal moveout correction in window;
TEC time error correction unit 300, for the seismic channel data after described normal moveout correction is carried out TEC time error correction, obtains the seismic channel data after TEC time error correction;
Second seismic channel data acquiring unit 400, for obtain from the seismic channel data after described TEC time error correction described default time window in the second seismic channel data;
Correlation determination unit 500, for determining the correlativity of described first seismic channel data and described second seismic channel data;
Bad track data determination unit 600, for according to the described correlativity determined, determines whether the seismic channel data of described collection is bad track data.
In the application one preferred embodiment, described TEC time error correction unit 300 specifically comprises:
Duration determination module, for determining TEC time error correction duration;
Degree of depth modular converter, for utilizing described TEC time error correction duration, is converted to second degree of depth by the seismic channel data after described normal moveout correction from first degree of depth, and described first degree of depth is greater than described second degree of depth;
Determination module, for being defined as the seismic channel data after TEC time error correction by the seismic channel data being positioned at described second depth.
Wherein, described duration determining unit determines TEC time error correction duration according to the following equation:
h≤L-2l-x
max/v
Wherein, h is described TEC time error correction duration, and L is the duration of the seismic channel data of described collection, l be described default time window duration, x
maxfor the maximum offset that the seismic channel data of described collection is corresponding, v is the speed of seismic wave propagation.
In another preferred embodiment of the application, described correlation determination unit 500 specifically comprises:
Root-mean-square amplitude computing module, for the second RMS amplitude of the first root-mean-square amplitude and described second seismic channel data that calculate described first seismic channel data;
Ratio calculation module, for calculating the ratio of described first root-mean-square amplitude and described second root-mean-square amplitude;
Correlation determining module, for being defined as the correlativity of described first seismic channel data and described second seismic channel data by the ratio of described first root-mean-square amplitude and described second root-mean-square amplitude.
It should be noted that, about in above-described embodiment of device, the implementation procedure of each functional module and the computing formula related to all identical with step S1 to S6, just repeat no more here.
Can being found by above-described embodiment, a kind of device determining bad track data that the embodiment of the present application provides, by carrying out normal moveout correction and TEC time error correction to the seismic channel data gathered, thus the first seismic channel data and the second seismic channel data can be obtained.Further, described first seismic channel data and described second seismic channel data are arranged in the stratum of different depth, by the difference between comparison first seismic channel data and the second seismic channel data, thus can determine whether the seismic channel data gathered is bad track data.A kind of method determining bad track data that the embodiment of the present application provides, can evade the error that artificial screening exists, and fundamentally can reject, bad track data to avoid again occurring bad track data in follow-up data processing procedure simultaneously.
In this manual, such as first and second such adjectives only may be used for an element or action and another element or action to distinguish, and without requiring or imply this relation or the order of any reality.When environment allows, should not be construed as one that is confined in only element, parts or step with reference to element or parts or step (s), and can be one or more etc. in element, parts or step.
With the object described, those skilled in the art are supplied to the description of the various embodiments of the application above.It is not intended to is exhaustive or is not intended to the present invention to be limited to single disclosed embodiment.As mentioned above, the various alternative and change of the application will be apparent for above-mentioned technology one of ordinary skill in the art.Therefore, although specifically discuss the embodiment of some alternatives, other embodiment will be apparent, or those skilled in the art relatively easily draw.The application is intended to be included in that of the present invention all that this had discussed substitute, amendment and change, and drops on other embodiment in the spirit and scope of above-mentioned application.
Each embodiment in this instructions all adopts the mode of going forward one by one to describe, between each embodiment identical similar part mutually see, what each embodiment stressed is the difference with other embodiments.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.
The application can be used in numerous general or special purpose computing system environments or configuration.Such as: personal computer, server computer, handheld device or portable set, laptop device, multicomputer system, system, set top box, programmable consumer-elcetronics devices, network PC, small-size computer, mainframe computer, the distributed computing environment comprising above any system or equipment etc. based on microprocessor.
The application can describe in the general context of computer executable instructions, such as program module.Usually, program module comprises the routine, program, object, assembly, data structure etc. that perform particular task or realize particular abstract data type.Also can put into practice the application in a distributed computing environment, in these distributed computing environment, be executed the task by the remote processing devices be connected by communication network.In a distributed computing environment, program module can be arranged in the local and remote computer-readable storage medium comprising memory device.
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 (10)
1. determine a method for bad track data, it is characterized in that, comprising:
Normal moveout correction process is carried out to the seismic channel data gathered, obtains the seismic channel data after normal moveout correction;
The first seismic channel data when presetting in window is obtained from the seismic channel data after described normal moveout correction;
Seismic channel data after described normal moveout correction is carried out TEC time error correction, obtains the seismic channel data after TEC time error correction;
Obtain from the seismic channel data after described TEC time error correction described default time window in the second seismic channel data;
Determine the correlativity of described first seismic channel data and described second seismic channel data;
According to the described correlativity determined, determine whether the seismic channel data of described collection is bad track data.
2. a kind of method determining bad track data as claimed in claim 1, is characterized in that, described seismic channel data after described normal moveout correction is carried out TEC time error correction, obtains the seismic channel data after TEC time error correction and specifically comprises:
Determine TEC time error correction duration;
According to described TEC time error correction duration, the seismic channel data after described normal moveout correction is converted to second degree of depth from first degree of depth, described first degree of depth is greater than described second degree of depth;
The seismic channel data being positioned at described second depth is defined as the seismic channel data after TEC time error correction.
3. a kind of method determining bad track data as claimed in claim 2, is characterized in that, determine TEC time error correction duration according to the following equation:
h≤L-2l-x
max/v
Wherein, h is described TEC time error correction duration, and L is the duration of the seismic channel data of described collection, l be described default time window duration, x
maxfor the maximum offset that the seismic channel data of described collection is corresponding, v is the speed of seismic wave propagation.
4. a kind of method determining bad track data as claimed in claim 1, is characterized in that, describedly determines that the correlativity of described first seismic channel data and described second seismic channel data specifically comprises:
Calculate the first root-mean-square amplitude of described first seismic channel data and the second RMS amplitude of described second seismic channel data;
Calculate the ratio of described first root-mean-square amplitude and described second root-mean-square amplitude;
The ratio of described first root-mean-square amplitude and described second root-mean-square amplitude is defined as the correlativity of described first seismic channel data and described second seismic channel data.
5. a kind of method determining bad track data as claimed in claim 4, is characterized in that, calculate the first root-mean-square amplitude of described first seismic channel data and the second RMS amplitude of described second seismic channel data according to the following equation:
Wherein, X
rmsqfor described first RMS amplitude, X
rmssfor described second root-mean-square amplitude, n is the sampled point number of described first seismic channel data or the second seismic channel data, X
qkt () is a kth sampling number certificate in described first seismic channel data, X
skt () is a kth sampling number certificate in described second seismic channel data.
6. a kind of method determining bad track data as described in claim 1 or 4, is characterized in that, the described described correlativity according to determining, determines whether the seismic channel data of described collection is that bad track data specifically comprise:
When described correlativity does not meet pre-conditioned, the seismic channel data of described collection is defined as bad track data;
When described correlativity meet described pre-conditioned time, the seismic channel data of described collection is defined as normal earthquake track data.
7. determine a device for bad track data, it is characterized in that, comprising:
Normal moveout correction unit, for carrying out normal moveout correction process to the seismic channel data gathered, obtains the seismic channel data after normal moveout correction;
First seismic channel data acquiring unit, the first seismic channel data when presetting for obtaining from the seismic channel data after described normal moveout correction in window;
TEC time error correction unit, for the seismic channel data after described normal moveout correction is carried out TEC time error correction, obtains the seismic channel data after TEC time error correction;
Second seismic channel data acquiring unit, for obtain from the seismic channel data after described TEC time error correction described default time window in the second seismic channel data;
Correlation determination unit, for determining the correlativity of described first seismic channel data and described second seismic channel data;
Bad track data determination unit, for according to the described correlativity determined, determines whether the seismic channel data of described collection is bad track data.
8. a kind of device determining bad track data as claimed in claim 7, it is characterized in that, described TEC time error correction unit specifically comprises:
Duration determination module, for determining TEC time error correction duration;
Degree of depth modular converter, for utilizing described TEC time error correction duration, is converted to second degree of depth by the seismic channel data after described normal moveout correction from first degree of depth, and described first degree of depth is greater than described second degree of depth;
Determination module, for being defined as the seismic channel data after TEC time error correction by the seismic channel data being positioned at described second depth.
9. a kind of device determining bad track data as claimed in claim 8, it is characterized in that, described duration determining unit determines TEC time error correction duration according to the following equation:
h≤L-2l-x
max/y
Wherein, h is described TEC time error correction duration, and L is the duration of the seismic channel data of described collection, l be described default time window duration, x
maxfor the maximum offset that the seismic channel data of described collection is corresponding, v is the speed of seismic wave propagation.
10. a kind of device determining bad track data as claimed in claim 7, it is characterized in that, described correlation determination unit specifically comprises:
Root-mean-square amplitude computing module, for the second RMS amplitude of the first root-mean-square amplitude and described second seismic channel data that calculate described first seismic channel data;
Ratio calculation module, for calculating the ratio of described first root-mean-square amplitude and described second root-mean-square amplitude;
Correlation determining module, for being defined as the correlativity of described first seismic channel data and described second seismic channel data by the ratio of described first root-mean-square amplitude and described second root-mean-square amplitude.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510852243.2A CN105445793B (en) | 2015-11-30 | 2015-11-30 | A kind of method and device for determining bad track data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510852243.2A CN105445793B (en) | 2015-11-30 | 2015-11-30 | A kind of method and device for determining bad track data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105445793A true CN105445793A (en) | 2016-03-30 |
| CN105445793B CN105445793B (en) | 2018-04-06 |
Family
ID=55556205
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510852243.2A Active CN105445793B (en) | 2015-11-30 | 2015-11-30 | A kind of method and device for determining bad track data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105445793B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106772576A (en) * | 2016-12-01 | 2017-05-31 | 中国石油天然气集团公司 | Seismic data processing technique and device |
| CN107193045A (en) * | 2017-07-07 | 2017-09-22 | 中国石油天然气集团公司 | A kind of seismic data processing technique and device |
| CN107884826A (en) * | 2017-12-13 | 2018-04-06 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | The method of dead trace recognition methods and optimization seismic data process in geological data |
| CN109557583A (en) * | 2017-09-26 | 2019-04-02 | 中国石油化工股份有限公司 | A kind of seismic attribute extracting method and system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030031091A1 (en) * | 2001-07-26 | 2003-02-13 | Ki-Young Kim | Method, system, and apparatus for processing seismic data by means of the time-varying optimum offset concept |
| CN101630019A (en) * | 2008-07-16 | 2010-01-20 | 中国石油集团东方地球物理勘探有限责任公司 | Method for correcting predicated large time difference of multiple benthal waves |
| CN104391324A (en) * | 2014-12-03 | 2015-03-04 | 成都理工大学 | Seismic trace set dynamic correction stretching correction pre-processing technology before AVO inversion depending on frequency |
-
2015
- 2015-11-30 CN CN201510852243.2A patent/CN105445793B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030031091A1 (en) * | 2001-07-26 | 2003-02-13 | Ki-Young Kim | Method, system, and apparatus for processing seismic data by means of the time-varying optimum offset concept |
| CN101630019A (en) * | 2008-07-16 | 2010-01-20 | 中国石油集团东方地球物理勘探有限责任公司 | Method for correcting predicated large time difference of multiple benthal waves |
| CN104391324A (en) * | 2014-12-03 | 2015-03-04 | 成都理工大学 | Seismic trace set dynamic correction stretching correction pre-processing technology before AVO inversion depending on frequency |
Non-Patent Citations (2)
| Title |
|---|
| 刘成斋 等: "地震数据体中非正常地震道的自动识别与评价方法", 《油气地质与采收率》 * |
| 张志让 等: "一种基于模型的海上地震资料中废炮废道的检测方法", 《石油地球物理勘探》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106772576A (en) * | 2016-12-01 | 2017-05-31 | 中国石油天然气集团公司 | Seismic data processing technique and device |
| CN107193045A (en) * | 2017-07-07 | 2017-09-22 | 中国石油天然气集团公司 | A kind of seismic data processing technique and device |
| CN107193045B (en) * | 2017-07-07 | 2019-02-15 | 中国石油天然气集团公司 | A kind of seismic data processing technique and device |
| CN109557583A (en) * | 2017-09-26 | 2019-04-02 | 中国石油化工股份有限公司 | A kind of seismic attribute extracting method and system |
| CN107884826A (en) * | 2017-12-13 | 2018-04-06 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | The method of dead trace recognition methods and optimization seismic data process in geological data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105445793B (en) | 2018-04-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Grigoli et al. | Automated seismic event location by travel‐time stacking: An application to mining induced seismicity | |
| Li et al. | Joint microseismic location and anisotropic tomography using differential arrival times and differential backazimuths | |
| Zhang et al. | Cross-well seismic waveform tomography for monitoring CO2 injection: a case study from the Ketzin Site, Germany | |
| Restivo et al. | Teleseismic shear wave splitting measurements in noisyenvironments | |
| EP2745146B1 (en) | System and method for subsurface characterization including uncertainty estimation | |
| RU2570221C2 (en) | Determining position of geologic horizon relative to manifestation of seismic pulse in seismic data | |
| CN104375188A (en) | Seismic wave transmission attenuation compensation method and device | |
| Song et al. | Microseismic event estimation and velocity analysis based on a source-focusing function | |
| CN105445793A (en) | Method and apparatus for determining bad trace data | |
| US20160161619A1 (en) | Systems and methods for aligning a monitor seismic survey with a baseline seismic survey | |
| CN104502996A (en) | Density curve correction method and system | |
| CN103592682A (en) | Seismic wave field energy balancing method and processing device | |
| CN104375178A (en) | Carbonate rock fracture-cave reservoir prediction method and device | |
| Li et al. | Chaotic system detection of weak seismic signals | |
| Dando et al. | Complexity in microseismic phase identification: full waveform modelling, traveltime computations and implications for event locations within the Groningen gas field | |
| Purcell et al. | Integrating velocity measurements in a reservoir rock sample from the SACROC unit with an AVO proxy for subsurface supercritical CO 2 | |
| CN106707345B (en) | Method and device for identifying lithology of angle elastic parameter | |
| CN106501851B (en) | A kind of optimum methods of seismic attributes and device | |
| CN104536046A (en) | Focus excitation signal consistency evaluating method based on earthquake records | |
| CN117233859A (en) | Method and system for monitoring carbon sequestration based on seismic exploration and electromagnetic exploration technology | |
| CN104570115A (en) | Surface wave attenuating method and device | |
| Hendriyana et al. | AIC-based diffraction stacking for local earthquake locations at the Sumatran Fault (Indonesia) | |
| Riazi et al. | Application of focal-time analysis for improved induced seismicity depth control: A case study from the Montney Formation, British Columbia, Canada | |
| US20140350857A1 (en) | Method Of Mapping A Subterranean Formation Based Upon Wellbore Position And Seismic Data And Related System | |
| Dai et al. | Study of an Automatic Picking Method for Multimode Dispersion Curves of Surface Waves Based on an Improved U-Net |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |