Abstract:
In a Doppler acoustic radar system, bistatic transmitting and receiving antennas are employed in a manner to significantly improve wind measurement performance, and to simplify equipment deployment.

Description:
BACKGROUND OF THE INVENTION 
     This invention relates generally to Doppler acoustic echo sounder (or acoustic radar) systems (also called Doppler sodar or Doppler acdar or Doppler echosonde systems), and more particularly concerns the employment of bistatic transmitting and receiving antennas in such manner as to significantly improve wind measurement performance of the system, and to simplify equipment deployment. The term &#34;acoustic radar&#34; as used herein refers to an acoustic echo sounder. 
     A Doppler bistatic acoustic wind monitoring system of the pulse type transmits an acoustic pulse upward along a first beam A (defined by the transmit antenna pattern), and then listens for and analyzes echoes returned from atomspheric scatterers. The listening is accomplished along a beam B defined by the receive antenna pattern. A and B beams must be coincident, at least at the scatter range being investigated; in fact, in a monostatic (MS) system, the beams are everywhere coincident because the same antenna is used first as the transmitter and then as the receiver. A bistatic (BS) system is defined as one in which A and B beams are not everywhere coincident. 
     The Doppler frequency shift of the returned signal is found to be proportional to the component of the velocity of the scatterers along a line bisecting the angle made by the lines extending from the scatterers to the receive antenna and to the transmit antenna. The scatterers are small scale temperature and velocity anomalies that move generally with the local wind. The location and size of the region of the scatterers are determined by the beam geometry (beam width at scatter zone, etc.), by the time interval between transmission of the acoustic pulse and receipt of the echoes, by the pulse length, and by the time constraints in the spectrum analyzer device. Thus, for one set of antennas and beams A and B, it is possible to obtain one velocity component for the local region; for a second antenna pair with different orientation a second velocity component can be found; and a third antenna system will yield a third velocity component. If three components can be found for the same height region, and if the atmosphere can be expected to be horizontally homogeneous, then even if the sensed regions are somewhat displaced horizontally, a representative wind vector V can be derived from the three components (or from two components if we assume the wind is horizontal) using geometrical relationships. 
     Temperature variations within the sensed region cause back scatter of sound generally at all angles toward the transmitted beam (although not at right angles to the beam); velocity variations do likewise, except for angles around 180° (backscatter direction). Thus, if there is turbulence but no temperature microstructure at the scatter region, a monostatic antenna system will not receive echoes, but a bistatic antenna system will. On the other hand, conventional bistatic systems require excessive spacing of the transmit and receive antennas for surveys of higher scatter regions; also, one of the antennas must be of sufficiently wide beam or have several beams in a fan orientation so as to permit surveys of scatter regions throughout a range of heights. 
     There are various complex relations between the scattering sensitivity, beam geometry, range, and ambient noise which enter the evaluation of what beam systems will give the best signal/noise ratio, and hence the best range, in a particular meteorological situation--and, for the deviation of a wind vector, there are further considerations of the accuracy consequences pertaining to the geometric relationships between the observed components and the total wind. See in this regard, &#34;Review of Geophysics and Space Physics&#34; Volume 16, No. 1, February 1978, and &#34;Journal of Geophysical Research&#34;, Volume 79, No. 36, Dec. 20, 1974, pages 5585-5591. These and other considerations demonostrate the need for simplification of the equipment needed to achieve high performance. 
     SUMMARY OF THE INVENTION 
     It is a major object of the invention to provide a bistatic system which offers features of simplicity, accuracy, and high signal/noise ratio. A key feature of the system is that it retains much of the simplicity of a monostatic system yet achieves performance of a bistatic system. Also, very little additional or different software is required for the analysis of the Doppler signals and the construction of the wind components in the desired coordinate system. The relation of the bistatic to the monostatic system is so close that it becomes possible to operate in monostatic mode at low altitudes and bistatic mode at high altitudes, and thus cover a large altitude range while utilizing each mode where it has particular advantages. 
     Basically the system is embodied in a combination that comprises: 
     (a) a first transmitting antenna or transducer T 1  oriented to transmit acoustic radar beams toward an atmospheric locus of wind moving with velocity characterized by vector V, 
     (b) a first receiving antenna or transducer R 1  oriented to receive acoustic radar echo beam signals from said locus for processing to enable detection of a component of the velocity vector V, 
     (c) the antennas or transducers having beam characteristics and orientations such that the center of the beam pattern from one (the beam center path) intersects the center of the beam pattern of the other with relatively small angularity. 
     As will be seen, the reduced angularity is typically less than 25°. With the two beams T 1  and R 1  being relatively narrow (called pencil beams) and with the beam patterns being symmetrical (mirror images), the region of beam coincidence is not wide and the wind component being sensed in the region of greatest sensitivity at each range is in the same direction at each range. In fact, the velocity component sensed most strongly is everywhere along the bisector line, equidistant between the centers of the beams. With conventional bistatic system having non-symmetrical transmit - receive beam configurations, the velocity component at different ranges has different directions. The wind component sensed by a monostatic system is in the same direction at all ranges. The bistatic system which utilizes the monostatic antenna for either transmit or receive can be given a configuration which senses the wind component along approximately the same elevation angle as the monostatic system sensed (although at a different azimuth). This feature lets the data-handling techniques of a monostatic system be used for a bistatic which incorporates the monostatic antennas. 
     The reduced angularity of the beams means that the region of coincidence of the pencil beams can be long. This is a desirable feature for a Doppler acoustic radar. 
     Each sensing plane (defined as the plane containing T 1 , R 1 , and the point of intersection of the two beams) is typically from vertical in a configuration which is sensing a component of the horizontal wind and which has the sensed component directions at approximately the same elevation angle of the beam from either T 1  or R 1 . The sensing plane can be vertical for a bistatic system. 
     Additional objects include the provision of a first set of transmit and receive antennas T 1  and R 1 , and a second set of such antennas T 2  and R 2 , each set oriented to define beam center paths with relatively small angularity therebetween, and the antennas T 2  and R 2  being azimuthally offset from T 1  and R 1  in such manner that wind velocity vector components may be determined to enable determination of a true wind velocity vector V, in a horizontal plane. 
     These and other objects and advantages of the invention, as well as the details of an illustrative embodiment, will be more fully understood from the following description and drawings, in which: 
    
    
     DRAWING DESCRIPTION 
     FIG. 1 is a frontal elevation showing transmit and receive antennas with a perpendicular bisector plane; 
     FIG. 2 is a side elevation taken on lines 2--2 of FIG. 1; 
     FIG. 3 is a plan view taken on lines 3--3 of FIG. 2; 
     FIG. 4 is a representation of derivation of a wind vector V; 
     FIG. 5 is a plan view of a system embodying the invention; and employing two sets of transmit and receive antennas; 
     FIG. 6 is like FIG. 5 but showing three sets of T and R antennas; 
     FIG. 7 is a side elevation showing the antenna; 
     FIG. 8 is a waveform diagram; and 
     FIG. 9 is a side view of a bistatic configuration for an antenna pair with the sensing plane vertical. 
    
    
     DETAILED DESCRIPTION 
     1. Technical Background 
     The velocity of the wind can be measured by transmitting sound into the atmosphere and then listening for echoes generated by the scattering of the sound by turbulence anomalies. The frequency of the detected echoes may be different from that of the transmitted signal if the turbulence field is moving in relation to the transmitter or receiver. Since the turbulence is carried along with the atmosphere, this movement is equivalent to the wind velocity. 
     The relation between the change in frequency and the geometry of the transmitting, receiving, and velocity directions can be expressed as: 
     
         V=cΔf/2f cos φ                                   (1) 
    
     where V is the velocity component along the bisector line (the line midway between the transmit and receive paths), f is the frequency of the transmitted signal, Δf is the Doppler shift, c is the speed of sound in the medium and 2φ is the angle between the transmit and receive paths. 
     The equation shows that the frequency shift is proportional to the magnitude of the velocity vector component lying along the angle bisector of the transmit and receive directions. 
     The operation of a Doppler acoustic radar is dependent upon receiving an echo with an amplitude sufficiently greater than the background noise to permit the accurate determination of its frequency content. The amplitude of the echo can be expressed as the received power, P r , which is given by the following equation: 
     
         P.sub.r ασ/R.sup.2,                            (2) 
    
     for a system operating at a given transmitted power level, a given transmitted pulse length, a given antenna size and efficiency, and under fixed conditions of atmospheric attenuation. In this equation, σ is the portion of the transmitted power scattered back into a unit solid angle from a unit volume, and R is the distance from the transmitting and receiving antennas to the scattering medium. 
     The scattering coefficient, σ, is related to the level of atmospheric turbulence and the geometry of the acoustic radar system in accordance with the expression: ##EQU1## where C v   2  /C 2  is the normalized velocity structure function, C T   2  /T 2  is the normalized temperature structure function, and φ is the scattering angle, the angle between the transmitted and received directions. T is temperature. C is the speed of sound. 
     The accuracy of the velocity component measurement is proportional to the accuracy of the measurements of the frequency shift, and the proportionality constant is a function of the geometry. In the presence of ambient noise the accuracy of the frequency measurement is related to the signal amplitude which is also a function of the geometry. It is clear that the effectiveness of a Doppler acoustic radar wind measuring system is closely tied to geometrical configuration of the transmitting and receiving antenna beams. 
     2. System Description 
     In FIGS. 1-3, transmit and receive antennas T 1  and R 1  are shown at ground level and oriented so that T 1  transmits acoustic radar beams, with beam centerline designated by A, toward an atmospheric locus L moving with velocity characterized by vector V. The latter is seen to have components V x  and V y , in FIG. 3. The receiving antenna R 1  is oriented with its centerline designated by B, intersecting A at point M. R 1  receives acoustic radar echo beam signals from locus L, for processing to enable detection of one component of V. As shown in FIGS. 3 and 4, T 1  and R 1  are used to determine component V x  which extends along a bisector line C. L is spaced on or near line C and need not coincide with M. (Another set of antennas T 2  and R 2  is used to determine V v , as will appear). The term &#34;antenna&#34; as used herein refers to a transducer. 
     The antennas T 1  and R 1  are oriented to define beam center paths A and B with relatively small angularity α therebetween. Typically, α is less than about 25°, for best results; thus, under these conditions, there is substantial range over which there is overlap between the &#34;pencils&#34; defined by the beams transmitted from T 1  and received by R 1 . The centers of such pencils are defined by A and B. Line C in FIGS. 1 and 3 bisects the angle α. A &#34;pencil&#34; exists when, for example, the beam cone transmitted by T 1  defines a cone angle ψ less than about 15° between the half power points. 
     It will also be noted that the sensing plane, which includes beams A and B and the bisector line C, is tilted from vertical, in order that V x  may be detected (see FIG. 2). For best results, in measuring a horizontal wind component, the tilt angularity β is greater than about 10°, from the vertical, is less than 45°, and may typically be about 30°. 
     FIG. 5 shows a bistatic system wherein antennas T 1  and R 1  are employed, as described above. A similar set of antennas T 2  and R 2  is also employed, to derive V y . T 1  and T 2  are shown as at the same location, and R 2  is azimuthally offset from R 1 . In FIG. 5, R 2  is offset 90°, azimuthally, from R 1 . 
     In a typical example the beams A 1  from antenna T 1  points upwards 60° at an azimuth of 000°. The B 1  beam antenna R 1  is moved to a location 210 m distant from T 1  at an azimuth of 75° and pointed up 60° in a direction so the beam center intersects the A beam center at 700 m height point M 1 . The component observed, approximately along the bisector line C 1 , is at an elevation angle of 60.9° and an azimuth of 345°. The centerlines of the beams intersect at an angle α which is 15°. If the pencil beams are about 10° wide between their associated half power points, there is a long range of some coincidence of the beams A 1  and B 1 , from which echoes can be expected. To complete the two-component system, beam A 2  from antenna T 2  points up 60° at an azimuth of 90°, and the B 2  beam antenna R 2  is set out 210 m at an azimuth of 165°, pointed up 60° to intersect the A 2  beam center a 700 m height point M 2 . The velocity component observed with the A 1  - T 1 , antenna pair is at an elevation angle of 60.9° and an azimuth of 345°, while that from the A 2  - T 2  pair has an azimuth of 75°. These two components are orthogonal in the horizontal plane, and together define V. T 2  may be offset from T 1  so long as the line between T 2  and R 2  is perpendicular to the line between T 1  and R 1 , so long as the range of R 2  from T 2  and of R 1  from T 1  are the same, and so long as winds at L 1  and L 2  are, as a practical matter traveling at the same velocities V. 
     If the same data reduction algorithms are used for deriving wind components from the echoes of this new system as for deriving them from the old MS system, the wind will be almost correctly ascertained, except for a constant 15° azimuth correction which is to be subtracted from every azimuth reading. There are several additional small errors, and corrections can be made for them; but they are so small that they are generally insignificant in comparison to other errors inherent in acoustic Doppler systems and can be ignored. 
     One such small error is due to the fact the bisector line C 1  (and also line C 2 ) is at an elevation angle at 60.86° rather than the 60.00° of the original A-B beams. If the wind geometry algorithm with 60° is used for the bistatic mode, the wind component will be indicated to be 1% weaker than actual. A logical correction technique would be to tilt all antennas 59.6° instead at 60°, and use the same 60° algorithm for data reduction in both the MS and BS modes; each would be then be in error by only about 1/2%. 
     Another small error arises because the Doppler shift obtained with bistatic antennas is less as the observed scattering angle increases. For the geometry described herein, at altitudes exceeding 700 m (scatter angles less than 15°) the decrease in Doppler shift is less than 1% and so the indicated velocity component is low by less than 1%. Even at a height of 350 m the effect is under 3.5%. 
     Another error lies in the fact that the length of the path from T 1  to the scatterer locus L to R 1  is altered when T 1  and R 1  are separated in the bistatic mode as contrasted with the situation when they are co-located in the monostatic mode; this effect, plus the difference in elevation angle for the acoustic beams between the modes, means that a time gate will represent a slightly different height for the two modes. There is no differences at a height of 700 m, and differences are only about 9 m at heights of 350 m and 1400 m. 
     3. Variations 
     For any system the role of the transmitters and receivers can be reversed without altering the principles involved. 
     The particular example described above is for a bistatic system with antennas T 1  and R 1  (and T 2  and R 2 ) separated by 210 m, set at an elevation angle of 60°, and pointed so the beam centerlines intersect at an angle of 15° at a height of 700 m. If the T-R separation distance &#34;d&#34; is decreased, the pencil beams will overlap more over a longer range of distances (altitudes), which is desirable, but the scattering angles for each height will be less than for the example, which decreases the amount of echo received from the scatterer. For instance, a lesser antenna separation will require the azimuth of antenna R 1  to be altered (to greater than 75° but less than 90°), and the azimuth of the bisector line C will be closer to the azimuth of antenna T 1 . Of course, the system can be adapted to shorter ranges by preserving all angles but decreasing all dimensions proportionally. 
     Accordingly, there are various compromise configurations which are practical and meet the main objective of having a bistatic system which invokes similar geometry and uses similar algorithms to those used in the monostatic system. For example, receiver R 1  can be put 183.3 m out at an azimuth of 83.5°, pointing up to intersect beam A 1  at 700 m and an angle of 13°. The bisector C 1  turns out to be at an elevation angle of 60°, the same as beam A 1  (and azimuth of 353.5°). Beam B will be at an elevation angle very close to 60°, and the range from R 2  to the intersection point M or the locus L will be very close to that from A to the intersection point or locus. Various other configurations can be used, with the elevation angles of T 1 , R 1  and the bisector C 1  being so similar that in a practical application the differences can be ignored, and with the acoustic leg lengths from T 1  and R 1  being so nearly similar that the differences have no practical consequence. 
     Ordinarily, the system will be operated in both MS and BS modes. For example, during a five-minute period, operate in a sequence giving fast repetitions of the monostatic pulses to handle low altitudes, and slower repititions of the bistatic pulses to handle the high altitude. See FIG. 4 in this regard, using computer 10. Print out winds in 33 m height gates from 67 m to 400 m from the monostatic mode, and winds from 400 to 1500 m in the bistatic mode. The redundant measurement at 400 m serves as a consistency check. Representative pulse frequencies are 1,000 to 4,000 Hertz; repetition rates are 1 second to 10 seconds; and pulse durations are 20 milliseconds to 400 milliseconds. 
     The discussion has assumed each beam to be a narrow, pencil-beam, symmetrical about its centerline, with a width of about 10° between the 3 db or 6 db power points. The beam can be made narrower or wider or be asymmetrical about its own centerline. 
     A three-component monostatic system can be made by using three beams at elevation angles of 60° and azimuths of 000, 120, and 240, from co-located antennas. See T 1  and R 1  ; T 2  and R 2  ; and T 3  and R 3  in FIG. 6. This can be merged into a combined monostatic-bistatic system by adding three more pencil-beam transmitters at 210 m distances and azimuths of 75, 195 and 315°. These geometries constitute a symmetrical system. Asymmetrical variations can also be employed. For example, to the basic-two-component system discussed herein, add another antenna pair handling the W sector (T 3  pointing at an azimuth of 270°; R 3  being at an azimuth of 345° at a range of 210 m). 
     It should be noted that with a three-component system, the elevation angles can be increased (say, to 70°) while still retaining adequate accuracy for horizontal wind observations. With different elevation angles, the azimuth orientations of the extra bistatic transmitters must be altered somewhat. 
     As noted, the remote antennas can be either transmitting or receiving. If they are receiving, they do not require much power--and in fact could be made battery powered with telemetering of the received echo signals. 
     In FIG. 5, item 11 represents circuitry to transmit pulses to T 1  and T 2  and to receive pulses from R 1  and R 2  to reduce the data. See for example the circuitry in U.S. Pat. No. 3,675,191 to McAllister. Bus connections appear at 16-18. 
     Also in FIG. 5, in MS mode, T 1  represents both a transmitter and receiver. 
     In the above, and generally speaking, C 1  is perpendicular to the lines connecting T 1  and R 1  , but C 1  can in some cases be tilted from such perpendicularity. FIG. 9 shows such a case where the sensing plane is vertical and C a , the bisector line for transmitter T a  and receiver R a  is not vertical. FIG. 9 can also represent the case where the ground is tilted and C a  is vertical when viewed frontally in the sensing plane, but the sensing plane is tilted from the vertical. All these bistatic configurations can be used to augment monostatic acoustic radar of a type producing facsimile charts, the bistatic mode augmenting echoes, over a long range of coincidence. 
     FIG. 7 diagrammatically shows a section through a transmitter, with housing 13, antenna dish 14, and sound wave emitter 15.