CN112558159A - Acoustic logging first arrival picking method - Google Patents
Acoustic logging first arrival picking method Download PDFInfo
- Publication number
- CN112558159A CN112558159A CN202011424881.1A CN202011424881A CN112558159A CN 112558159 A CN112558159 A CN 112558159A CN 202011424881 A CN202011424881 A CN 202011424881A CN 112558159 A CN112558159 A CN 112558159A
- Authority
- CN
- China
- Prior art keywords
- arrival
- wave
- waveform
- data
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000001914 filtration Methods 0.000 claims abstract description 30
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 4
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 239000004576 sand Substances 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 239000011229 interlayer Substances 0.000 abstract description 6
- 238000011161 development Methods 0.000 abstract description 2
- 238000012545 processing Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000001615 p wave Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
Abstract
The invention discloses a sound wave logging first arrival picking method, which is characterized in that based on collected logging wave train data, on the basis of band-pass filtering and fuzzy nested multistage median filtering denoising, a dynamic self-adaptive first arrival picking range is obtained through time difference obtained by an STC method, and a waveform enveloping energy ratio method is adopted in the range to obtain a sound wave first arrival. The method can accurately pick up the sound wave first arrival of the stratum with low signal-to-noise ratio and serious attenuation and the stratum sections with thin interbed and interlayer development.
Description
Technical Field
The invention belongs to the field of geophysical logging, and relates to a sound wave logging first arrival picking method.
Background
In the process of processing acoustic data, accurate pickup of first arrival time is the basis of establishment of a rock physical model and is also an important parameter of longitudinal wave radial velocity imaging based on ray tracing. In recent years, many research results are obtained for predecessors of first-arrival picking methods, and georgel et al propose a mesh subdivision slowness model to perform first-arrival travel time calculation from four directions (CN 110568497A); the method comprises the following steps that a schargent slope and the like obtain a plurality of first arrivals through an energy ratio method, calculate first arrival confidence coefficients on the basis, and select the first arrival path with high confidence coefficient to perform fitting to obtain a first arrival predicted value (CN 102313901A); xukixiang provides a first-motion wave travel time picking method and a device, and the device extracts first-motion wave travel time (CN103616722A) in a processing time window of each seismic channel on a given initial reference line according to a calculated characteristic function value. Tomb et al have proposed a wave envelope-based anisotropic medium ultrasonic velocity automatic calculation method, obtain possible envelope first-break time according to solving its minimum value, set up the first-break screening condition of the wave packet, pick up the first-break time of the wave packet automatically (CN 106525979A); the method comprises the steps of providing a full-automatic seismic wave first arrival picking algorithm, calculating a reliability factor of a first arrival by analyzing the energy ratio difference between an accurate first arrival and an inaccurate first arrival, quantifying the change trend of the energy ratio, and setting reliability threshold value constraint to obtain the first arrival (CN 108072896A); li inspired to propose a long-short time mean ratio method for automatically picking up p-wave arrival times (CN201710842194.3), Heile proposed to adopt an improved energy ratio first arrival picking algorithm to calculate the first arrival wave position (CN110045417A) of each seismic channel in a template-based range band, Liuxin et al proposed a window fractional dimension method, Zhang et proposed a fractal box dimension method, Yuelong et al proposed a first arrival picking method based on time-frequency analysis; in recent years, with the popularity of artificial intelligence and intelligent image processing technology, artificial neural network methods and first-arrival wave picking methods based on image processing have become hot of research.
Comprehensively analyzing the existing sound wave first arrival picking algorithm, and finding that the current first arrival picking technology at least has the following problems: 1. the influence of signal-to-noise ratio is very obvious; 2. the accuracy of first arrival picking up is not high in stratum sections with serious attenuation of acoustic signals such as coal seams, diameter expansion and soft strata; 3. the difficulty in identifying stratum sections with low signal resolution such as thin interbed and interlayer is high.
Disclosure of Invention
The invention aims to provide a sound wave logging first arrival picking method which can more accurately pick sound waves of weak signal stratum sections with serious attenuation and stratum sections with low signal resolution such as thin interbed and interlayer on the basis of eliminating random noise and reducing the influence of signal to noise ratio.
The invention is realized by the following technical scheme:
a sonic logging first arrival pickup method, comprising:
s1, performing acoustic logging to obtain wave train data;
s2, comprehensively analyzing waveform distribution characteristics and frequency spectrum distribution characteristics of the wave train data, and determining the distribution range of the wave train effective signals in the wave train data in a time domain and a frequency domain;
s3, sequentially performing band-pass filtering and fuzzy nested multilevel median filtering on the frequency section of the wave train data in the frequency domain distribution range to eliminate noise signals;
s4, solving a sound wave time difference curve of the wave train data after the noise signal is eliminated by adopting an STC method; calculating a self-adaptive first arrival picking range based on the obtained sound wave time difference curve;
and S5, calculating a dynamic first arrival picking range on the basis of S4, and calculating a sound wave first arrival by adopting a waveform envelope energy ratio method in the dynamic first arrival picking range.
Preferably, in S1, the wave train logging is digital sonic logging, variable density sonic logging or array sonic logging.
Preferably, S2 is specifically: carrying out waveform distribution characteristic analysis on the wavetrain data, and determining the distribution range of the wavetrain effective signals in the wavetrain data in a time domain and the distribution range of longitudinal waves; and carrying out Fourier transform on the wavetrain data to obtain a spectrogram, analyzing the spectral distribution characteristics of the wavetrain data according to the spectrogram, and determining the distribution range of the wavetrain effective signals in the wavetrain data in a frequency domain.
Preferably, the fuzzy nested multistage median filtering in S3 is specifically: setting the waveform signal of the wave train data as x (t), Y (t) as a signal sequence of x (t) which is ordered from small to large, and alpha (t) as the middle value of Y (t);
the long filter length filtering results are:
the short filter length filtering results are:
the setting method of the threshold value comprises the following steps:
in the formula (3), m isThe total number of (2), C, represents the average of the amplitudes of the output of the wave train data after being filtered by the long filter.
Then fuzzy nested multistage median filtering outputs signals: when in useWhen the temperature of the water is higher than the set temperature,when in useWhen the temperature of the water is higher than the set temperature,
preferably, S4 finds the adaptive first arrival picking range, and the specific method is as follows:
S=Δt×(LR+n×LP)/ts-WF+LP×Δt×n/ts (4)
E=Δt×(LR+n×LP)/ts+WR+LP×Δt×n/ts (5)
s is a first arrival picking starting position; e is a first arrival picking end position; delta t is the sound wave time difference obtained by STC; l isRIs the source distance; l isPIs the receiving track pitch; n is the track number (0, 1.. 7); t is tsIs the sampling rate; wFThe energy of the front and rear time windows of the waveform envelope is longer than that of the front time window; wRThe energy of the front and rear time windows of the waveform envelope is longer than that of the rear time window.
Preferably, S5 specifically includes:
s5.1, calculating a dynamic first arrival picking range on the basis of S4;
s5.2, solving a waveform envelope signal of the wave train data in the dynamic first arrival picking range;
s5.3, solving the energy ratio of the time windows before and after the waveform envelope based on the waveform envelope signal;
and S5.4, solving the first wave first arrival TT of the wave train data based on the energy ratio of the time windows before and after the waveform envelope.
Further, S5.1 specifically is: and repeating the step S4 to obtain the dynamic first arrival picking range on the basis of the sound wave time difference obtained by the STC.
Further, S5.2 specifically is: let the waveform signal of the wave train data be x (t),hilbert transform for x (t);
the waveform envelope signal is:
further, S5.3 specifically is: the energy ratio of the time window before and after the waveform envelope is
Wherein the content of the first and second substances,is a relative energy of the wave train, where N is the number of points of the wave train.
Further, S5.2 specifically is: the energy ratio R of the time window before and after the envelope of the waveform is maximummaxThe time point t corresponding to the time plus the energy of the time window before and after the waveform envelope is longer than the time window before, namely the first arrival point, and the first arrival point is multiplied by the sampling rate tsAnd obtaining the first wave first arrival TT.
Compared with the prior art, the invention has the following beneficial technical effects:
on the basis of band-pass filtering and fuzzy nested multistage median filtering denoising, the sound wave time difference is obtained by adopting an STC method, a dynamic self-adaptive first arrival picking range is obtained according to the obtained sound wave time difference, and a waveform envelope energy ratio is adopted in the obtained range to obtain a sound wave first arrival. The waveform envelope energy ratio is based on Hilbert transform, the waveform envelope energy ratio can play a role in amplifying sound waves, the effect of first arrival on weak signal stratums with serious attenuation is better, but the influence of random noise is large, and the stability is poor. The STC method has good stability but low identification precision in weak signal stratum sections with serious attenuation, thin interbed sections and interlayer sections. The invention fully combines the advantages of the two, obtains the dynamic self-adaptive first arrival picking range through the sound wave time difference obtained by the STC method, obtains the first arrival by adopting the waveform envelope energy ratio on the basis, and reduces the artificial influence while simplifying the operation; in addition, the invention also adopts band-pass filtering and fuzzy nested multistage median filtering to preprocess the original sound wave signal, and eliminates the noise signal under the condition of keeping the effective signal as much as possible. Compared with the traditional waveform envelope energy ratio method, the method has better stability and small influence of signal-to-noise ratio, and compared with the STC method widely applied to actual processing and interpretation software, the method can more accurately pick up the sound wave first arrivals of weak signal stratum sections with serious attenuation, thin interbed sections with low resolution and high identification difficulty and interlayer development sections, and ensure that the sound wave first arrivals capable of truly reflecting corresponding stratum information are obtained.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 illustrates a distribution characteristic and a spectrum diagram of an acoustic waveform according to an embodiment of the present invention;
FIG. 3 illustrates a first arrival acoustic log pick according to an embodiment of the present invention.
Detailed Description
The present invention will now be described in further detail with reference to specific examples, which are intended to be illustrative, but not limiting, of the invention.
Based on the collected logging wave train data, on the basis of band-pass filtering and fuzzy nested multistage median filtering denoising, a dynamic self-adaptive first arrival picking range is obtained through a time difference obtained by an STC method, and a waveform enveloping energy ratio method is adopted in the range to obtain a sound wave first arrival. The method comprises the following steps:
performing wave train logging, such as digital acoustic logging, variable density acoustic logging, array acoustic logging and the like to obtain wave train data;
step two, carrying out comprehensive analysis on waveform distribution characteristics and frequency spectrum distribution characteristics of the wave train data, and determining the distribution range of the wave train effective signals in a time domain and the distribution range of a frequency domain;
thirdly, selecting a proper frequency section for band-pass filtering according to the distribution range of the effective signal of the wave train obtained in the second step in the frequency domain, and eliminating noise caused by frequency; carrying out fuzzy nested multistage median filtering on the wave train data obtained by the band-pass filtering, and eliminating noise signals under the condition of keeping effective signals as much as possible;
step four, solving a sound wave time difference curve of the denoised wave train data by adopting an STC (slowness-time correlation method); calculating a self-adaptive first arrival picking range based on the obtained sound wave time difference curve;
and step five, calculating a dynamic range of signal first arrival pickup on the basis of the step four, and calculating the sound wave first arrival by adopting a waveform envelope energy ratio method in the dynamic range.
Examples
As shown in fig. 1, the first arrival pickup method for acoustic logging of the present invention has the following specific implementation processes:
performing acoustic logging, such as digital acoustic logging, variable density acoustic logging, array acoustic logging and the like, to obtain wave train data;
secondly, performing waveform distribution characteristic analysis on the wave train data to determine the approximate distribution range of the wave train effective signals in a time domain and the approximate distribution range of longitudinal waves; and carrying out Fourier transform on the wavetrain data to obtain a spectrogram, and analyzing the spectral distribution characteristics of the wavetrain data to determine the distribution range of the wavetrain effective signals in the frequency domain. See figure 2.
Thirdly, performing band-pass filtering according to the distribution range of the effective wave train signals obtained in the second step in the frequency domain to eliminate noise signals caused by frequency; and carrying out fuzzy nested multistage median filtering on the obtained wave train data. The fuzzy nested multistage median filtering principle is that a threshold value is designed as a judgment parameter, so that a long filter is adopted for filtering when random noise is eliminated, and a short filter is adopted for filtering when effective information is reserved, so that the random noise can be eliminated well, and the effective information can be not damaged to the maximum extent. Let the waveform signal of the wave train data be x (t), Y (t) be the signal sequence of x (t) ordered from small to large, and alpha (t) be the middle value of Y (t).
The long filter length filtering results are:
the short filter length filtering results are:
the setting method of the threshold value comprises the following steps:
in the formula (3), m isThe total number of the cells. C represents the average of the amplitudes of the output of the wave train data after being filtered by the long filter.
Outputting a signal: when in useWhen the temperature of the water is higher than the set temperature,when in useWhen the temperature of the water is higher than the set temperature,
and fourthly, solving the sound wave time difference delta t of the denoised wave train data through an STC (slowness-time correlation method), and solving the self-adaptive first arrival picking range on the basis of the sound wave time difference delta t obtained by the STC method. The specific method comprises the following steps:
S=Δt×(LR+n×LP)/ts-WF+LP×Δt×n/ts (4)
E=Δt×(LR+n×LP)/ts+WR+LP×Δt×n/ts (5)
s is a first arrival picking starting position; e is a first arrival picking end position; delta t is the sound wave time difference obtained by an STC method; l isRIs the source distance; l isPIs the receiving track pitch; n is the track number (0, 1.. 7); t is tsIs the sampling rate; wFThe energy of the front and rear time windows of the waveform envelope is longer than that of the front time window; wRThe energy of the front and rear time windows of the waveform envelope is longer than that of the rear time window;
and step five, changing the traditional static picking range into a dynamic picking range by taking the first arrival calculated by the STC as a constraint, and calculating the first arrival time of the first wave by adopting a waveform envelope energy ratio method on the basis. The specific method comprises the following steps:
1) finding dynamic first arrival picking range
And (3) repeating the step (4) on the basis of the sound wave time difference obtained by the STC to obtain a dynamic first arrival picking range, changing the traditional static picking range into the dynamic picking range to restrict the picking of the first arrival, and solving the advantage of good stability of the sound wave time difference based on the STC method, wherein the dynamic range can greatly improve the problem of poor stability of the energy ratio method.
2) Evaluating a waveform envelope signal
Let the waveform signal of the wave train data be x (t),the hilbert transform is x (t), and an analysis signal formed by hilbert only contains positive frequency components and is 2 times of the positive frequency components of the original signal, so that the waveform amplitude envelope energy of the acoustic wave signal is stronger than the amplitude energy, the seismic wave signal actually plays a role in amplifying seismic waves, and the weak first arrival recognition effect is better.
The waveform envelope signal is:
3) energy ratio of front time window and rear time window of waveform envelope is calculated
The energy ratio of the time window before and after the waveform envelope is
4) Finding first arrival TT of array acoustic wave
The energy ratio R of the time window before and after the envelope of the waveform is maximummaxThe time point t corresponding to the time plus the energy of the time window before and after the waveform envelope is longer than the time window before, namely the first arrival point, and the first arrival point is multiplied by the sampling rate tsThe first arrival TT of the first wave can be obtained. Namely:
TT=(t+WF)*ts (8)
the result is shown in fig. 3, compared with the characteristic that the conventional waveform envelope energy ratio method is greatly influenced by random noise and has poor first-arrival stability, the method has better stability in the first-arrival picking of the whole well section; compared with the first arrivals obtained by an STC method widely applied to actual processing and interpretation software, the method has the advantage that the sound wave first arrival pick-up is more accurate in low-resolution stratum sections such as severely attenuated weak signal stratum sections (3607m-3613m), thin interbedders (3591m-3599m), and interlayers (3573m-3579 m).
The present invention is not limited to the above-mentioned embodiments, and based on the technical solutions disclosed in the present invention, those skilled in the art can make some substitutions and modifications to some technical features without creative efforts according to the disclosed technical contents, and these substitutions and modifications are all within the protection scope of the present invention.
Claims (10)
1. A sonic logging first arrival pickup method, comprising:
s1, performing acoustic logging to obtain wave train data;
s2, comprehensively analyzing waveform distribution characteristics and frequency spectrum distribution characteristics of the wave train data, and determining the distribution range of the wave train effective signals in the wave train data in a time domain and a frequency domain;
s3, sequentially performing band-pass filtering and fuzzy nested multilevel median filtering on the frequency section of the wave train data in the frequency domain distribution range to eliminate noise signals;
s4, solving the acoustic wave time difference of the wavetrain data after the noise signal is eliminated by adopting an STC method; calculating a self-adaptive first arrival picking range based on the obtained sound wave time difference;
and S5, calculating a dynamic first arrival picking range on the basis of S4, and calculating a sound wave first arrival by adopting a waveform envelope energy ratio method in the dynamic first arrival picking range.
2. The sonic logging first arrival pickup method of claim 1 wherein at S1 the wavetrain logging is digital sonic logging, variable density sonic logging or array sonic logging.
3. The sonic logging first arrival pickup method of claim 1, wherein S2 is embodied as: carrying out waveform distribution characteristic analysis on the wavetrain data, and determining the distribution range of the wavetrain effective signals in the wavetrain data in a time domain and the distribution range of longitudinal waves; and carrying out Fourier transform on the wavetrain data to obtain a spectrogram, analyzing the spectral distribution characteristics of the wavetrain data according to the spectrogram, and determining the distribution range of the wavetrain effective signals in the wavetrain data in a frequency domain.
4. The sonic logging first arrival pickup method of claim 1 wherein the fuzzy nested multistage median filtering in S3 is embodied as: setting the waveform signal of the wave train data as x (t), Y (t) as a signal sequence of x (t) which is ordered from small to large, and alpha (t) as the middle value of Y (t);
the long filter length filtering results are:
the short filter length filtering results are:
the setting method of the threshold value comprises the following steps:
in the formula (3), m isC represents the average value of the amplitude output after the wavetrain data is filtered by the long filter;
5. the sonic logging first arrival pick-up method of claim 1 wherein S4 finds an adaptive first arrival pick-up range by:
S=Δt×(LR+n×LP)/ts-WF+LP×Δt×n/ts (4)
E=Δt×(LR+n×LP)/ts+WR+LP×Δt×n/ts (5)
s is a first arrival picking starting position; e is a first arrival picking end position; delta t is the sound wave time difference obtained by STC; l isRIs the source distance; l isPIs the receiving track pitch; n is the track number (0, 1.)..7);tsIs the sampling rate; wFThe energy of the front and rear time windows of the waveform envelope is longer than that of the front time window; wRThe energy of the front and rear time windows of the waveform envelope is longer than that of the rear time window.
6. The sonic logging first arrival pickup method of claim 1 wherein S5 specifically comprises:
s5.1, calculating a dynamic first arrival picking range on the basis of S4;
s5.2, solving a waveform envelope signal of the wave train data in the dynamic first arrival picking range;
s5.3, solving the energy ratio of the time windows before and after the waveform envelope based on the waveform envelope signal;
and S5.4, solving the first wave first arrival TT of the wave train data based on the energy ratio of the time windows before and after the waveform envelope.
7. The acoustic logging first arrival pickup method according to claim 6, wherein S5.1 is specifically: and repeating the step S4 to obtain the dynamic first arrival picking range on the basis of the sound wave time difference obtained by the STC.
9. the acoustic logging first arrival pickup method according to claim 6, wherein S5.3 is specifically: the energy ratio of the time window before and after the waveform envelope is
10. The acoustic logging first arrival pickup method according to claim 6, wherein S5.2 is specifically: the energy ratio R of the time window before and after the envelope of the waveform is maximummaxThe time point t corresponding to the time plus the energy of the time window before and after the waveform envelope is longer than the time window before, namely the first arrival point, and the first arrival point is multiplied by the sampling rate tsAnd obtaining the first wave first arrival TT.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011424881.1A CN112558159A (en) | 2020-12-08 | 2020-12-08 | Acoustic logging first arrival picking method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011424881.1A CN112558159A (en) | 2020-12-08 | 2020-12-08 | Acoustic logging first arrival picking method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112558159A true CN112558159A (en) | 2021-03-26 |
Family
ID=75059842
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011424881.1A Pending CN112558159A (en) | 2020-12-08 | 2020-12-08 | Acoustic logging first arrival picking method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112558159A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11467302B1 (en) | 2021-05-26 | 2022-10-11 | Saudi Arabian Oil Company | Seismic first breaks onset times determination by interface tracking using level-sets method |
WO2023184055A1 (en) * | 2022-03-28 | 2023-10-05 | Saudi Arabian Oil Company | A method and system for determining first breaks of sonic waveforms |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5594706A (en) * | 1993-12-20 | 1997-01-14 | Schlumberger Technology Corporation | Downhole processing of sonic waveform information |
CN102243320A (en) * | 2011-04-12 | 2011-11-16 | 中国石油天然气股份有限公司 | Method for picking seismic wave first break |
CN103064111A (en) * | 2012-12-12 | 2013-04-24 | 中国石油天然气集团公司 | Micro seismic event recognition method based on morphological filtering |
CN106353812A (en) * | 2016-08-15 | 2017-01-25 | 黄河勘测规划设计有限公司 | Method of extracting of shear wave velocity of dual-spaced full-wave train logging based on slowness-time |
WO2017078741A1 (en) * | 2015-11-06 | 2017-05-11 | Halliburton Energy Services, Inc. | Methods and systems employing windowed frequency spectra analysis to derive a slowness log |
CN107272066A (en) * | 2017-06-22 | 2017-10-20 | 东华理工大学 | A kind of noisy seismic signal first-arrival traveltime pick-up method and device |
CN107678064A (en) * | 2016-08-02 | 2018-02-09 | 中石化石油工程技术服务有限公司 | A kind of interval transit time real time extracting method |
US20180196156A1 (en) * | 2017-01-10 | 2018-07-12 | Reeves Wireline Technologies Limited | Method of and Apparatus for Carrying Out Acoustic Well Logging |
CN108426949A (en) * | 2018-02-14 | 2018-08-21 | 国家海洋局第二海洋研究所 | A kind of seabed sediment acoustics original position data first arrival identification pick-up method |
US20190086571A1 (en) * | 2016-05-25 | 2019-03-21 | Halliburton Energy Services, Inc. | An improved stoneley wave slowness and dispersion curve logging method |
-
2020
- 2020-12-08 CN CN202011424881.1A patent/CN112558159A/en active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5594706A (en) * | 1993-12-20 | 1997-01-14 | Schlumberger Technology Corporation | Downhole processing of sonic waveform information |
CN102243320A (en) * | 2011-04-12 | 2011-11-16 | 中国石油天然气股份有限公司 | Method for picking seismic wave first break |
CN103064111A (en) * | 2012-12-12 | 2013-04-24 | 中国石油天然气集团公司 | Micro seismic event recognition method based on morphological filtering |
WO2017078741A1 (en) * | 2015-11-06 | 2017-05-11 | Halliburton Energy Services, Inc. | Methods and systems employing windowed frequency spectra analysis to derive a slowness log |
US20190086571A1 (en) * | 2016-05-25 | 2019-03-21 | Halliburton Energy Services, Inc. | An improved stoneley wave slowness and dispersion curve logging method |
CN107678064A (en) * | 2016-08-02 | 2018-02-09 | 中石化石油工程技术服务有限公司 | A kind of interval transit time real time extracting method |
CN106353812A (en) * | 2016-08-15 | 2017-01-25 | 黄河勘测规划设计有限公司 | Method of extracting of shear wave velocity of dual-spaced full-wave train logging based on slowness-time |
US20180196156A1 (en) * | 2017-01-10 | 2018-07-12 | Reeves Wireline Technologies Limited | Method of and Apparatus for Carrying Out Acoustic Well Logging |
CN107272066A (en) * | 2017-06-22 | 2017-10-20 | 东华理工大学 | A kind of noisy seismic signal first-arrival traveltime pick-up method and device |
CN108426949A (en) * | 2018-02-14 | 2018-08-21 | 国家海洋局第二海洋研究所 | A kind of seabed sediment acoustics original position data first arrival identification pick-up method |
Non-Patent Citations (9)
Title |
---|
JUNXIAO LI;KRISTOPHER A. INNANEN: "Extraction of reflected events from sonic-log waveforms using the Karhunen-Loève transform", 《GEOPHYSICS》, vol. 82, no. 5, 12 July 2017 (2017-07-12), pages 265 - 277 * |
姜文龙;张晓予;谢向文;涂善波: "基于慢度―时间相关的双源距全波列测井横波速度提取", 《水利水电技术》, vol. 47, no. 02, 29 February 2016 (2016-02-29), pages 115 - 118 * |
宋延杰, 宋传星, 石颖, 许秀英: "慢速地层偶极横波测井横波初至检测新方法", 《大庆石油学院学报》, vol. 29, no. 02, 30 April 2005 (2005-04-30), pages 7 - 9 * |
张伟: "阵列声波测井仪研制及测井数据处理方法研究", 《中国博士学位论文全文数据库信息科技辑》, no. 12, 15 December 2010 (2010-12-15), pages 1 - 115 * |
张伟;王海;李洪臣;王彦春;段云卿: "用变换时窗统计能量比法拾取地震初至波", 《物探与化探》, vol. 33, no. 02, 30 April 2009 (2009-04-30), pages 178 - 180 * |
王雷;周继宏: "多分辨率地层纵横波时差提取", 《工程地球物理学报》, vol. 14, no. 06, 30 November 2017 (2017-11-30), pages 686 - 693 * |
秦晅;宋维琪: "基于时窗能量比与互信息量的微地震初至拾取方法", 《物探与化探》, vol. 40, no. 02, 30 April 2016 (2016-04-30), pages 374 - 379 * |
胡文祥: "声波测井资料弱初至波检测新方法", 《江汉石油学院学报》, vol. 16, no. 1, 31 December 1994 (1994-12-31), pages 23 - 25 * |
郑江龙;许江: "能量差法拾取直达波初至的应用研究", 《工程地球物理学报》, vol. 10, no. 06, 30 November 2013 (2013-11-30), pages 834 - 839 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11467302B1 (en) | 2021-05-26 | 2022-10-11 | Saudi Arabian Oil Company | Seismic first breaks onset times determination by interface tracking using level-sets method |
WO2022250801A1 (en) * | 2021-05-26 | 2022-12-01 | Saudi Arabian Oil Company | Seismic first breaks onset times determination by interface tracking using level-sets method |
WO2023184055A1 (en) * | 2022-03-28 | 2023-10-05 | Saudi Arabian Oil Company | A method and system for determining first breaks of sonic waveforms |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2017024702A1 (en) | Inversion system for ray elastic parameter | |
CN111723329B (en) | Seismic phase feature recognition waveform inversion method based on full convolution neural network | |
CN110058305A (en) | A kind of DAS seismic data noise-reduction method based on convolutional neural networks | |
CN112558159A (en) | Acoustic logging first arrival picking method | |
CN107678064B (en) | Real-time extraction method for sound wave time difference | |
CN109283576B (en) | Method for automatically picking up seismic phase of P wave by taking amplitude as characteristic function | |
CN112780259B (en) | Method and device for determining well cementation quality and storage medium | |
CN116828354B (en) | Radio quality optimization method and system for collar clamp wireless microphone | |
CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
CN112593922A (en) | Method and device for evaluating cementing quality of two well cementation interfaces through array acoustic logging | |
CN107605470A (en) | A kind of wave path in length and breadth is to velocity variations imaging method | |
CN108490488A (en) | It is a kind of based on when window energy ratio channel wave seismic scattering imaging method and device | |
CN114488293A (en) | High-resolution inversion method based on sensitive logging curve | |
CN111538086B (en) | First arrival automatic pickup method for improving seismic data first arrival wave quality | |
CN107679614B (en) | Particle swarm optimization-based real-time sound wave time difference extraction method | |
CN110007342B (en) | Time-frequency domain direct pickup first arrival method and system for seismic signals with low signal-to-noise ratio | |
CN110879413A (en) | Ray parameter domain converted wave static correction method and system | |
CN110632660B (en) | Thin sand body characterization method and device based on seismic data body | |
CN112882097A (en) | Calibration method for highly deviated well and horizontal well | |
CN105425297A (en) | Method and device of suppressing ghost signals | |
CN111045083A (en) | Reservoir gas-containing property detection method | |
CN114035229B (en) | Pre-stack seismic data wavelet threshold denoising optimal wavelet base selection method | |
CN112925020B (en) | Parametric array type shallow profile intelligent energy compensation method | |
KR102512123B1 (en) | An method for picking first arrival using machine learning | |
CN113589374B (en) | Effective speed extraction method based on ray density |
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 |