Drought-flood rapid alternation event recognition method, device, equipment and readable storage medium

Abstract

The present application discloses a method, apparatus, device, and readable storage medium for identifying dry-wet abrupt transition events. The method obtains a hydro-meteorological sequence and the change amount of groundwater resources corresponding to a continuous period of a region; for each period, based on the hydro-meteorological sequence and the change amount of groundwater resources of the period, calculates the available water flux and the atmospheric evaporation demand corresponding to the period; based on the available water flux and the atmospheric evaporation demand of the period, calculates the residual ratio of water-heat flux of the period; when the probability of a data point is less than N%, determines that a drought occurs in the region during the corresponding period; when the probability of the data point is greater than M%, determines that a flood occurs in the region during the corresponding period; determines whether the interval between the periods when the drought and the flood occur in the region is less than a period threshold, and if so, determines that a dry-wet abrupt transition event occurs in the region. It can be seen that the present application can improve the accuracy of identifying dry-wet abrupt transition events by calculating the residual ratio of water-heat flux reflecting the degree of land-air water-heat interaction in different periods.

Description

Technical Field

[0001] The present application relates to the field of hydrological and water resources applications, and more specifically, to a method, device, equipment and readable storage medium for identifying rapid alternation events between drought and flood. Background Art

[0002] A rapid alternation event between drought and flood refers to a compound disaster in which drought and flood disasters alternate in a short period of time in a basin, reflecting the coexistence of drought and flood disasters in a short period of time. In the context of continuously intensifying climate, rapid alternation events between drought and flood are increasing. Identifying rapid alternation events between drought and flood is beneficial for formulating relevant preventive measures, preventing rapid alternation disasters between drought and flood in the basin, providing a data basis for ensuring the safety of water resources in the basin, and is of great significance in the water resources planning of the basin.

[0003] In the prior art, rapid alternation events between drought and flood are mainly identified from the perspective of atmospheric precipitation. However, rapid alternation events between drought and flood involve the water-heat interaction process between land and atmosphere. Simply starting from the meteorological perspective, it is difficult to identify the actual rapid alternation changes between drought and flood in the basin, resulting in low accuracy in identifying rapid alternation events between drought and flood in the prior art. Summary of the Invention

[0004] In view of this, the present application provides a method, device, equipment and readable storage medium for identifying rapid alternation events between drought and flood, aiming to solve the problem of low accuracy in identifying rapid alternation events between drought and flood in the prior art.

[0005] To achieve the above object, the following solutions are proposed:

[0006] A method for identifying rapid alternation events between drought and flood, comprising:

[0007] Obtaining a hydro-meteorological sequence and the change amount of groundwater resources corresponding to a continuous period of a target area;

[0008] For each period, based on the hydro-meteorological sequence and the change amount of groundwater resources of the period, calculating the available water flux and the atmospheric evaporation demand corresponding to the period; based on the available water flux and the atmospheric evaporation demand of the period, calculating the water-heat flux residual ratio of the period;

[0009] Based on each water-heat flux residual ratio, constructing a normal distribution curve;

[0010] Determining the data points matched by each water-heat flux residual ratio on the normal distribution curve;

[0011] When the probability of the data point corresponding to any period is less than N%, it is determined that a drought occurs in the target area during the corresponding period; when the probability of the data point corresponding to any period is greater than M%, it is determined that a flood occurs in the target area during the corresponding period, where N is less than M;

[0012] Determine whether the interval between the occurrence of drought and flood in the target area is less than the cycle threshold. If so, determine that a rapid drought-flood transition event occurs in the target area; if not, determine that no rapid drought-flood transition event occurs in the target area.

[0013] Optionally, each hydro-meteorological sequence includes precipitation data and regional actual evapotranspiration corresponding to the period.

[0014] Calculate the available water flux corresponding to the period based on the hydro-meteorological sequence of the period and the change in groundwater resources, including:

[0015] Calculate the available water flux of the period based on a preset water flux expression, the change in groundwater resources of the period, the precipitation data in the hydro-meteorological sequence, and the regional actual evapotranspiration.

[0016] The water flux expression is as follows:

[0017] WA i =P i -ET i +ΔTWSC i

[0018] In the formula, WA i is the available water flux of the i-th period; P i is the precipitation data of the i-th period; ET i is the regional actual evapotranspiration of the i-th period; ΔTWSC i is the change in groundwater resources of the i-th period.

[0019] Optionally, calculate the atmospheric evaporation demand corresponding to the period based on the hydro-meteorological sequence of the period and the change in groundwater resources, including:

[0020] Calculate the saturation vapor pressure, actual vapor pressure, soil heat flux, average temperature at the target height, and slope of the saturation vapor pressure curve corresponding to the period based on the hydro-meteorological sequence of the period.

[0021] Calculate the atmospheric evaporation demand of the period based on a preset evaporation demand expression, comprehensively considering the saturation vapor pressure, actual vapor pressure, soil heat flux, average temperature at the target height, slope of the saturation vapor pressure curve of the period, surface net radiation amount in the hydro-meteorological sequence, and wind speed at the target height.

[0022] The evaporation demand expression is as follows:

[0023]

[0024] In the formula, ED i is the atmospheric evaporation demand corresponding to the i-th period; Ki is the slope of the saturated water vapor pressure curve corresponding to the i-th cycle; r is the psychrometer constant; U i is the wind speed at the target height corresponding to the i-th cycle; Rn i is the net surface radiation corresponding to the i-th cycle; G i is the soil heat flux corresponding to the i-th cycle; T i is the average temperature at the target height corresponding to the i-th cycle; es i is the saturated water vapor pressure corresponding to the i-th cycle; ea i is the actual water vapor pressure corresponding to the i-th cycle.

[0025] Optionally, each hydro-meteorological sequence further includes an atmospheric temperature sequence, a target height temperature sequence, and a dew point temperature corresponding to the cycle;

[0026] Based on the hydro-meteorological sequence of the cycle, calculating the saturated water vapor pressure, actual water vapor pressure, soil heat flux, average temperature at the target height, and slope of the saturated water vapor pressure curve corresponding to the cycle, including:

[0027] Using a preset function relationship group, based on the atmospheric temperature sequence, target height temperature sequence, and dew point temperature of the cycle, calculating the saturated water vapor pressure, actual water vapor pressure, soil heat flux, average temperature at the target height, and slope of the saturated water vapor pressure curve corresponding to the cycle;

[0028] The function relationship group is as follows:

[0029]

[0030]

[0031]

[0032]

[0033]

[0034]

[0035] In the formula, Δt i is the time step of the i-th cycle; Cs is the soil heat capacity; T iy is the atmospheric temperature at the end of the i-th cycle; T iy-1 is the atmospheric temperature at the start of the i-th cycle; T imax is the highest temperature at the target height of the i-th cycle; T imin is the lowest temperature at the target height of the i-th cycle; e0 i is the air temperature of the i-th cycle; T idew is the dew point temperature of the i-th cycle.

[0036] Optionally, each hydro-meteorological sequence includes the actual regional evapotranspiration corresponding to the period.

[0037] Calculating the residual ratio of water and heat fluxes for the period based on the available water flux and atmospheric evaporation demand for the period, including:

[0038] Using a preset residual ratio expression, the actual regional evapotranspiration, available water flux, and atmospheric evaporation demand for the period to calculate the residual ratio of water and heat fluxes for the period;

[0039] The residual ratio expression is as follows:

[0040]

[0041] In the formula, WER i is the residual ratio of water and heat fluxes for the i-th period; WA i is the available water flux for the i-th period; ED i is the atmospheric evaporation demand corresponding to the i-th period; ET i is the actual regional evapotranspiration for the i-th period.

[0042] Optionally, constructing a normal distribution curve based on each residual ratio of water and heat fluxes, including:

[0043] Detrending each residual ratio of water and heat fluxes and using the z-score method to process each detrended residual ratio of water and heat fluxes to form a normal residual ratio that follows a standardized normal distribution;

[0044] Drawing a normal distribution curve based on each normal residual ratio.

[0045] Optionally, detrending each residual ratio of water and heat fluxes and using the z-score method to process each detrended residual ratio of water and heat fluxes to form a normal residual ratio that follows a standardized normal distribution, including:

[0046] Calculating the mean and standard deviation of each detrended residual ratio of water and heat fluxes;

[0047] Substituting the mean, standard deviation, and each detrended residual ratio of water and heat fluxes into a preset z-score expression to calculate the normal residual ratio corresponding to each period;

[0048] The z-score expression is as follows:

[0049]

[0050] In the formula, SWER iis the normal residual ratio corresponding to the i-th cycle; WER′ i is the detrended hydrothermal flux residual ratio corresponding to the i-th cycle; μ is the mean of the detrended hydrothermal flux residual ratios; σ is the standard deviation of the detrended hydrothermal flux residual ratios.

[0051] A device for identifying drought-flood abrupt alternation events, comprising:

[0052] An acquisition module, configured to acquire a hydro-meteorological sequence and the change amount of groundwater resources corresponding to consecutive cycles in a target area;

[0053] A calculation module, configured to, for each cycle, calculate the available water flux and the atmospheric evaporation demand corresponding to the cycle based on the hydro-meteorological sequence and the change amount of groundwater resources in the cycle; and calculate the hydrothermal flux residual ratio corresponding to the cycle based on the available water flux and the atmospheric evaporation demand in the cycle;

[0054] A construction module, configured to construct a normal distribution curve based on each hydrothermal flux residual ratio;

[0055] A determination module, configured to determine the data point matched by each hydrothermal flux residual ratio on the normal distribution curve;

[0056] A comparison module, configured to determine that a drought occurs in the target area in the corresponding cycle when the probability of the data point corresponding to any cycle is less than N%; and determine that a flood occurs in the target area in the corresponding cycle when the probability of the data point corresponding to any cycle is greater than M%, where N is less than M;

[0057] A judgment module, configured to judge whether the interval between the cycles when a drought occurs and a flood occurs in the target area is less than a cycle threshold. If so, determine that a drought-flood abrupt alternation event occurs in the target area; if not, determine that no drought-flood abrupt alternation event occurs in the target area.

[0058] A device for identifying drought-flood abrupt alternation events, comprising a memory and a processor;

[0059] The memory is configured to store a program;

[0060] The processor is configured to execute the program to implement each step of the above-mentioned method for identifying drought-flood abrupt alternation events.

[0061] A readable storage medium, on which a computer program is stored, characterized in that when the computer program is executed by a processor, each step of the above-mentioned method for identifying drought-flood abrupt alternation events is implemented.

[0062] As can be seen from the above technical solution, the method for identifying drought-flood rapid alternation events provided by this application can obtain the hydrometeorological sequence and the change amount of groundwater resources corresponding to consecutive periods in the target area; for each period, based on the hydrometeorological sequence and the change amount of groundwater resources in the period, calculate the available water flux and the atmospheric evaporation demand corresponding to the period; based on the available water flux and the atmospheric evaporation demand in the period, calculate the residual ratio of water-heat flux in the period; on this basis, since the available water flux can usually be defined as the sum of the total available surface water resources and the change amount of groundwater resources reserves, the calculated available water flux reflects surface water information and groundwater information; since the atmospheric evaporation demand corresponds to the demand of the atmosphere for water, that is, the drying or evaporation ability, the calculated atmospheric evaporation demand can reflect the water-heat interaction process; based on this, the residual ratio of water-heat flux calculated using the available water flux and the atmospheric evaporation demand reflects the degree of water-heat interaction between the land and the atmosphere, not only considering the precipitation information of the atmosphere, but also comprehensively including the water information of the land and the water-heat interaction information between the atmosphere and the land, further improving the reliability of the residual ratio of water-heat flux; based on each residual ratio of water-heat flux, construct a normal distribution curve; determine the data points matched by each residual ratio of water-heat flux on the normal distribution curve; when the probability of the data point corresponding to any period is less than N%, determine that a drought occurs in the target area during the corresponding period; when the probability of the data point corresponding to any period is greater than M%, determine that a flood occurs in the target area during the corresponding period, where N is less than M; based on this, by comparing the probability size corresponding to the residual ratio of water-heat flux in each period with the preset threshold, the degree of water-heat interaction in each period can be evaluated, and further evaluate whether a drought or a flood occurs in the target area in each period; subsequently, it can be judged whether the interval between the periods of drought and flood occurrence in the target area is less than the period threshold, if so, determine that a drought-flood rapid alternation event occurs in the target area; if not, determine that no drought-flood rapid alternation event occurs in the target area; based on this, it can be determined whether a drought-flood rapid alternation event occurs in the target area in a short time; it can be seen that this application can complete drought identification and flood identification by calculating the residual ratio of water-heat flux reflecting the land-atmosphere water-heat interaction degree in different periods, and determine whether a drought-flood rapid alternation disaster occurs by judging whether drought and flood occur in a short time.

[0063] In the process of calculating the residual ratio of water-heat flux in this application, it comprehensively includes surface water information, groundwater information and the atmospheric evaporation demand reflecting the land-atmosphere water-heat interaction process, further ensuring that the residual ratio of water-heat flux can reflect the degree of land-atmosphere water-heat interaction and improving the identification reliability of this application. BRIEF DESCRIPTION OF THE DRAWINGS

[0064] To more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the accompanying drawings required for the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only the embodiments of the present application. For those of ordinary skill in the art, without creative efforts, other accompanying drawings can also be obtained based on the provided drawings.

[0065] Figure 1 Flowchart of a method for identifying dry-wet abrupt transition events disclosed in an embodiment of the present application;

[0066] Figure 2 Block diagram of a device for identifying dry-wet abrupt transition events disclosed in an embodiment of the present application;

[0067] Figure 3 Hardware block diagram of a device for identifying dry-wet abrupt transition events disclosed in an embodiment of the present application. Detailed implementation manners

[0068] The following will clearly and completely describe the technical solutions in the embodiments of the present application in conjunction with the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all of them. Based on the embodiments of the present application, all other embodiments obtained by those of ordinary skill in the art without creative efforts belong to the scope of protection of the present application.

[0069] The method for identifying dry-wet abrupt transition events of the present application can be used in many general or special computing device environments or configurations. For example: personal computers, server computers, handheld devices or portable devices, tablet devices, multi-processor devices, distributed computing environments including any of the above devices or equipment, and so on.

[0070] Next, in conjunction with Figure 1 A detailed introduction to the method for identifying dry-wet abrupt transition events of the present application is as follows:

[0071] Step S1: Obtain the hydrometeorological sequence corresponding to the continuous period of the target area and the change amount of groundwater resources.

[0072] Specifically, hydrological data and meteorological data of the target area from 1982 to 2015 can be obtained, and with one month as a period, based on the hydrological data and meteorological data of the target area from 1982 to 2015, multiple continuous periods of hydrometeorological data can be constructed.

[0073] Among them, the meteorological data can be sourced from the fifth-generation reanalysis global climate dataset of the European Centre for Medium-Range Weather Forecasts, the change in groundwater resources can be sourced from the terrestrial water storage change dataset of the Gravity Recovery and Climate Experiment (GRACE) launched by the National Aeronautics and Space Administration (NASA), and the actual regional evapotranspiration in the hydrological data can be sourced from the global terrestrial evapotranspiration dataset based on the PEW model of the National Tibetan Plateau Data Center of Sciences.

[0074] Each hydro-meteorological sequence may include the target high-altitude wind speed, precipitation data, actual regional evapotranspiration, atmospheric temperature sequence, surface net radiation, target high-altitude temperature sequence, and dew point temperature for the corresponding period.

[0075] Each change in groundwater resources may be the difference between the groundwater resource storage at the start of the corresponding period and the groundwater resource storage at the end of the period.

[0076] Step S2: For each period, based on the hydro-meteorological sequence and the change in groundwater resources of the period, calculate the available water flux and atmospheric evaporation demand corresponding to the period; based on the available water flux and atmospheric evaporation demand of the period, calculate the residual ratio of water-heat flux.

[0077] Specifically, the available water flux and atmospheric evaporation demand for each period can be calculated based on the hydro-meteorological sequence and the change in groundwater resources of each period.

[0078] Numerical operations can be performed on the available water flux and atmospheric evaporation demand for each period to calculate the residual ratio of water-heat flux for each period.

[0079] Among them, the available water flux of the target area in the water resources field is usually defined as the total available water resources, that is, within the foreseeable period, on the basis of overall consideration of domestic, industrial, and ecological water use, the maximum amount of water that can be used once in the local water resources through economically reasonable and technically feasible measures. It can be obtained by adding the non-repetitive amount of surface water resources and groundwater resources and the surface available water resources. In the natural environment, the water use of human activities does not need to be overall considered, so the available water flux can be further defined as the sum of the total surface available water resources and the change in groundwater resource storage.

[0080] Atmospheric evaporation demand is a complex concept. Different from the actual regional evapotranspiration that represents the upward water flux from the soil, free water surface, and / or vegetation to the atmosphere, atmospheric evaporation demand does not correspond to the land-atmosphere flux, but rather to the atmospheric demand for water, i.e., the drying or evaporation capacity. It is a function of the atmospheric boundary layer state and dynamics and is defined as the maximum amount of water that can enter the atmosphere through evaporation from a free water surface under certain meteorological conditions. Usually, the potential evapotranspiration (PET) can be used to represent it. The potential evapotranspiration cannot be directly measured and needs to be calculated through climate parameters, surface types, surface waters such as lakes and seas, soil types, vegetation cover, etc. In a natural environment, the Penman-Monteith formula of FAO 56 is usually used to calculate the potential evapotranspiration of the target area.

[0081] Each hydrothermal flux residual ratio can be used to determine the degree of land-atmosphere hydrothermal interaction in the corresponding period.

[0082] Step S3: Based on each hydrothermal flux residual ratio, construct a normal distribution curve.

[0083] Specifically, each hydrothermal flux residual ratio can be standardized, and a standardized normal distribution curve can be constructed based on the standardized hydrothermal flux residual ratios.

[0084] Step S4: Determine the data points matched by each hydrothermal flux residual ratio on the normal distribution curve.

[0085] Specifically, the data points corresponding to each hydrothermal flux residual ratio can be found on the normal distribution curve.

[0086] Step S5: When the probability of the data point corresponding to any period is less than N%, it is determined that a drought occurs in the target area during the corresponding period; when the probability of the data point corresponding to any period is greater than M%, it is determined that a flood occurs in the target area during the corresponding period.

[0087] Specifically, the hydrothermal flux residual ratio corresponding to the data point with a probability of N% is the drought critical value of the target area; the hydrothermal flux residual ratio corresponding to the data point with a probability of M% is the flood critical value of the target area. N is less than M.

[0088] Among them, N can be 10. M can be 90.

[0089] If the residual ratio of hydrothermal flux in any period is less than that corresponding to the data points with a probability of N%, the degree of land-atmosphere hydrothermal interaction in the corresponding period is poor and drought is likely to occur; if the residual ratio of hydrothermal flux in any period is greater than that corresponding to the data points with a probability of M%, the degree of land-atmosphere hydrothermal interaction in the corresponding period is high and waterlogging is likely to occur; if the probability of the data points corresponding to any period is greater than N% and less than M%, the degree of land-atmosphere hydrothermal interaction in the corresponding period is normal and neither drought nor waterlogging occurs in the target area during the corresponding period.

[0090] Therefore, the probability corresponding to the data points in each period can be compared with N% and M% to determine whether drought, waterlogging or normality occurs in the target area during the corresponding period.

[0091] Step S6: Determine whether the interval between the periods of drought and waterlogging occurring in the target area is less than the period threshold. If so, it is determined that a rapid alternation event of drought and waterlogging occurs in the target area; if not, it is determined that no rapid alternation event of drought and waterlogging occurs in the target area.

[0092] Specifically, the drought period corresponding to the drought and the waterlogging period corresponding to the waterlogging can be determined, the interval between the drought period and the waterlogging period can be calculated, and the interval between the periods is compared with the preset period threshold. If the interval between the periods is less than the period threshold, there is a rapid alternation event of drought and waterlogging in the target area; if it is greater than the period threshold, no rapid alternation event of drought and waterlogging occurs in the target area.

[0093] Among them, the period threshold can be 3.

[0094] As can be seen from the above technical solution, the method for identifying dry-wet abrupt alternation events provided by this application can obtain the hydrometeorological sequence and the change amount of groundwater resources corresponding to consecutive periods in the target area; for each period, based on the hydrometeorological sequence and the change amount of groundwater resources in the period, calculate the available water flux and the atmospheric evaporation demand corresponding to the period; based on the available water flux and the atmospheric evaporation demand in the period, calculate the residual ratio of water and heat fluxes in the period; on this basis, since the available water flux can usually be defined as the sum of the total available surface water resources and the change amount of groundwater resources reserves, thus, the calculated available water flux reflects surface water information and groundwater information; since the atmospheric evaporation demand corresponds to the demand of the atmosphere for water, that is, the drying or evaporation ability, thus, the calculated atmospheric evaporation demand can reflect the water and heat interaction process; based on this, the residual ratio of water and heat fluxes calculated using the available water flux and the atmospheric evaporation demand reflects the degree of water and heat interaction between the land and the atmosphere, not only considering the precipitation information of the atmosphere, but also comprehensively including the water information of the land and the water and heat interaction information between the atmosphere and the land, further improving the reliability of the residual ratio of water and heat fluxes; based on each residual ratio of water and heat fluxes, construct a normal distribution curve; determine the data points matched by each residual ratio of water and heat fluxes on the normal distribution curve; when the probability of the data point corresponding to any period is less than N%, determine that a drought occurs in the target area in the corresponding period; when the probability of the data point corresponding to any period is greater than M%, determine that a flood occurs in the target area in the corresponding period, where N is less than M; based on this, by comparing the probability size corresponding to the residual ratio of water and heat fluxes in each period with a preset threshold, the degree of water and heat interaction in each period can be evaluated, and further evaluate whether a drought or a flood occurs in the target area in each period; subsequently, it can be judged whether the interval between the periods of drought and flood occurrence in the target area is less than the period threshold, if so, determine that a dry-wet abrupt alternation event occurs in the target area; if not, determine that no dry-wet abrupt alternation event occurs in the target area; based on this, it can be determined whether a dry-wet abrupt alternation event occurs in the target area in a short time; it can be seen that this application can complete drought identification and flood identification by calculating the residual ratio of water and heat fluxes reflecting the land-atmosphere water and heat interaction degree in different periods, and determine whether a dry-wet abrupt alternation disaster occurs by judging whether drought and flood occur in a short time.

[0095] In the process of calculating the residual ratio of water and heat fluxes in this application, surface water information, groundwater information and the atmospheric evaporation demand reflecting the land-atmosphere water and heat interaction process are comprehensively included, further ensuring that the residual ratio of water and heat fluxes can reflect the degree of land-atmosphere water and heat interaction and improving the identification reliability of this application.

[0096] In addition, the accurate identification of dry-wet abrupt change events can comprehensively and effectively monitor or predict the changes in the dry-wet conditions of the target area, so as to more accurately identify dry-wet abrupt change events, which is of great significance for the analysis and prevention of drought and flood disasters in the target area.

[0097] In some embodiments of the present application, the process of calculating the available water flux corresponding to the period based on the hydrometeorological sequence and the change in groundwater resources in step S2 is described in detail as follows:

[0098] S20. Calculate the available water flux of the period based on a preset water flux expression, the change in groundwater resources in the period, the precipitation data in the hydrometeorological sequence, and the actual evapotranspiration of the region.

[0099] Specifically, the change in groundwater resources in each period, the precipitation data in the hydrometeorological sequence, and the actual evapotranspiration of the region can be substituted into the water flux expression to calculate the available water flux of each period.

[0100] Among them, the water flux expression can be as follows:

[0101] WA i =P i -ET i +ΔTWSC i

[0102] In the formula, WA i is the available water flux of the i-th period; P i is the precipitation data of the i-th period; ET i is the actual evapotranspiration of the region in the i-th period; ΔTWSC i is the change in groundwater resources in the i-th period.

[0103] It can be seen from the above technical solutions that this embodiment provides an optional method for calculating the available water flux. Through the above method, the available water flux can be further calculated by comprehensively considering the precipitation, the actual evapotranspiration of the region, and the change in groundwater resources, so that the available water flux can reflect the actual land-air interaction process of the target area.

[0104] In some embodiments of the present application, the process of calculating the atmospheric evaporation demand corresponding to the period based on the hydrometeorological sequence and the change in groundwater resources in step S2 is described in detail as follows:

[0105] S21. Calculate the saturation vapor pressure, actual vapor pressure, soil heat flux, average temperature at the target height, and slope of the saturation vapor pressure curve corresponding to the period based on the hydrometeorological sequence of the period.

[0106] Specifically, based on the target high-altitude wind speed, atmospheric temperature sequence, target high-altitude temperature sequence, and dew point temperature in each cycle of the hydrometeorological sequence, the saturation vapor pressure, actual vapor pressure, soil heat flux, average temperature at the target height, and slope of the saturation vapor pressure curve corresponding to each cycle can be calculated.

[0107] S22. Based on a preset evaporation demand expression, comprehensively calculate the atmospheric evaporation demand of the cycle by using the saturation vapor pressure, actual vapor pressure, soil heat flux, average temperature at the target height, slope of the saturation vapor pressure curve, and target high-altitude wind speed in the hydrometeorological sequence of the cycle.

[0108] Specifically, the saturation vapor pressure, actual vapor pressure, soil heat flux, average temperature at the target height, slope of the saturation vapor pressure curve, surface net radiation in the hydrometeorological sequence, and target high-altitude wind speed of each cycle can be substituted into the preset evaporation demand expression to calculate the atmospheric evaporation demand of each cycle.

[0109] Among them, the evaporation demand expression can be as follows:

[0110]

[0111] In the formula, ED i is the atmospheric evaporation demand corresponding to the i-th cycle; K i is the slope of the saturation vapor pressure curve corresponding to the i-th cycle; r is the psychrometer constant; U i is the target high-altitude wind speed corresponding to the i-th cycle; Rn i is the surface net radiation corresponding to the i-th cycle; G i is the soil heat flux corresponding to the i-th cycle; T i is the average temperature at the target height corresponding to the i-th cycle; es i is the saturation vapor pressure corresponding to the i-th cycle; ea i is the actual vapor pressure corresponding to the i-th cycle.

[0112] The height of the target height can be set according to actual needs. For example, it can be set to 2m.

[0113] It can be seen from the above technical solutions that this embodiment provides an optional method for calculating the atmospheric evaporation demand. Through the above method, the calculation of the atmospheric evaporation demand can be further completed from multiple perspectives, and the accuracy of the atmospheric evaporation demand can be further improved.

[0114] In some embodiments of the present application, the process of step S21, calculating the saturation vapor pressure, actual vapor pressure, soil heat flux, average temperature at the target height, and slope of the saturation vapor pressure curve corresponding to the period based on the hydrometeorological sequence of the period, is described in detail as follows:

[0115] S210. Using a preset function group, calculate the saturation vapor pressure, actual vapor pressure, soil heat flux, average temperature at the target height, and slope of the saturation vapor pressure curve corresponding to the period based on the atmospheric temperature sequence, target height temperature sequence, and dew point temperature of the period.

[0116] Specifically, the atmospheric temperature sequence, target height temperature sequence, and dew point temperature of each period can be substituted into the preset function group to calculate the saturation vapor pressure, actual vapor pressure, soil heat flux, average temperature at the target height, and slope of the saturation vapor pressure curve corresponding to each period.

[0117] Among them, the function group can be as follows:

[0118]

[0119]

[0120]

[0121]

[0122]

[0123]

[0124] In the formula, Δt i is the time step of the i-th period; Cs is the soil heat capacity; T iy is the atmospheric temperature at the end of the i-th period; T iy-1 is the atmospheric temperature at the start of the i-th period; T imax is the highest temperature at the target height of the i-th period; T imin is the lowest temperature at the target height of the i-th period; e0 i is the air temperature of the i-th period; T idew is the dew point temperature of the i-th period.

[0125] It can be seen from the above technical solutions that this embodiment provides an optional method for calculating the saturation vapor pressure, actual vapor pressure, soil heat flux, average temperature at the target height, and slope of the saturation vapor pressure curve, thereby better completing the acquisition of the atmospheric evaporation demand.

[0126] In some embodiments of the present application, the process of calculating the hydrothermal flux residual ratio of the cycle based on the available water flux and the atmospheric evaporation demand of the cycle in step S2 is described in detail as follows:

[0127] S223. Calculate the hydrothermal flux residual ratio of the cycle by using a preset residual ratio expression, the actual evapotranspiration of the region of the cycle, the available water flux, and the atmospheric evaporation demand.

[0128] Specifically, the actual evapotranspiration of the region, the available water flux, and the atmospheric evaporation demand of each cycle can be substituted into the preset residual ratio expression to calculate the hydrothermal flux residual ratio of each cycle.

[0129] Among them, the residual ratio expression can be as follows:

[0130]

[0131] In the formula, WER i is the hydrothermal flux residual ratio of the i-th cycle; WA i is the available water flux of the i-th cycle; ED i is the atmospheric evaporation demand corresponding to the i-th cycle; ET i is the actual evapotranspiration of the region of the i-th cycle.

[0132] It can be seen from the above technical solution that this embodiment provides an optional method for calculating the hydrothermal flux residual ratio. Through the above method, the hydrothermal flux residual ratio can be calculated by comprehensively considering the actual evapotranspiration of the region, the available water flux, and the atmospheric evaporation demand, further improving the characterization ability of the hydrothermal flux residual ratio for the land-atmosphere hydrothermal interaction process and further improving the accuracy of the present application.

[0133] In some embodiments of the present application, the process of constructing a normal distribution curve based on each hydrothermal flux residual ratio in step S3 is described in detail as follows:

[0134] S30. Detrend each hydrothermal flux residual ratio and use the z-score method to process each detrended hydrothermal flux residual ratio to form a normal residual ratio that follows a standardized normal distribution.

[0135] Specifically, the optimal fitting method can be used to construct a trend line for fitting the curve formed by each hydrothermal flux residual ratio;

[0136] Based on the trend line, the trend value of each hydrothermal flux residual ratio can be determined, and the corresponding trend value can be subtracted from each hydrothermal flux residual ratio to obtain the detrended hydrothermal flux residual ratio;

[0137] The residual ratio normalization expression constructed based on the z-score method can be used to process each detrended hydrothermal flux residual ratio to form a normal residual ratio that follows a standardized normal distribution.

[0138] S31. Draw a normal distribution curve based on each normal residual ratio.

[0139] Specifically, a normal distribution curve is drawn with the mean of each normal residual ratio as the central position of the normal distribution curve and the standard deviation of each normal residual ratio as the shape parameter of the normal distribution curve.

[0140] As can be seen from the above technical solution, this embodiment provides an optional method for constructing a normal distribution curve based on each hydrothermal flux residual ratio. Through the above method, the probability corresponding to each hydrothermal flux residual ratio can be better determined using the normal distribution curve, and the identification of drought and flood disasters can be better completed.

[0141] In some embodiments of the present application, the process of using the z-score method in step S30 to process each detrended hydrothermal flux residual ratio to form a normal residual ratio that follows a standardized normal distribution is described in detail as follows:

[0142] S300. Calculate the mean and standard deviation of each detrended hydrothermal flux residual ratio.

[0143] Specifically, the ratio between the sum of each detrended hydrothermal flux residual ratio and the number value of each hydrothermal flux residual ratio can be used as the mean of each detrended hydrothermal flux residual ratio.

[0144] The standard deviation of each detrended hydrothermal flux residual ratio can be calculated based on the mean of each detrended hydrothermal flux residual ratio.

[0145] S301. Substitute the mean, standard deviation, and each detrended hydrothermal flux residual ratio into a preset z-score expression to calculate the normal residual ratio corresponding to each period.

[0146] Specifically, the residual ratio normalization expression can be as follows:

[0147]

[0148] In the formula, μ is the mean of each hydrothermal flux residual ratio; σ is the standard deviation of each hydrothermal flux residual ratio; SWER i is the normal residual ratio of the i-th period.

[0149] As can be seen from the above technical solution, this embodiment provides an optional method for converting the hydrothermal flux residual ratio into a normal residual ratio, and the probability of each hydrothermal flux residual ratio can be further determined through the above method.

[0150] Next, the device for identifying dry-wet abrupt transition events provided by this application will be introduced in detail. The device for identifying dry-wet abrupt transition events provided below can be compared with the method for identifying dry-wet abrupt transition events provided above. Figure 2

[0151] See Figure 2 It can be found that the device for identifying dry-wet abrupt transition events may include:

[0152] An acquisition module 10, configured to acquire a hydro-meteorological sequence corresponding to a continuous period of a target area and the change amount of groundwater resources;

[0153] A calculation module 20, configured to, for each period, calculate the available water flux and the atmospheric evaporation demand corresponding to the period based on the hydro-meteorological sequence and the change amount of groundwater resources of the period; and calculate the hydrothermal flux residual ratio corresponding to the period based on the available water flux and the atmospheric evaporation demand of the period;

[0154] A construction module 30, configured to construct a normal distribution curve based on each hydrothermal flux residual ratio;

[0155] A determination module 40, configured to determine the data point matched by each hydrothermal flux residual ratio on the normal distribution curve;

[0156] A comparison module 50, configured to determine that a drought occurs in the target area during the corresponding period when the probability of the data point corresponding to any period is less than N%; and determine that a flood occurs in the target area during the corresponding period when the probability of the data point corresponding to any period is greater than M%, where N is less than M;

[0157] A judgment module 60, configured to judge whether the interval between the periods of drought and flood occurrence in the target area is less than a period threshold. If so, it is determined that a dry-wet abrupt transition event occurs in the target area; if not, it is determined that no dry-wet abrupt transition event occurs in the target area.

[0158] Furthermore, the calculation module may include:

[0159] An available water flux calculation unit, configured to calculate the available water flux of the period based on a preset water flux expression, the change amount of groundwater resources of the period, the precipitation data in the hydro-meteorological sequence, and the actual evapotranspiration of the area;

[0160] The water flux expression is as follows:

[0161] WA i =P​i -ET i +ΔTWSC i

[0162] In the formula, WA i is the available water flux in the i-th cycle; P i is the precipitation data in the i-th cycle; ET i is the actual evapotranspiration in the i-th cycle; ΔTWSC i is the change in groundwater resources in the i-th cycle.

[0163] Furthermore, the calculation module may further include:

[0164] A soil heat flux calculation unit for calculating the saturation vapor pressure, actual vapor pressure, soil heat flux, average temperature at the target height, and slope of the saturation vapor pressure curve corresponding to the cycle based on the hydrometeorological sequence of the cycle;

[0165] An atmospheric evaporation demand calculation unit for calculating the atmospheric evaporation demand of the cycle based on a preset evaporation demand expression, integrating the saturation vapor pressure, actual vapor pressure, soil heat flux, average temperature at the target height, slope of the saturation vapor pressure curve, surface net radiation amount, and target height wind speed in the hydrometeorological sequence of the cycle;

[0166] The evaporation demand expression is as follows:

[0167]

[0168] In the formula, ED i is the atmospheric evaporation demand corresponding to the i-th cycle; K i is the slope of the saturation vapor pressure curve corresponding to the i-th cycle; r is the psychrometer constant; U i is the target height wind speed corresponding to the i-th cycle; Rn i is the surface net radiation amount corresponding to the i-th cycle; G i is the soil heat flux corresponding to the i-th cycle; T i is the average temperature at the target height corresponding to the i-th cycle; es i is the saturation vapor pressure corresponding to the i-th cycle; ea i is the actual vapor pressure corresponding to the i-th cycle.

[0169] Furthermore, the soil heat flux calculation unit may include:

[0170] The function relation group uses components to calculate the saturated water vapor pressure, actual water vapor pressure, soil heat flux, average temperature at the target height, and slope of the saturated water vapor pressure curve corresponding to the period based on the periodic atmospheric temperature sequence, target height temperature sequence, and dew point temperature using a preset function relation group;

[0171] The function relation group is as follows:

[0172]

[0173]

[0174]

[0175]

[0176]

[0177]

[0178] In the formula, Δt i is the time step of the i-th period; Cs is the soil heat capacity; T iy is the atmospheric temperature at the end of the i-th period; T iy-1 is the atmospheric temperature at the start of the i-th period; T imax is the highest temperature at the target height in the i-th period; T imin is the lowest temperature at the target height in the i-th period; e0 i is the air temperature in the i-th period; T idew is the dew point temperature in the i-th period.

[0179] Furthermore, the calculation module may further include:

[0180] A water-heat flux residual ratio calculation unit for calculating the water-heat flux residual ratio of the period using a preset residual ratio expression, the regional actual evapotranspiration, available water flux, and atmospheric evaporation demand of the period;

[0181] The residual ratio expression is as follows:

[0182]

[0183] In the formula, WER i is the water-heat flux residual ratio of the i-th period; WA i is the available water flux of the i-th period; ED i is the atmospheric evaporation demand corresponding to the i-th period; ET i is the regional actual evapotranspiration of the i-th period.

[0184] Further, the building module may include:

[0185] A hydrothermal flux residual ratio processing unit, configured to detrend each hydrothermal flux residual ratio, and process each detrended hydrothermal flux residual ratio by using the z-score method to form a normal residual ratio that follows a standardized normal distribution;

[0186] A normal distribution curve drawing unit, configured to draw a normal distribution curve based on each normal residual ratio.

[0187] Further, the hydrothermal flux residual ratio processing unit may include:

[0188] A standard deviation calculation component, configured to calculate the mean and standard deviation of each detrended hydrothermal flux residual ratio;

[0189] A normal residual ratio calculation component, configured to substitute the mean, standard deviation, and each detrended hydrothermal flux residual ratio into a preset z-score expression to calculate the normal residual ratio corresponding to each period;

[0190] The z-score expression is as follows:

[0191]

[0192] In the formula, SWER i is the normal residual ratio corresponding to the i-th period; WER' i is the detrended hydrothermal flux residual ratio corresponding to the i-th period; μ is the mean of each detrended hydrothermal flux residual ratio; σ is the standard deviation of each detrended hydrothermal flux residual ratio.

[0193] The device for identifying dry-wet abrupt transition events provided by the embodiments of the present application can be applied to devices for identifying dry-wet abrupt transition events, such as PC terminals, cloud platforms, servers, server clusters, etc. Optionally, Figure 3 shows a hardware structure block diagram of a device for identifying dry-wet abrupt transition events. Referring to Figure 3 , the hardware structure of the device for identifying dry-wet abrupt transition events may include: at least one processor 1, at least one communication interface 2, at least one memory 3, and at least one communication bus 4;

[0194] In the embodiments of the present application, the number of the processor 1, the communication interface 2, the memory 3, and the communication bus 4 is at least one, and the processor 1, the communication interface 2, and the memory 3 complete mutual communication through the communication bus 4;

[0195] The processor 1 may be a central processing unit (CPU), or an application specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present invention, etc.;

[0196] The memory 3 may include a high-speed RAM memory, and may also include a non-volatile memory, such as at least one disk memory;

[0197] Wherein, the memory stores a program, and the processor can call the program stored in the memory, and the program is used for:

[0198] Obtain the hydro-meteorological sequence and the change amount of groundwater resources corresponding to consecutive periods of the target area;

[0199] For each period, based on the hydro-meteorological sequence and the change amount of groundwater resources in the period, calculate the available water flux and the atmospheric evaporation demand corresponding to the period; based on the available water flux and the atmospheric evaporation demand in the period, calculate the residual ratio of water-heat flux in the period;

[0200] Based on each residual ratio of water-heat flux, construct a normal distribution curve;

[0201] Determine the data points matched by each residual ratio of water-heat flux on the normal distribution curve;

[0202] When the probability of the data point corresponding to any period is less than N%, it is determined that a drought occurs in the target area in the corresponding period; when the probability of the data point corresponding to any period is greater than M%, it is determined that a waterlogging occurs in the target area in the corresponding period, where N is less than M;

[0203] Judge whether the interval between the periods of drought and waterlogging in the target area is less than the period threshold. If so, it is determined that a rapid alternation event of drought and waterlogging occurs in the target area; if not, it is determined that no rapid alternation event of drought and waterlogging occurs in the target area.

[0204] Optionally, the refined functions and extended functions of the program may refer to the above description.

[0205] The embodiment of the present application also provides a readable storage medium, which can store a program suitable for being executed by a processor, and the program is used for:

[0206] Obtain the hydro-meteorological sequence and the change amount of groundwater resources corresponding to consecutive periods of the target area;

[0207] For each period, based on the hydrometeorological sequence and the change in groundwater resources of the period, calculate the available water flux and the atmospheric evaporation demand corresponding to the period; based on the available water flux and the atmospheric evaporation demand of the period, calculate the residual ratio of water and heat fluxes for the period;

[0208] Based on the residual ratios of water and heat fluxes, construct a normal distribution curve;

[0209] Determine the data points corresponding to each residual ratio of water and heat fluxes on the normal distribution curve;

[0210] When the probability of the data point corresponding to any period is less than N%, it is determined that a drought occurs in the target area during the corresponding period; when the probability of the data point corresponding to any period is greater than M%, it is determined that a flood occurs in the target area during the corresponding period, where N is less than M;

[0211] Judge whether the interval between the periods of drought and flood in the target area is less than the period threshold. If so, it is determined that a rapid alternation event of drought and flood occurs in the target area; if not, it is determined that no rapid alternation event of drought and flood occurs in the target area.

[0212] Optionally, the refinement function and the extension function of the program can refer to the above description.

[0213] Finally, it should also be noted that in this article, relational terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the term "comprising", "including" or any other variant thereof is intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements not only includes those elements, but also includes other elements not expressly listed, or also includes elements inherent in such process, method, article or device. Without further limitation, an element defined by the statement "including a..." does not exclude the existence of additional identical elements in the process, method, article or device including the element.

[0214] The various embodiments in this specification are described in a progressive manner, and the key points of each embodiment are the differences from other embodiments. The same or similar parts between the various embodiments can be referred to each other.

[0215] The foregoing description of the disclosed embodiments enables those skilled in the art to implement or use the present application. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the present application. The various embodiments of the present application can be combined with each other. Therefore, the present application will not be limited to the embodiments shown herein, but rather is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for identifying drought-flood sudden change events, characterized in that: include: Obtain the hydrological and meteorological sequences and groundwater resource changes corresponding to the continuous periods in the target area; For each period, based on the hydrological and meteorological sequence of the period and the change in groundwater resources, the available water flux and atmospheric evaporation demand corresponding to the period are calculated, wherein the hydrological and meteorological sequence includes the actual regional evapotranspiration of the period; the water-heat flux residual ratio of the period is calculated using a preset residual ratio expression, the actual regional evapotranspiration of the period, the available water flux and the atmospheric evaporation demand; Based on the residual ratios of each water-heat flux, a normal distribution curve is constructed; determining a matching data point of each water heat flux residual ratio on the normal distribution curve; When the probability of the corresponding data point in any period is less than N%, it is determined that the target area has suffered a drought in the corresponding period; when the probability of the corresponding data point in any period is greater than M%, it is determined that the target area has suffered a waterlogging in the corresponding period, where N is less than M; Determine whether the period interval between the occurrence of drought and flood in the target area is less than the period threshold, if so, determine that the target area has experienced a drought-flood transition event; if not, determine that the target area has not experienced a drought-flood transition event; The residual ratio expression is as follows: In the formula, WER i is the residual ratio of water heat flux in the ith cycle; WA i is the available water flux in the ith cycle; ED i is the atmospheric evaporation demand corresponding to the i-th period; ET i is the actual regional evapotranspiration in the ith period.

2. The drought-flood rapid transition event identification method according to claim 1, characterized in that: Each hydro-meteorological series includes the precipitation data of the corresponding period and the actual evapotranspiration of the region; Based on the hydrological and meteorological sequence of the period and the change in groundwater resources, the available water flux corresponding to the period is calculated, including: Based on the preset water flux expression, the change of groundwater resources in the period, the precipitation data in the hydrometeorological sequence and the actual evapotranspiration of the region, the available water flux of the period is calculated; The water flux expression is as follows: WA i =P i -AND i +ΔTWSC i In the formula, WA i is the available water flux in the ith cycle; P i is the precipitation data of the ith period; ET i is the actual regional evapotranspiration in the i-th cycle; ΔTWSC i is the change of groundwater resources in the ith period.

3. The drought-flood rapid transition event identification method according to claim 1, characterized in that: Based on the hydrological and meteorological sequence of the period and the change in groundwater resources, the atmospheric evaporation demand corresponding to the period is calculated, including: Based on the hydrological and meteorological sequence of the period, the saturated water vapor pressure, the actual water vapor pressure, the soil heat flux, the average temperature at the target height and the slope of the saturated water vapor pressure curve corresponding to the period are calculated; Based on the preset evaporation demand expression, the atmospheric evaporation demand of the period is calculated by comprehensively considering the saturated water vapor pressure, actual water vapor pressure, soil heat flux, average temperature at the target height, slope of the saturated water vapor pressure curve, net surface radiation in the hydro-meteorological sequence, and wind speed at the target height; The evaporation demand expression is as follows: In the formula, ED i is the atmospheric evaporation demand corresponding to the i-th period; K i is the slope of the saturated water vapor pressure curve corresponding to the i-th cycle; r is the hygrometer constant; U i is the target high wind speed corresponding to the ith period; Rn i is the surface net radiation corresponding to the ith period; G i is the soil heat flux corresponding to the ith period; T i is the average temperature of the target height corresponding to the i-th cycle; es i is the saturated water vapor pressure corresponding to the i-th cycle; ea i is the actual water vapor pressure corresponding to the ith period.

4. The method for identifying drought-flood sudden transition events according to claim 3, characterized in that: Each hydro-meteorological series also includes the atmospheric temperature series, target high-altitude temperature series and dew-point temperature of the corresponding period; Based on the hydrological and meteorological sequence of the period, the saturated water vapor pressure, actual water vapor pressure, soil heat flux, target high average temperature and the slope of the saturated water vapor pressure curve corresponding to the period are calculated, including: Calculate the saturated water vapor pressure, actual water vapor pressure, soil heat flux, target high temperature average and the slope of the saturated water vapor pressure curve corresponding to the period based on the atmospheric temperature sequence, target high temperature sequence and dew point temperature of the period using a preset functional relationship group; The functional relationship group is as follows: In the formula, Δt i is the time step of the ith cycle; Cs is the soil heat capacity; T iy is the atmospheric temperature at the end of the ith period; T iy-1 is the atmospheric temperature at the beginning of the ith period; T imax is the target maximum temperature at high altitude in the i-th cycle; T imin is the target high minimum temperature of the i-th cycle; e0 i is the air temperature of the ith cycle; T idew is the dew point temperature of the ith cycle.

5. The drought-flood rapid transition event identification method according to claim 1, characterized in that: The normal distribution curve is constructed based on each water heat flux residual ratio, including: Each water-heat flux residual ratio is detrended, and the z-score method is used to process each detrended water-heat flux residual ratio to form a normal residual ratio that obeys a standardized normal distribution; Based on each normal residual ratio, a normal distribution curve is plotted.

6. The drought-flood rapid transition event identification method according to claim 5, characterized in that: The z-score method is used to process each detrended water and heat flux residual ratio to form a normal residual ratio that obeys a standardized normal distribution, including: Calculate the mean and standard deviation of each detrended water-heat flux residual ratio; Substitute the mean, standard deviation, and each detrended water and heat flux residual ratio into the preset z-score expression to calculate the normal residual ratio corresponding to each period; The z-score expression is as follows: In the formula, SWER i is the normal residual ratio corresponding to the i-th period; WER′ i is the detrended water heat flux residual ratio corresponding to the ith period; μ is the mean of each detrended water heat flux residual ratio; σ is the standard deviation of each detrended water heat flux residual ratio.

7. A drought-flood rapid transition event recognition device, characterized in that: include: An acquisition module is used to obtain the hydrological and meteorological sequences and groundwater resource changes corresponding to the continuous periods of the target area; A calculation module is used to calculate the available water flux and atmospheric evaporation demand corresponding to each period based on the hydrological and meteorological sequence and groundwater resource changes of the period, wherein the hydrological and meteorological sequence includes the actual regional evapotranspiration of the period; and calculate the water and heat flux residual ratio of the period by using a preset residual ratio expression, the actual regional evapotranspiration of the period, the available water flux and the atmospheric evaporation demand; wherein the residual ratio expression is as follows: In the formula, WER i is the residual ratio of water heat flux in the ith cycle; WA i is the available water flux in the ith cycle; ED i is the atmospheric evaporation demand corresponding to the i-th period; ET i is the actual regional evapotranspiration in the i-th cycle; A construction module for constructing a normal distribution curve based on the ratio of each water heat flux residual; A determination module, used for determining a data point matching each water heat flux residual ratio on the normal distribution curve; A comparison module, for determining that a drought disaster occurs in the target area in the corresponding period when the probability of the corresponding data point in any period is less than N%; for determining that a waterlogging disaster occurs in the target area in the corresponding period when the probability of the corresponding data point in any period is greater than M%, wherein N is less than M; The judgment module is used to judge whether the period interval between drought and flood in the target area is less than the period threshold. If so, it is determined that a drought-flood sudden transition event occurs in the target area; if not, it is determined that no drought-flood sudden transition event occurs in the target area.

8. A drought-flood rapid change event recognition device, characterized in that: including memory and processor; The memory is used to store programs; The processor is used to execute the program to implement each step of the drought-flood rapid transition event identification method as described in any one of claims 1-6.

9. A readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, each step of the method for identifying drought-flood transition events according to any one of claims 1 to 6 is implemented.

Images