CN104570116A  Geological marker bedbased time difference analyzing and correcting method  Google Patents
Geological marker bedbased time difference analyzing and correcting method Download PDFInfo
 Publication number
 CN104570116A CN104570116A CN201310520070.5A CN201310520070A CN104570116A CN 104570116 A CN104570116 A CN 104570116A CN 201310520070 A CN201310520070 A CN 201310520070A CN 104570116 A CN104570116 A CN 104570116A
 Authority
 CN
 China
 Prior art keywords
 time
 data
 layer
 time difference
 window
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Pending
Links
 239000003550 marker Substances 0.000 title abstract 6
 238000004364 calculation method Methods 0.000 claims abstract description 13
 238000000034 method Methods 0.000 claims description 37
 238000004458 analytical method Methods 0.000 claims description 24
 238000001914 filtration Methods 0.000 claims description 19
 230000003068 static Effects 0.000 claims description 15
 238000003384 imaging method Methods 0.000 claims description 7
 238000005314 correlation function Methods 0.000 claims description 4
 230000000875 corresponding Effects 0.000 claims description 4
 230000005012 migration Effects 0.000 claims description 4
 230000015572 biosynthetic process Effects 0.000 claims description 3
 238000004422 calculation algorithm Methods 0.000 claims description 3
 238000006243 chemical reaction Methods 0.000 claims description 3
 238000004587 chromatography analysis Methods 0.000 claims description 3
 238000005755 formation reaction Methods 0.000 claims description 3
 239000011435 rock Substances 0.000 claims description 3
 238000004088 simulation Methods 0.000 claims description 3
 230000000694 effects Effects 0.000 description 6
 238000001228 spectrum Methods 0.000 description 3
 230000003044 adaptive Effects 0.000 description 2
 230000001364 causal effect Effects 0.000 description 2
 238000005516 engineering process Methods 0.000 description 2
 239000004576 sand Substances 0.000 description 2
 230000003595 spectral Effects 0.000 description 2
 238000003325 tomography Methods 0.000 description 2
 210000003229 CMP Anatomy 0.000 description 1
 239000000470 constituent Substances 0.000 description 1
 238000010219 correlation analysis Methods 0.000 description 1
 230000002708 enhancing Effects 0.000 description 1
 238000007689 inspection Methods 0.000 description 1
 238000007781 preprocessing Methods 0.000 description 1
 230000001629 suppression Effects 0.000 description 1
 239000002699 waste material Substances 0.000 description 1
Abstract
The invention discloses a geological marker bedbased time difference analyzing and correcting method. The method comprises the following steps: acquiring geological singleshot data information of a seismic region through a singlepoint digital seismic method and an acquisition observation system, and transmitting the geological singleshot data information to a computer processing system; performing processing and analyzing on data by a computer processing system according to the explanation of an artificial geological marker bed, wherein the focal point is to determine the geological marker bed and perform quantitative calculation on crosscorrelation time difference and time difference correction and combination; finally, obtaining geological marker bedbased technical parameters and structural images, wherein the precision of the technical parameters is higher; the overall resolution is improved; the structural images are clearer. According to the method, the characteristics of relatively low frequency bandwidth and signaltonoise ratio of singlepoint digital data can be taken into consideration; the main geological marker bed level of the data and the signaltonoise ratio of the data are taken into the consideration when the time difference between tracks is calculated. The calculated time difference can better reflect the overall time difference of data tracks, so that the superposition precision and the resolution capability after combination can be improved better.
Description
Technical field
The present invention relates to oilgas exploration seism processing field, specifically a kind of stepout time analysis bearing calibration based on Geologic Indicators layer.
Background technology
Singlepoint digital landscape model belongs to the application initial stage, because the restriction being subject to equipment and investment is not promoted on a large scale.Be not a lot of for indoor combination method, majority directly combines after carrying out noise treatment to data, and the road combination TEC time error correction method retrieved has:
1, spectrum equalization adds non time difference combined method is a kind of seismic prospecting single detector recording channel indoor combination method.Mainly spectral equalization process is carried out to multiple seismic record signals of multiple single wave detector collection; Static TEC time error correction and dynamic float correction are carried out to the multiple seismic record signals through spectral equalization process, eliminates the dynamic float in described multiple seismic record signal and the static time difference; Combination superposition eliminates dynamic float, multiple seismic event trace record signals of the static time difference, effectively protects the frequency range of reflection wave, improves the resolution characteristic of reflection wave.
2, based on the consistance combination stacking method of adaptive matched filter.By building standard track, adopting adaptive matched filter method to carry out Concordance to the time difference between wavelet waveforms and road simultaneously, improving the effect of combination superposition.After correcting, the combination stack response of seismologic record can keep the original amplitude of reflection wave and waveform character preferably.
3, the static correcting method of the indoor combination of the wave detector of a kind of singlepoint high density earthquake.The feature of refraction time and waste mantle seismic event when utilizing the near surface of refraction wave to propagate, is expressed as the linear function of refraction wave slowness and the relative quiet school time difference first arrival time difference of adjacent geophone station, and adopts numerical statistic method to calculate the relatively quiet school time difference of each wave detector in geophone station combination range.
4, a kind of method of marine singlepoint acquired seismic data indoor dynamic combination.Utilize highdensity acquisition data, decay with selfadaptation noise pressing of surging after preservice, carry out random noise compacting with noncausal spatial prediction filter method, remove linear disturbance by FX territory prediction matching compacting linear disturbance method.Correct Dynamic combination method with phase spectrum and carry out phase shift TEC time error correction, after dynamic combined, export big gun collection record.The method is for the single inspection in sea highdensity acquisition geological data feature of noise, noise drawing method, noncausal spatial prediction filtering Attenuating Random Noise method, the FX territory prediction matching compacting linear disturbance method of selfadaptation being surged corrects dynamic combined with phase spectrum and effectively connects, various interference is effectively suppressed, reflected signal energy and continuity are obviously strengthened, and signal to noise ratio (S/N ratio) significantly improves.
These methods have respective advantage, and the TEC time error correction method based on Geologic Indicators layer of carrying out single track combination for indoor process has no corresponding Research Literature.
High density singlepoint digital geophone gathers data, has the advantages that single track receives, track pitch is little, resolution is high, signal to noise ratio (S/N ratio) is on the low side.Little track pitch receives and can carry out indoor combined treatment to strengthen received energy raising signal to noise ratio (S/N ratio).Owing to needing to there is the time difference between the seismic trace of combination, directly combination superposition weakens effective energy, have impact on the geology imaging effect after combination, reduces the effective frequency of data, be unfavorable for the highresolution imaging of data.During in order to combine, the effective constituent of seismic event is not lost, and ensures the resolution of data, needs to carry out the accurate Calculation of the time difference and correct before carrying out combination.
Summary of the invention:
The object of the present invention is to provide the stepout time analysis bearing calibration based on Geologic Indicators layer, carry out TEC time error correction by choosing Geologic Indicators layer thus protect effective highfrequency information in singlepoint Digital data processing.
In order to reach abovementioned purpose, present invention employs following technical scheme, based on the stepout time analysis bearing calibration of Geologic Indicators layer, comprising the following steps:
Gathering seismic region geology list big gun data message data by singlepoint digital seismics method and by gathering recording geometry, being then delivered to computer processing system;
Carry out data processing and analysis by computer processing system according to the explanation of artificial geology reference lamina, emphasis carries out the quantitative calculating of the crosscorrelation time difference and TEC time error correction and merging;
Finally obtain based on higher technical parameter, the better construct image clearly of whole resolution raising of the precision of Geologic Indicators layer.
The quantitative result of calculation of the described crosscorrelation time difference, depend on computation window position and time window in geological data signal to noise ratio (S/N ratio) height, the whether suitable degree of accuracy that can affect cross correlation results that height and window during crosscorrelation of geological data signal to noise ratio (S/N ratio) are chosen; So between crosscorrelation calculation road during the time difference, computation window is selected data signal to noise ratio (S/N ratio) higher and is had stable sign back layer position.
Explain at Geologic Indicators layer on the basis determined, the time window of design 250300ms, carries out the computing of crosscorrelation, obtains the time difference existed between two merging traffics at this time in window;
If x
_{1}(t) and x
_{2}t () is respectively two adjacent seismic trace signals, expect that earthquake is recorded as s (t), n
_{1}(t) and n
_{2}t average that () is respectively separate to be 0 variance be 1 white Gaussian noise, t
_{d}for time delay, α is attenuation coefficient, then seismic signal is: x
_{1}(t)=s (t)+n
_{1}(t); x
_{2}(t)=α s (tt
_{d})+n
_{2}(t);
X
_{1}(t) and x
_{2}t the cross correlation function of () is
${R}_{{x}_{1}{x}_{2}}\left(\mathrm{\τ}\right)=E\left[{x}_{1}\left(t\right){x}_{2}(t+\mathrm{\τ})\right]=\mathrm{\α}{R}_{\mathrm{ss}}(\mathrm{\τ}{t}_{d});$
Due to  R
_{ss}(τ) ≤R
_{ss}(0);
Therefore, R
_{ss}(τ) at τt
_{d}be maximum value when=0, successive relation tries to achieve t
_{d}, t now
_{d}be exactly seismic trace x
_{2}the correction time shift amount of (t); By seismic signal x
_{2}(t) time shift t
_{d}just achieve time shift calibration;
Utilize the time difference that crosscorrelation calculation goes out, TEC time error correction process is carried out to data and combination superposition in road after carrying out TEC time error correction.
The step that described computer processing system carries out the process of Geologic Indicators layer data and analysis comprises following:
Near surface process is carried out to data and eliminates static correction value;
Carry out the sound attenuation process of ground roll and random noise, improve data signal to noise ratio (S/N ratio);
Seismic section is carried out to explanation and the demarcation of Geologic Indicators layer in data work area;
Time window in carry out time difference between two neighboring traces and calculate;
In whole data, respectively a TEC time error correction is carried out to neighboring track;
Velocity analysis, normal moveout correction, superposition and migration imaging have been carried out to data after TEC time error correction.
Described to data carry out near surface process eliminate static correction value: by field survey gather earth's surface altitude figures, earth's surface elevation correction process is carried out to single big gun data; Adopt the formation velocity in the thickness of the low velocity layer of the near surface of field survey, low velocity layer, set up nearsurface model, utilize model constrained chromatography conversion method, calculate the correcting value of near surface, again consistency treatment is carried out to single big gun data, eliminate near surface and calculated to the impact produced the time difference.Model tuning computing formula:
Utilize earth's surface elevation E, datum elevation D, low reduction of speed layer thickness Z, low reduction of speed interval velocity Vw and basement rock speed Vb calculate correcting value.
The described sound attenuation process carrying out ground roll and random noise: need the process carrying out 2 aspects: a, the impact that ground roll is eliminated in region algorithm process is carried out on low frequency ground roll; Set up ground roll difference equation according to the physical model of ground roll, by the crossing prediction error simulation filtering in fx territory and tx territory, realize the object eliminating low frequency ground roll noise; B, adopt threedimensional useful signal to strengthen process to eliminate the impact of random noise; The method is the predicated error Wei Na filtering of an application frequency plane; This filtering supposes that predictable energy is useful signal, and dump energy is random noise; With correlation technique, input being calculated, determining to produce Wei Na filtering for predicting output, then by filtering for the treatment of data, when after the filtering of each window completing frequency plane, window carries out overlap, gets back to time domain export by antifouriertransform; Realize the random noise effective attenuation to geological data.
Described Geologic Indicators layer is determined: according to survey region seismotectonics stratigraphic interpretation data, analyze the Geologic Indicators stratum compared with shallowlayer, trend is moved towards, so that the design of computation window in time and the layer position of corresponding Seismic Stacked Section determination geology reference lamina.
Window during described designing and calculating: with the lineups time of seismic section Geologic Indicators layer, the center of window when being defined as, design one Geologic Indicators layer is included 250300ms time window, be decided to be and calculate the time window of the time difference.
Compared to prior art, the present invention has following beneficial effect:
This method can consider the bandwidth of singlepoint numerical data and the relatively low feature of signal to noise ratio (S/N ratio), considers the main geologic key horizon of data and the signal to noise ratio (S/N ratio) of data when calculating the time difference between road.Because it is very large that the signal to noise ratio (S/N ratio) of data calculates impact to the time difference, easily cause the error of calculation.The geology reflection coefficient being exactly Geologic Indicators layer position place is in addition strong, calculates the overall time difference that the time difference can reflect data track better, thus is more conducive to the superposition precision after improving combination and resolution characteristic.
Accompanying drawing explanation
Fig. 1 is the FB(flow block) of the treatment step of the stepout time analysis bearing calibration based on Geologic Indicators layer of the present invention;
Fig. 2 is for directly to combine superimposed image;
Fig. 3 combines superimposed image after conventional TEC time error correction;
Fig. 4 combines superimposed image after the inventive method TEC time error correction;
Fig. 5 be adopt static correction road, earth's surface combination and the inventive method correct combine after superpose contrast images;
Fig. 6 is that a specific embodiment carries out explanation and the uncalibrated image of Geologic Indicators layer to seismic section;
Fig. 7 is superimposed image after a specific embodiment earth's surface static correction combination;
Fig. 8 is that a specific embodiment adopts the inventive method to correct the rear superimposed image of combination;
Fig. 9 is migrated image after a specific embodiment earth's surface static correction combination;
Figure 10 is that a specific embodiment adopts the inventive method to correct combination to retrodeviate and move image.
Embodiment
Detailed description for the present invention and technology contents, coordinate accompanying drawing to be described as follows, but accompanying drawing only provides with reference to the use with explanation, is not used for being limited the present invention.
Fig. 1 is the treatment scheme adopted, and contains input singlepoint digital seismic data, the main processing steps of preprocessing of the information and this method.
Realization flow (Fig. 1):
1. geological data is inputted: input single big gun data that digital Singlereceiver seismics gathers.
2. near surface process: in order to improve the computational accuracy of the time difference, first will carry out the consistency treatment of near surface to data, the time difference elevation of near surface and low velocity layer caused eliminates.
By the earth's surface altitude figures that field survey gathers, earth's surface elevation correction process is carried out to single big gun data; Adopt the formation velocity in the thickness of the low velocity layer of the near surface of field survey, low velocity layer, set up nearsurface model, utilize model constrained chromatography conversion method, calculate the correcting value of near surface, again consistency treatment is carried out to single big gun data, eliminate near surface and calculated to the impact produced the time difference.Model tuning computing formula:
Utilize earth's surface elevation E, datum elevation D, low reduction of speed layer thickness Z, low reduction of speed interval velocity Vw and basement rock speed Vb calculate correcting value.
3. noise treatment: in order to improve the computational accuracy of the time difference, carries out Processing for removing to the noise jamming that data exist, and improves the effective energy of data.
Need the process carrying out 2 aspects: a, low frequency ground roll carried out to the impact that ground roll is eliminated in region algorithm process.Set up ground roll difference equation according to the physical model of ground roll, by the crossing prediction error simulation filtering in fx territory and tx territory, realize the object eliminating low frequency ground roll noise.B, adopt threedimensional useful signal to strengthen process to eliminate the impact of random noise.The method is the predicated error Wei Na filtering of an application frequency plane.This filtering supposes that predictable energy is useful signal, and dump energy is random noise.With correlation technique, input being calculated, determining to produce Wei Na filtering for predicting output, then by filtering for the treatment of data, when after the filtering of each window completing frequency plane, window carries out overlap, gets back to time domain export by antifouriertransform.Realize the random noise effective attenuation to geological data.
4. Geologic Indicators layer is determined: according to survey region seismotectonics stratigraphic interpretation data, analyze the Geologic Indicators stratum compared with shallowlayer, trend is moved towards, so that the design of computation window in time and the layer position of corresponding Seismic Stacked Section determination geology reference lamina.
5. window during designing and calculating: with the lineups time of seismic section Geologic Indicators layer, the center of window when being defined as, design one Geologic Indicators layer is included 250300ms time window, be decided to be and calculate the time window of the time difference.
6. crosscorrelation analysis: time window in adopt crosscorrelation calculation formula, calculate the time difference existed between adjacent twice.In the computation window of design, utilize cross correlation function
${R}_{{x}_{1}{x}_{2}}\left(\mathrm{\τ}\right)=E\left[{x}_{1}\left(t\right){x}_{2}(t+\mathrm{\τ})\right]=\mathrm{\α}{R}_{\mathrm{ss}}(\mathrm{\τ}{t}_{d})$ Calculate the crosscorrelation energy value between twice.Pass through  R
_{ss}(τ) ≤R
_{ss}(0) discriminant differentiates, at τt
_{d}r when=0
_{ss}(τ) be maximum value, the t tried to achieve
_{d}it is exactly the correction time shift amount between two seismic traces.
7. TEC time error correction between road: the correction time shift amount utilizing step 6. to calculate, the seismic trace of offset distance relatively far away carries out time shift calibration.Exactly the seismic trace of offset distance relatively far away is carried out t
_{d}time shift, obtain the new data as road combination after time shift.
Two adjacent seismic traces are carried out combination and superpose by 8. road combination superposition: the geological data after utilizing step 7. to carry out time shift.By the tracestacking combination after TEC time error correction, owing to eliminating the time difference existed between twice, thus improve the energy of inphase stacking.
9. export: the data after combination superposition are exported, carry out the process such as the velocity analysis in later stage, residual static correction, normal moveout correction, overall data superposition and migration imaging.Thus improve seismic imaging precision and strengthen seismic resolution ability.
Focus on:
The quantitative result of calculation of the described crosscorrelation time difference, depend on computation window position and time window in geological data signal to noise ratio (S/N ratio) height, the whether suitable degree of accuracy that can affect cross correlation results that height and window during crosscorrelation of geological data signal to noise ratio (S/N ratio) are chosen; So between crosscorrelation calculation road during the time difference, computation window is selected data signal to noise ratio (S/N ratio) higher and is had stable sign back layer position.Geologic Indicators reflection horizon has phase wiggles and stablizes, and energy stronger data signal to noise ratio (S/N ratio) is relatively high, and waveform can be followed the trail of continuously, and has good Hyperbolic Feature, wants high for time difference relative accuracy between the road calculated.See Fig. 2.
Explain at Geologic Indicators layer on the basis determined, the time window of design 250300ms, carries out the computing of crosscorrelation, obtains the time difference existed between two merging traffics at this time in window;
If x
_{1}(t) and x
_{2}t () is respectively two adjacent seismic trace signals, expect that earthquake is recorded as s (t), n
_{1}(t) and n
_{2}t average that () is respectively separate to be 0 variance be 1 white Gaussian noise, t
_{d}for time delay, α is attenuation coefficient, then seismic signal is: x
_{1}(t)=s (t)+n
_{1}(t); x
_{2}(t)=α s (tt
_{d})+n
_{2}(t);
X
_{1}(t) and x
_{2}t the cross correlation function of () is
${R}_{{x}_{1}{x}_{2}}\left(\mathrm{\τ}\right)=E\left[{x}_{1}\left(t\right){x}_{2}(t+\mathrm{\τ})\right]=\mathrm{\α}{R}_{\mathrm{ss}}(\mathrm{\τ}{t}_{d});$
Due to  R
_{ss}(τ) ≤R
_{ss}(0);
Therefore, R
_{ss}(τ) at τt
_{d}be maximum value when=0, successive relation tries to achieve t
_{d}, t now
_{d}be exactly seismic trace x
_{2}the correction time shift amount of (t); By seismic signal x
_{2}(t) time shift t
_{d}just achieve time shift calibration;
Utilize the time difference that crosscorrelation calculation goes out, TEC time error correction process is carried out to data and combination superposition in road after carrying out TEC time error correction.
Experimental study process: get a certain digital singlepoint acquiring data, track pitch is 10 meters.Superposition is combined after combination superposition and this method TEC time error correction after adopting adjacent 2 roads directly to combine superposition, earth's surface static correction respectively.
Fig. 2 directly combines superposition, and Fig. 3 is combination superposition after conventional TEC time error correction, and Fig. 4 is combination superposition after this method TEC time error correction.Be not difficult to find out, after Fig. 2 directly combines, lineups poor continuity, signal to noise ratio (S/N ratio) is low, and wave group feature is not obvious; Fig. 3 combined effect obtains certain improvement; Fig. 4 is owing to employing stepout time analysis bearing calibration between the road based on Geologic Indicators layer, and its combined effect is better than front 2, signal to noise ratio (S/N ratio), wave group feature and all increase with in phase continuity.Fig. 5 is superposing and contrast after adopting the combination of static correction road, earth's surface and this method to correct combine, is not difficult to find out: this method corrects the raising of the enhancing of the superposition signal to noise ratio (S/N ratio) after combining resolution, and tomography knows that wave group feature obviously, is beneficial to the Seismic Interpretation research in later stage.
Embodiment 1:
Area: certain oily district shlj area, singlepoint digital seismics method gathers data.
Acquisition parameter: gather recording geometry: 28L10S400T, blow out mode: 24956062495, single line receives number of channels: 400 roads, track pitch: 12.5m, CMP bin grid: 6.25m × 6.25m, degree of covering: 7(is horizontal) × 20(indulges)=140 times, receive linespacing: 125m, horizontal shotpoint spacing: 25m, longitudinal perpendicular offset: 125m, maximum offset: 3079m, minimum vertical misalignment distance: 6.25m.
Workload: gather data 29400 big gun, receive 110, line, bunch 42 is restrainted, full time data area 60km
^{2}.
Geological tasks: 1, improve data middle level stratigraphic resolution; 2, fracture system and breakpoint clear, clearcut; 3, wave group feature obviously, construct complete; 4, buried hill, deep structure and inside story clear.
Processing procedure: (carrying out according to treatment scheme) 1, data are carried out near surface process and eliminate static correction value; 2, carry out the sound attenuation process of ground roll and random noise, improve data signal to noise ratio (S/N ratio); 3, seismic section is carried out to explanation and the demarcation (Fig. 6) of Geologic Indicators layer in data work area; 4, time window in carry out time difference between two neighboring traces and calculate; 5, in whole data, respectively a TEC time error correction is carried out to neighboring track; 6, velocity analysis, normal moveout correction, superposition (Fig. 7, Fig. 8) and migration imaging (Fig. 9, Figure 10) have been carried out to data after TEC time error correction.
Effect analysis: Fig. 7 is that employing earth's surface corrects the rear superposition of combination, and Fig. 8 superposes after this method correction combination, and Fig. 9 adopts earth's surface correction combination to retrodeviate to move, and Figure 10 is the migrated section after drawing method corrects combination.
Analyze known: adopt the data after this technical finesse to have: 1, data signal to noise ratio (S/N ratio) is moderate, and wave group feature is obvious; 2, entirety resolution improves; 3, the weak layer signal waveform character between strongly reflecting layer is clear, differentiates and identifies that sand body ability strengthens; 4, tomography and breakpoint are simply clear; 5, buried hill, deep structure is clear; 6, data is beneficial to the precision that sand body is explained and improved Seismic Reservoir Prediction attributive analysis.
The foregoing is only preferred embodiment of the present invention, be not used to limit the scope of the claims of the present invention, other use the equivalence of patent spirit of the present invention to change, and all should all belong to the scope of the claims of the present invention.
Claims (8)
1., based on the stepout time analysis bearing calibration of Geologic Indicators layer, comprise the following steps:
Gathering seismic region geology list big gun data message data by singlepoint digital seismics method and by gathering recording geometry, being then delivered to computer processing system;
Carry out data processing and analysis by computer processing system according to the explanation of artificial geology reference lamina, emphasis carries out the quantitative calculating of the crosscorrelation time difference and TEC time error correction and merging;
Finally obtain based on higher technical parameter, the better construct image clearly of whole resolution raising of the precision of Geologic Indicators layer.
2. the stepout time analysis bearing calibration based on Geologic Indicators layer according to claim 1, it is characterized in that, the quantitative result of calculation of the described crosscorrelation time difference, depend on computation window position and time window in geological data signal to noise ratio (S/N ratio) height, the whether suitable degree of accuracy that can affect cross correlation results that height and window during crosscorrelation of geological data signal to noise ratio (S/N ratio) are chosen; So between crosscorrelation calculation road during the time difference, computation window is selected data signal to noise ratio (S/N ratio) higher and is had stable sign back layer position.
3. the stepout time analysis bearing calibration based on Geologic Indicators layer according to claim 1, it is characterized in that, explain at Geologic Indicators layer on the basis determined, the time window of design 250300ms, at this time, in window, carry out the computing of crosscorrelation, obtain the time difference existed between two merging traffics;
If x
_{1}(t) and x
_{2}t () is respectively two adjacent seismic trace signals, expect that earthquake is recorded as s (t), n
_{1}(t) and n
_{2}t average that () is respectively separate to be 0 variance be 1 white Gaussian noise, t
_{d}for time delay, α is attenuation coefficient, then seismic signal is: x
_{1}(t)=s (t)+n
_{1}(t); x
_{2}(t)=α s (tt
_{d})+n
_{2}(t);
X
_{1}(t) and x
_{2}t the cross correlation function of () is
Due to  R
_{ss}(τ) ≤R
_{ss}(0);
Therefore, R
_{ss}(τ) at τt
_{d}be maximum value when=0, successive relation tries to achieve t
_{d}, t now
_{d}be exactly seismic trace x
_{2}the correction time shift amount of (t); By seismic signal x
_{2}(t) time shift t
_{d}just achieve time shift calibration;
Utilize the time difference that crosscorrelation calculation goes out, TEC time error correction process is carried out to data and combination superposition in road after carrying out TEC time error correction.
4. the stepout time analysis bearing calibration based on Geologic Indicators layer according to claim 1, is characterized in that, the step that described computer processing system carries out the process of Geologic Indicators layer data and analysis comprises following:
Near surface process is carried out to data and eliminates static correction value;
Carry out the sound attenuation process of ground roll and random noise, improve data signal to noise ratio (S/N ratio);
Seismic section is carried out to explanation and the demarcation of Geologic Indicators layer in data work area;
Time window in carry out time difference between two neighboring traces and calculate;
In whole data, respectively a TEC time error correction is carried out to neighboring track;
Velocity analysis, normal moveout correction, superposition and migration imaging have been carried out to data after TEC time error correction.
5. the stepout time analysis bearing calibration based on Geologic Indicators layer according to claim 4, it is characterized in that, described to data carry out near surface process eliminate static correction value: by field survey gather earth's surface altitude figures, earth's surface elevation correction process is carried out to single big gun data; Adopt the formation velocity in the thickness of the low velocity layer of the near surface of field survey, low velocity layer, set up nearsurface model, utilize model constrained chromatography conversion method, calculate the correcting value of near surface, again consistency treatment is carried out to single big gun data, eliminate near surface and calculated to the impact produced the time difference.Model tuning computing formula:
Utilize earth's surface elevation E, datum elevation D, low reduction of speed layer thickness Z, low reduction of speed interval velocity Vw and basement rock speed Vb calculate correcting value.
6. the stepout time analysis bearing calibration based on Geologic Indicators layer according to claim 4, it is characterized in that, described in carry out the sound attenuation process of ground roll and random noise: the process needing to carry out 2 aspects: a, the impact that ground roll is eliminated in region algorithm process is carried out on low frequency ground roll; Set up ground roll difference equation according to the physical model of ground roll, by the crossing prediction error simulation filtering in fx territory and tx territory, realize the object eliminating low frequency ground roll noise; B, adopt threedimensional useful signal to strengthen process to eliminate the impact of random noise; The method is the predicated error Wei Na filtering of an application frequency plane; This filtering supposes that predictable energy is useful signal, and dump energy is random noise; With correlation technique, input being calculated, determining to produce Wei Na filtering for predicting output, then by filtering for the treatment of data, when after the filtering of each window completing frequency plane, window carries out overlap, gets back to time domain export by antifouriertransform; Realize the random noise effective attenuation to geological data.
7. the stepout time analysis bearing calibration based on Geologic Indicators layer according to claim 4, it is characterized in that, described Geologic Indicators layer is determined: according to survey region seismotectonics stratigraphic interpretation data, analyze the Geologic Indicators stratum compared with shallowlayer, trend is moved towards, so that the design of computation window in time and the layer position of corresponding Seismic Stacked Section determination geology reference lamina.
8. the stepout time analysis bearing calibration based on Geologic Indicators layer according to claim 4, it is characterized in that, window during described designing and calculating: with the lineups time of seismic section Geologic Indicators layer, the center of window when being defined as, design one Geologic Indicators layer is included 250300ms time window, be decided to be and calculate the time window of the time difference.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201310520070.5A CN104570116A (en)  20131029  20131029  Geological marker bedbased time difference analyzing and correcting method 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201310520070.5A CN104570116A (en)  20131029  20131029  Geological marker bedbased time difference analyzing and correcting method 
Publications (1)
Publication Number  Publication Date 

CN104570116A true CN104570116A (en)  20150429 
Family
ID=53086630
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201310520070.5A Pending CN104570116A (en)  20131029  20131029  Geological marker bedbased time difference analyzing and correcting method 
Country Status (1)
Country  Link 

CN (1)  CN104570116A (en) 
Cited By (3)
Publication number  Priority date  Publication date  Assignee  Title 

CN106353818A (en) *  20160812  20170125  中国石油天然气集团公司  Seismic data time difference correcting method and system 
CN107515422A (en) *  20170830  20171226  中国海洋大学  A kind of highprecision shallow seismic profile acquisition system 
CN110998369A (en) *  20170601  20200410  沙特阿拉伯石油公司  Detecting underground structures 
Citations (3)
Publication number  Priority date  Publication date  Assignee  Title 

US20060253256A1 (en) *  20020927  20061109  Robertsson Johan O A  Processing seismic data 
CN102221709A (en) *  20110601  20111019  中国石油集团川庆钻探工程有限公司地球物理勘探公司  Velocity analysis and dynamic correction method based on formation parameter information 
CN102759746A (en) *  20110428  20121031  中国石油天然气集团公司  Method for inverting anisotropy parameters using variable offset vertical seismic profile data 

2013
 20131029 CN CN201310520070.5A patent/CN104570116A/en active Pending
Patent Citations (3)
Publication number  Priority date  Publication date  Assignee  Title 

US20060253256A1 (en) *  20020927  20061109  Robertsson Johan O A  Processing seismic data 
CN102759746A (en) *  20110428  20121031  中国石油天然气集团公司  Method for inverting anisotropy parameters using variable offset vertical seismic profile data 
CN102221709A (en) *  20110601  20111019  中国石油集团川庆钻探工程有限公司地球物理勘探公司  Velocity analysis and dynamic correction method based on formation parameter information 
NonPatent Citations (5)
Title 

中国石油天然气总公司技术推广中心: "《油气勘探综合配套技术应用范例》", 31 December 1996, 石油工业出版社 * 
李辉峰 等: "《地震勘探新技术》", 31 December 2009, 石油工业出版社 * 
熊章强 等: "《地震勘探》", 30 September 2010, 中南大学出版社 * 
赵金洲: "《石油物探技术》", 31 August 2013, 中国石化出版社 * 
金龙 等: ""时移地震互均衡法混合因素处理技术"", 《石油地球物理勘探》 * 
Cited By (5)
Publication number  Priority date  Publication date  Assignee  Title 

CN106353818A (en) *  20160812  20170125  中国石油天然气集团公司  Seismic data time difference correcting method and system 
CN110998369A (en) *  20170601  20200410  沙特阿拉伯石油公司  Detecting underground structures 
CN110998369B (en) *  20170601  20210924  沙特阿拉伯石油公司  Detecting underground structures 
CN107515422A (en) *  20170830  20171226  中国海洋大学  A kind of highprecision shallow seismic profile acquisition system 
CN107515422B (en) *  20170830  20190419  中国海洋大学  A kind of highprecision shallow seismic profile acquisition system 
Similar Documents
Publication  Publication Date  Title 

US9075163B2 (en)  Interferometric seismic data processing  
US8737165B2 (en)  Interferometric seismic data processing for a towed marine survey  
CN102841379B (en)  Method for analyzing prestack time migration and speed based on common scatter point channel set  
CN104570125A (en)  Method for improving imaging velocity model precision by utilizing well data  
CN105388518A (en)  Centroid frequency and spectral ratio integrated borehole seismic quality factor inversion method  
CN106896409B (en)  A kind of varying depth cable ghost reflection drawing method based on wave equation boundary values inverting  
CN102176054B (en)  Nearsurface comprehensive information processing explanation method  
CN104216014A (en)  Seismic signal frequency division processing method  
CN103926623B (en)  Method for suppressing reverse time migration low frequency noise  
CN104820242B (en)  A kind of road collection amplitude towards prestack inversion divides compensation method  
CN104459794A (en)  Commonreflectionpoint gather timevariable time difference correction method and device  
CN101201409B (en)  Method for revising earthquake data phase  
CN103675900A (en)  Method for determining optimum velocity profile of convertedwave prestack time migration during seismic data processing process  
CN107526101B (en)  A kind of acquisition and processing method obtaining earthquake reflected wave  
CN102565852B (en)  Angle domain prestack offset data processing method aiming to detect oilgasbearing property of reservoir  
CN107807393A (en)  Separate unit station collection preliminary wave Enhancement Method based on seismic interference method  
CN104570116A (en)  Geological marker bedbased time difference analyzing and correcting method  
CN102053260B (en)  Method for acquiring azimuth velocity of primary wave and method for processing earthquake data  
CN105607119B (en)  Nearsurface model construction method and static correction value acquiring method  
CN109669212B (en)  Seismic data processing method, stratum quality factor estimation method and device  
CN105425300A (en)  Residual static correction method  
CN105093296A (en)  Method and device for optimizing observation system  
EA030770B1 (en)  System and method for seismic adaptive optics  
CN105093301B (en)  The generation method and device of common imaging point angle of reflection angle gathers  
CN104977615B (en)  A kind of multiple ripple drawing method of deep water OBC data based on modeling statistics pickup 
Legal Events
Date  Code  Title  Description 

C06  Publication  
PB01  Publication  
C10  Entry into substantive examination  
SE01  Entry into force of request for substantive examination  
WD01  Invention patent application deemed withdrawn after publication  
WD01  Invention patent application deemed withdrawn after publication 
Application publication date: 20150429 