CN109558617A - Acquisition methods, device, computer equipment and the readable storage medium storing program for executing of diameter flow data - Google Patents

Abstract

The present invention provides acquisition methods, device, computer equipment and the readable storage medium storing program for executing of a kind of diameter flow data.This method comprises: determining accumulated snow melt water amount, glacier melting amount and the snow cover maximum water-holding capacity at current time according to the driving data in region to be measured and preset hydrological simulation model；By the difference of the gross precipitation in region to be measured and the solid precipitation quantity in region to be measured, the liquid precipitation amount that is determined as in preset time period；According to the freezing water quantity again of the snow cover of the water-holding capacity of previous moment, previous moment and the water holding variable quantity of previous moment, the water-holding capacity at current time is obtained；According to accumulated snow melt water amount, glacier melting amount, liquid precipitation amount, snow cover maximum water-holding capacity and the water-holding capacity at current time, the water holding variable quantity at current time is obtained；By the difference between ideal meltwater runoff and the water holding variable quantity at current time, it is determined as the practical diameter flow data at current time.The above method improves the counting accuracy of the diameter flow data in region to be measured.

Description

Acquisition methods, device, computer equipment and the readable storage medium storing program for executing of diameter flow data

Technical field

The present invention relates to Water Resources Domain, more particularly to a kind of acquisition methods of diameter flow data, device, computer equipment and Readable storage medium storing program for executing.

Background technique

With the development of science and technology, people are more and more deep for the research of terrestrial climate.And cryosphere is as global climate The important component of system has important influence to climate change, therefore, just seems for the hydrologic research of cryosphere It is even more important.Diameter flow data has very important work to climate change as an important research object of hydrologic process With.

In traditional technology, the measuring and calculating for the diameter flow data of cryosphere is to be predicted to obtain by traditional hydrological model 's.

But the diameter flow data predicted by traditional hydrological model cannot accurately react frost collar region Diameter flow data, cause cryosphere diameter flow data calculate accuracy it is lower.

Summary of the invention

Based on this, the present invention provide the acquisition methods of diameter flow data for cryosphere a kind of, device, computer equipment and Readable storage medium storing program for executing, for accurately calculating the diameter flow data of cryosphere, to improve the precision of diameter flow data.

In a first aspect, the embodiment of the present invention provides a kind of acquisition methods of diameter flow data, comprising:

According to the driving data in region to be measured and preset hydrological simulation model, determine current time accumulated snow melt water amount, Glacier melting amount and snow cover maximum water-holding capacity；The driving data in the region to be measured includes: the area to be measured in preset time period The air themperature of the gross precipitation in domain, the surface temperature in region to be measured and region to be measured；In the preset hydrological simulation model Parameter comprising the melt water characteristic for characterizing cryosphere；

By the difference of the gross precipitation in the region to be measured and the solid precipitation quantity in the region to be measured, it is determined as described Liquid precipitation amount in preset time period；

Become according to the freezing water quantity again of the snow cover of the water-holding capacity of previous moment, previous moment and the water holding of previous moment Change amount obtains the water-holding capacity at current time；

According to the accumulated snow melt water amount, the glacier melting amount, the liquid precipitation amount, the snow cover maximum water-holding capacity And the water-holding capacity at the current time, obtain the water holding variable quantity at current time；

By the difference between ideal meltwater runoff and the water holding variable quantity at the current time, it is determined as the reality at current time Border diameter flow data.

The driving data according to region to be measured and preset hydrological simulation model in one of the embodiments, really Determine accumulated snow melt water amount, glacier melting amount and the snow cover maximum water-holding capacity at current time, comprising:

According to the snow melt parameter in the driving data in the region to be measured and preset hydrological simulation model, when determining current The possibility amount of snowmelt amount of the water equivalent of snow at quarter and the region to be measured；The snow melt parameter includes melting for characterizing precipitation or accumulated snow The environmental parameter of change and parameter for characterizing the amount of snowmelt amount under unit time unit temperature；

Minimum value between the possibility amount of snowmelt amount and the water equivalent of snow in the region to be measured, is determined as the accumulated snow melt water Amount；

According to the possibility amount of snowmelt amount in the region to be measured, the accumulated snow melt water amount in the region to be measured, the water equivalent of snow with And the ice-melt parameter in the hydrological simulation model, determine the glacier melting amount；The ice-melt parameter includes for characterizing ice The environmental parameter that river is melted and the parameter for characterizing the ice storage under unit time unit temperature；

By the product of the water equivalent of snow and the preset water holding factor, it is determined as the snow cover maximum water-holding capacity；Wherein, The liquid water that the water holding factor is held by unit water equivalent of snow.

Being used to characterize precipitation or the environmental parameter of snow melting described in one of the embodiments, includes: snow melt threshold It is worth temperature, snowfall threshold temperature and rainfall threshold temperature；

The parameter for characterizing the amount of snowmelt amount under unit time unit temperature includes: that the maximum snow melt of summer solstice is subsisted The minimum snow melt degree-day factor of the factor and winter solstice；

It is described for characterize the environmental parameter of Melting Glacierss to include: ice-melt threshold temperature；

The parameter for characterizing the ice storage under unit time unit temperature includes: ice-melt degree-day factor.

It is described according to the accumulated snow melt water amount, the glacier melting amount, the liquid precipitation in one of the embodiments, Amount, the snow cover maximum water-holding capacity and the water-holding capacity at the current time, obtain the water holding variable quantity at current time, wrap It includes:

Determine the accumulated snow melt water amount, the glacier melting amount and the liquid precipitation amount and value；

Determine the difference between the snow cover maximum water-holding capacity and the water-holding capacity at the current time；

Determine that the minimum value between value and the difference is the water holding variable quantity at the current time.

The water-holding capacity according to previous moment, the snow cover of previous moment freeze again in one of the embodiments, The water holding variable quantity of water and previous moment, obtains the water-holding capacity at current time, comprising:

According to formulaDetermine the possibility of previous moment freezing water quantity again W_t-1,ref；Wherein, the RFC is preset water holding freeze-out fraction again, the T_i,mltFor the ice-melt in the hydrological simulation model Threshold temperature, the d_i,mltFor the ice-melt degree-day factor in the hydrological simulation model, the Δ t is to calculate duration；

By the possibility of previous moment freezing water quantity W again_t-1,refMinimum between the water-holding capacity of the previous moment Value, is determined as the freezing water quantity again of the snow cover of previous moment；

By between the water-holding capacity of the previous moment and the water holding variable quantity of previous moment and value, previous moment is subtracted Obtained difference after the freezing water quantity again of snow cover, is determined as the water-holding capacity at the current time.

The driving data according to region to be measured and preset hydrological simulation model in one of the embodiments, really Before accumulated snow melt water amount, glacier melting amount and the snow cover maximum water-holding capacity of determining current time, the method also includes:

According to known model training data and preset CREST model, the actual value and ice-melt ginseng of snow melt parameter are determined Several actual values；Wherein, the initial setting of the initial set value and ice-melt parameter in the CREST model including snow melt parameter Value；

The initial set value of snow melt parameter in the CREST model is replaced with to the actual value of the snow melt parameter, with And the initial set value of the ice-melt parameter in the CREST model is replaced with to the actual value of the ice-melt parameter, it obtains described Hydrological simulation model；

Wherein, the model training data include: the gross precipitation in trained region, the trained region surface temperature, The air themperature in the trained region, the practical water equivalent of snow of actual measurement, actual measurement practical glacier melting amount.

The model training data according to known to and preset CREST model in one of the embodiments, determination are melted Avenge the actual value of parameter, comprising:

According to formulaDetermine the possibility amount of snowmelt amount in the trained region S_t-1,mlt'；

Wherein, snow melt degree-day factor Doy is number of days, T_s,mlt' it is snow melt threshold temperature, d_s,mlt6For the maximum snow melt degree-day factor of summer solstice, d_s,mlt12For the winter solstice Minimum snow melt degree-day factor, T_ls' be the trained region surface temperature, T_airmx' be the trained region Air Temperature Degree, Δ t' are to calculate duration；

According to formula S M_t-1'=min { S_t-1,mlt',SWE_t-1' determine the practical amount of snowmelt amount SM in the trained region_t-1'； Wherein, the SWE_t-1' initial value be setting value；According to formula Calculate solid precipitation quantity P_s'；Wherein, T_s' be the trained region snowfall threshold temperature, T_r' be the trained region rainfall Threshold temperature, P_tot' be the trained region gross precipitation；

According to formula S WE_t'=SWE_t-1'-SM_t-1'+P_s' be calculated to convergent water equivalent of snow SWE_t', and judge institute It states to convergent water equivalent of snow SWE_t' actual measurement practical water equivalent of snow between relationship whether meet the default condition of convergence；

If it is not, then adjust the value of the snow melt parameter, obtain it is new to convergent water equivalent of snow, until described new due-in Until relationship between the water equivalent of snow held back and the practical water equivalent of snow of actual measurement meets the default condition of convergence；

Using the value of snow melt parameter adjusted as the actual value of the snow melt parameter.

The model training data according to known to and preset CREST model in one of the embodiments, determination are melted The actual value of icing parameter, comprising:

According to formulaThe training center is calculated Domain to convergent glacier melting amount GM_t'；Wherein, the SWE_t'、SM_t'、S_t,mlt' to utilize the first hydrological simulation model and institute State what the model training data in trained region were calculated, the first hydrological simulation model is will be in the CREST model The initial set value of snow melt parameter replaces with the model that the actual value of the snow melt parameter obtains later；

Judgement is described to convergent glacier melting amount GM_t' whether the relationship between the practical glacier melting amount of actual measurement full The default condition of convergence of foot；

If it is not, then adjust the value of the ice-melt parameter, obtain it is new to convergent glacier melting amount, until it is described it is new to Until relationship between convergent glacier melting amount and the practical glacier melting amount of actual measurement meets the default condition of convergence；

Using the value of ice-melt parameter adjusted as the actual value of the ice-melt parameter.

Second aspect, the embodiment of the present invention provide a kind of acquisition device of diameter flow data, comprising:

First determining module, for according to region to be measured driving data and preset hydrological simulation model, determine current Accumulated snow melt water amount, glacier melting amount and the snow cover maximum water-holding capacity at moment；The driving data in the region to be measured includes: default The air themperature of the precipitation in the region to be measured in the period, the surface temperature in region to be measured and region to be measured；It is described preset It include the parameter that can be used in characterizing the melt water characteristic of cryosphere in hydrological simulation model；

Second determining module, for by the solid precipitation quantity in the precipitation in the region to be measured and the region to be measured Difference, the liquid precipitation amount being determined as in the preset time period；

First processing module, for according to the freezing water quantity again of the water-holding capacity of previous moment, the snow cover of previous moment with And the water holding variable quantity of previous moment, obtain the water-holding capacity at current time；

Second processing module, for according to the accumulated snow melt water amount, the glacier melting amount, the liquid precipitation amount, institute The water-holding capacity for stating snow cover maximum water-holding capacity and the current time obtains the water holding variable quantity at current time；

Third processing module, for by the difference between ideal meltwater runoff and the water holding variable quantity at the current time, It is determined as the practical diameter flow data at current time.

It the third aspect, a kind of computer equipment provided in an embodiment of the present invention, including memory, processor and is stored in On reservoir and the computer program that can run on a processor, realized when the processor is for executing the computer program The step of stating method described in any embodiment.

A kind of fourth aspect, readable storage medium storing program for executing provided in an embodiment of the present invention, is stored thereon with computer program, described The step of method described in any of the above-described embodiment is realized when computer program is executed by processor.

Acquisition methods, device, computer equipment and the readable storage medium storing program for executing of diameter flow data provided by the invention, by The driving data in survey region and preset hydrological simulation model, determine accumulated snow melt water amount, glacier melting amount and the product at current time Snow deposit maximum water-holding capacity, and by the difference of the gross precipitation in region to be measured and the solid precipitation quantity in the region to be measured, it determines For the liquid precipitation amount in preset time period；Then freezed again according to the snow cover of the water-holding capacity of previous moment, previous moment The water holding variable quantity of water and previous moment obtains the water-holding capacity at current time, and is melted according to above-mentioned accumulated snow melt water amount, glacier Water, liquid precipitation amount, snow cover maximum water-holding capacity and the water-holding capacity at current time obtain the water holding variation at current time Amount, so that the difference between ideal meltwater runoff and the water holding variable quantity at the current time is determined as the reality at current time Border diameter flow data, the i.e. method of the present embodiment, when calculating the practical diameter flow data in region to be measured, by accumulated snow melt water and glacier The production stream that melt water is formed is incorporated into diameter flow data, so that the measuring method of the diameter flow data of cryosphere is more accurate；Meanwhile In practical calculating process, also using water-holding capacity and again freezing water quantity as the runoff Consideration in region to be measured, it greatly improves The counting accuracy of the diameter flow data in region to be measured.

Detailed description of the invention

Fig. 1 is the acquisition methods flow diagram of diameter flow data provided by one embodiment of the present invention；

Fig. 2 is the flow diagram of the acquisition methods for the diameter flow data that another embodiment of the present invention provides；

Fig. 3 is the flow diagram of the acquisition methods for the diameter flow data that another embodiment of the invention provides；

Fig. 4 is the flow diagram of the acquisition methods for the diameter flow data that another embodiment of the invention provides；

Fig. 5 is the flow diagram of the acquisition methods for the diameter flow data that another embodiment of the invention provides；

Fig. 6 is the flow diagram of the acquisition methods for the diameter flow data that another embodiment of the invention provides；

Fig. 7 is the flow diagram of the acquisition methods for the diameter flow data that another embodiment of the invention provides；

Fig. 8 is a kind of structural schematic diagram of the acquisition device of diameter flow data provided in an embodiment of the present invention；

Fig. 9 is the structural schematic diagram of the acquisition device of another diameter flow data provided in an embodiment of the present invention；

Figure 10 is the structural schematic diagram of the acquisition device of another diameter flow data provided in an embodiment of the present invention；

Figure 11 is the structural schematic diagram of the acquisition device of another diameter flow data provided in an embodiment of the present invention；

Figure 12 is the structural schematic diagram of the acquisition device of another diameter flow data provided in an embodiment of the present invention；

Figure 13 is the structural schematic diagram of the acquisition device of another diameter flow data provided in an embodiment of the present invention；

Figure 14 is the structural schematic diagram of the acquisition device of another diameter flow data provided in an embodiment of the present invention.

Specific embodiment

With the progress of mankind's science and technology, people have further deep understanding for the earth environment locating for oneself.Pass through To astronomical and meteorological further investigation, people have grasped the rule of climate change to a certain extent, with such rule Rule carrys out Instructing manufacture activity, while so that people can adapt to the natural law for the transformation of the Nature, preferably serves People's lives.Wherein, daily life is directly related to for the research of Water Resources Domain, and hydrology field is as water conservancy The important branch in field, research effect value can not be underestimated.

In hydrology field, important component of the cryosphere as weather system has climate change important It influences, usual people calculate the hydrologic process (such as diameter flow data) in region by establishing hydrological model, thus right Weather in region is assessed.But it can not be accurate by the diameter flow data that traditional hydrological model is predicted Reaction frost collar region practical runoff situation, and then cryosphere can not accurately be estimated to the actual influence of weather.

The acquisition methods of the diameter flow data of cryosphere provided in an embodiment of the present invention, device, computer equipment and readable Storage medium calculates the meltwater runoff data of cryosphere by the accurate hydrological model established, substantially increases ice Freeze the precision of the measuring and calculating of circle diameter flow data, and then improves estimation cryosphere to the accuracy of the actual influence of weather.

It should be noted that the executing subject of following methods embodiment can be the acquisition device of diameter flow data, the device The some or all of of computer equipment can be implemented as by way of software, hardware or software and hardware combining.Optionally, The computer equipment can be the electronic equipments having data processing function such as PC, portable device, server, the present embodiment pair The concrete form of computer equipment and without limitation.The executing subject of following methods embodiment is carried out by taking computer equipment as an example Explanation.

It is following to the present embodiments relate to technical term or noun explain or illustrate:

Cryosphere: at the earth's surface, generally 3000 meters of height above sea level or more and high latitude area, there is ice sheet and glacier covering, It is frozen soil layer under it, these ice sheets, glacier and frozen soil layer etc. is known as cryosphere.Can occur in cryosphere snowfall, snow melt with And the weather circle of at least one of ice-melt natural phenomena.

Water holding: water holding refers to the liquid water for being maintained in ice and snow layer and not formed and produce stream.

Runoff: runoff refers to rainfall, snowmelt runoff or in flowing water stream under the effect of gravity.

The freezing water quantity again of snow cover: the liquid water freezed again can occur in snow cover.

In order to make the objectives, technical solutions, and advantages of the present invention clearer, pass through following embodiments and combine attached Figure, the further description of technical solution in the embodiment of the present invention.It should be appreciated that specific embodiment described herein Only to explain the present invention, it is not intended to limit the present invention.

Fig. 1 is the acquisition methods flow diagram of diameter flow data provided by one embodiment of the present invention, and the present embodiment is related to Be accumulated snow melt water amount, glacier melting amount and snow cover maximum water-holding capacity according to identified current time, and according to institute The water-holding capacity of the snow cover of liquid precipitation amount and current time in determining preset time period calculates the practical diameter at current time The detailed process of flow data.As shown in Figure 1, this method comprises:

S101, the driving data according to region to be measured and preset hydrological simulation model, determine that the accumulated snow at current time melts Water, glacier melting amount and snow cover maximum water-holding capacity；The driving data in the region to be measured include: in preset time period to Survey gross precipitation, the surface temperature in region to be measured and the air themperature in region to be measured in region；The preset hydrological simulation mould Include the parameter for characterizing the melt water characteristic of cryosphere in type.

Specifically, the region to be measured or cryosphere in the present embodiment can be to contain the area on accumulated snow, ice sheet, land etc. Domain needs to utilize the hydrological simulation in the present embodiment when for the carry out melt water Predicting Performance Characteristics in the region to be measured of this type Model contains the parameter of the melt water characteristic for characterizing cryosphere in the model, such as can be the parameter of characterization amount of snowmelt amount, The parameter that can also be characterization ice storage can also be the combination of the two.In the present embodiment, in the practical runoff of prediction cryosphere When data, computer equipment according to the driving data in above-mentioned region to be measured and preset hydrological simulation model, determines current first Accumulated snow melt water amount, glacier melting amount and the snow cover maximum water-holding capacity at moment.Wherein, the driving data in above-mentioned region to be measured includes Precipitation, the air themperature of the surface temperature in region to be measured and region to be measured of region to be measured within a preset period of time.In general, working as When time step is using day as unit, the surface temperature in region to be measured can be mean daily temperature in the measuring and calculating time and to be measured The air themperature in region can be the mean daily temperature in the measuring and calculating time.Optionally, the region to be measured in above-mentioned preset time period Precipitation may include solid precipitation and liquid precipitation, wherein solid precipitation is often referred to snowfall, and liquid precipitation is often referred to drop Rain；The surface temperature in above-mentioned region to be measured and the air themperature in region to be measured are acquired to the actual temperature in region to be measured After obtaining raw temperature data, raw temperature data is obtained later by rising scale or NO emissions reduction processing.It needs to illustrate , since the raw temperature data that computer equipment acquires is the number being distributed according to the grid dividing in region to be measured According to, however the precision of the usually grid dividing is not necessarily matched with the precision that preset hydrological simulation model needs, it is therefore desirable to Raw temperature data is subjected to liter scale or NO emissions reduction processing to adapt to the needs of above-mentioned hydrological simulation model.For example, when obtaining Region to be measured surface temperature be 2km × 2km, and above-mentioned hydrological simulation model need data precision be 10km × 10km When, then it needs the surface temperature in region to be measured carrying out a liter scale processing, making as grid precision is 10km × 10km Surface temperature；Likewise, when the air themperature in the region to be measured obtained is 20km × 20km, and what above-mentioned simulation model needed When data precision is 10km × 10km, then need to carry out the surface temperature in region to be measured that scale processing is made grid Precision is the air themperature of 10km × 10km.It should be noted that due to its spatial resolution of the surface temperature in region to be measured compared with It is high, it is therefore desirable to liter scale processing be carried out to the surface temperature in region to be measured to adapt to the processing of above-mentioned hydrological simulation model and need It asks；And the air themperature in region to be measured is unable to satisfy use demand since its spatial resolution is low, it is therefore desirable to region to be measured Air themperature carry out NO emissions reduction processing to adapt to the process demand of above-mentioned hydrological simulation model.In general, to the to be measured of acquisition When zone air temperature carries out liter scale processing, since the zone air temperature to be measured of acquisition and region surface temperature to be measured have Same variation tendency, therefore zone air temperature to be measured processing can be carried out according to the regularity of distribution of region surface temperature to be measured NO emissions reduction is handled to adapt to the process demand of above-mentioned hydrological model.

In general, the surface temperature in above-mentioned region to be measured can be obtained by satellite, the air themperature in above-mentioned region to be measured can To be obtained by ground observation website, the specific acquisition modes embodiment of the present invention and with no restrictions.Optionally, satellite can be with Hydrographic data is obtained by the MODI sensor of loading.When the surface temperature in the region to be measured of acquisition and the air in region to be measured Temperature is passed through after data processing, be can be used as a part of driving data, is input to above-mentioned preset hydrological simulation model, with It is exported accordingly.

Due to containing the parameter of the melt water characteristic for characterizing cryosphere in above-mentioned preset hydrological simulation model, When input of the driving data in region to be measured as the model, which can be according to the melt water characteristic of the characterization cryosphere Parameter predict that region to be measured, obtaining region current time to be measured obtains accumulated snow melt water amount, glacier melting amount and accumulated snow Layer maximum water-holding capacity, specific simulation and forecast process the present embodiment and without limitation, as long as this can be utilized to contain above-mentioned melt The hydrological simulation model of the parameter of water characteristic obtains the accumulated snow melt water amount, glacier melting amount and snow cover maximum water holding in region to be measured Amount.

S102, by the difference of the gross precipitation in the region to be measured and the solid precipitation quantity in the region to be measured, determine For the liquid precipitation amount in the preset time period.

Specifically, as described above, for the precipitation in a region usually there are two types of form, one kind is that solid-state drops Water, i.e. snowfall；Another is liquid precipitation, i.e. rainfall.Therefore, the gross precipitation in above-mentioned region to be measured is subtracted into region to be measured Solid precipitation quantity, the liquid precipitation amount of region to be measured within a preset period of time can be obtained.Wherein, solid precipitation quantity and liquid The accounting of precipitation is influenced and different by the surface temperature in region to be measured and the air themperature in region to be measured.

S103, according to the water-holding capacity of previous moment, the freezing water quantity again of the snow cover of previous moment and previous moment Water holding variable quantity obtains the water-holding capacity at current time.

Usually in cryosphere, a part of accumulated snow melt water, glacier melting and liquid precipitation can be retained in accumulated snow and glacier is worked as In, this operative liquid water is known as water holding by us.Since the water holding in accumulated snow and glacier will not generate runoff, and snowmelt runoff Also can there is a situation where freeze again, it is therefore desirable to using water-holding capacity and again freezing water quantity as the reference of the diameter flow data in region to be measured Factor.That is, ideally, if accumulated snow or glacier are there is no water holding and the case where freeze again, snow cover is practical The water of thawing adds the runoff water that the practical water melted in glacier can be region to be measured.But just because of accumulated snow or Person glacier device in actual scene can freeze there is a situation where water holding and again, and therefore, the present embodiment freezes by water-holding capacity and again Runoff Consideration of the water as region to be measured, it greatly improves the counting accuracies of the diameter flow data in region to be measured.

Therefore, in this step, computer equipment is melted in the accumulated snow melt water amount at current time, glacier obtaining region to be measured After liquid precipitation amount in water, snow cover maximum water-holding capacity and preset time period, computer equipment be will continue to before The water-holding capacity at one moment can be in conjunction with the freezing water quantity again of the snow cover of previous moment and the water holding variable quantity of previous moment Obtain the water-holding capacity at current time.

It is S104, maximum according to the accumulated snow melt water amount, the glacier melting amount, the liquid precipitation amount, the snow cover Water-holding capacity and the water-holding capacity at the current time obtain the water holding variable quantity at current time.

S105, by the difference between ideal meltwater runoff and the water holding variable quantity at the current time, when being determined as current The practical diameter flow data carved.

Specifically, comprehensively consider accumulated snow melt water amount, glacier melting amount, liquid precipitation amount and snow cover maximum water-holding capacity, In conjunction with the water-holding capacity at current time, then water holding variable quantity of the available current time relative to last moment utilizes ideal Meltwater runoff subtracts the water holding variable quantity, and obtained difference is determined as to the practical diameter flow data at current time.It needs Bright, above-mentioned ideal meltwater runoff is the meltwater runoff for not considering the cryosphere of ice and snow layer water holding variation.

The acquisition methods of diameter flow data provided in this embodiment, driving data and preset hydrology mould by region to be measured Analog model, determines accumulated snow melt water amount, glacier melting amount and the snow cover maximum water-holding capacity at current time, and by the total of region to be measured The difference of the solid precipitation quantity in precipitation and the region to be measured, the liquid precipitation amount being determined as in preset time period；Then According to the freezing water quantity again of the snow cover of the water-holding capacity of previous moment, previous moment and the water holding variable quantity of previous moment, obtain To the water-holding capacity at current time, and according to above-mentioned accumulated snow melt water amount, glacier melting amount, liquid precipitation amount, snow cover maximum water holding Amount and the water-holding capacity at current time, obtain the water holding variable quantity at current time, thus by ideal meltwater runoff and described current Difference between the water holding variable quantity at moment is determined as the practical diameter flow data at current time, the i.e. method of the present embodiment, is counting When calculating the practical diameter flow data in region to be measured, the production stream that accumulated snow melt water and glacier melting are formed is incorporated into diameter flow data, So that the measuring method of the diameter flow data of cryosphere is more accurate；Meanwhile in practical calculating process, also freeze by water-holding capacity and again Runoff Consideration of the amount of bearing water as region to be measured, it greatly improves the counting accuracies of the diameter flow data in region to be measured.

Fig. 2 is the flow diagram of the acquisition methods for the diameter flow data that another embodiment of the present invention provides.The present embodiment What is involved is according to region to be measured driving data and preset hydrological simulation model, determine current time accumulated snow melt water amount, The detailed process of glacier melting amount and snow cover maximum water-holding capacity.On the basis of the above embodiments, optionally, as shown in Fig. 2, Above-mentioned S101 may include steps of:

S201, according to the snow melt parameter in the driving data and preset hydrological simulation model in the region to be measured, determine The possibility amount of snowmelt amount of the water equivalent of snow at current time and the region to be measured.

Wherein, the snow melt parameter includes environmental parameter for characterizing precipitation or snow melting and for characterizing The parameter of amount of snowmelt amount under unit time unit temperature.

Specifically, being described by above-described embodiment it is found that comprising for characterizing frost in above-mentioned preset hydrological simulation model The parameter of the melt water characteristic of circle, the parameter of the melt water characteristic may include snow melt parameter and/or ice-melt parameter.Join for the snow melt Number, optionally, which may include: the environmental parameter for characterizing precipitation or snow melting, and for characterizing The parameter of amount of snowmelt amount under unit time unit temperature.Optionally, the above-mentioned environment ginseng for characterizing precipitation or snow melting Counting this for example can be snow melt threshold temperature, snowfall threshold temperature and rainfall threshold temperature；It is above-mentioned to be used to characterize the unit time The parameter of amount of snowmelt amount under unit temperature for example can be maximum snow melt degree-day factor and the minimum snow melt of winter solstice of summer solstice Degree-day factor.Wherein, snow melt threshold temperature is that accumulated snow starts the critical-temperature melted, the value range of usual snow melt threshold temperature It can be between -5 DEG C~+5 DEG C；Snowfall threshold temperature is the critical-temperature that snowfall takes place, usual snowfall threshold temperature Value range can be between -5 DEG C~+5 DEG C；Rainfall threshold temperature is the critical-temperature that rainfall takes place, usual rainfall threshold The value range for being worth temperature can be between -5 DEG C~+5 DEG C；The maximum snow melt degree-day factor of summer solstice is the earth's surface in region to be measured The average value of temperature and the air themperature in region to be measured had often been higher by snow melt threshold temperature once, in melting for every day time-division summer solstice Snowfall, the minimum snow melt degree-day factor of winter solstice are the average value of the surface temperature in region to be measured and the air themperature in region to be measured Snow melt threshold temperature was often higher by once, in the amount of snowmelt amount of every day time-division winter solstice, wherein since the temperature of summer solstice is high, because The snow melt degree-day factor of summer solstice is considered the maximum snow melt degree-day factor in whole year by this, and usually maximum snow melt degree-day factor takes Value may range from 0mm DEG C^-1d^-1~10mm DEG C^-1d^-1Between；And Winter Solstice degree/day is low, therefore by the snow melt of winter solstice subsist because Son is considered the minimum snow melt degree-day factor in whole year, and usually minimum snow melt degree-day factor value range can be 0mm DEG C^-1d^-1 ~10mm DEG C^-1d^-1Between.

Based on the snow melt parameter in above-mentioned hydrological simulation model, computer equipment can be according to the driving data in region to be measured With the snow melt parameter, the water equivalent of snow at current time and the possibility amount of snowmelt amount in region to be measured are determined.Current time is determined specific Possibility amount of snowmelt amount when, detailed process can be with are as follows:

Computer equipment first according to the maximum snow melt degree-day factor of summer solstice and winter solstice most Small snow melt degree-day factor obtains snow melt degree-day factor, and possible implementation can pass through formulaOr the modification of the formula acquires.Wherein, d_s,mltFor snow melt degree-day factor, d_s,mlt6For the maximum snow melt degree-day factor of summer solstice, d_s,mlt12For the minimum snow melt degree of winter solstice Day factor, and d_s,mlt6And d_s,mlt12For parameter preset in model, above-mentioned doy is number of days, can according to need measuring and calculating Time determines.

Secondly, after computer equipment has determined snow melt degree-day factor, computer equipment can be subsisted according to snow melt because What the surface temperature in region sub, to be measured, the air themperature in region to be measured and above-mentioned snow melt threshold temperature obtained region to be measured can Energy amount of snowmelt amount, possible implementation method can pass through formulaOr the change of the formula Type acquires.Wherein, S_s,mltFor possible amount of snowmelt amount, T_lsFor the surface temperature in region to be measured, T_airmxFor the Air Temperature in region to be measured Degree, T_s,mltFor snow melt threshold temperature, Δ t is to calculate duration, and T_lsAnd T_airmxIt can be by satellite or ground observation website Obtained temperature data is handled to obtain.

In addition, detailed process can be with are as follows: computer when computer equipment specifically determines the water equivalent of snow at current time Equipment by the difference of accumulated snow melt water amount known to water equivalent of snow known to previous moment and previous moment, with solid-state snowfall and Value is determined as the water equivalent of snow at current time, and possible implementation method can pass through formula S WE_t=SWE_t-1-SM_t-1+P_s,t-1 Or the modification of the formula acquires.Wherein, SWE_tFor the water equivalent of snow at current time, SWE_t-1Then work as the snow-broth of previous moment Amount, SM_t-1For the accumulated snow melt water amount of previous moment, P_s,t-1For the solid precipitation quantity of previous moment.Optionally, solid precipitation is obtained The possible implementation of amount can pass through formulaOr the formula Modification acquires, wherein P_sFor the solid precipitation quantity in the region to be measured of t moment, wherein T_rFor the rainfall threshold temperature in region to be measured； T_sFor the snowfall threshold temperature in region to be measured, and T_rAnd T_sFor data preset in model, usual T_rAnd T_sValue range be -5 DEG C -- between+5 DEG C；P_totFor the gross precipitation in region to be measured, and P_totThe reality observed for satellite or ground observation website Data.

It should be noted that the water equivalent of snow and accumulated snow melt water amount at current time all rely on the data of previous moment, It is that the data of t moment are obtained by the superposition at moment.Under normal conditions, it is calculated since June, since there is not accumulated snow in June, Therefore 0 can be set by the water equivalent of snow of initial time, with the variation of time, temperature is gradually decreased, along with rainfall, drop The generation of phenomena such as snow, is superimposed by the moment, obtains the water equivalent of snow and accumulated snow melt water amount at current time by above-mentioned formula.

S202, by the minimum value between the possibility amount of snowmelt amount and the water equivalent of snow in the region to be measured, be determined as described Accumulated snow melt water amount.

Specifically, above-mentioned accumulated snow melt water amount can pass through formula S M_t=min { S_t,mlt,SWE_tObtain, wherein SM_tFor product Snowmelt amount.

S203, worked as according to the possibility amount of snowmelt amount in the region to be measured, the accumulated snow melt water amount in the region to be measured, the snow-broth Ice-melt parameter in amount and the hydrological simulation model, determines the glacier melting amount.

Wherein, the ice-melt parameter includes environmental parameter for characterizing Melting Glacierss and for characterizing the unit time The parameter of ice storage under unit temperature.

Specifically, being described by above-described embodiment it is found that comprising for characterizing frost in above-mentioned preset hydrological simulation model The parameter of the melt water characteristic of circle, the parameter of the melt water characteristic may include snow melt parameter and/or ice-melt parameter.Join for the snow melt Number, optionally, the ice-melt parameter in the region to be measured may include: the environmental parameter for characterizing Melting Glacierss, and be used for table Levy the parameter of the ice storage under unit time unit temperature.Optionally, the above-mentioned environmental parameter for characterizing Melting Glacierss can be with It is ice-melt threshold temperature；The above-mentioned parameter for characterizing ice storage under unit time unit temperature can be ice-melt subsist because Son.Wherein, ice-melt threshold temperature is that glacier starts the critical-temperature melted, and ice-melt degree-day factor is the surface temperature in region to be measured Ice-melt threshold temperature had often been higher by once with the average value of the air themperature in region to be measured, the ice storage of the region every day to be measured. It should be noted that under normal conditions, the melt temperature of accumulated snow is high compared with the temperature of Melting Glacierss, it is therefore, when temperature is raised, first First start to melt is accumulated snow, and secondly glacier just starts to melt.

Specifically, when computer equipment obtain the possibility amount of snowmelt amount in region to be measured, region to be measured accumulated snow melt water amount and After the water equivalent of snow at region current time to be measured, computer equipment can be according to the possibility in identified region to be measured Ice-melt parameter in amount of snowmelt amount, the accumulated snow melt water amount in region to be measured, water equivalent of snow and above-mentioned hydrological simulation model, Determine the glacier melting amount at region current time to be measured, possible implementation can be to pass through formulaOr the modification of the formula acquires, wherein GM_tMelt for glacier Water, d_i,mltFor ice-melt degree-day factor, T_i,mltFor ice-melt threshold temperature.

S204, by the product of the water equivalent of snow and the preset water holding factor, be determined as the snow cover maximum water-holding capacity； Wherein, the liquid water that the water holding factor is held by unit water equivalent of snow.

Specifically, computer equipment will be worked as after computer equipment obtains the water equivalent of snow at region current time to be measured The water equivalent of snow SWE at preceding moment_tWith the product of preset water holding factor HWC, it is determined as snow cover maximum water-holding capacity.Optionally, Snow cover maximum water-holding capacity MHW_tIt can be according to formula MHW_t=HWCSWE_tOr the modification of the formula acquires.Wherein, HWC is The preset water holding factor, the liquid water that water holding factor characterization unit water equivalent of snow is held, in general, the water holding factor takes Value may range between 0.05~0.15.

Method provided in this embodiment, computer equipment pass through driving data according to region to be measured and preset hydrology mould Snow melt parameter in analog model determines the water equivalent of snow at current time and the possibility amount of snowmelt amount in region to be measured, then by area to be measured Minimum value between the possibility amount of snowmelt amount and water equivalent of snow in domain is determined as above-mentioned accumulated snow melt water amount, and on this basis, computer Equipment is according to the possibility amount of snowmelt amount in the region to be measured of above-mentioned determination, the accumulated snow melt water amount in region to be measured, water equivalent of snow and above-mentioned Ice-melt parameter in hydrological simulation model determines glacier melting amount, meanwhile, by identified water equivalent of snow and preset water holding because The product of son, is determined as snow cover maximum water-holding capacity.That is, the method for the present embodiment, based on identified region to be measured Accumulated snow melt water amount, glacier melting amount and the snow cover maximum water-holding capacity at current time calculate the practical diameter fluxion in region to be measured According to the production stream that accumulated snow melt water and glacier melting are formed is incorporated into diameter flow data by the present embodiment, so that the diameter fluxion of cryosphere According to measuring method it is more accurate.

Fig. 3 is the flow diagram of the acquisition methods for the diameter flow data that another embodiment of the invention provides.The present embodiment What is involved is computer equipments how according to the water-holding capacity of previous moment, the freezing water quantity again of the snow cover of previous moment and preceding The water holding variable quantity at one moment, obtains the detailed process of the water-holding capacity at current time.On the basis of the above embodiments, optional , as shown in figure 3, above-mentioned S103 may include steps of:

S301, according to formula

Determine the possibility of previous moment freezing water quantity W again_t-1,ref；Wherein, the RFC be preset water holding freeze again because Son, the T_i,mltFor the ice-melt threshold temperature in the hydrological simulation model, the d_i,mltFor in the hydrological simulation model Ice-melt degree-day factor, the Δ t are to calculate duration.

It should be noted that freeze-out fraction is the preset value of above-mentioned hydrological simulation model to preset water holding again, the water holding is again Freeze-out fraction characterizes the water holding in ice and snow layer and glaciated ratio occurs again.Under normal conditions, above-mentioned water holding freezes again The value range of the factor can be set between 0.05~0.15；The value range of above-mentioned ice-melt threshold temperature can be set -5 DEG C~+5 DEG C between；The value range of ice-melt degree-day factor can be set at 0mm DEG C^-1d^-1~10mm DEG C^-1d^-1Between, wherein it is single Mm DEG C of position^-1d^-1It indicates the ice storage of every degree daily, and calculates duration usually as unit of day.

Specifically, the possibility of previous moment freezing water quantity W again_t-1,refIt can be according to formulaAnd the modification of the formula acquires, each parameter in the formula is specifically It is bright referring to above, details are not described herein again.

S302, by the possibility of previous moment freezing water quantity W again_t-1,refBetween the water-holding capacity of the previous moment Minimum value is determined as the freezing water quantity again of the snow cover of previous moment.

Since snowmelt runoff produces stream and liquid water holding in addition to will form in ice and snow layer, some in actual conditions Can there is a phenomenon where freeze again, it is therefore desirable to by ice and snow layer it is total in liquid water freezing water quantity again in view of in diameter stream calculation.

Specifically, by the possibility of previous moment freezing water quantity W again_t-1,refWith the water-holding capacity HW of previous moment_t-1In minimum Value, is determined as the freezing water quantity RF again of the snow cover of previous moment_t-1, wherein the freezing water quantity again of the snow cover of previous moment with And the water-holding capacity of previous moment is for current time.At current time, it is generally recognized that the snow cover of previous moment Freezing water quantity again and the water-holding capacity of previous moment be the given data that can be calculated by model.Therefore, when previous The freezing water quantity RF again of the snow cover at quarter_t-1Formula RF can be passed through_t-1=min { HW_t-1,W_t-1,refIndicate.

S303, by it is between the water-holding capacity of the previous moment and the water holding variable quantity of previous moment and value, subtract previous Obtained difference after the freezing water quantity again of the snow cover at moment, is determined as the water-holding capacity at the current time.

Specifically, the water-holding capacity HW of previous moment_t-1With the water holding variation delta HW of previous moment_t-1Between and value, subtract Remove the freezing water quantity RF again of the snow cover of previous moment_t-1Obtained difference characterizes holding for the current time of snow cover later Water HW_t, wherein the water-holding capacity of previous moment and the water holding variable quantity of previous moment are for current time.? Current time, it is generally recognized that the water holding variable quantity of the water-holding capacity of previous moment and previous moment is that can be calculated by model The given data arrived.Therefore, the water-holding capacity HW at current time_tFormula HW can be passed through_t=HW_t-1+ΔHW_t-1-RF_t-1Or the public affairs The modification of formula acquires.

It should be noted that the water-holding capacity at current time depends on the water-holding capacity data of previous moment, to pass through the moment The obtained data of superposition.Under normal conditions, it is calculated since June, since there is not accumulated snow in June, there will not be product The water holding of snow deposit, so 0 is set by the water-holding capacity of initial time.With the variation of time, temperature is gradually decreased, along with drop The generation of phenomena such as rain, snowfall, by the superposition at moment, the water-holding capacity at available current time.

In the embodiment shown in fig. 3, by determining the possibility product of freezing water quantity and previous moment again of previous moment The freezing water quantity again of snow deposit, so that it is determined that the water-holding capacity at current time so that by the possibility of snow cover again freezing water quantity and The water-holding capacity of ice and snow layer is added in the calculating for producing stream, as the runoff Consideration in region to be measured, so that cryosphere The estimation precision of diameter flow data is more accurate, preferably people's production and life can be instructed to live.

Fig. 4 is the flow diagram of the acquisition methods for the diameter flow data that another embodiment of the invention provides.The present embodiment What is involved is computer equipments how according to above-mentioned identified accumulated snow melt water amount, glacier melting amount, liquid precipitation amount, snow cover Maximum water-holding capacity and the water-holding capacity at current time obtain the detailed process of the water holding variable quantity at current time.Any of the above-described On the basis of embodiment, optionally, above-mentioned S104 be may include steps of:

S401, determine the accumulated snow melt water amount, the glacier melting amount and the liquid precipitation amount and value.

Specifically, when computer equipment obtains the accumulated snow melt water amount SM in region to be measured by the process of above-described embodiment_t, ice River melt water amount GM_t, liquid precipitation amount P_t,l, snow cover maximum water-holding capacity MHW_tAnd the water-holding capacity HW at current time_tLater, it calculates Machine equipment calculates the accumulated snow melt water amount SM in region to be measured_t, glacier melting amount GM_tWith liquid precipitation amount P_t,lAnd value.By this and Value, characterizes the generation total amount of the liquid water in region to be measured.

S402, difference between the snow cover maximum water-holding capacity and the water-holding capacity at the current time is determined.

Specifically, the computer equipment snow cover maximum water-holding capacity MHW according to region to be measured_tWith the water holding at current time Measure HW_tDifference is calculated.By the difference, characterizes and reach what water holding saturation state also needed in the snow cover in region to be measured Liquid water.

Minimum value described in S403, determination between value and the difference is the water holding variable quantity at the current time.

Specifically, the water holding that the above-mentioned minimum value between value and difference is determined as the current time in region to be measured is changed Measure HW_t, the water holding variation delta HW at the region current time to be measured_tFormula Δ HW can be passed through_t=min { SM_t+GM_t+P_t,l, MHW_t-HW_tIndicate, wherein Δ HW_tCharacterize variable quantity of the current time relative to the water-holding capacity of the snow cover of previous moment.

It should be noted that two steps of S401 and S402 in the present embodiment, above-mentioned for determination and value and difference elder generation Afterwards sequence and without limitation.

In the present embodiment, computer equipment pass through the accumulated snow melt water amount that determines, glacier melting amount and liquid precipitation amount and value, And the difference between determining snow cover maximum water-holding capacity and the water-holding capacity at current time, and by above-mentioned between value and difference Minimum value is determined as the water holding variable quantity at current time, so that the water holding variable quantity of ice and snow layer is added in the calculating for producing stream, So that the measuring method of the diameter flow data of cryosphere is more accurate, preferably people's production and life can be instructed to live.

The various embodiments described above list respectively computer equipment how driving data according to region to be measured and preset water Literary simulation model obtains the diameter flow data in region to be measured, by the ginseng that cryosphere melt water characteristic is arranged in hydrological simulation model Number, so that the measuring method of the diameter flow data of cryosphere is more accurate.Next, describing the hydrology mould by following embodiments The training process of analog model.It should be noted that following methods limit for explaining only the invention rather than to the present invention It is fixed.

Fig. 5 is the flow diagram of the acquisition methods for the diameter flow data that another embodiment of the invention provides.The present embodiment What is involved is the detailed processes how computer equipment determines above-mentioned hydrological simulation model.It is above-mentioned apply example on the basis of, it is optional , can also include: before S101

S501, according to known model training data and preset CREST model, determine the actual value of snow melt parameter and melt The actual value of icing parameter；Wherein, the initial set value in the CREST model including snow melt parameter and the initial of ice-melt parameter set Definite value.Wherein, the model training data include: the gross precipitation in trained region, the surface temperature in the trained region, training The air themperature in region, the practical water equivalent of snow of actual measurement, actual measurement practical glacier melting amount.

It should be noted that the initial set value of above-mentioned snow melt parameter and the initial set value of ice-melt parameter are model specification Initial value, which can be the empirical value with reference to previous model, is also possible to use for reference other correlation models and obtains , without limitation to this present invention.

Specifically, computer equipment will be preset by satellite or the acquisition model training data input of ground observation website CREST model is calculated by the model, and is exported to convergent melt water index；Wherein, include to convergent melt water index To convergent snow melt parameter and to convergent ice-melt parameter.Further, computer equipment is due-in by exporting the model Correlation between the melt water index held back and the melt water index of actual measurement carries out convergence restriction, the actual value of available snow melt parameter With the actual value of ice-melt parameter.Optionally, the calculating process of the actual value of the actual value of snow melt parameter and ice-melt parameter, Ke Yican See below the description of Fig. 6 and embodiment illustrated in fig. 7.

In addition, the practical glacier melting amount of the practical water equivalent of snow and actual measurement surveyed in the training data of above-mentioned model can be with It surveys to obtain by satellite, can also survey to obtain by ground observation website, in this regard, the present invention is without limitation.In training data The air themperature of the gross precipitation in training region, the surface temperature in training region and training region is known data.

S502, the reality that the initial set value of the snow melt parameter in the CREST model is replaced with to the snow melt parameter It is worth, and the initial set value of the ice-melt parameter in the CREST model is replaced with to the actual value of the ice-melt parameter, obtains The hydrological simulation model.

In the present embodiment, according to known model training data and preset CREST model, the reality of snow melt parameter is determined The actual value of actual value and ice-melt parameter, and the initial set value of the snow melt parameter in CREST model is replaced with into snow melt parameter Actual value, and the initial set value of the ice-melt parameter in CREST model is replaced with into the actual value of ice-melt parameter, to obtain New hydrological simulation model.The actual value of snow melt parameter in the new hydrological model, the actual value of ice-melt parameter are with reality The glacier water capacity of the practical water equivalent of snow and actual measurement surveyed is data obtained from training objective, therefore in new hydrological model The actual value of snow melt parameter and the actual value of ice-melt parameter are more accurate.

The present embodiment is by determining the actual value of more accurate snow melt parameter and the actual value of ice-melt parameter, so that frost Natural phenomena in circle, such as: the factors such as freezing water quantity can pass through again for snowfall, snow melt, ice-melt, ice and snow water holding and water holding CREST model quantifies into the measuring and calculating of diameter flow data, so that the diameter flow data of the cryosphere of the CREST calculation using models is more quasi- Really, preferably people's production and life can be instructed to live.

Fig. 6 is the flow diagram of the acquisition methods for the diameter flow data that another embodiment of the invention provides.Fig. 6 is Fig. 5 Shown in embodiment according to known model training data and preset CREST model, determine the actual value of snow melt parameter A kind of possible implementation.Optionally, as shown in fig. 6, on the basis of above-mentioned embodiment shown in fig. 5, in S501 " according to Known model training data and preset CREST model, determine the actual value of snow melt parameter " the step of may include:

S601, according to formula

Determine the possibility amount of snowmelt amount S in the trained region_t-1,mlt'；

Wherein, snow melt degree-day factor Doy is number of days, T_s,mlt' it is snow melt threshold temperature, d_s,mlt6' be the summer solstice maximum snow melt degree-day factor, d_s,mlt12' it is Winter Solstice The minimum snow melt degree-day factor of day, T_ls' be the trained region earth's surface average daily temperature, T_airmx' be the trained region sky Gas average daily temperature, Δ t' are to calculate duration.

It should be noted that the snow melt degree-day factor can be according to formulaOr the modification of the formula acquires, and passes through The snow melt degree-day factor and above-mentioned preset CREST model acquired, further according to the change of above-mentioned formula (2) or the formula Type acquires may amount of snowmelt amount.In addition, the specific paraphrase of each parameter and value range are with reference to described previously in formula, herein not It repeats again.

S602, according to formula

SM_t-1'=min { S_t-1,mlt',SWE_t-1'} (3)

Determine the practical amount of snowmelt amount SM in the trained region_t-1'；Wherein, the SWE_t-1' initial value be setting value, T_ls' For the surface temperature in the trained region, T_airmx' be the trained region air themperature.It should be noted that due to snow-broth Equivalent is calculated since June, and snowfall does not occur June, therefore can be set the water equivalent of snow of initial time to 0。

S603, according to formula

Calculate solid precipitation quantity P_s'。

Specifically, computer equipment can calculate solid precipitation quantity according to the modification of above-mentioned formula (4) or the formula P_s', wherein P_tot' it is the gross precipitation for training region, P_tot' arrived for the preset value of above-mentioned model, value for moonscope Or ground observation website observes the actual value of the gross precipitation got.It should be noted that the possibility in training region is melted Minimum value in the water equivalent of snow in snowfall and training region is determined as training the practical amount of snowmelt amount in region, can pass through formula SM_t-1'=min { S_t-1,mlt',SWE_t-1' indicate.Optionally, the air of the surface temperature in above-mentioned trained region and training region The data that temperature can obtain for satellite or ground observation website, or the original number that satellite or ground observation website obtain According to obtaining by processing, for example, the data that satellite or ground observation website are obtained according to the demand of preset CREST model into Row rises scale or NO emissions reduction obtains；Alternatively, can also be the sequence of the more degree/days obtained to satellite or ground observation website It is averaged, i.e. T_ls' it can be the ground average daily temperature and T in training center domain_airmx' it is air average daily temperature.

S604, according to formula

SWE_t'=SWE_t-1'-SM_t-1'+P_s' (5)

It is calculated to convergent water equivalent of snow SWE_t', and judge described to convergent water equivalent of snow SWE_t' and actual measurement Whether the relationship between practical water equivalent of snow meets the default condition of convergence.

S605, if it is not, then adjust the value of the snow melt parameter, obtain it is new to convergent water equivalent of snow, until described new Meet the default condition of convergence to the relationship between convergent water equivalent of snow and the practical water equivalent of snow of actual measurement until, and will adjustment Actual value of the value of snow melt parameter afterwards as the snow melt parameter.

Under normal conditions, preset CREST model calculating is calculated since June, since June does not drop Snow, therefore 0 can be set by the water equivalent of snow of initial time.Over time, temperature gradually decreases, along with drop The generation of phenomena such as rain, snowfall.After gradually generating accumulated snow and freezing, then by the solid precipitation quantity of previous moment, previous The practical amount of snowmelt amount of moment water equivalent of snow and previous moment, determines the water equivalent of snow at current time, specifically: computer equipment According to formula S WE_t'=SWE_t-1'-SM_t-1'+P_s' or the modification of the formula acquire to convergent water equivalent of snow SWE_t', and sentence It is disconnected to be somebody's turn to do to convergent water equivalent of snow SWE_t' actual measurement practical water equivalent of snow between relationship whether meet the default condition of convergence. Wherein, each parameter indicates the parameter in training region in the formula, specific paraphrase and value range with reference to described previously, Details are not described herein again.

When the output of above-mentioned CREST model is to convergent water equivalent of snow SWE_t' between the practical water equivalent of snow of actual measurement When relationship is unsatisfactory for the preset condition of convergence, then the value of the snow melt parameter in above-mentioned CREST model is adjusted, thus based on adjustment CREST model afterwards be calculated it is new to convergent water equivalent of snow, and further judgement it is new to convergent water equivalent of snow with Whether the relationship between the practical water equivalent of snow of actual measurement meets the preset condition of convergence；If new to convergent water equivalent of snow and reality The relationship between practical water equivalent of snow surveyed still not satisfies the preset condition of convergence, then continues to adjust in above-mentioned CREST model The value of snow melt parameter, until what is obtained again new reaches pre- to convergent water equivalent of snow and the practical water equivalent of snow of actual measurement If the condition of convergence, then stop adjusting, and using the value of snow melt parameter adjusted as the actual value of snow melt parameter.Optionally, The initial set value of snow melt parameter in above-mentioned CREST model is replaced with obtained after the actual value of the snow melt parameter Model is the first hydrological simulation model, for the use when determining the actual value of ice-melt parameter of following embodiment illustrated in fig. 7.

Optionally, the above-mentioned preset condition of convergence can be so as to convergent water equivalent of snow and satellite or ground observation Objective function between the water equivalent of snow of website actual measurement reaches the preset condition of convergence.Optionally, above-mentioned objective function can wrap It includes: that assorted efficiency function NSE, relative error Bias and index of correlation CC；And the above-mentioned preset condition of convergence can be with are as follows: makes The value for obtaining above-mentioned objective function reaches ideal value or ideal convergence range.Wherein, receive assorted efficiency function NSE ideal value be 1, The ideal value of relative error function Bias is 0, and the ideal value of correlation of indices function CC is 1.Specifically, above-mentioned ideal convergence range It can be a range near the ideal value of each objective function, the setting of the range, which can according to need, to be selected, such as: The ideal convergence range of assorted efficiency function NSE received can be set between 0.8~1, the ideal convergence of relative error function Bias Range can be set between 0 ± 0.5, and the ideal convergence range of correlation of indices function CC can be set between 1 ± 0.3.Into One step, when needing the data precision calculated high, ideal convergence range can be arranged it is smaller, to meet high accuracy data Measuring and calculating demand；When needing the data precision calculated low, ideal convergence range can be arranged larger, meet data essence The efficiency of data measuring and calculating is improved in the case where degree.

It should be noted that above-mentioned objective function is embodied as: that assorted efficiency functionRelative error functionAnd the index of correlationWherein, OBSi is the reality of satellite or ground observation website actual observation Border data, S IMi are that preset CREST simulation model calculates exporting as a result, in the present embodiment, OBSi is satellite or ground The practical water equivalent of snow of website actual measurement is observed, SIMi works as what preset CREST simulation model was calculated to convergent snow-broth Measure SWE_t'。

The convergence process of above-mentioned model in order to facilitate understanding, it is following to cite a plain example: assuming that for the first time according to pre- If CREST model and known model training data obtain for the first time to convergent water equivalent of snow SWE_t' value be A, calculate SIMi in above three objective function is replaced with A by machine equipment, and OBSi is then replaced with satellite or ground observation website is real The practical water equivalent of snow surveyed, the value of three objective functions is calculated, while judging whether the value of three objective functions is distinguished Reach respective ideal value or ideal convergence range.If do not reached, the snow melt parameter in above-mentioned CREST model is adjusted Value, continue through CREST model adjusted and known model training data obtain second to convergent water equivalent of snow SWE_t', it is assumed that the value is B, and the SIMi in above three objective function is replaced with B, then replaces OBSi by computer equipment For the practical water equivalent of snow that satellite or ground observation website are surveyed, the value of three objective functions is calculated, while three are judged Whether the value of a objective function respectively reaches respective ideal value or ideal convergence range；If reaching, after above-mentioned adjustment Snow melt parameter actual value of the value as snow melt parameter.

In the present embodiment, by constantly adjusting the snow melt parameter in CREST model, so that the snow-broth of simulation model output Equivalent and the practical water equivalent of snow of actual measurement meet the preset condition of convergence, so that the snow melt parameter in CREST model is more accurate, To quantify the influence of the factors diameter flow data such as snowfall, snow melt and ice and snow water holding into the measuring and calculating of diameter flow data, so that The measuring and calculating of the diameter flow data of cryosphere is more accurate, so as to preferably instruct people's lives and production.

Fig. 7 is the flow diagram of the acquisition methods for the diameter flow data that another embodiment of the invention provides.Fig. 7 is Fig. 6 Shown in embodiment according to known model training data and preset CREST model, determine the actual value of ice-melt parameter A kind of possible embodiment.Optionally, as shown in fig. 7, on the basis of above-mentioned embodiment shown in fig. 6, after S605 also May include:

S701, according to formula

Be calculated the trained region to convergent glacier melting amount GM_t'；Wherein, the SWE_t'、SM_t'、S_t,mlt' For what is be calculated using the model training data of the first hydrological simulation model and the trained region, first hydrological simulation After actual value of the model for the initial set value of the snow melt parameter in the CREST model to be replaced with to the snow melt parameter The model arrived.

Specifically, the description based on above-mentioned embodiment illustrated in fig. 6, computer equipment is by the above-mentioned reality comprising snow melt parameter The hydrological model of value passes through the first hydrological simulation model, the model in the training region as the first hydrological simulation model The glacier melting amount GM in the trained region is calculated in the deformation of training data and above-mentioned formula (6) or the formula_t'。

Glacier melting amount GM described in S702, judgement to convergent trained region_t' with actual measurement practical glacier melting amount it Between relationship whether meet the default condition of convergence.

S703, if it is not, then adjust the value of the ice-melt parameter, obtain it is new to convergent glacier melting amount, until described New meets default convergence item to the relationship between convergent glacier melting amount and the practical glacier melting amount surveyed by satellite Until part, and using the value of ice-melt parameter adjusted as the actual value of the ice-melt parameter.

It is carried out on the basis of snow melt parameter determines it should be noted that the determination of ice-melt parameter can establish.Specifically, The glacier melting amount GM that judgement is determined by the first hydrological simulation model_t' the reality surveyed with satellite or ground observation website Whether the relationship between glacier melting amount meets the default condition of convergence.

When the output of above-mentioned CREST model is to convergent glacier melting amount GM_t' between the practical glacier melting amount of actual measurement Relationship when being unsatisfactory for the preset condition of convergence, then the value of the ice-melt parameter in above-mentioned CREST model is adjusted, thus based on adjusting CREST model after whole is calculated new to convergent glacier melting amount, and further judgement is new melts to convergent glacier Whether the relationship between water and the practical glacier melting equivalent of actual measurement meets the preset condition of convergence；If new to convergent ice Relationship between river melt water amount and the practical glacier melting amount of actual measurement still not satisfies the preset condition of convergence, then continues in adjustment The value for stating ice-melt parameter in CREST model, until the new reality to convergent glacier melting amount and actual measurement obtained again Border glacier melting amount reaches the preset condition of convergence, then stops adjusting, and joins the value of ice-melt parameter adjusted as ice-melt Several actual values.Optionally, the initial set value of the ice-melt parameter in above-mentioned CREST model is replaced with to the reality of the ice-melt parameter Obtained model (i.e. above-mentioned preset hydrological simulation model) after actual value, for calculating the diameter flow data of cryosphere.

Optionally, the above-mentioned preset condition of convergence can be so that the output of CREST model to convergent glacier melting amount and Objective function between satellite or the practical glacier melting amount of ground observation website actual measurement reaches the preset condition of convergence.It is optional , above-mentioned objective function may include: that assorted efficiency function NSE, relative error Bias and index of correlation CC；And it is above-mentioned pre- If the condition of convergence can be with are as follows: so that the value of above-mentioned objective function reaches ideal value or ideal convergence range.Wherein: receiving assorted effect The ideal value of rate function NSE is 1, and the ideal value of relative error function Bias is 0, and the ideal value of correlation of indices function CC is 1.Tool Body, above-mentioned ideal convergence range can be a range near the ideal value of each objective function, and the setting of the range can be with It is selected as needed, such as: the ideal convergence range of assorted efficiency function NSE received can be set between 0.8~1, relatively The ideal convergence range of error function Bias can be set between 0 ± 0.5, and the ideal convergence range of correlation of indices function CC can To be set as between 1 ± 0.3.Further, when the data precision for needing to calculate is high, ideal convergence range can be arranged It is smaller, to meet the measuring and calculating demand of high accuracy data；When the data precision for needing to calculate is low, ideal convergence range can be set That sets is larger, and the efficiency of data measuring and calculating is improved in the case where meeting data precision.

It should be noted that the OBSi in above-mentioned objective function is the actual number of satellite or ground observation website actual observation According to SI Mi is the result that preset CREST simulation model calculates output.In the present embodiment, OBSi is CREST simulation model Be calculated to convergent glacier melting amount, SIMi is the practical glacier Rong Shui of satellite or ground observation website actual observation Amount.

The convergence process of above-mentioned model in order to facilitate understanding, it is following to cite a plain example: assuming that for the first time according to pre- If CREST model and known model training data obtain for the first time to convergent glacier measure GM_t' value be C, computer sets It is standby that S IMi in above three objective function is replaced with into C, OBSi is then replaced with into satellite or the actual measurement of ground observation website Practical glacier melting amount the value of three objective functions is calculated, while judging whether the value of three objective functions reaches respectively To respective ideal value or ideal convergence range.If do not reached, the ice-melt parameter in above-mentioned CREST model is adjusted Value, continues through CREST model adjusted and known model training data obtain second to convergent glacier melting amount GM_t', it is assumed that the value is D, and the S IMi in above three objective function is replaced with D, then replaces OBSi by computer equipment The practical glacier melting amount surveyed for satellite or ground observation website is judged simultaneously with the value of three objective functions is calculated Whether the value of three objective functions respectively reaches respective ideal value or ideal convergence range；If reaching, by above-mentioned adjustment Actual value of the value of ice-melt parameter afterwards as ice-melt parameter.

In the present embodiment, by constantly adjusting the ice-melt parameter in CREST model, so that the ice-melt of simulation model output The practical ice-melt parameter of parameter actual measurement meets the preset condition of convergence, so that the ice-melt parameter in CREST model is more accurate, from And quantify ice-melt and the again influence of the factors such as freeze water diameter flow data into the measuring and calculating of diameter flow data, so that cryosphere The measuring and calculating of diameter flow data is more accurate, so as to preferably instruct people's lives and production.

Fig. 8 is a kind of structural schematic diagram of the acquisition device of diameter flow data provided in an embodiment of the present invention.As shown in figure 8, The device includes: the first determining module 11, the second determining module 12, first processing module 13, Second processing module 14 and third Processing module 15.

Specifically, the first determining module 11, for according to region to be measured driving data and preset hydrological simulation model, Determine accumulated snow melt water amount, glacier melting amount and the snow cover maximum water-holding capacity at current time；The driving data in the region to be measured It include: the air themperature of the precipitation in the region to be measured in preset time period, the surface temperature in region to be measured and region to be measured；Institute State the parameter in preset hydrological simulation model comprising can be used in characterizing the melt water characteristic of cryosphere；

Second determining module 12, for by the solid precipitation quantity of the precipitation in the region to be measured and the region to be measured Difference, the liquid precipitation amount being determined as in the preset time period；

First processing module 13, for the freezing water quantity again according to the water-holding capacity of previous moment, the snow cover of previous moment And the water holding variable quantity of previous moment, obtain the water-holding capacity at current time；

Second processing module 14, for according to the accumulated snow melt water amount, the glacier melting amount, the liquid precipitation amount, The snow cover maximum water-holding capacity and the water-holding capacity at the current time obtain the water holding variable quantity at current time；

Third processing module 15, for by the difference between ideal meltwater runoff and the water holding variable quantity at the current time Value, is determined as the practical diameter flow data at current time.

The acquisition device of diameter flow data provided in an embodiment of the present invention can execute above method embodiment, realize former Reason is similar with technical effect, and details are not described herein.

Fig. 9 is the structural schematic diagram of the acquisition device of another diameter flow data provided in an embodiment of the present invention.In above-mentioned Fig. 8 On the basis of illustrated embodiment, as shown in figure 9, above-mentioned first determining module 11, can specifically include: the first determination unit 111, Second determination unit 112, third determination unit 113 and the 4th determination unit 114.

Specifically, the first determination unit 111, for the driving data and preset hydrological simulation according to the region to be measured Snow melt parameter in model determines the water equivalent of snow at current time and the possibility amount of snowmelt amount in the region to be measured；The snow melt ginseng Number includes environmental parameter for characterizing precipitation or snow melting and for characterizing the snow melt under unit time unit temperature The parameter of amount；

Second determination unit 112, between the possibility amount of snowmelt amount and the water equivalent of snow by the region to be measured most Small value is determined as the accumulated snow melt water amount；

Third determination unit 113, for being melted according to the possibility amount of snowmelt amount in the region to be measured, the accumulated snow in the region to be measured Ice-melt parameter in water, the water equivalent of snow and the hydrological simulation model, determines the glacier melting amount；The ice-melt Parameter includes the environmental parameter for characterizing Melting Glacierss and the ginseng for characterizing the ice storage under unit time unit temperature Number；

4th determination unit 114, for being determined as the product for the product of the water equivalent of snow and the preset water holding factor Snow deposit maximum water-holding capacity；Wherein, the liquid water that the water holding factor is held by unit water equivalent of snow.

The acquisition device of diameter flow data provided in an embodiment of the present invention can execute above method embodiment, realize former Reason is similar with technical effect, and details are not described herein.

Being used to characterize precipitation or the environmental parameter of snow melting described in one of the embodiments, includes: snow melt threshold It is worth temperature, snowfall threshold temperature and rainfall threshold temperature；

The parameter for characterizing the amount of snowmelt amount under unit time unit temperature includes: that the maximum snow melt of summer solstice is subsisted The minimum snow melt degree-day factor of the factor and winter solstice；

It is described for characterize the environmental parameter of Melting Glacierss to include: ice-melt threshold temperature；

The parameter for characterizing the ice storage under unit time unit temperature includes: ice-melt degree-day factor.

Figure 10 is the structural schematic diagram of the acquisition device of another diameter flow data provided in an embodiment of the present invention.In above-mentioned reality On the basis of applying example, as shown in Figure 10, above-mentioned Second processing module 14 be can specifically include: first processing units 141, second Processing unit 142 and third processing unit 143.

Specifically, first processing units 141, for determining the accumulated snow melt water amount, the glacier melting amount and the liquid State precipitation and value；

The second processing unit 142, for determine the snow cover maximum water-holding capacity and the current time water-holding capacity it Between difference；

Third processing unit 143, for determining that the minimum value between value and the difference is the current time Water holding variable quantity.

The acquisition device of diameter flow data provided in an embodiment of the present invention can execute above method embodiment, realize former Reason is similar with technical effect, and details are not described herein.

Figure 11 is the structural schematic diagram of the acquisition device of another diameter flow data provided in an embodiment of the present invention.In above-mentioned reality On the basis of applying example, as shown in figure 11, above-mentioned first processing module 13 be can specifically include: fourth processing unit the 131, the 5th Processing unit 132 and the 6th processing unit 133.

Specifically, fourth processing unit 131, for according to formulaIt determines The possibility of previous moment freezing water quantity W again_t-1,ref；Wherein, the RFC is preset water holding freeze-out fraction again, the T_i,mltFor Ice-melt threshold temperature in the hydrological simulation model, the d_i,mltFor the ice-melt degree-day factor in the hydrological simulation model, The Δ t is to calculate duration；

5th processing unit 132, for by the possibility of previous moment freezing water quantity W again_t-1,refWith it is described previous when Minimum value between the water-holding capacity at quarter is determined as the freezing water quantity again of the snow cover of previous moment；

6th processing unit 133, for will be between the water-holding capacity of the previous moment and the water holding variable quantity of previous moment And value, subtract obtained difference after the freezing water quantity again of the snow cover of previous moment, be determined as the current time Water-holding capacity.

The acquisition device of diameter flow data provided in an embodiment of the present invention can execute above method embodiment, realize former Reason is similar with technical effect, and details are not described herein.

Figure 12 is the structural schematic diagram of the acquisition device of another diameter flow data provided in an embodiment of the present invention.In above-mentioned reality On the basis of applying example, as shown in figure 12, above-mentioned apparatus can also include: third determining module 16 and the 4th determining module 17.

Specifically, third determining module 16, for first determining module according to the driving data in region to be measured and Preset hydrological simulation model, before the accumulated snow melt water amount, glacier melting amount and the snow cover maximum water-holding capacity that determine current time, According to known model training data and preset CREST model, the actual value of snow melt parameter and the reality of ice-melt parameter are determined Value；Wherein, the initial set value of the initial set value in the CREST model including snow melt parameter and ice-melt parameter；

4th determining module 17, it is described for replacing with the initial set value of the snow melt parameter in the CREST model The actual value of snow melt parameter, and the initial set value of the ice-melt parameter in the CREST model is replaced with into the ice-melt and is joined Several actual values obtains the hydrological simulation model；

Wherein, the model training data include: the gross precipitation in trained region, the trained region surface temperature, The air themperature in the trained region, the practical water equivalent of snow surveyed by satellite, the practical glacier melting surveyed by satellite Amount.

The acquisition device of diameter flow data provided in an embodiment of the present invention can execute above method embodiment, realize former Reason is similar with technical effect, and details are not described herein.

Figure 13 is the structural schematic diagram of the acquisition device of another diameter flow data provided in an embodiment of the present invention.Figure 13 be On the basis of Figure 12 embodiment, as shown in figure 13, above-mentioned third determining module 16 be can specifically include: the 5th determination unit 161, the 6th determination unit 162, the first judging unit 163 and the first adjustment unit 164.

Specifically, the 5th determination unit 161, for according to formula Determine the possibility amount of snowmelt amount S in the trained region_t-1,mlt'；Wherein, snow melt degree-day factorDoy is number of days, T_s,mlt' it is snow melt threshold value Temperature, d_s,mlt6For the maximum snow melt degree-day factor of summer solstice, d_s,mlt12For the minimum snow melt degree-day factor of winter solstice, T_ls' for institute State the surface temperature in trained region, T_airmx' be the trained region air themperature, Δ t' be calculate duration；

6th determination unit 162, for according to formula S M_t-1'=min { S_t-1,mlt',SWE_t-1' determine the training center The practical amount of snowmelt amount SM in domain_t-1'；Wherein, the SWE_t-1' initial value be setting value；According to formulaCalculate solid precipitation quantity P_s,'；Wherein, T_s' it is the training center The snowfall threshold temperature in domain, T_r' be the trained region rainfall threshold temperature；

First judging unit 163, for according to formula S WE_t'=SWE_t-1'-SM_t-1'+P_s' be calculated to convergent snow Water equivalent SWE_t', and judge described to convergent water equivalent of snow SWE_t' and by satellite survey practical water equivalent of snow between Whether relationship meets the default condition of convergence；

The first adjustment unit 164, it is described to convergent water equivalent of snow SWE for being determined in first judging unit_t' with When the relationship between practical water equivalent of snow surveyed by satellite is unsatisfactory for the condition of convergence, the snow melt parameter is adjusted Value, obtain it is new to convergent water equivalent of snow, until described new to convergent water equivalent of snow and the reality that passes through satellite and survey Until relationship between water equivalent of snow meets the default condition of convergence；And the value of snow melt parameter adjusted is joined as the snow melt Several actual values.

The acquisition device of diameter flow data provided in an embodiment of the present invention can execute above method embodiment, realize former Reason is similar with technical effect, and details are not described herein.

Figure 14 is the structural schematic diagram of the acquisition device of another diameter flow data provided in an embodiment of the present invention.Figure 14 be On the basis of Figure 13 embodiment, as shown in figure 14, above-mentioned 4th determining module 17 be can specifically include: computing unit 171, Two judging units 172 and second adjustment unit 173.

Specifically, computing unit 171, for according to formula

The trained region is calculated To convergent glacier melting amount GM_t'；Wherein, described, SM_t'、S_t,mlt' it is to utilize the first hydrological simulation model and the training What the model training data in region were calculated, the first hydrological simulation model is to join the snow melt in the CREST model Several initial set values replaces with the model that the actual value of the snow melt parameter obtains later；

Second judgment unit 172, it is described to convergent glacier melting amount GM for judging_t' with pass through satellite survey reality Whether the relationship between the glacier melting amount of border meets the default condition of convergence；

Second adjustment unit 173, for judging the glacier melting amount in the trained region in the second judgment unit GM_t' and the relationship that passes through between the practical glacier melting amount that satellite is surveyed when being unsatisfactory for the default condition of convergence, adjust the ice-melt The value of parameter, obtains new to convergent glacier melting amount, new to convergent glacier melting amount and passes through satellite until described Until relationship between the practical glacier melting amount of actual measurement meets the default condition of convergence；And the value of ice-melt parameter adjusted is made For the actual value of the ice-melt parameter.

The acquisition device of diameter flow data provided in an embodiment of the present invention can execute above method embodiment, realize former Reason is similar with technical effect, and details are not described herein.

In one embodiment of the invention, a kind of computer equipment is also provided, including memory, processor and be stored in storage On device and the computer program that can run on a processor, the computer program are realized above-mentioned when being executed by the processor The acquisition methods of diameter flow data involved in embodiment of the method.

Computer equipment provided in an embodiment of the present invention, can be used for executing above method embodiment, realization principle and Technical effect is similar, and details are not described herein.

In one embodiment of the present of invention, a kind of readable storage medium storing program for executing is additionally provided, is stored in the readable storage medium storing program for executing There is computer program, when the computer program is run on computers, so that computer executes in above method embodiment The acquisition methods for the diameter flow data being related to.

Readable storage medium storing program for executing provided in an embodiment of the present invention can be used for executing above method embodiment, realization principle Similar with technical effect, details are not described herein.

Above-described embodiment all or part of can be realized by software, hardware, firmware or any combination thereof.When It when being realized using software program, can realize in the form of a computer program product in whole or in part, the computer journey Sequence product may include one or more computer instruction.Load and execute on computers the computer program instructions When, it is generated in whole or in part according to process or function described in the embodiment of the present invention.Wherein, the computer can be General computer, special purpose computer, computer network or other programmable devices；The computer instruction can store In computer readable storage medium, or from a computer readable storage medium to another computer readable storage medium pass It is defeated.For example, the computer instruction can pass through wired (example from a web-site, computer, server or data center Such as coaxial cable, optical digital subscriber lines DSL) or wireless (such as infrared, wireless, microwave) mode to another website station Point, computer, server or data center are transmitted；The computer readable storage medium can be computer and can deposit Any usable medium taken either includes the data storage such as one or more usable medium integrated server, data center Equipment.The usable medium can be magnetic medium, such as floppy disk, hard disk, tape, optical medium (such as DVD) or semiconductor Medium (such as solid-state hard disk SSD) etc..

It should be understood by those skilled in the art that, embodiments herein can provide for method, apparatus, computer equipment or The computer program products such as readable storage medium storing program for executing.Therefore, the application can be used complete hardware embodiment, complete software embodiment, Or the form of embodiment combining software and hardware aspects.Moreover, it wherein includes meter that the application, which can be used in one or more, Computer-usable storage medium (including but not limited to magnetic disk storage, CD-ROM, the optical memory of calculation machine usable program code Deng) on the form of computer program product implemented.

The application is referring to method, the process of equipment (system) and computer program product according to the embodiment of the present application Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.

These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates, Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or The function of being specified in multiple boxes.

These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one The step of function of being specified in a box or multiple boxes.

Each technical characteristic of embodiment described above can be combined arbitrarily, for simplicity of description, not to above-mentioned reality It applies all possible combination of each technical characteristic in example to be all described, as long as however, the combination of these technical characteristics is not deposited In contradiction, all should be considered as described in this specification.

In several embodiments provided herein, it should be understood that disclosed device, device and method are not having It has more than in scope of the present application, can realize in other way.For example, embodiments described above is only to illustrate Property, for example, the division of the module or unit, only a kind of logical function partition can have other in actual implementation Division mode, such as multiple units or components can be combined or can be integrated into another system or some features can neglect Slightly, it or does not execute.Wherein the unit as illustrated by the separation member may or may not be physically separated, and make It may or may not be physical unit for the component that unit is shown, it can it is in one place, or can also divide On cloth to multiple network units.Some or all of the modules therein can be selected to realize the present embodiment according to the actual needs The purpose of scheme.

The embodiments described above only express several embodiments of the present invention, and the description thereof is more specific and detailed, but simultaneously Limitations on the scope of the patent of the present invention therefore cannot be interpreted as.It should be pointed out that for those of ordinary skill in the art For, without departing from the inventive concept of the premise, various modifications and improvements can be made, these belong to guarantor of the invention Protect range.Therefore, the scope of protection of the patent of the invention shall be subject to the appended claims.

Claims (11)

1. a kind of acquisition methods of diameter flow data characterized by comprising

According to the driving data in region to be measured and preset hydrological simulation model, accumulated snow melt water amount, the glacier at current time are determined Melt water amount and snow cover maximum water-holding capacity；The driving data in the region to be measured includes: the region to be measured in preset time period The air themperature of gross precipitation, the surface temperature in region to be measured and region to be measured；Include in the preset hydrological simulation model For characterizing the parameter of the melt water characteristic of cryosphere；

By the difference of the gross precipitation in the region to be measured and the solid precipitation quantity in the region to be measured, it is determined as described default Liquid precipitation amount in period；

Changed according to the freezing water quantity again of the snow cover of the water-holding capacity of previous moment, previous moment and the water holding of previous moment Amount, obtains the water-holding capacity at current time；

According to the accumulated snow melt water amount, the glacier melting amount, the liquid precipitation amount, the snow cover maximum water-holding capacity and The water-holding capacity at the current time obtains the water holding variable quantity at current time；

By the difference between ideal meltwater runoff and the water holding variable quantity at the current time, it is determined as the practical diameter at current time Flow data.

2. the method according to claim 1, wherein the driving data according to region to be measured and preset water Literary simulation model determines accumulated snow melt water amount, glacier melting amount and the snow cover maximum water-holding capacity at current time, comprising:

According to the snow melt parameter in the driving data in the region to be measured and preset hydrological simulation model, current time is determined The possibility amount of snowmelt amount of water equivalent of snow and the region to be measured；The snow melt parameter includes for characterizing precipitation or snow melting Environmental parameter and parameter for characterizing the amount of snowmelt amount under unit time unit temperature；

By the minimum value between the possibility amount of snowmelt amount and the water equivalent of snow in the region to be measured, it is determined as the accumulated snow melt water Amount；

According to the possibility amount of snowmelt amount in the region to be measured, accumulated snow melt water amount, the water equivalent of snow and the institute in the region to be measured The ice-melt parameter in hydrological simulation model is stated, determines the glacier melting amount；The ice-melt parameter includes melting for characterizing glacier The environmental parameter of change and parameter for characterizing the ice storage under unit time unit temperature；

By the product of the water equivalent of snow and the preset water holding factor, it is determined as the snow cover maximum water-holding capacity；Wherein, described The liquid water that the water holding factor is held by unit water equivalent of snow.

3. according to the method described in claim 2, it is characterized in that, the environment ginseng for characterizing precipitation or snow melting Number includes: snow melt threshold temperature, snowfall threshold temperature and rainfall threshold temperature；

The parameter for characterizing the amount of snowmelt amount under unit time unit temperature includes: the maximum snow melt degree-day factor of summer solstice With the minimum snow melt degree-day factor of winter solstice；

It is described for characterize the environmental parameter of Melting Glacierss to include: ice-melt threshold temperature；

The parameter for characterizing the ice storage under unit time unit temperature includes: ice-melt degree-day factor.

4. the method according to claim 1, wherein described according to the accumulated snow melt water amount, the glacier melting Amount, the liquid precipitation amount, the snow cover maximum water-holding capacity and the water-holding capacity at the current time, obtain current time Water holding variable quantity, comprising:

Determine the accumulated snow melt water amount, the glacier melting amount and the liquid precipitation amount and value；

Determine the difference between the snow cover maximum water-holding capacity and the water-holding capacity at the current time；

Determine that the minimum value between value and the difference is the water holding variable quantity at the current time.

5. the method according to claim 1, wherein the water-holding capacity according to previous moment, previous moment The freezing water quantity again of snow cover and the water holding variable quantity of previous moment, obtain the water-holding capacity at current time, comprising:

According to formulaDetermine the possibility of previous moment freezing water quantity W again_t-1,ref； Wherein, the RFC is preset water holding freeze-out fraction again, the T_i,mltFor the ice-melt threshold value temperature in the hydrological simulation model Degree, the d_i,mltFor the ice-melt degree-day factor in the hydrological simulation model, the Δ t is to calculate duration；

By the possibility of previous moment freezing water quantity W again_t-1,refMinimum value between the water-holding capacity of the previous moment, really It is set to the freezing water quantity again of the snow cover of previous moment；

By between the water-holding capacity of the previous moment and the water holding variable quantity of previous moment and value, the accumulated snow of previous moment is subtracted Obtained difference after the freezing water quantity again of layer, is determined as the water-holding capacity at the current time.

6. according to the described in any item methods of claim 3-5, which is characterized in that the driving data according to region to be measured and Preset hydrological simulation model, before the accumulated snow melt water amount, glacier melting amount and the snow cover maximum water-holding capacity that determine current time, The method also includes:

According to known model training data and preset CREST model, the actual value and ice-melt parameter of snow melt parameter are determined Actual value；Wherein, the initial set value of the initial set value in the CREST model including snow melt parameter and ice-melt parameter；

The initial set value of snow melt parameter in the CREST model is replaced with to the actual value of the snow melt parameter, and will The initial set value of ice-melt parameter in the CREST model replaces with the actual value of the ice-melt parameter, obtains the hydrology Simulation model；

Wherein, the model training data include: the surface temperature, described of the gross precipitation in trained region, the trained region Training the air themperature in region, the practical water equivalent of snow of actual measurement, actual measurement practical glacier melting amount.

7. according to the method described in claim 6, it is characterized in that, model training data according to known to and preset CREST model determines the actual value of snow melt parameter, comprising:

According to formulaDetermine the possibility amount of snowmelt amount S in the trained region_t-1,mlt'； Wherein, snow melt degree-day factorDoy is day Number, T_s,mlt' it is snow melt threshold temperature, d_s,mlt6For the maximum snow melt degree-day factor of summer solstice, d_s,mlt12Melt for the minimum of winter solstice Avenge degree-day factor, T_ls' be the trained region surface temperature, T_airmx' be the trained region air themperature, Δ t' is Calculate duration；

According to formula S M_t-1'=min { S_t-1,mlt',SWE_t-1' determine the practical amount of snowmelt amount SM in the trained region_t-1'；Wherein, The SWE_t-1' initial value be setting value；

According to formulaCalculate solid precipitation quantity P_s'；Wherein, T_s' For the snowfall threshold temperature in the trained region, T_r' be the trained region rainfall threshold temperature, P_tot' it is the training center The gross precipitation in domain；

According to formula S WE_t'=SWE_t-1'-SM_t-1'+P_s' be calculated to convergent water equivalent of snow SWE_t', and judge it is described to Convergent water equivalent of snow SWE_t' actual measurement practical water equivalent of snow between relationship whether meet the default condition of convergence；

If it is not, then adjust the value of the snow melt parameter, obtain it is new to convergent water equivalent of snow, until described new to convergent Until relationship between water equivalent of snow and the practical water equivalent of snow of actual measurement meets the default condition of convergence；

Using the value of snow melt parameter adjusted as the actual value of the snow melt parameter.

8. the method according to the description of claim 7 is characterized in that model training data according to known to and preset CREST model determines the actual value of ice-melt parameter, comprising:

According to formulaThe trained region is calculated To convergent glacier melting amount GM_t'；Wherein, the SWE_t'、SM_t'、S_t,mlt' it is to utilize the first hydrological simulation model and the instruction What the model training data in white silk region were calculated, the first hydrological simulation model is by the snow melt in the CREST model The initial set value of parameter replaces with the model that the actual value of the snow melt parameter obtains later；

Judgement is described to convergent glacier melting amount GM_t' actual measurement practical glacier melting amount between relationship whether meet it is default The condition of convergence；

If it is not, then adjust the value of the ice-melt parameter, obtain it is new to convergent glacier melting amount, until described new wait restrain Glacier melting amount and actual measurement practical glacier melting amount between relationship meet the default condition of convergence until；

Using the value of ice-melt parameter adjusted as the actual value of the ice-melt parameter.

9. a kind of acquisition device of diameter flow data characterized by comprising

First determining module, for according to region to be measured driving data and preset hydrological simulation model, determine current time Accumulated snow melt water amount, glacier melting amount and snow cover maximum water-holding capacity；The driving data in the region to be measured includes: preset time The precipitation in region to be measured, the air themperature of the surface temperature in region to be measured and region to be measured in section；The preset hydrology It include the parameter that can be used in characterizing the melt water characteristic of cryosphere in simulation model；

Second determining module, for by the difference of the precipitation in the region to be measured and the solid precipitation quantity in the region to be measured Value, the liquid precipitation amount being determined as in the preset time period；

First processing module, for according to the water-holding capacity of previous moment, the freezing water quantity again of the snow cover of previous moment and preceding The water holding variable quantity at one moment, obtains the water-holding capacity at current time；

Second processing module, for according to the accumulated snow melt water amount, the glacier melting amount, the liquid precipitation amount, the product Snow deposit maximum water-holding capacity and the water-holding capacity at the current time obtain the water holding variable quantity at current time；

Third processing module, for determining the difference between ideal meltwater runoff and the water holding variable quantity at the current time For the practical diameter flow data at current time.

10. a kind of computer equipment including memory, processor and stores the meter that can be run on a memory and on a processor Calculation machine program, which is characterized in that any one of claim 1-8 is realized when the processor is for executing the computer program The step of method.

11. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer program The step of any one of claim 1-8 method is realized when being executed by processor.