AU2019100631A4

AU2019100631A4 - Self-correcting multi-model numerical rainfall ensemble forecasting method

Info

Publication number: AU2019100631A4
Application number: AU2019100631A
Authority: AU
Inventors: Chuanzhe LI; Jia Liu; Jiyang Tian; Yang Wang; Fuliang Yu
Original assignee: China Institute of Water Resources and Hydropower Research
Current assignee: China Institute of Water Resources and Hydropower Research
Priority date: 2019-06-12
Filing date: 2019-06-12
Publication date: 2019-07-18
Anticipated expiration: 2027-06-12

Abstract

Abstract Abstract The present application relates to a self-correcting multi-model numerical rainfall ensemble forecasting method, comprising the following steps: step 1, selecting various numerical weather prediction models; step 2, simulating forecasting and outputting rainfall data for every T hours; step 3, evaluating rainfall forecast results; step 4, determining a forecast weight coefficient of each model; and step 5, releasing a forecast result. The present application can more objectively evaluate the rainfall forecast results of all numerical weather prediction models on the basis of existing multi-model ensemble rainfall forecast, so that the final ensemble rainfall forecast result does not depend too much on man-made decisions and thus the released rainfall forecast result is more objective. Drawings First numerica l Second numerical M-th numeric a l weather prediction weather prediction weather prediction model model model Actuallymeasured 6hr.rainfall 6hr.rainfall 6hr.rainfall value of 6hr. rainfall forecast forecast forecast Probabilty ofdetection(POD) Maximum error (ME) Frequently bias Root mean Square erro (RMSB) False alarm ratiFAR Mean bi as error (MBE) Critical success index (CSI) a a deva on (S D) R ainfall foreca st Qualitative evaluati on:] Quantitativeevaluati on evaluation H Normalization Comprehensive model score [Determination of forecast weight c oefficient of each mode I [Output of forecast result Figure 1 Memory '32 Output device Processor Figure 2

Description

SELF-CORRECTING MULTI-MODEL NUMERICAL RAINFALL ENSEMBLE FORECASTING METHOD

TECHNICAL FIELD [0001] The present application relates to a self-correcting multi-model numerical rainfall ensemble forecasting method which is mainly used in multi-model ensemble rainfall forecast carried out by a meteorological department, a water conservancy department and other departments.

BACKGROUND [0002] For a long time, rainfall forecast, as an important part of numerical weather prediction, has been widely concerned by related scholars. As a formation process and occurrence of rainfall are affected by multiple aspects, such as large-scale atmospheric circulation, ocean current, land and sea location, topography, underlying surface and human activities, there are uncertainties in the spatial and temporal distribution of rainfall, which makes the rainfall forecast more difficult than other meteorological factors. In general, a rainfall forecast time step is 6h for numerical weather prediction model; the shorter the time step is, the more difficult the forecast is, resulting in lower forecast accuracy; while if the time step is too long, the forecast accuracy will decrease gradually due to the forecast period influence of the model itself. In recent years, with continuous development and improvement of numerical weather prediction models as well as progress of computer technology, ensemble rainfall forecast, having the advantages of reduced uncertainty of single numerical weather prediction model and improved rainfall forecast reliability, has become a main means of the rainfall forecast in meteorological and water conservancy departments and other departments.

[0003] At present, there are numerous numerical weather prediction models, such as the US WRF model, the UKMO model, the Canada MC2 model, the JRSM model, the Chinese GRAPES model and other models, whose application is wider. In a process of implementing the present application, the inventors at least found following problems in related arts: accuracies of the same type of rainfall forecast by various models at the same time are different; so, for the ensemble rainfall forecast, how to select and evaluate a numerical weather prediction model is the most important problem. However, if a man-made decision that one numerical weather prediction model is no longer selected due to its poor forecast accuracy of once or several times of rainfall is approved, the model forecast uncertainties will be increased. Although man-made decisions are

2019100631 12 Jun 2019

Description important to determine an ensemble rainfall forecast result, it depends on experience greatly, possibly causing wrong judgments or choices.

SUMMARY [0004] The present application designs a self-correcting multi-model numerical rainfall ensemble forecasting method, and solves the technical problems of different accuracies of the same type of rainfall forecast by various models at the same time and how to select and evaluate the optimal numerical rainfall ensemble forecast.

[0005] To solve one or more technical problems in the prior art, the present application provides a self-correcting multi-model numerical rainfall ensemble forecasting method.

[0006] According to a first aspect of the embodiments of the present application, provided is a self-correcting multi-model numerical rainfall ensemble forecasting method, comprising the following steps: step 1, selecting various numerical weather prediction models;

step 2, simulating forecasting, and outputting rainfall data for every T hours;

step 3, evaluating rainfall forecast results;

step 4, determining a forecast weight coefficient of each model; and step 5, releasing a forecast result.

[0007] Further, in the step 2, a rainfall output time step T is set as 6 hours.

[0008] Further, in the step 3, after the forecast results of 6 hours of rainfall through the selected various numerical models are output, based on an actually measured result of the rainfall, the rainfall forecast results are comprehensively evaluated in qualitative and quantitative manners respectively while considering time and space as well as point rainfall and areal rainfall, and the comprehensive evaluation results are scored.

[0009] Further, in the step 3, the comprehensive evaluation comprises qualitative evaluation in which first a forecasted rainfall value and an actually measured value are compared and assessed in a graded manner, and then classification evaluation indices are established according to an assessment result.

[0010] Specifically, when the classification indexes are used in spatial dimension evaluation, firstly forecasted values and actually measured values of a specific observation time step i at different observation locations are compared to acquire classification variables NA_t, NB_t and NC_t in a rainfall grade table, then the classification indices at all the time steps are statistically averaged according to equations (1)-(4), and

2019100631 12 Jun 2019

Description finally a classification evaluation result in the spatial dimension is obtained, wherein spatial scale evaluation indexes comprise:

POO,(probability of detection)— — V------ (1);

N NA_: + NC_:

A zvz?

FAR, (false alarm ratio)= — V------’--- (3); and

N NA. + NB_j

[0011] In the above equations, NA_t, NBi and NC_t respectively indicate whether the forecasted values and the actually measured values at the different observation locations within an i-th 6h observation time period are in corresponding rainfall grades in the rainfall grade table, N is the number of observation time periods, and the areal rainfall is a rainfall mean value at all rainfall stations.

[0012] For temporal dimension, firstly forecasted values and actually measured values of a specific observation location j at different observation time points are compared, the classification variables in the rainfall grade table are counted, then classification indices at all observation locations in a study area are statistically averaged according to equations (5)-(8), and finally a classification evaluation result in temporal dimension is obtained, wherein temporal dimension evaluation indices comprise:

M^NA_j+NC_j

.... 1 NA. + NB .

FS/tffrequently bias index)— — > ---------^L (6);

M fa NA.+NC.

[0013] NAj , NBj and NCj respectively indicate whether the forecasted values and the actually measured values of the observation location j at the different observation time points are in corresponding rainfall grades in the rainfall grade table, and M is the number of the observation locations.

[0014] Further, the rainfall grade table is shown hereinafter:

Description

2019100631 12 Jun 2019

Rainfall grades	Light rain	Moderate rain	Heavy rain	Torrential rain	Downpour	Excessively heavy rain
6hr rainfall (mm)	0.1-2.5	2.6-6	6.1-12	12.1-25	25.1-60	>60

[0015] The above variables NA_t, NB. and NC, are calculated as follows: for spatial dimension evaluation, within a specific observation time step i, if both a rainfall forecasted value and a rainfall observation value at an observation location are within any one of the above six rainfall grades, NA_f is marked as 1; if the rainfall observation value is within any one of the above six rainfall grades, but the rainfall forecasted value is not in any one of the above six rainfall grades, and is not equal to 0, NB_t is marked as 1; if the rainfall observation value is within any one of the above six rainfall grades, and the rainfall forecast value is 0mm, that is, the numerical weather model does not acquire rainfall information, JVC,. is marked as 1.

[0016] The above variables NA, , NB, and NC. are calculated as follows: for temporal dimension evaluation, within a specific observation location j, if both a rainfall forecasted value and a rainfall observation value at the observation location are within any one of the above six rainfall grades, NAj is marked as 1; if the rainfall observation value is within any one of the above six rainfall grades, but the rainfall forecasted value is not in any one of the above six rainfall grades, and is not equal to 0, NBj is marked as 1; if the rainfall observation value is within any one of the above six rainfall grades, and the rainfall forecast value is 0 mm, that is, the numerical weather model does not acquire rainfall information, NCj is marked as f.

[0017] Further, in the step 3, the comprehensive evaluation further comprises quantitative evaluation which adopts four quantitative evaluation indexes in error analysis. For temporal dimension evaluation, 7? and (9, respectively represent a forecast value and an actually measured value of the mean rainfall in the study area at an observation time point i, which are shown in equations (9)-(12):

Description

2019100631 12 Jun 2019

Λ-/Ε,(maximum error)—max|P,-0,|(9);

I ] λ· ~

RMSE/root mean square error)— J— ΣΖ ( / ^- β ) (1 θΧ V-Y i=l ^N

AffiE/(mean bias error)— — ^(/-//(ll); and

N ,=i l~l n 7

SD_t(standard deviation)= — Σ(/ “β “MBE) (12).

y ΛΜ ,-=1 [0018] where i represents one of different observation time periods, N is the number of the observation time periods, and MBE is the value of the mean deviation MBE_t.

[0019] For spatial dimension evaluation, P. and Qj respectively represent a forecasted value and an actually measured value of accumulated rainfall in the whole observation time period at the specific spatial location j, which are shown in equations (13)-(16):

AfEvtiTiiiximuin error) —max|P/-0/|(13);

RMSEs (root mean square error)—

SD_s (standard deviation)—

[0020] where j represents one of different observation locations, M is the number of observation locations, and MBE is the value of the mean deviation MBEs.

[0021] Further, in step 3, the above 8 classification evaluation indexes and 8 quantitative evaluation indices are used to establish an index system for rainfall forecast of each numerical weather prediction model, and thus a rainfall forecast result of each numerical weather prediction model is scored based on the above 16 evaluation indices.

[0022] Assuming that m numerical weather prediction models are adopted, each evaluation index is normalized. For example, with respect to indices of k numerical weather prediction models, PODtk: SPODtk=(PODtk-PODtmin)/(PODtmax-PODtmin)(17), wherein k is 1, ..., or m, and is the number of numerical weather prediction models; and PODtmax and PODtmin respectively represent the maximum and the minimum of m PODt corresponding to m numerical weather prediction models. The normalization of other evaluation indexes is calculated according to the above equation (17).

2019100631 12 Jun 2019

Description [0023] After normalization, each numerical weather prediction model is scored, a comprehensive score is represented by S, and Sk represents the comprehensive score of a k-th numerical weather prediction model, wherein ^x^.we_sj. ^x f*RMSE,_t ^X 8RMSE,t ^X MB E,t ^X ^MBE* ^X ^SD,t ^X ^SDst (18) [0024] Further, in the step 4, a coefficient, obtained by using a rainfall forecast score of any one of numerical weather prediction models to divide the sum of comprehensive scores of the all models, is used as a rainfall forecast weight coefficient of the numerical weather prediction model. As a solution for the next ensemble rainfall forecast, the weight coefficient a_k is calculated as follows: a_k = S_k (S) +...+ S_m ), wherein k is 1, ..., or m, and is the number of the numerical weather prediction models, and S_k represents the comprehensive score of the k-th numerical weather prediction model.

[0025] Further, in the step 4, after a previous rainfall forecast weight coefficient a_k is obtained and the next rainfall is completed, the previous rainfall forecast weight coefficient a_k is corrected based on a forecast value and an actually measured value of the next rainfall to be used as a solution for subsequent rainfall forecast.

[0026] Further, in the step 5, the forecast result of ensemble rainfall forecast is obtained by each model forecast result multiplied by its forecast weight coefficient:

P? = P_n x a, + P_n x a, + · + P_Pm x a_m [0027] P_Pm represents forecast rainfall of an m-th numerical weather prediction model at an observation location within a time period, and Gt_m represents a weight coefficient of the m-th numerical weather prediction model at the observation location within the time period.

[0028] According to a second aspect of the embodiments of the present application, provided is a non-transitory computer storage medium storing computer-executable instructions, which cause a computer to perform any above-mentioned self-correcting multi-model numerical rainfall ensemble forecasting method.

2019100631 12 Jun 2019

Description [0029] According to a third aspect of the embodiments of the present application, provided is a computer program product comprising computer programs stored in a non-transitory computer-readable storage medium and comprising program instructions, which cause a computer to perform the any above-mentioned self-correcting multi-model numerical rainfall ensemble forecasting method when the program instructions are executed by the computer.

[0030] According to a fourth aspect of the embodiments of the present application, provided is an electronic equipment, comprising at least one processor and a memory configured to store instructions, and when executed by the at least one processor, causing the at least one processor to perform any above-mentioned self-correcting multi-model numerical rainfall ensemble forecasting method.

[0031] The self-correcting multi-model numerical rainfall ensemble forecasting method provided by the embodiments of the present application has the following advantageous effects:

1) the embodiments of the present application can more objectively evaluate rainfall forecast results of all numerical weather prediction models on the basis of existing multi-model ensemble rainfall forecast, so that a final result of the ensemble rainfall forecast does not depend too much on man-made decisions and thus the released rainfall forecast result is more objective; and

2) the embodiments of the present application provide a comprehensive evaluation index system for analysis of the rainfall forecast results qualitatively and quantitatively in consideration of time and space as well as point rainfall and areal rainfall, each model is scored through a corresponding evaluation result by the index system, then the coefficient, obtained by using the rainfall forecast score of the numerical weather prediction model to divide the sum of the scores of the all models, is used as the rainfall forecast weight coefficient of the numerical weather prediction model, and finally, a next rainfall forecast result is determined.

BRIEF DESCRIPTION OF THE DRAWINGS [0032] One or more embodiments are illustrated by corresponding accompanying drawings which are not intended to limit the scope of the present invention. Components in the drawings with the same reference numbers in the accompanying drawings represent similar elements and there is no scale limitation in the drawings otherwise particularly represented.

[0033] Fig. 1 is a schematic flowchart of a self-correcting multi-model numerical rainfall ensemble forecasting method according to the first embodiment of the present application; and

2019100631 12 Jun 2019

Description [0034] Fig. 2 is a schematic diagram showing a hardware structure of an equipment for a self-correcting multi-model numerical rainfall ensemble forecasting method provided by the fourth embodiment of the present application.

DETAILED DESCRIPTION [0035] In order to illustrate purposes, technical solutions and advantages of the present application more clearly, the technical solutions will be clearly and completely described through implementations with reference to the accompanying drawings in the embodiments of the present application hereinafter. Obviously, the described embodiments below are merely for illustrating some embodiments of the present application.

Embodiment I [0036] The present application will be further illustrated with reference to Fig.las follows.

[0037] The technical solution adopted in the present application is a self-correcting multi-model numerical rainfall ensemble forecasting method based on a scoring method. The method mainly comprises two parts, namely, firstly, running of each numerical weather prediction model, and secondly, evaluation on each numerical weather prediction model running result and ensemble of the forecast results so that more objective ensemble rainfall forecast can be achieved and uncertainty of model forecast can be reduced. The method can be implemented by a self-correcting numerical rainfall ensemble forecasting device. For example, the device may be a weather forecasting platform, and can be configured in a smart terminal for use. The method is implemented by the steps as follows.

[0038] In step 1, numerical weather prediction models are selected. In this step, a plurality of currently popular numerical weather prediction models are selected to be installed on the same weather forecast platform.

[0039] In step 2, a forecast is simulated. In this step, on the weather forecast platform, an initial time, a boundary condition, a physical parameterization solution, terrain data and the like of each model are set and processed respectively. And, running is carried out according to an operation method of each model respectively to perform rainfall forecast and rainfall forecast results are output based on a time step of 6hours.

[0040] In step 3, the rainfall forecast results are evaluated, wherein, based on an actually measured result of the rainfall, the rainfall forecast results of all models are comprehensively evaluated in a qualitative and quantitative manner respectively considering time and space as well as point rainfall and areal rainfall, and

Description

2019100631 12 Jun 2019 the evaluated results are scored.

[0041] For qualitative evaluation, firstly, rainfall forecasted values and actually measured values are compared and assessed in a graded manner, in which grading standards are shown in below table 1. According to the assessment results, classification evaluation indices are established. Spatial dimension evaluation indices are shown in equations (1)-(4), and temporal dimension evaluation indexes are shown in equations (5)-(8).

Table 1 Rainfall Grade Table

Rainfall grades	Light rain	Moderate rain	Heavy rain	Torrential rain	Downpour	Excessively heavy rain
6hr. rainfall (mm)	0.1-2.5	2.6-6	6.1-12	12.1-25	25.1-60	>60

A NAPO£>_s(probability of detcction)= — V-----—--- (1);

N NA; + NC;

] * yVA + NB

FBI]frequently bias index)=— V----------(2);

N NA; + NC_t /ARJfalse alarm ratio)=—

Σ ,=1

NB;

NA; + NB;

(3); and

CSTfcritical success index)= —

NA;

NA; + NB; + NC.

(4).

[0042] NA;, NB; and NC_t respectively indicate whether the forecasted values and the actually measured values at different observation locations within an i-th 6h observation time step are in corresponding rainfall grades in the table 1, N is the number of observation time periods (6hr), and the areal rainfall is a rainfall average at all rainfall stations.

^M ΝΛ

PODt (probability of detection)— — V--------- (5);

M ΝΑ_} + NC.

..... 1 ™NA;+NB;

frequently bias index)— — V—^L----(6);

M NA;+NC;

NB.

FAR_t(false alarm ratio)= — V------I---(7); and

M NA. + NB_:

AL, NA.

CSI/(critical success index)=—y ---------------- (8).

M^NA_j+NB_j+NC_J [0043] NA., NB. and NC respectively indicate whether the forecast values and the actually measured values of the observation location j at the different observation time points are in corresponding rainfall

2019100631 12 Jun 2019

Description grades in the table 1, and M is the number of the observation locations.

[0044] The variables NA, NB and NC are calculated as follows: for example, during the spatial dimension evaluation, within a specific observation time step i, if both a rainfall forecasted value and a rainfall observation value at an observation location are in the range of 0.1-2.5 mm (light rain), NA_{ is marked as 1; if the rainfall observation value is in the range of 0.1-2.5 mm (light rain), but the rainfall forecasted value is not in the range, and is not equal to 0, ΝΒ_ί is marked as 1; if the rainfall observation value is in the range of 0.1-2.5 mm (light rain), and the rainfall forecasted value is 0 mm; that is, the numerical weather prediction model does not acquire rainfall information, NC_{ is marked as 1.

[0045] Assuming that there are six observation locations totally, if there is one observation location whose forecasted value and the observation value are in the range of 0.1-2.5 mm (light rain), NA, = 1; if there are two observation locations whose forecasted values and the observation values are in the range of 0.1-2.5 mm (light rain), NB_t = 2; if there are three observation locations whose forecasted values and the observation values are in the range of 0.1-2.5 mm (light rain), NCI =3. Therefore, within an i-th time period, PODs=l/(3+l)=l/4, FBIs=(l+2)/(l+3)=3/4, FARs=2/(l+2)=2/3andCSIs=l/(l+2+3)=l/6, which are statistical results within the i-th time period, and then all index values within N time periods are calculated to obtain a mean. For temporal dimension evaluation, the calculation method is the same as that of the spatial dimension evaluation.

[0046] Quantitative evaluation adopts four common quantitative evaluation indices in the error analysis. For temporal dimension evaluation, and O, respectively represent a forecast value and an actually measured value of a mean rainfall in the study area at the observation time i, which are shown in equations (9)-(12).

A/iXmaximum error)=max|/’_J-Oi|(9);

RMSEt(roo\. mean square error)=

SDt (standard deviation)=

Description

2019100631 12 Jun 2019 [0047] For spatial dimension evaluation, P. and O respectively represent a forecast value and an actually measured value of accumulated rainfall in the whole observation time period at a specific spatial location j, which are shown in equations (13)-(16):

ME_S (maximum error) =max|/^,/-O,|

RMSE_s (root mean square error)=

AffiA(mean

SA (standard deviation)=

-« , ,, <¹⁴);

[0048] The above 8 classification indices and 8 quantitative indices are combined to establish an index system for rainfall forecast of each numerical weather prediction model, and thus a rainfall forecast result of each numerical weather prediction model is scored based on the above 16 indices. Assuming that m numerical weather prediction models are used, each index is normalized.

[0049] With respect to indexes, POD_tk :

SPODtk=(PODtk-PODin,i_n)/(PODt_max-PODt_mtn) (17) wherein k is 1,..., or m.

[0050] After normalization, each numerical weather prediction model is scored, a comprehensive score is represented by S, and S_k represents the comprehensive score of a k-th numerical weather prediction model, wherein

'•St

Λ '^sk ^lk

RMSE.

(18) [0051] In step 4, a forecast weight coefficient of each model is determined, wherein a coefficient, obtained by using a rainfall forecast score of any one of numerical weather prediction models to divide the sum of the scores of all models, is used as a rainfall forecast weight coefficient of the numerical weather prediction model. As a solution for the next ensemble rainfall forecast, if and only if the actually measured rainfall of this rainfall is greater than 0.1mm, a rainfall forecast weight coefficient can be adjusted, and each weight coefficient is calculated as follows:

Description

2019100631 12 Jun 2019 ^k A/^1 + + ί',η) QQ) [0052] The larger the weight coefficient is, the greater the a_k is, which indicates that a forecast value of a k-th numerical weather prediction model is closer to its observation value.

[0053] In step 5, a forecast result is released, wherein the next rainfall is forecasted and the forecast result is released according to the determined forecast weight coefficients of the all models in this rainfall, and the forecast result is obtained by each model forecast result multiplied by its forecast weight coefficient:

A - Al ^{X ft}l ⁺ A2 ^{X a}2 + · + Pp_m ^X in which P_Pm represents forecast rainfall of an m-th numerical weather prediction model at an observation location within a time period.

Embodiment II [0054] It should be understood by those skilled in the art that, all or part of the steps of the above method provided by the embodiments may be implemented through programs that give instructions to respective hardware. The above programs may be stored in a computer-readable storage medium. During program implementation, the steps of the above method provided by the embodiments are implemented. The above storage medium may be an ROM, an RAM, a magnetic disk, an optical disk or other media capable of storing program codes.

[0055] The embodiments of the present application provide a non-transitory computer storage medium storing computer-executable instructions, which cause a computer to perform the self-correcting multi-model numerical rainfall ensemble forecasting method provided by any of the above embodiments.

[0056] As an embodiment, the non-transitory computer storage medium of the present application stores computer-executable instructions, and the computer-executable instructions is set as follows:

step 1, selecting various numerical weather prediction models;

step 2, simulating forecasting, and outputting rainfall data for every T hours;

step 3, evaluating rainfall forecast results;

[0057] As a non-transitory computer-readable storage medium, it can be used to store non-transitory software programs, non-transitory computer-executable programs and modules, and corresponding program

2019100631 12 Jun 2019

Description instructions/modules used in the self-correcting multi-model numerical rainfall ensemble forecasting method provided by the embodiments of the present application. When the one or more modules stored in the non-transitory computer-readable storage medium are executed by the processor, the self-correcting multi-model numerical rainfall ensemble forecasting method provided by any of the above embodiments is performed.

[0058] The non-transitory computer-readable storage medium may include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required by at least one function; the storage data area may store data and the like created during the operation of a self-correcting multi-model numerical rainfall ensemble forecasting device. In addition, the non-transitory computer-readable storage medium may include a high-speed random access memory and may also include a non-transitory memory. For example, the memory comprises at least one disk storage device, a flash memory device or other non-transitory solid state memory. In some embodiments, the non-transitory computer-readable storage medium may optionally include memories remotely configured with respect to the processor, and the memories may be connected to the self-correcting multi-model numerical rainfall ensemble forecasting device via networks. Examples of the networks include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

Embodiment III [0059] The embodiments of the present application provide a computer program product comprising a computer program stored on a non-transitory computer-readable storage medium, when program instructions included in the computer program are executed by a computer, the computer can perform any above-mentioned self-correcting multi-model numerical rainfall ensemble forecasting method.

Embodiment IV [0060] Fig. 2 is a schematic diagram showing a hardware structure of electronic equipment used for implementing the self-correcting multi-model numerical rainfall ensemble forecasting method and provided by the fourth embodiment of the present application. As shown in Fig.2, the equipment comprises:

[0061] One or more processors 210 and a memory 220, wherein in Fig. 2, one processor 210 is provided. [0062] The equipment for implementing the self-correcting multi-model numerical rainfall ensemble forecasting method may further include an input device 230 and an output device 240.

2019100631 12 Jun 2019

Description [0063] The processor 210, the memory 220, the input device 230, and the output device 240 may be connected via a bus or other means. As shown in Fig. 2, they are connected through a bus.

[0064] The input device 230 may receive input digital or character information and generate a key signal input related to user setting and function control of the self-correcting multi-model numerical rainfall ensemble forecasting device. The output device 240 may include display equipment, such as a display screen.

[0065] When the one or more modules stored in the memory 220 are executed by the one or processors 220, the self-correcting multi-model numerical rainfall ensemble forecasting method provided by any of the above embodiments is performed.

[0066] The above-described product can implement the method provided by the embodiments of the present invention, and has corresponding function modules for implementing the method and beneficial effects. Technical details which are not described in detail in the embodiments can refer to the method provided by the embodiments of the present application.

[0067] The electronic equipment of the embodiments of the present application exists in a variety of forms, and comprises, but is not limited to:

1) mobile communication equipment: this type of equipment is characterized by having mobile communication capabilities and mainly aims to provide voice and data communication, and these terminals include: smart phones (such as iPhone), multimedia phones, functional phones, low-end phones and the like;

2) ultra-mobile personal computer equipment: this type of equipment belongs to the field of personal computers, has computing and processing functions, and generally, also has a mobile Internet feature, and these terminals include: PDA, MID, UMPC and others, such as an iPad;

3) a server: this is equipment used for providing computing services, the server is composed of a processor, a hard disk, a memory, a system bus and the like, an architecture of the server is similar to that of a general computer, however, the server needs to provide highly reliable services, so it has high requirements on processing capacity, stability, reliability, security, scalability, manageability and other aspects; and

4) other electronic devices with data processing functions.

[0068] The above device embodiments are illustrative only. The units described as separate members may be or may not be physically separated. The members described as units may be or may not be physical units, may be located at the same place or may be distributed in multiple network units. The objectives of the

2019100631 12 Jun 2019

Description solutions of this application may be realized by selecting some or all of the modules according to the actual needs.

[0069] Through the description of the above embodiments, those skilled in the art can understand clearly that the all embodiments may be implemented through software and an indispensable universal hardware platform, of course, also be implemented through hardware. Based on such understanding, essentially, the above technical solutions or parts contributing to the related arts can be embodied in the form of a software product, the computer software product may be stored in a computer-readable storage medium, such as ROM/RAM, a magnetic disk, an optical disk, or the like, which includes a plurality of instructions to make computer equipment (which may be a personal computer, a server, network equipment, or the like) to perform all or part of the steps of the method of all embodiments of the application.

[0070] At last, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present application and not intended to limit them. Although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art can understand that the technical solutions described in the foregoing embodiments can be modified or some of the technical features thereof can be equivalently replaced, and these modifications or substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A self-correcting multi-model numerical rainfall ensemble forecasting method, comprising the following steps:

step 1, selecting various numerical weather prediction forecast models;

step 2, simulating forecasting and outputting rainfall data for every T hours;

step 3, evaluating rainfall forecast results;

step 4, determining a forecast weight coefficient of each model; and step 5, releasing a forecast result, wherein the step 3 further comprises: after outputting the forecast results of T hours of rainfall through the selected various numerical models, based on an actually measured result of the rainfall, comprehensively evaluating the rainfall forecast results in qualitative and quantitative manners respectively while considering time and space as well as point rainfall and areal rainfall, and scoring the comprehensive evaluation results;

in the step 3, the comprehensive evaluation comprises the qualitative evaluation in which first a forecasted rainfall value and an actually measured value are compared and assessed in a graded manner, and then classification evaluation indices are established according to an assessment result; specifically, when the classification evaluation indices are used in spatial dimension evaluation, firstly forecasted values and actually measured values of a specific observation time step i at different observation locations are compared to acquire classification variables NAi, NBi and NCi in a rainfall grade table, then the classification indices at all the time steps are statistically averaged according to equations (1)-(4), and finally a classification evaluation result in the spatial dimension is obtained; spatial scale evaluation indexes comprise:

1 * jyj

POD_S (probability of detection)= — / ’-----—---- (1);

N fc NA. + NC_i

FBI, (frequently bias index)=—Τ’---^!------^L (2);

N fc + PC,

FAR, (false alarm ratio)= — V------ (3); and

N fc NA. + NB,

1 * Ay

CS1, (critical success index)=— 2 -----------^!--------(4);

N fc Nf + NB_t + NQ in the above equations, NA,, NB, and NC_t respectively indicate whether the forecasted values and the actually measured values at the different observation locations within an i-th T observation time period are in corresponding rainfall grades in the rainfall grade table, N is the number of observation time periods,

2019100631 12 Jun 2019

Claims and the areal rainfall is a rainfall mean value at all rainfall stations;

when the classification evaluation indices are used in temporal dimension evaluation, firstly forecasted values and actually measured values of a specific observation location j at different observation time points are compared, the classification variables NAj , NB, and NC. in the rainfall grade table are counted, then classification indices at all observation locations in a study area are statistically averaged according to equations (5)-(8), and finally a classification evaluation result in temporal dimension is obtained; the temporal dimension evaluation indices comprise:

1 m na.

POD/probability of detection)=—V-------- (5);

M jzt NA. + NC_j „ ..... 1 “ΝΑ,+ΝΒ, .

FBI/frequently bias index)=— > ---------- (6);

M NA_j + NC.

1 ^w NB

FAR/false alarm ratio)= — V--------- (7); and

M ΝΑ. + NB ' J ,1 ] ^M NA

CSHcritical success index)= — V--------¹------ (8),

Mj^NA.+NB. + NC.

NAj , NBj and NC respectively indicate whether the forecasted values and the actually measured values of the observation location j at the different observation time points are in corresponding rainfall grades in the rainfall grade table, and M is the number of the observation locations;

the rainfall grade table is as follows:

for spatial scale dimension evaluation, the above variables NA_t, NBj and NCj are calculated as follows: within an observation time step i, if both a rainfall forecasted value and a rainfall observation value at an observation location are within any one of the above six rainfall grades, NAj is marked as 1; if the rainfall observation value is within any one of the above six rainfall grades but the rainfall forecasted value is not in any one of the above six rainfall grades and is not equal to 0, NBj is marked as 1; if the rainfall observation value is within any one of the above six rainfall grades, and the rainfall forecasted value is 0 mm; that is, the numerical weather prediction model does not acquire rainfall information, NC_t is marked as 1;

Claims

2019100631 12 Jun 2019 for temporal dimension evaluation, the above variables ΝΑ -, NB, and NC are calculated as follows: within a specific observation location j, if both a rainfall forecasted value and a rainfall observation value at the observation location are within any one of the above six rainfall grades, NA- is marked as 1; if the rainfall observation value is within any one of the above six rainfall grades but the rainfall forecasted value is not in any one of the above six rainfall grades and is not equal to 0, NB, is marked as 1; if the rainfall observation value is within any one of the above six rainfall grades, and the rainfall forecast value is 0 mm; that is, the numerical weather prediction model does not acquire rainfall information, NC j is marked as 1;

in the step 3, the comprehensive evaluation also includes quantitative evaluation adopting four quantitative evaluation indexes in error analysis;

for temporal dimension evaluation, P_t and O_t respectively represent a forecasted value and an actually measured value of the mean rainfall in the study area at the observation time point i, which are shown in equations (9)-(12):

(11); and

Substandard deviation)^ j—!— ^(Z’ ~O_t~MBE)² (12), where i represents one of different observation time steps, N is the number of the observation time steps, and MBE is the value of the mean deviation MBE_t;

for spatial dimension evaluation, Pj and Q- respectively represent a forecasted value and an actually measured value of accumulated rainfall in the whole observation time period at the specific spatial location j, which are shown in equations (13)-(16):

ME_S (maximum error) =max|F/-0_;| (13);

RMSEftoaX mean square etTor}= (15); and

2019100631 12 Jun 2019

Claims where j represents one of different observation locations, M is the number of the observation locations, and MBE is the value of the mean deviation MBE_s;

in the step 3, the above 8 classification evaluation indexes and 8 quantitative evaluation indices are used to establish an index system for rainfall forecast of each numerical weather prediction model, and thus a rainfall forecast result of each numerical weather prediction model is scored based on the above 16 evaluation indices;

assuming that m numerical weather prediction models are adopted, each evaluation index is normalized; for example, with respect to indices of k numerical weather prediction models, POD_tk :

SpoDtkAPODtk-POD^KPOD^-POD^C) (17) wherein k is 1,..., or m, and is the number of the numerical weather prediction models, POD_tmm and

POD_rmin respectively represent the maximum and the minimum of m POD_t corresponding to the m numerical weather predictionmodels, and the normalization of other evaluation indices is calculated according to the above equation (17); and after normalization, each numerical weather prediction model is scored, a comprehensive score is represented by S, and Sk represents the comprehensive score of a k-th numerical weather prediction model, which is shown in the followings:

/ ^{x x X} ^^EAk_t * ^{X X} / ^{X X X} SuBE^ ^{X X} j (18).

2. The self-correcting multi-model numerical rainfall ensemble forecasting method according to claim 1, wherein in the step 2, the rainfall output time step T is set as 6 hours.

3. The self-correcting multi-model numerical rainfall ensemble forecasting method according to claim 1, wherein in the step 4, a coefficient, obtained by using a rainfall forecast score of each of numerical weather prediction models to divide the sum of comprehensive scores of all models, is used as a rainfall forecast weight coefficient of each of the numerical weather prediction model; and as a solution for next ensemble rainfall forecast, the weight coefficient ak is calculated as follows:

^ak “ ^Sj(^St ⁺ ⁺

Claims

2019100631 12 Jun 2019 wherein k is 1,..., or m, and is the number of the numerical weather prediction models, and S_k represents the comprehensive score of the k-th numerical weather prediction model.

4. The self-correcting multi-model numerical rainfall ensemble forecasting method according to claim 3, wherein in the step 4, after a previous rainfall forecast weight coefficient ak is obtained and when the next rainfall is completed, the previous rainfall forecast weight coefficient ak is corrected based on a forecast value and an actually measured value of the next rainfall to be used as a solution for subsequent rainfall forecast.

5. The self-correcting multi-model numerical rainfall ensemble forecasting method according to claim 1, wherein in the step 5,the forecast result of ensemble rainfall forecast is obtained by each model forecast result multiplied by its forecast weight coefficient:

^PP = ^x «1 ^{+ p}P2 ^x «2 + + P_Pai * a_m (20) wherein P represents forecast rainfall of the m-th numerical weather prediction model at an observation location within a time period, and am represents a weight coefficient of the m-th numerical weather prediction model at the observation location within the time period.

6. Electronic equipment, comprising:

at least one processor, and a memory in communication with the at least one processor, wherein the memory stores instructions executable by the at least one processor, and when executed by the at least one processor, causing the at least one processor to perform the method according to claim 1.

7. A non-transitory computer-readable storage medium storing computer instructions, which cause a computer to perform the method according to claim 1 when the computer instructions are executed by the computer.