Bubble water depth measuring method and system thereof

Abstract

A water depth measuring method comprises steps of making a bundle of reference water column pressure measuring tubes of a U character form of the number of n≧2 at a same length with a water depth measuring tube in order to compensate for a total error ΣΔx of a water depth measurement; filling an amount of water into the reference water column pressure measuring tubes to form water columns ho1, ho2, . . . hon therein; measuring the water pressure in the reference water column pressure measuring tube by a pressure transducer before measuring a water depth; subtracting previously known water column pressure γhoi and γhoj derived from the measured water pressure values Pmi, and Pmj to calculate total errors ΣΔi and ΣΔj; calculating a total error ΣΔx based on the total errors ΣΔi and ΣΔj, in which the total error ΣΔx was generated upon measuring of the bubble generating pressure Pmx in the water column pressure measuring tubes; subtracting the total error ΣΔx from the bubble generating pressure Pmx, dividing the calculated water column pressures by the water column pressure γhoi(j) in the reference water column pressure measuring tube having most approximate value thereto and multiplying the re-calculated value by the water column pressure hoi(j) thereby to measure a water depth in a higher accuracy, in which the water depth measuring tube has a length corresponding to a water depth, the reference water column pressure measuring tubes has the same inner diameter as that of the water depth measuring tube coupled therewith, in which a diaphragm is connected through an electromagnetic valve to one tube of the reference water column pressure measuring tube to adjust the speed of the compressed gas from a buffer tank to be supplied to the reference water column pressure measuring tubes, a transparent container is connected to the other tube of the reference water column pressure measuring tubes to measure an amount of water therein and check it, periodically, and includes a nipper for supplementing water reduced due to the evaporation, a tube is directly mounted above the transparent container and has the same inner diameter as that of the reference water column pressure measuring tubes and a length portion, on the surface of which scales are formed, a thin film tube is made of rubber and directed to the upper portion of the tube to be expanded by a gas pressure, and a manual valve connected to the upper end of the thin film tube, and a system thereof.

Description

BACKGROUND OF THE INVENTION

[0001] The invention is related to providing a bubble water depth measuring method for compensating for a total error of a bubble water depth or water level measurement in a reservoir that a water level shifts in a wider range to enable the measurement of the water depth or water level in a higher accuracy, using an industrial pressure transducer, in which the pressure transducer is manufactured in an environment that a circumferential air temperature is changed in the range of −40° C. to +50° C., and a system thereof.

[0002] Particularly, the invention is related to providing a bubble water depth measuring method for measuring a water column pressure to enable the remote-measurement of a water depth or water level in a reservoir, a lake and a river. Herein, it is noted that only a water depth measurement now will be explained below, because a bubble water level meter is a device for measuring a water depth and calculating it into a water level.

PRIOR ARTS

[0003] There has been much interested in a bubble water level meter that can be used for a hydrology observatory, because the bubble water level meter can measure a water depth under an icy, if a water surface is frozen in a reservoir, a lake and a river. The water level meter for the hydrology observatory must be operated under the conditions as follows:

[0004] In a reservoir the changing width of a water level is usually up to 10 m.

[0005] There are many areas where a circumferential air temperature of a water level observatory post is in the range of −40° C. to +50° C. in seasons.

[0006] Most of the water level observatory posts are not equipped with not warming and cooling utilities but also a power source unit.

[0007] It is now tendency that a demand of a remote-water level measuring system is being increased because many water level observatory posts are established in uninhabited places.

[0008] In light of these facts, there is a problem in that a bubble type water level meter has a lower accuracy under the conditions of the severe weather change in seasons and a larger water level change width.

[0009] The bubble type water level meter has advantages in that its installing, maintenance and operating costs are relatively lower and can a water level even when a reservoir and a river are frozen in winter, but it is broadly not used as a hydrology observatory because its measuring error is larger.

[0010] Also, the bubble type can be operated in stability comparing with other water level meters in a river that the concentration of floating particles is higher and a sand and earth layer is swiftly changed and saves the maintainnence and operating costs, steeply, but its measuring accuracy is low.

[0011] Referring to FIG. 1, error factors that are caused upon measuring a water depth according a bubble generating method will be explained below, and an unit of a water column pressure will be represented as cmH2O or mmH2O for the purpose of consulting the convenience in calling the term “error” as an abbreviate of “absolute error”. 1 is a water column pressure tube, which will be called “a water depth measuring tube”, 10 is a compressed gas generating device and 3 is a pressure transducer. A pressure applied to the lower end portion of the water depth measuring tube 1 is as follows:

Pc=γhx+Pa,   (1)

[0012] Wherein, hx is an altitude difference of water filled in the water depth measuring tube 1, which is considered as a water depth. γ is a specific gravity (gm/cm3) of water. γhx is a water column pressure. Pa is an atmosphere pressure on a water surface.

[0013] The pressure transducer 3 measures a surplus pressure, not for an atmosphere pressure (P≈Pab−Pa; Pab-Absolute pressure). The water depth measuring tube 1 filled with water generates bubbles at the lower end portion with compressed gas being supplied to the upper end portion thereof. Upon generating of the bubbles, the pressure transducer 3 measures the compressed gas pressure, and its result is as follows:

P mx =γh x +Δpx+(Pa−Pao)−γgHo+Δpb;  (2)

[0014] Wherein, Δpx is an absolute error of the pressure transducer 3 at the time of measuring Pmx. ΔPa=(Pa−Pao) is a difference between an atmosphere Pa on the water surface and an atmosphere Pao applied to the pressure transducer 3. Generally, Pao≠Pa because a water level observatory post is placed on a much higher position than a water surface. Ho is an altitude difference between the lower end portion of the water depth measuring tube 1 and the mounting position of the pressure transducer 3. γg is a density of compressed gas to be supplied to the water depth measuring tube 1. Δpb is a pressure of a supplementary compressed gas changed according to the bubble pressure that is formed at the lower end portion of the water depth measuring tube 1. Δpb will be ignored because its reducing method is now developed. In expression (2), all items are a measuring error of a water column pressure γhx except for the water column pressure γhx. What the measuring errors are summed up is assumed as total absolute error ΣΔx of the water column pressure γhx as follows:

ΣΔx=Δpx+ΔPa−γgHo,  (3)

[0015] Wherein, ΔPa=Pa−Pao=γα(Ho−hx), γa is an air density, which can be seen as 1${\gamma }_g\approx {\gamma }_{\mathrm{go}}\frac{P_{\mathrm{mx}}+P_a}{P_a}.$

[0016] γgo is a density of a compressed gas to be used at an atmosphere. Un the condition that a water depth is changed over 10 m in a reservoir, the altitude difference between the lower end portion of the water depth measuring tube 1 and the mounting position of the pressure transducer 3 usually becomes Ho≧20 m.

[0017] Looking into ΔPa and γgHo, if Ho=20 m=2000 cm, hx is changed in the range of 200 to 1000 cm, an air density γα=1.2·10−3 gm/cm3 and a water depth is measured using a compressed gas, the change of ΔPa is as follows:

ΔPa=1.2·10−3(2000-200)=2.16 gm/cm2≈2.16 cmH2O(if hx=200 cm)

ΔPa=1.2·10−3(2000-200)=1.2 gm/cm2≈1.2 cmH2O(if hx=1000 cm)

ΔgHo=1.2·10−3·1.2·2000=2.88 gm/cm2≈2.9 cmH2O(if hx=200 cm)

ΔgHo=1.2·10−3·2·2000=4.8 gm/cm2≈4.8 cmH2O(if hx=1000 cm)

[0018] Herein, it is known that when hx=200˜1000 cm, the error changing range is 2.16-2.9≈−0.74 cm to 1.2-4.8≈−3.6 cm due to ΔPa−γgHo. If an allowance error of the water depth, water level is ±1 cm, the error component of ΔPa−γgHo cannot be ignored. Of course, γa was ignored though it is changed according to a temperature.

[0019] The pressure transducer 3 includes a temperature compensation circuit for correcting the property that a pressure sensor mounted therein is changed according to a temperature. But the temperature compensation circuit is lacking of compensating for an error changed according to a temperature t and a pressure P to be measured, perfectly. An error of a curve Δp=f(t;P) always happens.

[0020] For example, when the pressure transducer (Model PTX1000) that secures a pressure error δp (=±0.25% fs) at room temperature (t≈20-24° C.) measures a pressure of P=1000 cm H2O with being cooled at −10° C., its measuring error Δp≅−31 cm H2O, even though its use temperature range is introduced as −40° C.˜+90° C. More precisional pressure transducer (Model PDCR862) has a useable temperature range of −54° C.˜+125° C. and a measuring error δp of ±0.1% fs. Therefore, Δp=8 gm/cm2, and an absolute error hx is ˜8 cm. Of course, if only a gas temperature is changed without cooling or heating the pressure transducer on the whole, the error is reduced. But, when a bubble water level meter is installed outdoor, the pressure transducer is cooled or heated on the whole. Under the condition that a circumferential air temperature is changed in the range of −40° C.˜+50° C., there doesn't exist any pressure transducer for securing Δp<±1 gm/cm3 in full range to be measured.

[0021] Therefore, even though it is possible to compensate for the error ΣΔx caused at the time of measuring the water column pressure γhx, perfectly, when the water depth hx is calculated into 2$\mathrm{hx}=\frac{\gamma \mathrm{hx}}{{\gamma }_l},$

[0022] knowing nothing about an average density γ of water causes a corresponding error h′x as follows: 3$h_x^{\prime }=\frac{\gamma h_x}{{\gamma }^{\prime }};\delta h_x^{\prime }=\frac{h_x^{\prime }}{h_x}-1=\left(\frac{\gamma }{{\gamma }^{\prime }}-1\right)={\delta }_{\gamma }$

[0023] Therefore, a total error of the measurement of the water depth hx is as follows: 4$\begin{array}{cc}{h}_x^{\prime }=\frac{\gamma h_x+{\mathrm{\Sigma \Delta }}_x}{{\gamma }^{\prime }};{\delta }_{h^{\prime }x}=\frac{h_x^{\prime }}{h_x}-1={\delta }_{\gamma }+\frac{{\mathrm{\Sigma \Delta }}_x}{\gamma h_x}={\delta }_{\gamma }+{\mathrm{\Sigma \delta }}_x& \left(4\right)\end{array}$

[0024] Wherein, an error δh′x is a parameter that varies according to a circumferential temperature and a measuring range, but not integer. Due to it, it cannot be simply compensated. Therefore, the error inquired above is one that a bubble type water depth measuring method has. There has been known a bubble type water depth, water level meter for excluding or compensating for δγ and Σδx. Only, the bubble water measuring methods are disclosed in Publications as follows:

[0025] 1) U.S. Pat. No. 5,791,187 issued on Aug. 11, 1998, which is titled “Level Measurement Method using Measurements of Water Column Pressure Therefor”

[0026] 2) Canadian Patent No. 2,171,801 on Jan. 23, 2000, which is titled “Level Measurement Method using Measurements of Water Column Pressure Therefor”

[0027] 3) Gerrnany Patent DE19620656C2 issued on Mar. 2, 2000

[0028] 4) Japanese Patent No. 2855423 issued on Feb. 10, 1999

[0029] 5) South Korean Patent No. 185260 issued on Dec. 28, 1998

[0030] These known corresponding publications disclose a method of measuring a water depth hx, which will be explained with reference to FIG. 2. 2 is a water column pressure tube that has a shorter length than a water depth measuring tube 1 by Δh. A water pressure of the water column pressure tube 2 filled with water is as follows: γh2=γ(h2−Δh). The water depth measuring tube 1 and the water column pressure measuring tube 2 are connected through valves 51 and 52 to a buffer tank 4. Compressed gas is supplied through a valve 50 to the buffer tank 4. 3 is a pressure transducer.

[0031] As a part for enhancing the accuracy of the water depth measurement the arrangement uses two measuring tubes 1 and 2 to measure a water specific gravity γ that is changed according to a temperature and components of water and then divide a water column pressure γhx by the water specific gravity γ. But, it is interested in that the separate measurement of the water specific gravity γ is not necessary. If only the water depth hx is measured, it can be measured in a simple and effective expression as follows: 5$\begin{array}{cc}{h}_x=\frac{\Delta h}{1-\frac{P_{\mathrm{m2}}}{P_{\mathrm{mx}}}\left(1-{\gamma }_{\mathrm{ao}}\frac{T_o}{T}\Delta h\right)};& \left(5\right)\end{array}$

[0032] Wherein, Pm2 is a bubble pressure generated at the water column pressure measuring tube 2, γao is a specific gravity of air, and 6${\gamma }_{\mathrm{ao}}\frac{T_o}{T}\Delta h$

[0033] is ignored because of a very little value. Assuming that ΣΔx, ΣΔ2=0, the expression (5) is as follows: 7$h_x=\frac{P_{\mathrm{mx}}}{P_{\mathrm{mx}}-P_{\mathrm{m2}}}\Delta h=\frac{\gamma h_x·\Delta h}{\gamma h_x-\gamma \left(h_x-\Delta h\right)}=\frac{\gamma h_x\Delta h}{\mathrm{\gamma \Delta }h}=h_x$

[0034] Herein, it is noted that the water specific gravity γ is entirely excluded from the expression (5). It is because hx and h2 has a difference Δh therebetween, and it is seen as γx=γ2. But, because ΣΔx≠0, ΣΔ2≠0, a water depth h′x is as follows: 8$\begin{array}{cc}{h}_x^{\prime }=\frac{\gamma h_x+{\mathrm{\Sigma \Delta }}_x}{\gamma h_x+{\mathrm{\Sigma \Delta }}_x-\gamma h_2-{\mathrm{\Sigma \Delta }}_2}·\Delta h=\frac{\gamma h_x\left(1+\frac{{\mathrm{\Sigma \Delta }}_x}{\gamma h_x}\right)}{\mathrm{\gamma \Delta }h\left(1+\frac{{\mathrm{\Sigma \Delta }}_x-{\mathrm{\Sigma \Delta }}_2}{\mathrm{\gamma \Delta }h}\right)}\Delta h=\frac{h_x\left(1+\frac{{\mathrm{\Sigma \Delta }}_x}{\gamma h_x}\right)}{\left(1+\frac{{\mathrm{\Sigma \Delta }}_x-{\mathrm{\Sigma \Delta }}_2}{\mathrm{\gamma \Delta }h}\right)}& \left(6\right)\end{array}$

[0035] Herein, assuming that 9$\frac{{\mathrm{\Sigma \Delta }}_x-{\mathrm{\Sigma \Delta }}_2}{\mathrm{\gamma \Delta }h}<<1.0$

[0036] is very little value, it is ignorable. The expression (6) is represented as follows: 10$\begin{array}{cc}{h}_x^{\prime }=h_x\left(1+\frac{{\mathrm{\Sigma \Delta }}_x}{\gamma h_x}\right)\left(1-\frac{{\mathrm{\Sigma \Delta }}_x-{\mathrm{\Sigma \Delta }}_2}{\mathrm{\gamma \Delta }h}\right)\approx h_x\left(1+\frac{{\mathrm{\Sigma \Delta }}_x}{\gamma h_x}-\frac{{\mathrm{\Sigma \Delta }}_x-{\mathrm{\Sigma \Delta }}_2}{\mathrm{\gamma \Delta }h}\right)& \left(7\right)\end{array}$

[0037] Therefore, an absolute error Δh′x that is resulted from the measurement of the water is as follows: 11${\mathrm{\gamma \Delta }}_{h^{\prime }x}={\mathrm{\Sigma \Delta }}_x-\left({\mathrm{\Sigma \Delta }}_x-{\mathrm{\Sigma \Delta }}_2\right)\left(\frac{h_x}{\Delta h}\right)=\Delta \mathrm{px}+\Delta \mathrm{Pax}-{\gamma }_{\mathrm{gx}}{H}_o-\left[\Delta p_x+\Delta P_{\mathrm{ax}}-{\gamma }_{\mathrm{gx}}{H}_o-\Delta p_2-\Delta {\mathrm{Pa}}_2+{\gamma }_{\mathrm{g2}}{H}_o\right]\frac{h_x}{\Delta h}=\Delta \mathrm{px}+\Delta \mathrm{Pax}-{\gamma }_g{H}_o-\left[\Delta \mathrm{px}-\Delta p_2+\Delta \mathrm{pa}+\left({\gamma }_g-{\gamma }_{\mathrm{gx}}\right)H_o\right]\frac{h_x}{\Delta h};$

[0038] Herein, considering that ΔΔpa (=ΔPax−ΔPa2) and (γg2−γgx)Ho is a very little value, it is ignored. The resulting expression is as follows: 12${\mathrm{\gamma \Delta }}_{h^{\prime }x}\approx {\Delta }_{\mathrm{px}}-\left({\Delta }_{\mathrm{px}}-{\Delta }_{\mathrm{p2}}\right)\frac{h_x}{\Delta h}+\Delta P_{\mathrm{ax}}-{\gamma }_{\mathrm{gx}}{H}_o$

[0039] Therefore, Δh′x is as follows: 13$\begin{array}{cc}{\Delta }_{\mathrm{hx}}\cong \frac{\left[{\Delta }_{\mathrm{px}}+\left({\Delta }_{\mathrm{p2}}-{\Delta }_{\mathrm{px}}\right)\frac{\mathrm{hx}}{\Delta h}+\left(\Delta P_{\mathrm{ax}}-{\gamma }_{\mathrm{gx}}{H}_o\right)\right]}{\gamma }& \left(8\right)\end{array}$

[0040] The conventional method excludes the water specific gravity γ, but it doesn't compensate for errors γpx, γp2 of the pressure transducers and an error component due to an altitude difference. 14$\left(\Delta P_{\mathrm{ax}}-{\gamma }_{\mathrm{gx}}{H}_o\right)={\gamma }_{\mathrm{ao}}\left(H_o-h_x\right)-{\gamma }_{\mathrm{ao}}·\frac{P_{\mathrm{mx}}}{P_o}{H}_o={\gamma }_{\mathrm{ao}}{H}_o\left(1-\frac{P_{\mathrm{mx}}+P_o}{P_o}\right)-{\gamma }_{\mathrm{ao}}{h}_x$

[0041] For these reasons, the conventional method doesn't secure the higher accuracy of the water depth measurement. For example, if Ho=2000 Cm, hx=1000 Cm, Δh=100 Cm, t=−10° C., measuring the water depth hx (=1000 cm) is as follows: 15$\frac{{\Delta }_{\mathrm{p2}}}{\gamma }=-31\mathrm{cm},\frac{{\Delta }_{\mathrm{px}}}{\gamma }=-34\mathrm{cm},\frac{P_{\mathrm{ax}}-{\gamma }_{\mathrm{gx}}{H}_o}{\gamma }=-3.6\mathrm{cm}\therefore \Delta h_x^{\prime }=-34+\left(-31+34\right)\frac{10}{1}-36=-76\mathrm{cm}$

[0042] As a result, a water column pressure measuring error becomes larger, if a water depth in a reservoir, a lake, a river, etc. is measured by a water level meter, because an error of a pressure transducer that measures a bubble pressure is greatly changed dependent upon a circumferential air temperature t and a pressure P to be measured at an unattended observation post under the condition that an atmosphere temperature is in the range of −40° C. in Summer to +40° C. in Winter.

[0043] There often have occurred larger measuring errors of water depth or water level, due to it that a specific gravity of water gravity shifts at any time and an altitude difference between a water surface and a position of the transducer to be mounted, an altitude difference between lower and upper ends of a water column pressure measuring tube and a difference between a pressure of a bubble generating compressed gas to be measured and a pressure at a water depth to be measured based on a compressed density upon the bubble generating lies in a larger range.

[0044] As described above, the conventional technology has a disadvantage in that its water level measurement in a reservoir cannot secure the accuracy of ±1 cm contrary to the fact that the water level measurement of the reservoir requires the accuracy of less than ±1 cm.

[0045] For the purpose of explaining the invention in detail, firstly considering the property of a pressure transducer, the common characteristics of pressure transducers broadly well known are as follows:

[0046] The pressure transducer includes a temperature characteristic compensation device, because the property of a pressure sensor is changed according to a temperature. But, it generally doesn't compensate for the temperature. If a circumferential temperature is beside a normal one, for example the temperature 20˜25° C., the pressure measuring error is increased.

[0047] For example, the pressure transducers Model PTX1000 and PDCR882 having a relatively good performance has the following characteristics: when the rom temperature t=24° C., the Model PTX1000 represents that δ=±0.25% fs, the Model PDCR862 represents that δ=±0.1%fs. The usable temperature ranges of each of Models PTX1000 and PDCR862 are represented in table 1 as −40° C.˜+90° C. and −54° C.˜+125° C. Tables 1 and 2 represents results that a pressure error Δp is calculated into a water column pressure cmH2O, in which the pressure error Δp is a difference between a reference pressure based on an outputting signal, when the temperature t is 24° C., and a pressure based on an outputting signal, when the temperature t is changed in the range of −20° C.˜+50° C. the pressure difference Δp is equal to an absolute error of a water depth measurement, if a water specific gravity γ=1.0. Relative errors δ% are represented under the pressure error Δp. (Calculated data is rounded to the two decimals.)

TABLE 1
PTX1000 Δp cmH2O, δ%
	P: kg/cm
t° C.	0	0.1	0.2	0.4	0.6	0.8	1
(+50)	Δp	9	11	11.8	13.8	14.6	15.3	16
	δ%	—	11	5.9	3.45	2.43	1.9	1.6
(+40)	Δp	8	8.2	8.5	8.84	9.9	10.24	10.9
	δ	—	8.2	4.25	2.2	1.65	1.28	1.09
(+30)	Δp	3.5	5.1	5.43	6	4.26	3.8	3.4
		—	5.1	2.7	1.5	0.7	0.47	0.3
(+10)	Δp	−9	−6.7	−7.9	−9.6	−10.3	−11	−11.8
	δ	—	−6.7	−3.9	−2.4	−1.7	−1.4	−1.2
0	Δp	−17	−12.8	−16.3	−18	−19.3	−20.4	−21.2
	δ	—	−12.8	−8.2	−4.5	−3.2	−2.55	−2.1
−10	Δp	−26	−19.4	−23.5	−26.7	−28.7	−30.3	−31.7
	δ	—	−19.4	−11.7	−6.7	−4.7	−3.7	−3.17
−20	Δp	−34.5	−32.7	−35.2	−38.4	−42.2	−43.5	−45.4
	δ	—	−32	−17.6	−9.6	−7	−5.4	−4.5

[0048] Changing only the temperature of compressed air, but not cooling or freezing the pressure transducer on the whole, means that the pressure transducers are represented as better characteristics than those of Table 1. But, if the pressure transducer is used in a bubble water level meter, its whole temperature is changed according to a circumferential temperature. Table I represents characteristics when a pressure is increased from 0 to 1 kg/cm2, in which the characteristics are better than when a pressure is reduced from 1 to 0 kg/cm2. (It means Histeresis error or characteristics.)

[0049] Table 2 represents the characteristics of the pressure transducer Model PDCR826 that is better than those of Model PTX1000. Seeing Table 2, the temperature compensation error is little, when the temperature becomes +40, +30 or +10° C. But, the temperature compensation error is relatively larger, when the temperature t becomes 0, −10, −20 or +50° C. the Histeresis characteristics of the Model PDCR826 is inferior to that of Model PTX1000 the pressure error Δp is up to two times over Table 2, if the pressure is shifted from a high state to a low state. if the Histeriesis characteristics is bad, the pressure error becomes greater, when the water depth is small.

TABLE 2
PDCR826
	P: kg/cm2
t° C.	0	0.1	0.2	0.4	0.6	0.8	1
(+50)	Δp	1	1.08	1.15	1.32	1.59	1.83	1.98
	δ%	—	1.08	0.57	0.33	0.26	0.23	0.19
(+40)	Δp	1	0.3	0.31	0.43	1.11	1.01	0.78
	δ	—	0.3	0.155	0.108	0.186	0.127	0.08
(+30)	Δp	−1.5	0.017	0.097	0.044	0.035	0.06	0.34
	δ	—	0.017	0.048	0.011	0.006	0.0075	0.034
(+10)	Δp	−0.8	−0.92	−0.88	−0.0035	−0.085	−0.68	−0.94
	δ	—	−0.92	−0.44	˜0	−0.014	−0.085	−0.094
0	Δp	1.5	7.5	8.5	11.5	14.37	17.2	19.5
	δ	—	7.5	4.25	2.87	2.4	2.15	1.95
−10	Δp	1	1.5	1.03	1.56	1.9	2.2	2.6
	δ	—	1.5	0.52	0.4	0.32	0.275	0.26
−20	Δp	0.5	−1.84	−1.35	−0.5	˜0	1.04	2.07
	δ	—	−1.84	−0.67	−0.13	˜0	0.13	0.2

[0050] Seeing the temperature compensation characteristics, when a maximum pressure is applied, the pressure transducer is adjustable near to an error δ%fs at a normal temperature with a relative error being less. The less a measuring pressure, the greater an error δ is.

[0051] FIG. 3 shows a changing curve of a pressure error Δp when t=0° C. and −10° C., from which the following characteristics is known. When P=0, the output of Model PTX1000 becomes zero. The reason is because a temperature compensation circuit itself mounted in the pressure transducer represents a zero deviation according to the temperature change. If only the temperature of the compressed air is changed, but the temperature compensation circuit is not heated or cooled, the curve of Δp=f(t,P) is raised to about 26 cmH20, and the pressure error Δp becomes less, significantly.

[0052] But, the bubble type water level meter is heated or cooled on the whole in seasons, if it is mounted a water level observatory post that is not equipped with cooling and heating devices. Therefore, a water column pressure measuring error is changed due to the pressure error characteristics Δp=f(t;P) changed according to the temperature t and measuring pressure P in Tables 1 and 2. The curve of Δp=f(t;P) is divided into several sections to be represented as an approximate straight line in each section.

[0053] An object of the invention is to measure a water depth in an accuracy using water column pressure measuring tubes of n≧2 under the condition that a circumferential temperature t is changed in the range of −40° C.˜+50° C., even though the property of a pressure transducer is severely changed.

SUMMARY OF THE INVENTION

[0054] According to the invention, a water depth measuring method comprises steps of making a bundle of reference water column pressure measuring tubes of a U character form of the number of n≧2 at a same length with a water depth measuring tube in order to compensate for a total error ΣΔx of a water depth measurement; filling an amount of water into the reference water column pressure measuring tubes to form water columns ho1, ho2, . . . hon therein; measuring the water pressure in the reference water column pressure measuring tube by a pressure transducer before measuring a water depth; subtracting previously known water column pressure γhoi and γhoj derived from the measured water pressure values Pmi, and Pmj to calculate total errors ΣΔi and ΣΔj; calculating a total error ΣΔx based on the total errors ΣΔi and ΣΔj, in which the total error ΣΔx was generated upon measuring of the bubble generating pressure Pmx in the water column pressure measuring tubes; subtracting the total error ΣΔx from the bubble generating pressure Pmx, dividing the calculated water column pressures by the water column pressure γhoi(j) in the reference water column pressure measuring tube having most approximate value thereto and multiplying the re-calculated value by the water column pressure hoi(j) thereby to measure a water depth in a higher accuracy.

[0055] A bubble water depth measuring system comprises a compressed gas generator, a pressure transducer, a reference water column pressure measuring tube, a water depth measuring tube, an arithmetic logical transducer for calculating a water depth and a drive controller for supplying/interrupting compressed gas.

[0056] The water depth measuring tube has a length corresponding to a water depth. The reference water column pressure measuring tubes of n≧2 of a U shape has the same inner diameter as that of the water depth measuring tube coupled therewith. A diaphragm is connected through an electromagnetic valve to one tube of the reference water column pressure measuring tube to adjust the speed of the compressed gas from a buffer tank to be supplied to the reference water column pressure measuring tubes. A transparent container is connected to the other tube of the reference water column pressure measuring tubes to measure an amount of water therein and check it, periodically, and includes a nipper for supplementing water reduced due to the evaporation. A tube is directly mounted above the transparent container and has the same inner diameter as that of the reference water column pressure measuring tubes and a length portion, on the surface of which scales are formed. A thin film tube is made of rubber and directed to the upper portion of the tube to be expanded by a gas pressure, and a manual valve connected to the upper end of the thin film tube.

BRIEF DESCRIPTION OF THE DRAWINGS

[0057] The invention now will be described in detail with the reference to the accompanying, in which:

[0058] FIG. 1 is a view illustrating an error measured according to a bubble type water depth measuring method;

[0059] FIG. 2. is a view illustrating an apparatus for performing a conventional water depth measuring method;

[0060] FIG. 3 is a view illustrating error-changing curves dependent upon the change of a temperature and a measuring pressure of a pressure transducer;

[0061] FIG. 4 is a vies illustrating a bubble type water depth measuring method according to the invention;

[0062] FIG. 5 is a view illustrating the transformation of a curve of an error Δp of a pressure transducer into a straight line;

[0063] FIG. 6 is a view illustrating a configuration of an apparatus for performing a bubble type water depth, water level measuring method according to the invention;

[0064] FIG. 7 is a view illustrating the closure of a reference water depth measuring tube according to the invention; and,

[0065] FIG. 8 is a view illustrating an angle error of mounting a reference water depth measuring tube according to the invention.

DETAILED DESCRIPTION OF THE INVENTION

[0066] Referring to FIG. 4, only two water pressure transducers (n=2) are shown for simply explaining the invention. 1 is a water depth measuring tube, 21 and 22 are reference water pressure measuring tubes that are connected through valves 51, 52 and 53 to a buffer tank 4. The buffer tank 4 is connected to a pressure transducer 3. Pressure gas Pg is supplied through a valve 50 to the buffer tank 4. The reference water depth measuring tubes 21 and 22 are made as U shape, each of one portion of which is connected to the valves 52 and 53 and other portion of which are connected in turn to the buffer tank 4 and the pressure transducer 3 with the other end being opened. The water pressure measuring tube, 21 and 22 are filled with river water. As compressed gas is supplied to the water pressure measuring tubes for a predetermined short time, the water previously filled is pushed in a direction contrary to the gas supply. At this time, the reference water pressure measuring tubes each is kept at water column h01 and h02 (h01<h02). The water column h01 and h02 are selected according to a measuring range of a water depth hx and the property of a pressure transducer.

[0067] An algorithms of measuring the water depth hx is as follows:

[0068] When a bubble is generated at the water depth measuring tube 1 with the valves 52 and 53 being closed and the compressed gas being supplied to the buffer tank 4, the valve 50 is closed. At this time, the gas pressure Pmx is measured by the pressure transducer 3 and then memorized in a memory of a system.

Pmx=γhx+ΣΔx   (9)

[0069] As the valve 51 is closed and then the valve 53 is opened and closed for a short time, water filled in the reference water pressure tubes is pushed up along the other tube thereof. At the this time, the water column h02 is maintained, and the pressure Pm2 is measured and then memorized.

Pm 2 =γh 02 +ΣΔ2  (10)

[0070] Herein, the water column h02 is previously known. In the expression (10), ΣΔ′2 is memorized with γ′h02 being subtracted from Pm2.

ΣΔ′2=Pm2−γ′h02  (11)

[0071] ΣΔ′2 may be not equal to ΣΔ2, because an average specific gravity γ′ is used with the specific gravity γ of the water column h02 being unknown.

[0072] With the valve 53 being closed and the valve 52 being opened for a short time and closed, the water column pressure γh01 is measured.

Pm 1 =γh 01 +ΣΔ1  (12)

[0073] A result that γ′h01 is subtracted from Pm1 is memorized.

ΣΔ′1=Pm1−γ′h01  (13)

[0074] (Of course, the measuring turn of Pm1 and Pm2 is changeable.)

[0075] In the expression (9), ΣΔx is calculated as follows: 16$\begin{array}{cc}{\mathrm{\Sigma \Delta }}_x^{\prime }={\mathrm{\Sigma \Delta }}_1^{\prime }+\frac{{\Sigma }_2^{\prime }-{\mathrm{\Sigma \Delta }}_1^{\prime }}{P_{\mathrm{m2}}-P_{\mathrm{m1}}}\left(P_{\mathrm{mx}}-P_{\mathrm{m1}}\right)& \left(14\right)\end{array}$

[0076] The water depth hx is calculated as follows: 17$\begin{array}{cc}{h}_x^{\prime }=\frac{P_{\mathrm{mx}}-{\mathrm{\Sigma \Delta }}_x^{\prime }}{P_{\mathrm{m2}}-{\mathrm{\Sigma \Delta }}_2^{\prime }}·h_{02}=\frac{\gamma h_x+{\mathrm{\Sigma \Delta }}_x-{\mathrm{\Sigma \Delta }}_x^{\prime }}{\gamma h_{01}+{\mathrm{\Sigma \Delta }}_2-{\mathrm{\Sigma \Delta }}_2^{\prime }}{h}_{02}=h_x\left(\frac{1+\frac{{\mathrm{\Sigma \Delta }}_x-{\mathrm{\Sigma \Delta }}_x^{\prime }}{\gamma h_x}}{1+\frac{{\mathrm{\Sigma \Delta }}_2-{\mathrm{\Sigma \Delta }}_2^{\prime }}{\gamma h_{02}}}\right)\mathrm{or}& \left(15\right)\\ h_x^{\prime }=\frac{P_{\mathrm{mx}}-{\mathrm{\Sigma \Delta }}_x^{\prime }}{P_{\mathrm{m1}}-{\mathrm{\Sigma \Delta }}_1^{\prime }}{h}_{01}=h_x\left[\frac{1+\frac{{\mathrm{\Sigma \Delta }}_x-{\mathrm{\Sigma \Delta }}_x^{\prime }}{\gamma h_x}}{1+\frac{{\mathrm{\Sigma \Delta }}_1-{\mathrm{\Sigma \Delta }}_1^{\prime }}{\gamma h_{01}}}\right]& \left(16\right)\end{array}$

[0077] The expression (15) or (16) is selected as follows: if Pmx is close to Pm2, the expression (15) is selected, and if Pmx is closer to Pm1, the expression (16) is selected.

[0078] Assuming that 18$\frac{{\Delta }_2}{\gamma h_{02}}<<1.0,$

[0079] the expression (1) is as follows: 19$\begin{array}{cc}{h}_x^{\prime }=h_x\left(1+\frac{{\Delta }_x}{\gamma h_x}\right)\left(1-\frac{{\Delta }_2}{\gamma h_{02}}\right)\approx h_x\left(1+\frac{{\Delta }_x}{\gamma h_x}-\frac{{\Delta }_2}{\gamma h_{02}}\right)& \left(17\right)\end{array}$

[0080] Herein, 20$\frac{{\Delta }_x·{\Delta }_2}{{\gamma }^2·h_{02}·h_x}$

[0081] is ignored, because it is a very little value.

[0082] A relative error δhx of h′x resulted from the water depth measurement is as follows: 21$\begin{array}{cc}{\delta }_{\mathrm{hx}}=\frac{h_x^{\prime }}{h_x}-1=\frac{{\Delta }_x}{\gamma h_x}-\frac{{\Delta }_2}{\gamma h_{02}},& \left(18\right)\end{array}$

[0083] An absolute error Δhx is as follows: 22$\begin{array}{cc}{\Delta }_{\mathrm{hx}}={\delta }_{\mathrm{hx}}·h_x=\frac{{\Delta }_x}{\gamma }-\frac{{\Delta }_2{h}_x}{\gamma h_{02}}& \left(19\right)\end{array}$

[0084] Δx and Δ2 are as follows: Δx=ΣΔx−ΣΔ′x, in which ΣΔx is a measuring absolute error. The error is an approximate one that a curve of Δp=f(t;P) is converted into a straight line. Therefore, the larger the number of sections that the curve Δp=f(t;P) is formed into the straight line is, the less Δx is. 23${\Delta }_2={\mathrm{\Sigma \Delta }}_2-{\mathrm{\Sigma \Delta }}_2^{\prime }=\left(P_{\mathrm{m2}}-\gamma h_{02}\right)-\left(P_{\mathrm{m2}}-\gamma h_{02}\right)=\gamma h_{02}\left(\frac{{\gamma }^{\prime }}{\gamma }-1\right).$

[0085] Therefore, if γ′=γ, Δ2=0.

[0086] Assuming that the water specific gravity γ=0.0998 gm/cm3, but γ′=1.0 gm/cm3 in the section of h02 and h02=500 m, 24${\Delta }_2=0.9998·500\left(\frac{1}{0.9998}-1\right)=0.099gm\text{/}{\mathrm{cm}}^2$

[0087] At that time, the water depth error is about 0.1 cm (1 mm). When the expression (19) is derived, the water columns h02 and h02 were assumed to be not changed at the reference water column pressure measuring tubes 21 and 22. But, as time passed, water is escaped off in vapor. If the water columns h01 and h02 each becomes h01−Δh01 and h02−Δh02, Δh′x is increased more, to which will be explained below. Also, an evaporation-preventing device will be explained below.

[0088] First, according to a method of the invention, when the water depth hx were measured using an pressure transducer Model PTX1000, a water depth measuring absolute error Δh′x (cm) is represented in Table 3. The measuring conditions are as follows:

[0089] When t=−10° C. in a water depth of 200˜600 cm, the curve of Δp=f(t;P) is represented as a straight line converted. the water depth is selected so that h01=200 cm and h02=600 cm. The measuring range of the water depth hx is 200˜600 cm. Assuming that the water specific gravity γ=0.0998 gm/cm3, but γ′=1.0 gm/cm3, Ho=20m (2000 cm), in which the compressed gas is a compressed air

TABLE 3
hx cm	200	300	400	500	600	800	1000
h′x cm	199,999	299.65	399.28	499.79	599.94	800.8	1002.48
Δhx cm	˜0	−0.35	−0.72	−0.204	−0.06	+0.8	+2.48

[0090] When Table 3 is written, Δhx is calculated as follows: if the water depth hx of 500 cm is measured, as shown in Table 1, h01=200 and h02=600, when t=−10° C. And, when hx=500 cm, the errors of the pressure transducer Model PTX1000 are as follows: Δp1=−23.5, Δp2=−28.7, Δpx=−27.6 cm H20. Under these conditions, according to the algorithm of the invention, h′x and the measuring result of γhx are as follows:

P mx =γh x +ΣΔx=γhx+Δpx+ΔPα−γgHo=0.9998·500−27.6+1.2·10−3(2000-500)−1.2·10−3·1.5·2000=470.5 cm H2O

[0091] The measuring results of γh01 and γh02 are follows: 25$P_{\mathrm{m1}}=\gamma h_{01}+{\mathrm{\Sigma \Delta }}_1=0.9998·200-23.5+1.2·{10}^3\left(2000-200\right)-1.2·{10}^{-3}·1.2·2000=175.74\mathrm{cm}{H}_{2}O$$P_{\mathrm{m2}}=\gamma h_{02}+{\mathrm{\Sigma \Delta }}_2=0.9998·600-28.7+1.2·{10}^{-3}\left(2000-600\right)-1.2·{10}^{-3}·1.6·2000=569.02\mathrm{cm}{H}_{2}O$${\mathrm{\Sigma \Delta }}_1^{\prime }=P_{\mathrm{m1}}-\gamma h_{01}=175.74-1.0·200=-24.26$${\mathrm{\Sigma \Delta }}_2^{\prime }=P_{\mathrm{m2}}-\gamma h_{02}=569.02-1.0·600=-30.98$${\mathrm{\Sigma \Delta }}_x^{\prime }={\mathrm{\Sigma \Delta }}_1^{\prime }+\frac{{\mathrm{\Sigma \Delta }}_2^{\prime }-{\mathrm{\Sigma \Delta }}_1^{\prime }}{P_{\mathrm{m2}}-P_{\mathrm{m1}}}·\left(P_{\mathrm{mx}}-P_{\mathrm{m1}}\right)=-24.26+\frac{-30.98+24.26}{569.02-175.54}\left(470.5-175.54\right)=-29.296\therefore h_x^{\prime }=\frac{P_{\mathrm{mx}}-{\mathrm{\Sigma \Delta }}_x^{\prime }}{P_{\mathrm{m2}}-{\mathrm{\Sigma \Delta }}_2^{\prime }}{h}_{02}=\frac{470.5+29.296}{569.02+30.98}·600=499.796\mathrm{cm}$$\Delta h_x^{\prime }=499.796-500=-0.204\mathrm{cm}$

[0092] Herein, it is known that when hx=500 cm, a measuring error Δh′x=−2 mm.

[0093] In Table 3, the obsolute error Δh′x is represented as a negative “−” sign. When hx>600 cm, Δh′x becomes larger, which is represented as a positive “+” sign. When h=1000 cm, Δh′x=+2.48 cm. The cause is follows:

[0094] Referring to FIG. 5, when t=−10° C., a changing curve of Δpi is exaggeratedly shown. It is represented in a straight line as follows: 26${\mathrm{\Sigma \Delta }}_x^{\prime }={\mathrm{\Sigma \Delta }}_1^{\prime }+\frac{{\mathrm{\Sigma \Delta }}_1^{\prime }-{\mathrm{\Sigma \Delta }}_1^{\prime }}{P_{\mathrm{m2}}-P_{\mathrm{m1}}}\left(P_{\mathrm{mx}}-P_{\mathrm{m1}}\right)$

[0095] If the curve of Δp is approximated into a straight line II (Dotted line) in a section of hx that is a range of 600˜1000 cm, the error Δh′x may be significantly reduced. Unlike FIG. 4, the reference water column pressure measuring tubes of three (n=3) are selected to have water column pressures h01, h02 and h03, respectively. If Pmx>600, ΣΔ′x is calculated as follows: 27${\mathrm{\Sigma \Delta }}_x^{\prime }={\mathrm{\Sigma \Delta }}_2^{\prime }+\frac{{\mathrm{\Sigma \Delta }}_3^{\prime }-{\mathrm{\Sigma \Delta }}_2^{\prime }}{P_{\mathrm{m3}}-P_{\mathrm{m2}}}\left(P_{\mathrm{mx}}-P_{\mathrm{m2}}\right)$

[0096] If n≧3, the expressions (14) and (15) are represented as follows: 28$\begin{array}{cc}{\mathrm{\Sigma \Delta }}_x^{\prime }={\mathrm{\Sigma \Delta }}_i^{\prime }+\frac{{\mathrm{\Sigma \Delta }}_j^{\prime }-{\mathrm{\Sigma \Delta }}_i^{\prime }}{P_{\mathrm{mj}}-P_{\mathrm{mi}}}\left(P_{\mathrm{mx}}-P_{\mathrm{mi}}\right)& \left(20\right)\\ h_x^{\prime }=\frac{P_{\mathrm{mx}}-{\mathrm{\Sigma \Delta }}_x^{\prime }}{P_{\mathrm{mj}}-{\mathrm{\Sigma \Delta }}_j^{\prime }}{h}_{\mathrm{oj}}\mathrm{or}& \left(21\right)\\ h_x^{\prime }=\frac{P_{\mathrm{mx}}-{\mathrm{\Sigma \Delta }}_x^{\prime }}{P_{\mathrm{mi}}-{\mathrm{\Sigma \Delta }}_i^{\prime }}{h}_{\mathrm{oi}}& \left(22\right)\end{array}$

[0097] Like this, it is noted that the invention can secure a higher accuracy because the total error of the water column pressure is compensated at once.

[0098] The error Δh′x represented in Table 3 is a result of assuming that h01=const and h02=const. But, if water filled in a reference water column pressure measuring tube 2i is escaped off as a time passed, hoi is reduced and becomes hoi−Δh. In other words, h′oi=hoi−Δh. If hoi is substituted into the measurement expressions of ΣΔi and h′x, an error Δh′x is increased. In this case, an error Δh″x caused due to Δh will be obtained by the reference of the expressions (20) and (21).

[0099] All reference water column pressure-measuring tubes are in water, which means the temperature conditions are same. Therefore, water columns Δh reduced in all the reference measuring tubes are same to one another. ΣΔ″i, ΣΔ″j, ΣΔ″x and Δh″x are as follows: 29$\begin{array}{cc}{\mathrm{\Sigma \Delta }}_i^{″}={\mathrm{Pm}}_i-\gamma h_{\mathrm{oi}}=\gamma \left(h_{\mathrm{oi}}-\Delta h\right)+{\mathrm{\Sigma \Delta }}_i-\gamma h_{\mathrm{oi}}={\mathrm{\Sigma \Delta }}_i-\mathrm{\gamma \Delta }h{\mathrm{\Sigma \Delta }}_j^{″}={\mathrm{Pm}}_j-\gamma h_{\mathrm{oj}}={\mathrm{\Sigma \Delta }}_j-\mathrm{\gamma \Delta }h\begin{array}{c}{\mathrm{\Sigma \Delta }}_x^{″}={\mathrm{\Sigma \Delta }}_i^{″}+\frac{{\mathrm{\Sigma \Delta }}_j^{″}-{\mathrm{\Sigma \Delta }}_i^{″}}{P_{\mathrm{mj}}-P_{\mathrm{mi}}}\left(P_{\mathrm{mx}}-P_{\mathrm{mi}}\right)\\ =\left({\mathrm{\Sigma \Delta }}_i-\mathrm{\gamma \Delta }h\right)+\frac{{\mathrm{\Sigma \Delta }}_j-\mathrm{\gamma \Delta }h-{\mathrm{\Sigma \Delta }}_i+\mathrm{\gamma \Delta }h}{\gamma h_{\mathrm{oj}}+{\mathrm{\Sigma \Delta }}_j-\gamma h_{\mathrm{oi}}-{\mathrm{\Sigma \Delta }}_i}\left[P_{\mathrm{mx}}-\left(\gamma h_{\mathrm{oi}}+{\mathrm{\Sigma \Delta }}_i\right)\right]\\ =\left({\mathrm{\Sigma \Delta }}_i-\mathrm{\gamma \Delta }h\right)+\frac{{\mathrm{\Sigma \Delta }}_j-{\mathrm{\Sigma \Delta }}_i}{\gamma (h_{\mathrm{oj}}-h_{\mathrm{oi}})+{\mathrm{\Sigma \Delta }}_j-{\mathrm{\Sigma \Delta }}_i}\left[P_{\mathrm{mx}}-\left(\gamma h_{\mathrm{oi}}+{\mathrm{\Sigma \Delta }}_i\right)\right]\end{array}& \left(23\right)\end{array}$

[0100] Δ″x is changed only by −γΔh. 30$\begin{array}{cc}\begin{array}{c}{h}_x^{″}=\frac{P_{\mathrm{mx}}-{\mathrm{\Sigma \Delta }}_x^{″}}{P_{\mathrm{mj}}-{\mathrm{\Sigma \Delta }}_j^{″}}{h}_{\mathrm{oj}}=\frac{\gamma h_x+{\mathrm{\Sigma \Delta }}_x-{\mathrm{\Sigma \Delta }}_x^{″}}{\left(\gamma h_{\mathrm{oj}}-\mathrm{\gamma \Delta }h\right)+{\mathrm{\Sigma \Delta }}_j-{\mathrm{\Sigma \Delta }}_j^{″}}{h}_{\mathrm{oj}}\\ =\frac{\gamma h_x+{\mathrm{\Sigma \Delta }}_x-\left({\mathrm{\Sigma \Delta }}_x-\mathrm{\gamma \Delta }h\right)}{\gamma h_{\mathrm{oj}}-\mathrm{\gamma \Delta }h+{\mathrm{\Sigma \Delta }}_j-\left({\mathrm{\Sigma \Delta }}_j-\mathrm{\gamma \Delta }h\right)}{h}_{\mathrm{oj}}\\ =\frac{\gamma h_x\left(1+\frac{\Delta h}{h_x}\right)}{\gamma h_{\mathrm{oj}}}{h}_{\mathrm{oj}}=h_x\left(1+\frac{\Delta h}{h_x}\right)\end{array}& \left(24\right)\\ \delta h_x^{″}=\frac{\Delta h}{h_x}& \left(25\right)\end{array}$

[0101] Comparing the expression (25) with the expression (18), only 31$\frac{{\Delta }_x}{\gamma h_x}$

[0102] is replaced by 32$\frac{\Delta h}{\gamma h_x},$

[0103] the absolute error as follows:

Δhx″=Δh  (26)

[0104] If the water column hoi in the reference water column pressure measuring tube is reduced by Δh=5 mm due to water evaporation, Δh′x increases upto +5 mm before evaporation. If Δh′x is positive “+”, the error will be increased more +5 mm. On the contrary, if Δh′x is negative “−”, the error will be reduced less +5 mm. For these reasons, there need countermeasures for preventing the evaporation of water filled in the reference water column pressure-measuring tubes, checking its water column hoi and supplementing the evaporated water. The invention is designed to comply with the countermeasures.

[0105] Another problem is that upon the calculation of ΣΔI and ΣΔj γ′hoi and γ′hoj are subtracted from Pmi and Pmj using an average value γ′ because the water specific gravity γ is unknown. The water depth measuring error δhx(γ) happens because of the error of δγ=(γ′−γ)γ resulting from the water specific gravity γ. 33$\begin{array}{cc}{\delta }_{h^{\prime }x\left(\gamma \right)}\approx -{\delta }_{\gamma }\frac{h_{\mathrm{oi}\left(j\right)}}{h_x}& \left(27\right)\end{array}$

[0106] A water temperature of reservoir is usually 3˜4° C. after the water surface is frozen in winter. An average temperature of a water depth doesn't exceed 22° C. The temperature of the water surface is raised up to 24° C. If the water is clean, its density is changed in the range of 0.999992 to 0.997802 gm/cm3. But, the concentration of floating materials is usually 0.5 gm/l (=5·10−4 gm/cm3). When the water temperature reaches around 22° C., the water specific gravity γ=0.997802+5·10−4=0.998302 gm/cm3. In winter, the concentration of floating particles is reduced, significantly. If the average density is used irrelevant to seasons, δγ is changed in the range of −8.4·10−4 to +8.5·10−4.

[0107] If hoi=200 cm and hx=300 cm in the expression (27), δHx(γ)=|5.66·10−4|, and Δhx=0.169≈0.17 cm. If hx=600 cm, Δhx=0.11 cm. Therefore, the influence of γ′ is very small. For it, it is most important to take a measure to prevent the evaporation of water filled in the water column pressure measuring tube 2i.

[0108] An effective and simple method of the evaporation prevention is to store water in a closed space. The closed space is maintained at a relative humidity of 100%. Only, the relative humidity is changed by an absolute amount until reaching 100%. Another liquid can be filled in the water column pressure measuring tube instead of water. Herein, it is noted that the specific gravity change according to a temperature of the liquid is exactly confirmed.

[0109] Referring to FIGS. 6 and 7, an apparatus of performing a water depth, water level measuring method is shown according to the invention.

[0110] As shown in FIG. 6, 1 is a water depth measuring tube, 21, 22 and 23 are reference-measuring tubes of U shape that forms water columns of h01<h02<h03. 3 is a pressure transducer, and 4 is a buffer tank which supplies compressed gas to a plurality of measuring tubes. 10 is a compressed gas generator, for which a micro compressor or a compressed nitrogen tank can be used. 50, 51, 52, 53 and 54 is a diaphragm for adjusting an amount of compressed gas to be introduced into the reference measuring tubes, but instead of it a valve having a small inner diameter can be used. 71, 72 and 73 are transparent containers for measuring an amount of water filled in the reference measuring tubes 21, 22 and 23. 81, 82 and 83 are a corrugated tube made of a thin rubber film or another soft material. 91, 92 and 93 are manual valves. The detailed configuration of the transparent container 7i and the corrugated tube 8i is shown in FIG. 8. 11 is a digital converter for converting an outputting signal from the pressure transducer 3 into an unit of mmH2O. 12 is an arithmetic logical controller or microprocessor to calculate a water depth according to the invention, which includes devices for forming outputting signals such as code, analog signals, etc. 13 is a drive controller for controlling the operating of the valves Si and the compressed gas generators 10 and 11. Of course, the digital converter 11, the arithmetic logical controller 12 and the drive controller 13 can be integrated with each another as one integrated circuit chip.

[0111] 2 i is closed to prevent the evaporation of water filled therein. The valves 50 and 51 are opened in turn based on the outputting signal from the drive controller 13, and the other valves 52, 53 and 54 are closed. Then, as the compressed gas is supplied to the water depth measuring tube 1, the water filled in the water depth measuring tube 1 is gotten out downward and bubbles are generated. At this time, the valve 50 is closed. When the outputting signal from the pressure transducer 3 is stabilized, the water column Pmx is measured, which is stored in the arithmetic logical controller 12, and the valve 51 is closed.

[0112] The arithmetic logical controller 12 judges on which section of γh01˜γh02 and γh02˜γh03 the water column pressure Pmx is corresponded to and outputs its control signal to the drive controller 13. For example, if the water column pressure Pmx is corresponded to the section of γh02˜γh03, the valve 53 is opened for a predetermined short time, the compressed gas is supplied to the reference measuring tube 22 to push up the water therein along its left tube and form the water column h02. at this time, the water column pressure Pmx is measured. The arithmetic logical controller 12 calculates a water depth hx according to a water depth-measuring algorithm of the invention.

[0113] The characterized portion of the system is a configuration of the reference measuring tube 2i. Referring to FIG. 7, the transparent container 7i, the corrugated tube 8i and the manual valve 9i are enlarged. The transparent container 7i measures an amount of water to be injected into the reference measuring tube 2i and also is used in checking an amount of water evaporated and supplementing corresponding amount of water in the reference measuring tube 2i. 14i is a tube having the same inner diameter as that of the reference measuring tube 2i, on the outer wall of which scales are formed in an unit of mm. 15i is a nipper such as an injector which is used in injecting water into the transparent container 7i or supplementing water upon evaporating. The tube 14i has a length of l that ≈10 cm. The transparent container 7i has an inner volume vi as follows: 34$\begin{array}{cc}{v}_i=\frac{\pi d^2}{4}\left(h_{\mathrm{oi}}-\frac{l}{2}\right)& \left(28\right)\end{array}$

[0114] If the water column hoi is subject to being formed in the reference measuring tube 2i, the water volume 35$v_{\mathrm{io}}=\frac{\pi {\stackrel{\_}{d}}^2}{4}{h}_{\mathrm{oi}}.$

[0115] {overscore (d)} is an inner diameter of the tube. When such like an amount of water is filled in the transparent container 7i, the water is fully filled in the transparent container 7i, and the remaining water is raised up to a middle position (˜l/2) of the tube 14i. Therefore, an exactly water volume v, can be confirmed through the scales of the tube 14i.

[0116] After a certain time period passed, the water column hoi is reduced due to the evaporation. At this time, the water column error Δhoi is checked as follows: the compressed gas continues to be supplied through the valve 5i to the reference measuring tube 2i. The water in the reference measuring tube 2i is filled in the transparent container 7i, and the scale of the tube 14i is confirmed to check Δhoi. Then, additional water is supplemented using the nipper 15i. In this case, the manual valve 9i may be opened.

[0117] The tube 8i is expanded upon the air injection thereinto, the volume of which is sufficient. Even if air of 500 cm3 is injected, the tube 8i is easily inflated.

[0118] The tube 2i for a time period τ, the water is raised up along the left portion of the reference measuring tube 2i to form the water column hoi. The diaphragm 6i acts to adjusting an amount of compressed gas to be injected or injecting speed, a hole diameter of which is smaller that the inner one of the reference measuring tube 2i. The diaphragm 6i is not necessary to be separately made, but instead of it a manual adjustable valve is available.

[0119] At this time, air filled in the left portion of the reference measuring tube 2i is injected into the transparent container 7i and the tube 8i. The air volume is about vio. If the transparent container 7i is closed, the air pressure is increased. But, as the tube 8i is expanded, the air pressure becomes smaller. The water column pressure P′mi of γhoi is as follows:

P′ mi =γh oi +ΣΔ1+Δpc=Pm+Δpc;  (29)

[0120] Wherein, Δpc is a pressure required for expanding the tube 8i. The Δpc becomes smaller, if its material is soft and its size is larger. The Δpc is previously measured and stored in the arithmetic logical controller 12.

[0121] Next, the valve 9i is opened to measure the water column pressure Pmi, and then the valve 9i is closed to measure a water column pressure P′mi. The Δpc is as follows:

Δpc=P′m−Pmi  (30)

[0122] Of course, the tube expanding pressure Δpc is changed according to the temperature of air in the tube 8i and the position where the water column hoi is raised upward or fell downward an appointed position. In order to prevent the increasing of the water depth measuring error, the change of the tube expanding pressure Δpc must be secured so that Δpc<<ΣΔ. To it, the tube 8i is made of a very thin film and has a sufficient volume. For example, a balloon can be preferably used. Herein, most important condition is as follows:

Δpc=Δpc1=Δpc2=Δpc3= . . . =Δpcn  (31)

[0123] If the condition is secured, the errors of the tube expanding pressure Δpc are offset to one another. The condition of the expression (31) is secured by adjusting the size of the tube 8i. For example, the position where the valve 9i is mounted on the tube 8i is adjustable. Under the condition that the tube expanding pressure Δpc is measured and memorized, the water depth is measured according to the algorithms of the invention. 36$\begin{array}{cc}\Sigma {\Delta }_i=P_{\mathrm{mi}}^{\prime }-\left(\gamma h_{\mathrm{oi}}+\Delta \mathrm{pc}\right);\Sigma {\Delta }_j=P_{\mathrm{mj}}^{\prime }-\left(\gamma h_{\mathrm{oj}}+\Delta \mathrm{pc}\right);\Sigma {\Delta }_x^{\prime }\begin{array}{c}=\Sigma {\Delta }_i+\Delta \Delta \mathrm{pc}+\frac{\Sigma {\Delta }_j+\Delta \Delta \mathrm{pc}-\Sigma {\Delta }_i-\Delta \Delta \mathrm{pc}}{P_{\mathrm{mj}}-P_{\mathrm{mi}}}\left(P_{\mathrm{mx}}-P_{\mathrm{mi}}\right)\\ =\Sigma {\Delta }_i+\Delta \Delta \mathrm{pc}+\frac{\Sigma {\Delta }_j-\Sigma {\Delta }_i}{P_{\mathrm{mj}}-P_{\mathrm{mi}}}\left(P_{\mathrm{mx}}-P_{\mathrm{mx}}\right)\end{array}& \left(32\right)\end{array}$

[0124] Wherein, ΔΔpc is a difference between Ape previously memorized and Δ′pc caused upon measuring; ΔΔpc=Δ′pc−Δpc. Therefore, ΣΔ′x has the difference by Δpc compared with ΣΔx, but if ΔΔpc is smaller than the error of ΣΔx, it is ignoble.

[0125] As described above, the reference measuring tubes 21, 22 and 23 are almost closed not to evaporate the water filled therefrom. Especially, an amount of evaporated water is much more reduced in winter because the water temperature is lower. If the circumferential temperature is below −5° C., an amount of water can be checked without using the transparent tube 7i, because the possibility of freezing water is very high. Therefore, if an amount of water is checked and supplemented before the winter starts, it is enough. In summer, it is enough that an amount of water vi filled in the reference measuring tube 2i is checked.

[0126] The water depth measuring tube 1 and the reference measuring tube 2i are made of materials to which water is not adhesive. Cheapest thing is a polyurethane tube, and the more preferable thing is Teflon tube. The suitable inner diameter d of the measuring tubes is 2˜3 mm. All measuring tubes are coupled with each another in a bundle.

[0127] According to the invention, a system is a little complex, but the accuracy of the water depth measurement is very higher independent of the weather change. It is also not necessary to use an expensive pressure transducer. The invention has an advantage in that a pressure transducer not only having an inferior property but also being cheaper can be used.

[0128] On the other hand, when the bubble water level measuring system is mounted, there are cases that the measuring tubes 1 and 2i are mounted vertically and along the slope of a reservoir bank.

[0129] As shown in FIG. 8A, even though the measuring tube is vertically mounted, it is slanted at an angle α. The water column pressure γhoi is as follows:

γh′oi cos α  (33)

[0130] Wherein, h′oi is a length of a water column in the reference measuring tube 2i slanted at the angle α, and γhoi is a water column pressure of h′oi.

[0131] As shown in FIG. 8B, if a bundle of the measuring tubes are mounted along the slope of the reservoir bank, γhoi is as follows:

γhoi=γh′oi sin {overscore (β)}  (34)

[0132] Wherein, {overscore (β)} is an average gradient angle h′oi is exactly confirmable by an amount of water vi filled in the reference measuring tube 2i, but α and {overscore (β)} can't be exactly measured. If an angle error a of mounting the measuring tube is 1˜3°, a water depth measuring supplementary error is −0.03˜−0.3 cm. Unless the mounting angle error exceeds 2°, the supplementary error is ignorable. But, if the measuring tube is mounted at the gradient angle β and the gradient angle β is exactly not found, the measuring error of the water depth hx becomes larger. But, there happen cases that it is not possible to measure the gradient angle β. The result of measuring the water column pressure γhoi of the reference measuring tube 2i is as follows:

Pmi=γhoi+ΣΔi=γh′oi sin β+ΣΔi

[0133] Herein, what is exactly known is only h′oi. Therefore, in order to find the γhoi the ΣΔI gets calculated, exactly. For preparing these cases, the curve or table of Δp=f(t;P) is written in advance in a manner that a pressure transducer having a higher accuracy is corrected and checked. The water column pressure Pmi is measured by the pressure transducer, Δpa−γgHo is calculated and γhoi is calculated as follows:

γhoi=Pmi−(Δpi+ΔPa−γgHo)  (35)

[0134] As described above, according to the invention, reference water column pressure measuring tubes of n≧2 in an U shape are disposed at the same length along with a water depth measuring tube, an amount of water is filled in the reference water column pressure measuring tubes, compressed gas is supplied to the reference water column pressure measuring tubes from one side tube thereof for a short time period and the water is pushed up in the other tube to form the water column pressures γh01, γh02 . . . γh0n and measure the water depth.

[0135] Accordingly, the invention is mostly used in using the measurement of the water depth, water level in a reservoir, a lake and an underground water and can measure the water depth, water level in a higher accuracy compensating for errors caused due to the circumferential environment such as a temperature, an altitude difference, a water specific gravity, etc. at once.

Claims

1. A bubble water depth measuring method of finding a water column pressure of a water depth measuring tube with a bubble pressure to measure a water depth comprising steps of: disposing reference water column pressure measuring tubes of a U shape of the number of n≧2 at a same length with the water depth measuring tube; filling an amount of water into the reference water column pressure measuring tubes ho1, ho2, . . . hon to form a water column therein; supplying compressed gas to the reference water column pressure measuring tubes for a short time period from one tube thereof and raising up along the other tube thereof to form water column pressures γho1, γho2, . . . γhon; memorizing the water column pressure Pmx that the water column pressures of the water depth measuring tube are measured with the bubble generating pressure; measuring gas pressures Pm1 and Pm2 for maintaining the water column pressures γho1 and γho2 at the reference water column pressure measuring tubes corresponding to sections hoi˜hoj positioned under the water column pressure Pmx; calculating total errors ΣΔi and ΣΔj in the water column pressures; calculating a total error ΣΔx based on the total errors ΣΔi and ΣΔj as follows; and, 37Σ⁢ ⁢Δi=Pmi-γ⁢ ⁢hoi,Σ⁢ ⁢Δj=Pmj-γ⁢ ⁢hojΣ⁢ ⁢Δx=Σ⁢ ⁢Δi+Σ⁢ ⁢Δj-Σ⁢ ⁢ΔiPmj-Pmi⁢(Pmx-Pmi)calculating a water depth hx as follows; 38hx=Pmx-Σ⁢ ⁢ΔxPmi⁡(j)-Σ⁢ ⁢Δi⁡(j)·hoi⁡(j);wherein, Pmi(j) and Hoi(j) mean Pmi or Pmj and hoi or hoj. Among them the water column pressure closest to Pmx is selected.

2. The bubble water depth measuring method as claimed in claim 1, in which: the step of measuring the water column pressure Pmx with the gas pressures in the range of hxmin to hxmax that the water depth hx is measured based on the number n of the reference water column pressure measuring tube and the water column pressures γho1, γho2, . . . γhon further comprise steps of dividing a curve Δp=f(t;P) of an absolute error changed according to a pressure P that is used in measuring the water column pressure Pmx into sections of n≧2 to represent each section into a straight line; selecting the number n according to an allowance error of the straightened section; selecting the water column pressure γho1, γho2, γho3 corresponding to the pressure P1, P2, P3, between two points of each section divided and filling water of a reservoir or a river into the reference water column pressure measuring tubes.

3. A bubble water depth measuring system including a compressed gas generator, a pressure transducer, a reference water column pressure measuring tube, a water depth measuring tube, an arithmetic logical transducer for calculating a water depth and a drive controller for supplying/interrupting compressed gas comprising: the water depth measuring tube having a length corresponding to a water depth; the reference water column pressure measuring tubes of a U shape of n≧2 having the same inner diameter as that of the water depth measuring tube coupled therewith in a bundle; a diaphragm connected through an electromagnetic valve to one tube of the reference water column pressure measuring tube to adjust the speed of the compressed gas from a buffer tank to be supplied to the reference water column pressure measuring tubes; a transparent container connected to the other tube of the reference water column pressure measuring tubes to measure an amount of water therein and check it, periodically, and including a nipper for supplementing water reduced due to the evaporation; a tube directly mounted above the transparent container and haing the same inner diameter as that of the reference water column pressure measuring tubes and a length portion, on the surface of which scales are formed; a thin film tube made of rubber and directed to the upper portion of the tube to be expanded by a gas pressure; and a manual valve connected to the upper end of the thin film tube.

4. The bubble water depth measuring system as claimed in claim 3, in which: the reference water column pressure measuring tube is made of Teflon and Polyurethane.