SCANNING SYNTHETIC APERTURE SONAR SYSTEM

Abstract

A synthetic aperture sonar (SAS) system utilizes a novel timing and pointing method to illuminate and process data from multiple receive channels over one or more elevation swaths. The system further utilizes cross-track interferometry to improve accuracy of three-dimensional mapping.

Claims

1. A synthetic aperture sonar apparatus comprising: a multi-channel digital receiver comprising an array of transducers; and a processor operatively coupled to the multi-channel digital receiver, the processor configured to receive timeshared data from said array of transducers responsive to reflections from a first elevation swath illuminated with a first pulse rate interval and reflections from a second elevation swath illuminated using a second pulse rate interval, the processor further configured to generate image data based on the received timeshared data.

2. The synthetic aperture sonar apparatus of claim 1 further comprising: a second array of transducers spaced vertically from the first-mentioned array of transducers by a baseline distance B, said processor being further configured to receive additional timeshared data from said second array of transducers responsive to reflections from the first elevation swath illuminated with a first pulse rate interval and reflections from the second elevation swath illuminated using a second pulse rate interval, and to apply interferometric processing to generate 3D topographic image data from the timeshared data received from the first-mentioned and second arrays of transducers.

3. The synthetic aperture sonar apparatus of claim 1 wherein the first elevation swath and the second elevation swath have different elevation beam widths.

4. The synthetic aperture sonar apparatus of claim 1 wherein the first elevation swath and the second elevation swath have the same elevation beam width.

5. The synthetic aperture sonar apparatus of claim 1 wherein said transducer receive array is not longer than 25 cm.

6. The synthetic aperture sonar apparatus of claim 1 wherein said transducer receive array is not wider than 6 cm.

7. The synthetic aperture sonar array of claim 1 wherein said transducer receive array comprises a conformable material.

8. A method of operating a synthetic aperture sonar apparatus comprising the steps of: receiving timeshared data from a transducer array responsive to illumination of a first elevation swath using a first pulse interval rate and a second elevation swath using a second pulse interval rate; and applying interferometric processing to generate topographic image data from the timeshared data received from the array of transducers.

9. The method of claim 8 further comprising the steps of: further receiving timeshared data from a second transducer array across said first and second elevation swaths; and generating 3D topographic image data therefrom.

10. The method of claim 8 wherein the first elevation swath and the second elevation swath have different elevation beam widths.

11. The method of claim 8 wherein the first elevation swath and the second elevation swath have the same elevation beam width.

12. The method of claim 9 wherein the first elevation swath and the second elevation swath have different elevation beam widths.

13. The method of claim 9 wherein the first elevation swath and the second elevation swath have the same elevation beam width.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0020] While the specification concludes with claims particularly pointing out and distinctly claiming particular embodiments of the present invention, various embodiments of the invention can be more readily understood and appreciated by one of ordinary skill in the art from the following descriptions of various embodiments of the invention when read in conjunction with the accompanying drawings in which:

[0021] FIG. 1 is an illustration of a multi-channel digital beamforming process in an exemplary radar configuration;

[0022] FIG. 2 illustrates a novel timing concept for increasing SAS area coverage through the use of 2 pulse rate intervals over two elevation swaths in accordance with the present disclosure;

[0023] FIG. 3 is an illustration of across-track interferometry using a phase difference between two vertically spaced receiver arrays;

[0024] FIG. 4 illustrates a simulated forward-looking interferometer power (left) and height (right) data, showing a long flat surface in front of a moving platform;

[0025] FIG. 5 is a schematic illustration of a complete ScanSAS system disposed within a 15 cm diameter UUV; and

[0026] FIG. 6 is a schematic block diagram of the various exemplary system components.

DETAILED DESCRIPTION OF THE INVENTION

[0027] Certain exemplary embodiments will now be described to provide an overall understanding of the principles of the structure, function, manufacture, and use of the device and methods disclosed herein. One or more examples of these embodiments are illustrated in the accompanying drawings. Those skilled in the art will understand that the devices and methods specifically described herein and illustrated in the accompanying drawings are non-limiting exemplary embodiments and that the scope of the present invention is defined solely by the claims. The features illustrated or described in connection with one exemplary embodiment may be combined with the features of other embodiments. Such modifications and variations are intended to be included within the scope of the present disclosure. Further, in the present disclosure, like-numbered components of the embodiments generally have similar features, and thus within a particular embodiment each feature of each like-numbered component is not necessarily fully elaborated upon. Additionally, to the extent that linear or circular dimensions are used in the description of the disclosed systems, devices, and methods, such dimensions are not intended to limit the types of shapes that can be used in conjunction with such systems, devices, and methods. A person skilled in the art will recognize that an equivalent to such linear and circular dimensions can easily be determined for any geometric shape. Further, to the extent that directional terms like top, bottom, up, or down are used, they are not intended to limit the systems, devices, and methods disclosed herein. A person skilled in the art will recognize that these terms are merely relative to the system and device being discussed and are not universal.

[0028] Turing now to the drawing figures, there is described a novel timing and illumination methodology that forms the basis for the invention.

[0029] The theoretical along track (hereafter referred to as azimuth) resolution of a SAS system, .sub.y, is

[00004]${}_y=\frac{l}{2}$

[0030] While in theory one pays a coverage penalty for a shorter physical aperture, the maximum azimuth resolution improves linearly as well. To overcome this limitation, it is proposed to subdivide the receive array into N sections in the along track and using digital receivers achieve a N improvement in resolution, allowing fine resolution from a relatively large aperture.

[0031] From a 0.25 m aperture, the use of, for example, 8 subchannels (N=8), allows a theoretical maximum azimuth resolution of 1.6 cm. At the same time, the longer array reduces the Doppler spread allowing for larger swath coverage.

SAS Design

[0032] In this section, the practical capabilities of a single-sided 0.25 m aperture SAS is discussed.

[0033] As an example, we start with platform with velocity of 2.5 m/s (4.9 kts), which from a critical sampling perspective requires a pulse repetition frequency (PRF) of 20 Hz. For illustrative purposes, we first look at a 10 m mapping vehicle altitude (h) above the sea floor, although we examine the implications of altitude changes later in this section. From a signal-to-noise ratio (SNR) perspective, we also assume a 10 ms pulse, although the necessary pulse length may be relaxed given the significant array gain.

[0034] With a 50 ms pulse repetition interval and a desire to map larger look angles to maximize swath, our design requires 2 pulses in the water at any given time (i.e. the first range ambiguity is used for the swath window). The available swath then constitutes a range of 60-100 ms, or a slant range of 45-75 m. The equivalent ground range, for h=10 m, is 43.9-74.3 m, or an ACR=0.27 km.sup.2/hr. The other parameters are listed in Table 1 and calculated for several different altitudes.

TABLE-US-00001 TABLE 1 single-swath SAS design for a 0.25 m aperture and vehicle speed v.sub.s = 2.5 m/s. Parameter h = 5 m h = 10 m h = 15 m h = 20 m h = 25 m Slant 45.0-75.0 45.0-75.0 45.0-75.0 45.0-75.0 45.0-75.0 range [m] Ground 44.7-74.8 43.9-74.3 42.4-73.5 40.3-72.3 37.4-70.7 range [m] Ground 30.1 30.5 31.1 32.0 33.3 Swath [m] Angle 83.6-86.2 77.2-82.3 70.5-78.5 63.6-74.5 56.3-70.5 range [deg] Look 84.9 79.7 74.5 69.1 63.4 Angle [deg] Elevation 2.6 5.2 7.9 10.9 14.3 Beamwidth [deg] ACR 0.27 0.27 0.28 0.29 0.30

[0035] The required elevation beamwidth to cover the expected swath is 5.2 deg at h=10 m. However, because there is a guard band of the transmit pulse on either side, this requirement is not particularly tight; further, because there is an approximately linear need to grow the antenna beamwidth with increasing altitude (or to perform some level of adaptive beamforming), one must set a maximum antenna size and degrees of freedom based on the expected altitude range, as well as the vehicle orientation requirements. This is explored below, following the discussion of elevation scanning.

Elevation Scanning to Increase Swath: ScanSAS

[0036] The shorter antenna imposes a higher PRF and as such decreases the available timing for receive, reducing the coverage rate. To improve on the ACR values shown in Table 1 we utilize a unique and novel timing and pointing configuration to illuminate and process one or more additional swaths. Although two elevation swaths are described herein, it should be understood that the technique is applicable to multiple elevation swaths and the limitation is dependent on processing power only.

[0037] FIG. 2 illustrates the present timing scenario for the nominal height of 10 m. The upper pulse sequence represents the nominal burst or PRI listed in Table 1. With a minimum PRI of 50 ms, we can achieve 30 m of slant range resolution. The lower PRI must be higher than the minimum, and sogiven the incidence angle ranges at workwe choose to maintain the number of range ambiguities, allowing us to also acquire data in a closer swath, extending the slant range swath by 15 m, and the ground range swath, at h=10 m, by 15.6 m, a 50% increase in the swath and coverage rate. Across all the altitudes, the ACR thus improves to 0.41-0.49 km.sup.2/hr.

[0038] Nominally the penalty one pays for increasing coverage is in azimuth resolution. Achieving full azimuth resolution with SAS requires that a target be imaged with continuous pulses (at PRF=2v.sub.p/l) as long as the target is in the field of view, i.e. for a vehicle velocity v.sub.p, target range r.sub.0, and antenna beamwidth .sub.bw, the target must be imaged for time t=.sub.bwr.sub.0/v.sub.p. However, by timesharing among the two elevation swaths, we limit the number of pulses on a given target.

[0039] The reference design above provides an example. With a vehicle velocity of 2.5 m/s and a center frequency of 450 Hz, each element (not the full array) in an 8 element array sees a Doppler bandwidth of +/80 Hz. The full aperture length for the far swath (r.sub.0=75 m) can be estimated from

t.sub.a=.sub.bwr/v.sub.p [0040] or 3.2 s of synthesis time (8 m). If only 1.6 s of dwell time are available due to the addition of a second swath, the eventual aperture length is 4 m, yielding a final resolution of

[00005]${}_{\mathrm{az}}=\frac{}{2D}$ [0041] or 0.4 mrad. (Recall that the denominator factor of 2 comes from the SAS formation algorithm). This can be confirmed in the Doppler domain by considering that in 1.6 s we send 32 pulses, which when multiplied by the 8 receivers, yields a Doppler resolution of B.sub.D=160 Hz/256=0.625 Hz. In the azimuth direction, this is found from

[00006]${}_{\mathrm{az}}=\frac{B_D}{2v_p}$ [0042] or 0.4 mrad, again. Notably, at 75 m range, 0.4 mrad yields an azimuth resolution of 3 cm, which is still reasonable and in family with state-of-the-art systems.

Transducer Implications

[0043] In azimuth the implications of our design for the transducer array are relatively straightforward: we require 8 independent, equally spaced channels over a distance of 0.25 m. This design parameter is independent of center frequency. For reference we note that at 450 Hz center frequency, 0.25 m is 75 wavelengths, implying that each of the 8 channels is made of a minimum of 19 subchannels, assuming half-wavelength element spacing.

[0044] In elevation the situation is more complex. Table 2 lists the range of look-angles and elevation beamwidth to maintain the correct illumination on the swath:

TABLE-US-00002 TABLE 2 Look angles and elevation beamwidth for near- and far-swaths as a function of altitude. Parameter h = 5 m h = 10 m h = 15 m h = 20 m h = 25 m Look Angle 84.9 79.7 74.5 69.1 63.4 [deg] - Swath 1 Elevation BW 2.6 5.2 7.9 10.9 14.3 [deg] - Swath 1 Look Angle 82.0 73.8 65.3 55.9 44.9 [deg] - Swath 2 Elevation BW 3.2 6.6 10.5 15.4 22.7 [deg] - Swath 2

[0045] While in theory, one would utilize different elevation beamwidths for the two swaths, in practice it is acceptable to have the same beamwidth, as some amount of guard timing is available due to the transmit event (see FIG. 2).

[0046] If we assume [0047] Finest elevation beamwidth of 3 degrees [0048] Coarsest beamwidth of 15 degrees

[0049] For our same 450 kHz center frequency we arrive at an antenna 6.3 cm at its maximum extent in elevation dimension, or approximately 19 wavelengths (i.e. 38 elements). This dimension does scale with center frequency, so the higher center frequency that is chosen keeps the aperture to a reasonable size. Further, while there is no direct need to sample directly all channels, the ability to control and limit the array would arrive at the steerability required, which is up to 40 deg if the vehicle roll is not considered, or 20 deg if vehicle orientation can be used to make a coarse roll correction.

Transducer Design

[0050] Transducers made with piezocomposite are naturally broadband (up to 80%), have high sensitivity (up to 10 dB more than solid ceramic transducers), and can be made in cost effective large sheets. The material is conformable and can be shaped to tailor the directivity (beam patterns) or to mimic the shape of a vehicle to reduce hydrodynamic drag. Because a significant amount (60% to 85%) of the material is polymer, piezocomposite transducers are much lighter weight and couple better to the water. There is a large design space that allows a wider transducer optimization by adjusting the matrix material, active material (PZT in most cases), and the ratio of active to inactive material (volume fraction). Shading a piezocomposite transducer to reduce side lobes is straightforward and very cost effective. This is accomplished by screen printing a shaped electrode directly on the piezocomposite.

[0051] The transducer design for this system leverages the advantages of piezocomposite construction. Multi-element arrays can be constructed out of sheets of composite. The element configuration, count, and size can all be adjusted based on the desired parameters. Arrays can also be tiled together for longer vehicles, should that be desired. Should the frequency need to be modified, this is done by simply grinding the composite to a different thickness, or using a different ceramic preform.

[0052] Piezocomposite transducers have been used in a wide range of applications, including mine detection, mine classification, swimmer detection, harbor security, torpedo homing, parametric sonar for sub-bottom profiling, AUV sonar including side-scans, forward looking sonar, obstacle avoidance, bottom mapping, and acoustic communications (ACOMMS). They perform well in shallow water, at full ocean depth, low duty and high cycles (including 100%) and are amenable to creating low-cost, high element count arrays.

Interferometric Height Mapping

[0053] While side-looking SAS can provide two-dimensional imagery, an ambiguity still exists between distance and topography. That ambiguity can only be solved through the integration of another dimension of imaging. While stereo imaging is possible, better accuracy is obtained through the use of the received phase from two across-track (or vertically) separated antennas, i.e. cross-track interferometry, or XTI. The XTI concept is illustrated in FIG. 3.

[0054] The improvements by the use of height, versus contrast, for object discrimination is considerable.

[0055] FIG. 4 illustrates a simulated forward-looking interferometer power (left) and height (right) data, including a long flat surface in front of a moving system. A 50 cm obstacle is present at 25 m distance. The obstacle here has lower contrast than the background but is still highly visible in the height image (i.e. the elevated cellyellow color).

[0056] With a high SNR, favorable geometry, and reasonable separation (6 cm, or the elevation transducer size), we can achieve sub-centimetric-class precision in our height mapping providing true 3D mapping. With the present system parameters, the phase wrap height would be 0.5 m which allows for a significant amount of relief across-track (range pixel to range pixel).

Multichannel Digital Transceiver and Processor

[0057] Referring to FIGS. 5 and 6 there is illustrated an exemplary ScanSAS system 100 in accordance with the teachings described herein. The system 100 may be deployed in a UUV housing 900 which measures approximately 15 cm or 6 inches in diameter and may be of indeterminate length, although it must accommodate at least the 25 cm length of the transducer receive array.

[0058] The core multi-channel digital receiver and processing technologies that enable the present system are based on the Remote Sensing Solutions ARENA technology (ARENA is a registered trademark of Remote Sensing Solutions, Inc.). Representative ARENA systems, subsystems and the underlying digital architecture are fully explained and disclosed in US Patent Publication No. US20180321358A1, the entire contents of which are incorporated herein by reference.

[0059] The ARENA digital architecture allows the reception and scanSAS processing of the full 240 channels of the exemplary design. With a full processing throughput of less than 20 MHz, the design created herein is an order of magnitude less than the full ARENA capabilities, and fits comfortably in a smaller ARENA 103 series chassis, which may be re-packaged for the vehicle form factor (discussed below).

[0060] The exemplary system 100 as illustrated comprises the following elements: [0061] A small transmit element 102 and associated drive electronics 104; [0062] A 25 cm6 cm transducer receive array 106; [0063] An optional 25 cm6 cm second transducer receive array 108 for topographic mapping; [0064] A miniature 240 channel digitizer 110, 112 mounted to the back of each array (one per receiver) 106, 108; [0065] A single, miniature receiver and processor card 114, on which all motion compensation and InSAS processing is performed, as well as the system timing.

[0066] The overall system 100 can be extremely compact. A small picture of a comparable RSS ARENA 103 unit is shown in FIG. 5, which fits easily into the vehicle diameter as illustrated.

[0067] A system as illustrated in FIG. 5 can fit within a 100 W target power value (50 W for system operation, and 50 W for InSAS processing), some 30-50% lower in power than other competitive systems, and, with optimization, possibly achieve far lower values.

Collected System Parameters

[0068] All of the reference design system parameters are collected into the following table for reference.

TABLE-US-00003 Parameter Unit Value Platform Altitude Range m 5-25 Platform Velocity Range m/s 1-2.5 System Center Frequency kHz 450 System Bandwidth kHz 80 Receive Transducer Array m 0.25 Length Number of sub-arrays - 8 azimuth Receive Transducer Array m 0.063 Width Elevation Steerable Range deg 20 or 40 Theoretical Resolution cm cm 3 3 Baseline Separation m 0.063 Transmit Transducer m 0.03 Length Transmit Transducer Width m 0.01 Total ACR (single swath) km.sup.2/hr 0.41-0.49

[0069] Then on a per-swath basis, for a 10 m altitude:

TABLE-US-00004 Parameter Unit Swath 1 Swath 2 Nominal Look Angle deg 73.8 79.7 Swath Near Ground m 43.9 28.3 Range Swath Far Ground Range m 74.3 43.9 ACR Km{circumflex over ()}2/hr 0.27 0.13 Topography Precision cm 1.6 0.5 Topography Horz. Posting cm cm 10 10 10 10

[0070] The reference design chosen utilizes a center frequency of 450 kHz, and arrives at a basic mapping receive transducer of 25 cm in azimuth and 6.3 cm in elevation. Adding interferometic mapping doubles the width of the array to 25 cm12.6 cm.

[0071] The length of the aperture is independent of center frequency, and is instead directly related to the Area Coverage Rate (ACR). Any change in the frequency would not impact the length; however, if additional length on the vehicle is available, a longer aperture may provide increased ACR or a simplification of the ScanSAS timing scheme. This also enables this design to be deployed on larger UUVs as well.

[0072] The width of the aperture is directly related to the chosen center frequency. As such, a lower center frequency may provide less loss as a function of distance, but also suffers from lower antenna gain for a given aperture length, even if the aperture can grow in width. As such, reducing the center frequency may not provide significant if any gains, depending on the final area available for the transducer array.

[0073] Having thus described certain particular embodiments of the invention, it is understood that the invention defined by the appended claims is not to be limited by particular details set forth in the above description, as many apparent variations thereof are contemplated. Rather, the invention is limited only be the appended claims, which include within their scope all equivalent devices or methods which operate according to the principles of the invention as described.