CN109284724B - Filtering calculation method for vertical water temperature monitoring data - Google Patents
Filtering calculation method for vertical water temperature monitoring data Download PDFInfo
- Publication number
- CN109284724B CN109284724B CN201811154693.4A CN201811154693A CN109284724B CN 109284724 B CN109284724 B CN 109284724B CN 201811154693 A CN201811154693 A CN 201811154693A CN 109284724 B CN109284724 B CN 109284724B
- Authority
- CN
- China
- Prior art keywords
- water temperature
- layer
- envelope
- minimum
- curve
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Abstract
The invention discloses a filtering calculation method of vertical water temperature monitoring data, which comprises the following steps of: A. acquiring a minimum envelope curve of the vertical water temperature measured data; B. acquiring a maximum envelope curve of the vertical water temperature measured data; acquiring a water temperature correction curve according to the minimum envelope line and the maximum envelope line; D. acquiring the minimum envelope line of the water temperature correction curve in the step C as the updated minimum envelope line; E. c, acquiring the maximum envelope curve of the water temperature correction curve in the step C as the updated maximum envelope curve; F. circularly executing the steps C to E at least once; G. and calculating the water temperature monitoring data after filtering correction according to the updated minimum envelope line and the maximum envelope line. The invention can eliminate abnormal fluctuation of isothermal layer and other data, ensure that the vertical water temperature monitoring data after filtering calculation is monotonically increased, retain thermocline characteristics and restore the actual water temperature distribution characteristics.
Description
Technical Field
The invention belongs to the field of correction of vertical water temperature monitoring data in the field of engineering, and particularly relates to a filtering calculation method of vertical water temperature monitoring data.
Background
The vertical water temperature monitoring is an important link in the water temperature monitoring work, and the measured vertical water temperature distribution data is also important data for analyzing the water temperature structure of the reservoir.
The existing vertical water temperature monitoring instruments can be divided into a plurality of types such as temperature chains, optical fibers and the like. Limited by the stability and the measurement precision of the measurement data of the monitoring instrument, unreasonable data fluctuation often occurs in the vertical water temperature monitoring data, even under the condition that the water temperature exceeds 4 ℃, an inverse temperature layer with the height of 10m appears, which is not consistent with the actual situation, so that the originally obtained vertical water temperature monitoring data needs to be filtered and calculated to be close to the actual situation.
The commonly used filtering algorithm is a finite amplitude filtering algorithm, a median filtering algorithm, an arithmetic mean filtering algorithm, a weighted mean filtering algorithm and the like. However, these filtering algorithms are not suitable for correcting the vertical water temperature monitoring data. Because under normal conditions, when the water temperature is higher than 4 ℃, the temperature inversion layer does not appear in the reservoir water temperature, namely, the water temperature distribution is an increasing function. Meanwhile, the distribution curve of the water temperature along the elevation is complex in form, and has isothermal distribution, single thermocline distribution, double thermocline distribution and the like, the characteristics of the water temperature distribution curve are difficult to describe by common mathematical equations, and common filtering algorithms cannot ensure that the corrected water temperature monitoring data is monotonically increased and the distribution characteristics of the thermocline are reserved.
However, with the arrival of the big data era, more and more reservoirs start to implement vertical water temperature monitoring, and filtering correction on vertical water temperature monitoring data correctly and efficiently has a vital significance on data analysis work such as reservoir water temperature structure evaluation, reservoir water temperature structure change research and the like.
Disclosure of Invention
The invention aims to provide a filtering calculation method for vertical water temperature monitoring data, which can eliminate abnormal fluctuation of data such as an inverse temperature layer and the like, ensure that the vertical water temperature monitoring data after filtering calculation is monotonically increased, retain thermocline characteristics and restore actual water temperature distribution characteristics.
In order to solve the technical problems, the technical scheme adopted by the invention is as follows:
a filtering calculation method for vertical water temperature monitoring data is characterized by comprising the following steps:
a, acquiring a minimum envelope line of vertical water temperature measured data;
b, acquiring a maximum envelope line of the vertical water temperature measured data;
c, acquiring a water temperature correction curve according to the minimum envelope line and the maximum envelope line;
step D, acquiring the minimum envelope of the water temperature correction curve in the step C as the updated minimum envelope;
step E, acquiring the maximum envelope line of the water temperature correction curve in the step C as the updated maximum envelope line;
f, circularly executing the steps C to E at least once;
and G, calculating the water temperature monitoring data after filtering correction according to the updated minimum envelope line and the maximum envelope line.
By the method, two monotonically increasing minimum envelope lines and maximum envelope lines on the vertical water temperature measured data are found, so that the data after filtering correction are ensured to be monotonically increased, the envelope lines and the fluctuation of the flattened data are further shrunk on the basis, a filtering calculation result of the vertical water temperature monitoring data is obtained, abnormal fluctuation of data such as an inverse temperature layer can be eliminated, the vertically increasing water temperature monitoring data after filtering calculation are ensured, the thermocline characteristic is reserved, and the actual water temperature distribution characteristic is restored.
As a preferable mode, in step a, the minimum envelope of the measured water temperature data in the vertical direction is obtained as the minimum envelope t of the layer 1min1;
In the step B, the maximum envelope curve of the vertical water temperature measured data is obtained and used as the maximum envelope curve t of the 1 st layermax1;
In step C, according to tmin1And tmax1Acquiring a water temperature correction curve as a 1 st water temperature correction curve tavg1;
In step D, t is obtainedavg1Is taken as the layer 2 minimum envelope tmin2;
In step E, t is obtainedavg1As the maximum envelope t of layer 2max2;
In the step F, the steps C to E are executed circularly once, and the method comprises the following steps:
step F1, according to tmin2And tmax2Acquiring a water temperature correction curve as a 2 nd water temperature correction curve tavg2;
Step F2, obtaining tavg2Is taken as the layer 3 minimum envelope tmin3;
Step F3, obtaining tavg2As the maximum envelope t of layer 3max3;
In step G, according to tmin3And tmax3And calculating the filtered and corrected water temperature monitoring data.
Preferably, in step a, the layer 1 minimum envelope t ismin1Points onIs determined by the following formula:
in the formula, n is the number of vertical water temperature sampling points; i is the serial number of a vertical water temperature sampling point, and i is 1,2, …, n-1, n; t is tiIs the measured value of the water temperature of the ith sampling point, tiCorresponding water level height of hi,hi<hi+1;
In the step B, the maximum envelope t of the 1 st layermax1Points onIs determined by the following formula:
in the step C, the 1 st water temperature correction curve tavg1Points onIs determined by the following formula:
in the step D, the minimum envelope t of the layer 2min2Upper each nodeThe determination method comprises the following steps: if it isThenIf it isThen
In the step E, the maximum envelope t of the layer 2max2Upper each nodeThe determination method comprises the following steps: if it isThenIf it isThen
In the step F1, the 2 nd time water temperature correction curve tavg2Points onIs determined by the following formula:
in the step F2, the 2 nd time water temperature correction koji is obtainedLine tavg2Second derivative of (2)Will be provided withThe corresponding point is taken as the minimum envelope t of the layer 3min3Node (a) of
In the step F3, the 2 nd time water temperature correction curve t is obtainedavg2Second derivative of (2)Will be provided withThe corresponding point is taken as the maximum envelope t of the layer 3max3Node (a) of
In the step G, the water temperature monitoring data after the filtering correctionThe following equation is obtained:
preferably, the layer 2 minimum envelope tmin2Maximum envelope t of layer 2 at the value between adjacent nodesmax2Value between upper adjacent nodes, layer 3 minimum envelope tmin3Maximum envelope t of layer 3, the value between upper adjacent nodesmax3The values between the upper adjacent nodes are determined by linear interpolation.
Compared with the prior art, the invention finds the two monotonically increasing minimum envelope lines and the maximum envelope lines on the vertical water temperature measured data, thereby ensuring that the data after filtering correction is monotonically increased, further shrinking the envelope lines and flattening the data fluctuation on the basis, obtaining the filtering calculation result of the vertical water temperature monitoring data, eliminating abnormal fluctuation of data such as an adverse temperature layer and the like, ensuring that the vertical water temperature monitoring data after filtering calculation is monotonically increased and keeps the thermocline characteristic, reducing the actual water temperature distribution characteristic, reasonably and objectively correcting the vertical water temperature monitoring data, and laying a solid data foundation for deep analysis of the vertical water temperature monitoring data.
Drawings
FIG. 1 is a flow chart of an embodiment of the present invention.
FIG. 2 shows measured data of water temperature in vertical direction and the minimum envelope t of layer 1min1Layer 1 maximum envelope tmax1A graph of (a).
FIG. 3 shows the minimum envelope t of layer 1min1Layer 1 maximum envelope tmax1The 1 st water temperature correction curve tavg1A graph of (a).
FIG. 4 shows the minimum envelope t of layer 1min1Layer 2 minimum envelope tmin2Layer 1 maximum envelope tmax1Layer 2 maximum envelope tmax2The 1 st water temperature correction curve tavg1The 2 nd time water temperature correction curve tavg2A graph of (a).
FIG. 5 shows the minimum envelope t of layer 2min2Layer 2 maximum envelope tmax2The 2 nd time water temperature correction curve tavg2A graph of (a).
FIG. 6 shows the minimum envelope t of layer 3min3Layer 3 maximum envelope tmax3The 2 nd time water temperature correction curve tavg2A graph of (a).
FIG. 7 shows the minimum envelope t of layer 3min3Layer 3 maximum envelope tmax3And the water temperature monitoring data t after filtering correctionavg3A graph of (a).
FIG. 8 is a comparison graph of measured vertical water temperature data and filtered and corrected water temperature monitoring data.
Detailed Description
The specific implementation mode of the invention is introduced by taking the fiber monitoring data filtering calculation of the vertical water temperature in front of the Xiluodi dam as an example.
As shown in fig. 1, the present invention comprises the steps of:
step A, acquiring a minimum envelope of vertical water temperature measured data as a layer 1 minimum envelope tmin1。
Specifically, n is the number of vertical water temperature sampling points; i is the serial number of a vertical water temperature sampling point, and i is 1,2, …, n-1, n; t is tiIs the measured value of the water temperature of the ith sampling point, tiCorresponding water level height of hi,hi<hi+1. Calculating the minimum envelope t of the layer 1 from the (n-1) th measuring point to the 1 st measuring point along the descending elevation sequencemin1Layer 1 minimum envelope tmin1Points onIs determined by the following formula:
b, acquiring the maximum envelope curve of the vertical water temperature measured data as the maximum envelope curve t of the 1 st layermax1。
Specifically, the maximum envelope t of the layer 1 from the 2 nd measuring point to the nth measuring point is calculated along the ascending order of elevationmax1Maximum envelope t of layer 1max1Points onIs determined by the following formula:
layer 1 minimum envelope tmin1Layer 1 maximum envelope tmax1The calculation results are shown in fig. 2.
C, calculating each elevation h according to the step A and the step BiCorresponding toAndcalculating the mean value as the 1 st water temperature correction curve tavg11 st corrected value of water temperature at each elevationThe calculation results are shown in FIG. 3, the 1 st water temperature correction curve tavg1Points onIs determined by the following formula:
the 1 st water temperature correction curve tavg1The calculation results are shown in fig. 3.
Step D, obtaining tavg1Is taken as the layer 2 minimum envelope tmin2。
In particular, layer 2 minimum envelope tmin2Upper each nodeThe determination method comprises the following steps: for theDescending along the elevation from the n-1 measuring point to the 1 measuring point ifThenIf it isThen
Step E, obtaining tavg1As the maximum envelope t of layer 2max2。
In particular, the layer 2 maximum envelope tmax2Upper each nodeThe determination method comprises the following steps: for theAscending along the elevation, from the 2 nd measuring point to the nth measuring point, ifThenIf it isThen
Layer 2 minimum envelope tmin2Layer 2 maximum envelope tmax2The calculation results are shown in fig. 4.
Step F1, according to tmin2And tmax2Acquiring a water temperature correction curve as a 2 nd water temperature correction curve tavg2。
Specifically, calculatingAndcorresponding mean valueAndand calculateEach elevation hiCorresponding toAndthe weighted mean value of (2) is used as the second water temperature correction curve
2 nd water temperature correction curve t of each elevationavg2The calculation results are shown in fig. 5.
Step F2, obtaining tavg2Is taken as the layer 3 minimum envelope tmin3. Specifically, the 2 nd water temperature correction curve t is obtainedavg2Second derivative of (2)Will be provided withThe corresponding point is taken as the minimum envelope t of the layer 3min3Node t ofi min3;
Step F3, obtaining tavg2As the maximum envelope t of layer 3max3. Specifically, the 2 nd water temperature correction curve t is obtainedavg2Second derivative of (2)Will be provided withThe corresponding point is taken as the maximum envelope t of the layer 3max3Node (a) of
Layer 3 minimum envelope tmin3Layer 3 maximum envelope tmax3The calculation results are shown in fig. 6.
Step G. according to tmin3And tmax3And calculating the filtered and corrected water temperature monitoring data.
In particular, from each elevation hiCorresponding toCalculating the mean value as the 3 rd water temperature correction curve tavg3And correcting the 3 rd water temperature by the curve tavg3And the water temperature is used as a water temperature monitoring data curve after final filtering correction.
Water temperature monitoring data corrected by filtering each elevationThe following equation is obtained:
water temperature monitoring data t after filtering correctionavg3The calculation results are shown in fig. 7.
In this embodiment, the layer 2 minimum envelope tmin2Maximum envelope t of layer 2 at the value between adjacent nodesmax2Value between upper adjacent nodes, layer 3 minimum envelope tmin3Maximum envelope t of layer 3, the value between upper adjacent nodesmax3The values between the upper adjacent nodes are determined by linear interpolation.
As shown in FIG. 8, the invention eliminates abnormal fluctuation of data such as an inverse temperature layer, and the like, the vertical water temperature monitoring data after filtering calculation is monotonically increased, and the thermocline characteristic is kept, so that the actual water temperature distribution characteristic is basically restored.
While the present invention has been described with reference to the embodiments shown in the drawings, the present invention is not limited to the embodiments, which are illustrative and not restrictive, and it will be apparent to those skilled in the art that various changes and modifications can be made therein without departing from the spirit and scope of the invention as defined in the appended claims.
Claims (4)
1. A filtering calculation method for vertical water temperature monitoring data is characterized by comprising the following steps:
a, acquiring a minimum envelope line of vertical water temperature measured data;
b, acquiring a maximum envelope line of the vertical water temperature measured data;
c, acquiring a water temperature correction curve according to the minimum envelope line and the maximum envelope line;
step D, acquiring the minimum envelope of the water temperature correction curve in the step C as the updated minimum envelope;
step E, acquiring the maximum envelope line of the water temperature correction curve in the step C as the updated maximum envelope line;
f, circularly executing the steps C to E at least once;
and G, calculating the water temperature monitoring data after filtering correction according to the updated minimum envelope line and the maximum envelope line.
2. The method of claim 1, wherein the vertical water temperature monitoring data is filtered,
in the step A, the minimum envelope curve of the vertical water temperature measured data is obtained and used as the minimum envelope curve t of the layer 1min1;
In the step B, the maximum envelope curve of the vertical water temperature measured data is obtained and used as the maximum envelope curve t of the 1 st layermax1;
In step C, according to tmin1And tmax1Acquiring a water temperature correction curve as a 1 st water temperature correction curve tavg1;
In step D, t is obtainedavg1Is taken as the layer 2 minimum envelope tmin2;
In step E, t is obtainedavg1As the maximum envelope t of layer 2max2;
In the step F, the steps C to E are executed circularly once, and the method comprises the following steps:
step F1, according to tmin2And tmax2Acquiring a water temperature correction curve as a 2 nd water temperature correction curve tavg2;
Step F2, obtaining tavg2Is taken as the layer 3 minimum envelope tmin3;
Step F3, obtaining tavg2As the maximum envelope t of layer 3max3;
In step G, according to tmin3And tmax3And calculating the filtered and corrected water temperature monitoring data.
3. The method of claim 2, wherein the vertical water temperature monitoring data is filtered,
in the formula, n is the number of vertical water temperature sampling points; i is the serial number of a vertical water temperature sampling point, and i is 1,2, …, n-1, n; t is tiIs the measured value of the water temperature of the ith sampling point, tiCorresponding water level height of hi,hi<hi+1;
In the step B, the maximum envelope t of the 1 st layermax1Points onIs determined by the following formula:
in the step C, the 1 st water temperature correction curve tavg1Points onIs determined by the following formula:
in the step D, the minimum envelope t of the layer 2min2Upper each nodeThe determination method comprises the following steps: if it isAnd isThenIf it isThen
In the step E, the maximum envelope t of the layer 2max2Upper each nodeThe determination method comprises the following steps: if it isAnd isThenIf it isThen
In the step F1, the 2 nd time water temperature correction curve tavg2Points onIs determined by the following formula:
in the step F2, the 2 nd time water temperature correction curve t is obtainedavg2Second derivative of (2)Will be provided withThe corresponding point is taken as the minimum envelope t of the layer 3min3Node (a) of
In the step F3, the 2 nd time water temperature correction curve t is obtainedavg2Second derivative of (2)Will be provided withThe corresponding point is taken as the maximum envelope t of the layer 3max3Node (a) of
In the step G, the water temperature monitoring data after the filtering correctionThe following equation is obtained:
4. the method of claim 3, wherein the layer 2 minimum envelope t is calculated by filtering the vertical water temperature monitoring datamin2Maximum envelope t of layer 2 at the value between adjacent nodesmax2Value between upper adjacent nodes, layer 3 minimum envelope tmin3Maximum envelope t of layer 3, the value between upper adjacent nodesmax3The values between the upper adjacent nodes are determined by linear interpolation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811154693.4A CN109284724B (en) | 2018-09-30 | 2018-09-30 | Filtering calculation method for vertical water temperature monitoring data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811154693.4A CN109284724B (en) | 2018-09-30 | 2018-09-30 | Filtering calculation method for vertical water temperature monitoring data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109284724A CN109284724A (en) | 2019-01-29 |
CN109284724B true CN109284724B (en) | 2021-09-10 |
Family
ID=65182508
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811154693.4A Active CN109284724B (en) | 2018-09-30 | 2018-09-30 | Filtering calculation method for vertical water temperature monitoring data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109284724B (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5516276A (en) * | 1978-07-21 | 1980-02-04 | Matsushita Electric Ind Co Ltd | Vertical water temperature distribution measuring method |
CN103162869A (en) * | 2013-02-05 | 2013-06-19 | 中国长江三峡集团公司 | Measuring method of deepwater reservoir vertical direction water temperature distribution |
-
2018
- 2018-09-30 CN CN201811154693.4A patent/CN109284724B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5516276A (en) * | 1978-07-21 | 1980-02-04 | Matsushita Electric Ind Co Ltd | Vertical water temperature distribution measuring method |
CN103162869A (en) * | 2013-02-05 | 2013-06-19 | 中国长江三峡集团公司 | Measuring method of deepwater reservoir vertical direction water temperature distribution |
Also Published As
Publication number | Publication date |
---|---|
CN109284724A (en) | 2019-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108108530B (en) | Fatigue life calibration method suitable for structural connecting piece | |
CN112100713B (en) | Cable force automatic identification method based on variable-pitch grid | |
CN109284724B (en) | Filtering calculation method for vertical water temperature monitoring data | |
JP2012159414A (en) | Battery status estimation device | |
CN114116689A (en) | Big data cleaning method based on building structure safety monitoring | |
WO2015097795A1 (en) | Wind-speed evaluation method, power generation amount estimation method, output and voltage estimation method, placement method for wind power generator, and wind-speed evaluation device | |
CN110866682B (en) | Underground cable early warning method and device based on historical data | |
CN106897789B (en) | Wind power active power climbing event detection method and device | |
CN111667550A (en) | Test method for judging quality of fitting curve | |
CN103116877B (en) | A kind of reservoir level process smoothing processing method | |
CN117332205A (en) | High-precision automatic optimization method and device for temperature compensation resistance of piezoresistor | |
CN111475972B (en) | High-temperature fatigue-creep interaction damage assessment method, system and storage medium | |
CN111238667B (en) | Temperature compensation method, printed circuit board, compressor and vehicle | |
CN111368459A (en) | Vibration fatigue evaluation method and device for wind power generation supporting tower, medium and terminal | |
CN116150991A (en) | Transmission line tower unbalance tension calculation method based on laser point cloud | |
CN113625104B (en) | Line length correction method for traveling wave fault location | |
CN112685449A (en) | Thermodynamic system performance calculation method, device and equipment | |
CN114235275A (en) | Calibration method and device of pressure sensor and electronic equipment | |
CN108345573B (en) | Differential expansion determining function calculation method for differential expansion measuring probe of high-pressure cylinder of steam turbine | |
JP6011650B2 (en) | Wind noise evaluation method for overhead power transmission lines | |
CN104101507B (en) | The data processing equipment of a kind of steam turbine traditional performance parameter and method | |
CN108599976B (en) | Network quality evaluation model construction method and device and computer equipment | |
CN109960818B (en) | Method and device for generating simulated wind speed data of wind power plant | |
CN117313014B (en) | Real-time monitoring and early warning method for abnormal energy consumption data of kiln operation | |
CN111950123A (en) | Gyroscope error coefficient curve fitting prediction method and system |
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 |