Method and apparatus for measuring water velocities and discharges in rivers and the like

Abstract: An improved method of measuring the flow velocity of water bodies, characterised in that, in addition to a principal direct path along which acoustic signals are propagated across and through the direction of flow of water, an auxiliary indirect path for acoustic signals is utilised comprising two segments of an acoustic signal reflected at a water-air or water-ice boundary of said body of water.

Description: BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates generally to flowmeters, and particularly to acoustic flowmeters. More particularly still, it relates to acoustic flowmeters based on the measurement of propagation times across the whole river, and which are suitable for measuring the average flow velocity in rivers, or the like bodies of water, in which the water level varies, and also under ice. 2. Prior Art of the Invention Supersonic flow meters are being used successfully to determine the flow velocity of water volume in rivers and canals. The mean velocity of water in the reference plane of the acoustic path is determined by measuring the propagation times of acoustic pulses in opposite directions. In order to determine the mean velocity for the entire flow cross-section, a correction factor (k) must be applied. k depends on the position of the acoustic transducers, the flow cross-section profile, the relative heights of the acoustic path, and the bottom roughness. It is known that if the acoustic path is located about 40% above the bottom, then the correction factor becomes independent of the bottom roughness and equal to unity. 40% above bottom is therefore a preferred depth of installation. However, as the water level rises, and the relative installation level above ground falls to 20% or less, the effects of the bottom contour become more significant and the uncertainty in the value of the proper correction factor k increases. The latter may be estimated by computation using uncertain assumptions about bottom contour variations upstream as well as about river width, which is not satisfactory. The problem is further aggravated in winter, when k is also influenced by the surface contour of the ice cover, and more uncertain assumptions about the same must be made in order to determine the new k. Errors in flow velocity determination may, therefore, reach or exceed 40% in winter, because the logarithmic vertical velocity distributions, which are stable in the summer with good accuracy, are no longer stable in the winter under ice. In order to mitigate the above mentioned deficiencies, it has been suggested that acoustic flowmeters be stacked vertically to measure at several planes, thereby apprehending the different vertical distribution of flow under ice. Such an arrangement is also likely to improve accuracy in the summer. Practical experience, however, shows that the complexity of such installation is significant, and that also as a result the susceptibility to failure is increased, the complexity of apparatus increases proportionately with the number of planes, while the improvement in accuracy remains marginal. A further problem with multi-level installations is that a top level installation close to the ice layer (for an accurate assessment of the vertical velocity distribution) is incompatible with other theoretical requirements, for in order to avoid multi-path echoes, it is necessary to measure at minimum depth, which grows as the acoustic path length does. Thus, the top installation is constrained downwardly from above, while the bottom installation is constrained upwardly from below. These constraints sometimes cause both top and bottom installations to coincide, which means there is no room for multi-level measurement. A novel solution is therefore indicated. The use of acoustic energy to measure flow velocity is the subject of three United States patents of general interest. Unites States patent 3,633,415 granted January 11, 1972 to Lu is entitled Flowmeter. It teaches how flow velocity is measured in any direction in the plane by placing transducers at the corners of a triangle in the plane and measuring time differences of travel of acoustic energy between transducers. United States patent 3,861,211 granted January 21, 1975 to Dewan is entitled Ultra-Low Flow Velocity Current Meter. This invention relates to the use of an acoustic signal to determine current velocity of water. The device makes use of a centrally single positioned transducer with four equispaced receivers equidistant from the transducer. A carrier frequency with two different modulator frequencies is directed to the transducers and the receivers receive the transmitted pulses through the water current. The velocity of the water is determined by measuring the phase angle of the different signals that pass through the water and which are received by the receivers. The different signals permits one to determine the velocity of the water. In United States patent 3,949,605 granted April 13, 1975 to Stallworth et al., both acoustic and electrical signals are used to measure flow in large bodies of water by comparing the delay between acoustic and electrical propagation. SUMMARY OF THE INVENTION The present invention provides a solution to the problem of the prior art, which avoids multi-level installations by operating a single level installation in two modes. In the first, principal mode, the principal acoustic signal path is direct and crosses the flow direction at an angle of approximately between 5 and 60 , with the acoustic transducers place about 20% to 40% above the bottom. In addition to the principal, direct, path, an auxiliary composite path is also utilized. The auxiliary measurement path uses the same transducers as are used for the principal path, but is composed of the two straight paths comprising the total internal reflection path of the acoustic signal at the water-air or water-ice boundary. Thus, according to the present invention, a broad aspect thereof is the improvement of the method of measuring the flow velocity of water bodies characterised in that, in addition to a principal direct path along which acoustic signals are propagated across and through the direction of flow of water, an auxiliary indirect path for acoustic signals is utilised comprising two segments of an acoustic signal reflected at a water-air or water-ice boundary of said body of water. While it is of course possible to utilise two separate acoustic transducers for the auxiliary path, this is generally not necessary, since the acoustic transducers are not too narrowly directional such that the use of an auxiliary pair becomes necessary. The use of only one pair of transducers in fact has distinct advantages: first, it is economically the optimum solution, but more importantly, it allows prediction of the arrival time of the' reflected signal with better precision, because only one distance between transducer, the path length, is involved. The velocities obtained from the propagation times along the principal and the auxiliary paths are divided and the obtained ratio permits determination of the equivalent bottom roughness in the summer and, assuming the letter remains constant in winter, the equivalent roughness of the ice-water boundary in the winter. The installation depth of 40% corresponds often to 20% or less at high water level or 50% at low water level. 50% is the position at which an acoustic flowmeter could work at minimum water level. The acoustic energy components reflected by the surface and by the bottom then have the same delay and the delay must be at least 1.5 wavelengths in order to allow discrimination from the direct signal. It is then not possible to use the auxiliary path as described, because the signals of the path reflected by the surface and those reflected by the bottom arrive at the same time and cannot be distinguished from each other. Fortunately, this does not represent a significant limitation to the utilization of the present invention, since in every body of water there are high water phases during which the low water level is significantly exceeded, often several times per year. It is then easy to conduct a measurement via the auxiliary path, but moreover such measurement is particularly desirable, first because at such high water phases it is highly probable that the river flow roughness is significantly altered by the increased water flow velocity, and second, because the k-factor deviates significantly from unity. In the present invention the preferred transducer depth is between 10% and 40% average above the bottom, because then the pulses through the auxiliary path reflected by the surface are delayed sufficiently as compared with the direct principal path pulses and those reflected at the bottom. Therefore, three pulses arrive in sequence at the receiving transducer. The first pulse coming through the direct path, the second (sometimes more than one appearing as a group reflected from different boulders) coming through the reflection off the bottom contour, and the third and last (single) pulse coming via the auxiliary path reflected by the water-air or water-ice boundary. The three pulses are influenced differently by the flow velocity, because all use different parts of the vertical distribution profile in the river cross-section. Therefore. all three pulses carry different information, and it is an object of the invention to extract this information and decipher the unknown parameters which govern the flow in the river. Changes of the velocity distribution near the surface are sensed better through the auxiliary path than through the direct path. It is for this reason that the ratio of the flow velocity measured by means of the auxiliary path to that measured with the direct path is a direct measure of both the bottom roughness in summer as well as the ice roughness in winter. In naturally flowing streams, also in those where free flow is obstructed by water backup, the width is often much larger than the depth. In such cases the friction of the banks can be ignored as compared to the friction of the bottom. BRIEF DESCRIPTION OF THE DRAWINGS The preferred embodiment of the present invention will now be described in detail in conjunction with the annexed drawings, in which: Figure 1 is a plan view of the arrangement according to the present invention for measuring the flow velocity of river water; Figure 2 is a front elevation of the arrangement shown in Figure 1 as seen from the river bank; Figure 3 shows typical vertical distribution profiles of the velocity in the summer; Figure 4 shows similar conditions but with ice cover; Figure 5 is a schematic diagram of the apparatus of the present invention; Figure 6 illustrates the observed propagation times relative to each other; and Figure 7a and 7b shows a high-level flowchart for controlling the apparatus of Figure 5. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 1. Theoretical Considerations It was found empirically that not only in rough pipes and artificial conduits, but also in natural streams at any place along the vertical, the logarithmic law of velocity distribution (1) holds. u(y) = v* (A + 5.75 1g y/ks') u(y) = flow at the height y v* = average shear stress velocity y = vertical distance from the bottom ks' = equivalent sand roughness in cm A = Roughness constant according to Nikoradse = 8.50 The ratio y/ks' can be interpreted as admittance, which is decreased in winter by "connecting" the admittance of the ice friction parallel: u(y) = v* (A + 5.75 1g y ks' + keis y h - y ) If the water depth h and the height of the acoustic path y1 between the acoustic transducers is about the same over the width of the river, then the two velocities-along the direct path and along the reflected path - can be calculated. Along the direct path according to equation (2): u1(y1) = v* (A + 5.75 1g y1 ks' + keis y1 h - y1 ) with y1 = average height of the acoustic direct path above the bottom. Along the reflected path the velocity must be found by integration starting from the height of the acoustic transducer up to the surface: Forming the ratio u1/u2 will eliminate the still unknown shear stress velocity v* and it remains equation (5) with the measurable ratio u1/u2 on the left side and a term with the bottom roughness and the ice roughness on the right side. In the summer, keis is known to be 0 and the only variable left is ks'. In the winter, a time with low velocities, the bottom roughness is supposed to be the same as in the summer and the summer value is still valid. Then keis is the only variable. Equation (5) can be solved by trial and error, most straightforward by applying a procedure known as the Newton approximation method and by handling the integral as a finite sum. Ks' in the summer and keis in the winter will be varied by the microprocessor in the flowmeter according to the invention until it fits the measured ratio u1/u2. The given time when to keep keis = 0 or when to keep ks' constant can be determined by the acoustic flowmeter. The velocity of sound C is found by dividing the known (constant) path length by the measured propagation time. In fresh water the velocity of sound is also a reliable indicator of the temperature. If the temperature drops below 0 degrees centigrade with no major changes of the velocity (which means the bottom roughness remains stable) it is most likely that ice buildup takes place. A similar, but different pattern can be identified in the spring, when the ice will break up. As soon as the equivalent average bottom roughness and the ice roughness in the winter is known, the shear stress velocity can be derived using equation (3). The average velocity in the whole cross-section is found by using the following integral: Vm(ks,keis,h) = v* h (A + 5.75 1g y ks' + keis y h - y )dy The integral is again solved pragmatically by numerically adding a finite sum. The wanted discharge is finally calculated by applying the well-known equation (7): Q = Vm A(h) with Q = discharge [m 3 /s] Vm = calculated average velocity according to elu.(6) A(h) = cross-section area as a function of water level h The acoustic flowmeter determines the average velocity V of the water and the velocity of sound along the acoustic path by solving the following equations V = ( 1 T12 - 1 T21 ) Sp 2cos@ + dC dt t C = ( 1 T12 + 1 T21 ) Sp 2 + dV dt t with T12 = propagation time transducer 1 to transducer 2 T21 = propagation time transducer 2 to transducer 1 Sp = programmed pathlength which must be identical with the real pathlength @ = angle of acoustic path to the mainstream dC/dt = change of velocity of sound with time dV/dt = change of velocity of sound with time t = time between transmitting in either direction As equation (5) shows there is a mutual interference between C and V, if one or both change drastically. Because C is in the order of 1000 times larger than V a stochastic change of C will cause noise in the measured V value. To reduce this error the pulses should be transmitted at the same time. Unfortunately, this would require two separate transmitters and receivers. As described later, the proposed flowmeter will use a time offset between the transmitter pulses for each direction which is small compared to the fluctuation of the water surface but large enough to discriminate the signal coming from transducer 1 from those coming from transducer 2 but using the same digitally gain controlled receiver. For the composite path similar equations hold, but with slightly higher propagation times due to the longer path length: Vr = ( 1 T12r - 1 T21r ) Spr 2cos@ Cr = ( 1 T12r + 1 T21r ) Spr 2 with T12r = propagation time transd. 1-surface-transd. 2 T21r = propagation time transd. 2-surface-transd. 1 Spr = pathlength transd.1-surface-transd.2 @ = angle of acoustic path to mainstream The relative water level above the transducers can be calculated using the measured propagation times: Wre1 = H1 + H2 2 mit with H1 = P 2 T12 T12r 2 - T12 2 H2 = P 2 T21 T21r 2 - T21 2 The error increases with larger pathlength because ray bending due to temperature layers becomes more serious. In this case separate measurement with an upward looking vertical transducer will provide the relative water level Wrel more reliable. The auxiliary composite pathlength is calculated with equation (9) and the known direct pathlength Sp: Spr = Sp 2 + 4Wre1 2 The auxiliary path length is required to solve equations (7) and (8) and to determine the effective cross-section area A by interpolating a stored array with the cross-section A(h) as function of the water level h, A(h) = f(h). The purpose of this invention is to create an acoustic flowmeter which can handle the above mentioned propagation time measurements with high accuracy and which can solve the complete set of above mentioned equations. The average velocities along the direct path and along the auxiliary path are measured sequentially: first the direct pulse in one direction and a delayed time later, while the first pulse is still under way to the opposite transducer, the opposite transducer is fired in the opposite direction. Then, second, the reflected pulse is measured by setting the time window that late that the direct pulse and the bottom reflected pulse will be blocked and only the pulse reflected at the surface will be detected. The signals running the direct way has been observed to be very stable in terms of propagation time and amplitude, only affected by the velocity of sound, velocity of the water and the attenuation of the medium. The signals of the auxiliary path are additionally affected by the conditions of the reflecting surface; the surface can be smooth, rough or even swell modulated, depending on the direction and the strength of the wind. This will change the acoustic conditions from second to second, even changing the position, when the reflection occurs. The propagation time as well as the reflection coefficient will be modulated by waves. It is most important that both pulses, running along the auxiliary path in opposite directions, will find almost identical conditions at the reflecting boundary to avoid jumping from one reflection point to another. The probability for such dramatic event obviously will diminish with increasing simultaneousness of both pulse reflections. This definitely calls for the time coincidence of both transducers to be fired. Another intention of this invention is to solve the problem in the most economical way. Transmission at the same time would require to receive the opposite signals at the same time. Therefore two equal receivers would be needed. This is not only expensive but additional problems will arise. Bandlimited receivers have a certain group propagation time. Both receivers must have the same time delays to meet the accuracy requirements of the acoustic flowmeter and with two receivers additional calibration would be necessary. The invention solves the problem by firing the two transducers shortly one after another with a digital preprogrammed time delay. The delay is long enough to distinguish both groups of signals running in opposite directions, even when received only with one receiver, but short enough to provide simultaneous conditions for the reflection. After having received the first group coming from one direction the receiver is switched with a wideband switch to the opposite direction. This switch does not introduce any phase shift between the two directions and the group propagation time is kept within the tolerance of a few nanoseconds. 2. Structure and Operation When referring to the drawing figures in the following description, the designation numerals 1 - 54a,b are defined as follows: 1 Transducer 1. 2 Transducer 2. 3 Direct path. 4 Path with bottom reflection. 5 Auxiliary path, upper boundary reflection. 6 Acoustic Flowmeter. 7 Cable to transducer 1. 8 Cable to transducer 2. 9 Vertical velocity distribution, smooth bottom, deep water. 10 Vertical velocity distribution, rough bottom contour, shallow water. 11 Velocity u1, measured at height y1. 12 Average velocity u'2 measured along the auxiliary path, smooth bottom. 13 Average velocity u''2 measured along the auxiliary path rough bottom. 14 Vertical velocity distribution, smooth bottom, shallow water, ice cover. 15 Vertical velocity distribution, rough bottom, shallow water, ice cover. 16 Average velocity along the auxiliary path, rough bottom, shallow water, ice cover. 17 Average velocity along the auxiliary path, smooth bottom, deep water, ice cover. 18 Oscillator. 19 Timer Unit: 16-Bit free-running counter. 20 16-Bit Comparator. 21 16-Bit Output Compare Register. 22 Microprocessor. 23 Interrupt lines. 24 Data-Bus. 25 8-Bit Register, 2 Bit Direction Bits, 6 Bits Gain word. 26 controlled direction switch. to fire either transmitter 1 or transmitter 2. 27 Control input of direction switch 26 28 Controlled direction switch. to connect receiver either to transducer 2 or transducer 1 29 Control input of direction switch 28 30 Transmitter TX1. 31 Transmitter TX2. 32 Digitally gain controlled amplifier. 33 Pulse shaping network to convert analog signal pulse to digital pulse. The designations for signals sent by transducer 2 are the same as the following designations for transducer 1, but with a "b" subscript instead of an "a" subscript. 34a Shutter enable D-Flipflop with processor controlled Set-input and signal-controlled Clock-input. 34b Shutter enable D-Flipflop with processor controlled Set-input and signal-controlled Clock-input. 35 Hardware Output OC2 from Compare Register to trigger either transmitter TX1, setting D-Flipflop 34a and opening shutter 39a or transmitter TX2, setting D-flipflop 34b opening the shutter 39b. 36 Hardware Output OC3 from Compare Register to enable the D-input of the D-flipflop 34a. The incoming signal is clocking the flipflop thus closing shutter 39a. 37 Hardware Output OC4 from Compare Register to enable the D-input of the D-flipflop 34b. The incoming signal is clocking the flipflop thus closing shutter 39b. 38 HF-clock (about 100MHZ). 39a Shutter gate, controlled by D-Flipflop 34a. 39b Shutter gate, controlled by D-flipflop 34b. 40a High resolution timer, 32 bit wide. 40b High resolution timer, 32 bit wide. 41a Transmission pulse 1, fired by transducer 1, undelayed. 41b Transmission pulse 2, fired by transducer 2, delayed by tau. 42a Direct signal, sent from transducer 1 appearing at transducer 2. 42b Direct signal, sent delayed from transducer 2 received from transducer 1. 43a Bottom reflected signals appearing at transducer 2. 43b Bottom reflected signals appearing at transducer 1. 44a Surface reflected signal using the auxiliary path and appearing at transducer 2. 44b Surface reflected signal using the auxiliary path and appearing at transducer 1. 45 Direction bit signal, switching after transmission of transducer 1 to prepare for transmission of transducer 2 and reception from transducer 2. 46 Direction bit signal, switching after reception of all signals sent by transducer 1 to prepare reception of signals coming from transducer 2. 47a Time window set by OC3, calculated by Compare/Capture Unit, to enable D-flipflop 34a to clock, when first direct signal 47a arrives, closing shutter gate 39a. 47b Time window set by OC4, calculated by Compare/Capture Unit, to enable D-flipflop 34b to clock, when second direct signal 47b arrives, closing shutter gate 39b. 48a Time window reset by clearing interrupt OC3. 48b Time window reset by clearing interrupt OC4. 49a Open shutter gate 39a by setting D-flipflop 34a when transmitter TX1 fires. 49b Open shutter gate 39b by setting D-flipflop 34b when transmitter TX2 fires. 50a Close shutter gate 39a by clocking D-flipflop 34a with signal 42a. 50b Close shutter gate 39b by clocking D-flipflop 34b with signal 42b. 51a Counting pulses to determine propagation time T12. 51b Counting pulses to determine propagation time T21. 52a Equivalent to 49a. 52b Equivalent to 49b. 53a This time the window 47a, 48a is delayed (not shown) to catch signal 44a and the shutter gate closes at 53a. 53b This time the window 47b, 48b is delayed (not shown) to catch signal 44b and the shutter gate closes at 53b. 54a Counting pulses to determine propagation time T12r. 54b Counting pulses to determine propagation time T21r. Referring now to figures 1 and 2, the line connecting transducers 1 and 2 represents the acoustic direct path 3 , the path with the signal reflected at the bottom 4 , and the path with the signal reflected at the surface 5 ; all of which coincide in the plan-view of figure 1. Both the transducers 1 and 2 are connected to acoustic flowmeter apparatus 6 via cables 7 and 8 . In figure 3, the curve 9 represents the distribution when the ratio h/ks' is large (deep water and smooth bottom), while curve 10 shows the distribution when this ratio is small (shallow water and rough bottom). The auxiliary composite path runs from depth h1, which is the transducer position, to the surface h and down to the second transducer again at height h1. The acoustic pulse therefore is influenced by different amplitudes of flow, in case of smooth bottom according to curve 9 and in case of rough bottom according to curve 10 all along the way from h1 to h and back to h1. 11 is the average velocity for a given vertical, the cross-over point of at 40% of the height h, and 12 is the average velocity between h1 and h for smooth bottom and 13 is the average velocity between h1 and h for rough bottom. Figure 4 shows conditions similar to those in figure 3, but with ice cover. The acoustic flowmeter measures now either by calculation according to equation (9) or by direct vertical measurement the clear depth under the ice, which is the effective depth for the flow. Profiles 4 and 15 are the vertical distributions of velocity with different roughness of the bottom. The average velocities measured with the auxiliary path are the velocities u2' 16 and u2'' 17 . With reference to Figures 5, 6, 7a and 7b the apparatus and its operation will now be described. The controlling unit is a Motorola HC16 microprocessor 22, which comprises the following integral and auxiliary components: (i) A highly stable, phaselock loop type master clock 18 (ii) A 16-Bit free-running counter 19 (iii) A 16-Bit Comparator 20 (iv) A 16-Bit Output Compare Register 21 (v) A Processor Core 22 ; (vi) Interrupt-lines 23 from the timer unit 19 , 20 to the processor core 22 ; The processor is complemented with the following parts: A Data-Bus 24 , connecting the internal and external registers 21 , 31 , 40 and the memory 41 with the processor 22 ; A 8-Bit D-flipflop type register 25 , providing 2 Bits to disable both transmitters or to enable either transmitter TX1 30 or TX2 31 via direction switch 26 and associated control input 27 and to connect the receiver 32 via direction switch 28 and associated control input 28 with the transducers TD1 1 or TD2 2 ; The remaining 6 Bits are used to control the gain of the receiver-amplifier 32 ; Two transmitters 30 and 31 , the first TX1 connected to transducer 1 via cable 7 and the second TX2 connected via cable 8 to transducer 2 ; A pulse shaping network 33 , transforming the analog signal into a digital signal; Two enable D-flipflops 34a,b , with D-inputs controlled by Output Compare lines 36 and 37 . They are connected to shutter gates 39a and 39a , to control the pulses counted by the high resolution counters 40a and 40b . The D-Flipflops 34a and 34b are set by hardware output OC2 35 , from the Compare/Capture Unit 20 , controlled by switch 26 when at the same time either TX1 30 or TX2 31 is triggered. The enable flipflops close the corresponding shuttergate, when they get a clock pulse generated by the arriving direct or reflected signal/ The microprocessor 22 calculates when to switch the appropriate direction switches 26 and 28 to control the direction of travel and when to fire transducers 1 or 2; when to set and when to reset the time windows controlled by the Comparator-output OC3 36 to select the direct signal; or when to set and when to reset the time window to suppress the direct signal and select the signal reflected at the surface instead. All other circuits which are conventionally used with a microprocessor, such as RAM, EPROM, clock, auxiliary registers, etc. are shown as "SRAM/EPROM". We now describe the operation of the unit with specific reference to Figures 5, 6, 7a and 7b, which shows the timing and propagation times in relative juxtaposition, and the flowchart of the microprocessor software. 41a is the time of transmission for transducer 1, 41b the time for transmission of transducer 2. Both transmitters are edgetriggered using the COMPARE/CAPTURE UNIT of the microprocessor, when the relative timer output OC2 goes from high to low state. The time difference between both pulses is set by software to allow ample time between the first set of signals at transducer 2 and the following set of signals arriving at transducer 1. This time difference will be a couple of hundred microseconds. The pulses travel across the river in opposite directions and are picked up by transducer 2 at time 42a as direct signal, at time 43a as bottom returns and at 44a as surface reflected signal. The second pulse transmitted a little later arrives at transducer 1 at time 42b as direct signal, 43b as bottom returns and at time 44b as surface reflected signal. 45 shows the direction bit, causing the transmitter to transmit with transducer 1 and to receive with transducer 2 on high state and with the opposite transducer in low state. 47a and 47b show the action of the enable-flipflops 34a and 34b , which are enabled (allowed to clock when signal arrives) again by the COMPARE/CAPTURE UNIT processor controlled at different possible times: 1. either just before the direct signal at 47a or 47b with automatic reset 48a or 48b , respectively. 2. after the direct signal and the bottom returns have diminished and just before the surface reflected signal arrives (not explicitly shown). When the first transmitter 30 TX1 fires 41a , the shutter gate 49a is opened for the counting pulses coming from the clock 38 and which appear as 51a at the output of the shutter gate, which is connected to the high resolution counter 40a . The same description holds for the second transmitter 31 TX2, opening shutter gate 49b with counting pulses 51b applying to high resolution counter 40b . When the direct signals 42a or 42b arrive the processor had already noticed how many times (i) an overflow occurred. With a master clock of 4 MHz each overflow is 16 348 microseconds apart. Then the processor also reads the time comparison register 21 , which was set each time when one of both signals has arrived. By adding "i" times the overflow interval plus the time comparison readings the total propagation times T12 and T21 are calculated. The timing procedure for the surface reflected signals are precisely the same, except that the D-flipflops are enabled later, leaving the shutter gates open from 52a till 53a , and from 52b till 53b , respectively. The corresponding counting pulses are 54a and 54b , respectively. The embodiment of the inventions will be explained more by giving two numerical examples. The acoustic direct path will be 20.000 m long. To maintain sufficient time delay between the direct path, bottom reflecting path and the auxiliary path with the reflection at the surface has to meet certain depth requirements. To obtain a delay of at least 1.5 lamda (wavelength in the water), the minimum depth below and above the acoustic path must be according to Pythagoras with lamda = 0,0075 m and f = 200 kHz (transmitter frequency): h1 = 1.5 lamda P = 0.474 m The minimum depth of the acoustic transducers according to the invention: h2 = 5 lamda P = 0.866 m or about 2 h1 The actual depth is in fact 162 cm and the cross-section area A at the valid depth is 124 m 2 . The transducers are mounted 56 cm above the bottom. The following propagation times have been measured: T12 = 13,786904 ms T12r = 13,863416 ms T21 = 13,799308 ms T21r = 13,877344 ms The following velocities have been calculated: u1 = 0,652 m/s u2 = 0,728 m/s C = 1450,000 m/s Further: Wre1 = H1 + H2 2 = (1,0549894 + 1,0649935) = 1,05999m Spr = 20,112 m Then: u1/u2 = 0,896 ks' = 5,5 cm keis = 0,0 (known) v* = 4,6 cm vm = 0,668 m/s A = 123,99 m 2 interpolated from array at relative water level Wre1 = 1,06 Q = 82,83 m 3 /s The second example uses similar numbers however with ice cover. The following velocities have been calculated: u1 = 0,625 m/s u2 = 0,614 m/s u1/u2 = 1,0178 ks' = 5,5 cm (known from summer) keis = 3,0 cm v* = 4,6 cm vm = 0,588 m/s Q = 72,87 m 3 /s Due to the invention it is possible for the first time to determine the roughness of the bottom and the roughness of the ice. This leads to a significant improvement in the accuracy of the discharge measurement. Further, the real effective depth and cross-section area is used for the calculation of the discharge, which allows accurate calculation of the discharge and overcoming existing uncertainties and eliminating elaborate calibration with mechanical current meters. So far the invention was described with an example using a single path only. To overcome the sensitivity to changes of the direction of the stream often more complex arrangements are used. Popular configurations are V-configurations with a passive corner reflector on the opposite river bank and X-configurations with the active electronic responders. The invention can be applied in those configurations as well using the same basic principles as described.

Claims: A method for determining the flow velocity of a body of water flowing in a first direction, the method comprising: (a) emitting acoustic signals from an emitter to a receiver; (b) measuring propagation times of the acoustic signals propagated along a principal direct path (3) across and through the body of water, and (c) measuring propagation times of acoustic signals propagated along a path opposite to the principal direct path (3), characterised in that said method further comprises: (d) measuring the propagation times of acoustic signals propagated along an auxiliary indirect path (5) across and through the body of water, the auxiliary indirect path (5) comprising two segments of a reflection path of acoustic signals at a water-air or water-ice boundary of the body of water; and (e) utilising the measured propagation times of acoustic signals propagated along the principal path, opposite path and indirect path, to determine the flow velocity of the body of water in a first direction. A method according to Claim 1, characterised in that the principal direct path is inclined with respect to the first direction of flow of the body of water with an acute azimuth angle. A method according to Claim 2, characterised in that the acute azimuth angle is at least 5 degrees. A method according to Claim 2, characterised in that the acute azimuth angle is at most 60 degrees. A method according to any preceding Claim, characterised in that the acoustic signals propagated along and opposite to the principal direct path (3) and the acoustic signals propagated along the auxiliary indirect path (4) are propagated by the same signal propagation means (1,2). A method according to Claim 5, characterised in that the propagation means comprises a signal transducer (1,2). A method according to any preceding Claim, characterised in that the two segments of the auxiliary indirect path (5) comprise a first segment from a first transducer (1) to a reflection point at the water-air or water-ice boundary, and a second segment from the reflection point to a second transducer (2). A method according to any preceding Claim, characterised in that the auxiliary indirect path (5) is above the principal direct path (3). A method according to Claim 8, characterised in that the principal direct path (3) is arranged such that it is at a sufficient depth such that a reflected acoustic signal from the above auxiliary indirect path (5) is received later than a reflected acoustic signal from the principal direct path (3) or the path opposite the principal direct path (3), except when the body of water is lower than a predetermined level. A method according to any preceding Claim, characterised in that when measuring the propagation times of acoustic signals propagated along the principal direct path (3) only a direct signal triggers a signal detector and when measuring the propagation times of reflected acoustic signals propagated along the path opposite to the principal direct path (3) or the auxiliary path (5) only a reflected signal triggers the signal detector. A method according to Claim 10, wherein the approximate arrival time of the reflected signal is precalculated from the known mechanical length of the principal direct path (3), the measured velocity of sound in the principal direct path (3), the measured velocity of the water in the principal direct path (3), and the measured height of the water surface above the horizontal transducers (1,2); and that during the second measuring sequence the signal detector is enabled shortly before the expected approximate arrival time of the reflected signal. A method according to Claim 11, characterised in that velocities ul and u2, measured through the principal direct path (3) and the reflected path, respectively, are used to solve the following equation for ksm in the summer and for kice in the winter. A method according to Claim 12, characterised in that the following values are computed: V* = u1 A + 5,75 1g y1 ksm + kice(y1/(h-y1)) Q = Vm A(h)

Inventors: STEDNITZ WOLFGANG, STEDNITZ, WOLFGANG
Applicants: STEDNITZ WOLFGANG, STEDNITZ, WOLFGANG
Application year: 1994
Main classification: G01P5
IPC codes: G01P5