Water resource space-time balance evaluation method and system, electronic equipment and storage medium

Abstract

The application discloses a kind of water resource space-time balanced evaluation method, system, electronic equipment and storage medium, method includes: obtaining the water resource supply-demand basic attribute data of the region to be evaluated, to obtain two-dimensional matrix data of space-time according to the water resource supply-demand basic attribute data;According to the two-dimensional matrix data of space-time, obtain the preset Gini coefficient;According to the two-dimensional matrix data of space-time, obtain the preset Theil index;According to the preset Gini coefficient and the preset Theil index, balanced degree is evaluated, and the space-time balanced degree evaluation result is obtained.The embodiment of the application can effectively reflect the difference between different regions of water resources, and can reflect the time characteristics of water resource changes, so as to realize more comprehensive water resource space-time balanced evaluation.The application can be widely applied to the field of water resource evaluation technology.

Description

Technical Field

[0001] This application relates to the field of water resources assessment technology, and in particular to a method, system, electronic device and storage medium for assessing the spatiotemporal equilibrium of water resources. Background Technology

[0002] Water resources are a crucial foundational and strategic resource for economic and social development. While related technologies can analyze static spatial distribution in water resource spatial equilibrium analysis, they struggle to capture dynamic changes within or between years. For instance, the abundance of total water resources in southern regions can mask seasonal and interannual shortages, making a comprehensive assessment of water resource equilibrium difficult.

[0003] In summary, the technical problems existing in the relevant technologies need to be improved. Summary of the Invention

[0004] The main objective of this application is to propose a method, system, electronic device, and storage medium for assessing the spatiotemporal balance of water resources, which can effectively reflect the differences in water resources between different regions and reflect the temporal characteristics of water resource changes, thereby achieving a more comprehensive assessment of the spatiotemporal balance of water resources.

[0005] To achieve the above objectives, one aspect of this application proposes a method for assessing the spatiotemporal equilibrium of water resources, the method comprising the following steps:

[0006] Obtain basic water resource supply and demand attribute data for the area to be evaluated, and construct a spatiotemporal two-dimensional matrix data based on the basic water resource supply and demand attribute data;

[0007] The preset Gini coefficient is calculated based on the spatiotemporal two-dimensional matrix data;

[0008] The preset Theil index is calculated based on the spatiotemporal two-dimensional matrix data;

[0009] The spatiotemporal equilibrium evaluation result is obtained by evaluating the equilibrium degree based on the preset Gini coefficient and the preset Theil index.

[0010] In some embodiments, acquiring basic water resource supply and demand attribute data for the area to be evaluated, and constructing a spatiotemporal two-dimensional matrix data based on the basic water resource supply and demand attribute data, includes:

[0011] The region to be evaluated is divided according to a preset region division rule to obtain several preset sub-regions;

[0012] Obtain the basic water resource supply and demand attribute data for each of the preset sub-regions within a preset time period; wherein, the basic water resource supply and demand attribute data includes water resource quantity data and water demand basic attribute data;

[0013] The spatiotemporal two-dimensional matrix data is constructed based on the preset sub-regions and the basic attribute data of water resource supply and demand.

[0014] In some embodiments, calculating the preset Gini coefficient based on the spatiotemporal two-dimensional matrix data includes:

[0015] Construct a spatial distance matrix based on the spatial distance data of the preset sub-regions;

[0016] Construct a time distance matrix based on the time distance data of the basic attributes of water resource supply and demand;

[0017] The spatiotemporal distance weight matrix is ​​calculated based on the spatial distance matrix and the temporal distance matrix, and the distance-weighted two-dimensional Gini coefficient is calculated based on the spatiotemporal distance weight matrix and the spatiotemporal two-dimensional matrix data.

[0018] In some embodiments, calculating the preset Gini coefficient based on the spatiotemporal two-dimensional matrix data further includes:

[0019] A preset block window is constructed to traverse the spatiotemporal two-dimensional matrix data according to the preset block window to obtain several sub-blocks; wherein, the preset block window includes a spatial dimension window and a temporal dimension window;

[0020] Calculate the preset Gini coefficient corresponding to each sub-block;

[0021] Data visualization is performed based on the preset Gini coefficient to generate a Gini coefficient surface plot.

[0022] In some embodiments, calculating the preset Theil index based on the spatiotemporal two-dimensional matrix data includes:

[0023] The spatiotemporal two-dimensional matrix data is grouped according to preset partitioning conditions to obtain several groups of preset feature data; wherein, the preset partitioning conditions include time conditions or spatial conditions.

[0024] Preset inter-group imbalance data are calculated based on the preset feature data;

[0025] The imbalance data within the preset group is calculated based on the preset feature data;

[0026] The spatiotemporal Theil index is calculated based on the preset inter-group imbalance data and the preset intra-group imbalance data.

[0027] In some embodiments, the step of evaluating the spatiotemporal equilibrium based on the preset Gini coefficient and the preset Theil index to obtain a spatiotemporal equilibrium evaluation result includes:

[0028] The first balance evaluation data is determined by querying the preset Gini value evaluation table based on the preset Gini coefficient.

[0029] Based on the preset Theil index, query the preset Theil index evaluation table to determine the second equilibrium evaluation data;

[0030] The spatiotemporal balance evaluation result is obtained by analyzing the first balance evaluation data and the second balance evaluation data.

[0031] In some embodiments, the preset Theil index includes a first Theil index and a second Theil index; the first Theil index is calculated by grouping the spatiotemporal two-dimensional matrix data according to time conditions, and the second Theil index is calculated by grouping the spatiotemporal two-dimensional matrix data according to spatial conditions.

[0032] The step of querying the preset Theil index evaluation table according to the preset Theil index to determine the second equilibrium evaluation data includes:

[0033] The overall balance is assessed based on the first Theil index and the second Theil index to obtain the third balance evaluation data.

[0034] Based on the first Telegraph index, the preset Telegraph index evaluation table is queried, and the balance analysis is performed based on the first evaluation level data obtained from the query to obtain the fourth balance evaluation data.

[0035] The preset Theil index evaluation table is queried according to the second Theil index, and the balance analysis is performed based on the second evaluation level data obtained from the query to obtain the fifth balance evaluation data.

[0036] The second balance evaluation data is determined based on the third balance evaluation data, the fourth balance evaluation data, and the fifth balance evaluation data.

[0037] To achieve the above objectives, another aspect of this application proposes a water resource spatiotemporal equilibrium assessment system, the system comprising:

[0038] The first module is used to acquire basic water resource supply and demand attribute data of the area to be evaluated, so as to construct a spatiotemporal two-dimensional matrix data based on the basic water resource supply and demand attribute data;

[0039] The second module is used to calculate a preset Gini coefficient based on the spatiotemporal two-dimensional matrix data;

[0040] The third module is used to calculate a preset Theil index based on the spatiotemporal two-dimensional matrix data;

[0041] The fourth module is used to perform a balance evaluation based on the preset Gini coefficient and the preset Theil index to obtain a spatiotemporal balance evaluation result.

[0042] To achieve the above objectives, another aspect of this application provides an electronic device, the electronic device comprising:

[0043] At least one processor;

[0044] At least one memory for storing at least one program;

[0045] When the at least one program is executed by the at least one processor, the at least one processor performs the method described above.

[0046] To achieve the above objectives, another aspect of the embodiments of this application proposes a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method.

[0047] The embodiments of this application include at least the following beneficial effects: This application provides a method, system, electronic device, and storage medium for assessing the spatiotemporal equilibrium of water resources. This scheme acquires basic water resource supply and demand attribute data for the area to be assessed, constructs a spatiotemporal two-dimensional matrix data based on this data, calculates a preset Gini coefficient and a preset Theil index based on the spatiotemporal two-dimensional matrix data, and then assesses the equilibrium degree based on the preset Gini coefficient and the preset Theil index to obtain the spatiotemporal equilibrium degree assessment result, thus achieving a spatiotemporal equilibrium assessment of water resources. It is readily understood that by constructing a spatiotemporal two-dimensional matrix data and then assessing the equilibrium degree based on the preset Gini coefficient and the preset Theil index calculated from the spatiotemporal two-dimensional matrix data, the dimensions of water resource equilibrium assessment are expanded from one-dimensional space to two-dimensional space-time. This can reflect the spatial and temporal imbalance of water resources, effectively reflect the differences in water resources between different regions, and reflect the temporal characteristics of water resource changes, thereby achieving a more comprehensive spatiotemporal equilibrium assessment of water resources. Attached Figure Description

[0048] Figure 1 This is a flowchart of the spatiotemporal equilibrium assessment method for water resources provided in an embodiment of the present invention;

[0049] Figure 2 This is a Gini coefficient surface plot provided in an embodiment of the present invention;

[0050] Figure 3 This is a histogram of the spatiotemporal decomposition of the Theil index provided in this embodiment of the invention;

[0051] Figure 4 This is a schematic diagram of the structure of the water resource spatiotemporal equilibrium assessment system provided in an embodiment of the present invention;

[0052] Figure 5 This is a schematic diagram of the hardware structure of the electronic device provided in an embodiment of the present invention. Detailed Implementation

[0053] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of this application and are not intended to limit it. In the following description, when referring to the accompanying drawings, unless otherwise indicated, the same numbers in different drawings represent the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with those of this application; they are merely examples of apparatuses and methods consistent with some aspects of the embodiments of this application as detailed in the appended claims.

[0054] It is understood that the terms “first,” “second,” etc., used in this application may be used herein to describe various concepts, but unless otherwise stated, these concepts are not limited by these terms. These terms are only used to distinguish one concept from another. For example, without departing from the scope of the embodiments of this application, first information may also be referred to as second information, and similarly, second information may also be referred to as first information. Depending on the context, the words “if,” “when,” or “in response to a determination” as used herein may be interpreted as “when…” or “when…” or “in response to a determination.”

[0055] As used in this application, the terms "at least one", "multiple", "each", "any", etc., "at least one" includes one, two or more, "multiple" includes two or more, "each" refers to each of the corresponding multiples, and "any" refers to any one of the multiples.

[0056] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used herein is for the purpose of describing embodiments of this application only and is not intended to limit this application.

[0057] Before providing a detailed description of the embodiments of this application, some of the nouns and terms involved in the embodiments of this application will be explained first. The nouns and terms involved in the embodiments of this application are subject to the following interpretations.

[0058] Water resources are a crucial foundational and strategic resource for economic and social development. In related technologies, spatial equilibrium analysis of water resources is typically conducted based on the total regional water resources (such as total natural water resources and total available water resources) or on indicators calculated from the total water resources and socio-economic output (such as per capita water resources and water resources per unit of GDP). For example, while traditional indicators like the Gini coefficient are often used to analyze spatial imbalances, they only reflect static spatial distribution and cannot capture dynamic changes within or between years. This is especially true in southern regions where abundant total water resources can mask seasonal and interannual shortages.

[0059] In view of this, this application provides a method, system, electronic device, and storage medium for assessing the spatiotemporal equilibrium of water resources. This solution obtains basic water resource supply and demand attribute data of the area to be assessed, constructs a spatiotemporal two-dimensional matrix data based on the basic water resource supply and demand attribute data, calculates a preset Gini coefficient and a preset Theil index based on the spatiotemporal two-dimensional matrix data, and conducts an equilibrium assessment based on the preset Gini coefficient and the preset Theil index to obtain the spatiotemporal equilibrium assessment result. This can effectively reflect the differences in water resources between different areas and reflect the temporal characteristics of water resource changes, thereby achieving a more comprehensive assessment of the spatiotemporal equilibrium of water resources.

[0060] The water resource spatiotemporal equilibrium assessment method provided in this application relates to the field of water resource assessment technology. This method can be applied to a terminal, a server, or software running on either a terminal or a server. In some embodiments, the terminal can be a smartphone, tablet, laptop, desktop computer, smart speaker, smartwatch, or vehicle terminal, but is not limited to these. The server can be configured as an independent physical server, a server cluster or distributed system composed of multiple physical servers, or a cloud server providing basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communication, middleware services, domain name services, security services, CDN, and big data and artificial intelligence platforms. The server can also be a node server in a blockchain network. The software can be an application implementing the water resource spatiotemporal equilibrium assessment method, but is not limited to the above forms.

[0061] This application can be used in a wide variety of general-purpose or special-purpose computer system environments or configurations. Examples include: personal computers, server computers, handheld or portable devices, tablet devices, multiprocessor systems, microprocessor-based systems, set-top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, and distributed computing environments including any of the above systems or devices. This application can be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, etc., that perform specific tasks or implement specific abstract data types. This application can also be practiced in distributed computing environments where tasks are performed by remote processing devices connected via a communication network. In distributed computing environments, program modules can reside in local and remote computer storage media, including storage devices.

[0062] Figure 1 This is an optional flowchart of the water resource spatiotemporal equilibrium assessment method provided in the embodiments of this application. Figure 1 The method may include, but is not limited to, steps S110 to S140.

[0063] Step S110: Obtain basic water resource supply and demand attribute data for the area to be evaluated, and construct a spatiotemporal two-dimensional matrix data based on the basic water resource supply and demand attribute data.

[0064] Step S120: Calculate the preset Gini coefficient based on the spatiotemporal two-dimensional matrix data.

[0065] Step S130: Calculate the preset Theil index based on the spatiotemporal two-dimensional matrix data.

[0066] Step S140: Equilibrium assessment is performed based on the preset Gini coefficient and preset Theil index to obtain the spatiotemporal equilibrium assessment result.

[0067] In this specific embodiment, the present invention first acquires basic water resource supply and demand attribute data for the area to be evaluated, and constructs a spatiotemporal two-dimensional matrix data based on this data. Specifically, the area to be evaluated in this embodiment refers to an area requiring water resource balance assessment, such as a province or river basin. Accordingly, the present invention acquires the preset basic water resource supply and demand attribute data for the area to be evaluated, such as the basic water resource supply and demand attribute data for the area in different time periods (e.g., month, quarter, or year). Here, the basic water resource supply and demand attribute data in this embodiment refers to water resource supply and demand related data for the relevant area, such as hydrological and meteorological data and socio-economic data. Accordingly, the present invention constructs a spatiotemporal two-dimensional matrix data by statistically analyzing the collected basic water resource supply and demand attribute data for the area to be evaluated in different time periods. Next, the present invention calculates a preset Gini coefficient and a preset Theil index based on the spatiotemporal two-dimensional matrix data. Specifically, the present invention uses the Gini coefficient and the Theil index as indicators for assessing the spatiotemporal balance of water resources. Accordingly, the present invention calculates the preset Gini coefficient and the preset Theil index based on the eigenvalues ​​in the spatiotemporal two-dimensional matrix data. Finally, this embodiment of the invention conducts a balance assessment based on a preset Gini coefficient and a preset Theil index to obtain a spatiotemporal balance assessment result. Specifically, by analyzing the preset Gini coefficient and the preset Theil index, this embodiment of the invention can effectively reflect the differences in water resources between different regions, and can also reflect the temporal characteristics of seasonal and interannual variations. This allows for a more scientific and comprehensive quantitative analysis of the water resource supply and demand pattern, providing strong decision-making support for the optimal allocation and precise regulation of regional water resources.

[0068] In some embodiments of the present invention, basic water resource supply and demand attribute data of the area to be evaluated are obtained, and spatiotemporal two-dimensional matrix data are constructed based on the basic water resource supply and demand attribute data, including but not limited to the following steps:

[0069] The area to be evaluated is divided according to the preset area division rules, resulting in several preset sub-areas.

[0070] Acquire basic water resource supply and demand attribute data for each preset sub-region within a preset time period. This data includes water resource quantity data and water demand attribute data.

[0071] Spatiotemporal two-dimensional matrix data is constructed based on preset sub-regions and basic attribute data of water resource supply and demand.

[0072] In this specific embodiment, the present invention first divides the area to be evaluated according to a preset regional division rule to obtain several preset sub-regions. Specifically, the preset regional division rule in the present invention refers to pre-set regional division conditions, such as according to administrative regions (e.g., cities, counties) or hydrological units (e.g., sub-basins). Accordingly, the present invention divides the area to be evaluated into several preset sub-regions according to the determined division method. Next, the present invention obtains basic water resource supply and demand attribute data for each preset sub-region within a preset time period. Specifically, the basic water resource supply and demand attribute data in the present invention includes water resource quantity data and water demand basic attribute data. The present invention collects the natural water resource quantity or usable water resource quantity of each sub-region within a certain time period T (e.g., month, quarter, or year) to obtain water resource quantity data, and simultaneously collects data reflecting water demand (water demand basic attribute data), such as the population, economic scale (GDP), and agricultural irrigation area of ​​each sub-region. Further, the present invention constructs a spatiotemporal two-dimensional matrix data based on the preset sub-regions and the basic water resource supply and demand attribute data. Specifically, in this embodiment of the invention, a spatiotemporal two-dimensional matrix is ​​obtained by statistically analyzing the water resources and water demand of each sub-region at different time periods, as shown in the following equation (1):

[0073] W i,t i∈{1,2,…,N},t∈{1,2,…,T}(1)

[0074] Wherein, W i,t W represents the characteristic value of the i-th region in the t-th time period. Accordingly, this characteristic value can be selected from different indicators depending on the assessment object; for example, if the assessment object is water resources, then W... i,t It can be characterized by natural water resources or usable water resources (such as water resources after deducting ecological baseflow or water resources whose quality does not meet utilization requirements). For example, if the assessment object is water pressure or water demand, then W i,t Water resources can be characterized by per capita water resources, available water resources per unit of GDP, and available water resources per mu of farmland irrigation. Accordingly, in the spatiotemporal two-dimensional matrix constructed in this embodiment of the invention, rows represent sub-regions, and columns represent time series data of each sub-region.

[0075] It should be noted that the embodiments of the present invention construct the spatiotemporal balance evaluation index based on a holistic approach: both "sub-regions" and "time periods" are regarded as basic units, and each unit corresponds to one of the aforementioned feature values ​​W. i,t Theoretically speaking, the more "average" the eigenvalues ​​are among these units, the more balanced they are in both space and time; conversely, the greater the difference, the more significant the space-time concentration.

[0076] In some embodiments of the present invention, a preset Gini coefficient is calculated based on spatiotemporal two-dimensional matrix data, including but not limited to the following steps:

[0077] A spatial distance matrix is ​​constructed based on the spatial distance data of the preset sub-regions.

[0078] A time distance matrix is ​​constructed based on the time distance data of basic attributes of water resource supply and demand.

[0079] The spatiotemporal distance weight matrix is ​​calculated based on the spatial distance matrix and the temporal distance matrix, and the distance-weighted two-dimensional Gini coefficient is calculated based on the spatiotemporal distance weight matrix and the spatiotemporal two-dimensional matrix data.

[0080] In this specific embodiment, the present invention first constructs a spatial distance matrix based on the spatial distance data of a preset sub-region, and constructs a temporal distance matrix based on the temporal distance data of the basic attributes of water resource supply and demand. Then, a spatiotemporal distance weight matrix is ​​calculated based on the spatial distance matrix and the temporal distance matrix. Specifically, in this embodiment, the spatial distance function and the temporal distance function are selected as linear functions, as shown in equations (2) and (3) below:

[0081] f space (d ij )=1-d ij / d max (2)

[0082] f time (d ts )=1-d ts / t max (3)

[0083] Wherein, d max The distance between two sub-regions within all study areas; t max The difference between two furthest moments in time; d ij It is the spatial distance between subregion i and subregion j; d ts It is the time distance (monthly difference, quarterly difference, etc.) between time period t and time period s.

[0084] Accordingly, in this embodiment of the invention, the corresponding spatial distance matrix and temporal distance matrix are constructed using the aforementioned spatial distance function and temporal distance function, respectively. Further, in this embodiment of the invention, the spatiotemporal distance weight matrix is ​​defined as the spatial distance function f. space (d ij ) and time distance function f time (d ts The product of ) is shown in equation (4) below:

[0085] δ(i,j,t,s)=f space(d ij )×f time (d ts (4)

[0086] Therefore, in this embodiment of the invention, the spatiotemporal distance weight matrix is ​​calculated by multiplying the spatial distance matrix with the temporal distance matrix.

[0087] Furthermore, in this embodiment of the invention, a distance-weighted two-dimensional Gini coefficient is calculated based on the spatiotemporal distance weight matrix and the spatiotemporal two-dimensional matrix data. Specifically, in this embodiment of the invention, the preset Gini coefficient includes a distance-weighted two-dimensional Gini coefficient, as shown in the following equation (5):

[0088]

[0089] Wherein, δ(i,j,t,s) is the spatiotemporal distance weighting function; is the overall mean; j and i are both sub-region labels, i,j∈{1,2,…,N}; t and s are both time period labels, t,s∈{1,2,…,T}.

[0090] In some embodiments of the present invention, the calculation of a preset Gini coefficient based on spatiotemporal two-dimensional matrix data further includes, but is not limited to, the following steps:

[0091] A preset block window is constructed to traverse the spatiotemporal two-dimensional matrix data according to the preset block window, resulting in several sub-blocks. The preset block window includes a spatial dimension window and a temporal dimension window.

[0092] Calculate the preset Gini coefficient for each sub-block.

[0093] Data visualization is performed based on the preset Gini coefficient, generating a Gini coefficient surface plot.

[0094] In this specific embodiment, the present invention first constructs a preset block window, and then traverses the spatiotemporal two-dimensional matrix data according to the preset block window to obtain several sub-blocks, and then calculates the preset Gini coefficient corresponding to each sub-block. Specifically, in order to more intuitively display the changing trend of water resource imbalance at different times and spatial locations, the present invention first sets a preset block window, including a spatial dimension window and a time dimension window. For example, the present invention sets a spatial dimension window Δi and a time dimension window Δt for the block, then each block is a matrix of size Δi×Δt. Then, the present invention traverses all possible sub-blocks W on the entire data matrix. (i,t) =[W i:i+Δi W t:t+Δt ], and calculate W for each sub-block (i,t) Gini coefficient G (i,t)Furthermore, in this embodiment of the invention, data visualization is performed based on the calculated preset Gini coefficient to generate a Gini coefficient surface plot. Specifically, in this embodiment of the invention, three-dimensional data (i, t, G) is plotted as a surface plot to obtain the Gini coefficient surface plot, as shown below. Figure 2 As shown. Among them, Figure 2 The time axis represents the alternating characteristics of water resources in terms of abundance and scarcity over an annual cycle; the spatial axis represents the structural differences in water resource endowments among different sub-regions; the color changes from blue to yellow, corresponding to Gini coefficient values ​​from low to high, with a higher Gini coefficient indicating a lower spatiotemporal water resource balance in that area (i.e., a sub-region within a certain time period). It is easy to understand that this embodiment of the invention, through this method, can further identify hotspots of imbalance within specific time periods and spatial regions, based on the analysis of the overall spatiotemporal water resource balance, and is more conducive to revealing the structural characteristics and evolutionary patterns of spatiotemporal imbalances.

[0095] In some embodiments of the present invention, a preset Theil index is calculated based on spatiotemporal two-dimensional matrix data, including but not limited to the following steps:

[0096] The spatiotemporal two-dimensional matrix data is grouped according to preset partitioning conditions to obtain several sets of preset feature data. The preset partitioning conditions include temporal conditions or spatial conditions.

[0097] Preset inter-group imbalance data are calculated based on preset feature data.

[0098] The unbalance data within the preset group is calculated based on the preset feature data.

[0099] The spatiotemporal Theil index is calculated based on preset inter-group imbalance data and preset intra-group imbalance data.

[0100] In this specific embodiment, the present invention first groups the spatiotemporal two-dimensional matrix data according to preset partitioning conditions to obtain several groups of preset feature data. Then, based on the preset feature data, preset inter-group imbalance data and preset intra-group imbalance data are calculated, and then the spatiotemporal Theil index is calculated. Specifically, the preset partitioning conditions in the present invention include time conditions and spatial conditions. For example, the present invention first groups the spatiotemporal two-dimensional matrix data according to time conditions, and then partitions the data within each group according to spatial conditions, that is, distinguishing spaces within each group. For example, the present invention groups the data by time, that is, each time period t is regarded as a group, then there are T groups in total. Each time period corresponds to a 1×N matrix, and T time periods will result in T 1×N matrices. For example, based on the water resource quantity data of 21 cities in a certain province from 1960 to 2024, the spatiotemporal balance of water resources is analyzed, then there is a group of water resource quantity data for each year of the 21 cities. Correspondingly, the t-th group contains all regions i = 1, 2, ..., N, then the total number of feature values ​​W of the group is... t with the mean As shown in equations (6) and (7) respectively:

[0101]

[0102]

[0103] Accordingly, embodiments of the present invention calculate the inter-group imbalance. That is, it measures the imbalance in the time dimension, reflecting the contribution of the difference in average water volume in different time periods to the overall imbalance, as shown in the following formula (8):

[0104]

[0105] Next, the embodiment of the present invention calculates the intra-group imbalance. That is, it measures the unevenness of spatial dimensions, reflecting the contribution of the distribution differences among regions within the same time period, as shown in the following equation (9):

[0106]

[0107] Accordingly, in this embodiment of the invention, the spatiotemporal Theil index, i.e., the overall imbalance T, is calculated based on the calculated preset inter-group imbalance data and preset intra-group imbalance data. total As shown in equation (10):

[0108]

[0109] Alternatively, embodiments of the present invention can first group the spatiotemporal two-dimensional matrix data according to spatial conditions, and then further divide the data within each group according to temporal conditions, i.e., distinguish time periods within each group. For example, in embodiments of the present invention, spatial grouping is used, where each sub-region i is considered as a group, resulting in a total of N groups. The i-th group contains all time periods t = 1, 2, ..., T, and the total number of eigenvalues ​​W for this group is... i with the mean As shown in equations (11) and (12) respectively:

[0110]

[0111]

[0112] Next, the embodiments of the present invention calculate the inter-group imbalance. That is, it measures the imbalance in spatial dimensions, reflecting the contribution of the difference in the overall average value of different regions to the imbalance, as shown in the following equation (13):

[0113]

[0114] Then, the embodiment of the present invention calculates the intra-group imbalance. That is, it measures the imbalance in the time dimension, reflecting the unequal contribution of time allocation within each region, as shown in the following equation (14):

[0115]

[0116] Accordingly, embodiments of the present invention calculate the overall imbalance T. total As shown in equation (15):

[0117]

[0118] In some embodiments of the present invention, a spatiotemporal equilibrium evaluation result is obtained by performing an equilibrium assessment based on a preset Gini coefficient and a preset Theil index, including but not limited to the following steps:

[0119] The first equilibrium evaluation data is determined by querying the preset Gini value evaluation table based on the preset Gini coefficient.

[0120] Based on the preset Theil index, consult the preset Theil index evaluation table to determine the second equilibrium evaluation data.

[0121] Based on the analysis of the first and second equilibrium evaluation data, the spatiotemporal equilibrium evaluation results are obtained.

[0122] In this specific embodiment, the present invention evaluates and analyzes the spatiotemporal balance of water resources based on two water resource spatiotemporal balance assessment indicators (preset Gini coefficient and preset Theil index). First, the present invention determines the first balance assessment data by querying the preset Gini value evaluation table. Specifically, the present invention first sets a Gini value range to assess the spatiotemporal balance, as shown in Table 1 below:

[0123] Table 1

[0124] Gini value range rating level 0.00–0.15 Highly balanced 0.15–0.30 higher equilibrium 0.30–0.50 Moderate equilibrium 0.50–0.70 Obviously unbalanced >0.70 Severe imbalance

[0125] Accordingly, in the embodiments of the present invention, a preset Gini coefficient (distance-weighted two-dimensional Gini coefficient G) is used. 2D,δ The value of ) is within the range of [0, 1]. The closer the value is to 0, the more balanced the spatiotemporal distribution of water resources; the closer the value is to 1, the more uneven the spatiotemporal distribution of water resources. Among them, the distance-weighted two-dimensional Gini coefficient in Table 1 is mainly used to measure the degree of unevenness in the distribution of water resources in the entire spatiotemporal matrix.

[0126] Next, in this embodiment of the invention, the preset Theil index is queried from a preset Theil index evaluation table to determine the second equilibrium evaluation data. Specifically, in this embodiment of the invention, the preset Theil index T... totalThe value range is within [0, ln(N×T)]. The closer the value is to 0, the more balanced the spatial and temporal distribution of water resources; the larger the value, the more uneven the spatial and temporal distribution of water resources. In order to facilitate decision analysis, the embodiment of the present invention first standardizes the preset Theil index to the range of [0, 1], as shown in the following formula (16):

[0127]

[0128] Where, T total,norm For the preset Theil index T total The standardized value.

[0129] Similarly, embodiments of the present invention can also standardize the indices after spatiotemporal decomposition, so as to... and For example, see equations (17) and (18) below:

[0130]

[0131]

[0132] Accordingly, the standardized preset Theil index is calculated as shown in equation (19):

[0133]

[0134] Where, in the formula and They are respectively and The standardized value.

[0135] Meanwhile, this embodiment of the invention constructs a preset Theil index evaluation table and sets a Theil index value range to evaluate the spatiotemporal equilibrium, as shown in Table 2 below:

[0136] Table 2

[0137]

[0138]

[0139] Table 2 above measures the degree of uneven distribution of water resources across the entire spatiotemporal matrix. Accordingly, embodiments of the present invention can assess the degree of unevenness of water resources in time and space based on the calculation methods of the two preset Theil indices described above. For example, if the unevenness between time groups... or regional imbalance A larger value indicates a significant difference in water resources between different months / quarters / years; if the imbalance is within a time group... or regional group imbalance A larger value indicates a significant difference in the total amount of water resources between regions.

[0140] Furthermore, this embodiment of the invention analyzes the first and second equilibrium evaluation data obtained from the assessment to obtain the spatiotemporal equilibrium evaluation result. Specifically, this embodiment of the invention combines the first equilibrium evaluation data obtained from the assessment analysis using a preset Gini coefficient and the second equilibrium evaluation data obtained from the assessment analysis using a preset Theil index to perform equilibrium analysis, thereby enabling a more scientific and comprehensive quantitative analysis of the water resource supply and demand pattern, providing strong decision-making support for the optimal allocation and precise regulation of regional water resources.

[0141] In some embodiments of the present invention, the preset Theil index includes a first Theil index and a second Theil index. The first Theil index is calculated by grouping spatiotemporal two-dimensional matrix data according to time conditions, i.e., by first dividing the time conditions and then spatially dividing within each group. The second Theil index is calculated by grouping spatiotemporal two-dimensional matrix data according to spatial conditions, i.e., by first dividing the spatial conditions and then temporally dividing within each group. Accordingly, in embodiments of the present invention, the second equilibrium evaluation data is determined by querying a preset Theil index evaluation table, including but not limited to the following steps:

[0142] The overall balance is assessed based on the first and second Theil indices, resulting in the third balance evaluation data.

[0143] Based on the first Theil index, a preset Theil index evaluation table is queried. The balance analysis is performed based on the first evaluation level data obtained from the query to obtain the fourth balance evaluation data.

[0144] The second Theil index is used to query the preset Theil index evaluation table, and the balance analysis is performed based on the second evaluation level data obtained from the query to obtain the fifth balance evaluation data.

[0145] The second balance evaluation data is determined based on the third, fourth, and fifth balance evaluation data.

[0146] In this specific embodiment, the present invention first performs an overall equilibrium assessment based on a first Theil index and a second Theil index to obtain third equilibrium evaluation data. Specifically, the present invention first performs an overall equilibrium assessment by combining preset Theil indices calculated using two spatiotemporal Theil index calculation methods. For example, as... Figure 3 As shown, Figure 3 This demonstrates the temporal (in this case, yearly) and spatial (in this case, city-level) imbalances obtained after further decomposing and calculating the spatiotemporal Theil index for a certain region. Among these, Figure 3The overall height of the two pillars represents the overall spatiotemporal balance of water resources in the region. Group (a) has a value of 0.18, and group (b) has a value of 0.23, both falling within the range of 0.15–0.30. Therefore, the overall spatiotemporal balance of water resources in this region is moderate. Next, in this embodiment of the invention, the first and second Theil indices are queried from a preset Theil index evaluation table to determine the corresponding evaluation levels, namely, the first evaluation level data and the second evaluation level data. Then, a balance analysis is performed to determine the corresponding fourth and fifth balance evaluation data. Specifically, in this embodiment of the invention, the first and second Theil indices are spatiotemporal Theil indices. This embodiment of the invention performs spatiotemporal balance analysis on the first and second Theil indices. For example, as... Figure 3 As shown, in group (a), the degree of imbalance between groups is... A value of 0.06 indicates relatively small differences in water resources across different years in this region; intra-group imbalance. A value of 0.12 indicates significant differences in water resources among different cities within the same year. In group (b), the inter-group imbalance... The value is 0.18, indicating significant differences in the multi-year average water resources among different cities; the degree of imbalance within the group. A value of 0.05 indicates a small difference in water resources across different years within the same city. Finally, this embodiment of the invention determines the second equilibrium evaluation data based on the third, fourth, and fifth equilibrium evaluation data. Specifically, this embodiment of the invention analyzes the corresponding overall equilibrium assessment results and spatiotemporal equilibrium analysis results to determine the equilibrium assessment result based on the spatiotemporal Theil index. For example, through the above-described... Figure 3 Analysis of the central region indicates that the imbalance in water resources in this area mainly stems from structural differences in water resource distribution among different cities. Accordingly, the assessment results of the spatiotemporal balance of water resources obtained from this analysis can generate corresponding decision support for water resource management. For example, a suggested water network engineering layout is: prioritize strengthening the interconnection of the internal water network within the region, and then consider the construction of reservoirs and other water storage projects.

[0147] It should be noted that, in some embodiments of the present invention, after calculating the distance-weighted two-dimensional Gini coefficient and the spatiotemporal Theil index, water resource spatiotemporal equilibrium assessment can be performed using at least one of these two indices. This reflects the differences in water resources between different regions and the temporal characteristics of water resource changes, thus enabling a more comprehensive assessment of water resource spatiotemporal equilibrium. It is easy to understand that related technologies focus on analyzing water resource equilibrium from a spatial perspective, but this is insufficient to comprehensively reflect the characteristics of water resources. For example, it cannot accurately reflect the abundant total water resources in southern regions but their uneven distribution across years and within the year. To address this problem, the embodiments of the present invention propose a water resource equilibrium assessment method that simultaneously considers both spatial and temporal distribution dimensions. This method not only effectively reflects the differences in water resources between different regions but also reflects the temporal characteristics of seasonal and interannual changes. It can more scientifically and comprehensively quantify and analyze the water resource supply and demand pattern in southern regions, providing strong decision-making support for the optimal allocation and precise regulation of regional water resources.

[0148] Please see Figure 4 This application also provides a water resources spatiotemporal equilibrium assessment system, which can implement the above-mentioned water resources spatiotemporal equilibrium assessment method. The system includes:

[0149] The first module 210 is used to acquire basic water resource supply and demand attribute data of the area to be evaluated, so as to construct a spatiotemporal two-dimensional matrix data based on the basic water resource supply and demand attribute data.

[0150] The second module 220 is used to calculate the preset Gini coefficient based on the spatiotemporal two-dimensional matrix data.

[0151] The third module 230 is used to calculate the preset Theil index based on the spatiotemporal two-dimensional matrix data.

[0152] The fourth module 240 is used to evaluate the balance based on the preset Gini coefficient and the preset Theil index, and obtain the spatiotemporal balance evaluation result.

[0153] It is understood that the content of the above method embodiments is applicable to this system embodiment. The specific functions implemented in this system embodiment are the same as those in the above method embodiments, and the beneficial effects achieved are also the same as those achieved in the above method embodiments.

[0154] This application also provides an electronic device, which includes a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the above-described water resource spatiotemporal equilibrium assessment method. This electronic device can be any smart terminal, including tablet computers, in-vehicle computers, etc.

[0155] It is understood that the content of the above method embodiments is applicable to this device embodiment. The specific functions implemented by this device embodiment are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0156] Please see Figure 5 , Figure 5 The hardware structure of an electronic device according to another embodiment is illustrated. The electronic device includes:

[0157] The processor 310 can be implemented using a general-purpose CPU (Central Processing Unit), microprocessor, application-specific integrated circuit (ASIC), or one or more integrated circuits, and is used to execute relevant programs to implement the technical solutions provided in the embodiments of this application.

[0158] The memory 320 can be implemented as a read-only memory (ROM), static storage device, dynamic storage device, or random access memory (RAM). The memory 320 can store the operating system and other applications. When the technical solutions provided in the embodiments of this specification are implemented through software or firmware, the relevant program code is stored in the memory 320 and is called and executed by the processor 310 using the water resource spatiotemporal equilibrium assessment method of the embodiments of this application.

[0159] Input / output interface 330 is used to realize information input and output;

[0160] The communication interface 340 is used to enable communication and interaction between this device and other devices. Communication can be achieved through wired means (such as USB, network cable, etc.) or wireless means (such as mobile network, WIFI, Bluetooth, etc.).

[0161] Bus 350 transmits information between various components of the device (e.g., processor 310, memory 320, input / output interface 330, and communication interface 340);

[0162] The processor 310, memory 320, input / output interface 330 and communication interface 340 are connected to each other within the device via bus 350.

[0163] This application also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for assessing the spatiotemporal equilibrium of water resources.

[0164] It is understood that the content of the above method embodiments is applicable to this storage medium embodiment. The specific functions implemented in this storage medium embodiment are the same as those in the above method embodiments, and the beneficial effects achieved are also the same as those achieved in the above method embodiments.

[0165] Memory, as a non-transitory computer-readable storage medium, can be used to store non-transitory software programs and non-transitory computer-executable programs. Furthermore, memory may include high-speed random access memory, and may also include non-transitory memory, such as at least one disk storage device, flash memory device, or other non-transitory solid-state storage device. In some embodiments, memory may optionally include memory remotely located relative to the processor, and these remote memories can be connected to the processor via a network. Examples of such networks include, but are not limited to, the Internet, intranets, local area networks, mobile communication networks, and combinations thereof.

[0166] The embodiments described in this application are for the purpose of more clearly illustrating the technical solutions of the embodiments of this application, and do not constitute a limitation on the technical solutions provided by the embodiments of this application. As those skilled in the art will know, with the evolution of technology and the emergence of new application scenarios, the technical solutions provided by the embodiments of this application are also applicable to similar technical problems.

[0167] Those skilled in the art will understand that the technical solutions shown in the figures do not constitute a limitation on the embodiments of this application, and may include more or fewer steps than shown, or combine certain steps, or different steps.

[0168] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs.

[0169] Those skilled in the art will understand that all or some of the steps in the methods disclosed above, as well as the functional modules / units in the systems and devices, can be implemented as software, firmware, hardware, or suitable combinations thereof.

[0170] The terms “first,” “second,” “third,” “fourth,” etc. (if present) in the specification and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of this application described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms “comprising” and “having,” and any variations thereof, are intended to cover non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0171] It should be understood that in this application, "at least one (item)" means one or more, and "more than" means two or more. "And / or" is used to describe the relationship between related objects, indicating that three relationships can exist. For example, "A and / or B" can represent three cases: only A exists, only B exists, and both A and B exist simultaneously, where A and B can be singular or plural. The character " / " generally indicates that the preceding and following related objects are in an "or" relationship. "At least one (item) of the following" or similar expressions refer to any combination of these items, including any combination of single or plural items. For example, at least one (item) of a, b, or c can represent: a, b, c, "a and b", "a and c", "b and c", or "a and b and c", where a, b, and c can be single or multiple.

[0172] In the several embodiments provided in this application, it should be understood that the disclosed apparatus and methods can be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative; for instance, the division of the units described above is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be through some interfaces; the indirect coupling or communication connection between apparatuses or units may be electrical, mechanical, or other forms.

[0173] The units described above as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0174] Furthermore, the functional units in the various embodiments of this application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit.

[0175] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes multiple instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of the various embodiments of this application. The aforementioned storage medium includes various media capable of storing programs, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0176] The preferred embodiments of the present application have been described above with reference to the accompanying drawings, but this does not limit the scope of the claims of the present application. Any modifications, equivalent substitutions, and improvements made by those skilled in the art without departing from the scope and substance of the embodiments of the present application shall be within the scope of the claims of the present application.

Claims

1. A method for assessing the spatiotemporal equilibrium of water resources, characterized in that, The method includes the following steps: Obtain basic water resource supply and demand attribute data for the area to be evaluated, and construct a spatiotemporal two-dimensional matrix data based on the basic water resource supply and demand attribute data; The preset Gini coefficient is calculated based on the spatiotemporal two-dimensional matrix data; The preset Theil index is calculated based on the spatiotemporal two-dimensional matrix data; The spatiotemporal equilibrium evaluation result is obtained by performing an equilibrium evaluation based on the preset Gini coefficient and the preset Theil index. The step of calculating the preset Gini coefficient based on the spatiotemporal two-dimensional matrix data includes: Construct a spatial distance matrix based on the spatial distance data of the preset sub-regions; Construct a time distance matrix based on the time distance data of the basic attributes of water resource supply and demand; A spatiotemporal distance weight matrix is ​​calculated based on the spatial distance matrix and the temporal distance matrix, and a distance-weighted two-dimensional Gini coefficient is calculated based on the spatiotemporal distance weight matrix and the spatiotemporal two-dimensional matrix data. The formula for calculating the distance-weighted two-dimensional Gini coefficient is as follows: ； Where, in the formula The spatiotemporal distance weighting function; This is the overall mean; Representing the The region in the first Feature values ​​for each time period; and Both are sub-region markers, ; and Both are time period markers, ; The step of calculating the preset Theil index based on the spatiotemporal two-dimensional matrix data includes: The spatiotemporal two-dimensional matrix data is grouped according to preset partitioning conditions to obtain several groups of preset feature data; wherein, the preset partitioning conditions include time conditions or spatial conditions. Preset inter-group imbalance data are calculated based on the preset feature data; The imbalance data within the preset group is calculated based on the preset feature data; The spatiotemporal Theil index is calculated based on the preset inter-group imbalance data and the preset intra-group imbalance data.

2. The method according to claim 1, characterized in that, The process of acquiring basic water resource supply and demand attribute data for the area to be evaluated, and constructing a spatiotemporal two-dimensional matrix data based on the basic water resource supply and demand attribute data, includes: The region to be evaluated is divided according to a preset region division rule to obtain several preset sub-regions; Obtain the basic water resource supply and demand attribute data for each of the preset sub-regions within a preset time period; wherein, the basic water resource supply and demand attribute data includes water resource quantity data and water demand basic attribute data; The spatiotemporal two-dimensional matrix data is constructed based on the preset sub-regions and the basic attribute data of water resource supply and demand.

3. The method according to claim 1, characterized in that, The step of calculating the preset Gini coefficient based on the spatiotemporal two-dimensional matrix data further includes: A preset block window is constructed to traverse the spatiotemporal two-dimensional matrix data according to the preset block window to obtain several sub-blocks; wherein, the preset block window includes a spatial dimension window and a temporal dimension window; Calculate the preset Gini coefficient corresponding to each sub-block; Data visualization is performed based on the preset Gini coefficient to generate a Gini coefficient surface plot.

4. The method according to claim 1, characterized in that, The process of evaluating the spatiotemporal equilibrium based on the preset Gini coefficient and the preset Theil index to obtain the spatiotemporal equilibrium evaluation result includes: The first balance evaluation data is determined by querying the preset Gini value evaluation table based on the preset Gini coefficient. Based on the preset Theil index, query the preset Theil index evaluation table to determine the second equilibrium evaluation data; The spatiotemporal balance evaluation result is obtained by analyzing the first balance evaluation data and the second balance evaluation data.

5. The method according to claim 4, characterized in that, The preset Theil index includes a first Theil index and a second Theil index; the first Theil index is calculated by grouping the spatiotemporal two-dimensional matrix data according to time conditions, and the second Theil index is calculated by grouping the spatiotemporal two-dimensional matrix data according to spatial conditions. The step of querying the preset Theil index evaluation table according to the preset Theil index to determine the second equilibrium evaluation data includes: The overall balance is assessed based on the first Theil index and the second Theil index to obtain the third balance evaluation data. Based on the first Telegraph index, the preset Telegraph index evaluation table is queried, and the balance analysis is performed based on the first evaluation level data obtained from the query to obtain the fourth balance evaluation data. The preset Theil index evaluation table is queried according to the second Theil index, and the balance analysis is performed based on the second evaluation level data obtained from the query to obtain the fifth balance evaluation data. The second balance evaluation data is determined based on the third balance evaluation data, the fourth balance evaluation data, and the fifth balance evaluation data.

6. A water resources spatiotemporal equilibrium assessment system, characterized in that, The system includes: The first module is used to acquire basic water resource supply and demand attribute data of the area to be evaluated, so as to construct a spatiotemporal two-dimensional matrix data based on the basic water resource supply and demand attribute data; The second module is used to calculate a preset Gini coefficient based on the spatiotemporal two-dimensional matrix data; The third module is used to calculate a preset Theil index based on the spatiotemporal two-dimensional matrix data; The fourth module is used to evaluate the balance based on the preset Gini coefficient and the preset Theil index, and obtain the spatiotemporal balance evaluation result. The step of calculating the preset Gini coefficient based on the spatiotemporal two-dimensional matrix data includes: Construct a spatial distance matrix based on the spatial distance data of the preset sub-regions; Construct a time distance matrix based on the time distance data of the basic attributes of water resource supply and demand; A spatiotemporal distance weight matrix is ​​calculated based on the spatial distance matrix and the temporal distance matrix, and a distance-weighted two-dimensional Gini coefficient is calculated based on the spatiotemporal distance weight matrix and the spatiotemporal two-dimensional matrix data. The formula for calculating the distance-weighted two-dimensional Gini coefficient is as follows: ； Where, in the formula The spatiotemporal distance weighting function; This is the overall mean; Representing the The region in the first Feature values ​​for each time period; With Both are sub-region markers, ; and Both are time period markers, ; The step of calculating the preset Theil index based on the spatiotemporal two-dimensional matrix data includes: The spatiotemporal two-dimensional matrix data is grouped according to preset partitioning conditions to obtain several groups of preset feature data; wherein, the preset partitioning conditions include time conditions or spatial conditions. Preset inter-group imbalance data are calculated based on the preset feature data; The imbalance data within the preset group is calculated based on the preset feature data; The spatiotemporal Theil index is calculated based on the preset inter-group imbalance data and the preset intra-group imbalance data.

7. An electronic device, characterized in that, include: At least one processor; At least one memory for storing at least one program; When the at least one program is executed by the at least one processor, the at least one processor implements the method as described in any one of claims 1-5.

8. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by a processor, it implements the method of any one of claims 1 to 5.

Images