CN106761715B - LWF memory-type acoustic logging first arrival detection method - Google Patents
LWF memory-type acoustic logging first arrival detection method Download PDFInfo
- Publication number
- CN106761715B CN106761715B CN201611239317.6A CN201611239317A CN106761715B CN 106761715 B CN106761715 B CN 106761715B CN 201611239317 A CN201611239317 A CN 201611239317A CN 106761715 B CN106761715 B CN 106761715B
- Authority
- CN
- China
- Prior art keywords
- mintrop wave
- wave
- negative peak
- mintrop
- point
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 65
- 230000002159 abnormal effect Effects 0.000 claims abstract description 188
- 238000000034 method Methods 0.000 claims abstract description 36
- 238000005070 sampling Methods 0.000 claims description 106
- 238000005553 drilling Methods 0.000 claims description 30
- 230000007257 malfunction Effects 0.000 claims description 20
- 238000012545 processing Methods 0.000 claims description 13
- 230000015572 biosynthetic process Effects 0.000 claims description 10
- 238000004891 communication Methods 0.000 claims description 8
- 239000000284 extract Substances 0.000 claims description 7
- 239000011435 rock Substances 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000003672 processing method Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims 1
- 238000004321 preservation Methods 0.000 claims 1
- 238000013461 design Methods 0.000 abstract description 4
- 230000035772 mutation Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000001737 promoting effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
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
- E21B47/12—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
- E21B47/14—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling using acoustic waves
-
- 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
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Mining & Mineral Resources (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a kind of LWF memory-type acoustic logging first arrival detection methods, comprising steps of one, acquisition well logging sound wave Wave data；Two, acoustic waveform is handled: 201, determining first arrival detection starting sample point sxRange, 202, determine first arrival detection starting sample point sx, 203, first arrival detection starting sample point sxRear negative peak position is extracted, and 204, determine Mintrop wave negative peak position, 205, complete the first arrival detection of well logging sound wave Wave data in Depth Domain, form Mintrop wave negative peak sequence ti；Three, on percentage regulation domain the well logging sound wave Wave data Mintrop wave abnormal point: 301, determine Mintrop wave sequence section, 302, the replacement of Mintrop wave abnormal point, 303, complete Mintrop wave abnormal point in Depth Domain and adjust process.Of the invention novel in design, the adjustment of time difference mutation carries out first arrival detection in the adjustment of wave distortion and longitudinal Depth Domain on binding time domain, method can in real time, acquisition stratum first wave of sound fast and automatically then, it is practical.
Description
Technical field
The present invention relates to sonic detection technical fields, specially LWF memory-type acoustic logging first arrival detection method.
Background technique
As the quickening of In Oil Field Exploration And Development process promotes, horizontal well, extended reach well, high angle hole and Holes of Complicated Wells increasingly increase
The problems such as more, bit speed also greatly promotes, hole condition complexity and hole anomaly seriously restricts traditional cable well logging and takes and take entirely
Quasi- data.Drill-pipe conveyed is the new technology of developed recently without wireline logging (LWF), it has abandoned the cable in traditional well logging, has had
Horizontal well, extended reach well, high angle hole and Holes of Complicated Wells are solved to effect using at high cost, time-consuming brought by drilling tool conveying well measuring
Long, the problems such as risk is big, it has also become the main logging method of horizontal well, extended reach well, high angle hole and Holes of Complicated Wells, but in reality
In the production of border, on the one hand, the hole condition of horizontal well, extended reach well and high angle hole is complicated, and drilling rod will shake left and right, logger without
Method is completely in center condition, will lead to the wave-shape amplitude that acoustic logging instrument receives and differs greatly, and noise amplitude is greater than
Longitudinal wave amplitude.On the other hand, since wellbore is big, instrument diameters are small, and emitted energy is small, cause in combination with factors such as reservoirs
The received wave distortion of acoustic logging instrument is serious, and influence of noise is very big, jump is serious.Meanwhile conventional arrangement Threshold detection is first
The method of wave, as a result serious distortion, most of situation must be using first arrival time be manually chosen, and processing speed is slow, efficiency
Low, accuracy rate is affected by human factors larger.Therefore, it is badly in need of studying adaptable first arrival detection method, improves production efficiency
And data reliability.
Summary of the invention
Aiming at the problems existing in the prior art, the present invention provides a kind of LWF memory-type acoustic logging first arrival detection method,
It is novel in design rationally, on binding time domain in the adjustment of wave distortion and longitudinal Depth Domain the adjustment of time difference mutation carry out it is first
Wave detection, method can in real time, acquisition stratum first wave of sound fast and automatically then, it is practical, convenient for promote the use of.
The present invention is to be achieved through the following technical solutions:
LWF memory-type acoustic logging first arrival detection method of the present invention, comprising the following steps:
Step 1 obtains well logging sound wave Wave data: in advance setting using acoustic logging instrument (1) and according to timer (2)
Surely sampling number N in the record time T and each record time T of waveform is received, each of well logging different depth position is obtained
The Wave data of time T is recorded, and is transmitted to microcontroller (3), the waveform for each record time T that microcontroller (3) will acquire
Data are stored in memory (4), and communication module (5) are connected on microcontroller (3), and N is positive integer；
The acoustic logging instrument (1) is mounted on the front end for the drilling rod extending into well logging wellbore, by being mounted on well
Computer (6) driving drilling rod dragging control device (7), drilling rod dragging control device (7) drives the drilling rod to move up and realizes institute
State the acquisition of well logging sound wave Wave data in Depth Domain in well logging wellbore；
Step 2, acoustic waveform processing: computer (6) obtains the waveform saved in memory (4) by communication module (5)
Data are respectively processed acoustic logging instrument in step 1 (1) in the Wave data that each record time T is obtained, each to record
The processing method for the Wave data that time T obtains is all the same；Respectively, it is first determined first arrival detection starting sample point sxModel
It encloses, then determines first arrival detection starting sample point sx, first arrival detection starting sample point s againxNegative peak position afterwards is extracted, and the 4th
It determines Mintrop wave negative peak position, finally completes the first arrival detection of well logging sound wave Wave data in Depth Domain, form Mintrop wave negative peak sequence
ti；
Step 3, the Mintrop wave abnormal point of the well logging sound wave Wave data on percentage regulation domain, it is first determined Mintrop wave sequence
Then section carries out the replacement of Mintrop wave abnormal point, finally complete Mintrop wave abnormal point in Depth Domain and adjust process.
Preferably, acoustic logging instrument described in step 1 (1) includes pinger and sends out for receiving the sound wave
The acoustic receiver of emitter is wireless acoustic signals, the pinger and the acoustic receiver are arranged at the drilling rod
On the connecting rod of front end, the acoustic receiver, which is located above the pinger, to be laid.
Preferably, in step 2, when handling any record time T Wave data obtained, detailed process is as follows:
Step 201, first arrival detection starting sample point s is determinedxRange: s1≤sx≤s2, wherein s1For in sampling number N
S1A sampled point, s2For s in sampling number N2A sampled point, sxFor between s1A sampled point and s2Between a sampled point
SxA sampled point, s1、sxAnd s2It is positive integer and s1≤s2≤ N, s1The corresponding sampling instant of a sampled point is
S2The corresponding sampling instant of a sampled point isSxThe corresponding sampling instant of a sampled point is
In the time-domain of record time T, the through waveform data obtained first in record time T, s are excluded1It is a
The sampled point acquired before sampled point is through waveform data, and waveform data sampling number amount of going directly is s1- 1；
Step 202, first arrival detection starting sample point s is determinedx, detailed process is as follows:
Step A1 selects sampled point sliding window n to carry out y (k)=a+bx (k) linear fit, wherein n is positive integer and 0
< n≤s2-s1+ 1, x (k) is k-th of sampled point corresponding time in sliding window n, and y (k) is the numerical value at x (k) sampled point, a
With b be sliding window n in linear fitting coefficient, k=1,2 ..., n；
Step A2, according to formulaCalculate fitting coefficient a in sliding window n
And b；
Step A3, according to formulaAll sampled points in sliding window n are calculated to deviate
Fitting a straight line y (k)=a+bx (k) variance var；
Step A1 is repeated several times in step B, until completing s1A sampled point is to s2All sliding window n on a sampled point
The variance var of interior sampled point；
Step C determines minimum variance varmin, computer (6) is compared variance var all in step B, minimum side
Poor varminThe corresponding sampled point of value be first arrival detection starting sample point sx；
Step 203, first arrival detection starting sample point sxExtract negative peak position afterwards: computer (6) originates first arrival detection
Sampled point sxFormer and later two sampling point values of n-th sampled point are compared in the record time T, when any sampled point
Numerical value when being respectively less than the numerical value of former and later two sampled points adjacent thereto, judge that this sampled point is negative peak dot, extracts all negative
Peak position is simultaneously stored in memory (4)；
Step 204, it determines Mintrop wave negative peak position: selecting the deviation threshold value Vs of Mintrop wave negative peak, step by computer (6)
First is more than that the negative peak position for deviateing threshold value Vs is determined as Mintrop wave negative peak position in all negative peak positions extracted in 203
It sets, the corresponding sampling instant in Mintrop wave negative peak position determined is t,
Step 205, step 201 is repeated several times, by drilling rod dragging control device (7) drive the drilling rod it is equidistant on
It moves, the next record time T Wave data obtained is handled, until completing the well logging sound wave Wave data in Depth Domain
First arrival detection, formed Mintrop wave negative peak sequence ti, whereinI=1,2 ..., m, m be Depth Domain on survey
Well position sequence and m are positive integer,L is depth-logger, and Δ l moves up equidistant for the drilling rod.
Further, the point of first arrival detection starting sample described in step 201 sxRange empirical value meet: 10≤sx≤45。
Further, the value of sampled point sliding window n described in step A1 is 9.
Further, the deviation threshold value Vs value of Mintrop wave negative peak described in step 204 is 180.
Further, detailed process is as follows for step 3:
Step 301, Mintrop wave sequence section is determined: by computer (6) to the sequence of Mintrop wave negative peak described in step 205 tiMiddle phase
The sampling instant of adjacent two Mintrop wave negative peaks carries out ratio calculation, the sampling instant ratio of at least three Mintrop wave negative peaks in continuous Depth Domain
Value γiMeet: 0.9≤γiWhen≤1.08, the correct Mintrop wave sequence section of sampling instant of each Mintrop wave is obtained, wherein
The quantity of the Mintrop wave sequence section is multiple, the number of the Mintrop wave abnormal point between the two Mintrop wave sequence sections that are connected
Amount is one or more, the sampling instant ratio γ of two adjacent Mintrop wave negative peaks in continuous Depth DomainiIt is unsatisfactory for: 0.9≤γi
When≤1.08, there are the Mintrop wave abnormal points under malfunction, and the Mintrop wave abnormal point is discrete in two adjacent Mintrop waves
Between sequence section；
The malfunction includes that expanding failure in wellbore, undergauge failure and stratum are collapsed failure, when expanding failure occurs
Or stratum collapse failure when, the transmission range of sound wave rock stratum into well logging wellbore is elongated, the delay of the Mintrop wave negative peak of Mintrop wave abnormal point
Occur, the sampling instant ratio γ of two Mintrop wave negative peaksi> 1.08；When undergauge failure occurs, sound wave rock into well logging wellbore
The transmission range of layer shortens, and the Mintrop wave negative peak of Mintrop wave abnormal point occurs in advance, the sampling instant ratio γ of two Mintrop wave negative peaksi<
0.9；
Step 302, the replacement of Mintrop wave abnormal point, process are as follows:
Step I judges the parity of the quantity of the Mintrop wave abnormal point between two Mintrop wave sequence sections: working as computer
(6) according to the quantity of the Mintrop wave abnormal point between two counted in the step 301 Mintrop wave sequence sections be even number when, hold
Row step II；When computer (6) is according to the Mintrop wave abnormal point between two counted in the step 301 Mintrop wave sequence sections
Quantity be odd number when, execute step III；
Step II, the equal part and classification processing of even number Mintrop wave abnormal point, process are as follows:
Step II a, the formation of two sequences and failure are sorted out in even number Mintrop wave abnormal point: by even number Mintrop wave abnormal point
Quantity is divided into two parts, and a part is continuous with a Mintrop wave sequence section in two Mintrop wave sequence sections and is included into two
First ray is formed in a Mintrop wave sequence section in the Mintrop wave sequence section；In another part and two Mintrop wave sequence sections
Another Mintrop wave sequence section is continuous and another Mintrop wave sequence section for being included into two Mintrop wave sequence sections in form second
Sequence, using the sampling instant ratio γ of two adjacent Mintrop wave negative peaksiCarry out failure classification；
The replacement of Mintrop wave negative peak: step II b when expanding failure occurs, is saved according in memory in step 203 (4)
All negative peak positions, the new negative peak position of Look-ahead, the new negative peak position in the time-domain locating for Mintrop wave abnormal point
Set the sampling instant ratio γ in Depth Domain with adjacent correct Mintrop wave negative peakiMeet: 0.9≤γiIt is described when≤1.08
New negative peak position is replaced Mintrop wave negative peak；When undergauge failure occurs, saved according in memory in step 203 (4)
All negative peak positions, search new negative peak position, the new negative peak position backward in the time-domain locating for Mintrop wave abnormal point
Set the sampling instant ratio γ in Depth Domain with adjacent correct Mintrop wave negative peakiMeet: 0.9≤γiIt is described when≤1.08
New negative peak position is replaced Mintrop wave negative peak；
Step II c, one or many to repeat step II a, the even number Mintrop wave between two Mintrop wave sequence sections
Abnormal point replacement finishes；
Step III, the classification processing of odd number Mintrop wave abnormal point, process are as follows:
Step III a, the formation and failure of Mintrop wave center abnormal point and two queues are sorted out in odd number Mintrop wave abnormal point:
Computer (6) is located in the Mintrop wave by searching for odd number Mintrop wave abnormal point median formation Mintrop wave center abnormal point is determined
The Mintrop wave abnormal point of heart abnormal point side is continuous with a Mintrop wave sequence section in two Mintrop wave sequence sections and is included into two
First queue is formed in a Mintrop wave sequence section in the Mintrop wave sequence section；Positioned at the Mintrop wave center abnormal point other side
Mintrop wave abnormal point is continuous with another Mintrop wave sequence section in two Mintrop wave sequence sections and is included into two Mintrop wave sequences
Second queue is formed in another Mintrop wave sequence section in section, using the sampling instant ratio γ of two adjacent Mintrop wave negative peaksi
Carry out failure classification；
Step III b, the replacement of Mintrop wave negative peak in two queues: when expanding failure occurs, according to memory in step 203
(4) all negative peak positions saved in, the new negative peak position of Look-ahead in the time-domain locating for Mintrop wave abnormal point are described new negative
Sampling instant ratio γ of the peak position in Depth Domain with adjacent correct Mintrop wave negative peakiMeet: 0.9≤γiWhen≤1.08,
The new negative peak position is replaced Mintrop wave negative peak；When undergauge failure occurs, protected according in memory in step 203 (4)
New negative peak position, the new negative peak position are searched in the time-domain locating for Mintrop wave abnormal point in all negative peak positions deposited backward
In Depth Domain with the sampling instant ratio γ of adjacent correct Mintrop wave negative peakiMeet: 0.9≤γiIt is described new when≤1.08
Negative peak position is replaced Mintrop wave negative peak；
Step III c, one or many to repeat step III a, the replacement up to completing Mintrop wave negative peak in two queues；
Step III d, the replacement of Mintrop wave center abnormal point calculate separately Mintrop wave center abnormal point by computer (6)
With the sampling instant ratio γ between two new negative peak positions adjacent theretoi, according in memory in step 203 (4)
The new negative peak moment is searched in the time-domain locating for the abnormal point of the Mintrop wave center, when described in all negative peak positions saved
The sampling instant ratio γ at new negative peak momentiMeet: 0.9≤γiWhen≤1.08, the new negative peak moment is the Mintrop wave
The replaced Mintrop wave negative peak of center abnormal point, the odd number Mintrop wave abnormal point replacement between two Mintrop wave sequence sections finish；
Step 303, one or many to repeat step 302, until completing Mintrop wave abnormal point in Depth Domain adjusts process.
Further, the quantity of the Mintrop wave abnormal point in step 301 between adjacent two Mintrop wave sequence sections is 1~5
It is a.
Preferably, record time T described in step 1 is that sampling number N is in 1792 μ s, each record time T
224.
Compared with prior art, the invention has the following beneficial technical effects:
1, first arrival detection device structure of the present invention is simple, and input cost is low, and it is convenient that installation is laid, convenient for promoting
It uses.
2, the sampled point that the present invention records in the time one in time-domain carries out first arrival detection starting sample point really
Recognize and all negative peak positions after the first arrival detection starting sample point are extracted, realizes the Mintrop wave negative peak in the record time
It extracts, high reliablity, effect is good.
3, the present invention carries out first arrival detection to the sampled point in each record time in well logging wellbore in Depth Domain, obtains
To the Mintrop wave negative peak sequence being initially formed, the initial Mintrop wave negative peak sequence of acquisition back can intuitively be checked by computer, it is real
When it is quick.
4, the present invention carries out the adjustment of Mintrop wave abnormal point to initial Mintrop wave negative peak sequence in Depth Domain, rejects due in wellbore
Expanding failure, undergauge failure and stratum collapse Mintrop wave negative peak caused by the influence of failure jump, improve first arrival detection accuracy rate, make
Obtained Mintrop wave negative peak sequence is smooth-out.
5, the method applied in the present invention step is simple, and design is rationally and realization is convenient, and precision is high, binding time domain upper ripple
The formation of Mintrop wave negative peak sequence is realized in the adjustment of time difference mutation in the adjustment of shape distortion and longitudinal Depth Domain, as a result accurately, just
In popularization and use.
In conclusion the present invention is novel in design rationally, the adjustment of wave distortion and longitudinal Depth Domain on binding time domain
The adjustment of upper time difference mutation carries out first arrival detection, method can in real time, acquisition stratum first wave of sound fast and automatically then, it is real
It is strong with property, convenient for promoting the use of.
Detailed description of the invention
Fig. 1 is the schematic block circuit diagram for the first arrival detection equipment that the present invention uses.
Fig. 2 is the method flow block diagram of first arrival detection method of the present invention.
Fig. 3 is sampled point and the branch for occurring to record sampled point in the time under expanding failure in the normal recordings time of the present invention
Schematic diagram.
Fig. 4 is sampled point and the branch for occurring to record sampled point in the time under undergauge failure in the normal recordings time of the present invention
Schematic diagram.
Fig. 5 is branch's schematic diagram that Mintrop wave negative peak sequence is initially formed in Depth Domain of the present invention.
Fig. 6 is branch's schematic diagram of Mintrop wave abnormal point Mintrop wave negative peak sequence adjusted in Depth Domain of the present invention.
In figure: 1 is acoustic logging instrument；2 be timer；3 be microcontroller；4 be memory；5 be communication module；6 are
Computer；7 drag control device for drilling rod.
Specific embodiment
Below with reference to specific embodiment, the present invention is described in further detail, it is described be explanation of the invention and
It is not to limit.
As depicted in figs. 1 and 2, LWF memory-type acoustic logging first arrival detection method of the present invention, comprising the following steps:
Step 1 obtains well logging sound wave Wave data: using acoustic logging instrument 1 and preset according to timer 2
Sampling number N in time T and each record time T is recorded, the wave of each record time T of well logging different depth position is obtained
Graphic data, and it is transmitted to microcontroller 3, the Wave data for each record time T that microcontroller 3 will acquire is stored in memory 4
In, communication module 5 is connected on microcontroller 3, N is positive integer；
In the present embodiment, record time T described in step 1 is sampling number N in 1792 μ s, each record time T
It is 224, the time interval between every two sampled point is 8 μ s；
The acoustic logging instrument 1 is mounted on the front end for the drilling rod extending into well logging wellbore, by being mounted on well
Computer 6 drives drilling rod to drag control device 7, and drilling rod dragging control device 7 drives the drilling rod to move up and realizes the well logging well
The acquisition of well logging sound wave Wave data in intraocular Depth Domain；
In actual use, acoustic logging instrument 1 described in step 1 is including pinger and for receiving the sound wave
The acoustic receiver of transmitter is wireless acoustic signals, the pinger and the acoustic receiver are arranged at the drilling rod
Front end connecting rod on, the acoustic receiver, which is located above the pinger, to be laid；
Step 2, acoustic waveform processing: computer 6 obtains the Wave data saved in memory 4 by communication module 5,
Acoustic logging instrument 1 in step 1 is respectively processed in the Wave data that each record time T is obtained, each time T that records is obtained
The processing method of the Wave data taken is all the same；When handling any record time T Wave data obtained, process is such as
Under:
Step 201, first arrival detection starting sample point s is determinedxRange: s1≤sx≤s2, wherein s1For in sampling number N
S1A sampled point, s2For s in sampling number N2A sampled point, sxFor between s1A sampled point and s2Between a sampled point
SxA sampled point, s1、sxAnd s2It is positive integer and s1≤s2≤ N, s1The corresponding sampling instant of a sampled point is
S2The corresponding sampling instant of a sampled point isSxThe corresponding sampling instant of a sampled point is
When practical operation, acoustic logging instrument 1 is set by timer 2 and is located at each record time T on depth-logger domain
Sampling time and every two sampled point between time interval, in conjunction with computer 6 periodically control drilling rod dragging control dress
Setting 7 moves up drilling rod, and acoustic logging instrument 1, which is located at the data acquired when underground and passes through microcontroller 3, is stored in memory 4
In, after DATA REASONING, the data in memory 4 are transmitted to by computer 6 by communication module 5 and carry out data processing；
In the present embodiment, the point of first arrival detection starting sample described in step 201 sxRange empirical value meet: 10≤sx
≤45；
In the time-domain of record time T, the through waveform data obtained first in record time T, s are excluded1It is a
The sampled point acquired before sampled point is through waveform data, and waveform data sampling number amount of going directly is s1- 1；
In the present embodiment, preceding 80 μ s is direct wave in the 1792 μ s of a record time T based on experience value, according to every two
Time interval between a sampled point is 8 μ s calculating, and before the 10th sampled point arrives, i.e., preceding 9 sampled points are direct wave, and
Occur between 80 μ of μ s~360 s in the 1792 μ s that longitudinal wave can record time T at one, i.e. first arrival detection starting sample point sxModel
It encloses for the 10th sampled point~45th sampled point；
Step 202, first arrival detection starting sample point s is determinedx, detailed process is as follows:
Step A1 selects sampled point sliding window n to carry out y (k)=a+bx (k) linear fit, wherein n is positive integer and 0
< n≤s2-s1+ 1, x (k) is k-th of sampled point corresponding time in sliding window n, and y (k) is the numerical value at x (k) sampled point, a
With b be sliding window n in linear fitting coefficient, k=1,2 ..., n；
In the present embodiment, the value of sampled point sliding window n described in step A1 is 9, for the number of the time interval of 8 μ s
Close to a wave period and at least will appear 2 zero points according to, 9 sampled points, the average value of a wave period close to zero curve,
Convenient for calculating variance；
Step A2, according to formulaCalculate fitting coefficient a in sliding window n
And b；
Step A3, according to formulaAll sampled points in sliding window n are calculated to deviate
Fitting a straight line y (k)=a+bx (k) variance var；
Step A1 is repeated several times in step B, until completing s1A sampled point is to s2All sliding window n on a sampled point
The variance var of interior sampled point；
Step C determines minimum variance varmin, computer 6 is compared variance var all in step B, minimum variance
varminThe corresponding sampled point of value be first arrival detection starting sample point sx；
In practical operation, sampled point in all sliding window n between the 10th sampled point~45th sampled point of calculating
Variance var is recycled 28 times and is found minimum variance varminThe corresponding sampled point of value be first arrival detection starting sample point sx；
Step 203, first arrival detection starting sample point sxExtract negative peak position afterwards: computer 6 adopts first arrival detection starting
Sampling point sxFormer and later two sampling point values of n-th sampled point are compared in the record time T, when any sampled point
When numerical value is respectively less than the numerical value of former and later two sampled points adjacent thereto, judge that this sampled point is negative peak dot, extracts all negative peaks
Position is simultaneously stored in memory 4；
Step 204, it determines Mintrop wave negative peak position: selecting the deviation threshold value Vs of Mintrop wave negative peak, step by computer 6
First is more than that the negative peak position for deviateing threshold value Vs is determined as Mintrop wave negative peak position in all negative peak positions extracted in 203
It sets, the corresponding sampling instant in Mintrop wave negative peak position determined is t,
In the present embodiment, the selection of the deviation threshold value Vs of Mintrop wave negative peak described in step 204 is determining Mintrop wave negative peak
Important parameter, if deviateed, threshold value Vs is too big, and fluctuating small Mintrop wave negative peak can be missed, if deviation threshold value Vs is too small,
It fluctuates big noise and is taken as Mintrop wave, the Mintrop wave negative peak error caused is excessive, big according to the first wave amplitude for counting whole mouth well
It is small, noise mean value and variance size are analyzed, choosing deviation threshold value Vs value is 180；
Step 205, step 201 is repeated several times, control device 7 is dragged by drilling rod, the drilling rod is driven equidistantly to move up,
The next record time T Wave data obtained is handled, until completing the well logging sound wave Wave data in Depth Domain
First arrival detection forms Mintrop wave negative peak sequence ti, whereinI=1,2 ..., m, m be Depth Domain on well logging
Position sequence and m are positive integer,L is depth-logger, and Δ l moves up equidistant for the drilling rod；
In actual use, depth-logger L be 2000m~3000m, the equidistant Δ l for using drilling rod to move up for 0.1m, m's
Value is 20001~30001, that is, the Mintrop wave negative peak sequence t formediIn have 20001~30001 Mintrop wave negative peaks；
Step 3, the Mintrop wave abnormal point of the well logging sound wave Wave data on percentage regulation domain, detailed process is as follows:
Step 301, Mintrop wave sequence section is determined: by computer 6 to the sequence of Mintrop wave negative peak described in step 205 tiIn it is adjacent
The sampling instant of two Mintrop wave negative peaks carries out ratio calculation, the sampling instant ratio of at least three Mintrop wave negative peaks in continuous Depth Domain
γiMeet: 0.9≤γiWhen≤1.08, the correct Mintrop wave sequence section of sampling instant of each Mintrop wave is obtained, wherein
It should be noted that when the Mintrop wave sequence section is the corresponding sampling of at least four Mintrop wave negative peaks in continuous Depth Domain
Carve the sampling instant ratio γ of the three Mintrop wave negative peaks formediIt is all satisfied 0.9≤γi≤ 1.08 sampling point set,Table
When showing that the sampling instant ratio of Mintrop wave negative peak is corresponding with the previous sampled point sampling of the corresponding sampling instant of latter sampled point
The ratio at quarter；
The quantity of the Mintrop wave sequence section is multiple, the number of the Mintrop wave abnormal point between the two Mintrop wave sequence sections that are connected
Amount is one or more, the sampling instant ratio γ of two adjacent Mintrop wave negative peaks in continuous Depth DomainiIt is unsatisfactory for: 0.9≤γi
When≤1.08, there are the Mintrop wave abnormal points under malfunction, and the Mintrop wave abnormal point is discrete in two adjacent Mintrop waves
Between sequence section；
In the present embodiment, the quantity of the Mintrop wave abnormal point in step 301 between the two neighboring Mintrop wave sequence section is 1~
5；
In conjunction with Fig. 3 and Fig. 4, in the present embodiment, the malfunction includes expanding failure in wellbore, undergauge failure and soft Jie
Matter failure, when expanding failure or soft-medium failure occurs, the transmission range of sound wave rock stratum into well logging wellbore is elongated, and Mintrop wave is different
The Mintrop wave negative peak delay often put occurs, the sampling instant ratio γ of two Mintrop wave negative peaksi> 1.08；When undergauge failure occurs,
Sound wave transmission range of rock stratum into well logging wellbore shortens, and the Mintrop wave negative peak of Mintrop wave abnormal point occurs in advance, two Mintrop wave negative peaks
Sampling instant ratio γi< 0.9；
Step 302, the replacement of Mintrop wave abnormal point, process are as follows:
Step I judges the parity of the quantity of the Mintrop wave abnormal point between two Mintrop wave sequence sections: when computer 6
According to the quantity of the Mintrop wave abnormal point between two counted in the step 301 Mintrop wave sequence sections be even number when, execute
Step II；When computer 6 is according to the number of the Mintrop wave abnormal point between two counted in the step 301 Mintrop wave sequence sections
When amount is odd number, step III is executed；
Step II, the equal part and classification processing of even number Mintrop wave abnormal point, process are as follows:
Step II a, the formation of two sequences and failure are sorted out in even number Mintrop wave abnormal point: by even number Mintrop wave abnormal point
Quantity is divided into two parts, and a part is continuous with a Mintrop wave sequence section in two Mintrop wave sequence sections and is included into two
First ray is formed in a Mintrop wave sequence section in the Mintrop wave sequence section；In another part and two Mintrop wave sequence sections
Another Mintrop wave sequence section is continuous and another Mintrop wave sequence section for being included into two Mintrop wave sequence sections in form second
Sequence, using the sampling instant ratio γ of two adjacent Mintrop wave negative peaksiCarry out failure classification；
Step II b, the replacement of Mintrop wave negative peak: when expanding failure occurs, according to what is saved in memory 4 in step 203
All negative peak positions, the new negative peak position of Look-ahead, the new negative peak position in the time-domain locating for Mintrop wave abnormal point
In Depth Domain with the sampling instant ratio γ of adjacent correct Mintrop wave negative peakiMeet: 0.9≤γiIt is described new when≤1.08
Negative peak position be replaced Mintrop wave negative peak；When undergauge failure occurs, according to the institute saved in memory 4 in step 203
There is negative peak position, search new negative peak position in the time-domain locating for Mintrop wave abnormal point backward, the new negative peak position exists
In Depth Domain with the sampling instant ratio γ of adjacent correct Mintrop wave negative peakiMeet: 0.9≤γiIt is described new when≤1.08
Negative peak position is replaced Mintrop wave negative peak；
Step II c, one or many to repeat step II a, the even number Mintrop wave between two Mintrop wave sequence sections
Abnormal point replacement finishes；
In the present embodiment, when the quantity of the Mintrop wave abnormal point between the Mintrop wave sequence section two neighboring in step 301 is 2
When a, 2 Mintrop wave abnormal points are respectively α1And α2, by Mintrop wave abnormal point α1It is included into two be contiguous therewith the Mintrop wave sequence
First ray is formed in a Mintrop wave sequence section in column section, by Mintrop wave abnormal point α2It is included into two Mintrop waves being contiguous therewith
The second sequence, Mintrop wave abnormal point α are formed in another Mintrop wave sequence section in sequence section1Corresponding sampling instant and First ray
In with Mintrop wave abnormal point α1The sampling instant of continuous correctly sampled point is compared, according to γiRatio result judge failure
State, in Mintrop wave abnormal point α1New Mintrop wave negative peak α is searched in locating time-domain1' replacement Mintrop wave abnormal point α1；Mintrop wave is abnormal
Point α2In corresponding sampling instant and the second sequence with Mintrop wave abnormal point α2The sampling instant of continuous correctly sampled point is compared
Compared with according to γiRatio result judge malfunction, in Mintrop wave abnormal point α2New Mintrop wave negative peak is searched in locating time-domain
α2' replacement Mintrop wave abnormal point α2, 2 Mintrop wave abnormal points, which are replaced, to be finished；
When the quantity of the Mintrop wave abnormal point between the Mintrop wave sequence section two neighboring in step 301 is 4,4 companies
Continuous Mintrop wave abnormal point is followed successively by β1、β2、β3And β4, by Mintrop wave abnormal point β1And β2It is included into two be contiguous therewith the Mintrop wave sequence
First ray is formed in a Mintrop wave sequence section in column section, by Mintrop wave abnormal point β3And β4It is included into described in two be contiguous therewith
The second sequence, Mintrop wave abnormal point β are formed in another Mintrop wave sequence section in Mintrop wave sequence section1Corresponding sampling instant and first
In sequence with Mintrop wave abnormal point β1The sampling instant of continuous correctly sampled point is compared, according to γiRatio result judgement
Malfunction, in Mintrop wave abnormal point β1New Mintrop wave negative peak β is searched in locating time-domain1' replacement Mintrop wave abnormal point β1, at this point,
New Mintrop wave negative peak β1' correct Mintrop wave negative peak is classified as in First ray；Mintrop wave abnormal point β2Corresponding sampling instant and the first sequence
New Mintrop wave negative peak β in column1' sampling instant be compared, according to γiRatio result judge malfunction, Mintrop wave exception
Point β2New Mintrop wave negative peak β is searched in locating time-domain2' replacement Mintrop wave abnormal point β2；Mintrop wave abnormal point β4When corresponding sampling
Carve and the second sequence in Mintrop wave abnormal point β4The sampling instant of continuous correctly sampled point is compared, according to γiRatio
As a result malfunction is judged, in Mintrop wave abnormal point β4New Mintrop wave negative peak β is searched in locating time-domain4' replacement Mintrop wave abnormal point
β4, at this point, new Mintrop wave negative peak β4' correct Mintrop wave negative peak is classified as in the second sequence；Mintrop wave abnormal point β3Corresponding sampling instant
With Mintrop wave negative peak β new in the second sequence4' sampling instant be compared, according to γiRatio result judge malfunction,
Mintrop wave abnormal point β3New Mintrop wave negative peak β is searched in locating time-domain3' replacement Mintrop wave abnormal point β3, 4 Mintrop wave abnormal points replace
It changes and finishes；
Step III, the classification processing of odd number Mintrop wave abnormal point, process are as follows:
Step III a, the formation and failure of Mintrop wave center abnormal point and two queues are sorted out in odd number Mintrop wave abnormal point:
Computer 6 is located at the Mintrop wave center by searching for odd number Mintrop wave abnormal point median formation Mintrop wave center abnormal point is determined
The Mintrop wave abnormal point of abnormal point side is continuous with a Mintrop wave sequence section in two Mintrop wave sequence sections and is included into two institutes
It states in a Mintrop wave sequence section in Mintrop wave sequence section and forms first queue；Head positioned at the Mintrop wave center abnormal point other side
Wave abnormal point is continuous with another Mintrop wave sequence section in two Mintrop wave sequence sections and is included into two Mintrop wave sequence sections
In another Mintrop wave sequence section in form second queue, using the sampling instant ratio γ of two adjacent Mintrop wave negative peaksiInto
Row failure is sorted out；
Step III b, the replacement of Mintrop wave negative peak in two queues: when expanding failure occurs, according to memory in step 203
All negative peak positions saved in 4, the new negative peak position of Look-ahead in the time-domain locating for Mintrop wave abnormal point, the new negative peak
Sampling instant ratio γ of the position in Depth Domain with adjacent correct Mintrop wave negative peakiMeet: 0.9≤γiWhen≤1.08, institute
Stating new negative peak position is replaced Mintrop wave negative peak；When undergauge failure occurs, according to what is saved in memory 4 in step 203
New negative peak position is searched in the time-domain locating for Mintrop wave abnormal point in all negative peak positions backward, and the new negative peak position is in depth
Spend the sampling instant ratio γ on domain with adjacent correct Mintrop wave negative peakiMeet: 0.9≤γiWhen≤1.08, the new negative peak
Position is replaced Mintrop wave negative peak；
Step III c, one or many to repeat step III a, the replacement up to completing Mintrop wave negative peak in two queues；
Step III d, the replacement of Mintrop wave center abnormal point, by computer 6 calculate separately Mintrop wave center abnormal point with
Sampling instant ratio γ between the new negative peak position of adjacent thereto twoi, according to being saved in memory 4 in step 203
All negative peak positions, the new negative peak moment is searched in the time-domain locating for the abnormal point of the Mintrop wave center, when described new
The sampling instant ratio γ at negative peak momentiMeet: 0.9≤γiWhen≤1.08, the new negative peak moment is the Mintrop wave center
The replaced Mintrop wave negative peak of abnormal point, the odd number Mintrop wave abnormal point replacement between two Mintrop wave sequence sections finish；
In the present embodiment, when the quantity of the Mintrop wave abnormal point between the Mintrop wave sequence section two neighboring in step 301 is 1
When a, 1 Mintrop wave abnormal point is δ, using two Mintrop wave sequence sections neutralize Mintrop wave abnormal point δ adjacent two it is correct
Mintrop wave negative peak respectively compared with Mintrop wave abnormal point δ, judge malfunction, adjust the front and back sampling instant ratio of Mintrop wave abnormal point δ
Value γi, new Mintrop wave negative peak δ ' replacement Mintrop wave abnormal point δ, 1 Mintrop wave exception are searched in the time-domain locating for Mintrop wave abnormal point δ
Point replacement finishes；
When the quantity of the Mintrop wave abnormal point between the Mintrop wave sequence section two neighboring in step 301 is 3,3 companies
Continuous Mintrop wave abnormal point is followed successively by ε1、ε2And ε3, Mintrop wave center abnormal point ε is searched by computer 62, by Mintrop wave abnormal point ε1Return
Enter in a Mintrop wave sequence section in two be contiguous therewith the Mintrop wave sequence section and form first queue, by Mintrop wave abnormal point ε3
Second queue, Mintrop wave abnormal point are formed in another Mintrop wave sequence section being included into two be contiguous therewith the Mintrop wave sequence section
ε1In corresponding sampling instant and first queue with Mintrop wave abnormal point ε1The sampling instant of continuous correctly sampled point is compared
Compared with according to γiRatio result judge malfunction, in Mintrop wave abnormal point ε1New Mintrop wave negative peak ε is searched in locating time-domain1'
Replace Mintrop wave abnormal point ε1, at this point, new Mintrop wave negative peak ε1' correct Mintrop wave negative peak is classified as in first queue；Mintrop wave abnormal point ε3It is right
In the sampling instant and second queue answered with Mintrop wave abnormal point ε3The sampling instant of continuous correctly sampled point is compared, root
According to γiRatio result judge malfunction, in Mintrop wave abnormal point ε3New Mintrop wave negative peak ε is searched in locating time-domain3' replacement
Mintrop wave abnormal point ε3, at this point, new Mintrop wave negative peak ε3' correct Mintrop wave negative peak is classified as in second queue；New Mintrop wave negative peak ε1' and it is new
Mintrop wave negative peak ε3' respectively with Mintrop wave center abnormal point ε2It is compared, judges malfunction, meet γiRatio range when,
Mintrop wave center abnormal point ε2New Mintrop wave negative peak ε is searched in locating time-domain2' replacement Mintrop wave abnormal point ε2, 3 Mintrop wave abnormal points
Replacement finishes；
When the quantity of the Mintrop wave abnormal point between the Mintrop wave sequence section two neighboring in step 301 is 5,5 companies
Continuous Mintrop wave abnormal point is followed successively by φ1、φ2、φ3、φ4And φ5, φ1'φ2'φ3'φ4'φ5' pass through the lookup Mintrop wave of computer 6
Center abnormal point φ3, by Mintrop wave abnormal point φ1And φ2A Mintrop wave being included into two be contiguous therewith the Mintrop wave sequence section
First queue is formed in sequence section, by Mintrop wave abnormal point φ4And φ5It is included into two be contiguous therewith the Mintrop wave sequence section
Second queue, Mintrop wave abnormal point φ are formed in another Mintrop wave sequence section1In corresponding sampling instant and first queue with Mintrop wave
Abnormal point φ1The sampling instant of continuous correctly sampled point is compared, according to γiRatio result judge malfunction,
Mintrop wave abnormal point φ1New Mintrop wave negative peak φ is searched in locating time-domain1' replacement Mintrop wave abnormal point φ1, at this point, new Mintrop wave negative peak
φ1' correct Mintrop wave negative peak is classified as in first queue；Mintrop wave abnormal point φ2New Mintrop wave in corresponding sampling instant and first queue
Negative peak φ1' sampling instant be compared, according to γiRatio result judge malfunction, in Mintrop wave abnormal point φ2Locating
New Mintrop wave negative peak φ is searched in time-domain2' replacement Mintrop wave abnormal point φ2；Mintrop wave abnormal point φ5Corresponding sampling instant and second
In queue with Mintrop wave abnormal point φ5The sampling instant of continuous correctly sampled point is compared, according to γiRatio result sentence
Disconnected malfunction, in Mintrop wave abnormal point φ5New Mintrop wave negative peak φ is searched in locating time-domain5' replacement Mintrop wave abnormal point φ5, this
When, new Mintrop wave negative peak φ5' correct Mintrop wave negative peak is classified as in second queue；Mintrop wave abnormal point φ4Corresponding sampling instant and
New Mintrop wave negative peak φ in two queues5' sampling instant be compared, according to γiRatio result judge malfunction, in Mintrop wave
Abnormal point φ4New Mintrop wave negative peak φ is searched in locating time-domain4' replacement Mintrop wave abnormal point φ4；New Mintrop wave negative peak φ2' and it is new
Mintrop wave negative peak φ4' respectively with Mintrop wave center abnormal point φ3It is compared, judges malfunction, meet γiRatio range when,
In Mintrop wave center abnormal point φ3New Mintrop wave negative peak φ is searched in locating time-domain3' replacement Mintrop wave abnormal point φ3, 5 Mintrop waves are different
Often point replacement finishes；
Step 303, one or many to repeat step 302, until completing Mintrop wave abnormal point in Depth Domain adjusts process.
In conjunction with Fig. 5 and Fig. 6, in the present embodiment, since the quantity of the Mintrop wave sequence section is multiple, Mintrop wave described in every two
The discrete Mintrop wave abnormal point processing method occurred between sequence section with the process phase of the replacement of Mintrop wave abnormal point in step 302
Together, this will not be repeated here, and for this method first arrival detection processing result as shown in fig. 6, processing is smooth, effect is good.
The above is only presently preferred embodiments of the present invention, is not intended to limit the invention in any way, it is all according to the present invention
Technical spirit any simple modification to the above embodiments, change and equivalent structural changes, still fall within skill of the present invention
In the protection scope of art scheme.
Claims (9)
1.LWF memory-type acoustic logging first arrival detection method, which comprises the following steps:
Step 1 obtains well logging sound wave Wave data: presetting using acoustic logging instrument (1) and according to timer (2) and connects
Sampling number N in the record time T and each record time T of waveform is received, each record of well logging different depth position is obtained
The Wave data of time T, and microcontroller (3) are transmitted to, the Wave data for each record time T that microcontroller (3) will acquire
It is stored in memory (4), communication module (5) is connected on microcontroller (3), N is positive integer；
The acoustic logging instrument (1) is mounted on the front end for the drilling rod extending into well logging wellbore, passes through the meter being mounted on well
Calculation machine (6) drives drilling rod dragging control device (7), and drilling rod dragging control device (7) drives the drilling rod to move up and realizes the survey
In well wellbore in Depth Domain well logging sound wave Wave data acquisition；
Step 2, acoustic waveform processing: computer (6) obtains the waveform number saved in memory (4) by communication module (5)
According to the Wave data to acoustic logging instrument in step 1 (1) in each record time T acquisition is respectively processed, respectively when record
Between T obtain Wave data processing method it is all the same；Respectively, it is first determined first arrival detection starting sample point sxRange,
Then first arrival detection starting sample point s is determinedx, first arrival detection starting sample point s againxNegative peak position afterwards is extracted, and the 4th determines
Mintrop wave negative peak position finally completes the first arrival detection of well logging sound wave Wave data in Depth Domain, forms Mintrop wave negative peak sequence ti；
Step 3, the Mintrop wave abnormal point of the well logging sound wave Wave data on percentage regulation domain, it is first determined Mintrop wave sequence section, so
The replacement for carrying out Mintrop wave abnormal point afterwards finally completes Mintrop wave abnormal point in Depth Domain and adjusts process.
2. LWF memory-type acoustic logging first arrival detection method according to claim 1, it is characterised in that: institute in step 1
Stating acoustic logging instrument (1) includes pinger and the acoustic receiver for receiving the wireless acoustic signals of the pinger
On the connecting rod for the front end that device, the pinger and the acoustic receiver are arranged at the drilling rod, the acoustic receiver
Device, which is located above the pinger, to be laid.
3. LWF memory-type acoustic logging first arrival detection method according to claim 1, it is characterised in that: right in step 2
When the Wave data that any record time T is obtained is handled, detailed process is as follows:
Step 201, first arrival detection starting sample point s is determinedxRange: s1≤sx≤s2, wherein s1For s in sampling number N1
A sampled point, s2For s in sampling number N2A sampled point, sxFor between s1A sampled point and s2Between a sampled point
sxA sampled point, s1、sxAnd s2It is positive integer and s1≤s2≤ N, s1The corresponding sampling instant of a sampled point isS2
The corresponding sampling instant of a sampled point isSxThe corresponding sampling instant of a sampled point is
In the time-domain of record time T, the through waveform data obtained first in record time T, s are excluded1A sampled point
The sampled point acquired before is through waveform data, and waveform data sampling number amount of going directly is s1- 1；
Step 202, first arrival detection starting sample point s is determinedx, detailed process is as follows:
Step A1 selects sampled point sliding window n to carry out y (k)=a+bx (k) linear fit, wherein n is positive integer and 0 < n
≤s2-s1+ 1, x (k) is k-th of sampled point corresponding time in sliding window n, and y (k) is the numerical value at x (k) sampled point, a and b
For linear fitting coefficient in sliding window n, k=1,2 ..., n；
Step A2, according to formulaCalculate fitting coefficient a and b in sliding window n；
Step A3, according to formulaIt calculates all sampled points in sliding window n and deviates fitting
Straight line y (k)=a+bx (k) variance var；
Step A1 is repeated several times in step B, until completing s1A sampled point is to s2It is adopted in all sliding window n on a sampled point
The variance var of sampling point；
Step C determines minimum variance varmin, computer (6) is compared variance var all in step B, minimum variance
varminThe corresponding sampled point of value be first arrival detection starting sample point sx；
Step 203, first arrival detection starting sample point sxExtract negative peak position afterwards: computer (6) is to first arrival detection starting sample point
sxFormer and later two sampling point values of n-th sampled point are compared in the record time T, when the numerical value of any sampled point
When the numerical value of former and later two sampled points respectively less than adjacent thereto, judge that this sampled point is negative peak dot, extracts all negative peak positions
And it is stored in memory (4)；
Step 204, it determines Mintrop wave negative peak position: selecting the deviation threshold value Vs of Mintrop wave negative peak, step 203 by computer (6)
First is more than that the negative peak position for deviateing threshold value Vs is determined as Mintrop wave negative peak position in all negative peak positions of middle extraction,
The determining corresponding sampling instant in Mintrop wave negative peak position is t,
Step 205, step 201 is repeated several times, drives the drilling rod equidistantly to move up by drilling rod dragging control device (7), it is right
The Wave data that next record time T is obtained is handled, until completing the head of the well logging sound wave Wave data in Depth Domain
Wave detection, forms Mintrop wave negative peak sequence ti, whereinI=1,2 ..., m, m be Depth Domain on well logging position
It sets sequence and m is positive integer,L is depth-logger, and Δ l moves up equidistant for the drilling rod.
4. LWF memory-type acoustic logging first arrival detection method described in accordance with the claim 3, it is characterised in that: institute in step 201
State first arrival detection starting sample point sxRange empirical value meet: 10≤sx≤45。
5. LWF memory-type acoustic logging first arrival detection method described in accordance with the claim 3, it is characterised in that: institute in step A1
The value for stating sampled point sliding window n is 9.
6. LWF memory-type acoustic logging first arrival detection method described in accordance with the claim 3, it is characterised in that: institute in step 204
The deviation threshold value Vs value for stating Mintrop wave negative peak is 180.
7. LWF memory-type acoustic logging first arrival detection method according to claim 3, it is characterised in that: the tool of step 3
Body process is as follows:
Step 301, Mintrop wave sequence section is determined: by computer (6) to the sequence of Mintrop wave negative peak described in step 205 tiIn adjacent two
The sampling instant of a Mintrop wave negative peak carries out ratio calculation, the sampling instant ratio of at least three Mintrop wave negative peaks in continuous Depth Domain
γiMeet: 0.9≤γiWhen≤1.08, the correct Mintrop wave sequence section of sampling instant of each Mintrop wave is obtained, wherein
The quantity of the Mintrop wave sequence section be it is multiple, the quantity of Mintrop wave abnormal point between the two Mintrop wave sequence sections that are connected is
One or more, the sampling instant ratio γ of two adjacent Mintrop wave negative peaks in continuous Depth DomainiIt is unsatisfactory for: 0.9≤γi≤
When 1.08, there are the Mintrop wave abnormal points under malfunction, and the Mintrop wave abnormal point is discrete in two adjacent Mintrop wave sequences
Between column section；
The malfunction includes that expanding failure in wellbore, undergauge failure and stratum are collapsed failure, when expanding failure or ground occurs
Layer collapse failure when, the transmission range of sound wave rock stratum into well logging wellbore is elongated, and the Mintrop wave negative peak of Mintrop wave abnormal point postpones to occur,
The sampling instant ratio γ of two Mintrop wave negative peaksi> 1.08；When undergauge failure occurs, the biography of sound wave rock stratum into well logging wellbore
Defeated distance shortens, and the Mintrop wave negative peak of Mintrop wave abnormal point occurs in advance, the sampling instant ratio γ of two Mintrop wave negative peaksi< 0.9；
Step 302, the replacement of Mintrop wave abnormal point, process are as follows:
Step I judges the parity of the quantity of the Mintrop wave abnormal point between two Mintrop wave sequence sections: when computer (6) root
According to the Mintrop wave abnormal point between two counted in the step 301 Mintrop wave sequence sections quantity be even number when, execute step
Rapid II；When computer (6) is according to the number of the Mintrop wave abnormal point between two counted in the step 301 Mintrop wave sequence sections
When amount is odd number, step III is executed；
Step II, the equal part and classification processing of even number Mintrop wave abnormal point, process are as follows:
Step II a, the formation of two sequences and failure are sorted out in even number Mintrop wave abnormal point: by even number Mintrop wave abnormal point quantity
Be divided into two parts, a part is continuous with a Mintrop wave sequence sections in two Mintrop wave sequence sections and is included into described in two
First ray is formed in a Mintrop wave sequence section in Mintrop wave sequence section；Another part with it is another in two Mintrop wave sequence sections
One Mintrop wave sequence section is continuous and another Mintrop wave sequence section for being included into two Mintrop wave sequence sections in form the second sequence,
Using the sampling instant ratio γ of two adjacent Mintrop wave negative peaksiCarry out failure classification；
Step II b, the replacement of Mintrop wave negative peak: when expanding failure occurs, according to the institute saved in memory in step 203 (4)
There is negative peak position, the new negative peak position of Look-ahead in the time-domain locating for Mintrop wave abnormal point, the new negative peak position exists
In Depth Domain with the sampling instant ratio γ of adjacent correct Mintrop wave negative peakiMeet: 0.9≤γiIt is described new when≤1.08
Negative peak position is replaced Mintrop wave negative peak；When undergauge failure occurs, according to the institute saved in memory in step 203 (4)
There is negative peak position, search new negative peak position in the time-domain locating for Mintrop wave abnormal point backward, the new negative peak position exists
In Depth Domain with the sampling instant ratio γ of adjacent correct Mintrop wave negative peakiMeet: 0.9≤γiIt is described new when≤1.08
Negative peak position is replaced Mintrop wave negative peak；
Step II c, one or many repetition step II a, until the even number Mintrop wave between two Mintrop wave sequence sections is abnormal
Point replacement finishes；
Step III, the classification processing of odd number Mintrop wave abnormal point, process are as follows:
Step III a, the formation and failure of Mintrop wave center abnormal point and two queues are sorted out in odd number Mintrop wave abnormal point: calculating
It is different to be located at the Mintrop wave center by searching for odd number Mintrop wave abnormal point median formation Mintrop wave center abnormal point is determined for machine (6)
Often the Mintrop wave abnormal point of point side is continuous with a Mintrop wave sequence section in two Mintrop wave sequence sections and is included into described in two
First queue is formed in a Mintrop wave sequence section in Mintrop wave sequence section；Mintrop wave positioned at the Mintrop wave center abnormal point other side
Abnormal point is continuous with another Mintrop wave sequence section in two Mintrop wave sequence sections and is included into two Mintrop wave sequence sections
Another Mintrop wave sequence section in form second queue, using the sampling instant ratio γ of two adjacent Mintrop wave negative peaksiIt carries out
Failure is sorted out；
Step III b, the replacement of Mintrop wave negative peak in two queues: when expanding failure occurs, according to memory in step 203 (4)
All negative peak positions of middle preservation, the new negative peak position of Look-ahead in the time-domain locating for Mintrop wave abnormal point, the new negative peak
Sampling instant ratio γ of the position in Depth Domain with adjacent correct Mintrop wave negative peakiMeet: 0.9≤γiWhen≤1.08, institute
Stating new negative peak position is replaced Mintrop wave negative peak；When undergauge failure occurs, saved according in memory in step 203 (4)
All negative peak positions, search new negative peak position backward in the time-domain locating for Mintrop wave abnormal point, the new negative peak position exists
In Depth Domain with the sampling instant ratio γ of adjacent correct Mintrop wave negative peakiMeet: 0.9≤γiIt is described new negative when≤1.08
Peak position is set to replaced Mintrop wave negative peak；
Step III c, one or many to repeat step III a, the replacement up to completing Mintrop wave negative peak in two queues；
Step III d, the replacement of Mintrop wave center abnormal point, by computer (6) calculate separately Mintrop wave center abnormal point with
Sampling instant ratio γ between its two adjacent new negative peak positioni, saved according in memory in step 203 (4)
All negative peak positions, the new negative peak moment is searched in the time-domain locating for the abnormal point of the Mintrop wave center, when described new
The sampling instant ratio γ at negative peak momentiMeet: 0.9≤γiWhen≤1.08, the new negative peak moment is the Mintrop wave center
The replaced Mintrop wave negative peak of abnormal point, the odd number Mintrop wave abnormal point replacement between two Mintrop wave sequence sections finish；
Step 303, one or many to repeat step 302, until completing Mintrop wave abnormal point in Depth Domain adjusts process.
8. LWF memory-type acoustic logging first arrival detection method according to claim 7, it is characterised in that: phase in step 301
The quantity for facing the Mintrop wave abnormal point between two Mintrop wave sequence sections is 1~5.
9. LWF memory-type acoustic logging first arrival detection method according to claim 1, it is characterised in that: institute in step 1
Stating record time T is that sampling number N is 224 in 1792 μ s, each record time T.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611239317.6A CN106761715B (en) | 2016-12-28 | 2016-12-28 | LWF memory-type acoustic logging first arrival detection method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611239317.6A CN106761715B (en) | 2016-12-28 | 2016-12-28 | LWF memory-type acoustic logging first arrival detection method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106761715A CN106761715A (en) | 2017-05-31 |
| CN106761715B true CN106761715B (en) | 2019-07-05 |
Family
ID=58923304
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611239317.6A Active CN106761715B (en) | 2016-12-28 | 2016-12-28 | LWF memory-type acoustic logging first arrival detection method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106761715B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110617859B (en) * | 2019-10-28 | 2021-04-06 | 贵州航天凯山石油仪器有限公司 | Transducer head wave threshold calculation method based on direct-reading ultrasonic flowmeter |
| CN113568044B (en) * | 2020-04-28 | 2023-09-26 | 中国石油天然气股份有限公司 | Method and device for determining head arrival time of array acoustic logging |
| CN115163052B (en) * | 2022-06-18 | 2023-07-28 | 杭州丰禾石油科技有限公司 | Parameter measurement method of ultrasonic well diameter and ultrasonic well diameter logging while drilling device |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU765772A1 (en) * | 1978-06-12 | 1980-09-23 | Всесоюзный научно-исследовательский институт нефтепромысловой геофизики | Apparatus for acoustic well-logging |
| CN101320017A (en) * | 2007-06-08 | 2008-12-10 | 武汉中科智创岩土技术有限公司 | Method for detecting head wave of sound wave transmission method |
| CN101793955A (en) * | 2010-03-30 | 2010-08-04 | 山东科技大学 | Method and device for detecting and locating head wave of ultrasonic wave |
| CN102094637A (en) * | 2009-12-14 | 2011-06-15 | 西安威尔罗根能源科技有限公司 | Method for detecting digital acoustic head waves |
| CN103388473A (en) * | 2012-05-10 | 2013-11-13 | 中国石油集团长城钻探工程有限公司 | Method and system for eliminating direct wave interference in acoustic logging, and acoustic logging instrument |
-
2016
- 2016-12-28 CN CN201611239317.6A patent/CN106761715B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU765772A1 (en) * | 1978-06-12 | 1980-09-23 | Всесоюзный научно-исследовательский институт нефтепромысловой геофизики | Apparatus for acoustic well-logging |
| CN101320017A (en) * | 2007-06-08 | 2008-12-10 | 武汉中科智创岩土技术有限公司 | Method for detecting head wave of sound wave transmission method |
| CN102094637A (en) * | 2009-12-14 | 2011-06-15 | 西安威尔罗根能源科技有限公司 | Method for detecting digital acoustic head waves |
| CN102094637B (en) * | 2009-12-14 | 2014-07-09 | 西安威尔罗根能源科技有限公司 | Method for detecting digital acoustic head waves |
| CN101793955A (en) * | 2010-03-30 | 2010-08-04 | 山东科技大学 | Method and device for detecting and locating head wave of ultrasonic wave |
| CN103388473A (en) * | 2012-05-10 | 2013-11-13 | 中国石油集团长城钻探工程有限公司 | Method and system for eliminating direct wave interference in acoustic logging, and acoustic logging instrument |
Non-Patent Citations (5)
| Title |
|---|
| 一种新的声波测井信号的首波检测方法;法林等;《石油管材与仪器》;20150831;第1卷(第4期);第33-36页 |
| 交叉偶极阵列声波测井仪控制和处理电路的设计与实现;陈磊等;《测控技术》;20100131;第29卷(第1期);第47-50页 |
| 基于FPGA 的短窗_长窗能量比首波检测技术;叶顶胜等;《重庆科技学院学报（自然科学版）》;20120229;第14卷(第1期);第137-138页 |
| 基于FPGA的存储式声波测井仪井下自动增益控制技术;叶顶胜等;《石油工业计算机应用》;20120930;第32-33页 |
| 阵列声波测井信号调理与首波提取技术研究;张嘉伟等;《中国测试技术》;20060731;第32卷(第4期);第4-5,100页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106761715A (en) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106353792B (en) | Method suitable for positioning micro-seismic source of hydraulic fracturing | |
| US11209565B2 (en) | High precision acoustic logging processing for compressional and shear slowness | |
| CN104216008B (en) | Downhole fracturing microseismic event identification method | |
| CN106761715B (en) | LWF memory-type acoustic logging first arrival detection method | |
| US11112519B2 (en) | Automatic slowness-frequency range determination for advanced borehole sonic data processing | |
| CN104570076B (en) | Automatic seismic wave first-arrival picking method based on dichotomy | |
| CN102073067B (en) | Method for increasing automatic pickup efficiency of geological data first-motion waves | |
| CN108442882B (en) | Shale gas large-displacement horizontal well while-drilling geosteering method based on two-dimensional seismic data | |
| CN107678064B (en) | Real-time extraction method for sound wave time difference | |
| CN109188520B (en) | Thin reservoir thickness prediction method and device | |
| US11353616B2 (en) | Methods and systems for processing slowness values from borehole sonic data | |
| CN105093306A (en) | Method for automatic interpretation and thickness calculation of reservoir in geophysical exploration | |
| CN106837313B (en) | LWF memory-type acoustic logging slowness extracting method | |
| CN103777236A (en) | Vibroseis synchronous scanning wave field forward numerical simulation method | |
| CN103995972A (en) | Stratigraphic correlation method | |
| CN105242311A (en) | Method and apparatus for predicting time window of first arrival pickup | |
| WO2018063156A1 (en) | Iterative migration velocity optimization for a vsp survey using semblance | |
| CN110927817A (en) | Method for optimally designing and controlling thin-layer superposed riverway horizontal well track | |
| CN110568493A (en) | Identification method of complex fault block basin hidden fault | |
| CN107121700B (en) | The production method of structural map under a kind of thrust-nappe belt in front of the mountains | |
| CN107679614B (en) | Particle swarm optimization-based real-time sound wave time difference extraction method | |
| CN105445795A (en) | Method and system for analyzing velocity of tilted interface | |
| CN104199107B (en) | Depth prediction approach and system before brill based on the many wave datum of vertical seismic | |
| US11927711B2 (en) | Enhanced-resolution sonic data processing for formation body wave slowness with full offset waveform data | |
| CN104714248A (en) | Formation overlap-denudation line recognition method based on iteration earthquake DNA detection |
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 |