SYSTEM FOR MARINE SEISMIC EXPLORATION

Abstract: System for seismic profiling of structures (2) beneath the sea bottom, comprising an acoustic transducer array (12) which synchronously and coherently generates an acoustic wave adapted to be reflected from said structures (2). When the reflected acoustic waves propagate towards the water surface (1), there is caused a perturbation of the electromagnetic reflection coefficient thereof. Receiver means for detection and characterization of the received acoustic field at the surface (1), comprises a radar system (10) illuminating the region concerned of the surface (1) and imposes thereat a number of interferograms each corresponding to the acoustic field from a scattering point given by the position thereof and the frequency of the field. By recording the degree of coincidence between the interferograms and the acoustic perturbation pattern on the water surface (1), the seismic properties of said structures (2) are surveyed.

Description: System for marine seismic exploration. Known methods for marine seismic exploration are based on the emission of acoustic waves in the water from acoustic sources in the form of explosive charges, air guns and the like, whereby the corresponding acoustic signals from the depth are received by means of hydrophones, usually arranged i a seismic streamer which is towed behind a vessel. This invention provides a completely new concept for such seismic explorations based on the merging of several more or less well-established fields of technology and science today, which gives a high degree of flexibility and makes it possible to design an optimum system for seismic explorations. Some main aspects of this new system are as follows: New principles in radar technology involving matched radiation (illumination) make it possible to take advantage of apriori information regarding the object, i.e. the structures or subsurface layers under investigation. Based on such information the frequency bandwidth of the system can be narrowed and the sensitivity correspondingly increased. The aim of such a "self-learning" system is to have the appropriate seismic signatures stored in the transmitter unit processor and to structure this in such a way that only two bits of information is required: Is the particular structure present or is it not. A specific aspect of the invention involves that optimum illumination is achieved by coupling a set of sinusoidal acoustic waves to the wave propagation medium (i.e. water). By using matching networks in the form of resonant cavities, like what is known from microwave circuitry, there can be achieved a very high degree of matching over a narrow bandwidth. Interaction mechanisms between acoustic and electromagnetic waves are well known in solid state physics and are utilized within electronics as a means of achieving high speed signal processing methods. As examples can be mentioned surface acoustic wave devices for implementation of filters and convolvers. Further can be mentioned opto-acoustic spectrum analysers based on surface acoustic waves in solids. Thus, the present system is primarily designed for offshore oil and gas explorations and is based on the modern high- technology concept mentioned briefly above. More specifically the invention relates to a system for marine seismic exploration of structures under the sea bottom, comprising transmitter means with an array of acoustic transducers for synchronously and coherently generating one or more monochromatic acoustic waves in the water, and receiver means for detecting acoustic waves reflected from these structures. What is novel and specific to the system according to the invention in the first place, consists therein that the receiver means comprises a radar system adapted to illuminate that area of the sea surface in which the' reflected acoustic waves appear, in order to therein establish a number of inter- ferograms each corresponding to the acoustic wave field from a scattering point in the structure given by the position of the point and the frequency of said field, and that said receiver means is adapted to determine the degree of coincidence between said interferograms and acoustic perturbation patterns affecting the electromagnetic reflection coefficient of the sea surface, thereby to obtain a recording of the seismic properties of said structures. As will appear from the following description, much of the theory and explanations given apply both to electromagnetic waves (the radar part) and acoustic waves (the acoustic part) as these are incorporated in. the new system described here. The invention shall be explained more closely in the following with reference to the drawings, of which: Fig. 1 illustrates the basic concept of the system according to the invention, Fig. 2 schematically shows an acoustic source arrangement for providing a narrow multifrequency beam in the water, to be used in the system according to the invention, Fig. 3 schematically shows the principle of radar detection of acoustic wave patterns on the water surface, resulting from reflection and/or refraction of an emitted acoustic beam, Fig. 4 is a diagram showing the formation of a surface wave pattern resulting from acoustic waves from the depth of the water mass, Fig. 5 shows the scattering of an electromagnetic wave, which illuminates a water surface with an acoustic wave pattern of the kind illustrated in fig. 4, Fig. 6 illustrates coherent scattering from an area of the water surface having an acoustic wave pattern and being illuminated with electromagnetic waves, Fig. 7 shows more in general as a spatial diagram, which magnitudes are involved in and determine the scattering of electromagnetic waves caused by an acoustic wave, Fig. 8 shows determination of scattering angles by geometrical construction, Fig. 9 shows in two diagrams the spatial resolution in a pulse system (time domain) , respectively in a multi¬ frequency system (frequency domain) , Fig. 10 illustrates wave number matching based on the carrier frequency itself and a beat frequency pattern respectively, Fig. 11 shows matched illumination making use of apriori information about the wave propagation medium, Fig. 12 schematically illustrates how any scattering object characterized by a distribution in space of scattering elements, can be decomposed into its Fourier- components, Fig. 13 shows matching to periodic wave structures (patterns) , Fig. 14 shows matching to a frequency modulated wave structure, Fig. 15 illustrates the use of and the advantages of having multiple apertures in the system, Fig. 16 shows the effect of focusing the acoustic and/or the electromagnetic system, Fig. 17 illustrates the scanning of an area with movable focal points, Fig. 18 shows how a transverse spatial interference pattern gives information about transverse shape and size of a scattering object, Fig. 19 illustrates the geometry of the back scattering process, Fig. 20 serves to explain the relationship between the extent of an antenna or transducer array and the resolution obtained, Fig. 21 shows a bistatic synchronized acoustic system with several controlled acoustic transducer arrays, Fig. 22 gives the conditions for equal resolution transversally and longitudinally, Fig. 23 shows a form of coherent acoustic point source. THE BASIC CONCEPT This is shown in figure 1. Beneath a water mass the surface of which is denoted 1 there is an underground structure 2 which for example can consist of several layers as indicated.. A number of acoustic sources or generators 18 for acoustic sinusoidal waves with frequencies chosen so as to provide a "matched" radiation or illumination for the geological structure 2, feed one or more transducers 12 so that acoustic waves propagate through the water and further on to the geological structure 2. The different frequencies can be transmitted simultaneously or in succession. By diffraction and/or refraction the acoustic signal is scattered back from irregularities in the geological structure 2. These acoustic return signals are sensed either by a set of conventional microphones (hydrophones) or by means of a radar system generally designated 10, the radar wave length of which is tuned so as to get a constructive coherent interaction with the acoustic waves caused on the water surface 1 by the return signals from the depth of the water. (The socalled electro- acoustic or electro-optic effect). In order to achieve this, the acoustic generators are controlled by means of a radio telemetry system which provides for phase-locking of the acoustic part to the radar part or the radar wave pattern. More specifically the acoustic part comprises a control unit 11 with a receiver antenna 11a which from a radio telemetry unit 13 with transmitter antenna 13a receives control signals as indicated with arrow 9. On the other hand the unit 13 is connected to radar transmitter means 14, which inter alia comprises phase shifters and which through one or more antennas 16 illuminates the water surface 1 with radio waves for sensing the acoustic waves as mentioned above. The returning radio waves are received by a radar antenna 17, which is connected to a signal processor 15 having functions to be explained more closely in the following. The acoustic generators 18 with associated transducer(s) , which advantageously constitute a transducer array, are thereby all controlled by the radar system 10, which illuminates the water surface 1 so that there is provided a phase-locked coherent radar "hydrophone" system. From the receiver antenna - possibly antennas - 17 the signal processor 15 receives information with respect to amplitude and phase of each composite acoustic wave and thereby can provide information regarding the transfer function of the geophysical medium (mass of water and underground structures) , by which the acoustic waves have been affected. The Fourier transform of this transfer function directly gives the desired delay function, i.e. the "echogram". Such a system can operate in a self-learning mode, whereby the radar system estimates the condition of the water surface 1 and selects frequencies for the acoustic field and for the radar or interference field, which gives a maximum of contrast for the acoustic field with respect to the water surface. More particularly there can be implemented algorithms, which enhance the surface wave pattern (footprint) of the acoustic fields in relation to the background represented by the water surface. These algorithms can for example be based on Doppler-filtering and rigidity filtering. Specific acoustic sources In the principle known forms of acoustic sources or transducers can be used in connection with the present new system, but these known designs have a poor efficiency and are not particularly well suited. As mentioned, there is proposed an optimum acoustic coupling in connection with the new system, this coupling being achieved by resonant cavities tuned to the individual matched frequencies. These acoustic frequencies are phase-locked to a corresponding set of radar frequencies, as will appear from figure 1. At the end of this description the particular form of acoustic sources referred to here, shall be discussed more closely. Another important aspect is that concentration of the acoustic beam is achieved by a cooperating group or array of acoustic sources. Such a source array can be composed by means of a number of remotely controlled boats having acoustic sources in their bottom, possibly by having the boat bottom itself made in the form of "tuned drums", i.e. acoustic sources with resonant cavities as mentioned above. The relative position of these source array elements are measured accurately by means of radio phase methods. Figure 2 illustrates schematically such a group or array of acoustic sources for providing a narrow multifrequent acoustic beam in the water. There are shown four boats 24, 25, 26 and 27, of which boat 24 is provided with an antenna 24a for-its control. Boat 24 carries 2 acoustic transducers 21 and 22 for transmitting frequencies Fl to F5 under control by radio from a main station. Boat 27 is provided with an acoustic transducer 23. The boats are shown with different spacings so adapted that from the total array of acoustic transducers on the boats there will be emitted a single acoustic beam 29 having a beam angle A /D where D is the maximum extent of the source array elements. In such an arrangement, as shown in figure 2, there is obtained a beam sharpening factor of D/d, where d is the source or transducer aperture, as shown for transducer 23 on boat 27. In view of the above it will be understood that the acoustic transducers can be arranged in an array which is one- dimensional, two-dimensional or three-dimensional for the purpose of providing harmonic, multifrequent signals. This arrangement of acoustic transducers on independent boats or vessels, the positions of which are determined by using radio navigation methods of high accuracy, makes it possible to compensate for deviations from an ideal position with respect to the other transducers in the array. The principles of a radar "hydrophone": Figure 3 additionally shows relationships as a supplement to the principle illustrated in figure 1, for explaining the so¬ called radar "hydrophone". A vehicle 30, for example an airplane, carries radar antennas 31 and 32 for transmission and reception respectively, whereby the radar illumination is determined by generators 33 for a number of frequencies Fl to F5. When transmitting acoustic waves from an acoustic source array established by means of 4 boats 35, 36, 37 and 38, there is provided at the water surface a resulting acoustic illumination or an acoustic surface wave pattern 39 (footprint) , resulting from reflection and/or refraction of the acoustic beam transmitted towards underground structures of interest beneath the sea bed. As will appear from figure 3, the radar hydrophone is matched to the acoustic waves so that there is provided a maximum of contrast between the acoustic field and the natural wave phenomena or patterns, which might be present on the water surface. Thereby information about amplitude and relative phase for each frequency involved is obtained. Computing the transfer functions V(F) as indicated in diagram 3B on figure 3, gives information about the spatial distribution of the scatterers constituting the scattering medium, by a Fourier transformation. An alternative to this description is that the radar hydrophone transmitter antenna has a fixed position on shore. The resulting delay function is shown in diagram 3C in figure 3. Such a system will have an inherent extreme sensitivity as a result of the following 4 factors: - The acoustic illumination from the source array on boats 35 - 38 is matched to the geological structure of interest, thus reducing the required wave number(K) spectrum to a minimum. - The multifrequency radar hydrophone is matched coherently to the surface wave footprint 39 of the set of coherent sinusoidal acoustic frequencies employed, over a large surface area of the sea, which gives a very large "sensing area" being coherent in space and time. The Doppler shift (diagram 3A in figure 3) , which the scattered radio wave is subjected to as a result of the velocity of the acoustic surface wave pattern or footprint, is apriorily known. The width of the Doppler filters, i.e. one filter for each of the acoustic frequencies used, can thus be made very small, which results in an extremely low noise bandwidth. — The velocity of the acoustic wave (VA ) is larger by several order of magnitude compared to the velocity of the gravity waves on the sea surface and the incoherent motion caused by the local wind field. Therefore these will not or will only to a very small degree disturb the detection process. INTERACTION MECHANISMS BETWEEN ACOUSTIC AND ELECTROMAGNETIC WAVES Acoustic waves (sound waves) As already mentioned previously and as illustrated in figure 4, acoustic waves 40 reflected from the sea floor and from structures underneath, will produce a pattern (footprint) on the water surface 41. This surface wave pattern can be detected by means of radio waves. The pattern is periodic provided that the field is a plane wave of a single frequency. The wave length is denoted by λ a . It appears from figure 4 that the periodJLof the surface wave pattern on the water surface 41 depends on the angle . It is evident from the figure that the period Λ of the pattern is larger than the wave length. Thus, the pattern on the surface travels with a velocity that is larger than the velocity of the wave itself, which is roughly 1500 meters/second. This acoustic wave, which comes from below, causes the water surface to deviate from its "ideal" plane level so as to take the form of a sine curve 44 as shown dotted in figure 4. The amplitude h of this curve depends on the power carried by the acoustic wave. If, for instance, this power is 1 milliwatt per m 2 , the amplitude is about 1 micrometer. As mentioned, there will normally exist ordinary ocean waves (gravitation waves and capillary waves) on the water surface 41. Some of these waves may easily have amplitudes of several meters. It may still be possible to detect the acoustic surface waves or sound waves 44, because they have a behaviour quite different from that of ordinary ocean waves. If the acoustic frequency employed is 30 Hz, which is an attractive frequency for seismic explorations, the wave length λa will be 50 meters. The resulting surface wave pattern 44 may then have a period depending on the angle α, of about 70 meters and travels at a velocity of 2100 meters/second. In contrast a gravitational wave of wave length 100 meters has a frequency of 0,12 Hz and travels with a velocity of 12 meters/second. Radio waves By assuming that the water surface (as already mentioned) is illuminated from above by a radio wave 42 (see figure 5) , which is also a pl-ane wave of a single frequency and with a polarization as shown, the interaction between the acoustic surface wave pattern and the electromagnetic waves can be studied more closely. From figure 5 it appears that the angle of incidence of the radio waves 42 is γi . Moreover, reference is made to figure 6, in which there is indicated an area A = i x w of the water surface being illuminated from above with the radio wave 42 and from below by the sound wave 40 (figure 4) . It is assumed that the dimensions I and w of this area are large compared with the wave lengths λ and Λ • Most of the incident electromagnetic power will be reflected in ordinary manner. However, a small part of the incident radio wave power will be scattered coherently by the tiny perturbations or acoustic surface wave patterns 44 being generated on the surface 41 by the acoustic sound field 40. In figures 5 and 6 this back-scattered portion of the incident radio wave 42 is designated 43. This has a scattering angle designated γ s and a frequency ω '. The amplitude of the back-scattered wave at a distance R from the scattering area A is proportional to the quantity: where Pa is the acoustic power density (power per square meter) and Pen is the electromagnetic power density. It is worth noting that the scattered field is proportional to the square root of the total acoustic power directed into the area A. Likewise, the field is proportional to the square root of the total electromagnetic power. The sensitivity of the radar hydrophone will thus increase with increasing power in the radio wave employed. Figure 7 with associated mathematical expressions show additional relationships in the scattering of an electromagnetic wave caused by an acoustic wave 49. This represents a generalization with respect to the preceding discussion, which for simplicity was restricted to the case when the plane of incidence of the radio wave also contains the wave vector of the acoustic wave. In figure 8 the projections of the wave vectors are —> -* denoted ki and k s respectively. The scattered wave k s has a frequency equal to the sum of the incident wave and the acoustic wave's frequencies. Another scattered wave, which is not shown in the figure, may also be present, with a frequency equal to the difference between the same frequencies. Figure 9 shows how- the scattering angle in the horizontal plane, Φ s , can be determined by a simple geometrical construct¬ ion. The angle Υ s . formed by the wave in relation to the vertical direction, is then determined by the length of ks being known and equal to ω/c . In practice the acoustic wave is not always a plane wave. The wave fronts may be curved and the apparent wave length may change along the direction of propagation. Modification and adaption in view of these problems may be made by experts on the basis of the general theory presented here, without any necessity for discussing this more closely here. From the above it also appears that on the basis of knowledge of the frequency of the acoustic field and by using only one interference pattern, there can be provided information simultaneously with respect to the scattering coefficients for a number of depth layers or underground structures. In this are involved the local angle of incidence and the phase velocity of the acoustic waves and requirements with respect to matching of wave number and velocity. MULTIFREQUENT MATCHED ILLUMINATION Given all degrees of freedom one can visualize many schemes for enhancing the total scattering from an object of interest (see surface,- underground structure) at the expense of the contribution from the background against which the object is viewed. In connection with this invention, however, one is facing the rather general problem of matching a radiation or an illumination to four dimensions (3 in space and 1 in time) , and it is of importance to offer a unified method in the sense that all four dimensions are addressed in the same manner. Seeking a system where the four dimensional tailoring of the illuminatio is individual for each dimension, will at best complicate the issue. It is also likely that such particular solutions would prove to be suboptimal. Accordingly, the problems shall here be solved in general terms, considering the scattering object (surface) to be characterized by a four-dimensional irregularity spectrum (wave numbers spectrum) . This approach is based on well-founded generalized Fourier theory, such as the classical theory presented by Wiener (1930) . This involves describing the scattering surface as decomposed into its Fourier components to which individual continuous and sinusoidal electromagnetical waves are matched. It will easily be seen that this involves a superposition of a set of individual phase-locked sine-waves (Fourier components) from individual electromagnetic oscillators, so as to produce a wave-form that is matched to the structure (Fourier components) of the scattering body or object. For this approach to be meaningful in relation to practical objects, one must confine oneself to an object characterization based on the diffuse back-scattered component. This means that the scattering surface must be illuminated with waves, whose wave length is small compared with the roughness scale of the scattering object (-surface.-structure) . For readers who are not quite familiar with Fourier domain operations and corresponding matched illumination wave number matching, it may be of interest to present an approach from a different angle. For this purpose figure 10 is presented. In this figure diagram 10a shows certain relationships in the time domain and diagram 10b corresponding relationships in the frequency domain. The spatial resolution of a pulsed system is determined by the pulse width (ΔZ = Cτ) . in order to obtain the same spatial resolution by means of a multifrequency "interferometer" system, a bandwidth ΔF = C /Δ Z is required. By resolving the rectangular pulse into its Fourier components (coefficients in the Fourier expansion) , the basis for an equivalent multifrequency system is acquired. The rectangular pulse is composed from a set of continuous wave oscillators with frequency and amplitude determined by the coefficients of the Fourier expansion. It is to be noted that the receiver noise bandwidth of such a multifrequency system is very much reduced i relation to the corresponding pulsed system. 12 a What is described here has some relation to the method proposed by present inventors Gjessing and Hjelmstad in International Patent Application PCT/NO86/00003. With advantage the method described therein can be employed also in connection with the present invention. The wave number matching mentioned above shall now first be discussed in a qualitative, semi-intuitiv manner based on simple principles of physics. Afterwards a short mathematical approach shall be presented. Wave number matching For detecting or identifying object structures by means of wave number matching (spatial resolution) , the object is illuminated with a set of electromagnetic or acoustic waves having different but mutually correlated frequencies. The properties of the back-scattered signal are observed and from 13 this there can be drawn very definite conclusions about the distribution of the scatterers constituting the scattering object or body along the direction of radio wave propagation. With reference to figure 10 a simple object is first considered, for example a ship 52 being represented by two identical scattering centers or reflection points having a spacing of Δ z longitudinally along the direction of wave propagation. From an antenna 51 this object is illuminated with a wide spectrum of frequencies, and the set of interfering waves scattered back from the object is analyzed (inversion techniques) . Two illuminating waves can be considered. As the difference in frequency between the two waves increases, a situation arises where the waves reflected from the two reflection points arrive in antiphase at the receiving site, causing a minimum in the resulting interferogram. If the carrier wave frequency itself is taken as a basis (diagram 10A in figure 10) , a carrier frequency given by F = c/2 Δ z gives constructive interference from two scattering centers spaced Δz apart. In exactly the same way a beat frequency pattern (diagram 10B) can be matched to the reflecting structure by choosing two illuminating sine-waves with a spacing A F = c/24 z, i.e. by choosing a Δ K of the illuminating wave, which is matched to a wave number of the scattering structure. If a complicated reflecting object is concerned, being characterized by several pairs of scattering centers (several Fourier components) , the same number of frequency pairs are needed for the ΔK-matching or several frequencies for the K-matching. If there exists apriori information about the sub-surface structure to be explored and it is desired to ignore certain information, the illumination can be matched to the known part of the propagation medium. This is illustrated in figure 11. In this figure there is shown an example of the use of three frequencies from generators 53 with corresponding radiation apertures 53a (antennas or acoustic transducers) illuminating a structure 55 consisting of several layers (model) . The number (M) of different frequencies employed is equal to the number of layers in the case when a model is concerned. If there is no model, the number of resolution cells required is equal to 14 M = Δ Z/ δZ. Figure 12 schematically illustrates how any scattering object 56 characterized by a distribution in space of the scatterers ( σε (r) ) can be decomposed into its Fourier components A(K) . There will be constructive interference from all these scatterers if there is chosen a combination of illuminating sine-waves ET (K) that couple to each of the Fourier components A(K) of the scattering object. There will be matched -. -r illumination when E(K) = A*(K) . A very important aspect of the matched illumination concept is the radar's ability to optimize the contrast between the useful signal provided by the acoustic field, and the background clutter (sea waves, turbulence etc.) in a self- learning mode. The radar measures the sea state and the detailed behaviour of the scatterers on the sea surface and avoids or compensates for the components in the 4-dimensional spectrum that are due to natural wind and wave fields. Focusing Improved detection and identification of structures are possible by employing focused acoustic and electromagnetic radiation arrays. As mentioned previously, it is possible simultaneously to obtain wave number and wave velocity matching by using suitably selected discrete frequency lines. Figure 13 shows matching to periodic wave structures obtained by simply selecting frequencies (or frequency spacings) for matching to the acoustic wave length. In applications where the structure of interest is moving with a velocity V with respect to the radar (e.g. when the structure is the acoustic wave field at the air/sea interface migrating with a velocity V = C s /sinus where C 3 is the velocity of sound) , a 4-dimensional matched filter can be realized through the use of phase shifters. By spinning the phase of the transmitted waves with a rate f = 2V/cF, the resulting hologram (interference pattern) will move or migrate at a velocity equal to the velocity of sound. 15 By using time-varying frequency components it is furthermore obvious that matched filters can be realized for structures with varying wave number versus position. This is explained by means of figure 14 where a linearly changing frequency versus time provides an interferogram with continuously varying local wave number. In the same way, using multiaperture, multifrequency techniques, any orientation of the resulting holograms can be obtained. This is illustrated in figure 15, wherein as an example there are shown two apertures 61, 62 (antennas or acoustic transducers) driven by a separate generator or signal source. This figure also contains a numerical example. Thus, more particularly the radar part of a system according to the invention can have a number of apertures located side by side for generating wave number vectors in any direction. From the above it appears that the focused matched filters can be realized in connection with the invention and that this gives a much improved sensitivity of the system. Figure 16 indicates how the focusing can be effected for the acoustic and the electromagnetic part of the system respectively. In the separate diagrams A, B, C and D in figure 16 water surface is designated 71 and a sub-surface object 72, whereas the sea bed is designated 75. The position of the radar is shown at 70, whereas corresponding acoustic sources or source arrays respectively are shown at 77a, 77c as well as 77b and 77d respectively. The various combinations of focusing serve to improve the detection and mapping performance. More particularly diagram A indicates a situation in which none of the system parts are focused. In diagram B the acoustic system beam energy is focused at specific points in the sub-surface structure 72. The radar part here operates without focusing. In diagram D the radar part is adapted for synchronously scanning with the acoustic source array. Diagram C is based on a focused radar part scanning the area illuminated with the acoustic part 77 c. Figure 17 shows in the horizontal plane 18A and in vertical section 18B respectively, how a focal point 88 for the combined acoustic-electromagnetic system can be scanned with 16 high speed. The focal point 88 can be located at an arbitrary distance from the radar antenna array 81 r - 82 r or the acoustic source array 81 a - 82 a. The focal point 88 can be directed to any structure of interest in the sub-surface formations beneath the sea bed 85. The sea surface is shown at 80. From the above it appears that one of the system parts (the radar part) can be stationary and remotely located, whereas the other part (the acoustic part - source array) can be movable. Moreover, it has appeared that the illuminated area (footprint) has large dimensions compared to the wave length of the acoustic interference pattern, so that there is obtained a focusing effect, and that there is employed a number of frequencies for matching to acoustic waves having their phase centers directly beneath the position of the radar. Advantageously the radar can be provided with more than one aperture and can operate with time varying frequencies thereby to generate interference patterns corresponding to a focused acoustic wave, the phase center of this acoustic wave being able to have an arbitrary position in three dimensions in relation to the radar platform or position. SPATIAL CORRELATION OF ACOUSTIC WAVES Earlier in this description there has been a discussion of the correlation properties of acoustic waves having different frequencies. By measuring the degree of correlation between waves having different frequencies, there is obtained informatio about the longitudinal distribution of the scatterers. When dealing with a thin reflector (zero distribution in depth or length) , then the bandwidth of the reflector is very large. Conversely, if the scatterers are distributed over a large region in space, the bandwidth is small. The attention shall now be focused on the transverse distribution of the scattering elements constituting the scattering object. In order to reveal this transverse structure use is made of another characteristic property of acoustic 17 waves, namely their spatial correlation properties. This forms the basis for a synthetic aperture method. The scattering object shall be illuminated with one single frequency, and at the receiving site there is used a set of antenna elements distributed along a base-line or in a plane which is perpendicular to the line joining the receiving antenna array and the scattering object. At each element the amplitude and the phase of the impinging wave are measured. On the basis of these point observations of field strength, information about the object or structure can be extracted. This shall be explained more closely in the following. Figure 18 illustrates the scattering process. Provided that the position of an object 92 is known, the transverse spatial interference pattern will give information about the transverse shape and size of the scattering object. The surface of the object is illuminated with a single wave, for example from an antenna 91 having a drive or illuminating unit 90. Depending on the properties of the surface, there is obtained a certain field strength distribution of the back- scattered field E s (x) . 94a and 94b designate the phase front from 92a and 92b respectively.. The basic principle of spatial interferometry (based on a synthetic aperture antenna) further appears from figure 18. Measuring the transverse spatial distribution of intensity takes place in a processor unit 95 with an accompanying antenna array. The processor unit 95 gives information about the transverse shape and size of the scattering object if the range thereof is known. Now follows a calculation of the angular distribution P(Θ) scattered from a surface. It has been deduced earlier (Gjessing 1978a) and besides is intuitively rather evident that the field strength of the back-scattered wave from an object characterized ■ → by the delay function f (z) is given by: v<|» E(K) ( 5z) . -e-JK- 2 d 3 £- (1) 18 From the above equation the following expression is found for the back-scattered field expressed by the intensity distribution over the scattering surface. E S (K) E(x) e (2) -^ 'X dx since K = . -- sm 2ττ (3) (for small angles] this equation shows how the back-scattered field is distributed in the direction Θ. From this the angular power distribution -=7 P(K) is deduced as follows (Gjessing and Boerresen 1968) : P ( K ) • ** E ( ) E S *(K) E*(x) ECx+r) e^ K'x β -JK(x+r ) d 3 3 d 3 ? (4) The second integral above is immediately recognized as the —> spatial autocorrelation RE (r) of the field strength distribution. Hence P( ) -v. P(sin a) R-(r) r^ κ ' r d ? (5) This equation shows that the angular power spectrum (radiation pattern.__ of the scattered wave is the Fourier transform of the intensity distribution over the scattering region when this distribution is expressed statistically in terms of its spatial autocorrelation, as illustrated in figure 19. Thus, figure 19 shows the geometry of the back-scattering process. The illuminating field gives rise to a scattered intensity distribution Es (x) . The above relationship is well known from antenna theory: The radiation pattern or diagram (the angular power distribution) of an antenna with aperture A is obtained by the Fourier 19 transform of the intensity distribution over this aperture. Thus, if the above function E(x) is rectangular, which implies that the intensity is evenly distributed over the antenna aperture, then the angular power distribution is of the form (sinus Θ/Θ ) 2 and the beam-width β = λ/D where D is the aperture size. Conversely, if the intensity distribution over the aperture is of the form sinus x / x, then the angular distri¬ bution of the back-scattered wave will be rectangular. Equation (5) above shall now be considered again and used as a basis for studying another important property of the scattered field, namely the spatial correlation of intensity. In the preceding discussion the case of a "transmitting antenna" was taken as a starting point. Now the situation shall be considered at the receiving side. The receiving antenna consists of a set of antenna array elements which makes possible measurements of amplitude and phase at each element in the antenna array. The power reaching this array antenna is distributed as P(Θ) over an angular region. Applying the inverse Fourier transform and besides by intuitive consideration, it can be found that information about the spatial correlation properties of the field strength is obtained as follows: jK-r -(r) P(K) dK (6) This equation says that the spatial correlation of the scattered field is the Fourier transform of the angular power distribution. From figure 18 it is seen that if an object of transverse extent Δ x shall be resolved by means of a receiving antenna array at a distance R from the object, it is necessary to measure the intensity distribution over a spatial region Lx = R \/Δx. Summing up the above, it can be remarked that by measuring the intensity distribution across a broadside array (measuring amplitude and phase at each antenna array point) , direct information is obtained about the transverse scattering properties of the scattering object. 20 This is illustrated in figure 20. Therein is shown (figure 21A) an antenna array consisting of antenna 96, 97, 98 and 99 with an accompanying interferometer processor unit 95. If the antenna array 96 to 99 consists of n elements with unequal spacing, there will be N = n(n-1) points on the spatial 2 covariance function (angular power spectrum) . If the total extent of the antenna array is D, the resolution of this system will beΔθ = λ/D. If phase can be measured with a resolution of ø degrees, the angular resolving power becomes: Δθ = φ /360 λ/D (7) Within the sector Θ there will be no ambiguities Θ nonambigous = λ/d (d is the smallest spacing between elements . ) It will be realized that there is analogy between this spatial autocorrelation function and the autocorrelation function in the frequency domain as discussed above. In the multifrequency case there can be "filtered out" certain longitudinal spatial distributions of the scattering object by using frequency filters. With a broadside antenna array it is possible to "filter out" certain transverse properties of the scattering elements by providing a "spatial filter". With a multifrequency multi-aperture system the two- dimensional Fourier transform (the hologram) of the scattering surface will be obtained directly. If by matched illumination there is produced a hologram of the structure of interest, there is obtained a system with optimum sensitivity and identification capability. In connection with the present discussion of the angular power spectrum of back-scattered acoustic waves, both the transmitting side and the receiving side have been considered, in order to best illustrate the physics of the problem. This treatment is based on the common and known assumption of reciprocity. 21 In connection with the above it may be of interest to consider a bistatic seismic system which is very much analogous to a bistatic radio wave scatter system. Here no use shall be made of any radar hydrophone as discussed earlier, but rather a set of tuned receiver hydrophones for continuous waves, these being almost identical to the transmitter arrangement based on "tuned drums" mentioned previously. There are employed a number of small remotely controlled boats like in the transmitting case. The receivers are synchronized with the transmitter generators by radio telemetry. This makes it possible to use receiver elements with extremely small bandwidth and correspondingly high sensitivity. Furthermore, by keeping track of absolute phase for all the frequencies at all the array element points, a dedicated high resolution mapping of the geological structure can be accomplished. This is illustrated in figure 21. The bistatic synchronized acoustic system in figure 21 comprises two groups 101 and 102 of remotely controlled element arrays communicating with each other through a radio telemetry channel indicated at 105. The group 102 is driven by a number of generators 103 having frequencies indicated at Fl, F2 and FM, whereas the receiver side has corresponding detectors etc. as indicated at 106. In addition to high sensitivity as mentioned, this system gives a large spatial resolving power both with respect to range (depth) and azimuth. In summing up what is explained above on transverse resolution and longitudinal resolution, respectively, reference is made to figure 22. In this figure it is seen that for equal transverse and longitudinal resolution, i.e. Δx / z =1, the following conditions have to be satisfied: F = D Δ F 2R Here F/F is the relative bandwidth, whereas R is the object range and D is the antenna aperture. The upper part of the figure (22A) shows that spaced receiving antennas give transverse 22 resolution, whereas the lower part (figure 22B) shows that spaced frequencies give resolution in depth. The fact that the acoustic sources shall be coherent, derived from a set of phase-locked sinusoidal generators and beam-shaped so as to appear as an array of point generators, makes it important to have transducers with a maximum of efficiency for coherent signals. This can be implemented with "matching circuits" as shown in figure 23. In this figure there is shown a sinus generator 110 energizing a number of resonant transducers 111-112 of high efficiency. Each of these transducers moreover is coupled to a matching unit in a group 113, in which there is effected matching to the acoustic propagation medium, by a form of coupled circuits in analogy to that which is employed in microwave technology. Any desired phase front can thus be produced, in particular in order that the arrangement shall operate like a point source. In this manner a number of phase-controlled resonant transducer cavities will constitute the acoustic sources. This phase front matching by distributed source elements gives a good impedance matching and reduced caviation. A design being advantageous in practice consists therein that the acoustic transducers in the form of tuned resonance elements 111-112 are coupled to an impedance matching network 113 in the form of quarter-wave transformers from an air column to the water mass or volume.

Claims: 23 C L I M S 1. System for marine seismic exploration of structures (2) under the sea bottom, comprising transmitter means (11) with an array of acoustic transducers (12, 21-23) for synchronously and coherently generating an acoustic wave in the water, and receiver means for detecting acoustic waves being reflected from said structures (2) , c h a r a c t e r i z e d i n that the receiver means comprises a radar system (10) adapted to illuminate (KR) an area of the water surface (1) in which the reflected acoustic waves (KA) appear, and to establish therein a number of interferograms each corresponding to the acoustic field from a scattering point in the structure (2) given by the position of the point and the frequency of the field, said receiver means being adapted to determine (15) the degree of coincidence between said interferograms and acoustic perturbation patterns affecting the electromagnetic reflection coefficient of the water surface (1) , thereby to obtain a recording of seismic properties of said structures (2) (figure 1, 2) . 2. System according to claim 1, c h a r a c t e r i z e d i n that the acoustic transducers (21-23, 53, 101, 102) are arranged in a one-dimensional, two-dimensional or three- dimensional array adapted to generate harmonic, multifrequent signals (figures 2, 11, 21) . 3. System according to claim 1, c h a r a c t e r i z e d i n that the acoustic transducers (77b, 77d) are adapted to produce acoustic wave fields being focused at a point beneath the surface so as to increase the energy density in strata of interest in said structures (72) (figure 15, 16) . 4. System according to claim 1, c h a r a c t e r i z e d i n that the operation of the acoustic transducers (12, 35-38) is synchronized with the operation of the radar system (10, 31- 33) (Figure 1,2.) . 24 5. System according to claim 1, c h a r a c t e r i z e d i n that the radar system (10, 31-33) is adapted to operate at a number of frequencies generating a large number of interference patterns with different wave numbers and having plane wave fronts (fig. 1, 2) . 6. System according to claim 5, c h a r a c t e r i z e d i n that the radar system has a number of apertures (61, 62) located side by side for generating wave number vectors of any direction (figure 15) . 7. System according to claim 1, c h a r a c t e r i z e d i n that the illuminated area (39) is wide compared to the wave length of the acoustic interference patterns so as to obtain a focusing effect, and that a number of frequencies (33) are used for matching to acoustic waves having phase centers directly beneath the radar position (figure 3) . 8. System according to claim 7, c h a r a c t e r i z e d i n that the radar system is provided with more than one aperture and operates with time varying frequencies for generating interference patterns corresponding to a focused acoustic wave, and that the phase center of the acoustic wave is adapted to have an arbitrary location in three dimensions with respect to the radar position (figure 17) . 9. System according to any one of claims 5, 6, 7 or 8 , c h a r a c t e r i z e d i n that on the basis of knowledge with respect to the frequency of the acoustic field and by using only one interference pattern, there is provided information about the scattering coefficients of a number of depth or structure strata simultaneously, taking into account the local incident angle and phase velocity of the acoustic waves and the requirements as to matching of wave numbers and velocity (figure 4) . 25 10. System according to claim 9, c h a r a c t e r i z e d i n that phase shifters are applied to all frequency components, thereby to give the interference patterns a migration (propagation) corresponding to the possibly projected velocity of the acoustic fields (figure 1) . 11. System according to claim 1, c h a r a c t e r i z e d i n that it operates in a self-learning mode, in which the radar system estimates the condition of the water surface and selects frequencies of the acoustic interference field and the radar interference field respectively, which gives a maximum of contrast of the acoustic field against the water surface. 12. System according to claim 11, c h a r a c t e r i z e d i n that there are implemented algorithms which enhance the surface or interference patterns (footprints) of the acoustic fields in relation to the background constituted by the water surface, for example by means of Doppler filtering and rigidity filtering (figure 3) . 13. System according to claim 1, c h a r a c t e r i z e d i n that the acoustic transducers are tuned resonant elements (111, 112) with impedance matching networks in the form of quarter-wave transformers (113) from an air column to the water mass or volume (figure 23) . 14. System according to claim 1, c h a r a c t e r i z e d, i n that the acoustic transducers (21-23) are arranged on vessels (24-27) having independent positions, and that there is employed radio navigation of high accuracy in order to compensate for deviations from an ideal position in relation to the other transducers in the transducer array (figure 2) .

Inventors: GJESSING DAG T, HJELMSTAD JENS, TONNING ANDREAS, GJESSING, DAG, T., HJELMSTAD, JENS, TONNING, ANDREAS
Applicants: GJESSING DAG T, HJELMSTAD JENS, TONNING ANDREAS, GJESSING, DAG, T., HJELMSTAD, JENS, TONNING, ANDREAS
Application year: 1987
Main classification: G01V1
IPC codes: G01H9, G01S15, G01V1