WO1985001586A1 - Noise subtraction filter - Google Patents

Noise subtraction filter Download PDF

Info

Publication number
WO1985001586A1
WO1985001586A1 PCT/US1983/001494 US8301494W WO8501586A1 WO 1985001586 A1 WO1985001586 A1 WO 1985001586A1 US 8301494 W US8301494 W US 8301494W WO 8501586 A1 WO8501586 A1 WO 8501586A1
Authority
WO
WIPO (PCT)
Prior art keywords
pump
current
average
signature
signals
Prior art date
Application number
PCT/US1983/001494
Other languages
French (fr)
Inventor
T. Umeda
Original Assignee
Exploration Logging, Inc.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Exploration Logging, Inc. filed Critical Exploration Logging, Inc.
Priority to PCT/US1983/001494 priority Critical patent/WO1985001586A1/en
Priority to AU21246/83A priority patent/AU555850B2/en
Priority to EP19830903340 priority patent/EP0159314A4/en
Publication of WO1985001586A1 publication Critical patent/WO1985001586A1/en
Priority to NO852096A priority patent/NO852096L/en

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/12Means 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/14Means 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
    • E21B47/18Means 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 through the well fluid, e.g. mud pressure pulse telemetry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H21/00Adaptive networks
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/30Noise handling
    • G01V2210/32Noise reduction
    • G01V2210/324Filtering

Definitions

  • This invention relates to methods and apparatus for filtering noise signals that occur on a periodic basis.
  • the invention is particularly suitable for use in filtering pump noise in telemetry systems in which data is transmitted in the form of pressure pulses through the fluid that is employed in drilling a well.
  • the oscillating component of the pump pressure signal that is generated by the pump or pumps is often much larger than the pulses that are telemetered through the well drilling fluid, and hence it is essential that the pump noise be removed so that the telemetered pulses can be recognized and detected.
  • the prior art arrangements for filtering pump noise include frequency responsive electric networks, and subtraction circuitry in which pump noise is largely cancelled out by subtracting a preceding pump noise signal from a current pump noise signal.
  • Subtraction circuitry in which a preceding pump noise signal is subtracted from a current pump noise signal has been limited to the use of a single pump because the noise signals of plural pumps interact and beat with one another to produce combined signals which change with time.
  • the aforesaid difficulties are overcome in the present invention by continuously developing and updating an "average pump signature" over a certain number of pump cycles so that the average pump signature is representative of the pressure signals that are produced by the pump during a predetermined interval of time, e.g. one pump cycle.
  • the time scale of the average pump signature is continuously compared witn the time scale of the current pump noise signal, and if necessary, the average pump signature is compressed or expanded in time by linear interpolation to correspond with the time scal of the current pump noise signal.
  • the resultant average pump signature is then subtracted from the current pressure noise signal to leave a "residual" signal which contains the telemetered data pulses.
  • the updating of the average pump signature combines a portion of the last average pump signature witn a portion of the last residual signal.
  • the data pulses have little effect upon the average pump signature since that signature is obtained by averaging over several pump cycles and hence the probability of several data pulses occurring in synchronism with the pump pressure signals is very low.
  • This method and apparatus functions satisfactorily when plural pumps (such as three) are used simultaneously because the pump noise signals are represented by a sum of the average pump signatures.
  • the invention may be carried out by programming a computer or by circuit hardware.
  • Fig. 1 illustrates one form of a pulse code for producing pressure pulses that may be employed in telemetering data from the interior of a well to the surface;
  • Fig. 2 illustrates how encoded pressure pulses may be processed at the surface to provide a read-out of the parameters that are measured downnole;
  • Fig. 3 illustrates four of the signals that are produced by the apparatus of this invention
  • Fig. 4 illustrates the sequence of operations during data acquisition while operating with one pump
  • Fig. 5 illustrates the sequence of operations during the subtraction procedure while operating with one pump
  • Fig. 6 illustrates the sequence of operations during data acquisition while operating with multiple pumps
  • Fig. 7 illustrates the sequence of operations during the subtraction procedure while operating with multiple pumps
  • Fig. 8 shows circuit nardware that may be employed to carry out the steps of Figs. 4 and 5 with a single pump;
  • Fig. 9 snows circuit hardware that may be employed to carry out the steps of Figs. 6 and 7 with multiple pumps.
  • Copending application Serial No. 405,063 discloses and claims methods and apparatus for encoding and synchronizing data in a pulse telemetry system.
  • Fig. 1 illustrates one of the encoding arrangements that is disclosed in that application.
  • a pulse is considered to be a drop in pressure for a fixed period of time (e.g. 0.5 second) followed by one second of normal pressure.
  • Synchronization pulses are produced at the beginning of each synchronization interval.
  • the time intervals between the successive data pulses (D1, D2, D3 ) in a sync interval provides an analog measurement of the data parameters that are to be telemetered to the surface, and three data words may be produced during each synchronization interval.
  • Redundant pulses (R1, R2) are provided midway in time betweeen adjacent pairs of data pulses.
  • the redundant pulses provide a means for discriminating between data pulses and noise pulses.
  • data pulses can be considered to be data pulses (D1, D2, D3) only if a redundant pulse exists midway in time between adjacent pairs of data pulses.
  • Each sync interval is a sub-frame of information, and a selected number of sub-frames (e.g. 11) comprise a frame which constitutes one complete set of data.
  • Fig. 2 illustrates one arrangement for processing the encoded data.
  • the drilling rig includes the usual rotary table 10, kelly 12, swivel 14, traveling block 16, mud pumps 18, mud pit 20, and a drill string made up of drill pipe sections 22 secured to the lower end of the kelly 12 and to the upper end of a drill collar 24 and terminating in the drill bit 26.
  • the downhole pulse encoding apparatus and a valve for producing negative pressure pulses in the drilling fluid may be located in a drill collar 28 located above the drill bit 26.
  • a pressure transducer 30 is coupled to the conduit for the drilling fluid and it senses the negative pressure pulses that are produced downhole.
  • the mud pumps 18 produce pressure signals that have certain characteristics which constitute noise with respect to the data pulses.
  • a strobe generator 32 produces strobe signals that represent the time required for each cycle of each mud pump.
  • the signals from tne pressure transducer 30 and the strobe generator 32 are applied to an input module 34 where the pressure signals are passed througn a low-pass filter and then are averaged in pairs and the average is stored 50 or 100 times per second in a ouffer until required by the filter module 36.
  • the filter module 36 is the subject of this invention. It serves to remove or reduce the noise signals produced by the mud pumps 18.
  • the output of the filter 36 is applied to a detector 38 which employs a matched fiiter to enhance the detectability of the encoded pulses.
  • the pulse signature used in the matched filter is a first order approximation to a rectangular pulse which has been high-pass filtered at a frequency corresponding to 1/(4*pulse width) hz.
  • the use of a high-pass filter allows a simple level detector to be used for pulse identification.
  • a high speed FFT convolution technique is used to implement the matched filter convolution.
  • the convolution is performed approximately once every 20 seconds whenever data is made available by the filter module 36.
  • the output of the matched filter is stored in a "Detection Buffer" which can hold data for more than one complete sub-frame. (If the sub-frame is 60 seconds long, tne buffer is set up to hold data for a 64-second interval.)
  • a simple level detection scheme is performed on the (overlapping) 64-second intervals.
  • the detection level is a user-controlled parameter.
  • Pulse position (time) is assigned at the pulse minimum; pulse height and width are also determined.
  • the decoded pulses are stored in a "Pulse Table" for use by the decoder module 40. If pulse width does not fit witnin certain predefined limits, then the pulse is considered to be "noise” and discarded.
  • Sub-frame sync pulses are found by a search of invariant pulses of one-minute periods over the last three minutes. When these pulses are found, the detection buffer is shifted in time to ensure that the start of the sub-frame will be located at the start of the buffer. The time of the sync pulse is recorded for use by the decoder module 40.
  • a search for frame ident pulses is enabled within the detector module 38 before the sub-frame sync pulses are identified. This enables the apparatus to synchronize on the first sub-frame that is transmitted from downhole.
  • the output of the detector 38 is applied to a decoder 4 ⁇ which identifies tne pulses and calculates the values.
  • This module operates once every sub-frame. Every time the detection buffer is filled, the decoder module 40 uses the information within the Pulse Table set-up by the detector module 38 to identify pulse types.
  • Sync pulses are first identified by their position in time. Frame ident pulses are then searched for, again at predefined time positions. All pulses are allowed a small error in time setup as a user setpoint called "Pulse Position Variance" (typically 0.1 second). All remaining pulses are searcned to identify the "groups of three equally spaced pulses" which identify data and redundant pulses.
  • the output buffer 42 transmits tne data within the decoder 40 to a computer 46 or a printer 44 for setting forth the data parameters that are transmitted to the surface.
  • Fig. 3A illustrates the pump strobe signals when one pump is operating. Each strobe signal is produced at the time of one pump cycle.
  • Fig. 3B illustrates the pump pressure signals when one pump is operating. These signals are much larger than the data pulses, and hence the pump noise signals must be removed or reduced in order to recognize and detect the data pulses.
  • Fig. 3C illustrates an "average pump signature" 50 that is obtained over a certain number of pump cycles. It is shown inverted with respect to tne pressure signals of Fig. 3B to illustrate the subtraction process.
  • the currently updated average pump signature of Fig. 3C is subtracted from eacn cycle of the pump pressure signals of Fig. 3B to expose a residual signal constituting the data signals, as illustrated in Fig. 3D.
  • the acquisition process will typically operate on "interrupts" generated by an external clock, with the subtraction process being the "main line” program.
  • Fig. 4 shows pressure data sampled at 100 times per second being stored into a n input buffer.
  • the current pressure data is transferred periodically to the current pressure buffer. Every time the leading edge of a Pump Strobe signal is recognized an entry is made in the Strobe Table.
  • the pump noise subtraction process of Fig. 5 starts every time a new Pump Strobe entry is made in the Pump Strobe Table. Since the pump speed can vary with time, the period of one cycle will change and therefore the total number of samples in one cycle will change.
  • the current pump signature 50 may require to be compressed or expanded to match the latest pump period.
  • the pump signature 52 illustrates the circumstance where the current pump signature is expanded so that the time base 54 is equal to the time of the current pump cycle. After this process there will be the same number of samples in the Pump Signature as there are in the current pump cycle.
  • the output of the pump noise subtraction filter is the "residual" 56 which is the difference between the pump pressure signals of the latest pump cycle and the current average pump signature.
  • the Pump Signature is updated by adding a fraction of the residual (known as the Weighting factor) to a fraction of the current signature. (Typically i/15 and 14/15 respectively.) This produces a new average pump signature 58. As the portions of the fractions are cnanged the filter can be made to distribute its "weight" over many cycles (and reduce the effect of any noise, including data pulses, on the signature) or fewer cycles, to make the signature more responsive to changes in pump signature.
  • the Strobe Table and Pressure Buffer are next reset to make space for new data.
  • the average Pump Signature is assumed to be zero for all points in the cycle. Normal operation of the filter will cause the signature to gradually take on the characteristics of the actual pump signature. For a weighting factor of W (0 ⁇ W ⁇ 1) the Calculated Pulse Signature will be the following if the last
  • Figs. 6 and 7 show the changes required to operate the model using multiple pumps (simultaneously).
  • Fig. 6 shows that the Pump Strobe Table must now identify the pump number as well as the time of arrival of the Pump Strobe signal.
  • Fig. 7 shows that the subtraction process must wait until a pump strobe signal has been obtained from each operating mud pump.
  • the "predicted" pressure is the sum of all the signatures of the operating pumps, taking into account the relative pnase position of each pump.
  • the residual pressure is now the difference between the current pressure data and the predicted pressure. As before, the residual pressure is the desired output from the filter.
  • the signatures are updated as before except that only i/(total pumps operating) of the residual is assigned to each signature.
  • the Strobe Table and Pressure Buffer are reset to end the process.
  • Fig. 8 illustrates how the input module 34 and filter 36 for removing pump noise can be implemented for one pump with circuit hardware.
  • Signals from the pressure transducer 30 are applied through a bandpass filter 60 to an analog-to-digital converter 62 where the pressure signals are sampled at 100 cycles per second.
  • the pressure signals are averaged in pairs and the average is stored 50 times a second in a current pressure buffer 64.
  • the current pressure buffer 64 accepts data from the analog-to-digital converter 62 during the interval of time between successive pump strobe signals (Fig. 3A).
  • a certain number of memory locations for data words in the current pressure buffer 64 are filled with the digital data when the pump strobe signals are one second apart.
  • the pump speed can vary with time, and tne period of one pump cycle will change. Therefore the number of samples in one pump cycle will change.
  • the current pressure buffer 64 will accept up to five seconds of digital data, so as to encompass all practical pump speeds.
  • the pump strobe signals of Fig. 3A are produced by the strobe generator 32 and applied to a strobe table register 66 whicn enters the successive strobe times.
  • the register 66 is coupled to a microprocessor 67 which controls the sequence of operations.
  • the microprocessor 67 controls the transfer of data from the buffer 64 and the transfer of signatures from buffers 74 and 76.
  • the data in the current pressure buffer 64 and the current average pump signature are shifted to the subtractor circuit 68 wnich subtracts one from tne other to produce at its output 70 the residual signal of Fig. 3D.
  • the residual signal is stored in a residual buffer 72 until it is called for by the detector 38.
  • the average pump signature is obtained by averaging the pressure signals that are produced over a selected number of the current pump cycles in an update pump signa ture averaging buffer 74.
  • the average pump signature is also updated periodically by the addition of a fraction of the current residual signal in the residual buffer 72 to a fraction of the current average pump signature.
  • the fraction of the current residual signal that is added to the current average pump signature is called the "weighting factor".
  • the proportions of the fractions can be cnanged to acnieve particular results. For example, a large number of pump cycles may be employed in developing a signature wherein the effect of any noise, including mud pulses, is reduced. If fewer pump cycles are employed in developing the signature, the average pump signature is more responsive to changes in the pump signature.
  • Typical operating parameters for producing an updated average pump signature employ 14/15 of the old average pump signature plus 1/15 of the current residual signal.
  • the averaging of the pressure signals that are produced during a predetermined number of pump cycles and the addition of a fractional part (e.g. 14/15) to a fractional part (e.g. 1/15) of the residual signal is performed by the update pump signature averaging buffer 74.
  • the current residual signal .from the residual buffer 72 is shifted into the averaging buffer during each pump cycle where the 14/15 and 1/15 addition is effected.
  • the resultant current average pump signature is snifted into the average pump signature buffer 76.
  • a compressor/expander 78 serves to compress or expand the time scale of the current average pump signature in buffer 76 so that it corresponds to tne time scale of the current pressure in buffer 64. Subtraction occurs on a data word by data word basis.
  • the filter operation commences the average pump signature is zero for all points in the cycle.
  • tney are shifted through the subtractor 68 and the residual buffer 72 into tne buffer 74 wnere tne pump signature gradually takes on the characteristics of the actual average pump signature.
  • the average pump signatur is continually updated by the Duffer 74 and the current average pump signature is stored in the buffer 76.
  • the proper data points in the buffer 76 are selected by the compressor/expander 73 and shifted into the subtractor 68 where the current average pump signature and the current pump pressure are subtracted from one another to produce the residual signal that is shifted into the buffer 72.
  • Fig. 9 illustrates how the filter 36 for removing pump noise can be implemented for two or more pumps with circuit hardware.
  • the current pressure data input circuitry is the same as that shown on Fig. 8.
  • An average pump signature is obtained for eacn pump, with the circuitry for the respective pumps being identified by the letters A and B.
  • Tne outputs of the compressor/expanders 78A and 78B are applied to an adder 80 which produces the sura of all of the signatures ot the operating pumps, taking into account the relative phase position of each pump.
  • the signatures are updated as before except that only 1/(number of pumps operating) of the residual is assigned to each signature.
  • the residual pressure is then the difference between the current pressure data in buffer 64 and the predicted pressure signature in adder 80.
  • noise subtraction filter 36 may be expanded to function with any desired number of pumps.
  • P.in(i) Input pressure data at the i-th sample.
  • Is(j,k) Sample index for k-th strobe from pump j
  • P.sig(j,m) Pump j's AC pressure profile at its m-th stroke phase position. 0 ⁇ m ⁇ Signature size - 1.
  • the signature size is 123.
  • P.mean(i) Predicted mean pressure for the i-th sample.
  • P.ave(j) Pump j's average signature value.
  • the filter section receives pressure data, P.in(i), and pump strobe indices, Is (j,k), from the input section.
  • Strobe data, S(j,i) has been evaluated for strobe sample index (strobe occurrence time) by the Input Section.
  • the kth index has been assigned for pump j at data sample i
  • the strode height is a user controlled parameter.
  • This particular algorithm searcnes for the leading edge of a strooe, and it is somewhat noise tolerant due to the two tescs described above.
  • the q-th pump's contribution to the pressure for the i-th sample is referenced in its pump signature table, P.sig, by linearly interpolating the index, m(q), for the sample index bracketed by the indices of two nearest strobes. That is,
  • m(q) Signature size * (i-Is(q,kq-l))/(Is(q,kq)-Is (q,kq-1))
  • the predicted pressure for the i-th sample is computed by summing the individual pump pressure contributions of each pump that is "in operation" to the mean pressure, P.mean(i),
  • the j-th pump is considered "in operation” if its last strobe occurred less than five seconds from the current input sample, i.e., if at the i-tn sample and for last strobe index, Is(j,k),
  • the difference between the predicted pressure, P.pred, and the pressure input, P. in, composes the residual pressure.
  • P.ave(q) is the mean value of the pump signature for q-th pump.
  • this average signature value is subtracted from the elements of the pump's signature table.
  • signature update weight is incremented until its maximum is reached
  • the filter algorithm Whenever there are no pumps in operation and for five seconds after this con ⁇ itions ceases to exist, the filter algorithm generated random pressure data for the residual pressure so that sample timing can be maintained,

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • Geophysics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Acoustics & Sound (AREA)
  • Mining & Mineral Resources (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Fluid Mechanics (AREA)
  • Geochemistry & Mineralogy (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Positive-Displacement Pumps (AREA)
  • Transplanting Machines (AREA)
  • Filters And Equalizers (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Inorganic Insulating Materials (AREA)

Abstract

A filter (36) for use in filtering noise from cyclical signals have a data component and a noise component in which a predetermined number of the cyclical signals are averaged to produce an average signature signal (76). The average signature signal (76) is updated each cycle (74) to produce a current updated average signature signal, and the current cyclical signals and the current updated average signals are substracted (68) from one another to remove the noise component and produce a residual signal (72) which contains the data component. The time frames of the current cyclical signals and the current updated average signals are adjusted (78) to be equal so as to enable point by point subtraction. As an additional feature, the updating of the average signature combines a portion of the last average pump signature with a portion of the last residual signal (74).

Description

NOISE SUBTRACTION FILTER
Field of the Invention
This invention relates to methods and apparatus for filtering noise signals that occur on a periodic basis. The invention is particularly suitable for use in filtering pump noise in telemetry systems in which data is transmitted in the form of pressure pulses through the fluid that is employed in drilling a well. In well drilling the oscillating component of the pump pressure signal that is generated by the pump or pumps is often much larger than the pulses that are telemetered through the well drilling fluid, and hence it is essential that the pump noise be removed so that the telemetered pulses can be recognized and detected.
The Prior Art
The prior art arrangements for filtering pump noise include frequency responsive electric networks, and subtraction circuitry in which pump noise is largely cancelled out by subtracting a preceding pump noise signal from a current pump noise signal.
In order for a frequency responsive electric network to function satisfactorily the frequencies of the pump noise signals must be significantly different from the frequencies of the telemetered pulse signals. Hence, such filters are satisfactory only when the frequency characteristics of the telemetered data pulse signals are different from the frequency characteristics of the pump noise signals.
Subtraction circuitry in which a preceding pump noise signal is subtracted from a current pump noise signal has been limited to the use of a single pump because the noise signals of plural pumps interact and beat with one another to produce combined signals which change with time.
Summary of the Invention
The aforesaid difficulties are overcome in the present invention by continuously developing and updating an "average pump signature" over a certain number of pump cycles so that the average pump signature is representative of the pressure signals that are produced by the pump during a predetermined interval of time, e.g. one pump cycle. The time scale of the average pump signature is continuously compared witn the time scale of the current pump noise signal, and if necessary, the average pump signature is compressed or expanded in time by linear interpolation to correspond with the time scal of the current pump noise signal. The resultant average pump signature is then subtracted from the current pressure noise signal to leave a "residual" signal which contains the telemetered data pulses.
The updating of the average pump signature combines a portion of the last average pump signature witn a portion of the last residual signal. The data pulses have little effect upon the average pump signature since that signature is obtained by averaging over several pump cycles and hence the probability of several data pulses occurring in synchronism with the pump pressure signals is very low. This method and apparatus functions satisfactorily when plural pumps (such as three) are used simultaneously because the pump noise signals are represented by a sum of the average pump signatures.
The invention may be carried out by programming a computer or by circuit hardware.
Brief Description of the Drawings
Fig. 1 illustrates one form of a pulse code for producing pressure pulses that may be employed in telemetering data from the interior of a well to the surface;
Fig. 2 illustrates how encoded pressure pulses may be processed at the surface to provide a read-out of the parameters that are measured downnole;
Fig. 3 illustrates four of the signals that are produced by the apparatus of this invention;
Fig. 4 illustrates the sequence of operations during data acquisition while operating with one pump;
Fig. 5 illustrates the sequence of operations during the subtraction procedure while operating with one pump;
Fig. 6 illustrates the sequence of operations during data acquisition while operating with multiple pumps;
Fig. 7 illustrates the sequence of operations during the subtraction procedure while operating with multiple pumps;
Fig. 8 shows circuit nardware that may be employed to carry out the steps of Figs. 4 and 5 with a single pump; and
Fig. 9 snows circuit hardware that may be employed to carry out the steps of Figs. 6 and 7 with multiple pumps.
Description of the Preferred Embodiments
Pulse Encoding
Copending application Serial No. 405,063 discloses and claims methods and apparatus for encoding and synchronizing data in a pulse telemetry system.
Fig. 1 illustrates one of the encoding arrangements that is disclosed in that application. A pulse is considered to be a drop in pressure for a fixed period of time (e.g. 0.5 second) followed by one second of normal pressure. Synchronization pulses are produced at the beginning of each synchronization interval. The time intervals between the successive data pulses (D1, D2, D3 ) in a sync interval provides an analog measurement of the data parameters that are to be telemetered to the surface, and three data words may be produced during each synchronization interval.
Redundant pulses (R1, R2) are provided midway in time betweeen adjacent pairs of data pulses. The redundant pulses provide a means for discriminating between data pulses and noise pulses. Thus, data pulses can be considered to be data pulses (D1, D2, D3) only if a redundant pulse exists midway in time between adjacent pairs of data pulses.
Each sync interval is a sub-frame of information, and a selected number of sub-frames (e.g. 11) comprise a frame which constitutes one complete set of data.
During certain sub-frames a particular pattern of frame identification pulses is produced and telemetered to the surface to enable the apparatus at the surface to identify the sequence of data transmission. Processing the Encoded Data
Fig. 2 illustrates one arrangement for processing the encoded data.
The drilling rig includes the usual rotary table 10, kelly 12, swivel 14, traveling block 16, mud pumps 18, mud pit 20, and a drill string made up of drill pipe sections 22 secured to the lower end of the kelly 12 and to the upper end of a drill collar 24 and terminating in the drill bit 26. The downhole pulse encoding apparatus and a valve for producing negative pressure pulses in the drilling fluid may be located in a drill collar 28 located above the drill bit 26.
A pressure transducer 30 is coupled to the conduit for the drilling fluid and it senses the negative pressure pulses that are produced downhole.
The mud pumps 18 produce pressure signals that have certain characteristics which constitute noise with respect to the data pulses.
A strobe generator 32 produces strobe signals that represent the time required for each cycle of each mud pump.
The signals from tne pressure transducer 30 and the strobe generator 32 are applied to an input module 34 where the pressure signals are passed througn a low-pass filter and then are averaged in pairs and the average is stored 50 or 100 times per second in a ouffer until required by the filter module 36.
The filter module 36 is the subject of this invention. It serves to remove or reduce the noise signals produced by the mud pumps 18.
The output of the filter 36 is applied to a detector 38 which employs a matched fiiter to enhance the detectability of the encoded pulses. The pulse signature used in the matched filter is a first order approximation to a rectangular pulse which has been high-pass filtered at a frequency corresponding to 1/(4*pulse width) hz. The use of a high-pass filter allows a simple level detector to be used for pulse identification.
"Ideal" mud pulse:
Actual matched filter pulse signature:
Figure imgf000009_0001
Output from matched filter:
Figure imgf000009_0002
The lowest sharp peak corresponds to the best match between actual psipulse and filter signature. seconds
A high speed FFT convolution technique is used to implement the matched filter convolution. In operation, the convolution is performed approximately once every 20 seconds whenever data is made available by the filter module 36. The output of the matched filter is stored in a "Detection Buffer" which can hold data for more than one complete sub-frame. (If the sub-frame is 60 seconds long, tne buffer is set up to hold data for a 64-second interval.)
A simple level detection scheme is performed on the (overlapping) 64-second intervals. (The detection level is a user-controlled parameter.) Pulse position (time) is assigned at the pulse minimum; pulse height and width are also determined. The decoded pulses are stored in a "Pulse Table" for use by the decoder module 40. If pulse width does not fit witnin certain predefined limits, then the pulse is considered to be "noise" and discarded.
Sub-frame sync pulses are found by a search of invariant pulses of one-minute periods over the last three minutes. When these pulses are found, the detection buffer is shifted in time to ensure that the start of the sub-frame will be located at the start of the buffer. The time of the sync pulse is recorded for use by the decoder module 40.
The first time a downhole Tool starts up downhole, a search for frame ident pulses is enabled within the detector module 38 before the sub-frame sync pulses are identified. This enables the apparatus to synchronize on the first sub-frame that is transmitted from downhole.
The output of the detector 38 is applied to a decoder 4ϋ which identifies tne pulses and calculates the values.
This module operates once every sub-frame. Every time the detection buffer is filled, the decoder module 40 uses the information within the Pulse Table set-up by the detector module 38 to identify pulse types.
Sync pulses are first identified by their position in time. Frame ident pulses are then searched for, again at predefined time positions. All pulses are allowed a small error in time setup as a user setpoint called "Pulse Position Variance" (typically 0.1 second). All remaining pulses are searcned to identify the "groups of three equally spaced pulses" which identify data and redundant pulses.
Once the pulses are identified sub-frame numbers and data values can be calculated and supplied to the output buffer 42 via the Decoder Buffer. The output buffer 42 transmits tne data within the decoder 40 to a computer 46 or a printer 44 for setting forth the data parameters that are transmitted to the surface.
The Pump Noise Subtraction Filter
Fig. 3A illustrates the pump strobe signals when one pump is operating. Each strobe signal is produced at the time of one pump cycle.
Fig. 3B illustrates the pump pressure signals when one pump is operating. These signals are much larger than the data pulses, and hence the pump noise signals must be removed or reduced in order to recognize and detect the data pulses.
Fig. 3C illustrates an "average pump signature" 50 that is obtained over a certain number of pump cycles. It is shown inverted with respect to tne pressure signals of Fig. 3B to illustrate the subtraction process.
The currently updated average pump signature of Fig. 3C is subtracted from eacn cycle of the pump pressure signals of Fig. 3B to expose a residual signal constituting the data signals, as illustrated in Fig. 3D.
There are two processes running simultaneously during tne pump noise filter operation. The first is tne data acquisition, shown m Figure 4 (this is the i.nput module 34 of Figure 2); the second is the subtraction process, snown in Figure 5 (this is trie filter module 36 of Figure 2).
In a computer base-d process the acquisition process will typically operate on "interrupts" generated by an external clock, with the subtraction process being the "main line" program.
Fig. 4 shows pressure data sampled at 100 times per second being stored into a n input buffer. The current pressure data is transferred periodically to the current pressure buffer. Every time the leading edge of a Pump Strobe signal is recognized an entry is made in the Strobe Table.
The pump noise subtraction process of Fig. 5 starts every time a new Pump Strobe entry is made in the Pump Strobe Table. Since the pump speed can vary with time, the period of one cycle will change and therefore the total number of samples in one cycle will change. The current pump signature 50 may require to be compressed or expanded to match the latest pump period. The pump signature 52 illustrates the circumstance where the current pump signature is expanded so that the time base 54 is equal to the time of the current pump cycle. After this process there will be the same number of samples in the Pump Signature as there are in the current pump cycle.
The output of the pump noise subtraction filter is the "residual" 56 which is the difference between the pump pressure signals of the latest pump cycle and the current average pump signature.
The Pump Signature is updated by adding a fraction of the residual (known as the Weighting factor) to a fraction of the current signature. (Typically i/15 and 14/15 respectively.) This produces a new average pump signature 58. As the portions of the fractions are cnanged the filter can be made to distribute its "weight" over many cycles (and reduce the effect of any noise, including data pulses, on the signature) or fewer cycles, to make the signature more responsive to changes in pump signature.
The Strobe Table and Pressure Buffer are next reset to make space for new data. When the filter is started (from scratch) the average Pump Signature is assumed to be zero for all points in the cycle. Normal operation of the filter will cause the signature to gradually take on the characteristics of the actual pump signature. For a weighting factor of W (0 < W < 1) the Calculated Pulse Signature will be the following if the last
Dump cycle was cycle "i" :
Signature=W(Residual of Pump Cycle(i))+W2(Residual of Cycle(i-l)) +W3(Residual of Cycle (i-2))+...
Figs. 6 and 7 show the changes required to operate the model using multiple pumps (simultaneously).
Fig. 6 shows that the Pump Strobe Table must now identify the pump number as well as the time of arrival of the Pump Strobe signal.
Fig. 7 shows that the subtraction process must wait until a pump strobe signal has been obtained from each operating mud pump. The "predicted" pressure is the sum of all the signatures of the operating pumps, taking into account the relative pnase position of each pump.
The residual pressure is now the difference between the current pressure data and the predicted pressure. As before, the residual pressure is the desired output from the filter.
The signatures are updated as before except that only i/(total pumps operating) of the residual is assigned to each signature.
The Strobe Table and Pressure Buffer are reset to end the process.
The algorithm for tne pump noise subtraction filter described above is set forth in detail in Appendix I hereto. Fig. 8 illustrates how the input module 34 and filter 36 for removing pump noise can be implemented for one pump with circuit hardware.
Signals from the pressure transducer 30 are applied through a bandpass filter 60 to an analog-to-digital converter 62 where the pressure signals are sampled at 100 cycles per second. The pressure signals are averaged in pairs and the average is stored 50 times a second in a current pressure buffer 64. The current pressure buffer 64 accepts data from the analog-to-digital converter 62 during the interval of time between successive pump strobe signals (Fig. 3A).
A certain number of memory locations for data words in the current pressure buffer 64 are filled with the digital data when the pump strobe signals are one second apart. However, the pump speed can vary with time, and tne period of one pump cycle will change. Therefore the number of samples in one pump cycle will change. The current pressure buffer 64 will accept up to five seconds of digital data, so as to encompass all practical pump speeds.
The pump strobe signals of Fig. 3A are produced by the strobe generator 32 and applied to a strobe table register 66 whicn enters the successive strobe times. The register 66 is coupled to a microprocessor 67 which controls the sequence of operations. The microprocessor 67 controls the transfer of data from the buffer 64 and the transfer of signatures from buffers 74 and 76.
One subtraction process is initiated each time that a new entry is made in tne strobe table register 66.
Upon the occurrence of a new entry in the strobe table register 66, the data in the current pressure buffer 64 and the current average pump signature are shifted to the subtractor circuit 68 wnich subtracts one from tne other to produce at its output 70 the residual signal of Fig. 3D. The residual signal is stored in a residual buffer 72 until it is called for by the detector 38.
The average pump signature is obtained by averaging the pressure signals that are produced over a selected number of the current pump cycles in an update pump signa ture averaging buffer 74.
Preferably, the average pump signature is also updated periodically by the addition of a fraction of the current residual signal in the residual buffer 72 to a fraction of the current average pump signature. The fraction of the current residual signal that is added to the current average pump signature is called the "weighting factor". The proportions of the fractions can be cnanged to acnieve particular results. For example, a large number of pump cycles may be employed in developing a signature wherein the effect of any noise, including mud pulses, is reduced. If fewer pump cycles are employed in developing the signature, the average pump signature is more responsive to changes in the pump signature.
Typical operating parameters for producing an updated average pump signature employ 14/15 of the old average pump signature plus 1/15 of the current residual signal.
The averaging of the pressure signals that are produced during a predetermined number of pump cycles and the addition of a fractional part (e.g. 14/15) to a fractional part (e.g. 1/15) of the residual signal is performed by the update pump signature averaging buffer 74. The current residual signal .from the residual buffer 72 is shifted into the averaging buffer during each pump cycle where the 14/15 and 1/15 addition is effected. The resultant current average pump signature is snifted into the average pump signature buffer 76. A compressor/expander 78 serves to compress or expand the time scale of the current average pump signature in buffer 76 so that it corresponds to tne time scale of the current pressure in buffer 64. Subtraction occurs on a data word by data word basis.
Thus the appropriate word from the pump signature buffer is subtracted from the current word in the data buffer. The word selected from the average pump signature buffer 76 is selected by the compressor/expander 78 under the control of microprocessor 67 according to the following formula: m = s
Figure imgf000016_0001
where m is the index of the word selected from the pump signature buffer, s is the total number of words in the pump signature buffer, i is the index of the current word m the data buffer 64 and n is the length, in words of the pressure data stored in tne data buffer 64.
When the filter operation commences the average pump signature is zero for all points in the cycle. As the initial current pressure cycles are stored in the buffer 64, tney are shifted through the subtractor 68 and the residual buffer 72 into tne buffer 74 wnere tne pump signature gradually takes on the characteristics of the actual average pump signature. The average pump signatur is continually updated by the Duffer 74 and the current average pump signature is stored in the buffer 76. The proper data points in the buffer 76 are selected by the compressor/expander 73 and shifted into the subtractor 68 where the current average pump signature and the current pump pressure are subtracted from one another to produce the residual signal that is shifted into the buffer 72. Fig. 9 illustrates how the filter 36 for removing pump noise can be implemented for two or more pumps with circuit hardware.
The current pressure data input circuitry is the same as that shown on Fig. 8.
An average pump signature is obtained for eacn pump, with the circuitry for the respective pumps being identified by the letters A and B.
Tne outputs of the compressor/expanders 78A and 78B are applied to an adder 80 which produces the sura of all of the signatures ot the operating pumps, taking into account the relative phase position of each pump.
The signatures are updated as before except that only 1/(number of pumps operating) of the residual is assigned to each signature.
The residual pressure is then the difference between the current pressure data in buffer 64 and the predicted pressure signature in adder 80.
It will be apparent that the noise subtraction filter 36 may be expanded to function with any desired number of pumps.
APPENDIX 1
Detailed Description of the Filter Section Algorithm
The filtering algorithm employed in this section is the pump phase partition model, and the assumptions upon which this model is based have been stated in the description of the pump noise subtraction filter. The implementation details are described in the following: SYMBOL DEFINITIONS
P.in(i) = Input pressure data at the i-th sample.
S(j,i) = Input strobe data at i-th sample from pump j (l<j<3)
Is(j,k) = Sample index for k-th strobe from pump j
P.pred(i) = Predicted pressure for the i-th sample
P.sig(j,m) = Pump j's AC pressure profile at its m-th stroke phase position. 0 < m < Signature size - 1. The signature size is 123.
P.mean(i) = Predicted mean pressure for the i-th sample.
P.res(i) = Residual pressure at the i-th sample =
P.in(i) - P.pred(i). This is the filter sections output to the detector.
P.ave(j) = Pump j's average signature value.
P.ran(i) = Computer generated random pressure data. = Pump j ' s signature update weight ; 0 <U< Update weight l imit .
The filter section receives pressure data, P.in(i), and pump strobe indices, Is (j,k), from the input section. Strobe data, S(j,i), has been evaluated for strobe sample index (strobe occurrence time) by the Input Section. The kth index has been assigned for pump j at data sample i
Is(j,k) = i whenever
S(j,i)+S(j,i-1)+S(j,i-2)-S(j,i-3)-S(j,1-4)-S(j,i -5) > (3/4) * Strobe height (j) and
Minimum (S(j,i),3(j,i-1),S(j,i-2))
> Maximum (S(j,i-3),S(j,i-4),S(j,i-5))
(The strode height is a user controlled parameter.) This particular algorithm searcnes for the leading edge of a strooe, and it is somewhat noise tolerant due to the two tescs described above.
The q-th pump's contribution to the pressure for the i-th sample is referenced in its pump signature table, P.sig, by linearly interpolating the index, m(q), for the sample index bracketed by the indices of two nearest strobes. That is,
m(q) = Signature size * (i-Is(q,kq-l))/(Is(q,kq)-Is (q,kq-1))
where kq is selected on the basis that
Is(q,kq) > i > Is (q,kq-1) This linear interpolation of the index estimates the phase of the pump.
The predicted pressure for the i-th sample is computed by summing the individual pump pressure contributions of each pump that is "in operation" to the mean pressure, P.mean(i),
P.pred(i)=P.mean (i)-t-[sum q over all pumps in operation:
P.sig(q,m(q))]
The j-th pump is considered "in operation" if its last strobe occurred less than five seconds from the current input sample, i.e., if at the i-tn sample and for last strobe index, Is(j,k),
i-Is(j,k) < 500 (= 5 seconds)
If the pump is considered "off", its signature and update weight are initilized to zero.
The difference between the predicted pressure, P.pred, and the pressure input, P. in, composes the residual pressure.
P.res(i) = P.in(i) - P.pred(i)
which is passed to the detector section.
pump signatures, P.sig(q,m(q) ), are updated at every sample point by the relation
P.sig(q,m(q)) =
(P.sig(q,m(q))*U(q)+P.res(i))/(number of pumps in operation)
Maximum (l,U(q)+1)
After a pump stroke is complete, i.e., i=Is(q,kq), the predicted mean pressure is updated by
P.mean(i)=(P.mean(i)+P.ave(q))/2
where P.ave(q) is the mean value of the pump signature for q-th pump.
After the predicted mean pressure is updated, this average signature value is subtracted from the elements of the pump's signature table.
Also, the signature update weight is incremented until its maximum is reached,
U(q) =Minimum(U(q)+1, Signature Update Weight Limit)
Whenever there are no pumps in operation and for five seconds after this conαitions ceases to exist, the filter algorithm generated random pressure data for the residual pressure so that sample timing can be maintained,
P.res(i)=P.ran(i).

Claims

WHAT IS CLAIMED IS:
I. A method of filtering noise from signals having a data component and a cyclical noise component comprising averaging a predetermined number of the cyclical signals to produce an average signature signal, updating the average signature signal each cycle to produce a current average signature signal, and subtracting the current cyclical signals and the current average signature signals from one another to produce a residual signal which contains the data component.
2. A method ot noise subtraction for use in processing data that is telemetered oy pressure pulses that are transmitted througn well drilling fluid tnat is circulated by pumping in repetitive cycles, comprising producing composite signals that are representative of the pressure changes in the fluid having a data component that is representative of the telemetered data and naving a cyclical noise component that is produced by the pumping of the fluid, producing an average signature signal that is representative of tne average of saic composite signals over a selected number or current pumping cycles, mid subtracting the current composite signals and tne current average signature signal from one another to cause the cyclical noise components that are produced PV pumping to cancel leaving a residual signal that contains the data component.
3. The method of Claim 4 whhrein the time frame of the signals is adjusted so that the composite signals and the average signature signais have the same time frame witn respect to the pumping cycles.
4. The method of Claim 2 wnerein the average signature signal is produced by adding a portion of the current pump signature to a portion of tne current residual signal.
5. The metnod of Claim 2 wnerein the average signature signal is produced by adding 14/15 or the current pump signature to 1/15 of the current residual signal.
6. A noise subtraction filter comprising means for producing a repetitive signal having a data component and a noise component, means for continuously averaging a predetermined number of the current repetitive signals to produce an average signature signal, means for updating the average signature signal periodically to produce a current average signature signal, and means for subtracting the current repetitive signals and tne current average signature signal, whereby the noise component in the current repetitive signal is substantially cancelled leaving a residual signal which contains the data component.
7. A noise subtraction filter for use in processing data that is telemetered through well drilling fluid that is circulated υy pump means that operates in repetitive cycles, comprising means adapted to oe coupled to tne well drilling fluid for producing composite signals representative of pressure changes in the fluid having a data component that is representative of the data that is telemetered through tne well drilling fluid and having a cyclical noise component that results from the changes in pressure that are produced by tne pump means, means for averaging a predetermined number υf said composite signals during a predetermined number of pump cycles to produce an average signature signal, and means for subtracting the current composite signals and the current average signature signal from one anotner to cause the cyclical noise components to cancel leaving a residual signal wnicn contains the data component.
8. A noise subtraction filter for use in processing data that is telemetered through well drilling fluid that is circulated by pump means that operates in repetitive cycles, comprising means coupled to tne well drilling fluid tor producing composite signals representative of pressure changes in the fluid having a data component that is representative of the data tnat is telemetered through the well drilling fluid and having a cyclical noise component that results from the cnanges in pressure that are produced by tne pump means, means for producing an average pump signature signal that is representative of the average of said composite signals during a predetermined number of pump cycles, and means for subtracting the current composite signals and the current average signature signal from one anotner to cause the cyclical noise components tnat are produced by the pump means to cancel leaving a residual signal at the output of tne subtracting means that contains tne data component.
9. The apparatus of claim 3 wherein the means for subtracting tne current composite signals and the current average signature signals includes means for compressing or expanding tne time rrame ot one of the signals to cause both of the signals to have the same time frame.
10. The apparatus of Claim 3 wnerein the average pump signature signal is produced by adding a portion of the current pump signature to a portion of tne current residual signal.
11. The apparatus of Claim 8 wherein the average pump signature signal is produced by adding 14/15 of the current pump signature to 1/15 of the current residual signal.
PCT/US1983/001494 1983-09-26 1983-09-26 Noise subtraction filter WO1985001586A1 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
PCT/US1983/001494 WO1985001586A1 (en) 1983-09-26 1983-09-26 Noise subtraction filter
AU21246/83A AU555850B2 (en) 1983-09-26 1983-09-26 Noise subtraction filter
EP19830903340 EP0159314A4 (en) 1983-09-26 1983-09-26 Noise subtraction filter.
NO852096A NO852096L (en) 1983-09-26 1985-05-24 PROCEDURE AND APPARATUS FOR NOISE FILTERING

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/US1983/001494 WO1985001586A1 (en) 1983-09-26 1983-09-26 Noise subtraction filter

Publications (1)

Publication Number Publication Date
WO1985001586A1 true WO1985001586A1 (en) 1985-04-11

Family

ID=22175458

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US1983/001494 WO1985001586A1 (en) 1983-09-26 1983-09-26 Noise subtraction filter

Country Status (4)

Country Link
EP (1) EP0159314A4 (en)
AU (1) AU555850B2 (en)
NO (1) NO852096L (en)
WO (1) WO1985001586A1 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU633673B2 (en) * 1990-01-18 1993-02-04 Matsushita Electric Industrial Co., Ltd. Signal processing device
EP0991958A4 (en) * 1996-01-16 2000-04-12 Mobil Oil Corp Method for removing borehole rugosity noise from well log data
GB2424707A (en) * 2005-03-31 2006-10-04 Schlumberger Holdings Well logging while a mud pump is off
EP3035083A1 (en) * 2014-12-17 2016-06-22 Services Pétroliers Schlumberger Systems and methods for removing noise from acoustic impedance logs

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3484591A (en) * 1966-07-18 1969-12-16 Hewlett Packard Co Extended bandwidth signal-to-noise ratio enhancement methods and means
US3714590A (en) * 1971-11-08 1973-01-30 Amoco Prod Co Averaging circuit
US4025724A (en) * 1975-08-12 1977-05-24 Westinghouse Electric Corporation Noise cancellation apparatus
US4144578A (en) * 1976-02-12 1979-03-13 Carl Schenck Ag Method and apparatus for monitoring and/or correcting of a physical characteristic
US4193118A (en) * 1978-07-18 1980-03-11 Motorola, Inc. Low pass digital averaging filter
US4357673A (en) * 1980-04-18 1982-11-02 Hewlett-Packard Company Apparatus for performing measurements and error analysis of the measurements

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA1189442A (en) * 1981-11-09 1985-06-25 Gary D. Berkenkamp Pump noise filtering apparatus for a borehole measurement while drilling system utilizing drilling fluid pressure sensing

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3484591A (en) * 1966-07-18 1969-12-16 Hewlett Packard Co Extended bandwidth signal-to-noise ratio enhancement methods and means
US3714590A (en) * 1971-11-08 1973-01-30 Amoco Prod Co Averaging circuit
US4025724A (en) * 1975-08-12 1977-05-24 Westinghouse Electric Corporation Noise cancellation apparatus
US4144578A (en) * 1976-02-12 1979-03-13 Carl Schenck Ag Method and apparatus for monitoring and/or correcting of a physical characteristic
US4193118A (en) * 1978-07-18 1980-03-11 Motorola, Inc. Low pass digital averaging filter
US4357673A (en) * 1980-04-18 1982-11-02 Hewlett-Packard Company Apparatus for performing measurements and error analysis of the measurements

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP0159314A4 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU633673B2 (en) * 1990-01-18 1993-02-04 Matsushita Electric Industrial Co., Ltd. Signal processing device
EP0991958A4 (en) * 1996-01-16 2000-04-12 Mobil Oil Corp Method for removing borehole rugosity noise from well log data
EP0991958A1 (en) * 1996-01-16 2000-04-12 Mobil Oil Corporation Method for removing borehole rugosity noise from well log data
GB2424707A (en) * 2005-03-31 2006-10-04 Schlumberger Holdings Well logging while a mud pump is off
GB2424707B (en) * 2005-03-31 2007-06-06 Schlumberger Holdings Pump-off measurements for quality control and wellbore stability prediction
US7251566B2 (en) 2005-03-31 2007-07-31 Schlumberger Technology Corporation Pump off measurements for quality control and wellbore stability prediction
EP3035083A1 (en) * 2014-12-17 2016-06-22 Services Pétroliers Schlumberger Systems and methods for removing noise from acoustic impedance logs
US10871589B2 (en) 2014-12-17 2020-12-22 Schlumberger Technology Corporation System and methods for removing noise from acoustic impedance logs

Also Published As

Publication number Publication date
EP0159314A1 (en) 1985-10-30
AU555850B2 (en) 1986-10-09
NO852096L (en) 1985-05-24
AU2124683A (en) 1985-04-23
EP0159314A4 (en) 1986-07-23

Similar Documents

Publication Publication Date Title
US4642800A (en) Noise subtraction filter
US5955966A (en) Signal recognition system for wellbore telemetry
US4513403A (en) Data encoding and synchronization for pulse telemetry
US4912764A (en) Digital speech coder with different excitation types
CA2021508C (en) Digital speech coder having improved long term lag parameter determination
US8433533B2 (en) High resolution sensor with scalable sample rate
AU555850B2 (en) Noise subtraction filter
EP0235180B1 (en) Voice synthesis utilizing multi-level filter excitation
JP4336745B2 (en) Lossless compression method for signals with wide dynamic range
EP0078907A2 (en) Pump noise filtering apparatus for a borehole measurement while drilling system utilizing drilling fluid pressure sensing
CN114640695A (en) High-frequency time sequence data effective transmission method based on long sequence bi-prediction and informar for intelligent factory
JPS5886296A (en) Pump noise filter for measuring device among bored well in oil well
EP0545403B1 (en) Speech signal encoding system capable of transmitting a speech signal at a low bit rate
US5202953A (en) Multi-pulse type coding system with correlation calculation by backward-filtering operation for multi-pulse searching
CN1085367A (en) Obtain the method for signal bursts quality numerical value
AU2121083A (en) Data encoding and synchronization for pulse telemetry
GB2206768A (en) Data transmission by concurrent relative phase coding in multiple spectral carriers
US5224061A (en) Method of and arrangement for restoring invalid samples of an equidistantly sampled signal
CA2201541C (en) Signal recognition system for wellbore telemetry
CA1257002A (en) Pole-zero analyzer
EP0341271A1 (en) Coded communications system
SU1054828A1 (en) System for transmitting telemetric data
CA2214582A1 (en) Speech signal encoding system capable of transmitting a speech signal at a low bit rate without carrying out a large volume of calculation
JPS6074744A (en) Synchronizing circuit
JPS6367200B2 (en)

Legal Events

Date Code Title Description
AK Designated states

Designated state(s): AU NO

AL Designated countries for regional patents

Designated state(s): DE FR GB

WWE Wipo information: entry into national phase

Ref document number: 1983903340

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 1983903340

Country of ref document: EP

WWW Wipo information: withdrawn in national office

Ref document number: 1983903340

Country of ref document: EP