CN105445793B - A kind of method and device for determining bad track data - Google Patents
A kind of method and device for determining bad track data Download PDFInfo
- Publication number
- CN105445793B CN105445793B CN201510852243.2A CN201510852243A CN105445793B CN 105445793 B CN105445793 B CN 105445793B CN 201510852243 A CN201510852243 A CN 201510852243A CN 105445793 B CN105445793 B CN 105445793B
- Authority
- CN
- China
- Prior art keywords
- seismic channel
- channel data
- data
- mrow
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000012937 correction Methods 0.000 claims abstract description 115
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000001568 sexual effect Effects 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 14
- 238000012216 screening Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 239000004744 fabric Substances 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
- 230000000750 progressive effect Effects 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 application provides a kind of method and device for determining bad track data, wherein, methods described includes：Correction process of taking action is entered to the seismic channel data of collection, obtains the seismic channel data after dynamic correction；The first seismic channel data when obtaining default from the seismic channel data after the dynamic correction in window；Seismic channel data after the dynamic correction is subjected to TEC time error correction, obtains the seismic channel data after TEC time error correction；The second seismic channel data when obtaining described default from the seismic channel data after the TEC time error correction in window；Determine the correlation of first seismic channel data and second seismic channel data；According to the correlation of determination, whether the seismic channel data for determining the collection is bad track data.A kind of method and device for determination bad track data that the embodiment of the present application provides, error that artificial screening is brought can be avoided and ensure no longer occur bad track data in follow-up data processing procedure.
Description
Technical field
The application is related to seismic data processing technology field, more particularly to a kind of method and device for determining bad track data.
Background technology
The process handled geological data can often be divided into following three phases：
1) first stage, earthquake data acquisition stage, the working contents in this stage are：In survey area, cloth
Put two-dimentional or three-dimensional survey line；Using dynamite source or controlled source earthquake-wave-exciting, wherein, where dynamite source or controlled source
Point is properly termed as shot point；Multiple wave detectors are equally spaced along survey line to receive seismic signal, wave detector or geophone group
Quantity can be set as needed, arrangement mode can be different between each geophone group, such as split shooting, end points
Arrangement etc. is blown out, can also be set as needed；After wave detector receives seismic signal, adopted so that constant duration is discrete
Sample geological data, and be recorded in digital form on tape.
2) second stage, seismic data processing stage, the working contents in this stage are：Managed with seimic wave propagation
Based on, using computer and Seismic data processing software, the geological data of field acquisition in the working process first stage will
Geological data becomes the seismic cross-section that can reflect subsurface structure and can reflect the seismic amplitude of subsurface lithologic change, frequency
The information profile such as rate and spread speed.
3) three phases, seismic data interpretation stage, the working contents in this stage are：Analysis interpretation second-order
The data such as the seismic cross-section in section obtained by seismic data process, it is theoretical according to petroleum geology principle and seimic wave propagation,
The construction of subterranean strata is determined, so as to find out the Favorable Zones of oily and propose drilling well position.
By the process of above-mentioned seismic prospecting it can be found that earthquake data acquisition is the basis of latter earthquake exploration steps,
The quality of earthquake data acquisition will directly affect the result of final seismic prospecting.In addition, earthquake data acquisition parameter is to determine
An important factor for earthquake data quality.Therefore, how rejecting bad track data from the geological data of collection will seem particularly significant.
The mode of the bad track data in identification geological data mainly has two kinds, respectively manual identified and machine pressure at present
System.However, the drawbacks of certain all be present in both approaches, wherein, manual identified often introduces excessive human error, so as to
So that the identification of bad track data is not thorough enough, eventually resulting in still has part bad track parameter to retain；Machine pressing can not be real
Bad track data are removed from geological data, but bad track data are suppressed in the seismic data, to reduce it in data
Inaccuracy in processing procedure, increase and data processing form however as subsequent processing steps it is changeable, it is downtrodden
Bad track data often reappear, so as to influence the accuracy of follow-up data processing.
It should be noted that the introduction to technical background above be intended merely to the convenient technical scheme to the application carry out it is clear,
Complete explanation, and facilitate the understanding of those skilled in the art and illustrate.Can not merely because these schemes the application's
Background section is set forth and thinks that above-mentioned technical proposal is known to those skilled in the art.
The content of the invention
The purpose of the embodiment of the present application is to provide a kind of method and device for determining bad track data, to avoid artificial screening
The error brought and ensure no longer occur bad track data in follow-up data processing procedure.
What a kind of method and device for determination bad track data that the embodiment of the present application provides was realized in：
A kind of method for determining bad track data, including：
Correction process of taking action is entered to the seismic channel data of collection, obtains the seismic channel data after dynamic correction；
The first seismic channel data when obtaining default from the seismic channel data after the dynamic correction in window；
Seismic channel data after the dynamic correction is subjected to TEC time error correction, obtains the seismic channel data after TEC time error correction；
The second seismic channel data when obtaining described default from the seismic channel data after the TEC time error correction in window；
Determine the correlation of first seismic channel data and second seismic channel data；
According to the correlation of determination, whether the seismic channel data for determining the collection is bad track data.
A kind of device for determining bad track data, including：
Dynamic correction unit, for entering correction process of taking action to the seismic channel data of collection, obtain the seismic channel after dynamic correction
Data；
First seismic channel data acquiring unit, when being preset for being obtained from the seismic channel data after the dynamic correction in window
The first seismic channel data；
TEC time error correction unit, for the seismic channel data after the dynamic correction to be carried out into TEC time error correction, obtain TEC time error correction
Seismic channel data afterwards；
Second seismic channel data acquiring unit, it is described default for being obtained from the seismic channel data after the TEC time error correction
When window in the second seismic channel data；
Correlation determination unit, for determining that first seismic channel data is related to second seismic channel data
Property；
Bad track data determination unit, for the correlation according to determination, determining the seismic channel data of the collection is
No is bad track data.
A kind of method and device for determination bad track data that the embodiment of the present application provides, passes through the seismic channel data to collection
Enter action correction and TEC time error correction, so as to obtain the first seismic channel data and the second seismic channel data.Further, institute
State the first seismic channel data and second seismic channel data is located in the stratum of different depth, by comparing the first seismic channel
Correlation between data and the second seismic channel data, whether the seismic channel data that may thereby determine that collection is bad track number
According to.A kind of method for determination bad track data that the embodiment of the present application provides, can evade error existing for artificial screening, while energy
It is enough fundamentally to reject bad track data, to avoid occurring bad track data again in follow-up data processing procedure.
With reference to following explanation and accompanying drawing, the particular implementation of the application is disclose in detail, specifies the original of the application
Reason can be in a manner of adopted.It should be understood that presently filed embodiment is not so limited in scope.In appended power
In the range of the spirit and terms that profit requires, presently filed embodiment includes many changes, modifications and is equal.
The feature for describing and/or showing for a kind of embodiment can be in a manner of same or similar one or more
Used in individual other embodiment, it is combined with the feature in other embodiment, or substitute the feature in other embodiment.
It should be emphasized that term "comprises/comprising" refers to the presence of feature, one integral piece, step or component when being used herein, but simultaneously
It is not excluded for the presence or additional of one or more further features, one integral piece, step or component.
Brief description of the drawings
Included accompanying drawing is used for providing being further understood from the embodiment of the present application, which constitutes one of specification
Point, for illustrating presently filed embodiment, and come together with word description to explain the principle of the application.Under it should be evident that
Accompanying drawing in the description of face is only some embodiments of the present application, for those of ordinary skill in the art, is not paying wound
On the premise of the property made is laborious, other accompanying drawings can also be obtained according to these accompanying drawings.In the accompanying drawings：
Fig. 1 is a kind of method flow diagram for determination bad track data that the embodiment of the present application provides；
Fig. 2 is a kind of principle schematic for dynamic correction that the embodiment of the present application provides；
Fig. 3 is a kind of functional block diagram of the device for determination bad track data that the embodiment of the present application provides.
Embodiment
In order that those skilled in the art more fully understand the technical scheme in the application, it is real below in conjunction with the application
The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described implementation
Example only some embodiments of the present application, rather than whole embodiments.It is common based on the embodiment in the application, this area
All other embodiment that technical staff is obtained under the premise of creative work is not made, it should all belong to the application protection
Scope.
Fig. 1 is a kind of method flow diagram for determination bad track data that the embodiment of the present application provides.Although flow is described below
Including the multiple operations occurred with particular order, but it should be clearly understood that these processes can include more or less behaviour
Make, these operations sequentially can be performed or performed parallel (such as using parallel processor or multi-thread environment).As shown in figure 1,
Methods described includes：
S1：Correction process of taking action is entered to the seismic channel data of collection, obtains the seismic channel data after dynamic correction.
When in the embodiment of the present application, in order to correct in the seismic channel data of collection due to caused by offset distance difference
Difference, correction process of taking action can be entered to the seismic channel data of collection.The seismic channel data to collection enter to take action correction process it
Before, the embodiment of the present application can use the conventional static correction means in this area, the seismic channel data of the collection first be carried out quiet
Correction.After by static correction, shot point and geophone station in seismic channel data can be moved into reference plane.With anti-altogether
Exemplified by exit point trace gather, in the case of formation velocity identical, the seismic channel of the back wave of same shot point in different offset distances
In, the time of its arrival geophone station is different.Offset distance is bigger, and the time of arrival then can be longer.
For the ease of subsequent treatment, the seismic channel data Jing Guo static correction can be entered school of taking action in the embodiment of the present application
Positive processing, a unified time is corrected to by the reflection interval of same shot point.Specifically, Fig. 2 is implemented for the application
A kind of principle schematic for dynamic correction that example provides.Fig. 2 show CRP gather excite and detection process.Such as Fig. 2 institutes
Show, S points are shot point, and P points are pip, and R points are geophone station, and M points are upright projection of the P points on stratum, and pip P is apart from ground
The vertical range in face is h, and the offset distance between shot point S and geophone station R is x, then the firing time t from P points to R points is：
Wherein, v represents formation velocity.
In the embodiment of the present application, the purpose for moving correction is the reflection interval t for the seismic channel that offset distance is x being corrected to
Offset distance be 0 seismic channel reflection interval, i.e. vertical reflection interval of the M points to P points：
Wherein, T0Vertical reflection interval for M points to P points.
So, the dynamic correction value Δ t of the seismic channel can be expressed as：
Δ t=t-T0
Eventually passing through abbreviation can obtain：
As can be seen here, dynamic correction value corresponding to each seismic channel is related to the offset distance x of the seismic channel, and offset distance is bigger,
Dynamic correction value is just bigger.By entering the processing of action correction to the seismic channel data of collection, so as to by same shot point
Reflection interval be corrected to a unified time, the seismic channel data after dynamic correction is obtained, so as to subsequent treatment.
S2：The first seismic channel data when obtaining default from the seismic channel data after the dynamic correction in window.
In the embodiment of the present application, the analysis for bad track data in seismic channel data can use the side analyzed paragraph by paragraph
Method.Specifically, the embodiment of the present application can select one it is default when window, the duration of window can be according to the seismic channel when described default
The durations of data determines.For example, the duration of window could be arranged to the 10% of the seismic channel data duration when described default.
In the embodiment of the present application, the duration of the seismic channel data can be 2000ms, then window can think 0ms extremely when described default
200ms, the duration of window can think 200ms when this is default., can be with from the dynamic correction when being determined described default after window
The first seismic channel data when obtaining described default in seismic channel data afterwards in window.First seismic channel data can be true
Reflect the attribute of the seismic channel data of the collection, i.e., whether be bad track data.
S3：Seismic channel data after the dynamic correction is subjected to TEC time error correction, obtains the seismic channel data after TEC time error correction.
In the embodiment of the present application, can be by collection in order to determine whether the seismic channel data of collection is bad track data
After seismic channel data is from deep layer data conversion to shallow-layer data, then the shallow-layer data being converted to and original deep layer data are entered
Row contrast, so that it is determined that whether the seismic channel data of the collection is bad track data.Enter in the seismic channel data to the collection
Take action after correction process, the embodiment of the present application can carry out TEC time error correction to the seismic channel data after the dynamic correction again, obtain
Seismic channel data after TEC time error correction.Specifically, the embodiment of the present application can predefine the duration of TEC time error correction.The duration can
The duration moved with being considered as the seismic channel data after the dynamic correction to need.It is described in the preferred embodiment of the application one
TEC time error correction duration can be determined according to the following equation：
h≤L-2l-xmax/v
Wherein, h is the TEC time error correction duration, and L is the duration of the seismic channel data of the collection, when l is described default
The duration of window, xmaxFor maximum offset corresponding to the seismic channel data of the collection, v is the speed of seimic wave propagation.
What above formula provided is a value range, and according to the demand in practical application scene, the embodiment of the present application can be
A specific value, such as 3000ms are determined in the value range, the duration so as to be determined using this corrects duration, by institute
The seismic channel data after correction is stated to change to the second depth from the first depth.In the embodiment of the present application, to the dynamic correction
The TEC time error correction that seismic channel data afterwards is carried out can be corrected upwards, so as to be changed positioned at the seismic channel data of deep layer
To the seismic channel data for being located at shallow-layer.Therefore, first depth is more than second depth, so, can be will be located at institute
The seismic channel data for stating the second depth is defined as the seismic channel data after TEC time error correction.Seismic channel data after the TEC time error correction
Reference data can be used as, can be to weigh the validity of the seismic channel data of collection, with described in determination using the reference data
Whether the seismic channel data of collection is bad track data.
S4：The second seismic channel data when obtaining described default from the seismic channel data after the TEC time error correction in window.
, can be with the acquisition from the data and first earthquake after the seismic channel data after obtaining the TEC time error correction
The corresponding seismic channel data of track data, so as to carry out follow-up comparison process.Specifically, the embodiment of the present application equally may be used
Seismic channel data in window during obtaining described default from the seismic channel data after the TEC time error correction, so as to obtain for reference
Second seismic channel data.Window should be consistent with window during default in step S2 when described default, such as be selected in step s 2
Take it is default when window be 200ms to 400ms, then now should also be as obtaining in the seismic channel data after the TEC time error correction
200ms to 400ms seismic channel data, so, first seismic channel data and second seismic channel data just can one
One correspondence, so as to facilitate follow-up comparison process.
S5：Determine the correlation of first seismic channel data and second seismic channel data.
In the embodiment of the present application, can after first seismic channel data and second seismic channel data is obtained
So that the two data to be analyzed.The purpose of analysis is to weigh the similarity degree between described two data.When this
Two data are more consistent, then illustrate that the seismic channel data of the collection becomes during TEC time error correction without generation is obvious
Change, so as to prove that the seismic channel data is normal seismic channel data；If however, first seismic channel data and described the
Two seismic channel datas difference is larger, then explanation is during TEC time error correction, and there occurs obvious for the seismic channel data of the collection
Change, then can illustrate that the seismic channel data of the collection has unstability, be bad track data.
In the embodiment of the present application, the phase of first seismic channel data and second seismic channel data can be calculated
Pass is worth.Specifically, a sampling parameter can be preset, to first seismic channel data and second seismic channel data
Carry out discrete., so, can be to be obtained from this two groups of data every 10ms for example, the sampling parameter could be arranged to 10ms
A sampled point is taken, can be with respectively from first seismic channel data and second ground for 200ms duration
Shake and 20 sampled points are obtained in track data.But data corresponding to this 20 sampled points can be calculated, it is described to ask for
Correlation between first seismic channel data and second seismic channel data.Specifically, calculating the formula of correlation can be：
Wherein, correlations of the σ between first seismic channel data and second seismic channel data, n are sampled point
Number, xiFor data, y corresponding to ith sample point in first seismic channel dataiFor in second seismic channel data
Data corresponding to ith sample point,For the average value of each sample point data in first seismic channel data,To be described
The average value of each sample point data in second seismic channel data.
In another preferred embodiment of the application, in order to simplify the contrast flow of two groups of seismic channel datas and improve contrast
Precision, can calculate first seismic channel data the first root-mean-square amplitude and second seismic channel data second
RMS amplitude.Specifically, first root-mean-square amplitude and second root-mean-square amplitude can be counted by following formula
Calculate：
Wherein, XrmsqFor first RMS amplitude, XrmssFor second root-mean-square amplitude, n is first earthquake
The sampled point number of track data or the second seismic channel data, Xqk(t) it is k-th of sampling number in first seismic channel data
According to Xsk(t) it is k-th of sample point data in second seismic channel data.
After first RMS amplitude and second RMS amplitude is calculated, described first can be calculated
The ratio of root-mean-square amplitude and second root-mean-square amplitude, and by first root-mean-square amplitude and the second root mean square width
The ratio of value is defined as the correlation of first seismic channel data and second seismic channel data.
S6：According to the correlation of determination, whether the seismic channel data for determining the collection is bad track data.
In the embodiment of the present application, when determining between first seismic channel data and second seismic channel data
, can be so that based on the correlation, whether the seismic channel data for determining the collection be bad track data after correlation.Specifically, originally
Application embodiment can pre-set a condition, and the preparatory condition can be used as screening the critical condition of bad track, when described
When correlation is unsatisfactory for preparatory condition, the seismic channel data of the collection is defined as bad track data；When the related sexual satisfaction
During the preparatory condition, the seismic channel data of the collection is defined as normal earthquake track data.
After bad track data are determined, number record corresponding to the seismic channel can be got off, and the ground can be inquired about
The numbering of the seismic data belonging to road is shaken, so as to know which of which seismic data seismic channel data is bad track number
According to so as to remove bad track data from the seismic data, to ensure the mistake subsequently analyzed based on seismic data
The accuracy of journey.
By above-described embodiment it can be found that the embodiment of the present application provide a kind of determination bad track data method, by right
The seismic channel data of collection enters action correction and TEC time error correction, so as to obtain the first seismic channel data and the second seismic channel
Data.Further, first seismic channel data and second seismic channel data are located in the stratum of different depth, lead to
The difference crossed between the first seismic channel data of comparison and the second seismic channel data, may thereby determine that the seismic channel data of collection
Whether it is bad track data.A kind of method for determination bad track data that the embodiment of the present application provides, artificial screening presence can be evaded
Error, while fundamentally bad track data can be rejected, to avoid occurring bad track again in follow-up data processing procedure
Data.
The embodiment of the present application also provides a kind of device for determining bad track data.Fig. 3 is one kind that the embodiment of the present application provides
Determine the functional block diagram of the device of bad track data.As shown in figure 3, described device includes：
Dynamic correction unit 100, for entering correction process of taking action to the seismic channel data of collection, obtain the earthquake after dynamic correction
Track data；
First seismic channel data acquiring unit 200, when being preset for being obtained from the seismic channel data after the dynamic correction
The first seismic channel data in window；
TEC time error correction unit 300, for the seismic channel data after the dynamic correction to be carried out into TEC time error correction, obtain arrival time difference school
Seismic channel data after just；
Second seismic channel data acquiring unit 400, described in being obtained from the seismic channel data after the TEC time error correction
The second seismic channel data when default in window；
Correlation determination unit 500, for determining the phase of first seismic channel data and second seismic channel data
Guan Xing；
Bad track data determination unit 600, for the correlation according to determination, determine the seismic channel data of the collection
Whether it is bad track data.
In the preferred embodiment of the application one, the TEC time error correction unit 300 specifically includes：
Duration determining module, for determining TEC time error correction duration；
Depth modular converter, for utilizing the TEC time error correction duration, by the seismic channel data after the dynamic correction from the
One depth, which is changed to the second depth, first depth, is more than second depth；
Determining module, for by the seismic channel after being defined as TEC time error correction positioned at the seismic channel data of second depth
Data.
Wherein, the duration determining unit determines TEC time error correction duration according to the following equation：
h≤L-2l-xmax/v
Wherein, h is the TEC time error correction duration, and L is the duration of the seismic channel data of the collection, when l is described default
The duration of window, xmaxFor maximum offset corresponding to the seismic channel data of the collection, v is the speed of seimic wave propagation.
In another preferred embodiment of the application, the correlation determination unit 500 specifically includes：
Root-mean-square amplitude computing module, for calculating the first root-mean-square amplitude of first seismic channel data and described
Second RMS amplitude of the second seismic channel data；
Ratio calculation module, for calculating the ratio of first root-mean-square amplitude and second root-mean-square amplitude；
Correlation determining module, for the ratio of first root-mean-square amplitude and second root-mean-square amplitude to be determined
For first seismic channel data and the correlation of second seismic channel data.
It should be noted that in above-described embodiment of device, the implementation process of each functional module and it is related to
Calculation formula with step S1 into S6 it is identical, just repeat no more here.
By above-described embodiment it can be found that the embodiment of the present application provide a kind of determination bad track data device, by right
The seismic channel data of collection enters action correction and TEC time error correction, so as to obtain the first seismic channel data and the second seismic channel
Data.Further, first seismic channel data and second seismic channel data are located in the stratum of different depth, lead to
The difference crossed between the first seismic channel data of comparison and the second seismic channel data, may thereby determine that the seismic channel data of collection
Whether it is bad track data.A kind of method for determination bad track data that the embodiment of the present application provides, artificial screening presence can be evaded
Error, while fundamentally bad track data can be rejected, to avoid occurring bad track again in follow-up data processing procedure
Data.
In this manual, such as first and second adjective can be only used for by an element or action with it is another
One element or action make a distinction, without requiring or implying any actual this relation or order.In the feelings that environment allows
Under condition, be should not be interpreted as limited to reference to element or part or step (s) in only element, part or step one, and can
To be one or more of element, part or step etc..
The description to the various embodiments of the application is supplied to those skilled in the art with the purpose described above.It is not
Be intended to exhaustion or it is not intended to and limits the invention to single disclosed embodiment.As described above, the application's is various
Substitute and change will be apparent for above-mentioned technology one of ordinary skill in the art.Therefore, although specifically begging for
Some alternative embodiments have been discussed, but other embodiment will be apparent, or those skilled in the art are relative
Easily draw.The application is intended to be included in this all replacement of the invention discussed, modification and change, and falls
Other embodiment in the spirit and scope of above-mentioned application.
Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment
Divide mutually referring to what each embodiment stressed is the difference with other embodiment.It is real especially for system
For applying example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method
Part explanation.
The application can be used in numerous general or special purpose computing system environments or configuration.Such as：Personal computer, clothes
Business device computer, handheld device or portable set, laptop device, multicomputer system, the system based on microprocessor, put
Top box, programmable consumer-elcetronics devices, network PC, minicom, mainframe computer including any of the above system or equipment
DCE etc..
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these DCEs, by
Task is performed and connected remote processing devices by communication network.In a distributed computing environment, program module can be with
In the local and remote computer-readable storage medium including storage device.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and
Change is without departing from spirit herein, it is desirable to which appended claim includes these deformations and changed without departing from the application's
Spirit.
Claims (10)
- A kind of 1. method for determining bad track data, it is characterised in that including：Correction process of taking action is entered to the seismic channel data of collection, obtains the seismic channel data after dynamic correction；The first seismic channel data when obtaining default from the seismic channel data after the dynamic correction in window；Seismic channel data after the dynamic correction is subjected to TEC time error correction, obtains the seismic channel data after TEC time error correction；The second seismic channel data when obtaining described default from the seismic channel data after the TEC time error correction in window；Determine the correlation of first seismic channel data and second seismic channel data；According to the correlation of determination, whether the seismic channel data for determining the collection is bad track data.
- 2. a kind of method for determining bad track data as claimed in claim 1, it is characterised in that described by after the dynamic correction Seismic channel data carries out TEC time error correction, obtains the seismic channel data after TEC time error correction and specifically includes：Determine TEC time error correction duration；According to the TEC time error correction duration, the seismic channel data after the dynamic correction is changed to the second depth from the first depth, First depth is more than second depth；By the seismic channel data after being defined as TEC time error correction positioned at the seismic channel data of second depth.
- 3. a kind of method for determining bad track data as claimed in claim 2, it is characterised in that determine the time difference according to the following equation Correct duration：h≤L-2l-xmax/vWherein, h is the TEC time error correction duration, and L is the duration of the seismic channel data of the collection, window when l is described default Duration, xmaxFor maximum offset corresponding to the seismic channel data of the collection, v is the speed of seimic wave propagation.
- 4. a kind of method for determining bad track data as claimed in claim 1, it is characterised in that described to determine first earthquake The correlation of track data and second seismic channel data specifically includes：Calculate the first root-mean-square amplitude of first seismic channel data and the second root mean square of second seismic channel data Amplitude；Calculate the ratio of first root-mean-square amplitude and second root-mean-square amplitude；By the ratio of first root-mean-square amplitude and second root-mean-square amplitude be defined as first seismic channel data and The correlation of second seismic channel data.
- 5. a kind of method for determining bad track data as claimed in claim 4, it is characterised in that described in calculating according to the following equation Second RMS amplitude of the first root-mean-square amplitude of the first seismic channel data and second seismic channel data：<mrow> <msub> <mi>X</mi> <mrow> <mi>r</mi> <mi>m</mi> <mi>s</mi> <mi>q</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <mfrac> <mn>1</mn> <mi>n</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>X</mi> <mrow> <mi>q</mi> <mi>k</mi> </mrow> </msub> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow><mrow> <msub> <mi>X</mi> <mrow> <mi>r</mi> <mi>m</mi> <mi>s</mi> <mi>s</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <mfrac> <mn>1</mn> <mi>n</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>X</mi> <mrow> <mi>s</mi> <mi>k</mi> </mrow> </msub> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow>Wherein, XrmsqFor first RMS amplitude, XrmssFor second root-mean-square amplitude, n is the first seismic channel number According to or the second seismic channel data sampled point number, Xqk(t) it is k-th of sample point data in first seismic channel data, Xsk(t) it is k-th of sample point data in second seismic channel data, t represents the time.
- A kind of 6. method of determination bad track data as described in claim 1 or 4, it is characterised in that the institute according to determination Correlation is stated, whether the seismic channel data for determining the collection is that bad track data specifically include：When the correlation is unsatisfactory for preparatory condition, the seismic channel data of the collection is defined as bad track data；When preparatory condition described in the related sexual satisfaction, the seismic channel data of the collection is defined as normal earthquake road number According to.
- A kind of 7. device for determining bad track data, it is characterised in that including：Dynamic correction unit, for entering correction process of taking action to the seismic channel data of collection, obtain the seismic channel data after dynamic correction；First seismic channel data acquiring unit, during for obtaining default from the seismic channel data after the dynamic correction in window the One seismic channel data；TEC time error correction unit, for the seismic channel data after the dynamic correction to be carried out into TEC time error correction, after obtaining TEC time error correction Seismic channel data；Second seismic channel data acquiring unit, window during for obtaining described default from the seismic channel data after the TEC time error correction The second interior seismic channel data；Correlation determination unit, for determining the correlation of first seismic channel data and second seismic channel data；Bad track data determination unit, for the correlation according to determination, determine the collection seismic channel data whether be Bad track data.
- 8. a kind of device for determining bad track data as claimed in claim 7, it is characterised in that the TEC time error correction unit is specific Including：Duration determining module, for determining TEC time error correction duration；Depth modular converter, it is for utilizing the TEC time error correction duration, the seismic channel data after the dynamic correction is deep from first Degree conversion to the second depth, first depth is more than second depth；Determining module, for by the seismic channel number after being defined as TEC time error correction positioned at the seismic channel data of second depth According to.
- A kind of 9. device for determining bad track data as claimed in claim 8, it is characterised in that the duration determining unit according to Following formula determine TEC time error correction duration：h≤L-2l-xmax/vWherein, h is the TEC time error correction duration, and L is the duration of the seismic channel data of the collection, window when l is described default Duration, xmaxFor maximum offset corresponding to the seismic channel data of the collection, v is the speed of seimic wave propagation.
- A kind of 10. device for determining bad track data as claimed in claim 7, it is characterised in that the correlation determination unit Specifically include：Root-mean-square amplitude computing module, for calculating the first root-mean-square amplitude and described second of first seismic channel data Second RMS amplitude of seismic channel data；Ratio calculation module, for calculating the ratio of first root-mean-square amplitude and second root-mean-square amplitude；Correlation determining module, for the ratio of first root-mean-square amplitude and second root-mean-square amplitude to be defined as into institute State the correlation of the first seismic channel data and second seismic channel data.
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 CN105445793A (en) | 2016-03-30 |
| CN105445793B true 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) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106772576B (en) * | 2016-12-01 | 2018-11-16 | 中国石油天然气集团公司 | Seismic data processing technique and device |
| CN107193045B (en) * | 2017-07-07 | 2019-02-15 | 中国石油天然气集团公司 | A kind of seismic data processing technique and device |
| CN109557583B (en) * | 2017-09-26 | 2020-12-01 | 中国石油化工股份有限公司 | Seismic attribute extraction 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 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100412097B1 (en) * | 2001-07-26 | 2003-12-24 | 주식회사 산하이앤씨 | Method, system and apparatus for processing seismic data using the time-varying optimum offset concept |
| CN101630019B (en) * | 2008-07-16 | 2011-12-07 | 中国石油集团东方地球物理勘探有限责任公司 | 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
Also Published As
| Publication number | Publication date |
|---|---|
| CN105445793A (en) | 2016-03-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105445793B (en) | A kind of method and device for determining bad track data | |
| CN103842853B (en) | Time shift geological data is corrected for coating and record effect | |
| CN104237940B (en) | A kind of diffraction wave imaging method based on dynamic characteristic and device | |
| CN104155691B (en) | Converted wave anisotropic velocity analysis method and device | |
| CN104199103B (en) | Static correction method and static correction device for converted wave | |
| CN106226818A (en) | Seismic data processing technique and device | |
| CN107678063B (en) | A kind of multi-component converted wave crack prediction method based on Rank correlation | |
| CN107807382A (en) | A kind of high density surface wave exploration data acquiring and processing method | |
| CN104570110A (en) | Multi-component data joint speed analysis method based on longitudinal and horizontal wave matching | |
| CN106777060A (en) | Node type seismic detector data file management system and method | |
| CN103869362A (en) | Method and equipment for obtaining body curvature | |
| CN105093297A (en) | Quality control method for minor tremor positioning precision | |
| CN105005075B (en) | Multi-wave matching method based on seismic frequency information | |
| CN105986819A (en) | Method and device used for automatic processing and comprehensive interpretation of logging information | |
| CN104570115B (en) | A kind of surface wave attenuation method and device | |
| Wang et al. | Identification Protocols for Horizontal‐to‐Vertical Spectral Ratio Peaks | |
| CN106501851B (en) | A kind of optimum methods of seismic attributes and device | |
| CN105445788B (en) | A kind of automatic means of interpretation of normal-moveout spectrum based on model and global optimizing | |
| CN107229069B (en) | The method that a kind of pair of common receiver data trace gather carries out velocity analysis | |
| CN104536046A (en) | Focus excitation signal consistency evaluating method based on earthquake records | |
| CN106707345A (en) | Angle elasticity parameter lithology identification method and device | |
| CN107991705B (en) | Log curve correction method and apparatus based on Two-dimensional Statistical feature | |
| US20230072038A1 (en) | Method and apparatus for removing tube wave interference from optical fiber acoustic wave sensing seismic data | |
| Copley | Seismic characterization of Niobrara fluid and rock properties: A 4D study and multicomponent (3C) analysis | |
| CN104459788B (en) | Method and system for determining seismic data collecting parameters |
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 |