CN110005396B - Depth elastic correction method for logging wave train data - Google Patents
Depth elastic correction method for logging wave train data Download PDFInfo
- Publication number
- CN110005396B CN110005396B CN201811512074.8A CN201811512074A CN110005396B CN 110005396 B CN110005396 B CN 110005396B CN 201811512074 A CN201811512074 A CN 201811512074A CN 110005396 B CN110005396 B CN 110005396B
- Authority
- CN
- China
- Prior art keywords
- depth
- wave
- correction
- head
- logging
- 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 29
- 238000005070 sampling Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/04—Measuring depth or liquid level
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a depth elastic correction method for logging wave train data, which comprises the following steps: 1) Advanced traveling wave column logging; 2) Calculating a head wave arrival time curve; 3) Performing elastic depth correction on each head wave arrival time curve by taking the standard depth conventional logging curve as reference; 4) Respectively aligning the arrival times of the head waves of the wave trains received by each receiver by taking the arrival times of the head waves at any fixed depth as a reference; 5) Performing elastic depth correction on the arrival times of the wave train head waves received by the plurality of aligned receivers by taking a depth correction reference curve as a reference; 6) Depth point D 'of head wave to time curve after elastic depth correction' i0 For reference, repeating the steps 2) and 3) to obtain the point number delta P needing to be corrected 1j (ii) a 7) Number of points Δ P corrected as required 1j The head wave of the wave train data of each depth is horizontally moved by delta P after the depth is elastically corrected 1j According to the method, the correction of the depth position and the correction of the arrival time of the transverse waveform can be comprehensively considered, so that the elastic correction of the depth of the logging wave train data is realized.
Description
Technical Field
The invention belongs to the field of geophysical logging, and relates to a depth elastic correction method for logging wave train data.
Background
In the field logging operation process, due to the influence of various factors such as a logging instrument, a stratum, logging speed and the like at that time, the measured curve data usually cannot truly reflect the geological condition of the underground actual corresponding depth, so that the depth correction needs to be carried out on the logging curve data.
In the existing logging data depth correction technology, a mature correction method is provided for a conventional logging curve (namely a two-dimensional logging curve), and after depth correction, data which are missing in depth are obtained in a linear interpolation mode and the like under the condition that characteristic points of the logging curve are guaranteed. This type of method is applicable to data containing only longitudinal (depth) information. However, for some array data, which contains not only information in the longitudinal direction but also information in the transverse direction, the existing well log depth correction method is not applicable.
The wave train data is a special array logging data, and not only has corresponding information on the depth, but also the arrival time of the wave train is corresponding to the waveform in the transverse direction. The existing logging curve depth correction method only considers the correction of a depth position and does not consider the correction of the arrival time of a transverse waveform.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a depth elastic correction method for logging wave train data, which can comprehensively consider the correction of a depth position and the correction of the arrival time of a transverse waveform to realize the elastic correction of the depth of the logging wave train data.
In order to achieve the purpose, the depth elasticity correction method of the logging wave train data comprises the following steps:
1) Carrying out traveling wave train logging to obtain a group of logging data and a plurality of groups of wave train data;
2) Calculating a head wave arrival time curve according to wave train data before depth correction;
3) Performing elastic depth correction on each head wave arrival time curve by taking the standard depth conventional logging curve as reference;
4) Respectively aligning the arrival times of the head waves of the wave trains received by each receiver by taking the arrival times of the head waves at any fixed depth as a reference;
5) Performing elastic depth correction on the arrival times of the wave train head waves received by the plurality of aligned receivers by taking a depth correction reference curve as a reference;
6) Depth point D 'of head wave to time curve after elastic depth correction' i0 For reference, repeating the steps 2) and 3) to obtain the point number delta P needing to be corrected 1j ;
7) The point number delta P needing to be corrected is obtained according to the step 6) 1j Horizontally moving the head wave of the wave train data of each depth by delta P after elastically correcting the depth 1j And obtaining the final elastic correction depth wave train data by using the data points.
The wave train logging comprises digital acoustic logging, variable density acoustic logging and array acoustic logging.
And 2) calculating a head wave arrival time curve by adopting a threshold method, a speed ratio method or a slowness time correlation method according to wave train data before depth correction.
The specific operation of the step 4) is as follows:
41 Setting logging instrument to receive eight-channel wave train data and determine head wave arrival time reference depth D i0 And calculating therefrom a head wave arrival time reference depth D i0 Corresponding arrival time TT of each first wave i0j0 ,j=1,…,8;
42 Calculate the time-of-arrival correction Δ TT for each of the eight channel-train head waves at each depth ij Then, the head wave arrival time correction amount Δ TT corresponding to each depth 8-channel wave train ij Comprises the following steps:
ΔTT ij =TT ij -TT i0j0
where i denotes the ith depth point, i =1, \ 8230;, i max ;
43 ) head wave arrival time correction quantity Δ TT obtained according to 42) ij Converting the amount of correction for the amount of time to the number of points Δ P that need to be corrected ij Wherein, in the step (A),
ΔP ij =ΔTT ij /ts
wherein ts represents the time interval of the wave train sampling;
44 Root of Chinese YamNumber of points Δ P corrected as needed ij At depth point D i0 Taking the arrival time of the wave train head wave as a reference, and aligning the wave train data at other depths with the wave train data at the depth.
The invention has the following beneficial effects:
when the depth elastic correction method for the logging wave train data is specifically operated, elastic depth correction is carried out on each head wave arrival time curve by taking a standard depth conventional logging curve as a reference, the arrival time of the head wave at any fixed depth is taken as a reference to respectively carry out dead center elimination on the head wave arrival time of the wave train received by each receiver, elastic depth correction is carried out on the head wave arrival times of the wave train received by a plurality of aligned receivers by taking the depth correction reference curve as a reference to obtain the number of points to be corrected, then horizontal movement is carried out on the head wave of each depth wave train data according to the number of points to be corrected to obtain the wave train data of the elastic correction depth, and the elastic correction of the logging wave train data depth is realized by considering the correction of the depth position and the correction of the arrival time of the transverse waveform.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a schematic diagram of a plurality of wave trains in accordance with the present invention;
FIG. 3 is a schematic diagram of multiple depth wavetrains after alignment of the head wave.
Detailed Description
The invention is described in further detail below with reference to the accompanying drawings:
referring to fig. 1, the method for correcting the depth elasticity of the log wave train data comprises the following steps:
1) Carrying out travelling wave train logging to obtain a group of logging data and a plurality of groups of wave train data, wherein the wave train logging comprises digital acoustic logging, variable density acoustic logging and array acoustic logging, and referring to figure 2;
2) Calculating a head wave arrival time curve according to the wave train data before depth correction, wherein the head wave arrival time curve is calculated according to the wave train data before depth correction by adopting a threshold method, a speed ratio method or a slowness time correlation method;
3) Performing elastic depth correction on each head wave arrival time curve by taking the standard depth conventional logging curve as reference;
4) Respectively aligning the arrival times of the head waves of the wave trains received by each receiver by taking the arrival times of the head waves at any fixed depth as a reference;
the specific operation of the step 4) is as follows:
41 Setting logging instrument to receive eight-channel wave train data and determine head wave arrival time reference depth D i0 And calculating therefrom a head wave arrival time reference depth D i0 Corresponding first wave arrival time TT i0j0 ,j=1,…,8;
42 Calculates the first-wave arrival time correction amount Δ TT of eight-channel wave trains at each depth ij Then, the head wave arrival time correction amount Δ TT corresponding to each 8-channel wave train of depth ij Comprises the following steps:
ΔTT ij =TT ij -TT i0j0
where i denotes the ith depth point, i =1, \8230;, i max ;
43 ) head wave arrival time correction amount Δ TT obtained according to 42) ij Converting the amount of correction for the amount of time to the number of points Δ P to be corrected ij Wherein, in the process,
ΔP ij =ΔTT ij /ts
wherein ts represents the time interval of the wave train sampling;
44 Point number Δ P corrected as needed ij By depth point D i0 Taking the arrival time of the wave train head wave as a reference, and aligning the wave train data at other depths with the wave train data at the depth.
5) Elastic depth correction is carried out on the arrival times of the wave train head waves received by the plurality of aligned receivers by taking a depth correction reference curve as a reference, and a corrected depth point D is recorded i0 Corresponding depth point D' i0 ;
And (3) performing elastic depth correction on each aligned depth wave train data by taking the depth correction reference curve as a reference, wherein the elastic depth correction mode is synchronous to the step 3). In the elastic depth correction process, a processing mode similar to a one-dimensional curve and linear interpolation are adopted for missing depth point data, so that the method is only suitable for non-dispersive wave trains, and the speed is not changed for the same receiving array.
6) The head wave arrival time curve after elastic depth correction is defined as a depth point D' i0 For reference, repeating the steps 2) and 3) to obtain the point number delta P needing to be corrected 1j ;
7) The point number delta P needing to be corrected is obtained according to the step 6) 1j The head wave of the wave train data of each depth is horizontally moved by delta P after the depth is elastically corrected 1j And obtaining the final elastic correction depth wave train data by using the data points.
Claims (3)
1. A depth elasticity correction method for logging wave train data is characterized by comprising the following steps:
1) Firstly, performing travelling wave train logging to obtain a group of logging data and a plurality of groups of wave train data;
2) Calculating a head wave arrival time curve according to wave train data before depth correction;
3) Performing elastic depth correction on each head wave arrival time curve by taking the standard depth conventional logging curve as reference;
4) Respectively aligning the arrival times of the head waves of the wave trains received by each receiver by taking the arrival times of the head waves at any fixed depth as a reference;
5) Performing elastic depth correction on the arrival times of the wave train head waves received by the plurality of aligned receivers by taking a depth correction reference curve as a reference;
6) The head wave arrival time curve after elastic depth correction is divided into a depth point D i0 For reference, repeating the steps 2) and 3) to obtain the point number delta P needing to be corrected 1j ;
7) The point number delta P needing to be corrected is obtained according to the step 6) 1j Horizontally moving the head wave of the wave train data of each depth by delta P after elastically correcting the depth 1j Obtaining the final elastic correction depth wave train data by using the data points;
the specific operation of the step 4) is as follows:
41 Setting logging instrument to receive eight wave trains data and determining reference depth point D of arrival time of head wave i0 And calculating the reference depth point D of the head wave in time i0 Corresponding arrival time TT of each first wave i0j0 ,j=1,...,8;
42 Calculate the time-of-arrival correction Δ TT for each of the eight channel-train head waves at each depth ij Then, the head wave arrival time correction amount Δ TT corresponding to each 8-channel wave train of depth ij Comprises the following steps:
ΔTT ij =TT ij -TT i0j0
where i denotes the ith depth point, i =1 max (ii) a j denotes the serial number of the receiver, j =1, \ 8230;, 8
43 ) head wave arrival time correction quantity Δ TT obtained according to 42) ij Converting the amount of time correction to the number of points Δ P to be corrected ij Wherein, in the step (A),
ΔP ij =ΔTT ij /ts
wherein ts represents the time interval of the wave train sampling;
44 Point number Δ P corrected as needed ij By reference to depth point D i0 Taking the arrival time of the wave train head wave as a reference, and aligning the wave train data at other depths with the wave train data at the depth.
2. The method of depth-elastic correction of logging wavetrain data according to claim 1, wherein the wavetrain logging comprises digital sonic logging, variable density sonic logging or array sonic logging.
3. The method for depth-elastic correction of logging wavetrain data according to claim 1, wherein in step 2), the head-wave arrival-time curve is calculated by using a threshold method, a velocity ratio method or a slowness time correlation method according to the wavetrain data before depth correction.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811512074.8A CN110005396B (en) | 2018-12-11 | 2018-12-11 | Depth elastic correction method for logging wave train data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811512074.8A CN110005396B (en) | 2018-12-11 | 2018-12-11 | Depth elastic correction method for logging wave train data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110005396A CN110005396A (en) | 2019-07-12 |
CN110005396B true CN110005396B (en) | 2023-02-10 |
Family
ID=67165145
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811512074.8A Active CN110005396B (en) | 2018-12-11 | 2018-12-11 | Depth elastic correction method for logging wave train data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110005396B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113568044B (en) * | 2020-04-28 | 2023-09-26 | 中国石油天然气股份有限公司 | Method and device for determining head arrival time of array acoustic logging |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014185870A (en) * | 2013-03-22 | 2014-10-02 | Hazama Ando Corp | Elastic wave survey method |
CN104160299A (en) * | 2012-03-09 | 2014-11-19 | 雪佛龙美国公司 | Correction of shear log for elastic anisotropy |
CN106054254A (en) * | 2016-08-01 | 2016-10-26 | 中国石油天然气集团公司 | Well logging curve correction method and device |
CN106353813A (en) * | 2016-10-25 | 2017-01-25 | 中国海洋石油总公司 | Method for identifying fluid properties based on array acoustic logging |
CN106595834A (en) * | 2016-11-10 | 2017-04-26 | 西北工业大学 | Method of acquiring deep sea great depth sound field horizontal and longitudinal correlation |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6845325B2 (en) * | 2001-11-08 | 2005-01-18 | Schlumberger Technology Corporation | Global classification of sonic logs |
-
2018
- 2018-12-11 CN CN201811512074.8A patent/CN110005396B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104160299A (en) * | 2012-03-09 | 2014-11-19 | 雪佛龙美国公司 | Correction of shear log for elastic anisotropy |
JP2014185870A (en) * | 2013-03-22 | 2014-10-02 | Hazama Ando Corp | Elastic wave survey method |
CN106054254A (en) * | 2016-08-01 | 2016-10-26 | 中国石油天然气集团公司 | Well logging curve correction method and device |
CN106353813A (en) * | 2016-10-25 | 2017-01-25 | 中国海洋石油总公司 | Method for identifying fluid properties based on array acoustic logging |
CN106595834A (en) * | 2016-11-10 | 2017-04-26 | 西北工业大学 | Method of acquiring deep sea great depth sound field horizontal and longitudinal correlation |
Also Published As
Publication number | Publication date |
---|---|
CN110005396A (en) | 2019-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104730579B (en) | A kind of joint static correcting method of ripple in length and breadth based on calculation of near surface shear velocity inverting | |
CN109669212B (en) | Seismic data processing method, stratum quality factor estimation method and device | |
CN102928873B (en) | Method for positioning ground micro-seismic based on four-dimensional energy focusing | |
CN102565857B (en) | Automatic remaining dynamic correction method | |
CN102841373A (en) | Microseism positioning method based on azimuth angle constraint | |
CN101598808B (en) | Method for improving image quality of seismic data | |
CN101630016A (en) | Method for improving imaging quality of vertical seismic profile | |
CN107605470B (en) | Longitudinal and transverse wave radial velocity change imaging method | |
CN103645497B (en) | Emergence angle based multi-component wave field separation method | |
CN103792579B (en) | A kind of dynamic(al) correction method of suppressing nmo stretching | |
CN105182420A (en) | Dynamic-matching dynamic-correction method | |
CN104932016A (en) | Viterbi-BMC algorithm automatic speed analysis method | |
CN110005396B (en) | Depth elastic correction method for logging wave train data | |
CN105093319A (en) | Ground micro-seismic static correction method based on three-dimensional seismic data | |
CN105116448A (en) | Conversion wave azimuthal anisotropy correction method and device thereof | |
CN103399346A (en) | Well-to-seismic integration initial wave impedance modeling method | |
CN103576197A (en) | Method for extracting converted wave angle channel set | |
CN104182651B (en) | For the automatic quality control method in micro-seismic event azimuth that three-component geophone is received | |
CN102830170A (en) | Control method and control device for obtaining coal sample transverse wave signal based on ultrasonic test | |
CN102338885B (en) | Three-component VSP data first arrival time automatic pick method | |
CN107015273A (en) | A kind of coal seam Love type groove wave pressure compression methods and device based on frequency domain | |
CN102998700B (en) | A kind of earth's surface-consistent energy compensation process of vertical seismic profile data | |
CN104536048B (en) | A kind of automatic excision seismic wave ripple method extremely | |
CN112901146A (en) | Method and device for identifying bad track in acoustic logging detection | |
CN107797140B (en) | Beading seismic response features focus optimization method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |