Abstract: Method for controlling the phase spectrum of seismic data to match assumptions inherent in subsequent processing steps. The source signature, after processing with the same initial processing steps used on the data, is used to design a phase control filter that shapes the seismic data to have a minimum phase spectrum or whatever other phase spectrum the subsequent processing algorithms may assume. The processed data is then filtered with a second phase-control filter, also designed using the parallel-processed signature, to shape the data to zero phase or whatever other phase may be desired for interpretation of the data.

Abstract: A method enables the time required to complete a seismic survey and the noises recorded in the seismic data to be reduced. One aspect of the method includes actuating a first vibrator group to start a first sweep at time T0; and a second vibrator group to start a second sweep at time T1, wherein T0<T1<T0+S1+L, where S1 is a sweep time of the first vibrator group and L is a listening time; wherein the time between the first sweep and the second sweep is T1?T0?((n?1)*f1*S1)/(n*(f1?f0)) where n is a natural number, f0 is the lower frequency limit of the vibrator sweep and f1 is the upper frequency limit of the vibrator sweep.

Abstract: An arrangement for conducting a seismic survey and methods for operating a vibrator array in a seismic survey. A signal source coupled to an array of vibrators drives the vibrators with a set of signals to generate seismic signals into a survey area. Each of the signals of the set has an autocorrelation function and a cross correlation function with the other signals of the set that decays faster with time than reflections of the seismic signals from the survey area. One such set of signals includes a linear FM sweep phase-modulated by a maximal-length sequence code shifted a different number of chips for each signal.

Abstract: A method and apparatus for minimizing interference between seismic systems. The method may generally include: (a) actuating a plurality of first sources associated with a first seismic system to generate a first plurality of phase-encoded vibratory signals; (b) actuating a plurality of second sources associated with a second seismic system to generate a second plurality of phase-encoded vibratory signals that are at least partially frequency separated from the first plurality of phase-encoded vibratory signals generated in (a); and (c) detecting the first and second plurality of phase-encoded vibratory signals utilizing a first receiver positioned at a location remote from the first and second sources.

Abstract: The force exerted by a reciprocating actuator is estimated as a weighted sum of three quantities: a reaction-mass acceleration, a baseplate acceleration, and a relative acceleration of the reaction-mass relative to the baseplate. In one embodiment, acceleration sensors measure the reaction-mass acceleration and the baseplate acceleration; and a displacement transducer measures a relative displacement of the reaction-mass relative to the baseplate. A double-differentiator generates the value of the relative acceleration from the second derivative of the relative displacement. In the sum, the relative acceleration is weighted by factor representative of the ratio of the difference between the masses of the reaction-mass and the baseplate to the sum of the masses of the reaction-mass and the baseplate.

Abstract: A method for simultaneously operating multiple seismic vibrators using continuous sweeps (little or no “listening” time between sweeps) for each vibrator, and recovering the separated seismic responses for each vibrator with the earth signature removed. Each vibrator is given a unique, continuous pilot signal. The earth response to the motion of each vibrator is measured or estimated. The vibrator motion records for each vibrator and the combined seismic data record for all the vibrators are parsed into separate shorter records. The shorter records are then used to form a system of simultaneous linear equations in the Fourier transform domain, following the HFVS method of Sallas and Allen. The equations are then solved for the separated earth responses.

Abstract: A method of processing seismic data obtained using a seismic vibrator that includes estimating travel times of direct wave arrivals between the seismic vibrator and an array of seismic receiver locations and wavefield deconvolving the seismic data using the direct wave arrival travel times. Also a related method of producing a high-frequency geological subsurface image that includes acquiring seismic data having significant harmonic energy using a seismic vibrator, estimating direct wave arrival travel times between the seismic vibrator and an array of seismic receiver locations, wavefield deconvolving the seismic data using the estimated direct wave arrival travel times, and using the wavefield deconvolved seismic data to produce a high-frequency geological subsurface image. A further related computer useable medium having computer readable program code means embodied therein practicing the inventive method.

Abstract: A method for the simultaneous operation of multiple seismic vibrators using unique modified pseudorandom sweeps and recovery of the transmission path response from each vibrator is disclosed. The vibrator sweeps are derived from pseudorandom binary sequences modified to be weakly correlated over a time window of interest, spectrally shaped and amplitude level compressed. Cross-correlation with each pilot signal is used to perform an initial separation of the composite received signal data set. Recordings of the motion of each vibrator are also cross-correlated with each pilot, windowed, and transformed to form a source cross-spectral density matrix in the frequency domain useful for source signature removal and for additional crosstalk-suppression between the separated records. After source signature removal in the frequency domain an inverse transform is applied to produce an estimate of each source-to-receiver earth response in the time domain.

Abstract: Embodiments of the present invention provide systems and methods for deriving a source signature for an array of seismic sources for marine seismic analysis, wherein the systems and methods include factors for determining the source signature taking into consideration both sea-surface and sea-floor reflections of signals produced by the seismic sources. In certain aspects, reflection coefficients of the sea-surface and the sea-floor and relative lengths of paths between a seismic source image and a detector are applied to a series of simultaneous equations that are solved using measurements of the pressure field produced by the array at a plurality of known locations to determine individual source signatures for each seismic source in the array and these individual source signatures may then be superposed to provide a source signature for the array that accounts for sea-floor reflections of outputs from the seismic sources.

Abstract: The method for exploring desired characteristics of a subsurface sector, having at least one resonant frequency, is based on selectively transmitting suitable narrowband energy waves into the subsurface sector, thereby producing narrowband signals reflected off the subsurface sector. The transmitted narrowband energy waves can be selectively and optimally adjusted in real time so as to provide optimum illumination of the desired characteristics from the explored sector.

Abstract: The present invention provides a seismic surveying method in which vibrators in a first vibrator group are actuated at time T0 and vibrators in a second vibrator group are actuated at time T1, where T0<T1<T0+S1+L and S1 is the sweep time of the first group and L is the listening time. This method enables the time required to complete a seismic survey to be reduced. In other embodiments, vibrators in the first vibrator group may be actuated at time T2, where T1<T2<T1+S2+L and S2 is the sweep time of the first group, and vibrators in the second vibrator group at time T3 where T2<T3<T2+S1+L and T3?T2?T1?T0. Such embodiments providing that by appropriately combining the shot records noise may be eliminated.

Abstract: A method and related apparatus are described for generating acoustic signals for use in a vibratory seismic survey, including the step of combining into a drive signal a high frequency sweep signal, which sweeps upwardly through a high frequency band during a first time interval, and a low frequency sweep signal which is of lower amplitude than the high frequency sweep signal and which sweeps upwardly through a low frequency band during a second time interval, wherein the second time interval starts during the first time interval but after the beginning thereof; and applying the drive signal to a mechanical drive system for a vibratable element. The method improves the utilization of a single vibratory source.

Abstract: In accordance with the present invention, there is provided a method of determining whether a particular high fidelity vibratory seismic survey phase encoding is likely to be a good one based on an analysis of the eigenvalue structure (i.e., eigenvalues, eigenvalue separation, condition number, and model resolution matrix) of a matrix formed from the Fourier transforms of the sweep signals. Preferably, a singular value decomposition will be used to calculate the eigenvalues. Using this same approach, the condition number and eigenvalues of matrices that are associated with multiple proposed designs can be compared with each other to determine which is likely to be yield the best seismic data. This approach is preferably used either as a component of the advanced planning for a survey or in the field during pre-survey testing. The use of the instant approach to determine an optimal phase encoding scheme is also taught.

Abstract: The present invention provides a method for harmonic noise attenuation in correlated sweep data in seismic exploration. The method includes forming a plurality of correlation data subsets using a plurality of sweep data sets and a correlation reference sequence. The method further includes estimating a noise level in a correlation data set using the correlation data subsets and subtracting the estimated noise level from the correlation data set.

Abstract: A method of seismic surveying comprising the steps of actuating the or each vibrator in a first vibrator group at time T0, and subsequently actuating the or each vibrator in a second vibrator group at time T1 that satisfies T0<T1<T0+S1+L where S1 is the sweep time of the first vibrator group and L is the listening time. At least one of the first vibrator group and the second vibrator group comprises at least two vibrators. The first group and the second group of vibrators may be the same group, or they may be different groups. This method enables the time required to complete a seismic survey to be reduced compared to the prior art “simultaneous shooting” and “slip-sweep shooting” techniques.

Abstract: Seismic prospecting method and device using simultaneous emission of seismic signals obtained by coding a signal by pseudo-random sequences, and notably a periodic signal phase modulated by such sequences. Signals such as, for example, periodic signals phase modulated according to a pseudo-random code (binary for example), produced by a control unit, are for example emitted in the ground by means of vibrators, the seismic signals reflected by the subsoil discontinuities in response to the periodic signals emitted are picked up by receivers coupled with the formation and recorded in an acquisition and recording system. The periodic signals are emitted simultaneously by the vibrators.

Abstract: A method of seismic surveying using one or more vibrational seismic energy sources activated by sweep signals. The highest order harmonic that has sufficient strength to cause significant harmonic distortion of a sweep segment is determined. A number of sweep segments in excess of the number of sources is selected. Initial phase angles are selected for each sweep segment of each seismic energy source so that substantially all harmonics up to and including the highest order harmonic are suppressed. Using cascaded sweeps, seismic data are acquired and processed to substantially attenuate harmonics upto the selected order.

Abstract: This invention relates to mapping the open fractures in the reservoir subsurface formations. To map these open fractures, the elastically nonlinear interaction of the two surface generated seismic waves, as they propagate through the open fractures and are reflected back from the subsurface boundaries, is recorded using conventional 3-D seismic receiver configurations. Recorded data are processed for standard 3-D reflection imaging and additionally for mapping the location and orientation of the subsurface open fractures.

Abstract: This invention relates to mapping the hydrocarbon reservoir characteristics by mapping the reservoir formations that display dynamic elastic nonlinearity responses to the seismic signals. The main reason of this nonlinear behavior in the reservoir rocks is their bulk rock property: the porosity, fractures, differential pressure and pore saturation. To map these bulk rock properties, the interaction of the two seismic waves as they propagate through elastically nonlinear rocks is recorded. Two compressional seismic signals are transmitted from the surface. Seismic reflection data using surface or borehole detectors are recorded. One of the transmitted signals is a conventional swept frequency and the other is a mono-frequency signal. During the propagation of these two signals through the elastically nonlinear reservoir rocks, there is an interaction between the two signals. The sum and difference frequencies of the two primary seismic signals that were transmitted from the surface are created.

Abstract: Near-offset and far-offset seismic data volumes are time-aligned by first selecting a plurality of time shifts. The near-offset and far-offset seismic data volumes are cross-correlated at the plurality of time shifts. An initial time-shift volume and a maximum correlation volume are created from the maximal cross-correlations at the plurality of time shifts. Areas of high time shift from the initial time-shift volume and areas of low cross-correlation from the maximum correlation volume are determined. The determined areas of high time shift and low cross-correlation are filtered from the initial time-shift volume, generating a filtered time-shift volume. The filtered time-shift volume is applied to the far-offset seismic volume to generate a time-aligned far-offset volume.

Abstract: Segmenting the spectral distribution of overlapped vibratory signals improves the efficiency of data acquisition while providing reduced harmonic distortion in the time zones of interest. Two identical sweep segments are used. Each sweep segment includes an earlier low frequency sweep and a later high frequency sweep, the individual sweeps having substantially no overlap in frequency except for tapering. The high frequency sweep in each pair starts before the end of the low frequency sweep with an overlap in time that is selected to avoid harmonics from the low frequency sweep. Correlation of the recorded signal separately with the low frequency sweep and the high frequency sweep gives data sets in which individual portions of the desired data are recoverable with the harmonic distortion largely separated from the desired data. The individual portions of the desired data are then spliced to give a broadband response of the earth.

Abstract: Methods are provided for automatic detection of linear seismic events and the detected events are used in correction of statics and geometry error.

Abstract: An apparatus and method are provided for detecting buries objects, geological formations or sediment properties by transmitting an acoustic signal and receiving a scattered signal. An acoustic source and an array of receivers are used to correlate the scattered signal with the transmitted signal.

Abstract: A method for identifying faults and stratigraphic features within seismic data without interpreter bias by processing seismic data to identify the minimum difference between seismic traces. Large values of difference are plotted as display attributes for seismic reflection data interpretation for two-dimensional and three-dimensional seismic data. The large values of difference represent faults and stratigraphic features within the seismic data. Dip azimuth and dip magnitude attributes can be generated and displayed as well.

Abstract: Harmonics, which appear with seismic reflection signals and are responsive to actual vibrations of a land-vibrator for petroleum exploration, are effectively utilized in pulse compression techniques by applying a correlation operator (CO) signal representative of the actual vibrations injected into the earth. In the invention, the CO signal may be acquired by detecting the motion of the baseplate of the land vibrator. Signal correlation procedures are then used to define an inverse filter whose output is pulse-like in response to the CO signal, and the reflected signal is applied to the inverse filter to yield a pulse compressed signal including harmonic energy.

Abstract: A reverse VSP system is provided and uses a frequency-controllable modulator of an MWD tool as a downhole seismic source. The modulator has a rotor and stator arranged to provide substantially sinusoidal outputs, and is controlled to provide a "Vibroseis"-type frequency sweep (preferably from 1 Hz to 50 Hz). With the bit of the MWD tool in contact with the bottom of the borehole, the force generated by the frequency-controllable modulator will generate measurable axial compressional (P-) waves. These waves which move through the formation are detected by an array of geophones on the formation surface, and a reference signal is obtained by a kelly-mounted accelerometer or strain gauge. The reference signal and signals obtained by the geophones are processed using known signal processing techniques.

Abstract: A method for separating and pre-processing vibratory source data includes varying the phase of vibratory sources according to two patterns. A first pattern is used for odd numbered sources while a second pattern is used for even numbered sources. Each pattern begins with a zero phase shift, the zero phase shift occurring at the sweep number corresponding to the position number of the source. The first pattern is one of zero phase shift, ninety degree phase shift, ninety degree phase shift and one hundred eighty degree phase shift. The second pattern is one of zero phase shift, one hundred eighty degree phase shift, ninety degree phase shift and ninety degree phase shift. The patterns are alternated for each source, the beginning of the pattern corresponds to the position of the source in a line. The first source begins with the first pattern. The second source begins with the second pattern, with zero starting at the second sweep.

Abstract: A method for recording and pre-processing vertical seismic profile high fidelity vibratory seismic data includes the steps of measuring the motion of the vibrator which is related to the vibrator applied force times a transfer function of minimum phase, causal, linear system relating the actual vibrator output with the measured vibrator motion, separating signals according to the generating source of each, determining a ratio by dividing the vibratory seismic data by the measured motion of the vibrator to remove the unknown applied force leaving the earth reflectivity times a time derivative divided by a minimum phase function, minimum phase band pass filtering the resulting ratio and performing minimum phase deconvolution to remove the time derivative divided by the transfer function of minimum phase.

Abstract: An instant acoustic logging method for use in oil drilling is disclosed in which a drill bit arranged at one end of a drill string and rotated in a well is used as the acoustic source. The method comprises the steps of measuring, at the top of the string, the vibration generated therein by the rotating drill bit, sensing acoustic signals from the drill bit at least one point on the ground away from said drill bit, and processing the signals to provide data on the rock around the well. The method comprises an additional preliminary step of processing the vibration measured at the top of the string to determine the nature of the pilot signals from the drill bit, the data on the rock around the well being obtained by comparing the acoustic signals sensed on the ground with the pilot signals from the drill bit.

Abstract: A method of separating the effects of the earth response on vibratory energy from individual ones of multiple vibrators to geophones, in generating a seismic survey, is disclosed. According to the disclosed method, measurements of the actual vibrations generated at the source are made, in addition to the geophone measurements of the refracted and reflected vibrations, both over a number of frequency sweeps. The source and geophone vibrations are first transformed into frequency domain representations, by way of Discrete Fourier Transforms. An inverse matrix of the recorded frequency-domain source vibrations is then generated for each frequency. When the number of sweeps equals the number of vibrators, the inverse matrix is merely the multiplicative inverse of the source measurements; if the number of sweeps exceeds the number of vibrators, the inverse matrix is the generalized inverse of the original matrix.

Abstract: A method of separating the effects of the earth response on vibratory energy from individual ones of multiple vibrators to geophones, in generating a seismic survey, is disclosed. According to the disclosed method, measurements of the actual vibrations generated at the source are made, in addition to the geophone measurements of the refracted and reflected vibrations, both over a number of frequency sweeps. The source and geophone vibrations are first transformed into frequency domain representations, by way of Discrete Fourier Transforms. An inverse matrix of the recorded frequency-domain source vibrations is then generated for each frequency. When the number of sweeps equals the number of vibrators, the inverse matrix is merely the multiplicative inverse of the source measurements; if the number of sweeps exceeds the number of vibrators, the inverse matrix is the generalized inverse of the original matrix.

Abstract: A process for preparing vertical seismic profiles (VSP) by bore drilling using, as the seismic source, the vibrations generated by the drilling bit. A cross-correlation between the bit signal and the signals collected by seismic geophone lines, is determined in which the bit signal is processed by arranging a plurality of several types of detectors on the bore drilling facility. A plurality of traces are recorded, in which the signal/noise components have different characteristics. These can be separated by means introducing filter functions by operating either on their Fourier Transforms or on the basis of their statistic distribution.

Abstract: A method for cascading or linking seismic vibrator sweeps to form a cascaded sweep sequence. According to the method, a first cascaded sweep sequence is generated containing N sweep segments linked end-to-end. The N sweep segments are substantially identical, except that the initial phase angle of each sweep segment within the sweep sequence is progressively rotated by a constant phase increment of about 360/N degrees. A second cascaded sweep sequence is generated consisting of (i) N consecutive sweep segments linked end-to-end which correspond to said first cascaded sweep sequence and (ii) an additional sweep segment linked to the N consecutive sweep segments which is positioned and phased so as to substantially suppress harmonic ghosts during correlation. One of these cascaded sweep sequences is used for the vibrator sweep sequence and the other is used for the correlation reference sequence.

Abstract: The present invention provides a method for on-line real-time processing of processing navigational data for determining the location of sensor and receiver points in a navigational network having a number of different types of devices. Observations from these devices are obtained using a coordinate system that follows appropriate nominal sailing lines. Outlying observations are discarded using w-statistics for the observations. Any correlated observations such as compass azimuths are uncorrelated. The uncorrelated observations are then sequentially processed in an extended sequential Kalman filter, which provides the best estimate of the station coordinates. These estimated coordinates are then used to determine the location of the source and receiver points.

Abstract: Seismic impulses irregularly distributed in time (according to a random or a pseudorandom code for example) are transmitted into the soil, and the direct waves and those reflected by the subsoil discontinuities are picked-up by a reception system (G). The frequency spectrum of the emitted and the received signals is divided into several bands, and the components of the signals in each one of the bands are reduced to elementary signals indicating the change of sign thereof (sign bit) and are recorded. Cross-correlation of the elementary signals respectively associated with the emitted and the received signals and stacking of the cross-correlation products associated with each sensor are carried out thereafter for each sensor of the reception system, and seismic sections very comparable to what would be obtained by digitizing the signals received with full precision, as in conventional methods, are achieved with the obtained stacks.

Abstract: A geophysical prospecting method comprises injecting signal energy in the form of at least two pseudo-random pulse sequences of the complementary sequence type into an underground formation. Next, the signal energy is detected after it has been attenuated by passing through the underground formations in the form of response signals. The response signals are correlated with reference signals comprising at least two pseudo-random binary sequences to yield correlations. The correlations are summed so that the resultant signal which corresponds to the subsurface response signal does not have secondary lobes. The method is applicable to seismic prospecting and more particularly to cross-well seismic surveying. It is also applicable to acoustic prospecting for logging measurements. In a preferred embodiment, the pseudo-random sequences comprise a pair of Golay sequences.

Abstract: An apparatus (10) for providing while drilling information on a subterranean geologic formation (28) includes a drilling rig (12) and a rotary drill bit (18) attached to the drilling rig (12) for providing seismic waves as it drills in the earth (16). Geophones (20) are spaced from the rotary drill bit (18) in the earth (16) and receive in direct seismic wave paths (26) and seismic wave paths (30) reflected from the subterranean geologic formation (28) the seismic waves provided by the drill bit (18). A reference sensor (24) is located on the drilling rig (12). The seismic signals sensed by the reference sensor (24) and by the geophones (20) are cross-correlated to separate the drill bit generated signals from interference signals by combining the reference signals and the signals received by the geophones (20).

Abstract: A method of deconcolving the non-ideal frequency response from acoustic vibrations transmitted along a structure such as a drill string is disclosed. The deconvolution retains the values of the transmission time between the signal source and the receiver, in the form of an exponential phase term, and is multiplied by the amplitude frequency response of the structure. The input data time series, after transformation into the frequency domain, is then divided by the deconvolution operator. The deconvolution method may be used in a noise reduction method where both axial and torsional vibrations are generated from the same location, where one of the time series is shifted by the amount of the time delay, so that the vibrations generated from the same location coincide, providing reinforcement of the desired signal.

Abstract: A method of correcting a seismic signal s(t) provided by a vibration sensor (100) in response to at least one emission signal F(t) produced by a least one seismic vibration generator (200). The real emission signal F(t) differs in phase and amplitude relative to a reference emission signal F(t). The correction method of the invention consists in determining the phase difference .DELTA..phi.(f) and the amplitude ratio A(f) between the real emission signal F(t) and the reference emission signal F(t) as a function of frequency f, and in applying the seismic signal s(t) from the vibration sensor (100) to a filter (304) whose response is defined in phase and amplitude respectively by -.DELTA..phi.(f) and by 1/A(f).

Abstract: An apparatus (10) for providing while drilling information on a subterranean geologic formation (28) includes a drilling rig (12) and a rotary drill bit (18) attached to the drilling rig (12) for providing seismic waves as it drills in the earth (16). Geophones (20) are spaced from the rotary drill bit (18) in the earth (16) and receive indirect seismic wave paths (26) and seismic wave paths (30) reflected from the subterranean geologic formation (28) the seismic waves provided by the drill bit (18). A reference sensor (24) is located on the drilling rig (12). The seismic signals sensed by the reference sensor (24) and by the geophones (20) are cross-correlated to separate the drill bit generated signals from interference signals by combining the reference signals and the signals received by the geophones (20).

Abstract: A method is provided for significantly reducing the distortion and crossfeed from any selected order harmonic for any number of vibratory seismic sources operated concurrently, at the same time providing for separation of the signals from the different sources and for improving the signal-to-noise ratio. After determining the highest order harmonic likely to cause distortion, a number of sweeps of each source in each position is selected. This number depends upon the number of sources and the highest order harmonic to be suppressed. Initial phase angles for each sweep of each source are then selected to permit signal separation while suppressing harmonics up to and including that highest order harmonic.

Abstract: One or more seismic sources impart seismic pulses into the earth in a coded sequence. The coded sequence is composed of a minimum of three cycles of a Galois sequence of order N, where N is a positive integer. The response of subsurface formations to the coded sequence is obtained and correlated with one cycle or the outgoing Galois sequence of the same order. Correlation properties of the Galois sequences provide a way to sum a large number of seismograms that are superimposed in time, improving the signal to noise ratio while keeping the field time in which the seismograms are obtained to a minimum. Using this encoding technique, responses of subsurface formations to concurrent activation of each of several sources can be obtained and separated from each other.

Abstract: A filtering system is provided for acquiring and processing signals using a sampled filter for signal separation and signal enhancement; such as for determing locations, distances, and times in a geophysical exploration system. Signature signals permit mixing and seperation with sampled filters, such as for sharing common circuitry and for increasing the amount of acquired information. A sampled filter, such as a digital correlator, is provided for generating high resolution output data in response to low resolution input data processed with low resolution computation circuits. A real-time time-domain correlator is provided with single-bit resolution computational elements to provide improved correclation filtering. A high speed real-time correlator is provided to enhance signals with copositing-after-correlation and with correlation using a plurality of correlation operators.

Abstract: An apparatus (10) for providing while drilling information on a subterranean geologic formation (28) includes a drilling rig (12) and a rotary drill bit (18) attached to the drilling rig (12) for providing seismic waves as it drills in the earth (16). Geophones (20) are spaced from the rotary drill bit (18) in the earth (16) and receive in direct seismic wave paths (26) and seismic wave paths (30) reflected from the subterranean geologic formation (28) the seismic waves provided by the drill bit (18). A reference sensor (24) is located on the drilling rig (12). The seismic signals sensed by the reference sensor (24) and by the geophones (20) are cross-correlated to separate the drill bit generated signals from interference signals by combining the reference signals and the signals received by the geophones (20).

Abstract: A method and apparatus for obtaining seismic vibrator reflection data is disclosed with a significant signal to noise ratio. The earth is sinusoidally vibrated with source vibrations linearly varying over time. The amplitudes of the reflections from the subsurface strata interfaces are sensed and multiplied by the amplitude of the source vibrations, resulting in a first representation of the product as a function of time. A second representation of the product is calculated in association with negative values of time, such that phase consistency between the two representations is maintained. The two representations are then mathematically combined, with values of zero associated with all other values of time. The two representations are then linearly transformed into a function of frequency, and loaded into an array for computer operation.

Abstract: Methods and apparatus for processing vibrative seismic data. Baseplate phaselock vibrative energy, reaction mass phaselock vibrative energy and ground force phaselock vibrative energy are transmitted into the earth using a land vibrator and the reflections of the transmitted waves recorded. The source signature of the land vibrator is computed for each type of vibrative energy and the source signature is correlated with the corresponding recorded data to produce an image of the subsurface. Regardless of the type of phaselock vibrative energy generated into the subsurface, the images of the subsurface produced by the correlation of the source signature and the recorded data are substantially the same. The correlated data may be then analyzed to yield useful information related to the characteristics of the subsurface.

Abstract: A seismic data processing system utilizes vibration type shock waves which are delivered by a reactive mass (52) acting in conjunction with a base plate (50) to output low frequency shock waves to the earth. The force signals for the reactive mass (52) and (50) are measured and input to a transmitter (42) for transmission to a receiver (48) at a control center (28) which also receives resultant measured data. The measured force signal is input to a data correlator (74) as a model signal and correlated with the resultant data to provide processed data out.

Abstract: An apparatus (10) for determining while drilling in the earth with a drill bit (12) the positions of seismic wave reflecting geologic formations (15, 17) has seismic wave sensors (A, B) positioned in the earth near its surface (42). The sensors (A, B) are disposed at known positions with respect to the borehole. The drill bit (12) generates the seismic waves. The sensors (A, B) receive the seismic waves both directly and by reflection from the geologic formations (15, 17). Adaptive filters (20, 26) are connected to each of the sensors (A, B). A summing circuit (34) receives outputs from the adaptive filters (20, 26). The converged impulse responses of the adaptive filters provide information about the amplitude of the reflected waves and their arrival times at the surface, thus allowing detection of the reflecting horizons and their location and physical characteristics.

Abstract: A continuous time-variant digital filtering method and apparatus. In particular, a method for continuous time-variant digital filtering of a digitized signal comprrises the steps of windowing the input digitized signal into discrete windows. For each discrete window of the input digitized signal, an autocorrelation matrix is determined. Thereafter, a filter coefficient matrix for each window of the input digitized signal is determined from the autocorrelation matrix. Corresponding filter coefficients from each filter coefficient matrix and window of the input digitized signals are sequentially collected and interpolated for each sample point of the input digital signal to form continuous filter operators. Each continuous filter operator is sequentially convolved with the input digitized signal to form a partial output signal and the partial output signals from the convolution of each continuous filter operator with the input digital signal are summed to form a final output signal.

Abstract: Vibrator seismic sources are operated concurrently to reduce the time required for acquiring the data for a seismic survey. Separate seismic records are produced for each of the concurrently operated vibrator seismic sources. Frequency sweeping pilot signals are generated for each of the separate vibrator seismic sources concurrently in a sequence of at least four sweeps. The pilot signals for each vibrator source have a plurality of separate phase angles during the sweeps. The seismic signals produced by the concurrently operated vibrator seismic sources are detected for each of the sweeps. The detected seismic signals are correlated separately with each of the corresponding pilot signals which were used to drive the vibrator seismic sources that produced the detected seismic signal. The correlating process produces first and second separate correlated records for each of the sweeps.