Water measurement apparatus and methods

Abstract

Water measurement apparatus and methods of operation. An exemplary implementation of water measurement apparatus comprises a measurement circuit outputting a voltage. A processor is operatively associated with the measurement circuit to receive the voltage and determine a capacitance value proportional to a water level. Computer readable program code is provided in computer readable storage and executable by the processor. The computer readable program code includes program code for determining the water level based on the capacitance value. Water measurement apparatus may also be implemented in an auto-network of water measurement stations.

Description

TECHNICAL FIELD

This invention relates generally to water measurement, and more particularly to water measurement apparatus and methods.

BACKGROUND

Water management is increasingly important as water supplies continue to be a limited resource for municipal, agricultural, and recreational purposes. Quantifying water supplies remains at the core of water management. However, quantifying water supplies often requires somebody visit the reservoirs and/or feeder streams in the water supply system and physically measure the water level. Of course this can be a time consuming process, particularly when the water supplies are spread out over a large geographic area, or in rough terrain. The water level may vary substantially between visits, making effective water management more difficult.

Accordingly, a number of water measurement devices are available to automatically measure and record water levels. One such device includes a floatation device slidably mounted around a pipe in the water. Fluctuations in the water level cause the floatation device to move up and down along the length of the pipe. The floatation device is connected to a strip chart recorder which produces markings that correspond to the water level indicated by the floatation device. Over time, however, the pipe may become corroded and impair movement of the floatation device. In addition, floatation devices are susceptible to damage and may need to be replaced. A chart house is also needed to protect the strip chart recorder from the environment (e.g., rain, snow, and wind). Furthermore, somebody still needs to visit the chart house periodically to retrieve the strip chart recordings and replace the ink cartridges and strip chart paper.

Ultrasonic measurement devices are also available. However, stray reflections increase the signal noise and decrease the reliability of the reading. Ultrasonic measurement devices depend on the velocity of sound and therefore are also sensitive to air temperature, humidity, and altitude. Measurement devices are also available that use microwaves. However, the electrical power requirements of these devices limit their use to areas with suitable power sources. Other devices may produce inaccurate measurements if the effects of water temperature and conduction (salinity) are not addressed.

SUMMARY

An exemplary implementation of a water measurement apparatus comprises a measurement circuit outputting a voltage. A processor is operatively associated with the measurement circuit to receive the voltage and determine a capacitance value proportional to a water level. Computer readable program code is provided in computer readable storage and executable by the processor. The computer readable program code includes program code for determining the water level based on the capacitance value.

An exemplary implementation of an auto-network of water measurement stations comprises at least a first and second water measurement station. A capacitive water measurement device at the first water measurement station determines a water level from a capacitance value based on a measured voltage. A transmitter at the first water measurement station is operatively associated with a receiver at the second water measurement station. The transmitter delivers water measurement data including at least the water level to the receiver at the second measurement station.

An exemplary method comprises: measuring a voltage between an inner conductor and an outer conductor, separating a real component from an imaginary component of the measured voltage, converting the measured voltage to a capacitance value using the real and imaginary component of the measured voltage, and determining the water level based on the capacitance value.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a high-level schematic diagram illustrating an exemplary implementation to auto-network a plurality of water measurement apparatus;

FIG. 2 is a graphical representation of a timing chart for operation of a water measurement apparatus in an exemplary auto-network;

FIG. 3 a is a perspective view of an exemplary water measurement apparatus as it may be installed on a post;

FIG. 3 b is a diagrammatic illustration of an exemplary water measurement apparatus;

FIG. 3 c is a cut-away perspective view of an exemplary water measurement apparatus;

FIG. 4 is a functional block diagram illustrating an exemplary implementation of control circuitry;

FIGS. 5-8 are circuit diagrams that may be used to implement an exemplary water measurement apparatus; and

FIG. 9 is a plot of a waveform showing sample points during operation of an exemplary water measurement apparatus.

DETAILED DESCRIPTION

Described herein are capacitive water measurement apparatus and methods to mitigate the effects of water salinity and the change in dielectric constant with temperature. Water measurement apparatus can be provided at low cost, is readily installed, and requires minimum maintenance. In addition, water measurement apparatus has low power requirements for operation. Water measurement apparatus may be implemented in an auto-networking environment to remotely transmit water measurement data (e.g., water level and corresponding time of the measurements) to a technician or other user. The water measurement apparatus may be implemented as a solid state device with no moving parts, eliminating mechanical failures, and is not readily corroded by water hardness. This and other implementations are described in more detail below with reference to the figures.

Implementations of the water measurement apparatus are accurate to at least ±0.25% accuracy for measurements of 0 to 4ft (i.e., ±0.01 ft, ±0.12 in) and have a resolution of at least 0.01 ft (0.12 in). In addition, implementations of the water measurement apparatus are insensitive to air temperature, water temperature, water salinity, and residual deposits (reducing the need for regular cleaning).

Data can be gathered on-site by read it directly from an optional LCD, or remotely by automatically transmitting data to a remote base station (e.g., accessible via the Internet). In addition, self-testing diagnostics may also be provided and may be implemented to remotely alert a user for maintenance or service, reducing the need for routing maintenance trips. The circuitry may also be field-programmable via RF commands.

Exemplary Auto-Network:

FIG. 1 is high-level schematic diagram illustrating an exemplary implementation to auto-network a plurality of water measurement apparatus. The auto-network may be implemented as a “call, talk, and hang-up” model. Data is forwarded in an optimal or near optimal path with some data being sent directly to an endpoint and other data passing through other measurement stations. Data can also be passed delivered to an accumulation point and modified (e.g., sending height, flow or volume data). The auto-network may also be auto-configured, so that if a gage is station is within a predetermined distance of another station, it is detected and added to the auto-network.

A water storage system 100 may include one or more water supplies 110a, 110b, 110c (also referred to generally by reference 110), such as, e.g., reservoirs, rivers, ditches, and/or streams. A plurality of stations 120a, 120b, 120c, 120d, 120e (also referred to generally by reference 120) may implement the water measurement apparatus, discussed in more detail herein, to record water levels of the water supplies 110.

One or more technicians (illustrated by vehicle 130) may be deployed to the stations 120 at various times to retrieve the water level data. In exemplary implementations, the recorded water levels are remotely transmitted as a data signal to the technician (e.g., to a laptop, PDA or other electronic device) without the technician having to visit each of the stations 120 individually. Such an implementation may be used, e.g., if the water storage system 100 is spread out over a large geographic area, or where the terrain is such that the technician cannot readily visit each of the stations 120 on a regular or semi-regular basis.

The stations may be auto-networked. Such an implementation may be used to transmit the recorded water level data to the technician if barriers 140, such as, e.g., mountains and/or buildings, exist in and around the water storage system 100 that block the data signal. A blocked data signal is illustrated in FIG. 1 with an “X” through signal path 150a, 150b.

Instead, the data signals are transmitted to other stations 120 in the auto-network using spread spectrum/frequency hopping to auto correct for RF congestion. For purposes of illustration, station 120a may transmit a data signal (e.g., including water level data for station 120a) to a technician at vehicle 130 via stations 120b-e over data path 160a-e. For example, station 120b may receive the data signal from station 120a and retransmit the data signal to station 120c, and so forth. Optionally, one or more of the stations 120b-e may include additional data in the data signal (e.g., water measurement data collected at these stations). Alternatively, communication can occur simultaneously between different stations.

FIG. 2 is a graphical representation of a timing chart 200 for operation of an exemplary auto-network. According to this implementation, the stations are in an active state to perform measurement and transmission operations, and are otherwise in an inactive state to reduce power consumption.

The timing chart 200 shows active states for stations 1-5 (e.g., stations 120a-e in FIG. 1). For purposes of illustration, Station 1 is in an active state from time t₁ to t₃, and is otherwise in an inactive state. During time t₁ to t₃, Station 1 makes a water level measurement and records it as a data signal. Transmission occurs during an overlap in the operation of Station 1 and Station 2. For example, in FIG. 2 transmission from Station 1 to Station 2 occurs on or after time t₂ when Station 2 becomes active, and prior to Station 1 becoming inactive at t₃.

In an exemplary implementation, the stations may be inactive or in a “sleep mode” every 1 to 2 seconds, and only need to be active during water level readings (e.g., <100 msec) and transmission operations (e.g., 15 msec), increasing the battery life and reducing maintenance. In other implementations, a plurality of water level measurements may be stored in memory and transmitted together. According to such an implementation, the active states of the stations do not have to overlap except during transmission from one station to another.

Exemplary Water Measurement Apparatus:

FIG. 3 a is a perspective view of an exemplary water measurement apparatus 300. Water measurement apparatus includes a measurement gage 310 that may be readily installed by attaching it to a post 320 (e.g., a T-post), e.g., using hose clamps or cable clamps 330 or other suitable fasteners. The post 320 may be driven into the ground so that the measurement gage 310 is in the water (e.g., as illustrated in FIG. 3b). Measurement gage 310 can be readily positioned up or down along the post 320 and rotated (e.g., to aim a directional antenna for an RF link).

Water measurement apparatus 300 may be located upstream from a weir and used to determine volumetric water flow (e.g., measured in cubic feet per second, CFS). Water law often refers to water in terms of acre feet (AF) available to users. A water manager may be obligated to deliver water (measured in CFS) for a defined period of time (often 24 hours).

Having said this, it is also noted that water measurement apparatus 300 does not need to be implemented with a weir. For example, the water measurement apparatus 300 may be implemented as a staff gage. Reservoir capacities may be predetermined and calibrated against the staff gage. The water height measured by the staff gage can be used to determine the stored water capacity (e.g., in AF) in the reservoir. With this data available at the beginning of the water season, the shareholders can determine in advance how much water they should expect for the season.

Measurement gage 310 is shown according to an exemplary implementation in more detail in FIG. 3b. Measurement gage 310 may be configured as co-axial cylinders, although other geometries are also possible. In one such implementation, measurement gage 310 includes an inner conductor 340 surrounded by insulating sheath 350, and an outer conductor 360. An optional opening (not shown) may be formed in the measurement gage 310 (e.g., in the top surface of the outer conductor) to where damping is a concern. The size of the opening may be increased for a faster response to changes in water height.

As water fills the region between the insulating sheath 350 and the outer conductor 360, the capacitance between the inner conductor 340 and outer conductor 360 increases because water has a significantly higher dielectric constant than air. The increase in capacitance is proportional to the water height 370 in the measurement gage 310. The reference (or zero) level 375 of the water produces a reference level of capacitance. When this capacitance is determined and subtracted from the measured capacitance (and the capacitance per unit of height is known), the water height can be computed.

In an exemplary implementation, the insulating sheath 350 is chosen for its lack of conduction current (i.e., it is a good insulator). Suitable materials for the insulation sheath may include, but are not limited to Teflon, PVC, or other insulating and low water-absorption materials (e.g., plastics). The resistive component of the insulating sheath (i.e., part of the sheath material) is determined by its displacement current (dielectric loss). Air is also a good insulator and for most practical purposes has no displacement current.

When more than one dielectric material is stacked between conductors, the capacitance is the same as the capacitance of each material connected in series. The capacitance of the insulating sheath 350 is in series with the capacitance of the water. When two materials are placed side-by-side between conductors, the capacitance is the same as the capacitance of each material connected in parallel. In the embodiment of the measurement gage 310 shown in FIG. 3b, the capacitance of the section filled with water is in parallel with the capacitance of the section filled with air (e.g., above the water level).

FIG. 3 c is a cutaway perspective view of an exemplary water height measurement apparatus showing the measurement gage 310 in more detail. The measurement gage may include a base portion 380 and a top portion 390. Base portion 380 may house the measurement capacitor 381. Measurement capacitor 381 may include an outer gage electrode 382 surrounding insulating sheath 350 and inner conductor 340. One or more openings 383 are formed in the outer conductor 360 to allow water to enter an integral stilling well 384 formed as part of the base portion 380 between the measurement capacitor 381 and the outer conductor 360.

Top portion 390 may house the circuitry, such as, e.g., a transceiver and battery board 391 and a measurement board 392. The circuit boards 391, 392 may be mounted in the top portion 390, e.g., by fasteners 393. Electrical connections (not shown) may also be provided between the circuit boards 391, 392 and the measurement capacitor 381 in base portion 380. The top portion 390 may be mounted to the bottom portion 380, e.g., by a threaded rod 395 using nuts 396 and compression fitting 397.

FIG. 4 is a functional block diagram illustrating exemplary control circuitry 400 that may be used to implement a water measurement apparatus. Control circuitry 400 may be provided, e.g., on a computer board 401 mounted in a protective housing to the water measurement apparatus (e.g., as illustrated in an exemplary implementation in FIG. 3c).

Control circuitry 400 is operatively associated with a gage or measurement circuit 410, discussed in more detail below. Control circuitry 400 may also be operatively associated with an optional transmitter 420 (e.g., RF transmitter) for remote data and program management (e.g., in the auto-network of FIG. 1). An optional user interface 430 (e.g., an interactive display implemented with an LCD) may also be provided.

Control circuitry 400 includes one or more processor 440 (or processing units), and computer-readable storage or memory 450 (e.g., Flash memory). Memory 450 may be used, for example, to store water level data (e.g., water height measurements). Processor is operatively associated with the measurement gauge 410, e.g., via converter 460 including a digital to analog (D/A) converter 461 and analog to digital (A/D) converter 462.

Control circuitry 400 also includes computer-readable program code for implementing operations on the processor 440. In an exemplary implementation, program code is provided for managing the application of the signal (e.g., sin wave) from the D/A converter 461, the measurement with the A/D converter 462, and storing water level data.

Program code is also provided for operations on the water level data using a Discrete Fourier Transform (DFT) algorithm, computing real and imaginary values of the gage measurement with respect to the phase of the input, and determining the capacitance based on the voltage measurement. These operations for determining the water height based at least in part on the measured voltage are described in more detail below. Optional program code may also be provided for determining the flow rate (e.g., over a weir) and/or the total volume of water (e.g., in a reservoir) based on one or more water height measurements.

In one embodiment, the program code may be implemented as scripts. Embodiments for controlling a device using scripts are described in co-pending, co-owned U.S. patent application entitled “DISTRIBUTED CONTROL SYSTEMS AND METHODS FOR BUILDING AUTOMATION” of Hesse, et al., filed on Apr. 24, 2003 (Ser. No. 10/422,525), which is hereby incorporated herein by reference for all that it discloses. The scripts may be defined based on various parameters, such as the needs and desires of those responsible for monitoring the water. The scripts can also be reconfigured based on the changing needs and/or desires of those responsible for monitoring the water.

Control circuitry 400 may also include a clock 470 (e.g., a battery-backed, real-time clock). Clock 470 may be used to record the time of a water height measurement. Clock 470 may also be used to initiate measurements and/or implement a timing schedule, such as the timing diagram illustrated in FIG. 2 for auto-networking.

Control circuitry 400 may also include a power management module 480 operatively associated with one or more power sources, such as, e.g., battery 481 and solar cell 482.

Exemplary Operations:

Briefly, the complex measurement method involves solving simultaneous equations for capacitance and resistance using the real and imaginary voltages. The real and imaginary voltage measurements are passed through a discrete Fourier Transform (DFT). The sine wave is smoothed by gage capacitance. The measurement shows real and imaginary values that are about equal. The result is independent of salinity and other factors affecting conduction current.

Temperature dependence may be reduced using a Teflon sleeve. The sensitivity of capacitance (C_t) may be determined as follows, where the thickness of the Teflon sleeve is chosen to determine S. $\begin{array}{c}{C}_t=\frac{C_{\mathrm{water}}{C}_{\mathrm{sheath}}}{C_{\mathrm{water}}+C_{\mathrm{sheath}}}\\ \mathrm{and}\\ S=\frac{C_{\mathrm{sheath}}}{C_{\mathrm{water}}+C_{\mathrm{sheath}}}\end{array}$

Where C_sheath for Teflon is 43.9 pF/in and C_water is 322.2 pF/in. Accordingly, S is 0.1199. The remaining variation may be compensated with an optional temperature sensor inside the inner tube.

The effects of residual deposits may also be reduced. Under the water, residual deposits have a much smaller dielectric constant than that of water and the effects on accuracy are negligible. Above water, the volume of water contributes to capacitance. Therefore, surface wetting effects on accuracy are negligible.

The following equations illustrate exemplary operations to implement a water measurement apparatus (such as the water measurement apparatus 300 shown in FIGS. 3a and 3b). In one embodiment, the operations may be implemented by control circuitry, such as the control circuitry 400 shown in FIG. 4.

A measurement circuit (e.g., gage 410 in FIG. 4) may be modeled according to one implementation by circuit 500 shown in FIG. 5 The total gage capacitance for circuit 500 is defined by equation (1): $\begin{array}{cc}{C}_{\mathrm{gage}}=\frac{C_{\mathrm{water}}{C}_{\mathrm{sheath}}}{C_{\mathrm{water}}+C_{\mathrm{sheath}}}+\frac{C_{\mathrm{air}}{C}_{\mathrm{sheath}}}{C_{\mathrm{air}}+C_{\mathrm{sheath}}}& \left(\mathrm{EQ}1\right)\end{array}$

When C_water is much larger than C_sheath the series capacitance is dominated by C_sheath. Also note that the capacitance of the empty section is dominated by C_air.

The gage capacitance may be corrected for temperature and salinity, as modeled in FIG. 6 by circuit 600. With regard to temperature, the dielectric constant of water is about 78.5 at room temperature and varies from about 85 near freezing to about 56 near the boiling point. The measured capacitance depends on the temperature of the water. However, if the geometries of the insulating sheath and the space that fills with water are chosen so that C_water is much greater than C_sheath and the empty space such that C_air is much less than C_sheath, the measured value of capacitance depends primarily upon the characteristics of C_sheath. It should be noted that the temperature need not be measured to compensate for its effect when a stable material (e.g., Teflon) is used for the sheath.

The conductivity of the water also varies with salinity. The measurement method is preferably insensitive to changes in R_water. Measurement techniques may be employed to separate resistance and capacitance in parallel, such as, e.g., those known in analog in-circuit testing.

The gage capacitance may also be corrected for conductivity and salinity of the water. The water also has a loss component in parallel (whether by conduction current, displacement current, or both) with the dielectric.

The gage may be modeled using the equivalent capacitance and dielectric loss, as shown by circuit 700 in FIG. 7. The circuit 700 is connected into a circuit 800 shown in FIG. 8 to measure voltage. The voltage, V_b, is (by direct circuit analysis): $\begin{array}{cc}\begin{array}{c}{V}_b=\frac{R_p+\frac{1}{j\omega C_{\mathrm{gage}}}}{R_p+\frac{1}{j\omega C_{\mathrm{gage}}}}\left(I_s\right)\\ \mathrm{Where}\\ I_s=\frac{V_a-V_b}{R_s}\end{array}& \left(\mathrm{EQ}2\right)\end{array}$

Expanding equation 2, the real component and the imaginary component of this equation form another system of equations which can be directly solved for the gage capacitance and parallel resistance. The solution for gage capacitance is: $\begin{array}{cc}{C}_{\mathrm{gage}}=\frac{-\mathrm{Im}\left\{V_b\right\}{V}_a}{\omega R_s\left({\left(\mathrm{Re}\left\{V_b\right\}\right)}^2+{\left(\mathrm{Im}\left\{V_b\right\}\right)}^2\right)}& \left(\mathrm{EQ}3\right)\end{array}$

The solution for parallel resistance is: $\begin{array}{cc}{R}_p=\frac{-\left({\left(\mathrm{Re}\left\{V_b\right\}\right)}^2+{\left(\mathrm{Im}\left\{V_b\right\}\right)}^2\right)R_s}{{\left(\mathrm{Re}\left\{V_b\right\}\right)}^2-\mathrm{Re}\left\{V_b\right\}{V}_a+{\left(\mathrm{Im}\left\{V_b\right\}\right)}^2}& \left(\mathrm{EQ}4\right)\end{array}$

The results are readily checked in Matlab as follows: $\%\mathrm{Script}\mathrm{to}\mathrm{check}\mathrm{derivation}$$\%\mathrm{Define}\mathrm{component}\mathrm{values}$$\mathrm{Rs}=10000;$$\mathrm{Rp}=1000;$$\mathrm{Cgage}=22\text{-}e\text{-}12;$$Va=5;$$F=10000;$$\%\mathrm{Compute}\mathrm{the}\mathrm{impedance}\mathrm{of}R\mathrm{in}\mathrm{parallel}\mathrm{with}\mathrm{Cgage}$$W=2*\mathrm{pi}*f;$$\mathrm{Xcp}=1/\left(j*w*\mathrm{Cgage}\right);$$\mathrm{Zp}=\mathrm{Rp}*\mathrm{Xcp}/\left(\mathrm{Rp}+\mathrm{Xcp}\right);$$\%\mathrm{Compute}\mathrm{the}\mathrm{voltage}Vb$$Vb=\left(\mathrm{Zp}/\left(\mathrm{Rs}+\mathrm{Zp}\right)\right)*Va$$\%\mathrm{Find}\mathrm{the}\mathrm{real}\mathrm{and}\mathrm{imag}\mathrm{values}\mathrm{of}Vb$$v1=\mathrm{real}\left(Vb\right);$$v2=\mathrm{imag}\left(Vb\right);$$\%\mathrm{Compute}\mathrm{Cgage}\mathrm{from}\mathrm{measured}\mathrm{voltages}$$\mathrm{Cgage\_meas}=-Va*v2/(w*\mathrm{Rs}*\left(v1^2+v2^2\right)\%\mathrm{Compute}\mathrm{Rp}\mathrm{from}\mathrm{measured}\mathrm{voltages}\mathrm{Rp\_meas}=-\left(v1^2+v2^2\right)*\mathrm{Rs}/\left(v1^2-Va*v1+v2^2\right)$

For a gage capacitance of 220 pF, a source resistance of 10K, and a parallel resistance of 10K, and a frequency of 10 KHz checks as: $Vb=0.4545-0.0057i$$\mathrm{Cgage\_meas}=2.2000e^{-010}$$\mathrm{Rp\_meas}=1.0000e^{+003}$

Changing the frequency to 1 KHz results in a more reasonable value of V_b but produces no other changes. Similar results were produced by varying the other parameters.

The discrete Fourier transform of a signal (DFT) is defined as: $\begin{array}{cc}V\left(k\Delta \omega \right)=\frac{j\sqrt{2}}{N}\sum _{n=0}^{N-1}V\left(n\right)e^{\frac{-j2\pi nk}{N}}& \left(\mathrm{EQ}5\right)\end{array}$

Where N is the number of samples of the waveform and delta co is w the change in radian frequency. If the time between samples is delta t, the change in radian frequency is defined as: $\begin{array}{cc}\Delta \omega =\frac{1}{N\Delta t}& \left(\mathrm{EQ}6\right)\end{array}$

The actual radian frequency is: $\begin{array}{cc}\omega =k\Delta \omega =\frac{k}{N\Delta t}& \left(\mathrm{EQ}7\right)\end{array}$

The integer k is chosen so that omega (w) is the desired frequency.

Exemplary Calculation

The parameters for the check solution are defined as: $\%\mathrm{Define}\mathrm{the}\mathrm{parameters}$$N=100;$$\mathrm{delta}T=1.12e^{-4}$$\mathrm{amplitude}=1;$$\mathrm{phaseAngle}=45;$$\mathrm{approx\_freq}=1000;$

Here N is the number of samples, deltaT is the time between samples, amplitude is the RMS value of the voltage, phaseAngle is the phase angle of the signal with respect to a sine wave, and approx_frequency is the target frequency. The actual frequency is the nearest integer frequency (because of deltaT) to the target frequency.

The vector of times and the actual frequency is determined by: $\%\mathrm{Computer}\mathrm{the}\mathrm{time}\mathrm{vector}$$T=0;\mathrm{delta}T;\left(\mathrm{Nsamples}-1\right)*\mathrm{delta}T;$$K=\mathrm{round}\left(F/\mathrm{deltaF}\right);$$\mathrm{Freq}=k*\mathrm{deltaF};$

The signal (as a function of time) is computed by: $\%\mathrm{Compute}\mathrm{the}\mathrm{time}\mathrm{varying}\mathrm{voltage}\mathrm{vector}$$v=\mathrm{sqrt}\left(2\right)*\mathrm{amplitude}*\sin \left(2*\mathrm{pi}*\mathrm{Freq}*t+\mathrm{pi}*\mathrm{phaseAngle}/180\right);$

The discrete fourier transform (DFT) of the signal is computed by: $\%\mathrm{Computes}\mathrm{the}\mathrm{discrete}\mathrm{fourier}\mathrm{transform}\mathrm{at}\mathrm{one}\mathrm{frequency}$$n=0;N\text{-}1;$$V\mathrm{ac}=\mathrm{sqrt}\left(2\right)*j*v*\exp \left(j*2*\mathrm{pi}*k*n*/N\right)N;$

Note that in Matlab, v is a row vector and exp(j*2*pi*k*n/N)' is a column vector. When these two quantities are multiplied together, the result is the sum of the product of the elements.

The real and imaginary values are equal and positive (as we expect them to be at a 45 degree angle). The magnitude of the voltage (sqrt(Vreal{circumflex over ( )}2+Vimag{circumflex over ( )}2)) is equal to 1 which was the value set.

FIG. 9 is a waveform showing sample points during operation of an exemplary water measurement apparatus. Note that there are multiple cycles of the input waveform sampled. The data for the D/A contains all the points shown as circles.

For computational convenience, the real and imaginary values can be separately computed by noting that: $\begin{array}{cc}V\left(k\Delta \omega \right)=\frac{j\sqrt{2}}{N}\sum _{n=0}^{N-1}V\left(n\right)e^{\frac{-j2\pi nk}{N}}& \left(\mathrm{EQ}8\right)\\ =\frac{j\sqrt{2}}{N}\sum _{n=0}^{N-1}V\left(n\right)\left(\cos \left(\frac{2\pi nk}{N}\right)+j\sin \left(\frac{2\pi nk}{N}\right)\right)& \\ =\frac{\sqrt{2}}{N}\sum _{n=0}^{N-1}V\left(n\right)\left(j\cos \left(\frac{2\pi nk}{N}\right)-\sin \left(\frac{2\pi nk}{N}\right)\right)& \\ \mathrm{OR}& \\ V\left(k\Delta \omega \right)=\frac{-\sqrt{2}}{N}\sum _{n=0}^{N-1}V\left(n\right)\sin \left(\frac{2\pi nk}{N}\right)+& \left(\mathrm{EQ}9\right)\\ j\frac{\sqrt{2}}{N}\sum _{n=0}^{N-1}V\left(n\right)\cos \left(\frac{2\pi nk}{N}\right)& \end{array}$ Exemplary Design Parameters

The empty value of capacitance is: $\mathrm{Design}\mathrm{Inputs}\text{:}$$\mathrm{Dielectric}\mathrm{constants}\text{:}$${\varepsilon }_0:=8.854*{10}^{-12}\mathrm{farad}\text{/}m$${\varepsilon }_{\mathrm{water}}:=78.54$${\varepsilon }_{\mathrm{pvc}}:=4.5$${\varepsilon }_{\mathrm{air}}:=1$$\mathrm{Outside}\mathrm{diameter}\mathrm{of}\mathrm{the}\mathrm{inner}\mathrm{conductor}\text{:}$$r_{\mathrm{ic}}:=\frac{0.875}{2}\mathrm{in}$$\mathrm{Inside}\mathrm{diameter}\mathrm{of}\mathrm{the}\mathrm{outer}\mathrm{conductor}\text{:}$$r_{\mathrm{oc}}:=\frac{1.500}{2}\mathrm{in}$$\mathrm{Inside}\mathrm{diameter}\mathrm{of}\mathrm{the}\mathrm{PVC}\mathrm{pipe}\text{:}$$r_{\mathrm{ipvc}}:=\frac{0.930}{2}\mathrm{in}$$\mathrm{Outside}\mathrm{diameter}\mathrm{of}\mathrm{the}\mathrm{PVC}\mathrm{pipe}\text{:}$$r_{\mathrm{opvc}}:=\frac{1.050}{2}\mathrm{in}$$\mathrm{Length}\mathrm{of}\mathrm{the}\mathrm{gage}\text{:}$$L:=22\mathrm{in}$

The capacitance for a gage height H is: $\mathrm{Height}\mathrm{of}\mathrm{water}\mathrm{in}\mathrm{gage}\text{:}H\equiv 0.0001\mathrm{inch}$$C_{\mathrm{gage}}=6.992\times {10}^{-11}F$

The full scale value of capacitance is: $\mathrm{Design}\mathrm{Inputs}\text{:}$$\mathrm{Dielectric}\mathrm{constants}\text{:}$$\begin{array}{c}{\varepsilon }_0:=8.854*{10}^{-12}\mathrm{farad}/m\\ {\varepsilon }_{\mathrm{water}}:=78.54\\ {\varepsilon }_{\mathrm{pvc}}:=4.5\\ {\varepsilon }_{\mathrm{air}}:=1\end{array}\mathrm{Outside}\mathrm{diameter}\mathrm{of}\mathrm{the}\mathrm{inner}\mathrm{conductor}\text{:}$$r_{\mathrm{ic}}:=\frac{0.875}{2}\mathrm{in}\mathrm{Inside}\mathrm{diameter}\mathrm{of}\mathrm{the}\mathrm{outer}\mathrm{conductor}\text{:}$$r_{\mathrm{oc}}:=\frac{1.500}{2}\mathrm{in}$$\mathrm{Inside}\mathrm{diameter}\mathrm{of}\mathrm{the}\mathrm{PVC}\mathrm{pipe}\text{:}$$r_{\mathrm{ipvc}}:=\frac{0.930}{2}\mathrm{in}$$\mathrm{Outside}\mathrm{diameter}\mathrm{of}\mathrm{the}\mathrm{PVC}\mathrm{pipe}\text{:}$$r_{\mathrm{opvc}}:=\frac{1.050}{2}\mathrm{in}$$\mathrm{Length}\mathrm{of}\mathrm{the}\mathrm{gage}\text{:}$$L:=22\mathrm{in}$

The capacitance for a gage height H is: $\mathrm{Height}\mathrm{of}\mathrm{water}\mathrm{in}\mathrm{gage}\text{:}H\equiv 18\mathrm{inches}$$C_{\mathrm{gage}}=2.878\times {10}^{-10}F$

It is readily apparent that the water measurement apparatus and methods of the present invention represent important developments in the field of water monitoring. Having herein set forth exemplary implementations, it is anticipated that suitable modifications can be made thereto which will nonetheless remain within the scope of the invention.

Claims

1. A water measurement apparatus comprising: a measurement circuit outputting a voltage; a processor operatively associated with the measurement circuit to receive the voltage and determine a capacitance value proportional to a water level; and computer readable program code provided in computer readable storage and executable by the processor, the computer readable program code including program code for determining the water level based on the capacitance value.

2. The water measurement apparatus of claim 1 wherein the measurement circuit includes an inner conductor and an outer conductor.

3. The water measurement apparatus of claim 2 wherein the inner conductor is surrounded by an insulating sheath.

4. The water measurement apparatus of claim 2 wherein the voltage is measured between the inner conductor and the outer conductor.

5. The water measurement apparatus of claim 2 wherein the capacitance value changes in proportion to the water height between the inner conductor and the outer conductor as indicated by the voltage.

6. The water measurement apparatus of claim 1 wherein the computer readable program code includes program code for determining a real component and an imaginary component of the voltage.

7. The water measurement apparatus of claim 1 further comprising program code for correcting the capacitance value for water conductivity.

8. The water measurement apparatus of claim 1 further comprising program code for correcting the capacitance value for water temperature.

9. The water measurement apparatus of claim 1 further comprising program code for correcting the capacitance value for water salinity.

10. The water measurement apparatus of claim 1 further comprising recording the water level and a corresponding sample time in the computer readable storage.

11. The water measurement apparatus of claim 1 further comprising a transmitter at a first station operatively associated with a receiver at a second station to deliver water measurement data including at least the water level from the first station to the second station.

12. The water measurement apparatus of claim 11 wherein the transmitter and the receiver are auto-networked with one another.

13. An auto-network of water measurement stations comprising: at least a first and second water measurement station; a capacitive water measurement device at the first water measurement station, the capacitive water measurement device determining a water level at the first water measurement station from a capacitance value based on a measured voltage; and a transmitter at the first water measurement station operatively associated with a receiver at the second water measurement station, the transmitter delivering water measurement data based on the water level to the receiver at the second measurement station.

14. The auto-network of water measurement stations of claim 13 wherein the water measurement data includes the water level and a corresponding sample time.

15. The auto-network of water measurement stations of claim 13 wherein at least the first measurement station toggles between an active state and an inactive state to conserve power.

16. The auto-network of water measurement stations of claim 13 wherein the first and second measurement stations toggle between an active state and an inactive state, the active state of the first measurement station overlapping in time with the active state of the second measurement station.

17. The auto-network of water measurement stations of claim 13 wherein the capacitive measurement device includes a measurement circuit outputting a voltage corresponding to the water level.

18. The auto-network of water measurement stations of claim 13 wherein the capacitive measurement device includes a processor executing computer readable program code including program code for determining the water level based on a measured voltage.

19. The auto-network of water measurement stations of claim 13 wherein the capacitive measurement device includes a measurement circuit outputting a voltage corresponding to the water level.

20. A method comprising: measuring a voltage between an inner conductor and an outer conductor; separating a real component and an imaginary component of the measured voltage; converting the measured voltage to a capacitance value using the real and imaginary component of the measured voltage; and determining the water level based on the capacitance value.

21. The method of claim 20 further comprising determining a flow rate based on the capacitance value.

22. The method of claim 20 further comprising determining a volume of water based on the capacitance value.

23. The method of claim 20 wherein the capacitance value changes in proportion to changes in the water height between the inner conductor and the outer conductor.

24. The method of claim 20 further comprising recording the water level and a corresponding sample time.

25. The method of claim 20 further comprising transmitting water measurement data including at least the water level in an auto-network of water measurement stations.