CN107679021A - A kind of river Ru Hu storehouses method of calculating flux - Google Patents
A kind of river Ru Hu storehouses method of calculating flux Download PDFInfo
- Publication number
- CN107679021A CN107679021A CN201710892965.XA CN201710892965A CN107679021A CN 107679021 A CN107679021 A CN 107679021A CN 201710892965 A CN201710892965 A CN 201710892965A CN 107679021 A CN107679021 A CN 107679021A
- Authority
- CN
- China
- Prior art keywords
- day
- river
- water
- lake
- rainfall
- 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.)
- Granted
Links
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Life Sciences & Earth Sciences (AREA)
- Algebra (AREA)
- Remote Sensing (AREA)
- Software Systems (AREA)
- Sewage (AREA)
Abstract
The present invention relates to a kind of river Ru Hu storehouses method of calculating flux, mainly enter lake storage hydraulic condition of river data and meteorological data by collecting, calculate the coefficient correlation of cross section of river flow and different precipitation duration, select the precipitation duration T significantly correlated with section flow, establish representative cross section of river flow and select total precipitation relational expression in precipitation duration T, the aquifer yield of the total precipitation of unit area in precipitation duration is calculated, water collection basin domain area is recycled, calculates each river course and enter lake storehouse flow day by day;The unit area Rainfall-runoff coefficient that the method for calculating flux of the present invention is also related to based on principle of water balance to calculating is adjusted, lifting flow, the harmony of waterlevel data.River Ru Hu storehouses flow rate calculation is carried out using the method for the present invention, improves the time into lake storehouse water and spatial resolution, enters the calculating of lake pollution flux for the main river courses of Hu Ku and estimation is laid a good foundation.
Description
Technical field
The invention belongs to lake storehouse water environment technical field, and in particular to a kind of Hu Ku rivers Ru Hu storehouses method of calculating flux.
Background technology
Data on flows is to calculate lake storehouse water balance, flood control, water supply and the Important Parameters of environmental quality management, accurate estimation
The main river waters of Hu Ku are significant with calculating.Acquisition river course enters lake Incoming water quantity and usually requires to enter lake storage in river course
Mouth door lays water quantity monitoring section, is obtained by observation procedure.Because river water is by upstream rainfall and lake reservoir level height
Influence, enter the lake storage discharge of river and quick time change be present, therefore, it is accurate observe river course enter lake Incoming water quantity need into
Row high frequency monitors.The discharge of river is so obtained using the method for monitoring, not only needs the fixed infrastructure of construction, instrument to set
It is standby, it is also necessary to be equipped with enough professionals.Enter lake storage river course it is numerous in the case of, the acquisition of lake storehouse calculated flow into lake at present
It is general that using setting measuring frequency section to carry out cross section of river flow in river mouth, once and several times observation procedure obtains.This method obtains
Lake storehouse water must be entered day needs to carry out big discharge observation, it is necessary to expend a large amount of manpower and financial resources.Therefore, carry out to be put in storage into lake at present and lead
The Hu Ku for wanting river water all standing to monitor is seldom, and most of Hu Ku only sets a representative station in main stream entering lake and carries out water
Test, then using represent estimation moon yardstick of standing, the whole Hu Ku collection waters of year yardstick enters water, not only lack day yardstick and enter lake
Incoming water quantity, and more lack each mainly to enter lake storage river course and enter lake Incoming water quantity.This is unfavorable for recognizing lake storehouse change of water quality
Reason and grasp pollutant sources.Particularly current, lake storehouse water environment protection and administer under the situation being highly valued, lack
The weary main river courses of Hu Ku enter water, allow for related work and lack scientific basis, it is difficult to concentrate on crucial points, work effect needs to be carried
It is high.
The content of the invention
It is an object of the invention to provide a kind of river Ru Hu storehouses stream for fully excavating Hu Ku collection waters hydrographic information data
Measure computational methods.
To realize above-mentioned technical purpose, the present invention adopts the following technical scheme that:
Step 1:Lack flow monitoring data problem for entering Hu Ku rivers river mouth section, collect Hu Ku and its River Basin Hydrology
Data and meteorological data, the hydrographic data include Hu Ku day by day water level, in basin control hydrometric station day by day flow, enter Hu Ku master
Want river basins area features, using Hu Ku as sewage treatment plant tail water discharge in water factory's water withdrawal at water source, Hu Ku collection waters and
Lake storehouse water outlet data on flows day by day;The meteorological data includes Hu Ku collection waters rainfall data, lake storehouse water surface evaporation number day by day
According to;
Wherein, it is described enter lake storehouse trunk river refer to the collection waters gross area in collected river and account for all to enter Hu Ku rivers and catchment
More than the 90% of the domain gross area;
Step 2:According to Hu Ku rivers collection waters internal control hydrometric station data on flows and rainfall data day by day day by day are entered, calculate
Hydrometric station control section day by day flow it is corresponding collection waters before 1~first 10 days different duration T the total precipitation between
Pearson correlation coefficient, the size identification according to coefficient go out and the section significantly correlated crucial duration T of flow day by day;
Step 3:Using the crucial duration T identified in step 2, the flow day by day surveyed according to hydrometric station control section,
Corresponding crucial duration T quantum of rainfall, upstream collection water surface area, using linear regression least square method, establish the control hydrology
Total rainfall linear relation in measured discharge of standing crucial duration T corresponding with collection waters;With the catchment area of control of section remove flow-
Total precipitation relational expression in T durations, tries to achieve basin unit area Rainfall-runoff coefficient;
Step 4:According to hydrometric station control section flow~rainfall of foundation, collection water surface area, runoff coefficient relation, meter
The flow day by day for respectively entering Hu Ku rivers is calculated, rainfall is obtained and enters Hu Ku total Waters day by day;According to lake reservoir level, lake reservoir level area
Curve calculates water level diurnal variation caused by rainfall runoff enters lake storehouse, obtains lake reservoir level diurnal variation sequence 1 caused by rainfall;
Step 5:According to the tail water for the discharge of sewage disposal day by day being collected into, water factory's water withdrawal, lake storehouse water surface evaporation, lake
Storehouse water surface rainfall, lake storehouse water yield and lake reservoir level data, based on principle of water balance, calculate is caused by rainfall runoff
The Day-to-day variability of lake reservoir level, obtain lake reservoir level diurnal variation sequence 2 caused by rainfall;
Step 6:Carry out the comparison of lake reservoir level diurnal variation sequence 2 and the change curve of lake reservoir level diurnal variation sequence 1, adjustment
Participate in respectively entering the runoff coefficient of the flow rate calculation day by day in Hu Ku rivers so that the change curve of lake reservoir level diurnal variation sequence 1 approaches
Lake reservoir level diurnal variation sequence 2, each day runoff coefficient after being corrected;
Step 7:The flow in Hu Ku river courses is respectively entered according to day runoff coefficient calculating after amendment, then by sewage treatment plant tail water
Day by day discharge capacity is added to day by day on flow of corresponding river, calculates and obtains the flow that each river enters Hu Ku day by day.
After coefficient adjustment is confluxed in production, enter lake enter river flow data realize it is consistent with lake waterlevel data day by day,
The harmony of flow, waterlevel data is improved, minimum flow data are provided for the research of lake storehouse Hydrodynamic Process.
The method of the present invention, in the step 2, according to hydrometric station measured discharge, section upstream collection waters rainfall, calculate
Section daily average water discharge with three kinds of situations comprising the same day, not comprising the same day, not comprising the same day and proxima luce (prox. luc), time span 1
To the linearly related Pearson correlation coefficient of the total rainfall in gathering ground on the 10th and hydrometric station measured discharge, maximum Pierre is chosen
Inferior coefficient correlation corresponds to time span T, the duration of raining used for section flow rate calculation.
In the step 2, the collection larger reservoir situation of water surface area be present for entering lake storehouse upper river collection waters, calculate river
Reservoir catchment domain area is deducted during adfluxion water surface area from river collection water surface area.
The method of the present invention has the advantages that:
(1) time into lake storehouse water data and spatial resolution are improved, main river course data on flows day by day is given, is
The main river courses of Hu Ku enter lake contaminant capacity calculating and estimation is laid a good foundation;
(2) temporal resolution significantly improves, and current most of lake storehouse calculated flow into lake is only capable of being estimated to water month by month, without day by day
Data on flows;
(3) cause all kinds of hydrographic informations in basin are worth fully to be excavated so that the data on flows day by day of acquisition, water level number
According to the compatibility and harmony for maintaining height.
Brief description of the drawings
Fig. 1 is Chaohu Huan Hu river courses of embodiment 1 method of calculating flux flow chart day by day;
Fig. 2 is that the rainfall coefficient correlation change per day with 1-10 before gathering ground of 2014-2015 peaches small stream diurnal courses amount is rolled over
Line chart;
Fig. 3 is peach small stream daily average water discharge in 2014 and 1-10 days average rainfall scatter diagrams before gathering ground;
Fig. 4 is peach small stream daily average water discharge in 2015 and 1-10 days average rainfall scatter diagrams before gathering ground;
Fig. 5 is to calculate water level and measured water level process comparison diagram according to the discharge of river of primary Calculation;
Fig. 6 is that each river course inquired into using the runoff coefficient after correction enters Chaohu that lake total amount calculates 2014 and 2015
Water level process and the contrast situation schematic diagram of measured water level process;
Fig. 7 be finally inquire into obtained 12, Chaohu main stream entering lake 2014 with 2015 day by day runoff process illustrate
Figure.
Embodiment
By taking Chaohu as an example, the technical scheme for calculating inventive flow method is further described embodiment.
Chaohu is located in the middle part of Anhui Province, is the fourth-largest fresh water lake in China Yangtze river basin.Its drainage area 13486km2, Chaohu
Drainage area 9153km on lock2, category THE LOWER YANGTZE VALLEY left bank water system, Chaohu water surface area 769.6km2(ordinary water level 8.5m Wusong is high
Journey), mean depth about 2.7m, the m of volume 20.7 hundred million3.Basin physical features Xi Gaodong is low, middle low-lying flat.Chaohu Basin tributary is many
More, Heavenly Stems and Earthly Branches stream is mostly tree-shaped, and hydraulic connection is clear.Main ingoing river has:Hang Buhe-Feng Le rivers, send river, the southern Fei river in Anhui Province-shop port
(tables 1) such as river, the white stone Milky Way, million rivers, Zhe Gao rivers, ten Wulihes, Chiung smelt rivers, wherein Hang Buhe-Feng Le rivers be calculated flow into lake most
Big river, it is secondly the southern Fei river in Anhui Province and the Bai Shi Milky Way, the Yu Xi rivers of lake east are that Chaohu uniquely goes out lake and enters Jiang Tongdao, million river sections
Check gate is had, section flow is had data on flows day by day by manual control.
The Hefei office Chaohu Basin fork situation table of table 1
Calculated using the river Ru Hu storehouses method of calculating flux of the present invention and Chaohu day enter lake current amount, as shown in figure 1, specific step
It is rapid as follows:
Step 1:The main ingoing river in Chaohu is chosen, obtains the hydrographic data and meteorological data for mainly entering Hu Ku rivers;
For Chaohu, according to upper table, the main ingoing river in Chaohu includes Hang Buhe, the southern Fei river in Anhui Province, sends river, the white stone Milky Way, Zhe Gao rivers, ten
Wulihe, the pool Xihe River, doube bridge river, Chiung smelt river, chicken Yu He, Jiang Kouhe and Hua Tanghe;
Collect and obtain the more complete hydrology of Chaohu 2014-2015, meteorological data.Including:
Hydrographic data:Feng Le rivers Tao Xi hydrometric stations intraday effect data, Chao Huzha, million river lock datas on flows;Chao Hutang
West, The Zhong Temple, the average daily waterlevel data of Chaohu lock;Drainage area characteristic;Chaohu water factory day water withdrawal, subordinate county-level city of Hefei City
Sewage treatment plant's day output, Tail water reuse related data;
Meteorological data:The weather station such as Chaohu Basin Hefei, Tao Xi, Bai Shan daily rainfall data;Chao Huzha, Dong Pu station day water
Face evaporation capacity data.
Step 2:Utilize the Shan Nan in the Tao Xi hydrometric stations daily mean data on flows of 2014~2015 years and gathering ground, peach
Small stream, mill pier reservoir weather station daily rain amount data, analysis calculate the time lag relation between peach streams amount and basin rainfall.It is first
Three precipitation station arithmetic mean of instantaneous values are first calculated, as the per day rainfall in Tao Xi gathering grounds;Then peach small stream is calculated day by day respectively
Flow with comprising the same day, not comprising the same day and not comprising under three kinds of situations of the same day and proxima luce (prox. luc), time span is to catchment for 1st~10
The linearly related Pearson correlation coefficient of total rainfall is slided in area, and the Pearson correlation coefficient for choosing maximum corresponds to time span T,
To influence the rainfall key duration T of peach small stream diurnal courses amount.
Step 3:Total rainfall in the crucial duration T for influenceing Tao Xi hydrometric station daily mean flows using being identified in step 2
Amount, flow day by day, corresponding crucial duration T quantum of rainfall, upstream collection water surface area according to the actual measurement of hydrometric station control section,
Using the least square method in linear regression method, establish in hydrometric station measured discharge crucial duration T corresponding with collection waters and always drop
The rainfall of rain~flow linear relation;Total precipitation relational expression in flow-T durations is removed with the catchment area of control section, is tried to achieve
Basin unit area Rainfall-runoff coefficient;
In the present embodiment, using the crucial duration T identified in step 2, peach small stream day in 2015 in 2014 are established respectively
Rainfall~runoff linear relation of equal runoff and gathering ground prophase programming:
QPeach small stream, t=α PPeach small stream, t-1·SPeach small stream (1)
In formula:QPeach small stream, tFor peach stream text stand control section t days the average daily runoff volume (m of actual measurement3);α is that Feng Le rivers are single
Plane accumulates Rainfall-runoff coefficient;PPeach small stream, t-1For average total rainfall (m) in the T durations of Tao Xi control catchments face;SPeach small streamFor peach stream text
Stand control drainage area (㎡).Wherein QPeach small stream, t、PPeach small stream, t-1And SPeach small streamAll it is known quantity, using least square method, 2014 can be tried to achieve
Year, Tao Xi control catchments unit area Rainfall-runoff factor alpha in 2015.
Fig. 2 is 2014, the Pearson came before peach small stream diurnal courses amounts in 2015 and gathering ground between 1~10 per day rainfall
Coefficient correlation situation of change.As can be seen that either 2014 or 2015, peach small stream diurnal courses amount with not comprising before the same day
Phase rainfall correlation is relatively preferable.Under the situation, when duration of raining length and when extending to first 4 days within 1 day before runoff, and day by day
Run-off correlation is in rising trend, and on a declining curve after extending to 5 days, and what this feature embodied for 2015 more fills
Point, i.e. total rainfall of first 4 days and peach small stream diurnal courses correlation are best, and the coefficient correlation of the two reached for 2014 respectively with 2015
0.817 and 0.847, in clear-cut correlation, total rain fall on the 4th is the most notable for peach small stream station flow effect before explanation.
Fig. 3 and Fig. 4 be respectively before peach small stream daily average water discharges in 2014,2015 and gathering ground 1-10 days moving average rainfalls dissipate
Point diagram and coefficient correlation.It can be seen that the point distribution formed with preceding 2 days rainfalls and peach small stream daily flow for first 1 day is more
Scattered, correlation is poor.And preceding 4 days rainfalls and peach small stream daily flow are totally linearly distributed, correlation is best.
It was on the basis of the crucial effect period at 4 days before identifying, using formula (1), is entered respectively using least square method
The average preceding four daily rain amounts total amount of row was fitted with, 2015 annual diurnal courses total amounts in 2014, can be obtained to produce per year with 2015 for 2014 and be flowed
Coefficient, respectively 0.0917 and 0.1046.
For the southern Fei river in Anhui Province, other 13, river etc. is sent without Streamflow Data river course, it will be assumed that its unit area rainfall in period t
Runoff coefficient is equal with Feng Le basins, then the discharge of river can be calculated by following formula:
QI, t=α Pi,t-1·Si (2)
In formula:QI, tFor river course i t days calculated flow into lake preresearch estimates value;Pi,t-1Dropped for sub-basin before t days in T durations
Rain total amount, its calculation interval are identical with peach small stream;SiFor sub-basin collection water surface area, due to southern Fei river in Anhui Province basin upstream have concurrently Dong Pu,
The big room capital of the ancient state of Chu and many emerging reservoirs for waiting collection water surface area larger, rainfall runoff use the southern Fei river in Anhui Province largely by reservoir-trap, during calculating
Drainage area is that southern Fei river in Anhui Province drainage area deducts above-mentioned reservoir catchment domain area.
For the southern Fei river in Anhui Province, other 13, river etc. is sent without footpath flow data river course, the preceding four dates Basin Rainfall total amount by corresponding to
And gathering ground area substitutes into formula (2), and the preliminary day calculated flow into lake in each bar river course is calculated.
Step 4:Rainfall~runoff linear relationship, collection water surface area, the relation of runoff coefficient established according to step 3, profit
Calculated with formula (2) and respectively enter the flow day by day in Chaohu river, obtained rainfall and enter Chaohu total Water day by day, according to Chaohu water level,
The lake level curve of areas calculates water level diurnal variation caused by rainfall runoff enters Chaohu, and Chaohu water level day caused by obtaining rainfall becomes
Change sequence 1;
Fig. 5 be using primary Calculation each bar river course Ru lakes total amount, and consider Chao Huzha, million river locks, rainfall, evaporation,
Water intaking, sewage treatment plant tail water emission behaviour, be calculated according to principle of water balance Chaohu day water level process with actual measurement water
Position contrast situation.Wherein dark solid is Chaohu water level calculated using the discharge of river inquired into, the shade area of coverage be Tang Xi, in
Mausoleum and the gaging station measured water level constant interval of Chaohu lock three, light dotted line are three station mean waters.As can be seen that using inquiring into
Chaohu water level process that river course calculated flow into lake calculates and actual measurement change in process trend are basically identical, illustrate to inquire into entering for each bar river course
Lake water amount change procedure matches with actual conditions, but certain deviation be present in different periods, and measured water level, such as
2~July in 2014, the water level for inquiring into flow rate calculation are less than actual measurement process apparently higher than measured water level, while 9~October, this master
If because basin soil moisture content is a dynamic changing process, such as in dry season, because Basin Rainfall amount is less, soil contains
Water will be significantly lower than the wet season, if using the runoff coefficient as the wet season, the river course for causing to inquire into is entered into lake current amount height
In actual conditions.Therefore, dynamic need to be made to runoff coefficient in different periods and is adjusted according to actual conditions.
Step 5:According to Chaohu water surface rainfall, evaporation from water surface water, Chaohu water factory water withdrawal, sewage treatment plant of Hefei City
Tail water reuse amount, Chaohu water yield data, based on principle of water balance, that is, enter Chaohu river and enter water sum and Chaohu water storage water
The change of amount is equal, calculates the Day-to-day variability for causing Chaohu water level by rainfall runoff, and Chaohu water level day caused by obtaining rainfall becomes
Change sequence 2;
Step 6:Though each river course inquired into using fixed runoff coefficient enters Chaohu water level of lake flow rate calculation with surveying water
Position situation of change is consistent, but obvious deviation be present in partial period.Therefore consider according to Chaohu water balance to runoff coefficient
Amendment is made, i.e., for period t, if calculating water level is higher than measured water level, it is bigger than normal that the river course for illustrating to inquire into enters lake current amount, then reduces
Runoff coefficient;Conversely, if calculating water level be less than measured water level, illustrate that the flow inquired into is less than normal, then accordingly increase runoff coefficient, directly
It is basically identical to calculating water level and measured water level.
Up to above-mentioned target, to carry out the ratio of Chaohu water level diurnal variation sequence 2 and the change curve of Chaohu water level diurnal variation sequence 1
Compared with adjustment participates in the runoff coefficient of the flow rate calculation day by day of each Chaohu ingoing river so that Chaohu water level diurnal variation sequence 1 changes
Curve approaches Chaohu water level diurnal variation sequence 2, obtains each day runoff coefficient;
According to principle of water balance, Chaohu lake surface rainfall and evaporated water, water factory's water withdrawal, sewage treatment plant tail water discharge
Amount, Chaohu lock day outflow water, the water of Chaohu Lake amount that came in and went out million river lock days and other need to inquire into river course calculated flow into lake sum should be with nest
The change of lake water storage water is equal, i.e. draws formula (3) accordingly.Here water storage water is approximately equal to Chaohu water surface area and water
The product of position change.
Δlt·SC, t=∑ Qi,t+(Pt-Et)·SC, t-At+Tt (3)
In formula:ΔltSEA LEVEL VARIATION in the t days obtained is calculated by measured water level value for Chaohu;SC, tPut down for Chaohu within t days
Equal water surface area, Δ lt·SC, tChaohu moisture storage capacity variable quantity as in t days;PtIt is Chaohu in t daily rainfalls, using in the week of lake
The mean value calculation at 6 stations such as mausoleum, Chiung smelt, Bai Shan, pool west;EtFor in t day Chaohu evaporation capacity, using the flat of Dong Pu stations and Chaohu lock
Average;AtIt is Chaohu water factory in t day water withdrawals;TtTailwater quantity, ∑ Q are discharged in t days for sewage treatment plantI, tTo enter Chaohu river t
Day enters water sum.
Formula (2) is substituted into formula (3), the water level variable by going out the estimation of calculated flow into lake preresearch estimates value can be extrapolated.Enter
Row estimated value and actual measurement Chaohu water level variable contrast, expand or reduce flow anomaly coefficient, ensureing estimated value and actual measurement nest
Synchronously increase and decrease simultaneously, makes the two as close possible to reduction error to Lake shoal amount.Formula (2) is substituted into formula (3) to be derived from adjusting
Shown in unit area Rainfall-runoff coefficient calculation such as formula (4) after whole:
α=[Δ lt·SC, t-(Pt-Et)·SC, t+A-Tt]/∑(Pi,t-1·Si) (4)
Chaohu measuring point water level can be impacted in view of wind gushes increase and decrease water, during to ensure that measuring point water potential truly reflects Chaohu
Between mean water change, survey here Chaohu water level variable take The Zhong Temple, on the lock of Chaohu, the western three gaging station water level variables in the pool it is equal
Moving averages on the seven of value, water level is influenceed with rejecting wind field.Calculated again with the unit area Rainfall-runoff coefficient after adjustment
Estimate the day calculated flow into lake in river.Sewage treatment plant tail water is added on corresponding river course again after this, 12 can be obtained out
Day by day the discharge process in river course.
The runoff coefficient that table 2 is 2014 and day part in 2015 uses, wherein on January 1st, 2014 to March 19, July
Adjustment was made that to two period runoff coefficients of September 1 day in 30th, other mosts of the time do not make a change.As can be seen that adjust
After whole, runoff coefficient and annual rainfall distribution space feature are more consistent, i.e., dry season due to soil moisture content it is relatively low, production stream system
Number is smaller, and wet season basin soil moisture content is in saturation or the state close to saturation, and Rainfall-runoff coefficient is with increase.2015
Year is because annual rainfall distribution spatio-temporal difference is larger, therefore runoff coefficient change is also relative complex, wherein on January 1 to 2 months 15
Day, March 16 to May 2, June 28 to July 29 and September 4 days to runoff coefficient on December 31 are made that certain tune
It is whole.
The Basin Rainfall runoff coefficient of table 2 corrects
Fig. 6 is that each river course for being inquired into using the runoff coefficient after correction enters Chaohu that lake total amount calculates 2014 and 2015
Water level process and the contrast situation of measured water level process.As can be seen that after corrected, calculate water level and measured water level is basic
It coincide, only less deviation is present in partial period, the discharge of river for illustrating to inquire into is closer to actual conditions.
Fig. 7 be finally inquire into obtain Chaohu doube bridge river, Zhe Gao rivers, Chiung smelt river, 12 main stream entering lakes such as the southern Fei river in Anhui Province
2014 with runoff process day by day in 2015.As can be seen that in 12 river courses, Hangzhoupro port river flow is maximum, and the southern Fei river in Anhui Province is taken second place, the flower pool
The river course daily average water discharges such as river, Jiang Kou rivers are smaller, and million ruisseaus are mainly controlled by gate, and for flow by manual control, itself has flow
Data, it is not necessary to calculate.
Claims (7)
1. a kind of river Ru Hu storehouses method of calculating flux, it is characterised in that comprise the following steps:
Step 1:Hu Ku and its River Basin Hydrology data and meteorological data are collected, the hydrographic data includes Hu Ku water level, stream day by day
Controlled in domain hydrometric station day by day flow, enter lake storehouse trunk river drainage area feature, water factory's water withdrawal using Hu Ku as water source, lake
Sewage treatment plant tail water discharge and lake storehouse water outlet data on flows day by day in Ku Ji waters;The meteorological data includes Hu Ku collection waters
Day by day rainfall data, lake storehouse water surface evaporation data;
Wherein, it is described enter lake storehouse trunk river refer to the collection waters gross area in collected river and account for that all to enter Hu Ku rivers collection waters total
More than the 90% of area;
Step 2:According to Hu Ku rivers collection waters internal control hydrometric station data on flows and rainfall data day by day day by day are entered, the hydrology is calculated
Stand control section day by day flow it is corresponding collection waters before 1~first 10 days different duration T the total precipitation between Pierre
Inferior coefficient correlation, the size identification according to coefficient go out and the section significantly correlated crucial duration T of flow day by day;
Step 3:Using the crucial duration T identified in step 2, flow day by day according to the actual measurement of hydrometric station control section, correspondingly
Crucial duration T quantum of rainfall, upstream collection water surface area, utilize linear regression least square method, it is real to establish control hydrometric station
Total rainfall linear relation in measurement of discharge crucial duration T corresponding with collection waters;When removing flow-T with the catchment area of control of section
Total precipitation relational expression in length, tries to achieve basin unit area Rainfall-runoff coefficient;
Step 4:According to hydrometric station control section flow~rainfall of foundation, collection water surface area, runoff coefficient relation, calculate each
Enter the flow day by day in Hu Ku rivers, obtain rainfall and enter Hu Ku total Waters day by day;According to lake reservoir level, the lake reservoir level curve of areas
Water level diurnal variation caused by rainfall runoff enters lake storehouse is calculated, obtains lake reservoir level diurnal variation sequence 1 caused by rainfall;
Step 5:According to the tailwater quantity for the discharge of sewage treatment plant day by day being collected into, water factory's water withdrawal, lake storehouse water surface evaporation, lake
Storehouse water surface rainfall, lake storehouse water yield and lake reservoir level data, based on principle of water balance, calculate is caused by rainfall runoff
The Day-to-day variability amount of lake reservoir level, obtain lake reservoir level diurnal variation sequence 2 caused by rainfall;
Step 6:The comparison of lake reservoir level diurnal variation sequence 2 and the change curve of lake reservoir level diurnal variation sequence 1 is carried out, adjustment participates in
The runoff coefficient of each flow rate calculation day by day for mainly entering Hu Ku rivers so that the change curve of lake reservoir level diurnal variation sequence 1 approaches
Lake reservoir level diurnal variation sequence 2, each day runoff coefficient after being corrected;
Step 7:The flow in Hu Ku river courses is respectively entered according to day runoff coefficient calculating after amendment, then by sewage treatment plant tail water day by day
Discharge capacity is added to day by day on flow of corresponding river, obtains the flow that each river enters Hu Ku day by day.
2. according to the method for claim 1, it is characterised in that in the step 2, according to hydrometric station measured discharge, section
Upstream collection waters rainfall, calculate section daily average water discharge and three comprising the same day, not comprising the same day, not comprising the same day and proxima luce (prox. luc)
Under kind situation, time span is the linearly related Pierre of 1 to 10 total rainfall and hydrometric station measured discharge in gathering ground
Inferior coefficient correlation, choose maximum Pearson correlation coefficient correspond to time span T, as section flow rate calculation use rainfall when
It is long.
3. according to the method for claim 1, it is characterised in that in the step 2, for entering lake storehouse upper river collection waters
In the presence of the collection larger reservoir situation of water surface area, reservoir catchment is deducted from river collection water surface area when calculating river collection water surface area
Domain area.
4. according to the method for claim 1, it is characterised in that in the step 3, basin unit area Rainfall-runoff is formed
Control section flow rate expression, obtained according to formula (1) using least square method:
QActual measurement, t=α PActual measurement, t-1·SActual measurement (1)
In formula, QActual measurement, tFor measured profile daily average water discharge of the river in t days, α is basin unit area Rainfall-runoff coefficient, PActual measurement, t-1
For quantum of rainfall in the T durations of section upstream collection waters, SActual measurementFor river course measured profile upstream collection water surface area.
5. according to the method for claim 4, it is characterised in that in the step 4, respectively entering Hu Ku rivers, flow passes through day by day
Following formula calculates,
QI, t=α Pi,t-1·Si (2)
In formula, QI, tEnter lake storehouse flow preresearch estimates value in t days for river course i;Pi,t-1For river course i collection waters before t days T durations
Interior quantum of rainfall;SiFor river course i collection water surface areas.
6. according to the method for claim 5, it is characterised in that in the step 6, based on lake storehouse principle of water balance, profit
Lake level Day-to-day variability rate is calculated with formula (3) measured data:
Δlt·SC, t=∑ Qi,t+(Pt-Et)·SC, t-At+Tt (3)
In formula:ΔltSEA LEVEL VARIATION in the t days obtained is calculated by measured water level value for Hu Ku;SC, tFor Hu Ku within t days average water
Face area, Δ lt·SC, tFor t day inland lakes reservoir storage variable quantity;PtIt is Hu Ku in t daily rainfalls;EtTo be evaporated in t Hu storehouses
Amount;AtIt is Hu Ku water factories in t day water withdrawals;TtTailwater quantity, ∑ Q are discharged in t days for sewage treatment plant in Hu Ku basinsI, tTo enter lake
Storage enters water sum river t days.
7. according to the method for claim 6, it is characterised in that be adjusted according to formula (3) watershed runoff coefficient, such as
Shown in formula (4):
α=[Δ lt·SC, t-(Pt-Et)·SC, t+A-Tt]/∑(PI, t-1·Si) (4)
Using the unit area Rainfall-runoff coefficient after adjustment, calculate and obtain the flow that each river enters Hu Ku day by day.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710892965.XA CN107679021B (en) | 2017-09-27 | 2017-09-27 | River lake-entering reservoir flow calculation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710892965.XA CN107679021B (en) | 2017-09-27 | 2017-09-27 | River lake-entering reservoir flow calculation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107679021A true CN107679021A (en) | 2018-02-09 |
CN107679021B CN107679021B (en) | 2020-10-30 |
Family
ID=61137759
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710892965.XA Active CN107679021B (en) | 2017-09-27 | 2017-09-27 | River lake-entering reservoir flow calculation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107679021B (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108428057A (en) * | 2018-03-16 | 2018-08-21 | 中国科学院南京地理与湖泊研究所 | A kind of computational methods of lake exchange water cycle |
CN109738019A (en) * | 2019-01-31 | 2019-05-10 | 山东大学 | A kind of small-sized shallow lake hydro-ecological safety red line pre-alarming control system and method |
CN109992904A (en) * | 2019-04-08 | 2019-07-09 | 中国电建集团贵阳勘测设计研究院有限公司 | Water balance calculation method for depression-spring system of karst area |
CN110541389A (en) * | 2019-09-06 | 2019-12-06 | 中国电建集团贵阳勘测设计研究院有限公司 | flood control and flood fighting emergency system in construction period and early warning method thereof |
CN111400655A (en) * | 2020-06-08 | 2020-07-10 | 中国水利水电科学研究院 | Correction optimization method and system for warehousing traffic |
CN112982273A (en) * | 2021-02-18 | 2021-06-18 | 中国水利水电科学研究院 | Ecological space optimization and function improvement method for river wetland |
CN113240171A (en) * | 2021-05-10 | 2021-08-10 | 福建省水利水电勘测设计研究院 | Method for determining optimal water diversion flow of estuary leaf vein-shaped river network main flow |
CN114819322A (en) * | 2022-04-18 | 2022-07-29 | 武汉大学 | Method for forecasting lake inflow flow of lake |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10204966A (en) * | 1997-01-22 | 1998-08-04 | Ebara Corp | Method for estimating tide level and method for simulating river level |
CN104021283A (en) * | 2014-05-23 | 2014-09-03 | 清华大学 | Prediction method and device of day runoff volume of snowmelt period |
CN104732073A (en) * | 2015-03-04 | 2015-06-24 | 河海大学 | Computing method for surface water and underground water coupling simulation |
CN106874625A (en) * | 2017-03-17 | 2017-06-20 | 河北省水利水电勘测设计研究院 | The determination method and determination system of a kind of front damming volume |
-
2017
- 2017-09-27 CN CN201710892965.XA patent/CN107679021B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10204966A (en) * | 1997-01-22 | 1998-08-04 | Ebara Corp | Method for estimating tide level and method for simulating river level |
CN104021283A (en) * | 2014-05-23 | 2014-09-03 | 清华大学 | Prediction method and device of day runoff volume of snowmelt period |
CN104732073A (en) * | 2015-03-04 | 2015-06-24 | 河海大学 | Computing method for surface water and underground water coupling simulation |
CN106874625A (en) * | 2017-03-17 | 2017-06-20 | 河北省水利水电勘测设计研究院 | The determination method and determination system of a kind of front damming volume |
Non-Patent Citations (4)
Title |
---|
WEIPING HU: "A review of the models for Lake Taihu and their application in lake environmental management", 《ECOLOGICAL MODELLING》 * |
李娜 等: "基于统计模型LOADEST的宝象河污染物通量估算", 《应用基础与工程科学学报》 * |
汪鹭: "博斯腾湖流域主要入湖河流径流序列分析", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
环境保护部环境工程评估中心: "《环境影响评价技术方法2017年版》", 28 February 2017 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108428057A (en) * | 2018-03-16 | 2018-08-21 | 中国科学院南京地理与湖泊研究所 | A kind of computational methods of lake exchange water cycle |
CN109738019A (en) * | 2019-01-31 | 2019-05-10 | 山东大学 | A kind of small-sized shallow lake hydro-ecological safety red line pre-alarming control system and method |
CN109992904A (en) * | 2019-04-08 | 2019-07-09 | 中国电建集团贵阳勘测设计研究院有限公司 | Water balance calculation method for depression-spring system of karst area |
CN110541389A (en) * | 2019-09-06 | 2019-12-06 | 中国电建集团贵阳勘测设计研究院有限公司 | flood control and flood fighting emergency system in construction period and early warning method thereof |
CN111400655A (en) * | 2020-06-08 | 2020-07-10 | 中国水利水电科学研究院 | Correction optimization method and system for warehousing traffic |
CN111400655B (en) * | 2020-06-08 | 2020-09-18 | 中国水利水电科学研究院 | Correction optimization method and system for warehousing traffic |
CN112982273A (en) * | 2021-02-18 | 2021-06-18 | 中国水利水电科学研究院 | Ecological space optimization and function improvement method for river wetland |
CN113240171A (en) * | 2021-05-10 | 2021-08-10 | 福建省水利水电勘测设计研究院 | Method for determining optimal water diversion flow of estuary leaf vein-shaped river network main flow |
CN113240171B (en) * | 2021-05-10 | 2022-06-03 | 福建省水利水电勘测设计研究院有限公司 | Method for determining optimal water diversion flow of estuary leaf vein-shaped river network main flow |
CN114819322A (en) * | 2022-04-18 | 2022-07-29 | 武汉大学 | Method for forecasting lake inflow flow of lake |
CN114819322B (en) * | 2022-04-18 | 2024-05-14 | 武汉大学 | Forecasting method for flow of lake entering lake |
Also Published As
Publication number | Publication date |
---|---|
CN107679021B (en) | 2020-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107679021A (en) | A kind of river Ru Hu storehouses method of calculating flux | |
US20180058932A1 (en) | Method for analyzing the types of water sources based on natural geographical features | |
Beaumont | Agricultural and environmental changes in the upper Euphrates catchment of Turkey and Syria and their political and economic implications | |
CN110838079A (en) | Intelligent city flood control and drainage remote monitoring system and method | |
CN111310976B (en) | Quick water quality forecasting and early warning method and system based on distributed computing | |
Ibrakhimov | Spatial and temporal dynamics of groundwater table and salinity in Khorezm (Aral Sea Basin), Uzbekistan | |
CN112116785A (en) | Tailing pond disaster early warning method and device based on strong rainfall weather forecast | |
Bonacci et al. | The water on a small karst island: the island of Korčula (Croatia) as an example | |
Wang et al. | Identifying the driving factors of sediment delivery ratio on individual flood events in a long-term monitoring headwater basin | |
CN109002946B (en) | Water resource optimization scheduling method for 'two-reservoir-two-station' system for replenishing water in rivers and lakes | |
CN110852525A (en) | Reservoir flood control limit water level dynamic adjustment method based on forecast error override defense | |
MacMillan | Predicting the general yield class of Sitka spruce on better quality land in Scotland | |
Tabatabaee et al. | Rainwater harvesting potentials for drought mitigation in Iran | |
CN104281783A (en) | Water pollutant net accumulation based water pollutant enrichment judging method and experimental device | |
CN113129318B (en) | Method for calculating water storage capacity of stagnant flood area by utilizing SAR (synthetic aperture radar) image | |
Ansarifard et al. | Monitoring drought by reconnaissance drought index (RDI) and standardized precipitation index (SPI) using DrinC software | |
CN114386796A (en) | Sponge city project construction effect and overall quality evaluation method | |
CN110084848B (en) | Method for determining water collection area threshold value during digital river network extraction | |
CN110400049A (en) | The ratification method of the ecological flow of rural hydropower station and appraise and decide device | |
Andersson et al. | Consequences of changed wetness on riverine nitrogen–human impact on retention vs. natural climatic variability | |
CN109002422B (en) | Method for acquiring influence specific gravity of recycled water reutilization and spatial variation on net inflow water production rate | |
Orphanopoulos et al. | Water use efficiency in an arid watershed: a case study | |
CN106368098A (en) | Rapid design method of road bioretention trough | |
Esther | Using GIS techniques to determine RUSLE’s ‘R’and ‘LS’Factors for kapingazi river catchment | |
Klepper | Model study of the Negara river basin to assess the regulating role of its wetlands |
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 |