CN109444995B - Method for obtaining accurate rainfall data based on bucket parameters and rainfall intensity - Google Patents
Method for obtaining accurate rainfall data based on bucket parameters and rainfall intensity Download PDFInfo
- Publication number
- CN109444995B CN109444995B CN201811606261.2A CN201811606261A CN109444995B CN 109444995 B CN109444995 B CN 109444995B CN 201811606261 A CN201811606261 A CN 201811606261A CN 109444995 B CN109444995 B CN 109444995B
- Authority
- CN
- China
- Prior art keywords
- rainfall
- bucket
- value
- tipping bucket
- sensor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/18—Testing or calibrating meteorological apparatus
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Environmental & Geological Engineering (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a method for obtaining accurate rainfall data based on bucket parameters and rainfall intensity, which is characterized by comprising the following steps of: step S1: setting a plurality of different average rainfall intensity values within a range under the total rainfall which can be determined in a test environment, testing the tipping bucket type rainfall sensor, and constructing a bucket parameter-rainfall intensity function according to a test result; and obtaining an instant rainfall function according to the bucket parameter-rain intensity function. The invention subverts the existing method for obtaining rainfall data by the tipping bucket type rainfall sensor, breaks the conventional effect that the error is smaller under the condition of higher precision, introduces the acquisition of the time interval of two times of tipping of the tipping bucket under the condition of not changing the main body component of the conventional tipping bucket type rainfall sensor, and constructs a brand-new solution scheme for obtaining the accurate rainfall by combining the characteristic attribute of the tipping bucket type rainfall sensor.
Description
Technical Field
The invention relates to the field of meteorological data acquisition, in particular to a method for obtaining accurate rainfall data based on bucket parameters and rainfall intensity.
Background
Precipitation data is used as an important meteorological data index, and has wide application and necessary significance. At present, precipitation data generally mainly take the total rainfall of statistics in a period, mainly realize the collection of real-time precipitation data and upload to meteorological monitoring terminal and summarize and count through devices such as rainfall sensor (rain gauge). The current mainstream rainfall sensor is a tipping bucket rainfall sensor.
The tipping bucket type rainfall sensor is a hydrological and meteorological instrument, and is used for measuring the rainfall in nature, and simultaneously converting the rainfall into digital information quantity expressed in a switching value form for output so as to meet the requirements of information transmission, processing, recording, display and the like. The traditional tipping bucket rainfall sensor can only achieve errors of +/-4% generally, and requirements on historical records and data analysis of weather are increasingly not met.
For the existing tipping bucket type rainfall sensor, the final total rainfall is basically obtained through preset precision (the quantity of water which can cause tipping bucket overturning is positively correlated) and the switching times (switching quantity) of a reed switch triggered by the tipping bucket overturning. If the calculation method has errors, the error is adjusted by an error correction method, so that the finally output switching value can be close to the ideal switching value as much as possible, and the errors are reduced.
However, with the current requirement for higher and higher precision of rainfall data, almost all of the adopted error correction schemes are not necessarily able to properly eliminate the influence of errors, and even new errors may be introduced, for example, although the rainfall sensor of the double-tipping bucket can eliminate dynamic loss (because rainfall is a continuous process, when the tipping bucket on one side receives rainwater and reaches a tipping bucket moment threshold, and when the tipping bucket on the other side receives rainwater, a part of rainwater is not received by the tipping bucket and is not included in the metering range of the tipping bucket type rainfall sensor, which may cause error loss), new random errors are introduced, or because of the limit of processing precision, the elimination of errors cannot be realized within the range of cost and technical capability, and the operation process of correcting errors is tedious and time-consuming and labor-consuming.
Disclosure of Invention
In order to solve the problems of defects and deficiencies in the prior art, the invention adopts a brand-new problem analysis and scheme design idea, and specifically adopts the following technical scheme without changing the main body component of the tipping bucket type rainfall sensor:
a method for obtaining accurate rainfall data based on bucket parameters and rainfall intensity is characterized in that: the rain sensor based on the tipping bucket comprises the following steps:
step S1: total rainfall H determinable in a test environment0Next, a plurality of different values are set to be (0, u)max]Mean rain intensity value of interval
Testing the tipping bucket type rainfall sensor, and constructing a bucket parameter-rain intensity function according to a test result:
F(u,y,t)＝0；
wherein u is the equivalent rain intensity of the bucket duration t, and the unit is mm/min; the bucket number y is the total overturning times of the tipping bucket, the bucket duration t is the interval time of two overturning of the tipping bucket, and the unit is s;
and according to the bucket parameter-rain intensity function, a relation established when the rain intensity is constant:
and definition of instant rain intensity
Obtaining an instant rainfall function: j (h, t) ═ 0; wherein the instant rainfall h is the rainfall generated by the equivalent rainfall intensity of the bucket duration t within the bucket duration t, and the unit is mm; total rainfall H0In units of mm; the bucket duration t is in units of s.
Preferably, in step S1, the method for testing the dump bucket rainfall sensor and constructing the bucket parameter-rain intensity function according to the test result includes:
total rainfall H given in the test environment0Next, a plurality of different values are set to be (0, u)max]The interval rainfall value u is used for measuring the total tipping bucket overturning times y value of the tipping bucket type rainfall sensor in the corresponding rainfall process under each rainfall value; and establishing a mapping corresponding relation between u and y;
wherein u ismaxIs an extreme value of rain intensity.
Preferably, in step S1, the method for testing the dump bucket rainfall sensor and constructing the bucket parameter-rain intensity function according to the test result includes:
total rainfall H determinable in a test environment0And then, providing a plurality of rainfall tests with different rainfalls, testing the tipping bucket type rainfall sensor, and obtaining a parameter-rainfall function of the bucket to be fixed through a fitting mode according to a test result:
the parameter-rain intensity function to be fixed for fighting comprises n (n is more than or equal to 2) fitting parameters, and the values of the fitting parameters are obtained by the following method:
in trialTotal rainfall H given by test environment0Next, n different ones are set to be (0, u)max]Mean rain intensity value of interval
At each mean rain intensity value
Then, the total tipping bucket overturning times y value of the tipping bucket type rainfall sensor in the corresponding rainfall process is measured, and n groups of different numerical values (u) are obtained1,y1)、(u2,y2)、……(un,yn) And the value of the n fitting parameters is determined by substituting the bucket parameter-rain intensity function.
Preferably, the specific steps of step S1 are:
step S11: based on the total rainfall H determinable in the test environment0And (3) performing multiple rainfall tests with different intensities, fitting multiple groups of test value data obtained by testing the tipping bucket rainfall sensor or constructing a parameter-rain intensity function to be fixed for the bucket according to a mathematical model established by the tipping bucket rainfall sensor:
the parameter-rain intensity function to be fixed includes n (n is more than or equal to 2) fitting parameters;
step S12: total rainfall H determinable in a test environment0Next, n different ones are set to be (0, u)max]Mean rain intensity value of interval
Testing the tipping bucket type rainfall sensor, and measuring corresponding y values under each average rainfall intensity value to obtain n groups of test values;
step S13: substituting n groups of test values into a parameter-rain intensity function to be determined, calculating to obtain n fitting parameters, and determining an instant rainfall function according to the n fitting parameters: j (h, t) ═ 0;
step S14: based on the determined instant rainfall function, testing the tipping bucket type rainfall sensor by adopting the same test conditions as those in the step S12, obtaining switching value z output by the tipping bucket type rainfall sensor according to the measured instant rainfall h and the precision epsilon of the tipping bucket type rainfall sensor, and under each average rainfall intensity value, measuring the corresponding value of the switching value z to obtain n switching values z;
step S15: comparing the n switching values z with the ideal switching value D to obtain n correction values, wherein:
epsilon is the precision of the tipping bucket rainfall sensor, if a certain z value falls within the error range determined by the ideal switching value D, the corresponding correction value is 0; when all the correction values are 0, finishing construction; as long as one correction value is not zero, the execution continues to step S16;
step S16: correcting the y value in the n groups of test values by adopting the n corrected values respectively to obtain n groups of corrected test values;
step S17: substituting the corrected n groups of test values into a parameter-rain intensity function to be determined, calculating to obtain n fitting parameters, and determining the parameter-rain intensity function again according to the n fitting parameters: f (u, y, t) is 0 and the instant rainfall function: j (h, t) ═ 0; and returns to step S14.
Preferably, the method further comprises the following steps:
step S2: in the continuous period of rainfall, obtaining instant rainfall h according to the instant rainfall function and the value of the bucket duration t;
step S3: and obtaining the total rainfall H according to the values of all the instant rainfall H in the total rainfall statistical period.
Preferably, the dump bucket rainfall sensor comprises: the water bearing device comprises a water bearing device, a funnel arranged below the water bearing device, a tipping bucket arranged below the funnel, magnetic steel driven by the tipping bucket, and a counting circuit module comprising a reed switch; the bucket time length t is the time interval between two adjacent switch pulses generated by the reed switch.
Preferably, the counting circuit module is provided with a clock circuit; the bucket time t is obtained through the time interval of two adjacent times of trigger of the reed switch, and an accurate value is calculated through a clock signal provided by a clock circuit arranged in the counting circuit module.
Preferably, in step S3, the total rainfall H is obtained by calculating the number of switching values received by the monitoring terminal from the dump bucket rainfall sensor, and the switching values are obtained by the instantaneous rainfall H and the precision epsilon of the dump bucket rainfall sensor; the switching value is obtained by dividing the accumulated value of the instant rainfall h by the precision epsilon of the tipping bucket rainfall sensor or by dividing the accumulated value of the instant rainfall h by the precision epsilon of the tipping bucket rainfall sensor.
Preferably, in step S2: and electronic signals generated by opening and closing the reed switch each time are transmitted to the monitoring terminal in real time, and the value of the bucket time t is obtained by calculation at the monitoring terminal.
Preferably, after step S1, according to the instant rainfall function and
obtaining an equivalent rain intensity function:
G(u,t)＝0；
wherein the unit of the equivalent rain intensity u of the bucket time t is mm/min.
The invention and the optimized scheme thereof subvert the existing method for obtaining rainfall data by the tipping bucket type rainfall sensor, break the conventional effect that the error is smaller under the condition of higher precision, introduce the acquisition of the time interval of two times of tipping bucket overturning under the condition of not changing the main body component of the conventional tipping bucket type rainfall sensor, and construct a brand-new solution for obtaining the accurate rainfall by combining the characteristic attributes of the tipping bucket type rainfall sensor. The method completely surpasses the methodology of analyzing and eliminating the error of the tipping bucket type rainfall sensor in the prior art.
In the construction process of the optimal scheme, the method has low requirement on test conditions, can measure that the calculation model can be accurately constructed under the condition of not controlling accurate and constant rain intensity, and has good universality.
Meanwhile, the scheme provided by the invention also provides two feasible schemes without changing the conventional equipment or adjusting the conventional monitoring terminal, does not increase excessive cost, has strong compatibility and obvious effect, is suitable for large-scale popularization and use, and has very high market value and social effect.
Drawings
The invention is described in further detail below with reference to the following figures and detailed description:
FIG. 1 is a schematic diagram of basic constitutional units and modules required for realizing embodiment 1 of the present invention;
FIG. 2 is a schematic cross-sectional view of a detailed structure of embodiment 1 of the present invention;
FIG. 3 is a perspective exploded view of a concrete structure of example 1 of the present invention;
FIG. 4 is a perspective view of a counter swing mechanism according to embodiment 1 of the present invention;
FIG. 5 is a schematic circuit diagram of a counting circuit module 1 according to an embodiment of the present invention;
FIG. 6 is a schematic circuit diagram of a counting circuit module of embodiment 1 of the present invention 2;
FIG. 7 is a schematic circuit diagram of a counting circuit module according to embodiment 1 of the present invention;
FIG. 8 is a schematic diagram of basic constitutional units and modules required for realizing embodiment 2 of the present invention;
FIG. 9 is a schematic diagram of a first rain calibrator suitable for use in embodiments of the present invention;
FIG. 10 is a schematic diagram of a second rain calibrator suitable for use in embodiments of the present invention;
FIG. 11 is a graph of h-t correlation obtained by an embodiment of the present invention;
FIG. 12 is a graph of the y-u correlation obtained by an embodiment of the present invention;
in the figure: 1-rain bearing device; 2-a scaffold; 3-a funnel; 4-tipping bucket; 5-counting the swing mechanism; 6-a counting circuit module; 7-adjusting the screw; 8-magnetic steel; 21-skip position limiter; 22-counting oscillating mechanism limit pieces; 41-tipping bucket rotating shaft; 42-left bucket chamber; 43-right bucket chamber; 51-counting swing mechanism rotating shaft; 52-a projection; 53-relief portion; 54-pointer part; 61-reed pipe; 62-an output terminal; 63-counting circuit module support; 64-a battery; 65-a clock circuit; 66-a main control chip; 67-programming the debug interface; 68-indicator light circuit; 69-a photocoupler; 610-antistatic protection circuit; 611-USB interface; 612-a power supply circuit; 613-power decoupling filter circuit; 100-a housing; 200-base.
Detailed Description
In order to make the features and advantages of the present invention comprehensible, 2 embodiments accompanied with figures are described in detail as follows:
the basic gist of the present invention is: the time interval (bucket duration t) of two-time overturning of the tipping bucket is collected, and the accurate total rainfall H is directly obtained by combining the correlation characteristics of the tipping bucket type rainfall sensor and the bucket duration t. The correlation characteristic of the tipping bucket rainfall sensor and the bucket duration t is one of characteristic attributes of the tipping bucket rainfall sensor. Characteristic attributes of the dump bucket rain sensor exist based on a specific dump bucket rain sensor, and can be used for describing the characteristics of the dump bucket rain sensor which are different from other dump bucket rain sensors, such as: the characteristic that the left and right tipping buckets can turn over when receiving a certain amount of rainwater, the characteristic that how much time is needed for reaching a stable state from turning beginning to turning ending, the characteristic of the total bucket number under a specific total rainfall and a specific constant rainfall and the like, and the correlation characteristic of the tipping bucket type rainfall sensor and the bucket duration t is a general name of the characteristic attribute which can be obtained and determined by the tipping bucket type rainfall sensor according to the bucket duration t. That is, based on the scheme of the present invention, after a specific dump-bucket rainfall sensor (including the insight of its characteristic attributes) is determined, in a rainfall scene, only the duration t of the dump bucket needs to be measured, and the accurate total rainfall H can be obtained.
In the invention, the correlation characteristic of the tipping bucket rainfall sensor and the bucket duration t is characterized by an instant rainfall function J (H, t) being 0, and the total rainfall H can be obtained through the instant rainfall H.
In the present invention, the definition of the instantaneous rainfall h is: the equivalent rain intensity of the fighting time length t generates the rainfall within the fighting time length t. Wherein for "equivalent rain intensity":
we consider fromHowever, in the environment of rainfall, when the time point is x, the rainfall intensity is R (x), and then
Thus, the exact amount of rainfall over the bucket duration t can be expressed as:
based on the fundamental design principle of the tipping bucket type rainfall sensor, the rainfall change condition within the bucket time t cannot be obtained by the measurement of the tipping bucket type rainfall sensor, so that the integral type cannot be used for obtaining an accurate solution. However, due to the objective presence of accurate rainfall over the bucket duration t, according to the median theorem of integrals: there must be a certain rain intensity value
Satisfy the requirement of
Although the rain intensity value
Nor is it available, but we have thus obtained the original meaning of "equivalent rain intensity".
In the invention, the equivalent rain intensity is defined as the rain intensity value which is uniquely corresponding to each bucket time length value obtained by the determined tipping bucket type rain sensor under the condition of constant rain intensity, so that the fundamental basis of calculating the rainfall in the bucket time length t is obtained. This means that after the form of the instantaneous rainfall function is determined, the rainfall in each bucket time period t can be approximately obtained directly through the bucket time period t, so as to achieve the purpose of directly obtaining the total rainfall H.
It should be noted that the implementation of the solution of the present invention does not depend on a dump-bucket rain sensor of a specific structure, and the requirements of the dump-bucket rain sensor device as a material basis for data measurement are only that: the method can obtain a t value (the t value of a common tipping bucket type rainfall sensor can be obtained by measuring the time interval between two adjacent switch pulses generated by a reed switch), and has certain relativity between the instant rainfall h and the bucket duration t, namely h-t related characteristic attributes (for example, certain tipping bucket type rainfall sensors provided with an upper tipping bucket and a lower tipping bucket for eliminating errors generated by dynamic loss have no determined relativity with the instant rainfall h or equivalent rainfall intensity due to the structural characteristics of the device, so the method is not in the applicable range of the scheme of the invention). Therefore, in a general case, it can be understood that any dump bucket rainfall sensor (which can be generally considered as a dump bucket rainfall sensor with t decreasing trend along with the increase of real-time rainfall intensity and only one dump bucket) with different bucket time lengths t capable of reflecting the change of the instant rainfall h can be used as the basis of the method of the present invention. The specific structure or circuit configuration of the following embodiments provided by the present invention is not intended to limit the application scope of the method of the present invention, but is merely a case for facilitating those skilled in the art to better understand a specific implementation of the solution of the present invention.
As shown in fig. 1, in the first embodiment of the present invention, the basic arithmetic unit is located locally in the dump bucket rainfall sensor for the purpose of compatibility between the existing dump bucket rainfall sensor and the docking of the monitoring terminal generally installed in the weather station and the manner of data interaction. The embodiment improves the counting circuit module circuit which only has the function of generating counting pulses for opening and closing the reed pipe 61 in the conventional tipping bucket type rainfall sensor, and a clock circuit 65 serving as a local acquisition t value and a main control chip 66 serving as a data operation and storage core are added in the counting circuit module circuit. The main control chip 66 generates a switching value capable of accurately calculating the rainfall amount after performing operation according to the value t, and transmits the switching value to the monitoring terminal through the output terminal 62.
The following describes in detail the data processing and calculating method adopted by the main control chip 66 in this embodiment and the process of generating the same.
(1) Obtaining y-u correlation features
Based on the definition of the instantaneous rain quantity h, although the most straightforward design solution should strive to obtain an equivalent rain intensity for the bucket duration t, on the one hand a very high constancy of the rain intensity is difficult to obtain in a test environment, and on the other hand an exact value of the constant rain intensity is difficult to obtain. Therefore, in the aspect of constructing the characteristic parameters and the correlation between the parameters of the tipping bucket rainfall sensor, the bucket number y is a quantity which can be accurately measured through a general rainfall test, at this time, although an accurate constant rainfall intensity value cannot be given, so that a direct u-t correlation is difficult to obtain, the direct u-t correlation can be tried to be replaced by an average rainfall intensity value, and the correlation between the average rainfall intensity value and the bucket number y is obvious, so that the characteristic attribute of the tipping bucket rainfall sensor is obtained from the y-u correlation in the embodiment.
In this example, the y-u correlation was constructed by the following method:
total rainfall H determinable in a test environment0Next, a plurality of different values are set to be (0, u)max]Mean rain intensity value of interval
Testing the tipping bucket type rainfall sensor, and constructing a bucket parameter-rain intensity function according to a test result:
F(u,y,t)＝0；
wherein u is the equivalent rain intensity of the bucket duration t, and the unit is mm/min; the bucket number y is the total overturning times of the tipping bucket.
The goal of the above construction of the bucket parameter-rain intensity function is to determine the total rainfall H0Under the precondition of (1), the total rainfall is controlled to be kept unchanged, the test is carried out by changing the magnitude of the rain intensity, and the bucket number y and the average rain intensity value given by the test at the moment
There will be significant correlation, and the present embodiment will establish a mathematical expression of this correlation by the scoop parameter-rain intensity function. The bucket parameter-rain intensity function means, as a characteristic attribute of a bucket rain sensor, that it characterizes: when the total rainfall is a determined value H0Time, bucket number y, bucket duration t, etcThe effective rain intensity u has certain relevance. When expressing the correlation of y-u using the same, the average rain intensity value is obtained as long as the change rate of the rain intensity is not large (approximately constant)
The equivalent rain intensity u of the fighting duration t can be approximately regarded as follows:
and eliminating the bucket time length t, so as to construct and obtain the y-u correlation characteristic.
Wherein, in the present embodiment, the total rainfall amount H that can be determined is obtained0The total rainfall can be determined either by the volume of the standard ball in the rain calibrator as shown in fig. 9 or by a container with a weight scale for receiving the dump body effluent below the dump body rain sensor as shown in fig. 10.
The adjustment of the rainfall is generally performed by adjusting a valve (e.g., a water outlet solenoid valve). The specific value of the rainfall intensity is not needed to be known when the scheme of the embodiment is implemented, and generally, the provided rainfall intensity is not required to be completely constant in the process of simulating rainfall (but should be as constant as possible in order to control an error, in this point, it is enough to achieve the purpose that the state of the water outlet electromagnetic valve in the rainfall calibrator shown in fig. 9 is controlled to be constant in the process of simulating rainfall), so that the test condition meeting the scheme of the embodiment can be completed only by using a general rainfall calibrator and matched equipment thereof.
The design scheme of the rainfall calibrator shown in fig. 10 can be used as a test device for providing constant rainfall intensity in the embodiment, and the rainfall calibrator comprises a water storage container, when in test, a water pump always injects water into the water storage container, the flow of the injected water is greater than the flow of the discharged water, the redundant water flows out from an overflow port, so that the height difference between the water outlet and the highest water level of the water storage container is kept unchanged, the flow of the constant water outlet through hole area is stable, namely the rainfall intensity is constant, and the output rainfall intensity can be accurately controlled and adjusted through a stepping motor and a flow control valve.
More specifically, the present embodiment provides a path constructed by three bucket parameters, the rain intensity function.
One is a construction mode of a discrete function, and specifically includes:
total rainfall H determinable in a test environment0Next, a plurality of rainfall tests of different intensity are provided, at each average rainfall value
Then, the value of the number y of the buckets is measured, and various rainfall scenes such as light rain, medium rain, heavy rain and the like are covered when the rainfall intensity is controlled and adjusted, and the rainfall intensity is ensured to be constant as far as possible.
Wherein at each average rain intensity value
And then, a plurality of tests can be carried out, and the obtained plurality of y values are determined by adopting statistical methods such as averaging, median and the like.
And by (near) constant rain intensity
And establishing a mapping corresponding relation between u and y according to the approximate relation of the equivalent rain intensity u of the bucket duration t.
This construction has the advantage of being directly usable, when the rain strength is sufficiently constant (in this case)
Basically, the equivalent rain intensity u of the bucket duration t) and the test point density and distribution are reasonable enough, the accuracy of the correlation of y-u can be ensured, and the higher the constancy of the rain intensity provided by the test,
the smaller the error between the equivalent rain intensity u and the fighting duration t is, the more accurate the obtained fighting parameter-rain intensity function is.
The second is the construction mode of the continuous (fitting) function:
by passing
The distribution and the form of the discrete points, and the bucket number-bucket duration function which meets the following conditions is constructed by taking fitting as a target:
in that
The internal requirements are as follows: y is less than or equal to D, and when uα＜uβAt a time yα≥yβWherein u ismaxIs an extreme value of rain intensity and an ideal switching value
Is a function determined by the bucket parameter-rain intensity function.
The bucket parameter-rain intensity function can be pre-constructed by presetting a function analytic expression comprising n (n is more than or equal to 2) fitting parameters (the more the number of the fitting parameters is, the more the fitting effect is accurate, but the larger the general calculated amount is), the pre-constructed analytic expression form is generally obtained by fitting an image formed by the distribution of a large number of test values, and the distribution of the test values of the tipping bucket rain sensor adopting similar construction and principle generally has common characteristics and rules, so the pre-constructed function analytic expression is common to the same type of tipping bucket rain sensors. As shown in fig. 12, the curve is a coordinate diagram of the y-u correlation curve of the dump bucket rainfall sensor used in this embodiment, and the pre-constructed bucket parameter-rain intensity function may be based on fitting to the curve.
Thereafter, the total rainfall H determinable in the test environment0Next, n different ones are set to be (0, u)max]Mean rain intensity value of interval
Obtaining n different sets of test values (u)1,y1)、(u2,y2)、……(un,yn). The obtained n groups of test values are taken into a bucket parameter-rainIn a strongly functional equation, the values of the n fitting parameters are thus determined, the more accurately the values of the fitting parameters are determined as the test provides a constant rain intensity.
Thereby completing the determination of the analytic expression of the bucket parameter-rain intensity function.
The construction and calculation of the specific scoop parameter-rain intensity function analytic expression are provided below:
for example, from the test data obtained in a number of tests, shaped as shown in fig. 12, the scoop parameter-rain intensity function is pre-constructed to the following form (hyperbolic function):
wherein a, b and c are fitting parameters.
By means of a clamshell-type rain sensor as shown in FIG. 9 at H0The following three sets of data were obtained from experiments performed in a 10mm scenario:
light rain: y: 119;
rain: y: 113;
heavy rain: y: 108;
the above data is taken into:
obtaining by solution:
thereby determining the concrete form of the analytic expression of the bucket parameter-rain intensity function.
Compared with the construction mode of a discrete function, the mode of constructing the continuous function has the advantages that under the condition of the determination of the function analytic formula of the pre-construction, the required test times are generally less for obtaining the y-u correlation with the same precision, the requirement on the constancy of the rain intensity provided by the test is lower than that of the construction mode of the discrete function, and the requirement on the constancy of the rain intensity provided by the test is lower than that of the discrete function
Correlation data that are difficult to directly measure by experiment, such as approaching 0, can also be fitted to obtain satisfactory output values.
Furthermore, although the pre-constructed function analytic expressions are generally preferred to be in a functional form that is easy to solve to reduce the amount of computation. However, even if the constructed function analytic expression cannot solve the exact solution of the fitting parameter by the method, the function analytic expression can obtain an approximate value by a computer program, and therefore, the method also belongs to a feasible scheme.
For example, the scoop parameter-rain intensity function can be pre-constructed as the following form (polynomial function):
F(y,u)＝y-(a0+a1u+a2u2+a3u3+…+anun)＝0
in this case, it is difficult to directly find an exact solution, but an acceptable approximate solution can be obtained by a computer program.
Meanwhile, the accuracy of the constructed continuous function can be verified or corrected in a test mode, and the reasonability of the constructed function can be verified through the inherent characteristics of the correlation of y-u, such as
Thirdly, a function analytical formula is constructed for a mathematical model established by the tipping bucket type rainfall sensor:
although the above two construction methods do not depend on mathematical modeling and analysis of the tipping bucket type rainfall sensor, the construction method based on the y-u correlation provided by the embodiment does not exclude that a model of an available function analytic expression is derived through a mathematical model established by the tipping bucket type rainfall sensor, and the composition of the bucket parameter-rain intensity function is finally determined through a test method.
For example, when only the dynamic loss is considered, the constant rain intensity value is set to
The unit is mm/min, the ideal fighting number without dynamic loss is C, and according to the definition of rainfall H, the following steps are provided:
the amount of rainfall for which dynamic losses can be obtained is
Then, the duration of each bucket in the dynamic loss state within the duration t is set as tau, and the rainfall of the dynamic loss can be obtained as
The following y-u dynamic loss equation can be obtained by combining two equations, which is equivalent to a pre-constructed bucket parameter-rain intensity function:
(C-y)ut＝yuτ
where C and τ may correspond to fitting parameters in a fitting process;
and then according to the relation formula satisfied when the rain intensity is constant:
an equation characterizing the y-u correlation can be obtained:
by means of a clamshell-type rain sensor as shown in FIG. 9 at H0The following two sets of data were obtained from experiments performed in a 10mm scenario:
the first measurement: y is equal to 108, and y is equal to 108,
and (3) second measurement: y is equal to 113. the number of the first symbols,
the substitution can be solved to obtain C-123.7936913, tau-0.20310817
Thereby determining the concrete form of the analytic expression of the bucket parameter-rain intensity function.
(2) Selection of output mode
And (6) acquiring the total rainfall.
According to the determined bucket parameter-rain intensity function, a relation formula is established when the rain intensity is constant:
and definition of instant rain intensity
Obtaining an instant rainfall function: j (h, t) ═ 0; wherein the unit of the instant rainfall h is mm; total rainfall H0In units of mm; the bucket duration t is in units of s.
According to the instant rainfall function, under the actual use state of the tipping bucket type rainfall sensor, the instant rainfall h can be obtained directly according to the value of the bucket duration t. The constructed h-t correlation is generally shown in the graph of FIG. 11.
As a preferable scheme provided by this embodiment, there are:
wherein a, b and c are fitting parameters.
As another preferable scheme provided by this embodiment, there are:
wherein C is the ideal bucket number without dynamic loss, and tau is the time length of each bucket in the dynamic loss state within the time length t.
Since the rainfall intensity of actual rainfall in the natural world is basically impossible to be a constant value even within a bucket duration, the total rainfall H obtained by accumulating all instant rainfall H within the whole rainfall statistic period is substantially an approximate operation similar to integral, however, according to the characteristic, the bucket duration t under the same rainfall intensity is shorter as the precision of the tipping bucket type rainfall sensor is higher, and according to the common knowledge, the rainfall intensity change of natural rainfall is naturally less obvious within a shorter time, and the approximate calculation error of the total rainfall H obtained by accumulating all instant rainfall H within the whole rainfall statistic period is smaller. The characteristic is completely opposite to the existing scheme of calculating rainfall through a tipping bucket type rainfall sensor, is obviously better, and is equivalent to directly obtaining a technical path for realizing ultrahigh-precision rainfall data acquisition.
Specifically, as for how to directly obtain the total rainfall H in the total rainfall statistic period through the instant rainfall H, it is the most ideal way that the skip bucket type rainfall sensor directly outputs each calculated instant rainfall H or the total rainfall H in one rainfall period, and in order to be compatible with the existing rainfall acquisition system, this embodiment provides two technical implementation paths of redefined switching quantities calculated and obtained from the counting circuit module by dividing the instant rainfall H and the precision epsilon of the skip bucket type rainfall sensor:
one of them is that: the method is characterized in that the method is obtained by dividing an accumulated value of the instant rainfall h by the precision epsilon of the tipping bucket type rainfall sensor, namely in a total rainfall statistic period, a counting circuit module performs division operation on the accumulated value of the instant rainfall h and the precision epsilon of the tipping bucket type rainfall sensor in a preset period, the obtained divisor is summed with the mantissa of the previous period, then the integer part of the sum is converted into metering signals with the same number to be used as the switching value z to perform periodic batch output, and the decimal part is used as the mantissa of the period.
The second is that: the real-time rainfall h is obtained by accumulating after being divided by the precision epsilon of the tipping bucket rainfall sensor, namely in the total rainfall statistic period, the counting circuit module divides each real-time rainfall h by the precision epsilon of the tipping bucket rainfall sensor, takes an integer part as the switching value z output, accumulates a decimal part, converts the integer part of an accumulated value into the switching value z output with the same number, and simultaneously keeps the decimal part of the accumulated value to participate in the next accumulation.
The real-time rainfall H obtained through local calculation can be accurately obtained through the two paths in a manner of being compatible with the rainfall obtaining and calculating of the existing monitoring terminal, so that the monitoring terminal can obtain the accurate total rainfall H. Although the present embodiment provides the data signal sent locally to the monitoring terminal in the form of the switching value, the switching value generated by the prior art solution is completely based on the counting signal of the reed switch 61, and the switching value generation mechanism of the present embodiment is completely different as described above.
② obtaining the instant raininess
According to the constructed instant rainfall function and the scheme of the invention, the instant rainfall h is defined as follows:
according to the constructed instant rainfall function and the scheme of the invention, the instant rainfall h is defined as follows:
G(u,t)＝0；
wherein the unit of the equivalent rain intensity u of the bucket time t is mm/min.
Under the actual use state of the tipping bucket type rainfall sensor, the equivalent rainfall intensity value of the bucket duration t can be directly obtained according to the value of the bucket duration t.
If the constructed discrete function is a discrete function, a u value corresponding to a t value closest to an actually measured t value can be directly taken as an output value on a corresponding relation table of u and t; for the continuous function, the actually measured t value is directly substituted into the analytical expression to obtain the corresponding u value as an output value.
If all the obtained equivalent rain intensity values u are represented on the time axis, a quite accurate and real-time picture capable of reflecting the instant rain intensity situation is obtained, and the change situation of the actual rain intensity along with the time can be quantitatively known, which can not be realized by the prior art.
Furthermore, in the present embodiment, it can be found that although the theoretical basis of the correlation characteristic construction of the present embodiment is based on the constant rain intensity, since a test scheme of the non-constant rain intensity can be adopted, the test values are generally measured under the non-constant rain intensity, which inevitably results in inaccuracy of the finally determined instant rain function (in an actual scenario, even though the test is performed by the rain calibrator as shown in fig. 10, it is difficult to obtain the truly completely constant rain intensity). Therefore, in order to further improve the accuracy of the instantaneous rainfall function obtained in the above embodiment, in consideration of the correlation between the fitting parameters and the test values in the construction scheme of the fitting method, a preferred calibration scheme dedicated to the scheme of constructing the instantaneous rainfall function by the continuous function is provided:
the calibration scheme requires the use of a calibrator which provides a certain total rainfall H, having at least the rainfall calibrator shown in FIG. 90And the test device can reproduce the function of a specific rainfall scene, and the general commonly used rainfall calibrator meets the condition.
In this embodiment, no matter the fitting method or the method for establishing the mathematical model is used to pre-construct the obtained determination function model, the undetermined parameters (such as the fitting parameters or the parameters such as C and τ, etc. that need to be determined through experiments in the mathematical model) that need to be determined through experiments are involved, so that, under the condition that the coincidence between the test value and the correlation model has a deviation, the calibration method considers that a method for calibrating the bucket number y in the test value is used, and specifically includes the following steps:
step S11: based on the total rainfall H determinable in the test environment0And (3) performing multiple rainfall tests with different intensities, fitting multiple groups of test value data obtained by testing the tipping bucket rainfall sensor or constructing a parameter-rain intensity function to be fixed for the bucket according to a mathematical model established by the tipping bucket rainfall sensor:
the parameter-rain intensity function to be fixed includes n (n is more than or equal to 2) fitting parameters;
step S12: total rainfall H determinable in a test environment0Next, n different ones are set to be (0, u)max]Mean rain intensity value of interval
Testing the tipping bucket type rainfall sensor, and measuring corresponding y values under each average rainfall intensity value to obtain n groups of test values;
step S13: substituting n groups of test values into a parameter-rain intensity function to be determined, calculating to obtain n fitting parameters, and determining an instant rainfall function according to the n fitting parameters: j (h, t) ═ 0;
step S14: based on the determined instant rainfall function, testing the tipping bucket type rainfall sensor by adopting the same test conditions as those in the step S12, obtaining switching value z output by the tipping bucket type rainfall sensor according to the measured instant rainfall h and the precision epsilon of the tipping bucket type rainfall sensor, and under each average rainfall intensity value, measuring the corresponding value of the switching value z to obtain n switching values z;
step S15: comparing the n switching values z with the ideal switching value D to obtain n correction values, wherein:
epsilon is the precision of the tipping bucket rainfall sensor, if a certain z value falls within the error range determined by the ideal switching value D, the corresponding correction value is 0; when all the correction values are 0, finishing construction; only one repair is neededIf the positive value is not zero, the step S16 is continued;
step S16: correcting the y value in the n groups of test values by adopting the n corrected values respectively to obtain n groups of corrected test values;
step S17: substituting the corrected n groups of test values into a parameter-rain intensity function to be determined, calculating to obtain n fitting parameters, and determining the parameter-rain intensity function again according to the n fitting parameters: f (u, y, t) is 0 and the instant rainfall function: j (h, t) ═ 0; and returns to step S14.
Through the steps, the error of the output switching value is used for representing the non-conformity of the test value and the pre-constructed function model, the switching value error is used for generating a correction value, the correction value is used for correcting and adjusting the bucket number in the test value in turn, and the ideal state is gradually approached through cyclic correction.
In the method, the correction value may be generated by calculating a difference or a scaling factor, for example, in the case that the ideal switching value D is 100 and the allowable error thereof is 1%, that is, the determined error range is 99-101, if the switching value z outputted in a certain test is 97, the correction value may be, by a difference method: +3(97+3 equals 100), and the bucket number in the corresponding test value can be adjusted by + 3; by scaling, the correction value may be: 1.03(100/97 equals 1.03), the number of buckets in the corresponding test value can be adjusted by multiplying the number by 1.03.
In the subsequent operation, if the switching value z output in a certain test falls in the range of 99-101, namely | z-D | ≦ 1, the correction value is directly determined to be 0, when the correction values corresponding to all test groups are 0, the correction process is ended, and the parameter value corresponding to the last calculation and the instant rainfall function are finally determined.
The optimal calibration method can be used for calibrating the characteristic function model before the tipping bucket rainfall sensor is used, so that the accuracy of the instant rainfall function determined through a non-constant rainfall intensity test can be ensured to be sufficiently consistent with the real situation, and the method is also suitable for error calibration and adjustment of the tipping bucket rainfall sensor in the long-term use process.
To further illustrate a possible implementation mechanism of the method of this embodiment, a specific possible implementation mechanism provided in this embodiment is specifically described below with reference to fig. 1 to 7:
in terms of the structure of the device, as shown in fig. 2 and 3, the entire device of the present embodiment, as a typical dump bucket type rain sensor having a single dump bucket, includes: the water receiver comprises a water receiver 1 arranged above a support 2, a funnel 32 arranged at the top of the support 2, a skip bucket 4 arranged below the funnel 3, a counting swing mechanism 5 arranged below the skip bucket 4, and a counting circuit module 6 comprising a reed switch 61.
The tipping bucket 4 is hinged with the bracket 2 through a tipping bucket rotating shaft 41; the counting swing mechanism 5 is hinged with the bracket 2 through a counting swing mechanism rotating shaft 51; the bucket rotation shaft 41 and the counter swing mechanism rotation shaft 51 are provided on the symmetrical surfaces of the left bucket chamber 42 and the right bucket chamber 43 of the bucket 4.
As shown in fig. 2, the device of the present embodiment is integrally mounted in a housing formed by a case 100 and a base 200, and a counter circuit module 6 is disposed on the side of a holder 2, provided with an output terminal 62 and a battery 61, and integrally assembled on a counter circuit module holder 63.
As shown in fig. 4, the counting and swinging mechanism 5 of the present embodiment includes a protrusion 52 protruding along the extending direction of the left bucket chamber 42 and the right bucket chamber 43, a receding part 53 surrounding the skip rotating shaft 41, and a pointer part 54 fixed with a magnetic steel 8; the highest point of the pointer 54 is as high as the reed switch 61 arranged on the side wall of the support 2 (a scheme slightly lower than the reed switch 61 is also feasible, and only the normally open contact of the reed switch 61 is required to be closed when the pointer 54 is at the highest point of the swing, and the normally open contact of the reed switch 61 is required to be opened when the pointer 54 deviates from the highest point); the magnetic steel 8 is arranged on the symmetrical surface of the counting swing mechanism 5.
The swinging ranges of the two sides of the dump bucket 4 and the counting swinging mechanism 5 are respectively limited by two limiting parts, including two dump bucket limiting parts 21 shaped like cylindrical protrusions and two counting swinging mechanism limiting parts 2222 shaped like cylindrical protrusions; the limiting part is used for controlling the maximum swinging range of the tipping bucket 4 and the counting swinging mechanism 5, so that the path of each overturning is controllable.
Furthermore, the form of the counter swing mechanism 5 needs to satisfy: when the counting and swinging mechanism 5 is at any angle within the allowable swinging range, the path along which the dump bucket 4 swings towards two sides intersects with the protruding part 52 of the counting and swinging mechanism 5, so that each overturning and swinging of the dump bucket 4 can correspondingly drive the counting and swinging mechanism 5 to swing, and the magnetic steel 8 can trigger the swinging of the reed pipe 61 to complete one counting.
Meanwhile, the gravity center of the tipping bucket 4 is higher than the tipping bucket rotating shaft 41; the gravity center of the counting swing mechanism 5 should be higher than the rotating shaft 51 of the counting swing mechanism to avoid the situation that the state change of the reed pipe 61 cannot be triggered because the magnetic steel 8 swings repeatedly near the reed pipe 61.
The above components mainly adjust the counting assembly composed of the reed pipe 61 (magnetic reed switch) and the magnetic steel 8, except for the most basic necessary components of the conventional dump-bucket rainfall sensor. The scheme of the embodiment arranges the magnetic steel 8 on the counting swing mechanism 5 structurally separated from the tipping bucket 4, thereby eliminating the error caused by the influence of the magnetic force generated when the magnetic steel 8 is close to the reed pipe 61 on the whole moment of the device. When the skip bucket 4 turns over and the skip bucket 4 chamber on the rain-bearing side drops, the corresponding one-side protruding part 52 acts on the bottom of the skip bucket 4 chamber, and the kinetic energy of the drop bucket 4 chamber drives the counting swing mechanism 5 to swing, so that the magnetic steel 8 arranged on the pointer part 54 swings and passes through the near point of the reed switch 61, and the triggering of the magnetic reed switch can be completed under the condition of not influencing the moment of the skip bucket 4.
It should be noted that, in the conventional dump bucket type rain sensor, the magnetic steel 8 only needs to be set on the counting swing mechanism 5 and the pointer part 54 thereof as described in the device scheme of this embodiment, but only needs to be set on the magnetic steel 8 to swing to a certain characteristic point, which can trigger the normally open contact of the reed pipe 61 to close, and when the magnetic steel 8 deviates from the characteristic point, the normally open contact of the reed pipe 61 is opened, and of course, the specific setting mode of the reed pipe 61 can be flexibly adjusted to match the specific setting position of the magnetic steel 8.
Preferably, in the present embodiment, the dump body 4 is further provided with an adjusting screw 7 capable of adjusting the center of gravity perpendicular to the dump body rotating shaft 41, and the adjusting screw 7 is capable of linearly adjusting the center of gravity of the combined body of the dump body 4 and the adjusting screw 7 according to the disclosure of chinese patent publication No. CN 107765349A. And the rain bearing surface of the tipping bucket 4 is a cylindrical surface. In this preferred structure, the adjusting screw 7 should be disposed inside the pointer 54 to ensure that the distance between the magnetic steel 8 and the reed switch 61 is not too large to trigger the reed switch 61.
In terms of device circuits, as shown in fig. 5 to 7, the main control chip 66 as a control, storage and operation core adopts SMT32L051C8, which is connected to the clock circuit 65, the programming and debugging interface 67, the reed pipe 61, the indicator lamp circuit 68, the photocoupler 69, the power supply circuit 612, the anti-static protection circuit 613 and the USB interface 611, respectively.
The reed switch counting circuit formed by the reed switch counting circuit is matched with magnetic steel 8 driven by the swinging of the tipping bucket 4 in the tipping bucket type rainfall sensor, when the tipping bucket 4 is overturned to a preset position, the distance between the magnetic steel 8 and the reed switch 61 is nearest, so that a normally open switch of the reed switch 61 is closed and conducted under the action of a magnetic field of the magnetic steel 8, and when the magnetic steel 8 moves continuously and is slightly far away from the reed switch 61, the influence of the magnetic steel on the magnetic field of the reed switch 61 is weakened, and the normally open switch of the reed switch 61 is disconnected. The pulse generated by the opening and closing of the reed switch 61 can be used as the counting basis of the overturning times of the tipping bucket.
The clock circuit 65, which is an important component implemented in the present embodiment, is composed of a crystal oscillator, a crystal oscillation control chip, and a capacitor, and is provided for the purpose of measuring and measuring the time interval (bucket length) between two adjacent triggerings of the reed switch 61.
The main control chip 66 undertakes data storage and calculation in the process of acquiring more accurate rainfall data, on one hand, the instant rainfall can be calculated through a preset instant rainfall function and the bucket duration, and then the output switching value is obtained through the accumulation of the instant rainfall and the division calculation of the precision of the tipping bucket rainfall sensor.
In the present embodiment, the switching amount transmitted to the precipitation amount monitoring terminal is output from the output terminal 62 via the photocoupler 69. The output mode can be replaced by a wireless signal transmission mode, so that the switching value output by the tipping bucket rainfall sensor does not need to be transmitted through a specific line.
In order to better achieve the effect of accurate calculation of the main control chip 66, the embodiment provides the USB interface 611 capable of directly performing data interaction with the rainfall calibrator or the upper computer, so as to achieve acquisition of output data of the rainfall sensor or parameter update, and the programming and debugging interface 67 capable of downloading and obtaining the accurate calibration program, where the program language converted by the method of the embodiment and the corresponding parameter value thereof can be downloaded to the local or updated through the port. An anti-static protection circuit 610 is disposed between the USB interface 611 and the main control chip 66 to ensure the overall security of the circuit system.
The power circuit 612 of this embodiment is a 1.8V power supply voltage stabilizing circuit, and is powered by a 8000 ma battery 64, and is also provided with a power decoupling filter circuit 613 to enhance the stability of the power supply. The large-capacity battery can be matched with low-power consumption chips and devices, so that the device can work for a long enough time (years) without charging.
An indicator light circuit 68 is also provided for indicating the state of charge of the power supply.
In contrast, in the second embodiment of the present invention, since the design concept is based on the consideration of the existing dump-bucket rain sensor device, it is not required that the dump-bucket rain sensor device has a specially designed counting circuit module as shown in fig. 5-7 of the first embodiment, and the basic device conditions only need to satisfy the basic structure of the existing dump-bucket rain sensor as shown in fig. 8.
In this embodiment, the task of obtaining the value of the bucket duration t and calculating the total rainfall according to the instant rainfall function is executed at the monitoring terminal, and the electronic signal generated by the reed switch of the local tipping bucket type rainfall sensor is transmitted to the monitoring terminal in real time. At this time, the monitoring terminal can restore the bucket duration corresponding to the tipping bucket type rainfall sensor through the time interval between two adjacent electronic signals, and then the bucket duration is calculated through the method provided by the embodiment. (of course, the determination of the instant rainfall function must be performed locally, and only the determined characteristic function can be placed in the monitoring terminal to calculate the total rainfall according to the duration of the bucket)
The realization idea has the advantages that the construction of the instant rainfall function can be completed only by locally testing the tipping bucket type rainfall sensor without basically changing the hardware of the tipping bucket type rainfall sensor of the existing rainfall collection point (which is not suitable for the double-tipping bucket type rainfall sensor).
Accordingly, however, with this solution, the existing monitoring terminal needs to be modified, and the measuring and calculating precision time and the corresponding device or calculation program for performing the solution calculation need to be increased.
Finally, it should be specifically noted that although the 2 embodiments provided above all convert the acquisition scheme of the bucket duration (the time interval of two-time overturning of the dump bucket) into: and acquiring the time interval between two adjacent switching pulses generated by the reed switch. However, in the solution of the present invention, the obtaining manner of the bucket duration is not limited by the specific device combination of the reed switch 61 and the magnetic steel 8, and the bucket duration can be obtained by using a plurality of other proximity switches or contact switches similar in principle, such as an infrared sensor, a contact sensor, and the like, and the solution of using the reed switch 61 and the magnetic steel 8 in the above 2 embodiments is only because the solution is a common solution of the current tipping bucket type rain sensor, and the use effect is more superior because no additional power supply is required.
The present invention is not limited to the above preferred embodiments, and any other various methods for obtaining accurate rainfall data based on bucket parameters and rainfall intensity can be obtained according to the teaching of the present invention.
Claims (10)
1. A method for obtaining accurate rainfall data based on bucket parameters and rainfall intensity is characterized in that: the rain sensor based on the tipping bucket comprises the following steps:
step S1: total rainfall H determinable in a test environment0Next, a plurality of different values are set to be (0, u)max]Mean rain intensity value of interval
Testing the tipping bucket type rainfall sensor, and constructing a bucket parameter-rain intensity function according to a test result:
F(u,y,t)＝0；
wherein u is the equivalent rain intensity of the bucket duration t, and the unit is mm/min; the bucket number y is the total overturning times of the tipping bucket, the bucket duration t is the interval time of two overturning of the tipping bucket, and the unit is s;
and according to the bucket parameter-rain intensity function, a relation established when the rain intensity is constant:
and definition of instant rain intensity
Obtaining an instant rainfall function: j (h, t) ═ 0; wherein the instant rainfall h is the rainfall generated by the equivalent rainfall intensity of the bucket duration t within the bucket duration t, and the unit is mm; total rainfall H0In units of mm; the bucket duration t is in units of s.
2. The method for obtaining precise rainfall data based on scoop parameters and rainfall intensity of claim 1, wherein in step S1, the method for testing the scoop-type rainfall sensor and constructing the scoop parameter-rainfall intensity function according to the test result comprises:
total rainfall H given in the test environment0Next, a plurality of different values are set to be (0, u)max]The interval rainfall value u is used for measuring the total tipping bucket overturning times y value of the tipping bucket type rainfall sensor in the corresponding rainfall process under each rainfall value; and establishing a mapping corresponding relation between u and y;
wherein u ismaxIs an extreme value of rain intensity.
3. The method for obtaining accurate rainfall data based on scoop parameters and rain intensity of claim 1, wherein: in step S1, the method for testing the dump bucket rainfall sensor and constructing the bucket parameter-rain intensity function according to the test result includes:
total rainfall H determinable in a test environment0And then, providing a plurality of rainfall tests with different rainfalls, testing the tipping bucket type rainfall sensor, and obtaining a parameter-rainfall function of the bucket to be fixed through a fitting mode according to a test result:
the parameter-rain intensity function to be fixed for fighting comprises n (n is more than or equal to 2) fitting parameters, and the values of the fitting parameters are obtained by the following method:
total rainfall H given in the test environment0Next, n different ones are set to be (0, u)max]Mean rain intensity value of interval
At each mean rain intensity value
Then, the total tipping bucket overturning times y value of the tipping bucket type rainfall sensor in the corresponding rainfall process is measured, and n groups of different numerical values (u) are obtained1,y1)、(u2,y2)、……(un,yn) And brought into saidA bucket parameter-rain intensity function, from which the values of n fitting parameters are determined.
4. The method for obtaining accurate rainfall data based on scoop parameters and rain intensity of claim 1, wherein: the specific steps of step S1 are:
step S11: based on the total rainfall H determinable in the test environment0And (3) performing multiple rainfall tests with different intensities, fitting multiple groups of test value data obtained by testing the tipping bucket rainfall sensor or constructing a parameter-rain intensity function to be fixed for the bucket according to a mathematical model established by the tipping bucket rainfall sensor:
the parameter-rain intensity function to be fixed includes n (n is more than or equal to 2) fitting parameters;
step S12: total rainfall H determinable in a test environment0Next, n different ones are set to be (0, u)max]Mean rain intensity value of interval
Testing the tipping bucket type rainfall sensor, and measuring corresponding y values under each average rainfall intensity value to obtain n groups of test values;
step S13: substituting n groups of test values into a parameter-rain intensity function to be determined, calculating to obtain n fitting parameters, and determining an instant rainfall function according to the n fitting parameters: j (h, t) ═ 0;
step S14: based on the determined instant rainfall function, testing the tipping bucket type rainfall sensor by adopting the same test conditions as those in the step S12, obtaining switching value z output by the tipping bucket type rainfall sensor according to the measured instant rainfall h and the precision epsilon of the tipping bucket type rainfall sensor, and under each average rainfall intensity value, measuring the corresponding value of the switching value z to obtain n switching values z;
step S15: comparing the n switching values z with the ideal switching value D to obtain n correction values, wherein:
epsilon is the precision of the tipping bucket rainfall sensor, if a certain z value falls within the error range determined by the ideal switching value D, the corresponding correction value is 0; when all the correction values are 0, finishing construction; as long as one correction value is not zero, the execution continues to step S16;
step S16: correcting the y value in the n groups of test values by adopting the n corrected values respectively to obtain n groups of corrected test values;
step S17: substituting the corrected n groups of test values into a parameter-rain intensity function to be determined, calculating to obtain n fitting parameters, and determining the parameter-rain intensity function again according to the n fitting parameters: f (u, y, t) is 0 and the instant rainfall function: j (h, t) ═ 0; and returns to step S14.
5. The method for obtaining accurate rainfall data based on scoop parameters and rain intensity according to any one of claims 1-4, further comprising:
step S2: in the continuous period of rainfall, obtaining instant rainfall h according to the instant rainfall function and the value of the bucket duration t;
step S3: and obtaining the total rainfall H according to the values of all the instant rainfall H in the total rainfall statistical period.
6. The method for obtaining accurate rainfall data based on scoop parameters and rain intensity of claim 5, wherein: the tipping bucket rainfall sensor comprises: the water bearing device comprises a water bearing device, a funnel arranged below the water bearing device, a tipping bucket arranged below the funnel, magnetic steel driven by the tipping bucket, and a counting circuit module comprising a reed switch; the bucket time length t is the time interval between two adjacent switch pulses generated by the reed switch.
7. The method for obtaining accurate rainfall data based on scoop parameters and rain intensity of claim 6, wherein: the counting circuit module is provided with a clock circuit; the bucket time t is obtained through the time interval of two adjacent times of trigger of the reed switch, and an accurate value is calculated through a clock signal provided by a clock circuit arranged in the counting circuit module.
8. The method for obtaining precise rainfall data based on bucket parameters and rainfall intensity of claim 7 wherein in step S3, the total rainfall H is obtained by calculating the number of switching values outputted by the dump bucket rainfall sensor received by the monitor terminal, the switching values being obtained by the instantaneous rainfall H and the precision e of the dump bucket rainfall sensor; the switching value is obtained by dividing the accumulated value of the instant rainfall h by the precision epsilon of the tipping bucket rainfall sensor or by dividing the accumulated value of the instant rainfall h by the precision epsilon of the tipping bucket rainfall sensor.
9. The method for obtaining accurate rainfall data based on scoop parameters and rain intensity of claim 6, wherein: in step S2: and electronic signals generated by opening and closing the reed switch each time are transmitted to the monitoring terminal in real time, and the value of the bucket time t is obtained by calculation at the monitoring terminal.
10. The method for obtaining accurate rainfall data based on scoop parameters and rain intensity of claim 1, wherein: after step S1, according to the instant rainfall function and
obtaining an equivalent rain intensity function:
G(u,t)＝0；
wherein the unit of the equivalent rain intensity u of the bucket time t is mm/min.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811606261.2A CN109444995B (en) | 2018-12-26 | 2018-12-26 | Method for obtaining accurate rainfall data based on bucket parameters and rainfall intensity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811606261.2A CN109444995B (en) | 2018-12-26 | 2018-12-26 | Method for obtaining accurate rainfall data based on bucket parameters and rainfall intensity |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109444995A CN109444995A (en) | 2019-03-08 |
| CN109444995B true CN109444995B (en) | 2021-02-02 |
Family
ID=65537510
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811606261.2A Active CN109444995B (en) | 2018-12-26 | 2018-12-26 | Method for obtaining accurate rainfall data based on bucket parameters and rainfall intensity |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109444995B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116226606B (en) * | 2023-05-05 | 2023-07-04 | 安徽省大气探测技术保障中心 | Unmanned automatic weather station rainfall data quality control method and system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004198366A (en) * | 2002-12-20 | 2004-07-15 | Yokogawa Denshikiki Co Ltd | Rain gauge |
| CN108873108A (en) * | 2018-05-08 | 2018-11-23 | 河海大学 | One kind can trace to the source tilting bucket rain and its measurement method |
-
2018
- 2018-12-26 CN CN201811606261.2A patent/CN109444995B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004198366A (en) * | 2002-12-20 | 2004-07-15 | Yokogawa Denshikiki Co Ltd | Rain gauge |
| CN108873108A (en) * | 2018-05-08 | 2018-11-23 | 河海大学 | One kind can trace to the source tilting bucket rain and its measurement method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109444995A (en) | 2019-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109407184B (en) | Bucket-tipping type rainfall sensor correction method based on bucket parameters and rainfall intensity | |
| CN109507760B (en) | Method for obtaining accurate rainfall data based on bucket number and bucket duration | |
| CN109444995B (en) | Method for obtaining accurate rainfall data based on bucket parameters and rainfall intensity | |
| CN109444996B (en) | Method and device for obtaining accurate rainfall data | |
| CN102104259A (en) | Electric quantity detecting and controlling method and device for rechargeable battery | |
| CN104154987B (en) | A kind of vibrating string type sensor plate reading | |
| US4965756A (en) | Method and apparatus for calibration of electronic gas meters | |
| CN209432825U (en) | A kind of list barrel type runoff plots silt automated watch-keeping facility | |
| CN100568023C (en) | Intelligence testing apparatus for precipitation rain fall based on computer vision | |
| CN102176112A (en) | Method for achieving accurate clock timing by arranging RTC (real time clock) in MCU (micro control unit) | |
| CN109407182B (en) | High-precision tipping bucket type rainfall sensor based on bucket duration | |
| CN109507761B (en) | Method for obtaining accurate rainfall data based on equivalent rainfall intensity | |
| CN105223632A (en) | A kind of size double skip rainfall takes into account the method improving rainfall measurement accuracy | |
| CN110196395A (en) | Battery SOC estimation method | |
| US4953386A (en) | Method and apparatus for proving electronic gas meters | |
| US6339315B1 (en) | Measuring counter of the state of charge of the powering battery of an electronic appliance | |
| CN207528947U (en) | Rainfall titrates tester | |
| CN109407183B (en) | Tipping bucket type rainfall sensor for eliminating dynamic loss | |
| CN104809762A (en) | Multifunctional taximeter verification system | |
| CN113945774A (en) | Powder static accumulation amount testing equipment and testing method | |
| CN105954485B (en) | A kind of circuit with monitoring TDS water quality functions | |
| CN203271718U (en) | Tipping bucket flowmeter | |
| CN209446806U (en) | Rainfall gauge | |
| CN201397233Y (en) | Three-line pendulum oscillatory cycle measuring instrument | |
| CN103017847B (en) | The signal transmission of vortex electronic type water meter and heat energy meter and calculating device and method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |