Rain enhancement operation effect evaluation method and device for intelligently predicting natural rainfall

Abstract

The invention provides a rainfall enhancement operation effect evaluation method and device for intelligently predicting natural rainfall, and the method comprises the steps: determining an operation influence region and a comparison region based on the range of artificial rainfall enhancement operation; the contrast area is not influenced by a catalyst for artificial precipitation increasing operation; respectively carrying out rainfall extrapolation prediction on the contrast area and the work influence area through a rainfall prediction model to obtain rainfall prediction data of the contrast area and rainfall prediction data of the work influence area; respectively obtaining real rainfall data of the comparison area and real rainfall data of the operation influence area; based on the rainfall prediction data of the comparison area and the real rainfall data of the comparison area, obtaining a rainfall prediction error of the rainfall prediction model; and on the basis of the rainfall prediction data of the operation influence area, the real rainfall data of the operation influence area and the rainfall prediction error of the rainfall prediction model, obtaining the absolute rainfall increment of the artificial rainfall enhancement operation. And the accuracy of effect evaluation is improved.

Description

Technical Field

[0001] The present application relates to the technical field of rainmaking operation effect evaluation, and in particular to a rainmaking operation effect evaluation method, device, electronic device and storage medium for intelligently predicting natural rainfall. Background Art

[0002] The vitality of weather modification science depends on the effect of its application and its role in disaster reduction and mitigation. Therefore, effect evaluation is one of the important links that cannot be avoided in artificial rainfall operations and occupies an important position in the entire operation project. Scientifically and objectively measuring the actual effectiveness of artificial rainfall operations can help us understand the response mechanism of clouds under different operating conditions, analyze the changes in rainfall and its relationship with operating parameters, and thus optimize key factors such as operation timing, catalyst type and catalyst amount. This is crucial to improving cloud and precipitation theory and improving the efficiency and service level of weather modification operations. However, due to the uncertainty of the evaluation object and the huge natural variation of clouds and precipitation, it is very difficult to effectively test the rainfall effect in these variations. Therefore, effect testing is still a global problem so far.

[0003] Traditional effect verification methods (physical test, statistical test, numerical model test) have many problems and difficulties, which are mainly manifested in: limitations of observation means, imperfections of statistical methods, the need to improve the reliability of numerical models, incomplete understanding of the laws of clouds and precipitation and the mechanism of artificial weather modification, etc., which affect the accuracy of effect evaluation. Summary of the invention

[0004] The present application aims to solve one of the technical problems in the related art at least to some extent.

[0005] To this end, the first purpose of this application is to propose a method for evaluating the effectiveness of rainmaking operations by intelligently predicting natural rainfall, so as to improve the accuracy of the effect evaluation.

[0006] The second purpose of this application is to provide a rain-making operation effect evaluation device that can intelligently predict natural rainfall.

[0007] The third objective of the present application is to provide an electronic device.

[0008] A fourth objective of the present application is to provide a computer-readable storage medium.

[0009] A fifth object of the present application is to provide a computer program product.

[0010] To achieve the above-mentioned purpose, the first embodiment of the present application proposes a method for evaluating the effect of rainmaking operations by intelligently predicting natural rainfall, comprising:

[0011] Based on the scope of the artificial rainmaking operation, determining the operation impact area and the comparison area; the comparison area is not affected by the catalyst of the artificial rainmaking operation;

[0012] Using a rainfall prediction model, rainfall extrapolation prediction is performed on the comparison area and the operation impact area respectively, to obtain rainfall prediction data for the comparison area and rainfall prediction data for the operation impact area;

[0013] Respectively obtaining the actual rainfall data of the comparison area and the actual rainfall data of the operation impact area;

[0014] Based on the rainfall prediction data of the comparison area and the actual rainfall data of the comparison area, obtaining the rainfall prediction error of the rainfall prediction model;

[0015] Based on the rainfall prediction data of the operation influence area, the actual rainfall data of the operation influence area and the rainfall prediction error of the rainfall prediction model, the absolute rainfall increase of the artificial rainfall increase operation is obtained, and the absolute rainfall increase of the artificial rainfall increase operation is used as an evaluation indicator of the rainfall increase operation effect.

[0016] In some implementations, before the rainfall prediction model is used to perform rainfall extrapolation prediction on the comparison area and the operation impact area respectively, the method includes:

[0017] Determine the rainfall type within the scope of the artificial rainfall enhancement operation; wherein the rainfall type includes stratiform cloud rainfall and convective cloud rainfall, and the rainfall type is divided according to the maximum rainfall intensity;

[0018] The rainfall extrapolation model corresponding to the rainfall type is selected as the rainfall prediction model.

[0019] In some implementations, the training method of the rainfall prediction model includes:

[0020] A sample set is constructed based on minute grid precipitation data, wherein each sample in the sample set includes 2L consecutive hours of rainfall data, wherein the K-minute grid rainfall data of the first L hours of each sample is used as input parameters, and the K-minute grid rainfall data of the last L hours is used as label data;

[0021] The PredRNN++ model is trained by the sample set to obtain a rainfall extrapolation model.

[0022] In some implementations, the step of performing rainfall extrapolation prediction on the comparison area and the operation impact area respectively through a rainfall prediction model to obtain rainfall prediction data for the comparison area and rainfall prediction data for the operation impact area comprises:

[0023] Based on the time of the artificial rainfall enhancement operation, obtaining first rainfall data collected by a rainfall station in the comparison area in a first time period before the artificial rainfall enhancement operation;

[0024] Based on the first rainfall data, performing rainfall extrapolation prediction on the comparison area through the rainfall prediction model to obtain rainfall prediction data for the comparison area in a second time period;

[0025] Based on the time of the artificial rainfall enhancement operation, obtaining second rainfall data collected by a rainfall station in the operation impact area in the first time period before the artificial rainfall enhancement operation;

[0026] Based on the second rainfall data, rainfall extrapolation prediction is performed on the operation impact area through the rainfall prediction model to obtain rainfall prediction data for the operation impact area in the second time period.

[0027] In some implementations, respectively obtaining the real rainfall data of the comparison area and the real rainfall data of the operation impact area includes:

[0028] Obtaining real rainfall data collected by n rainfall stations in the comparison area;

[0029] Obtain actual rainfall data collected by m rainfall stations in the operation impact area.

[0030] In some implementations, obtaining the rainfall prediction error of the rainfall prediction model based on the rainfall prediction data of the comparison area and the actual rainfall data of the comparison area includes:

[0031] Calculate the difference between the actual rainfall data collected by n rainfall stations in the comparison area and the rainfall forecast data in the comparison area to obtain n differences;

[0032] The average value of the n differences is taken as the rainfall prediction error of the rainfall prediction model.

[0033] In some implementations, obtaining the absolute rainfall increase of the artificial rainfall increase operation based on the rainfall prediction data of the operation influence area, the actual rainfall data of the operation influence area and the rainfall prediction error of the rainfall prediction model comprises:

[0034] Based on the rainfall forecast data of the operation influence area and the actual rainfall data collected by the m rainfall stations in the operation influence area, the average absolute rainfall increase of the m rainfall stations in the operation influence area is obtained;

[0035] Calculate the difference between the average absolute rainfall increase of m rainfall stations in the operation influence area and the rainfall prediction error of the rainfall prediction model to obtain the absolute rainfall increase of a single rainfall station in the operation influence area;

[0036] The absolute rainfall increase of the artificial rainfall increase operation is obtained by calculating the product of the absolute rainfall increase of a single rainfall station in the operation influence area and the area of ​​the operation influence area.

[0037] To achieve the above-mentioned purpose, the second embodiment of the present application proposes a rainmaking operation effect evaluation device for intelligently predicting natural rainfall, comprising:

[0038] A region division module is used to determine an operation impact area and a comparison area based on the scope of the artificial rainfall enhancement operation; the comparison area is not affected by the catalyst of the artificial rainfall enhancement operation;

[0039] A rainfall prediction module, used to perform rainfall extrapolation prediction on the comparison area and the operation impact area respectively through a rainfall prediction model to obtain rainfall prediction data for the comparison area and rainfall prediction data for the operation impact area;

[0040] A rainfall acquisition module, used to respectively acquire the real rainfall data of the comparison area and the real rainfall data of the operation impact area;

[0041] An error calculation module, used for obtaining a rainfall prediction error of the rainfall prediction model based on the rainfall prediction data of the comparison area and the actual rainfall data of the comparison area;

[0042] An effect evaluation module is used to obtain the absolute rainfall increase of the artificial rainmaking operation based on the rainfall prediction data of the operation influence area, the actual rainfall data of the operation influence area and the rainfall prediction error of the rainfall prediction model, and use the absolute rainfall increase of the artificial rainmaking operation as an evaluation indicator for the effect of the rainmaking operation.

[0043] In some implementations, the rainfall prediction module is further used to:

[0044] Determine the rainfall type within the scope of the artificial rainfall enhancement operation; wherein the rainfall type includes stratiform cloud rainfall and convective cloud rainfall, and the rainfall type is divided according to the maximum rainfall intensity;

[0045] The rainfall extrapolation model corresponding to the rainfall type is selected as the rainfall prediction model.

[0046] In some implementations, the apparatus further includes a model training module, configured to:

[0047] A sample set is constructed based on minute grid precipitation data, wherein each sample in the sample set includes 2L consecutive hours of rainfall data, wherein the K-minute grid rainfall data of the first L hours of each sample is used as input parameters, and the K-minute grid rainfall data of the last L hours is used as label data;

[0048] The PredRNN++ model is trained by the sample set to obtain a rainfall extrapolation model.

[0049] In some implementations, the rainfall prediction module is specifically used to:

[0050] Based on the time of the artificial rainfall enhancement operation, obtaining first rainfall data collected by a rainfall station in the comparison area in a first time period before the artificial rainfall enhancement operation;

[0051] Based on the first rainfall data, performing rainfall extrapolation prediction on the comparison area through the rainfall prediction model to obtain rainfall prediction data for the comparison area in a second time period;

[0052] Based on the time of the artificial rainfall enhancement operation, obtaining second rainfall data collected by a rainfall station in the operation impact area in the first time period before the artificial rainfall enhancement operation;

[0053] Based on the second rainfall data, rainfall extrapolation prediction is performed on the operation impact area through the rainfall prediction model to obtain rainfall prediction data for the operation impact area in the second time period.

[0054] In some implementations, the rainfall acquisition module is specifically used to:

[0055] Obtaining real rainfall data collected by n rainfall stations in the comparison area;

[0056] Obtain actual rainfall data collected by m rainfall stations in the operation impact area.

[0057] In some implementations, the error calculation module is specifically used to:

[0058] Calculate the difference between the actual rainfall data collected by n rainfall stations in the comparison area and the rainfall forecast data in the comparison area to obtain n differences;

[0059] The average value of the n differences is taken as the rainfall prediction error of the rainfall prediction model.

[0060] In some implementations, the effect evaluation module is specifically used to:

[0061] Based on the rainfall forecast data of the operation influence area and the actual rainfall data collected by the m rainfall stations in the operation influence area, the average absolute rainfall increase of the m rainfall stations in the operation influence area is obtained;

[0062] Calculate the difference between the average absolute rainfall increase of m rainfall stations in the operation influence area and the rainfall prediction error of the rainfall prediction model to obtain the absolute rainfall increase of a single rainfall station in the operation influence area;

[0063] The absolute rainfall increase of the artificial rainfall increase operation is obtained by calculating the product of the absolute rainfall increase of a single rainfall station in the operation influence area and the area of ​​the operation influence area.

[0064] To achieve the above-mentioned purpose, the third aspect embodiment of the present application proposes an electronic device, comprising: a processor, and a memory communicatively connected to the processor; the memory stores computer-executable instructions; the processor executes the computer-executable instructions stored in the memory to implement the method described in the first aspect.

[0065] To achieve the above-mentioned purpose, the fourth aspect of the present application proposes a computer-readable storage medium, in which computer-readable storage medium is stored computer-executable instructions, and when the computer-executable instructions are executed by a processor, they are used to implement the method described in the first aspect.

[0066] To achieve the above-mentioned purpose, the fifth aspect of the present application proposes a computer program product, including a computer program, which implements the method described in the first aspect when executed by a processor.

[0067] The method, device, electronic device and storage medium for evaluating the effect of rain-enhancing operations for intelligently predicting natural rainfall provided by the present application, by setting up an operation influence area and a comparison area, the comparison area is used to calculate the rainfall extrapolation error of the rainfall prediction model, and the operation influence area eliminates the rainfall extrapolation error when calculating the rainfall enhancement effect, thereby significantly improving the accuracy of the evaluation of the effect of rain-enhancing operations. In addition, the corresponding rainfall prediction models are trained for different rainfall types, which effectively reduces the sample complexity, improves the model training efficiency, and the accuracy of rainfall extrapolation prediction. In addition, the present invention realizes the accurate extrapolation of 10-minute rainfall grid data for the next two hours through the rainfall prediction model, which is used as the natural rainfall benchmark for evaluating the rainfall enhancement effect. And the PredRNN++ extrapolation model, with its powerful spatiotemporal state transmission capability, not only grasps the development context of precipitation distribution in precipitation prediction, but also accurately captures the dynamic changes of precipitation intensity, laying a solid foundation for the inspection of the effect of rain-enhancing operations. Compared with traditional methods, the rainfall extrapolation forecast provided by the PredRNN++ extrapolation model has higher temporal and spatial resolution, especially in the field of nowcasting, which effectively makes up for the shortcomings of traditional numerical models in time effectiveness and resolution.

[0068] Additional aspects and advantages of the present application will be given in part in the description below, and in part will become apparent from the description below, or will be learned through the practice of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0069] The above and / or additional aspects and advantages of the present application will become apparent and easily understood from the following description of the embodiments in conjunction with the accompanying drawings, in which:

[0070] Figure 1 A schematic diagram of a flow chart of a method for evaluating the effect of rainmaking operations for intelligently predicting natural rainfall provided in an embodiment of the present application;

[0071] Figure 2 Schematic diagram of a flow chart made for the data set provided in this application example;

[0072] Figure 3 A schematic diagram of the structure of the PredRNN++ extrapolation model provided in the embodiment of the present application;

[0073] Figure 4 A block diagram of a rainmaking operation effect evaluation device for intelligently predicting natural rainfall provided in an embodiment of the present application;

[0074] Figure 5 A block diagram of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0075] Embodiments of the present application are described in detail below, and examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to be used to explain the present application, and should not be construed as limiting the present application.

[0076] The following describes the rain-enhancement operation effect evaluation method, device and equipment for intelligently predicting natural rainfall according to an embodiment of the present application with reference to the accompanying drawings.

[0077] Figure 1 A flow chart of a method for evaluating the effectiveness of rainmaking operations by intelligently predicting natural rainfall provided in an embodiment of the present application.

[0078] It should be noted that the executor of the method for evaluating the effect of rainmaking operations by intelligently predicting natural rainfall in the embodiment of the present application is the device for evaluating the effect of rainmaking operations by intelligently predicting natural rainfall in the embodiment of the present application, and the device for evaluating the effect of rainmaking operations by intelligently predicting natural rainfall can be configured in an electronic device so that the electronic device can perform the function of evaluating the effect of rainmaking operations by intelligently predicting natural rainfall.

[0079] like Figure 1 As shown, the method for evaluating the effect of rainmaking operation by intelligently predicting natural rainfall includes the following steps:

[0080] Step S101, based on the scope of the artificial rainfall enhancement operation, determine the operation impact area and the comparison area; the comparison area is not affected by the catalyst of the artificial rainfall enhancement operation.

[0081] As a way of implementation, the principles for determining the scope of the operation influence area and the comparison area in actual rainfall enhancement operations are as follows:

[0082] (1) The prevailing wind direction in the comparison area and the operation impact area should be similar.

[0083] It can be understood that if the prevailing wind directions are different, the correlation coefficient between the two areas will be poor and the analysis sensitivity will be low.

[0084] To count the prevailing wind direction at high altitude, the comparison area should be selected on the upwind side or on the side perpendicular to the wind direction (but this does not exclude the possibility of pollution caused by dynamic effects). If the prevailing wind direction of the comparison area and the operation impact area is inconsistent, more than one comparison area should be selected, and different comparison areas should be selected according to different prevailing wind directions.

[0085] (2) The comparison area and the operation impact area should be as close as possible, and the comparison area should not be affected by the catalyst of the rain-enhancement operation.

[0086] Generally, the comparison area and the operation impact area are about 30-50km apart.

[0087] (3) The topography and area of ​​the comparison area are similar to those of the operation impact area.

[0088] (4) The comparison area and the operation impact area are affected by the same weather system and have similar rainfall types. This is reflected in the fact that the rainfall in the comparison area is closely related to the rainfall in the operation impact area. The significance level of the sample correlation coefficient should be above 0.01.

[0089] (5) The comparison area should have a dense network of rainfall stations, just like the operation impact area.

[0090] Usually every 20-50km 2 One, depending on the nature of the precipitation and the degree to which the equipment conditions permit.

[0091] (6) The comparison area and the operation impact area have historical data of the same long time series, so as to conduct data comparison analysis between the comparison area and the operation impact area and train the rainfall prediction model.

[0092] In some embodiments, when determining the operation impact area and the comparison area, the rainfall stations of the operation impact area and the comparison area are determined at the same time.

[0093] For example, in the process of an actual artificial rainfall enhancement operation, a part of the rainfall stations in a city area is determined as the rainfall stations in the operation impact area, with a total of m; the other part of the rainfall stations is determined as the rainfall stations in the comparison area (non-operation impact stations), with a total of n. Based on the rainfall stations in the operation impact area and the rainfall stations in the comparison area, the scope of the operation impact area and the scope of the comparison area are delineated.

[0094] The principles for determining the rainfall stations in the operation impact area during rainfall enhancement operations are as follows:

[0095] The propagation distance of the rainfall stations in the operation impact area along the high-altitude wind direction within a certain period of time is taken as the impact range. For example, the impact time of the catalyst is three hours, so the rainfall stations in the operation impact area are based on the propagation distance of the high-altitude wind (generally 700hpa) for 3 hours as the operation impact area, and the number of rainfall stations M in the operation impact area is counted.

[0096] Step S102, using a rainfall prediction model to perform rainfall extrapolation prediction on the comparison area and the operation impact area respectively, to obtain rainfall prediction data for the comparison area and rainfall prediction data for the operation impact area.

[0097] The rainfall prediction model is first described in detail below.

[0098] In some embodiments, a training method for a rainfall prediction model includes: constructing a sample set based on minute grid precipitation data, each sample in the sample set includes 2L consecutive hours of rainfall data, wherein the K-minute grid precipitation data for the first L hours of each sample is an input parameter, and the K-minute grid precipitation data for the last L hours is a label data; training a PredRNN++ model with the sample set to obtain a rainfall extrapolation model.

[0099] It can be understood that the rainfall prediction model is a rainfall extrapolation model, which is used for rainfall extrapolation prediction.

[0100] Exemplarily, the precipitation data at minute grid points is used as the training sample set for training the rainfall prediction model, with a spatial resolution of 0.01 degrees, and the training sample set is stored in the form of a tensor. According to previous experiments, the rainfall prediction model uses the rainfall data of the past two hours to predict the rainfall data of the next two hours with better results; therefore, L is 2, that is, 2L is 4 hours, K is 10 minutes, and the actual rainfall data at minute grid points every 10 minutes is used as the input data and label data in the training sample set, that is, the continuous rainfall data every 4 hours is used as a sample set, in which the 12 (2 hours, 10 minutes collection frequency, 12 data) rainfall data of the first 2 hours are input data, and the 12 rainfall data of the last 2 hours are output data.

[0101] It should also be noted that since different intervention methods are used in rain enhancement operations according to different rainfall types, the samples participating in model training and model testing are classified according to rainfall types.

[0102] As an implementation method, rainfall types are differentiated according to the maximum rainfall intensity of 5 mm / h, at which the reflectivity is approximately equal to 35 dBZ, which is an important reference indicator for stratiform cloud precipitation and convective cloud precipitation.

[0103] In some embodiments, by distinguishing the maximum rainfall intensity, a convective cloud precipitation sample set and a stratiform cloud precipitation sample set are established respectively, and a training set, a validation set and a test set are allocated in proportion to facilitate classification training and optimization of the model.

[0104] For example, the dataset production process is as follows: Figure 2 As shown, it includes obtaining grid rainfall data, continuously checking the grid rainfall data, decoding and clipping the grid rainfall data after verification, and obtaining a data tensor; at the same time, according to the maximum rainfall intensity, the decoded and clipped grid rainfall data are classified to obtain a layered cloud precipitation data set and a convective cloud precipitation data set, and a sample set is obtained based on the data set.

[0105] The rainfall prediction model of the embodiment of the present invention uses PredRNN++ as the basic training model. The basic unit of PredRNN++ is Causal LSTM, which realizes the step-by-step connection of time state, space state and output results, enhances the spatiotemporal prediction capability, and adds a gradient fast transfer unit (GHU) to prevent the disappearance of gradient. The model structure is as follows: Figure 3 As shown, it is the PredRNN++ extrapolation model.

[0106] During model training, the convective cloud precipitation extrapolation model and the stratiform cloud precipitation extrapolation model are trained separately. Due to the different rainfall intensities of different precipitation types, different normalization parameters can be used for training. The model training uses the Adam optimizer for training, and according to the performance of the trained rainfall prediction model on the validation data set, an automatic stop strategy is formulated to prevent overfitting. Among them, in view of the uneven distribution of rainfall intensity, the loss function of the rainfall prediction model is designed by combining MAE and RMAE to increase the complexity of the loss function, thereby improving the generalization ability of the rainfall prediction model. As shown in the following formula (1), the values ​​of the proportional coefficients α and β in formula (1) are determined by the effect test of the rainfall prediction model, and finally the minute rainfall extrapolation model is obtained, that is, the rainfall prediction model of the present invention.

[0107] loss=α·MAE(y',y)+β·RMSE(y',y) (1)

[0108] In the model verification stage, three indicators, namely critical success index CSI, hit rate POD and false alarm rate FAR, are introduced to verify the rainfall forecast data for the next two hours. The indicator formula is shown in formula (2), where a represents the number of hits, b represents the number of false alarms, and c represents the number of missed alarms.

[0109]

[0110] Exemplarily, the model effect test of the rainfall prediction model is completed using the data set from 2019 to 2022 when no rain enhancement operations were carried out. According to the sample characteristics of convective cloud precipitation and stratiform cloud precipitation, and the impact of rain enhancement operations on rainfall of different magnitudes, the extrapolation effect test of the model is carried out at four rainfall intensities of 1mm, 2mm, 5mm and above 5mm.

[0111] In order to test the impact of the artificial rainmaking operation in the next two hours, the extrapolation effect of the model is tested every 10 minutes, so as to realize the extrapolation prediction of rainfall in different rainfall periods, and use the predicted rainfall as the rainfall of natural clouds (without considering the catalytic effect) in the process of artificial rainmaking operation, and then carry out the next step of rainmaking effect evaluation.

[0112] Therefore, in some embodiments, before performing rainfall extrapolation prediction on the comparison area and the operation impact area respectively through the rainfall prediction model, it includes: determining the rainfall type within the scope of the artificial rainmaking operation; wherein the rainfall types include layered cloud rainfall and convective cloud rainfall, and the rainfall types are divided according to the maximum rainfall intensity; and selecting the rainfall extrapolation model corresponding to the rainfall type as the rainfall prediction model.

[0113] Therefore, by modeling stratiform cloud precipitation and convective cloud precipitation separately during the model training stage, the sample complexity is effectively reduced and the training efficiency is improved.

[0114] After determining the rainfall prediction model, as an implementation method, rainfall extrapolation prediction is performed on the comparison area and the operation influence area respectively through the rainfall prediction model to obtain rainfall prediction data for the comparison area and rainfall prediction data for the operation influence area; including: based on the time of the artificial rainmaking operation, obtaining the first rainfall data collected by the rainfall station in the comparison area in the first time period before the artificial rainmaking operation; based on the first rainfall data, extrapolating the rainfall in the comparison area through the rainfall prediction model to obtain rainfall prediction data for the comparison area in the second time period; based on the time of the artificial rainmaking operation, obtaining the second rainfall data collected by the rainfall station in the operation influence area in the first time period before the artificial rainmaking operation; based on the second rainfall data, extrapolating the rainfall in the operation influence area through the rainfall prediction model to obtain rainfall prediction data for the operation influence area in the second time period.

[0115] For example, according to the actual rain enhancement operation time, through the rainfall prediction model, the 10-minute rainfall data of a single rainfall station in the comparison area within 2 hours before the actual rain enhancement operation starts are selected as the input parameters of the rainfall prediction model, and the rainfall prediction model is used to calculate the 10-minute rainfall forecast data of the rainfall station in the next 2 hours.

[0116] Step S103, respectively obtaining the actual rainfall data of the comparison area and the actual rainfall data of the operation impact area.

[0117] As an implementation method, a method for respectively obtaining the real rainfall data of the comparison area and the real rainfall data of the operation impact area includes: obtaining the real rainfall data collected by n rainfall stations in the comparison area; obtaining the real rainfall data collected by m rainfall stations in the operation impact area.

[0118] Step S104, obtaining a rainfall prediction error of a rainfall prediction model based on the rainfall prediction data of the comparison area and the actual rainfall data of the comparison area.

[0119] As an implementation method, a method for obtaining the rainfall prediction error of a rainfall prediction model based on the rainfall prediction data and the actual rainfall data in the comparison area comprises: respectively calculating the difference between the actual rainfall data collected by n rainfall stations in the comparison area and the rainfall prediction data in the comparison area to obtain n differences; and taking the average value of the n differences as the rainfall prediction error of the rainfall prediction model.

[0120] For example, according to the actual rainfall enhancement operation time, the rainfall prediction model is used to select the 10-minute rainfall data of a single rainfall station in the comparison area within 2 hours before the actual rainfall enhancement operation starts as the input parameter of the rainfall prediction model, and the rainfall prediction model is used to calculate the 10-minute rainfall prediction data of the rainfall station in the next 2 hours; the first 10-minute rainfall prediction data in the next 2 hours is recorded as X' i ; The actual rainfall data of the first 10 minutes of the rainfall station in the comparison area is recorded as X i , then the rainfall estimation error for the first 10 minutes at this rainfall station is X' i -X i Then, the rainfall estimation errors of n rainfall stations in the comparison area are counted, and the average rainfall estimation error of the first 10 minutes is calculated. ave for:

[0121]

[0122] The average value X ave is the 10-minute average rainfall estimation error of multiple rainfall stations in the comparison area.

[0123] Step S105, based on the rainfall prediction data of the operation influence area, the actual rainfall data of the operation influence area and the rainfall prediction error of the rainfall prediction model, the absolute rainfall increase of the artificial rainfall increase operation is obtained, and the absolute rainfall increase of the artificial rainfall increase operation is used as an evaluation indicator of the rainfall increase operation effect.

[0124] As an implementation method, a method for obtaining the absolute rainfall increase of an artificial rainmaking operation based on rainfall forecast data of an operation influence area, actual rainfall data of the operation influence area and rainfall forecast error of a rainfall forecast model; comprising: obtaining the average absolute rainfall increase of the m rainfall stations in the operation influence area based on rainfall forecast data of the operation influence area and actual rainfall data collected by m rainfall stations in the operation influence area; calculating the difference between the average absolute rainfall increase of the m rainfall stations in the operation influence area and the rainfall forecast error of the rainfall forecast model to obtain the absolute rainfall increase of a single rainfall station in the operation influence area; calculating the product of the absolute rainfall increase of a single rainfall station in the operation influence area and the area of ​​the operation influence area to obtain the absolute rainfall increase of the artificial rainmaking operation.

[0125] Taking the above example as an example, according to the actual rainfall enhancement operation time, the 10-minute rainfall data of a single rainfall station (for example, the first rainfall station) in the operation impact area within 2 hours before the actual rainfall enhancement operation is selected as the input parameter of the rainfall prediction model. The rainfall prediction model is used to calculate the 10-minute rainfall prediction data of the rainfall station in the next 2 hours, and the first 10-minute rainfall prediction data in the next 2 hours is recorded as Y' j The actual rainfall data of the first 10 minutes of the rainfall station is recorded as Y j , then the absolute rainfall increase in the first 10 minutes of the rainfall station is Y j -Y' j Then, the absolute rainfall increase of the m rainfall stations in the operation impact area is counted, and the average value Y of the absolute rainfall increase in the first 10 minutes is calculated. ave for:

[0126]

[0127] The average value Yave is the average value of the absolute rainfall increase within 10 minutes at multiple rainfall stations in the operation impact area.

[0128] According to the scope of the operation influence zone, the area G of the operation influence zone is estimated, and the absolute rainfall increase Ra of the 10-minute rainfall enhancement operation in the operation influence zone is:

[0129] R a =(Y ave -X ave )×G(5)

[0130] Relative rainfall rate R of 10-minute rainfall enhancement operationr for:

[0131]

[0132] By analogy, the absolute rainfall increase every 10 minutes in the next 2 hours after the rain-making operation can be obtained, and the absolute rainfall increase of the entire rain-making operation can be accumulated.

[0133] It should be noted that, combined with the rainfall enhancement operation time, the absolute rainfall enhancement of multiple rainfall stations in multiple periods of time can be obtained according to formula (5), thereby obtaining the absolute rainfall enhancement of the entire rainfall enhancement operation; the rainfall enhancement rate obtained according to formula (6) can also be combined with the actual rainfall data collected in the operation impact area to obtain the absolute rainfall enhancement of the entire rainfall enhancement operation.

[0134] The rain-enhancing operation effect evaluation method for intelligently predicting natural rainfall in the embodiment of the present application, by setting up an operation influence area and a comparison area, the comparison area is used to calculate the rainfall extrapolation error of the rainfall prediction model, and the operation influence area eliminates the rainfall extrapolation error when calculating the rainfall enhancement effect, thereby significantly improving the accuracy of the rain-enhancing operation effect evaluation. In addition, the corresponding rainfall prediction models are trained for different rainfall types, which effectively reduces the sample complexity, improves the model training efficiency, and the accuracy of rainfall extrapolation prediction. In addition, the present invention realizes the accurate extrapolation of K-minute rainfall grid data for the next two hours through the rainfall prediction model, which is used as the natural rainfall benchmark for evaluating the rainfall enhancement effect. And the PredRNN++ extrapolation model, with its powerful spatiotemporal state transmission ability, not only grasps the development context of precipitation distribution in precipitation prediction, but also accurately captures the dynamic changes of precipitation intensity, laying a solid foundation for the inspection of rain-enhancing operation effects. Compared with traditional methods, the rainfall extrapolation forecast provided by the PredRNN++ extrapolation model has higher temporal and spatial resolution, especially in the field of nowcasting, which effectively makes up for the shortcomings of traditional numerical models in time effectiveness and resolution.

[0135] In order to implement the above-mentioned embodiment, the present application also proposes a rain-making operation effect evaluation device for intelligently predicting natural rainfall. Figure 4 This is a schematic diagram of the structure of a rain-enhancement operation effect evaluation device for intelligently predicting natural rainfall provided in an embodiment of the present application. Figure 4 As shown, the rain enhancement operation effect evaluation device for intelligently predicting natural rainfall may include: a region division module 401, a rainfall prediction module 402, a rainfall prediction module 403, an error calculation module 404 and an effect evaluation module 405.

[0136] The area division module 401 is used to determine the operation impact area and the comparison area based on the scope of the artificial rainfall enhancement operation; the comparison area is not affected by the catalyst of the artificial rainfall enhancement operation;

[0137] A rainfall prediction module 402 is used to perform rainfall extrapolation prediction on the comparison area and the operation impact area respectively through a rainfall prediction model to obtain rainfall prediction data for the comparison area and rainfall prediction data for the operation impact area;

[0138] A rainfall prediction module 403 is used to obtain the actual rainfall data of the comparison area and the actual rainfall data of the operation impact area respectively;

[0139] An error calculation module 404 is used to obtain a rainfall prediction error of a rainfall prediction model based on the rainfall prediction data of the comparison area and the actual rainfall data of the comparison area;

[0140] The effect evaluation module 405 is used to obtain the absolute rainfall increase of the artificial rainmaking operation based on the rainfall prediction data of the operation influence area, the actual rainfall data of the operation influence area and the rainfall prediction error of the rainfall prediction model, and use the absolute rainfall increase of the artificial rainmaking operation as an evaluation indicator of the effect of the rainmaking operation.

[0141] In some implementations, the rainfall prediction module 402 is further configured to:

[0142] Determine the rainfall type within the scope of artificial rainfall enhancement operations; rainfall types include stratiform cloud rainfall and convective cloud rainfall, and rainfall types are divided according to the maximum rainfall intensity;

[0143] The rainfall extrapolation model corresponding to the rainfall type is selected as the rainfall prediction model.

[0144] In some implementations, the apparatus further includes a model training module 406, configured to:

[0145] A sample set is constructed based on minute grid precipitation data. Each sample in the sample set includes 2L consecutive hours of rainfall data. The K-minute grid rainfall data for the first L hours of each sample is used as input parameters, and the K-minute grid rainfall data for the last L hours is used as label data.

[0146] The PredRNN++ model is trained by the sample set to obtain the rainfall extrapolation model.

[0147] In some implementations, the rainfall prediction module 402 is specifically configured to:

[0148] Based on the time of the artificial rainfall enhancement operation, first rainfall data collected by a rainfall station in the comparison area in a first time period before the artificial rainfall enhancement operation is obtained;

[0149] Based on the first rainfall data, the rainfall prediction model is used to perform rainfall extrapolation prediction on the comparison area to obtain rainfall prediction data for the comparison area in the second time period;

[0150] Based on the time of the artificial rainfall enhancement operation, second rainfall data collected by a rainfall station in the operation impact area in a first time period before the artificial rainfall enhancement operation is obtained;

[0151] Based on the second rainfall data, the rainfall in the operation impact area is extrapolated and predicted using a rainfall prediction model to obtain rainfall prediction data for the operation impact area in the second time period.

[0152] In some implementations, the rainfall prediction module 403 is specifically configured to:

[0153] Obtain the actual rainfall data collected by n rainfall stations in the comparison area;

[0154] Obtain the actual rainfall data collected by m rainfall stations in the operation impact area.

[0155] In some implementations, the error calculation module 404 is specifically configured to:

[0156] Calculate the difference between the actual rainfall data collected by n rainfall stations in the comparison area and the rainfall forecast data in the comparison area to obtain n differences;

[0157] The average value of the n differences is taken as the rainfall prediction error of the rainfall prediction model.

[0158] In some implementations, the effect evaluation module 405 is specifically configured to:

[0159] Based on the rainfall forecast data in the operation impact area and the actual rainfall data collected by m rainfall stations in the operation impact area, the average absolute rainfall increase of the m rainfall stations in the operation impact area is obtained;

[0160] Calculate the difference between the average absolute rainfall increase of m rainfall stations in the operation influence area and the rainfall prediction error of the rainfall prediction model to obtain the absolute rainfall increase of a single rainfall station in the operation influence area;

[0161] The absolute rainfall increase of the artificial rainfall increase operation is obtained by calculating the product of the absolute rainfall increase of a single rainfall station in the operation influence area and the area of ​​the operation influence area.

[0162] It should be noted that the aforementioned explanation of the embodiment of the method for evaluating the effectiveness of rainmaking operations for intelligently predicting natural rainfall is also applicable to the device for evaluating the effectiveness of rainmaking operations for intelligently predicting natural rainfall of this embodiment, and will not be repeated here.

[0163] In order to implement the above embodiment, the present application also proposes an electronic device. Figure 5 , Figure 5 Schematic diagram of the structure of the electronic device provided in the embodiment of the present application. Figure 5As shown, the electronic device 500 includes: a processor 501, and a memory 502 communicatively connected to the processor 501; the memory 502 stores computer-executable instructions; the processor 501 executes the computer-executable instructions stored in the memory to implement the method provided in the aforementioned embodiment.

[0164] In order to implement the above embodiments, the present application also proposes a computer-readable storage medium, in which computer-executable instructions are stored. When the computer-executable instructions are executed by a processor, they are used to implement the methods provided by the above embodiments.

[0165] In order to implement the above embodiments, the present application also proposes a computer program product, including a computer program, which implements the methods provided by the above embodiments when executed by a processor.

[0166] The collection, storage, use, processing, transmission, provision and disclosure of user personal information involved in this application are in compliance with relevant laws and regulations and do not violate public order and good morals.

[0167] It should be noted that personal information from users should be collected for legitimate and reasonable purposes and should not be shared or sold outside of these legitimate uses. In addition, such collection / sharing should be carried out after receiving the user's informed consent, including but not limited to notifying the user to read the user agreement / user notice and sign the agreement / authorization including authorization of relevant user information before the user uses the function. In addition, any necessary steps should be taken to protect and safeguard access to such personal information data and ensure that others who have access to personal information data comply with its privacy policy and procedures.

[0168] The present application is expected to provide an implementation scheme for users to selectively block the use or access of personal information data. That is, the present disclosure is expected to provide hardware and / or software to prevent or block access to such personal information data. Once the personal information data is no longer needed, the risk can be minimized by limiting data collection and deleting the data. In addition, when applicable, such personal information is de-identified to protect the privacy of the user.

[0169] In the description of the aforementioned embodiments, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, without contradiction.

[0170] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features defined as "first" and "second" may explicitly or implicitly include at least one of the features. In the description of this application, the meaning of "plurality" is at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0171] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, fragment or portion of code comprising one or more executable instructions for implementing the steps of a custom logical function or process, and the scope of the preferred embodiments of the present application includes alternative implementations in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in the reverse order depending on the functions involved, which should be understood by technicians in the technical field to which the embodiments of the present application belong.

[0172] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can fetch instructions from an instruction execution system, device or apparatus and execute the instructions), or in combination with these instruction execution systems, devices or apparatuses. For the purpose of this specification, "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in combination with these instruction execution systems, devices or apparatuses. More specific examples of computer-readable media (a non-exhaustive list) include the following: an electrical connection with one or more wires (electronic device), a portable computer disk box (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), a fiber optic device, and a portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium and then editing, interpreting or processing in other suitable ways if necessary, and then stored in a computer memory.

[0173] It should be understood that the various parts of the present application can be implemented by hardware, software, firmware or a combination thereof. In the above-mentioned embodiments, multiple steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0174] A person skilled in the art may understand that all or part of the steps in the method for implementing the above-mentioned embodiment may be completed by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, which, when executed, includes one or a combination of the steps of the method embodiment.

[0175] In addition, each functional unit in each embodiment of the present application may be integrated into a processing module, or each unit may exist physically separately, or two or more units may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.

[0176] The storage medium mentioned above may be a read-only memory, a magnetic disk or an optical disk, etc. Although the embodiments of the present application have been shown and described above, it can be understood that the above embodiments are exemplary and cannot be understood as limiting the present application. A person of ordinary skill in the art may change, modify, replace and modify the above embodiments within the scope of the present application.

Claims

1. A method for evaluating the effect of rainmaking operations by intelligently predicting natural rainfall, characterized in that: The following steps are involved: Based on the scope of the artificial rainmaking operation, determining the operation impact area and the comparison area; the comparison area is not affected by the catalyst of the artificial rainmaking operation; Using a rainfall prediction model, rainfall extrapolation prediction is performed on the comparison area and the operation impact area respectively, to obtain rainfall prediction data for the comparison area and rainfall prediction data for the operation impact area; Respectively obtaining the actual rainfall data of the comparison area and the actual rainfall data of the operation impact area; Based on the rainfall prediction data of the comparison area and the actual rainfall data of the comparison area, obtaining the rainfall prediction error of the rainfall prediction model; Based on the rainfall prediction data of the operation influence area, the actual rainfall data of the operation influence area and the rainfall prediction error of the rainfall prediction model, the absolute rainfall increase of the artificial rainfall increase operation is obtained, and the absolute rainfall increase of the artificial rainfall increase operation is used as an evaluation indicator of the rainfall increase operation effect.

2. The method according to claim 1, characterized in that Before the rainfall prediction model is used to perform rainfall extrapolation prediction on the comparison area and the operation impact area respectively, the method includes: Determine the rainfall type within the scope of the artificial rainfall enhancement operation; wherein the rainfall type includes stratiform cloud rainfall and convective cloud rainfall, and the rainfall type is divided according to the maximum rainfall intensity; The rainfall extrapolation model corresponding to the rainfall type is selected as the rainfall prediction model.

3. The method according to claim 1, characterized in that The training method of the rainfall prediction model comprises: A sample set is constructed based on minute grid precipitation data, wherein each sample in the sample set includes 2L consecutive hours of rainfall data, wherein the K-minute grid rainfall data of the first L hours of each sample is used as input parameters, and the K-minute grid rainfall data of the last L hours is used as label data; The PredRNN++ model is trained by the sample set to obtain a rainfall extrapolation model.

4. The method according to claim 1, characterized in that: The method of performing rainfall extrapolation prediction on the comparison area and the operation impact area respectively through a rainfall prediction model to obtain rainfall prediction data for the comparison area and rainfall prediction data for the operation impact area comprises: Based on the time of the artificial rainfall enhancement operation, obtaining first rainfall data collected by a rainfall station in the comparison area in a first time period before the artificial rainfall enhancement operation; Based on the first rainfall data, performing rainfall extrapolation prediction on the comparison area through the rainfall prediction model to obtain rainfall prediction data for the comparison area in a second time period; Based on the time of the artificial rainfall enhancement operation, obtaining second rainfall data collected by a rainfall station in the operation impact area in the first time period before the artificial rainfall enhancement operation; Based on the second rainfall data, rainfall extrapolation prediction is performed on the operation impact area through the rainfall prediction model to obtain rainfall prediction data for the operation impact area in the second time period.

5. The method according to claim 1, characterized in that The step of respectively obtaining the real rainfall data of the comparison area and the real rainfall data of the operation impact area comprises: Obtaining real rainfall data collected by n rainfall stations in the comparison area; Obtain actual rainfall data collected by m rainfall stations in the operation impact area.

6. The method according to claim 5, characterized in that The step of obtaining the rainfall prediction error of the rainfall prediction model based on the rainfall prediction data of the comparison area and the actual rainfall data of the comparison area comprises: Calculate the difference between the actual rainfall data collected by n rainfall stations in the comparison area and the rainfall forecast data in the comparison area to obtain n differences; The average value of the n differences is taken as the rainfall prediction error of the rainfall prediction model.

7. The method according to claim 6, characterized in that The step of obtaining the absolute rainfall increase amount of the artificial rainfall increase operation based on the rainfall prediction data of the operation influence area, the actual rainfall data of the operation influence area and the rainfall prediction error of the rainfall prediction model comprises: Based on the rainfall forecast data of the operation influence area and the actual rainfall data collected by the m rainfall stations in the operation influence area, the average absolute rainfall increase of the m rainfall stations in the operation influence area is obtained; Calculate the difference between the average absolute rainfall increase of m rainfall stations in the operation influence area and the rainfall prediction error of the rainfall prediction model to obtain the absolute rainfall increase of a single rainfall station in the operation influence area; The absolute rainfall increase of the artificial rainfall increase operation is obtained by calculating the product of the absolute rainfall increase of a single rainfall station in the operation influence area and the area of ​​the operation influence area.

8. A rain-increasing operation effect evaluation device for intelligently predicting natural rainfall, characterized in that: include: The regional division module is used to determine the operation impact area and comparison area based on the scope of artificial rainfall enhancement operations; The comparison area is not affected by the catalyst of the artificial rainmaking operation; A rainfall prediction module, used to perform rainfall extrapolation prediction on the comparison area and the operation impact area respectively through a rainfall prediction model to obtain rainfall prediction data for the comparison area and rainfall prediction data for the operation impact area; A rainfall acquisition module, used to respectively acquire the real rainfall data of the comparison area and the real rainfall data of the operation impact area; An error calculation module, used for obtaining a rainfall prediction error of the rainfall prediction model based on the rainfall prediction data of the comparison area and the actual rainfall data of the comparison area; An effect evaluation module is used to obtain the absolute rainfall increase of the artificial rainmaking operation based on the rainfall prediction data of the operation influence area, the actual rainfall data of the operation influence area and the rainfall prediction error of the rainfall prediction model, and use the absolute rainfall increase of the artificial rainmaking operation as an evaluation indicator for the effect of the rainmaking operation.

9. An electronic device, characterized in that: include: A processor, and a memory communicatively connected to the processor; The memory stores computer-executable instructions; The processor executes the computer-executable instructions stored in the memory to implement the method according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer-executable instructions, which are used to implement the method according to any one of claims 1 to 7 when executed by a processor.

Images