CN104570116A - Geological marker bed-based time difference analyzing and correcting method - Google Patents

Geological marker bed-based time difference analyzing and correcting method Download PDF

Info

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
Application number
CN201310520070.5A
Other languages
Chinese (zh)
Inventor
郭树祥
梁鸿贤
李建明
孔省吾
刘媛
毕丽飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Geophysical Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Geophysical Research Institute of Sinopec Shengli Oilfield Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Geophysical Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN201310520070.5A priority Critical patent/CN104570116A/en
Publication of CN104570116A publication Critical patent/CN104570116A/en
Pending legal-status Critical Current

Links

Abstract

The invention discloses a geological marker bed-based time difference analyzing and correcting method. The method comprises the following steps: acquiring geological single-shot data information of a seismic region through a single-point digital seismic method and an acquisition observation system, and transmitting the geological single-shot 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 cross-correlation time difference and time difference correction and combination; finally, obtaining geological marker bed-based 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 signal-to-noise ratio of single-point digital data can be taken into consideration; the main geological marker bed level of the data and the signal-to-noise 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

Based on the step-out time analysis bearing calibration of Geologic Indicators layer
Technical field
The present invention relates to oil-gas exploration seism processing field, specifically a kind of step-out time analysis bearing calibration based on Geologic Indicators layer.
Background technology
Single-point 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 single-point 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 single-point acquired seismic data indoor dynamic combination.Utilize high-density acquisition data, decay with self-adaptation noise pressing of surging after pre-service, carry out random noise compacting with non-causal 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 high-density acquisition geological data feature of noise, noise drawing method, non-causal spatial prediction filtering Attenuating Random Noise method, the FX territory prediction matching compacting linear disturbance method of self-adaptation 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 single-point 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 high-resolution 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 step-out time analysis bearing calibration based on Geologic Indicators layer, carry out TEC time error correction by choosing Geologic Indicators layer thus protect effective high-frequency information in single-point Digital data processing.
In order to reach above-mentioned purpose, present invention employs following technical scheme, based on the step-out time analysis bearing calibration of Geologic Indicators layer, comprising the following steps:
Gathering seismic region geology list big gun data message data by single-point 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 cross-correlation 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 cross-correlation 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 cross-correlation of geological data signal to noise ratio (S/N ratio) are chosen; So between cross-correlation 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 250-300ms, carries out the computing of cross-correlation, obtains the time difference existed between two merging traffics at this time in window;
If x 1(t) and x 2t () is respectively two adjacent seismic trace signals, expect that earthquake is recorded as s (t), n 1(t) and n 2t average that () is respectively separate to be 0 variance be 1 white Gaussian noise, t dfor time delay, α is attenuation coefficient, then seismic signal is: x 1(t)=s (t)+n 1(t); x 2(t)=α s (t-t d)+n 2(t);
X 1(t) and x 2t the cross correlation function of () is R x 1 x 2 ( τ ) = E [ x 1 ( t ) x 2 ( t + τ ) ] = α R ss ( τ - t d ) ;
Due to | R ss(τ) |≤R ss(0);
Therefore, R ss(τ) at τ-t dbe maximum value when=0, successive relation tries to achieve t d, t now dbe exactly seismic trace x 2the correction time shift amount of (t); By seismic signal x 2(t) time shift t djust achieve time shift calibration;
Utilize the time difference that cross-correlation 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 near-surface 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:
Δ τ D = - Z v w + D - E + Z v b
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 f-x territory and t-x territory, realize the object eliminating low frequency ground roll noise; B, adopt three-dimensional 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 anti-fourier-transform; 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 shallow-layer, 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 250-300ms 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 single-point 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 step-out 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 single-point digital seismic data, the main processing steps of pre-processing of the information and this method.
Realization flow (Fig. 1):
1. geological data is inputted: input single big gun data that digital Single-receiver 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 near-surface 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:
Δ τ D = - Z v w + D - E + Z v b
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 f-x territory and t-x territory, realize the object eliminating low frequency ground roll noise.B, adopt three-dimensional 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 anti-fourier-transform.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 shallow-layer, 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 250-300ms time window, be decided to be and calculate the time window of the time difference.
6. cross-correlation analysis: time window in adopt cross-correlation 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 ( τ ) = E [ x 1 ( t ) x 2 ( t + τ ) ] = α R ss ( τ - t d ) Calculate the cross-correlation energy value between twice.Pass through | R ss(τ) |≤R ss(0) discriminant differentiates, at τ-t dr when=0 ss(τ) be maximum value, the t tried to achieve dit 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 dtime 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 trace-stacking combination after TEC time error correction, owing to eliminating the time difference existed between twice, thus improve the energy of in-phase 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 cross-correlation 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 cross-correlation of geological data signal to noise ratio (S/N ratio) are chosen; So between cross-correlation 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 250-300ms, carries out the computing of cross-correlation, obtains the time difference existed between two merging traffics at this time in window;
If x 1(t) and x 2t () is respectively two adjacent seismic trace signals, expect that earthquake is recorded as s (t), n 1(t) and n 2t average that () is respectively separate to be 0 variance be 1 white Gaussian noise, t dfor time delay, α is attenuation coefficient, then seismic signal is: x 1(t)=s (t)+n 1(t); x 2(t)=α s (t-t d)+n 2(t);
X 1(t) and x 2t the cross correlation function of () is R x 1 x 2 ( τ ) = E [ x 1 ( t ) x 2 ( t + τ ) ] = α R ss ( τ - t d ) ;
Due to | R ss(τ) |≤R ss(0);
Therefore, R ss(τ) at τ-t dbe maximum value when=0, successive relation tries to achieve t d, t now dbe exactly seismic trace x 2the correction time shift amount of (t); By seismic signal x 2(t) time shift t djust achieve time shift calibration;
Utilize the time difference that cross-correlation 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 single-point 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 step-out 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, single-point digital seismics method gathers data.
Acquisition parameter: gather recording geometry: 28L10S400T, blow out mode: 2495-6-0-6-2495, 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 line-spacing: 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, clear-cut; 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 step-out time analysis bearing calibration of Geologic Indicators layer, comprise the following steps:
Gathering seismic region geology list big gun data message data by single-point 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 cross-correlation 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 step-out 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 cross-correlation 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 cross-correlation of geological data signal to noise ratio (S/N ratio) are chosen; So between cross-correlation 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 step-out 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 250-300ms, at this time, in window, carry out the computing of cross-correlation, obtain the time difference existed between two merging traffics;
If x 1(t) and x 2t () is respectively two adjacent seismic trace signals, expect that earthquake is recorded as s (t), n 1(t) and n 2t average that () is respectively separate to be 0 variance be 1 white Gaussian noise, t dfor time delay, α is attenuation coefficient, then seismic signal is: x 1(t)=s (t)+n 1(t); x 2(t)=α s (t-t d)+n 2(t);
X 1(t) and x 2t the cross correlation function of () is
Due to | R ss(τ) |≤R ss(0);
Therefore, R ss(τ) at τ-t dbe maximum value when=0, successive relation tries to achieve t d, t now dbe exactly seismic trace x 2the correction time shift amount of (t); By seismic signal x 2(t) time shift t djust achieve time shift calibration;
Utilize the time difference that cross-correlation 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 step-out 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 step-out 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 near-surface 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 step-out 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 f-x territory and t-x territory, realize the object eliminating low frequency ground roll noise; B, adopt three-dimensional 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 anti-fourier-transform; Realize the random noise effective attenuation to geological data.
7. the step-out 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 shallow-layer, 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 step-out 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 250-300ms time window, be decided to be and calculate the time window of the time difference.
CN201310520070.5A 2013-10-29 2013-10-29 Geological marker bed-based time difference analyzing and correcting method Pending CN104570116A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310520070.5A CN104570116A (en) 2013-10-29 2013-10-29 Geological marker bed-based time difference analyzing and correcting method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310520070.5A CN104570116A (en) 2013-10-29 2013-10-29 Geological marker bed-based time difference analyzing and correcting method

Publications (1)

Publication Number Publication Date
CN104570116A true CN104570116A (en) 2015-04-29

Family

ID=53086630

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310520070.5A Pending CN104570116A (en) 2013-10-29 2013-10-29 Geological marker bed-based time difference analyzing and correcting method

Country Status (1)

Country Link
CN (1) CN104570116A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106353818A (en) * 2016-08-12 2017-01-25 中国石油天然气集团公司 Seismic data time difference correcting method and system
CN107515422A (en) * 2017-08-30 2017-12-26 中国海洋大学 A kind of high-precision shallow seismic profile acquisition system
CN110998369A (en) * 2017-06-01 2020-04-10 沙特阿拉伯石油公司 Detecting underground structures

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060253256A1 (en) * 2002-09-27 2006-11-09 Robertsson Johan O A Processing seismic data
CN102221709A (en) * 2011-06-01 2011-10-19 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Velocity analysis and dynamic correction method based on formation parameter information
CN102759746A (en) * 2011-04-28 2012-10-31 中国石油天然气集团公司 Method for inverting anisotropy parameters using variable offset vertical seismic profile data

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060253256A1 (en) * 2002-09-27 2006-11-09 Robertsson Johan O A Processing seismic data
CN102759746A (en) * 2011-04-28 2012-10-31 中国石油天然气集团公司 Method for inverting anisotropy parameters using variable offset vertical seismic profile data
CN102221709A (en) * 2011-06-01 2011-10-19 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Velocity analysis and dynamic correction method based on formation parameter information

Non-Patent Citations (5)

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

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106353818A (en) * 2016-08-12 2017-01-25 中国石油天然气集团公司 Seismic data time difference correcting method and system
CN110998369A (en) * 2017-06-01 2020-04-10 沙特阿拉伯石油公司 Detecting underground structures
CN110998369B (en) * 2017-06-01 2021-09-24 沙特阿拉伯石油公司 Detecting underground structures
CN107515422A (en) * 2017-08-30 2017-12-26 中国海洋大学 A kind of high-precision shallow seismic profile acquisition system
CN107515422B (en) * 2017-08-30 2019-04-19 中国海洋大学 A kind of high-precision 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 pre-stack 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) Near-surface 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) Common-reflection-point gather time-variable time difference correction method and device
CN101201409B (en) Method for revising earthquake data phase
CN103675900A (en) Method for determining optimum velocity profile of converted-wave pre-stack time migration during seismic data processing process
CN107526101B (en) A kind of acquisition and processing method obtaining earthquake reflected wave
CN102565852B (en) Angle domain pre-stack offset data processing method aiming to detect oil-gas-bearing property of reservoir
CN107807393A (en) Separate unit station collection preliminary wave Enhancement Method based on seismic interference method
CN104570116A (en) Geological marker bed-based time difference analyzing and correcting method
CN102053260B (en) Method for acquiring azimuth velocity of primary wave and method for processing earthquake data
CN105607119B (en) Near-surface 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