CN108416049A - A kind of high and cold mountain area basin sleet mixing Runoff calculation method - Google Patents
A kind of high and cold mountain area basin sleet mixing Runoff calculation method Download PDFInfo
- Publication number
- CN108416049A CN108416049A CN201810223041.5A CN201810223041A CN108416049A CN 108416049 A CN108416049 A CN 108416049A CN 201810223041 A CN201810223041 A CN 201810223041A CN 108416049 A CN108416049 A CN 108416049A
- Authority
- CN
- China
- Prior art keywords
- basin
- soil
- period
- calculated
- snow
- 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
- 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
Abstract
The present invention discloses a kind of high and cold mountain area basin sleet mixing Runoff calculation method, and step includes：Obtain the precipitation, temperature and evaporation data of period to be calculated；3 production fluid space inhomogeneities curves, i.e. Snow Cover Area attenuation curve, disengagement area distribution curve and tension water reservoir capacity curve are established using history meteorological model data；It establishes the storage determined according to soil freezing state completely super ooze and combines Runoff calculation mechanism, to precipitation, temperature and evaporation data based on the period to be calculated, with 3 production fluid space inhomogeneities curves, Runoff calculation mechanism is combined using full super ooze is stored, the basin sleet mixing runoff yield at moment to be calculated is calculated.The present invention can solve high and cold mountain area basin sleet mixing Runoff calculation problem, reduce in hydrological distribution model Runoff calculation to the degree of dependence of underlying surface data.
Description
Technical field
The present invention relates to River Basin Hydrology forecasting technique field, especially a kind of high and cold mountain area basin sleet mixing Runoff calculation
Method.
Background technology
Application distribution formula hydrological model is the important tool for carrying out River Basin Hydrology process simulation and forecast.High-cold regions hydrology gas
As and land surface condition there is significant vertical variability, season/day-night change of temperature ice and snow, freeze-thawing process of soil to be caused to be handed over
Cycle is knitted, extremely complex hydrologic process is formed.For this purpose, when building high and cold mountain area hydrological distribution model, need to establish one
The Runoff calculation method that kind can characterize freeze-thawing process of soil, adapt to high and cold mountain area data condition.
Currently, in cold highland area basin carry out Runoff calculation hydrological distribution model it is numerous, such as SWAT, VIC and
TOPKAPI etc..It would be appreciated that these models lacked to the considerations of runoff-generating model changes in soil freezing-thawing conversion, it is high and cold
The scarcity of mountain area underlying surface data also further limits model computational accuracy.
Explanation of nouns
Basin of Xin An Jiang hydrological model, basic principle are：Full basin is divided into many module unit basins, to each unit stream
Production runoff concentration calculation is made in domain, obtains the rate of discharge process in unit basin, then carries out exporting flood routing for river channel below, acquires stream
The discharge process of domain outlet；The outflow process in each unit basin is added, total outflow process in basin is obtained.
SWAT hydrological distribution models, basic principle are：Full basin is divided into several Hydrologic response units, to each water
Literary response unit makees production stream and slope concentration calculates, obtain water, sand, nutriment and chemical substance of Hydrologic response units etc. to
The input process of main stem, then carry out concentration of channel calculating acquires water, sand, nutriment and chemical substance etc. from the network of waterways to stream
The defeated of domain outlet moves past journey.
Invention content
The object of the present invention is to provide a kind of high and cold mountain area basin sleet mixing Runoff calculation methods, expired by establishing storage-
It is super to ooze combination runoff mechanism and 3 production fluid space inhomogeneities curves (Snow Cover Area attenuation curve, disengagement area distribution songs
Line, tension water reservoir capacity curve), high and cold mountain area basin sleet mixing Runoff calculation problem is solved, Distributed Hydrological mould is reduced
Degree of dependence of the Runoff calculation to underlying surface data in type.
The technical solution that the present invention takes is：A kind of high and cold mountain area basin sleet mixing Runoff calculation method, including：
S1 obtains the basin precipitation, temperature and evaporation data of period to be calculated, basin drop is calculated based on the data got
Rainfall and snowfall；
S2, system when establishing snow monitoring and snowpack water equivalent and accumulated snow all standing in basin between water equivalent three
Count relation function, i.e. Snow Cover Area attenuation curve；
S3 establishes the disengagement area distribution curve of the spatial non-uniformity of characterization basin evaporation；
S4 judges whether period soil to be calculated is in and freezes shape according to the relationship between the earth's surface soil moisture and temperature
State；
S5 is based on Snow Cover Area attenuation curve and disengagement area distribution curve, and it is mixed to calculate sleet according to the judging result of S4
Close runoff yield：
When soil is in frozen state, sleet mixing runoff yield X is calculated using runoff yield excess mechanism；
When soil is in non-frozen state, tension is characterized using tension water reservoir capacity curve in the hydrological model of the Xinanjiang River
Inhomogeneities of the water reservoir capacity on space scale calculates sleet mixing runoff yield Y according to runoff yield under saturated storage mechanism.
Curve in the present invention in Snow Cover Area attenuation curve, disengagement area distribution curve and tension water reservoir capacity curve
Relation Parameters between parameter and the soil moisture and temperature can be based on history meteorological model data, using parameter calibration method
It chooses and determines.History meteorological model data includes each curve variable such as history precipitation, temperature, runoff yield, evaporation recorded
Related data.
Preferably, in S1, basin precipitation and temperature record based on the period to be calculated carry out solid-liquid using temperature threshold method
State precipitation detaches, and then calculates rainfall and snowfall.I.e. be arranged a mushy stage distinguish temperature threshold value, by temperature threshold value with
On precipitation data be regarded as precipitation data, temperature threshold value precipitation data below is regarded as snowfall data, you can obtain basin it is each when
The rainfall and snowfall of section.Temperature threshold value may be configured as 0 degree Celsius.
If the basin precipitation that S1 is obtained, temperature are history measured datas with evaporation data, runoff yield is finally calculated
It is historical simulation value；If the basin precipitation that S1 is obtained, temperature are predictions for future data with evaporation data, production is finally calculated
Flow is predictions for future value.
Preferably, S2 is based on Snow Cover Area attenuation curve in SWAT hydrological distribution models, obtains snow cover in basin
Rate is：
In formula (1), C1And C2For profile shape parameter, ScovFor snow monitoring, S100%Basin 100% is covered for accumulated snow
When water equivalent, StFor the basin snowpack water equivalent of period t, and：
St=Ps,t+St-1′ (2)
In formula (2), Ps,tFor the snowfall of t-th of period, St-1' be the t-1 period remaining snowpack water equivalent.Product
Snow face accumulates the profile shape parameter C of attenuation curve1And C2, history meteorological model data are based on, are chosen using parameter calibration method true
It is fixed.
Preferably, the disengagement area distribution curve of S3 foundation is：
In formula (3), F (Ei) be basin in any point i evaporability EiProbability distribution on space scale, EmmFor sky
Interior single-point maximum evaporation ability, Emm=k1·Ew, EwFor the evaporation data of the S1 respective points obtained, k1For evaporability sky
Between conversion coefficient；
The basin evaporability that is averaged is：λ is the index for reflecting evaporating space inhomogeneities,
Then basin actual evaporation is：Ea=k2·Em, wherein k2For evaporator reduction coefficient.
Above-mentioned parameter λ, k1And k2History meteorological model data can be utilized respectively, are chosen and are determined using parameter calibration method.Ew
It is that the historical values (the Pan evaporation amount of actual measurement) that can be surveyed or the evaporation capacity data of future anticipation (are calculated
Potential evaporation ability).
Preferably, in S4, judge soil whether in frozen state include step：
S41 establishes surface soil temperature TsWith temperature TaBetween linear statistical relationship be：
Ts=aTa+b (4)
In formula (4), a and b are coefficient；
The soil moisture value T of period to be calculated is calculated using formula (4) by S42s,i；
S43 defines soil freezing threshold value Ts,base, moment soil moisture value T to be calculated that S42 is calculateds,iWith soil
Earth freeze threshold Ts,baseIt is compared：If Ts,i≤Ts,base, then soil be in frozen state；If Ts,i＞ Ts,base, then at soil
In non-frozen state.
Above-mentioned surface soil temperature TsWith temperature TaBetween linear statistical relationship in coefficient a and b, utilize history meteorological
Surface soil temperature data in hydrographic data and temperature record are chosen using parameter calibration method and are determined.
Preferably, in S5, when soil is in frozen state, sleet mixing runoff yield X is calculated using runoff yield excess mechanism
For：
X=I-Ea (5)
In formula (5), I is period water purification amount, and the water purification amount I of period tt=Pr,t+Mt, Pr,tFor the rainfall of period t, Mt
For the snow melt water of period t.
Preferably, the snow melt water M of the period ttCalculating include step：
S511, using the potential snow melt ability M of degree-day factor method calculation interval tt,max：
Mt,max=DDFScov·(Ta,t-Tmelt) (6)
In formula (6), DDF is degree-day factor, TmeltFor snow melting temperature, Ta,tFor the temperature of period t；
S512, the potential snow melt ability M that will be acquiredt,maxWith snowpack water equivalent StIt is compared：If Mt,max≤St, then melt
Snow-broth amount Mt=Mt,max；If Mt,max＞ St, then snow melt water Mt=St；
Then the remaining snowpack water equivalent of period t is：S′t=St-Mt。
The determination of above-mentioned degree-day factor DDF is also chosen using history meteorological model data using parameter calibration method true
It is fixed.
Preferably, in S5, when soil is in non-frozen state, the calculating step of sleet mixing runoff yield X is：
S521 obtains day part soil water-containing in curve based on tension water reservoir capacity curve in the hydrological model of the Xinanjiang River
Measure the corresponding curve ordinate reservoir capacity value A of W：
In formula (7), WMM=WM (1+B) is basin single-point maximum water-storage, and WM is that basin is averaged reservoir capacity,；
S522, calculating sleet mixing runoff yield Y according to runoff yield under saturated storage mechanism is：
Reservoir capacity exponent B in above-mentioned Xinanjiang River hydrological model in tension water reservoir capacity curve, it is available to go through
Basin reservoir storage related data in history meteorological model data is chosen using parameter calibration method and is determined.
Advantageous effect
By establishing, storage is full-super to ooze combination Runoff calculation mechanism and 3 production fluid space inhomogeneities curves to the present invention
(Snow Cover Area attenuation curve, disengagement area distribution curve, tension water reservoir capacity curve), and thus calculate on the grid of basin
Runoff yield under saturated storage amount or runoff yield excess amount can solve high and cold mountain area basin sleet mixing Runoff calculation problem, reduce Distributed Hydrological
Runoff calculation has stronger engineering significance to the degree of dependence of underlying surface data in model.
Description of the drawings
Fig. 1 show a kind of method flow schematic diagram of specific embodiment of the present invention.
Specific implementation mode
It is further described below in conjunction with the drawings and specific embodiments.
With reference to figure 1, high and cold mountain area basin sleet mixing Runoff calculation method of the invention, including：
S1 obtains the basin precipitation, temperature and evaporation data of period to be calculated, basin drop is calculated based on the data got
Rainfall and snowfall；
S2, system when establishing snow monitoring and snowpack water equivalent and accumulated snow all standing in basin between water equivalent three
Count relation function, i.e. Snow Cover Area attenuation curve；
S3 establishes the disengagement area distribution curve of the spatial non-uniformity of characterization basin evaporation；
S4 judges whether period soil to be calculated is in and freezes shape according to the relationship between the earth's surface soil moisture and temperature
State；
S5 is based on Snow Cover Area attenuation curve and disengagement area distribution curve, and it is mixed to calculate sleet according to the judging result of S4
Close runoff yield：
When soil is in frozen state, sleet mixing runoff yield X is calculated using runoff yield excess mechanism；
When soil is in non-frozen state, tension is characterized using tension water reservoir capacity curve in the hydrological model of the Xinanjiang River
Inhomogeneities of the water reservoir capacity on space scale calculates sleet mixing runoff yield Y according to runoff yield under saturated storage mechanism.
In step S1, basin precipitation and temperature record based on the period to be calculated carry out mushy stage using temperature threshold method
Precipitation detaches, and then calculates rainfall and snowfall.The temperature threshold value that one mushy stage is distinguished is set, more than temperature threshold value
Precipitation data be regarded as precipitation data, temperature threshold value precipitation data below is regarded as snowfall data, you can obtains basin day part
Rainfall and snowfall.Temperature threshold value may be configured as 0 degree Celsius.
If the basin precipitation that S1 is obtained, temperature are history measured datas with evaporation data, runoff yield is finally calculated
It is historical simulation value；If the basin precipitation that S1 is obtained, temperature are predictions for future data with evaporation data, production is finally calculated
Flow is predictions for future value.Calculating for history or forecast sleet mixing runoff yield, the present invention have universality.
Step S2 is based on Snow Cover Area attenuation curve in SWAT hydrological distribution models, obtains snow monitoring in basin
For：
In formula (1), C1And C2For profile shape parameter, ScovFor snow monitoring, S100%Basin 100% is covered for accumulated snow
When water equivalent, StFor the basin snowpack water equivalent of period t, and：
St=Ps,t+St-1′ (2)
In formula (2), Ps,tFor the snowfall of t-th of period, St-1' be the t-1 period remaining snowpack water equivalent.Product
Snow face accumulates the profile shape parameter C of attenuation curve1And C2, history meteorological model data are based on, are chosen using parameter calibration method true
It is fixed.
The fixed Snow Cover Area attenuation curve of parameter is utilized, can be calculated according to the snowpack water equivalent of t-th of period
Obtain the basin snow monitoring of t-th of period.
The disengagement area distribution curve that step S3 is established is：
In formula (3), F (Ei) be basin in any point i evaporability EiProbability distribution on space scale, EmmFor sky
Interior single-point maximum evaporation ability, Emm=k1·Ew, EwFor the evaporation data of the S1 respective points obtained, k1For evaporability sky
Between conversion coefficient；
Then the basin evaporability that is averaged is：λ is the index for reflecting evaporating space inhomogeneities,
Basin actual evaporation is：Ea=k2·Em, wherein k2For evaporator reduction coefficient.
In step S4, judge soil whether in frozen state include step：
S41 establishes surface soil temperature TsWith temperature TaBetween linear statistical relationship be：
Ts=aTa+b (4)
In formula (4), a and b are coefficient；
The soil moisture value T of period to be calculated is calculated using formula (4) by S42s,i；
S43 defines soil freezing threshold value Ts,base, moment soil moisture value T to be calculated that S42 is calculateds,iWith soil
Earth freeze threshold Ts,baseIt is compared：If Ts,i≤Ts,base, then soil be in frozen state；If Ts,i＞ Ts,base, then at soil
In non-frozen state.
Above-mentioned surface soil temperature TsWith temperature TaBetween linear statistical relationship in coefficient a and b, utilize history meteorological
Surface soil temperature data in hydrographic data and temperature record are chosen using parameter calibration method and are determined.
In step S5, when soil is in frozen state, calculating sleet mixing runoff yield X using runoff yield excess mechanism is：
X=I-Ea (5)
In formula (5), I is period water purification amount, and the water purification amount I of period tt=Pr,t+Mt, Pr,tFor the rainfall of period t, Mt
For the snow melt water of period t.
The snow melt water M of the period ttCalculating include step：
S511, using the potential snow melt ability M of degree-day factor method calculation interval tt,max：
Mt,max=DDFScov·(Ta,t-Tmelt) (6)
In formula (6), DDF is degree-day factor, TmeltFor snow melting temperature, Ta,tFor the temperature of period t；Basin accumulated snow covers
Lid rate ScovIt is calculated by Snow Cover Area attenuation curve.
S512, the potential snow melt ability M that will be acquiredt,maxWith snowpack water equivalent StIt is compared：If Mt,max≤St, then melt
Snow-broth amount Mt=Mt,max；If Mt,max＞ St, then snow melt water Mt=St；
Then the remaining snowpack water equivalent of period t is：S′t=St-Mt.And then can calculation interval t+1 snowpack water equivalent
St+1=Ps,t+1+S′tAnd the snow monitoring S of period t+1cov, and further obtain the snow melt water M of period t+1t+1, then
The present invention can be achieved to calculate the snow melt watermeter of all periods to be calculated, in terms of the sleet mixing runoff yield for supporting each corresponding period
It calculates.
The determination of above-mentioned degree-day factor DDF is also chosen using history meteorological model data using parameter calibration method true
It is fixed.
In step S5, when soil is in non-frozen state, the calculating step of sleet mixing runoff yield X is：
S521 obtains day part soil water-containing in curve based on tension water reservoir capacity curve in the hydrological model of the Xinanjiang River
Measure the corresponding curve ordinate reservoir capacity value A of W：
In formula (7), WMM=WM (1+B) is basin single-point maximum water-storage, and WM is that basin is averaged reservoir capacity；
S522, calculating sleet mixing runoff yield Y according to runoff yield under saturated storage mechanism is：
Reservoir capacity exponent B in above-mentioned Xinanjiang River hydrological model in tension water reservoir capacity curve, it is available to go through
Basin reservoir storage related data in history meteorological model data is chosen using parameter calibration method and is determined.
Embodiment
Certain existing high and cold mountain area small watershed 1 year by hour precipitation, temperature and evaporation data, method, grinds according to the present invention
Study carefully the basin in the phase is by the calculating process of hour runoff yield：
(1) definition is P by hour precipitation series1,P2,...,Pt..., temperature series isEvaporation
Train of energy isIt is respectively P to obtain rainfall and snowfall series with 0 DEG C for threshold calculationsr,1,
Pr,2,...,Pr,t... and Ps,1,Ps,2,...,Ps,t,...；
(2) moment is calculated in t, the snowfall P for this period of adding ups,tWork as with the remaining water of the moment basins t-1 snowpack
Measure St-1', obtain the water equivalent S of this period basin snowpackt=Ps,t+St-1′；
(3) snow surface in SWAT hydrological distribution models is determined using parameter calibration method based on history meteorological model data
Product attenuation curve, that is, determine snow monitoring S in basincovWith water equivalent of snow cover St, accumulated snow all standing when water equivalent S100%Letter
Number relationshipWherein C1, C2 are profile shape parameter；
Using fixed Snow Cover Area attenuation curve, the snow monitoring S of day part is calculatedcov；
(4) use degree-day factor method that the potential snow melt ability M of t moment is calculatedt,max=DDFScov(Ta,t-
Tmelt), wherein DDF is degree-day factor, Ta,tFor temperature, TmeltFor snow melting temperature；
The potential snow melt ability M that t moment is acquiredt,maxWith snowpack water equivalent StIt is compared：If Mt,max≤St, then
Snow melt water Mt=Mt,max；If Mt,max＞ St, then Mt=St；Then the water equivalent of the remaining snowpack of t moment is S 't=St-Mt；
(5) add up rainfall Pr,1,Pr,2,...,Pr,t... with snow melt water M1,M2,...,Mt..., it obtains for steaming
The water purification amount I of hair and Runoff calculation1,I2,...,It,...：It=Pr,t+Mt；
(6) it is based on history meteorological model data, establishes disengagement area distribution curve, and true using parameter calibration method
Determine parameter of curve, that is, establishes the evaporability E of any point i in basiniProbability-distribution function on space scale：F(Ei) be basin in any point i evaporability EiProbability distribution on space scale, Emm
For single-point maximum evaporation ability, E in spacemm=k1·Ew, EwFor the S1 respective points obtained evaporation data (i.e.), k1For evaporability space conversion coefficient；
The basin evaporability that is averaged is：λ is the index for reflecting evaporating space inhomogeneities；
Then basin actual evaporation is：Ea=k2·Em, wherein k2For evaporator reduction coefficient；
(7) it is based on history meteorological model data, surface soil temperature T is established using parameter calibration methodsWith temperature Ta
Linear statistical relationship Ts=aTa+ b, wherein a, b are coefficient, and soil moisture series is calculated
(8) judge the soil freezing state of each period to be calculated：If soil moisture TsLess than or equal to soil freezing threshold value
Ts,base, soil is in frozen state, if TsMore than Ts,base, then soil be in non-frozen state；
(9) according to soil freezing state computation sleet mixing runoff yield：
When soil is in frozen state, sleet mixing runoff yield X=I-E is calculated using runoff yield excess mechanisma, i.e. Xt=It-
Em,t；
When soil is in non-frozen state, sleet mixing runoff yield Y is calculated according to runoff yield under saturated storage mechanism：
Tension water reservoir capacity is characterized in space ruler first with tension water reservoir capacity curve in the hydrological model of the Xinanjiang River
Inhomogeneities on degree obtains reservoir capacity curve ordinate value corresponding with moment soil moisture content W is calculatedWherein WMM=WM (1+B) is basin single-point maximum water-storage, and WM is average for basin
Reservoir capacity, B are reservoir capacity exponent；
Then sleet mixing runoff yield Y is：
Using above-mentioned steps (2)~(9), the sleet mixing runoff yield of all periods can be calculated successively.
Parameter includes in the present invention：Snow Cover Area attenuation curve form parameter C1, C2, evaporator reduction coefficient k, evaporating space
Heterogeneity index λ, soil moisture parameter a, b, degree-day factor DDF, snow melting temperature Tmelt, basin is averaged reservoir capacity
WM and reservoir capacity exponent B can be based on history meteorological model data (precipitation, temperature, evaporation, runoff yield number
According to), it is chosen and is determined using parameter calibration method.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, without departing from the technical principles of the invention, several improvement and deformations can also be made, these improvement and deformations
Also it should be regarded as protection scope of the present invention.
Claims (8)
1. a kind of high and cold mountain area basin sleet mixing Runoff calculation method, characterized in that including：
S1 obtains the basin precipitation, temperature and evaporation data of period to be calculated, Basin Rainfall amount is calculated based on the data got
And snowfall；
S2, statistics when establishing snow monitoring and snowpack water equivalent and accumulated snow all standing in basin between water equivalent three are closed
Be function, i.e. Snow Cover Area attenuation curve；
S3 establishes the disengagement area distribution curve of the spatial non-uniformity of characterization basin evaporation；
S4 judges whether period soil to be calculated is in frozen state according to the relationship between the earth's surface soil moisture and temperature；
S5 is based on Snow Cover Area attenuation curve and disengagement area distribution curve, and sleet mixing production is calculated according to the judging result of S4
Flow：
When soil is in frozen state, sleet mixing runoff yield X is calculated using runoff yield excess mechanism；
When soil is in non-frozen state, stored using tension water reservoir capacity curve characterization tension water in the hydrological model of the Xinanjiang River
Inhomogeneities of the water capacity on space scale calculates sleet mixing runoff yield Y according to runoff yield under saturated storage mechanism.
2. according to the method described in claim 1, it is characterized in that, in S1, based on basin precipitation and temperature historical data, using gas
Warm threshold method carries out mushy stage precipitation separation, and then calculates rainfall and snowfall.
3. according to the method described in claim 1, it is characterized in that, S2 be based on SWAT hydrological distribution models in Snow Cover Area decay
Curve, obtaining snow monitoring in basin is：
In formula (1), C1And C2For profile shape parameter, ScovFor snow monitoring, S100%When covering basin 100% for accumulated snow
Water equivalent, StFor the basin snowpack water equivalent of period t, and：
St=Ps,t+St-1′ (2)
In formula (2), Ps,tFor the snowfall of t-th of period, St-1' be the t-1 period remaining snowpack water equivalent.
4. according to the method described in claim 1, it is characterized in that, S3 establish disengagement area distribution curve be：
In formula (3), F (Ei) be basin in any point i evaporability EiProbability distribution on space scale, EmmFor in space
Single-point maximum evaporation ability, Emm=k1·Ew, EwFor the evaporation data of the S1 respective points obtained, k1Turn for evaporability space
Change coefficient；
The basin evaporability that is averaged is：λ is the index for reflecting evaporating space inhomogeneities；
Then basin actual evaporation is：Ea=k2·Em, wherein k2For evaporator reduction coefficient.
5. according to the method described in claim 1, it is characterized in that, in S4, judge soil whether in frozen state include step：
S401 establishes surface soil temperature TsWith temperature TaBetween linear statistical relationship be：
Ts=aTa+b (4)
In formula (4), a and b are coefficient；
The soil moisture value T of period to be calculated is calculated using formula (4) by S402s,i；
S403 defines soil freezing threshold value Ts,base, moment soil moisture value T to be calculated that S402 is calculateds,iWith soil
Freeze threshold Ts,baseIt is compared：If Ts,i≤Ts,base, then soil be in frozen state；If Ts,i＞ Ts,base, then soil be in
Non- frozen state.
6. according to the method described in claim 1, it is characterized in that, in S4, when soil is in frozen state, utilize runoff yield excess
Mechanism calculates sleet mixing runoff yield X：
X=I-Ea (5)
In formula (5), I is period water purification amount, and the water purification amount I of period tt=Pr,t+Mt, Pr,tFor the rainfall of period t, MtFor when
The snow melt water of section t.
7. according to the method described in claim 6, it is characterized in that, the snow melt water M of the period ttCalculating include step：
S511, using the potential snow melt ability M of degree-day factor method calculation interval tt,max：
Mt,max=DDFScov·(Ta,t-Tmelt) (6)
In formula (6), DDF is degree-day factor, TmeltFor snow melting temperature, Ta,tFor the temperature of period t；
S512, the potential snow melt ability M that will be acquiredt,maxWith snowpack water equivalent StIt is compared：If Mt,max≤St, then snowmelt
Measure Mt=Mt,max；If Mt,max＞ St, then snow melt water Mt=St；
Then the remaining snowpack water equivalent of period t is：St'=St-Mt。
8. according to the method described in claim 1, it is characterized in that, in S5, when soil is in non-frozen state, sleet mixing production
The calculating step of flow X is：
S521 obtains W pairs of day part soil moisture content in curve based on tension water reservoir capacity curve in the hydrological model of the Xinanjiang River
The curve ordinate reservoir capacity value A answered：
In formula (7), WMM=WM (1+B) is basin single-point maximum water-storage, and WM is that basin is averaged reservoir capacity,；
S522, calculating sleet mixing runoff yield Y according to runoff yield under saturated storage mechanism is：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810223041.5A CN108416049B (en) | 2018-03-19 | 2018-03-19 | Method for calculating rain and snow mixed runoff yield of drainage basin in alpine mountain area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810223041.5A CN108416049B (en) | 2018-03-19 | 2018-03-19 | Method for calculating rain and snow mixed runoff yield of drainage basin in alpine mountain area |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108416049A true CN108416049A (en) | 2018-08-17 |
| CN108416049B CN108416049B (en) | 2020-07-17 |
Family
ID=63132055
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810223041.5A Active CN108416049B (en) | 2018-03-19 | 2018-03-19 | Method for calculating rain and snow mixed runoff yield of drainage basin in alpine mountain area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108416049B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109685236A (en) * | 2019-02-11 | 2019-04-26 | 国家电网有限公司 | A kind of northern reservoir spring runoff water source analysis method |
| CN109816154A (en) * | 2018-12-31 | 2019-05-28 | 浙江同川工程咨询有限公司 | A kind of small watershed mountain torrents forecasting procedure |
| CN110619109A (en) * | 2019-08-01 | 2019-12-27 | 中国林业科学研究院 | Precipitation phase state identification method in alpine region |
| CN112651189A (en) * | 2020-12-11 | 2021-04-13 | 中国电建集团昆明勘测设计研究院有限公司 | General basin water circulation simulation calculation method based on natural sub-basins |
| WO2021129680A1 (en) * | 2019-12-26 | 2021-07-01 | 清华大学 | Snowmelt flood prediction method and apparatus, electronic device, and storage medium |
| CN117556223A (en) * | 2024-01-12 | 2024-02-13 | 国能大渡河流域水电开发有限公司 | Multi-factor similarity-based snow melt runoff forecasting method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102034001A (en) * | 2010-12-16 | 2011-04-27 | 南京大学 | Design method for distributed hydrological model by using grid as analog unit |
| KR20160110035A (en) * | 2015-09-25 | 2016-09-21 | 유찬호 | System for predict accident |
| CN106033480A (en) * | 2015-03-13 | 2016-10-19 | 西安山脉科技发展有限公司 | Method for forecasting mountain torrents based on improved Xin'anjiang model |
| CN106951612A (en) * | 2017-03-06 | 2017-07-14 | 河海大学 | Dynamic water storage capacity Runoff calculation method in freeze-thawing process of soil |
-
2018
- 2018-03-19 CN CN201810223041.5A patent/CN108416049B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102034001A (en) * | 2010-12-16 | 2011-04-27 | 南京大学 | Design method for distributed hydrological model by using grid as analog unit |
| CN102419788A (en) * | 2010-12-16 | 2012-04-18 | 南京大学 | Method for designing distributed-type hydrographical model based on penetration-storage integrated dynamic runoff yield mechanism |
| CN106033480A (en) * | 2015-03-13 | 2016-10-19 | 西安山脉科技发展有限公司 | Method for forecasting mountain torrents based on improved Xin'anjiang model |
| KR20160110035A (en) * | 2015-09-25 | 2016-09-21 | 유찬호 | System for predict accident |
| CN106951612A (en) * | 2017-03-06 | 2017-07-14 | 河海大学 | Dynamic water storage capacity Runoff calculation method in freeze-thawing process of soil |
Non-Patent Citations (5)
| Title |
|---|
| JUN LI ET AL.: ""Two universal runoff yield models: SCS vs. LCM"", 《JOURNAL OF GEOGRAPHICAL SCIENCES》 * |
| 崔泰昌 等: ""试论蓄满产流模型与超渗产流模型"", 《山西水利科技》 * |
| 梁忠民 等: ""基于统计理论的产流模型"", 《水科学发展》 * |
| 梁忠民 等: ""考虑空间变异性的统计产流模型研究"", 《南京大学学报（自然科学）》 * |
| 芮孝芳 等: ""流域产流分析及计算"", 《水利发电学报》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109816154A (en) * | 2018-12-31 | 2019-05-28 | 浙江同川工程咨询有限公司 | A kind of small watershed mountain torrents forecasting procedure |
| CN109685236A (en) * | 2019-02-11 | 2019-04-26 | 国家电网有限公司 | A kind of northern reservoir spring runoff water source analysis method |
| CN110619109A (en) * | 2019-08-01 | 2019-12-27 | 中国林业科学研究院 | Precipitation phase state identification method in alpine region |
| WO2021129680A1 (en) * | 2019-12-26 | 2021-07-01 | 清华大学 | Snowmelt flood prediction method and apparatus, electronic device, and storage medium |
| CN112651189A (en) * | 2020-12-11 | 2021-04-13 | 中国电建集团昆明勘测设计研究院有限公司 | General basin water circulation simulation calculation method based on natural sub-basins |
| CN117556223A (en) * | 2024-01-12 | 2024-02-13 | 国能大渡河流域水电开发有限公司 | Multi-factor similarity-based snow melt runoff forecasting method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108416049B (en) | 2020-07-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108416049A (en) | A kind of high and cold mountain area basin sleet mixing Runoff calculation method | |
| CN106951612B (en) | Dynamic water storage capacity runoff yield calculation method in soil freezing and thawing process | |
| Muthanna et al. | Seasonal climatic effects on the hydrology of a rain garden | |
| Kuchment et al. | A distributed model of runoff generation in the permafrost regions | |
| Kane et al. | An extreme rainfall/runoff event in Arctic Alaska | |
| CN115391712A (en) | Urban flood risk prediction method | |
| Cary et al. | Is the soil frozen or not? An algorithm using weather records | |
| Sato et al. | Analysis of long‐term water balance in the source area of the Yellow River basin | |
| Hrissanthou et al. | Estimate of sediment inflow into Vistonis Lake, Greece | |
| CN110147620A (en) | A kind of cryosphere hydrologic process analysis method and system | |
| Lehning et al. | Operational use of a snowpack model for the avalanche warning service in Switzerland: model development and first experiences | |
| YOSHIKAWA | Water balance dynamics of three small catchments in a Sub-Arctic boreal forest | |
| Janowicz et al. | Wolf Creek Research basin water balance studies | |
| Bahremand et al. | Application of WetSpa model for assessing land use impacts on floods in the Margecany–Hornad watershed, Slovakia | |
| Kido et al. | Subglacial drainage system changes of the Gulkana Glacier, Alaska: discharge and sediment load observations and modelling | |
| Mudashiru et al. | Evaluation of Methods of Peak Runoff Determination using Catchment Characteristics for Jere Sub-basin, Gurara River Basin, North Central Nigeria | |
| Hejduk et al. | Recorded lag times of snowmelt events in a small catchment | |
| Rychak et al. | Mathematical modeling of drinking water availability in Kharkiv region (Ukraine) at different dynamics of global climate warming | |
| Gelfan | Prediction of runoff in poorly gauged basins using a physically based model | |
| Gurevich | Freezing types and winter runoff of small rivers | |
| Srinivasan et al. | Regional validity of the parameters of a distributed runoff-erosion model in the semi-arid region of Brazil | |
| Zhou et al. | Research on the Forecasting and Risk Analysis Method of Snowmelt Flood. | |
| Qiao et al. | Key algorithms and its realization about snowmelt flood disaster model based on remote sensing and gIS | |
| Medina et al. | Stormwater: Urban Modeling | |
| Vieux et al. | Continuous Distributed Modeling of LID/GI: Scaling from Site to Watershed |
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 |