Wave propagation and sampling theory; Part I, Complex signal and scattering in multilayered media

Abstract From experimental studies in digital processing of seismic reflection data, geophysicists know that a seismic signal does vary in amplitude, shape, frequency and phase, versus propagation time. To enhance the resolution of the seismic reflection method, we must investigate these variations in more detail.We present quantitative results of theoretical studies on propagation of plane waves for normal incidence, through perfectly elastic multilayered media.As wavelet shapes, we use zero-phase cosine wavelets modulated by a Gaussian envelope and the corresponding complex wavelets. A finite set of such wavelets, for an appropriate sampling of the frequency domain, may be taken as the basic wavelets for a Gabor expansion of any signal or trace in a two-dimensional (2-D) domain (time and frequency). We can then compute the wave propagation using complex functions and thereby obtain quantitative results including energy and phase of the propagating signals. These results appear as complex 2-D functions of time and frequency, i.e., as 'instantaneous frequency spectra.'Choosing a constant sampling rate on the logarithmic scale in the frequency domain leads to an appropriate sampling method for phase preservation of the complex signals or traces. For this purpose, we developed a Gabor expansion involving basic wavelets with a constant time duration/mean period ratio.For layered media, as found in sedimentary basins, we can distinguish two main types of series: (1) progressive series, and (2) cyclic or quasi-cyclic series. The second type is of high interest in hydrocarbon exploration.Progressive series do not involve noticeable distortions of the seismic signal. We studied, therefore, the wave propagation in cyclic series and, first, simple models made up of two components (binary media). Such periodic structures have a spatial period. We present synthetic traces computed in the time domain using the Goupillaud-Kunetz model of propagation for one-dimensional (1-D) synthetic seismograms.Three different cases appear for signal scattering, depending upon the value of the ratio wavelength of the signal/spatial period of the medium.(1) Large wavelengthsThe composite medium is fully transparent, but phase delaying. It acts like an homogeneous medium, with an 'effective velocity' and an 'effective impedance.'(2) Short wavelengthsFor wavelengths close to twice the spatial period of the medium, the composite medium strongly attenuates the transmission, and superreflectivity occurs as counterpart.(3) Intermediate wavelengthsFor intermediate values of the frequency, velocity dispersion versus frequency appears.All these phenomena are studied in the frequency domain, by analytic formulation of the transfer functions of the composite media for transmission and reflection. Such phenomena are similar to Bloch waves in crystal lattices as studied in solid state physics, with only a difference in scale, and we checked their conformity with laboratory measurements.Such models give us an easy way to introduce the use of effective velocities and impedances which are frequency dependent, i.e., complex. They will be helpful for further developments of 'complex deconvolution.'The above results can be extended to quasi-cyclic media made up of a random distribution of double layers. For signal transmission, quasi-cyclic series act as a high cut filter with possible time delay, velocity dispersion, and 'constant Q' type of law for attenuation. For signal reflection they act as a low cut filter, with possible superreflections.These studies could be extended to three-dimensional (3-D) binary models (grains and pores in a porous reservoir), in agreement with well-known acoustic properties of gas reservoirs (theory of bright spots).We present some applications to real well data. Velocity dispersion may explain: mistying between sonic logs and velocity surveys; mistying between synthetic seismograms and seismic sections.Effective velocities lower than mean velocities may explain: very low velocities for P- and S-waves, especially in weathered shallow layers; anomalies on the values of the ratio V s /V p ; mistying between P and S seismic sections.Finally, the Gabor expansion provides a tool to obtain sampled instantaneous frequency spectra, and to carry out a suitable recording and processing method in high-resolution seismic, especially to preserve the phase information. Such a processing will involve complex signals, complex traces, complex velocities and complex impedances.For practical purpose, this paper comprises two separate parts. Here, we present some interesting features on the scattering of seismic signals obtained by a simulation method using simple models. We show the usefulness of the notions of complex signals, complex velocities, and complex impedances, and overall of the Gabor expansion, by simulation on very simple models.Morlet et al, (1982, this issue) will be concerned with the development of fundamental notions useful to handle seismic data from the sampling method of recording to processing methods. There we give theoretical and practical tools to sample and handle these data in the time-frequency domain, using complex functions.