A method and apparatus for water resource security assessment based on water quantity-water quality coupling of supply and demand spatial matching

Abstract

This invention discloses a water resource security assessment method and apparatus based on water quantity-water quality coupling using spatial matching of supply and demand, belonging to the technical field of water resource security assessment and geographic information system spatial analysis. In this method, data preprocessing is performed on supply pixels and demand pixels to make water quantity, water quality, social water demand, and ecological water demand indicators dimensionless; a water resource supply attraction function is constructed to transform water quantity and water quality into unified variables; a distance decay function is constructed to quantify the impact of spatial distance on supply and demand interaction; the probability of supply and demand interaction, the unit service capacity of supply pixels, and the water resource availability of demand pixels are calculated; and a water resource security index is constructed based on the ratio of water resource availability to demand intensity. This invention, through a spatial matching mechanism of supply and demand, couples water quantity, water quality, and spatial distance within a unified framework, achieving refined evaluation at the pixel scale. It can accurately identify areas of water resource supply and demand contradictions, providing scientific support for the optimal allocation of water resources.

Description

Technical Field

[0001] This application relates to the field of water resources security assessment and geographic information system spatial analysis technology, and in particular to a water resources security assessment method and device based on water quantity-water quality coupling of supply and demand spatial matching. Background Technology

[0002] Water security generally refers to the quantity and quality of water resources that a region can sustainably obtain within a certain spatial and temporal scale to meet the needs of socio-economic development and ecosystems. With the acceleration of climate change and urbanization, water security assessment has gradually become an important technical foundation for water resource management and regional sustainable development.

[0003] Currently, various technical solutions have been developed for water resource security assessment. Early studies typically used simple, empirical, single indicators to judge the state of water resource security, such as using per capita water resources below 1700 m³ / a as the critical value for water resource pressure. While these methods are simple and easy to implement, they neglect the spatial and temporal distribution of water resources and water use efficiency, and lack attention to water quality factors. To take into account the multiple socio-economic attributes of water resources, water resource vulnerability indices and water resource scarcity indices have been developed, incorporating data such as population and GDP. Chinese scholar Xia Jun proposed a comprehensive approach to water resource carrying capacity, considering inter-basin water transfer, ecological water use, and socio-economic factors, and evaluating water resource security by comparing supply and demand. Internationally, the UK Eco-Hydrology Centre's Water Poverty Index comprehensively considers the entire process of water resource utilization from five dimensions: water resources themselves, water access, socio-economic capacity, usage, and the water environment. The Asian Development Bank's Water Security Outlook report constructs a water resource security index framework based on five dimensions: household water use, economic feasibility, urban services, river and ecosystem restoration, and flood resilience. Comprehensive indicators are often combined with multi-criteria decision analysis to form an evaluation framework. Methods such as entropy weight method, fuzzy comprehensive evaluation method, and analytic hierarchy process are used to assign weights to indicators. However, the selection of indicators and determination of weights are subjective and lack universal comparability. Moreover, the data are mostly based on statistical reports of administrative units, which cannot achieve fine description at the spatial scale.

[0004] Spatial accessibility analysis is widely used in the site selection research of medical facilities and public service facilities. The Huff model and its improved forms (such as the two-step movement search method and the three-step movement search method) provide technical means for supply and demand spatial matching analysis, but it has not yet been specifically constructed for the characteristics of water resource systems.

[0005] In summary, existing water resource security assessment methods have the following shortcomings: they do not uniformly integrate water quantity, water quality, and spatial distance factors; they do not consider the competitive allocation mechanism between supply and demand, and assume that the supply resources are evenly distributed or fully accessible within the region by using the total quantity comparison method; they are difficult to achieve pixel-scale refined evaluation and cannot reflect the spatial heterogeneity within a city or region.

[0006] Therefore, how to couple and express water quantity, water quality and spatial distance within a unified framework, achieve refined evaluation at the pixel scale, and accurately identify areas with water supply and demand contradictions is a problem that needs to be solved by those skilled in the art. Summary of the Invention

[0007] This application provides a water resource security assessment method and apparatus based on water quantity-water quality coupling for supply and demand spatial matching. It can solve the technical problems of existing water resource security assessment methods, which are unable to integrate water quantity, water quality and spatial distance factors under a unified calculation framework, cannot characterize the competitive distribution relationship between supply and demand, and are difficult to achieve refined assessment.

[0008] The first aspect of this application provides a water resource security assessment method based on water quantity-water quality coupling of supply and demand spatial matching, including:

[0009] Data preprocessing is performed on supply pixels and demand pixels respectively. The water quantity and water quality indicators of supply pixels are made dimensionless, and the social water demand and ecological water demand indicators of demand pixels are made dimensionless. Based on the dimensionless water quantity and water quality indicators, a water resource supply attraction function is constructed for each supply cell. The water resource supply attraction function is used to characterize the comprehensive supply capacity of the supply cell. A distance decay function is constructed to quantify the impact of the spatial distance between demand pixels and supply pixels on the intensity of supply-demand interaction. For each demand cell, the supply and demand interaction probability between the demand cell and each supply cell within the distance threshold is calculated based on the distance decay function and the water resource supply attractiveness. For each supply cell, the unit service capacity of the supply cell is calculated based on the supply-demand interaction probability and the demand intensity of the demand cell. For each demand cell, the water resource availability of the demand cell is calculated based on the supply-demand interaction probability and the unit service capacity of the supply cell. A water resource security index is constructed based on the ratio of water resource availability to the demand intensity of demand pixels.

[0010] Optionally, a distance decay function is constructed, including: A Gaussian distance decay function is adopted, which uses the Euclidean distance between the demand pixel and the supply pixel as the independent variable and a preset distance decay parameter as the decay factor to calculate the decay weight within a preset distance threshold.

[0011] Optionally, the probability of supply-demand interaction between demand pixels and each supply pixel within a distance threshold is calculated, including: Multiplying the attractiveness of water resource supply to the supplied pixel by the distance decay function yields the supply-demand interaction strength value; Divide the supply-demand interaction strength value by the sum of the supply-demand interaction strength values ​​of all supply pixels within the distance threshold to obtain the supply-demand interaction probability of the demand pixel to the supply pixel.

[0012] Optionally, calculating the unit service capacity supplied to a cell includes: The attractiveness of water resource supply to the supplied pixels is used as the molecule; The sum of the products of the demand intensity of demand pixels and the probability of supply-demand interaction within the distance threshold is used as the denominator; The ratio of the numerator to the denominator is used as the unit service capacity supplied to a pixel.

[0013] Optionally, the demand intensity of a demand cell is the sum of the dimensionless social water demand index and the dimensionless ecological water demand index.

[0014] Optionally, the availability of water resources for the required pixels is calculated, including: The supply contribution value is obtained by multiplying the unit service capacity of the supplied pixels by the probability of supply and demand interaction. The water resource availability of demand pixels is obtained by summing the supply contribution values ​​of all supply pixels within the distance threshold.

[0015] Optionally, a water resource security index may be constructed, including: The ratio of water availability to demand intensity of demand pixels is used as the water security index of demand pixels. When the demand intensity is zero, a very small positive number is added to the denominator to prevent division by zero error.

[0016] This application provides a water resource security assessment device based on water quantity-water quality coupling for spatial matching of supply and demand, comprising: The data preprocessing unit is used to preprocess the supply pixels and demand pixels respectively, making the water quantity and water quality indicators of the supply pixels dimensionless, and making the social water demand and ecological water demand indicators of the demand pixels dimensionless. The supply attraction construction unit is used to construct a water resource supply attraction function for each supply cell based on the dimensionless water quantity and water quality indicators. The water resource supply attraction function is used to characterize the comprehensive supply capacity of the supply cell. Distance decay building blocks are used to construct distance decay functions, which are used to quantify the impact of the spatial distance between demand pixels and supply pixels on the intensity of supply-demand interaction. The interaction probability calculation unit is used to calculate the supply and demand interaction probability between each demand pixel and each supply pixel within the distance threshold, based on the distance decay function and the attractiveness of water resource supply. The unit service capacity calculation unit is used to calculate the unit service capacity of each supply cell based on the supply-demand interaction probability and the demand intensity of the demand cell. The water availability calculation unit is used to calculate the water availability of each demand cell based on the supply-demand interaction probability and the unit service capacity of the supply cell. The security index construction unit is used to construct a water resource security index based on the ratio of water resource availability to the demand intensity of demand pixels.

[0017] A third aspect of this application provides a water resource security assessment device based on water quantity-water quality coupling for supply and demand spatial matching, comprising: One or more processors; A memory on which one or more programs are stored; When the one or more programs are executed by the one or more processors, the one or more processors implement the water resource security assessment method based on water quantity-water quality coupling for supply and demand spatial matching as described above.

[0018] The fourth aspect of this application provides a computer storage medium for storing a program, which, when executed, is used to implement the water resource security assessment method based on water quantity-water quality coupling according to any one of the preceding claims.

[0019] Compared with the prior art, the present invention has the following beneficial effects: This invention transforms water quantity and water quality into unified variables by constructing a supply attraction function and introduces a distance decay model, so that water quantity, water quality and water intake distance are coupled and expressed in the same mathematical structure. This solves the problem that traditional methods treat water quantity and water quality as independent evaluation factors and do not consider spatial distance factors.

[0020] This invention calculates the probability of supply and demand interaction and the unit service capacity, so that the supply capacity is affected by the actual demand allocation pressure and the demand acquisition is constrained by supply competition, forming a two-way coupled feedback mechanism between supply and demand. This solves the problem of traditional methods using total quantity comparison and ignoring the competitive allocation mechanism.

[0021] In this invention, all supply and demand variables are expressed in raster form, and the supply and demand interaction is calculated at the pixel level. The output result is a continuous spatial distribution field. Compared with the traditional method that uses administrative regions as statistical units, it can accurately identify differences within regions and provide technical support for the refined management of water resources.

[0022] This invention can be extended to water demand in sectors such as agriculture and industry, and can accommodate various types of demand. The distance attenuation parameter and weighting coefficient can be adjusted or optimized according to regional characteristics to adapt to applications at different watershed and urban scales. Dynamic spatiotemporal evaluation and analysis can be carried out with the support of multi-temporal data. Attached Figure Description

[0023] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments recorded in this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0024] Figure 1 A flowchart illustrating a water resource security assessment method based on water quantity-water quality coupling for spatial matching of supply and demand, provided for an embodiment of this application; Figure 2 A schematic diagram of the structure of a water resource security assessment device based on water quantity-water quality coupling for supply and demand spatial matching provided in this application embodiment; Figure 3 This is a schematic diagram of the structure of a device provided in an embodiment of this application. Detailed Implementation

[0025] This application provides a water resource security assessment method and apparatus based on water quantity-water quality coupling for spatial matching of supply and demand, which is used to identify areas with water supply and demand contradictions and to provide technical support for water resource optimization, urban planning and ecological restoration.

[0026] To address the technical challenges of existing water resource security assessment methods, such as the inability to simultaneously integrate water quantity, water quality, and spatial distance factors within a unified computational framework, the inability to characterize the competitive allocation relationship between supply and demand, and the difficulty in achieving refined assessments, this invention provides a water quantity-water quality coupled water resource security assessment method based on a spatial matching mechanism for supply and demand. By constructing a function of water resource supply attractiveness and demand intensity, and introducing a spatial distance attenuation and probabilistic competitive allocation model, dynamic matching calculation of water resource supply and demand relationships is achieved at the pixel scale. This establishes a unified water resource security index calculation framework, thereby enabling a refined quantitative assessment of the regional water resource security status.

[0027] See Figure 1 This figure is a flowchart illustrating a water resource security assessment method based on water quantity-water quality coupling with spatial matching of supply and demand, provided in an embodiment of this application. The water resource security assessment method based on water quantity-water quality coupling with spatial matching of supply and demand provided in this embodiment can be implemented, for example, through the following steps S101-107.

[0028] S101: Perform data preprocessing on supply pixels and demand pixels respectively.

[0029] In this embodiment of the application, the water quantity and water quality indicators of the supplied pixels are dimensionless, and the social water demand and ecological water demand indicators of the demand pixels are dimensionless.

[0030] Specifically, regarding the supply of pixels Water volume indicators involved Water quality indicators and demand pixels Social water demand indicators involved Ecological water demand indicators Dimensionless processing is performed to transform different physical quantity indices into dimensionless variables. Dimensionless processing methods can include linear normalization, standardization, or other equivalent methods; this invention does not impose specific limitations.

[0031] S102: Construct a water resource supply attraction function for each supply cell based on the dimensionless water quantity and water quality indicators.

[0032] In this embodiment of the application, for the supply point Constructing an attractiveness function for water resource supply : ; in, and These are water quantity indicators and water quality indicators, and These are the water quantity and water quality indicators after dimensionless processing, respectively. and These are the weighting coefficients.

[0033] S103: Construct the distance decay function.

[0034] In this embodiment, the distance decay function is used to quantify the impact of the spatial distance between demand pixels and supply pixels on the intensity of supply-demand interaction. A Gaussian distance decay function is employed, which uses the Euclidean distance between demand pixels and supply pixels as the independent variable and a preset distance decay parameter as the decay factor to calculate the decay weight within a preset distance threshold.

[0035] Specifically, a Gaussian distance decay function is introduced in the supply and demand matching process: ; in, For distance decay weights, For demand points to supply point distance, For distance decay parameters; only at the distance threshold Calculations are performed internally.

[0036] S104: For each demand cell, calculate the supply-demand interaction probability between the demand cell and each supply cell within the distance threshold based on the distance decay function and the water resource supply attractiveness.

[0037] In this embodiment, the water resource supply attractiveness of the supply pixel is multiplied by the distance decay function to obtain the supply-demand interaction strength value; the supply-demand interaction strength value is divided by the sum of the supply-demand interaction strength values ​​of all supply pixels within the distance threshold to obtain the supply-demand interaction probability of the demand pixel to the supply pixel.

[0038] Specifically, for any demand pixel Calculate its value for all points at the distance threshold. Supply points within Interaction probability : ; in, For distance decay weights, For demand points to supply point The distance; The attraction function for water resource supply.

[0039] S105: For each supply cell, calculate the unit service capacity of the supply cell based on the supply-demand interaction probability and the demand intensity of the demand cell.

[0040] In this embodiment, the water resource supply attractiveness of the supply pixel is used as the numerator; the sum of the product of the demand intensity of the demand pixel and the probability of supply-demand interaction within a distance threshold is used as the denominator; and the ratio of the numerator to the denominator is used as the unit service capacity of the supply pixel. The demand intensity of the demand pixel is the sum of the dimensionless social water demand index and the dimensionless ecological water demand index.

[0041] Specifically, regarding supply points Demand is allocated to supply points according to interaction probability, and the unit service capacity of each supply point is calculated. : ; ; in, For the probability of supply and demand interaction; For the attraction function of water resource supply; Distance threshold; and These are the social water demand index and the ecological water demand index after dimensionless processing; To represent the demand pixels Total demand intensity.

[0042] S106: For each demand cell, calculate the water resource availability of the demand cell based on the supply-demand interaction probability and the unit service capacity of the supply cell.

[0043] In this embodiment, the unit service capacity of the supplied pixel is multiplied by the supply-demand interaction probability to obtain the supply contribution value; the supply contribution values ​​of all supplied pixels within the distance threshold are summed to obtain the water resource availability of the demand pixel.

[0044] Specifically, the supply capacity of the supply points Return to the demand pixel Calculate the water resource availability of this pixel. : ; in, This represents the probability of supply and demand interaction.

[0045] S107: Construct a water resource security index based on the ratio of water resource availability to the demand intensity of demand pixels.

[0046] In this embodiment of the application, the ratio of water resource availability to demand intensity of demand pixels is used as the water resource security index of demand pixels. When the demand intensity is zero, a very small positive number is added to the denominator to prevent division by zero error.

[0047] Specifically, the water resource security index is defined based on the supply-demand ratio. : ； in, To prevent extremely small positive numbers with a denominator of zero.

[0048] This invention takes spatial matching of supply and demand as its starting point, abstracting water resource security into the interaction between supply and demand units, fully reflecting the multidimensional attributes of the water resource system. By introducing the concept of spatial accessibility model, demand units allocate resources to multiple supply units, while supply units form an attractiveness competition relationship, creating a two-way feedback mechanism between supply and demand.

[0049] This invention constructs a supply attraction function to transform water quantity and water quality into unified variables, and introduces water intake distance constraints through a distance decay function, so that water quantity and water quality jointly determine the supply attraction, and distance affects the intensity of supply and demand interaction through a continuous decay function. All three are expressed in a coupled mathematical structure.

[0050] This invention differs from traditional evaluation methods based on total statistics or administrative region averages. All supply and demand variables are expressed in raster form, and supply and demand interactions are calculated at the pixel level, resulting in a continuous spatial distribution field.

[0051] Based on the methods provided in the above embodiments, this application also provides a water resource safety assessment device based on water quantity-water quality coupling of supply and demand spatial matching. The following describes the water resource safety assessment device based on water quantity-water quality coupling of supply and demand spatial matching in conjunction with the accompanying drawings.

[0052] See Figure 2 The figure is a schematic diagram of a water resource security assessment device based on water quantity-water quality coupling with supply and demand spatial matching provided in an embodiment of this application.

[0053] The water resource security assessment device 200 based on water quantity-water quality coupling with supply and demand spatial matching provided in this application includes: a data preprocessing unit 201, a supply attraction construction unit 202, a distance attenuation construction unit 203, an interaction probability calculation unit 204, a unit service capacity calculation unit 205, a water resource availability calculation unit 206, and a security index construction unit 207.

[0054] The data preprocessing unit 201 is used to preprocess the supply pixels and demand pixels respectively, making the water quantity and water quality indicators of the supply pixels dimensionless, and making the social water demand and ecological water demand indicators of the demand pixels dimensionless. The supply attraction construction unit 202 is used to construct a water resource supply attraction function for each supply cell based on the dimensionless water quantity and water quality indicators. The water resource supply attraction function is used to characterize the comprehensive supply capacity of the supply cell. Distance attenuation building unit 203 is used to construct a distance attenuation function, which is used to quantify the impact of the spatial distance between demand pixels and supply pixels on the intensity of supply and demand interaction. The interaction probability calculation unit 204 is used to calculate the supply and demand interaction probability between the demand pixel and each supply pixel within the distance threshold for each demand pixel, based on the distance decay function and the water resource supply attraction. Unit service capacity calculation unit 205 is used to calculate the unit service capacity of each supply cell based on the supply-demand interaction probability and the demand intensity of the demand cell. Water availability calculation unit 206 is used to calculate the water availability of each demand cell based on the supply-demand interaction probability and the unit service capacity of the supply cell. The security index construction unit 207 is used to construct a water resource security index based on the ratio of water resource availability to the demand intensity of demand pixels.

[0055] In one possible implementation, the distance attenuation building unit 202 is specifically used for: A Gaussian distance decay function is adopted, which uses the Euclidean distance between the demand pixel and the supply pixel as the independent variable and a preset distance decay parameter as the decay factor to calculate the decay weight within a preset distance threshold.

[0056] In one possible implementation, the interaction probability calculation unit 204 is specifically used for: Multiplying the attractiveness of water resource supply to the supplied pixel by the distance decay function yields the supply-demand interaction strength value; Divide the supply-demand interaction strength value by the sum of the supply-demand interaction strength values ​​of all supply pixels within the distance threshold to obtain the supply-demand interaction probability of the demand pixel to the supply pixel.

[0057] In one possible implementation, the unit service capacity calculation unit 205 is specifically used for: The attractiveness of water resource supply to the supplied pixels is used as the molecule; The sum of the products of the demand intensity of demand pixels and the probability of supply-demand interaction within the distance threshold is used as the denominator; The ratio of the numerator to the denominator is used as the unit service capacity supplied to a pixel.

[0058] In one possible implementation, the demand intensity of a demand cell is the sum of the dimensionless social water demand index and the dimensionless ecological water demand index.

[0059] In one possible implementation, the water resource availability calculation unit 206 is specifically used for: The supply contribution value is obtained by multiplying the unit service capacity of the supplied pixels by the probability of supply and demand interaction. The water resource availability of demand pixels is obtained by summing the supply contribution values ​​of all supply pixels within the distance threshold.

[0060] In one possible implementation, the security index construction unit 207 is specifically used for: The ratio of water availability to demand intensity of demand pixels is used as the water security index of demand pixels. When the demand intensity is zero, a very small positive number is added to the denominator to prevent division by zero error.

[0061] Since the water resource safety assessment device 200 based on supply and demand spatial matching and water quality coupling is the same device as the water resource safety assessment method based on supply and demand spatial matching and water quality coupling provided in the above method embodiments, the specific implementation of each unit of the water resource safety assessment device 200 based on supply and demand spatial matching and water quality coupling is based on the same concept as in the above method embodiments. Therefore, for the specific implementation of each unit of the water resource safety assessment device 200 based on supply and demand spatial matching and water quality coupling, please refer to the description of the water resource safety assessment method based on supply and demand spatial matching and water quality coupling in the above method embodiments, and will not be repeated here.

[0062] This application embodiment also provides a water resource security assessment device based on water quantity-water quality coupling of supply and demand spatial matching, the device including: a processor and a memory; The memory is used to store instructions; The processor is used to execute the instructions in the memory to perform the water resource security assessment method based on water quantity-water quality coupling of supply and demand spatial matching mentioned in the above embodiments.

[0063] It should be noted that the water resource security assessment equipment based on water quantity-water quality coupling for supply and demand spatial matching provided in this application embodiment can all have the following hardware structure: Figure 3 The structure shown, Figure 3 This is a schematic diagram of the structure of a device provided in an embodiment of this application.

[0064] Please see Figure 3 As shown, device 300 includes: a processor 310, a communication interface 320, and a memory 330. The number of processors 310 in device 300 can be one or more. Figure 3 Taking a processor as an example, in this embodiment, the processor 310, communication interface 320, and memory 330 can be connected via a bus system or other means. Figure 3 Taking the connection between China and Israel via bus system 340 as an example.

[0065] Processor 310 may be a central processing unit (CPU), a network processor (NP), or a combination of a CPU and an NP. Processor 310 may further include hardware chips. These hardware chips may be application-specific integrated circuits (ASICs), programmable logic devices (PLDs), or combinations thereof. The PLD may be a complex programmable logic device (CPLD), a field-programmable gate array (FPGA), generic array logic (GAL), or any combination thereof.

[0066] The memory 330 may include volatile memory, such as random-access memory (RAM); the memory 330 may also include non-volatile memory, such as flash memory, hard disk drive (HDD) or solid-state drive (SSD); the memory 330 may also include a combination of the above types of memory.

[0067] Optionally, the memory 330 stores an operating system and programs, executable modules, or data structures, or subsets thereof, or extended sets thereof. The programs may include various operation instructions for implementing various operations. The operating system may include various system programs for implementing various basic services and handling hardware-based tasks. The processor 310 can read the programs in the memory 330 to implement the water resource security assessment method based on water quantity-water quality coupling for supply and demand spatial matching provided in this embodiment.

[0068] The bus system 340 can be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus, etc. The bus system 340 can be divided into address bus, data bus, control bus, etc. For ease of representation, Figure 3The bus is represented by a single thick line, but this does not mean that there is only one bus or one type of bus.

[0069] This application also provides a computer-readable storage medium including instructions that, when executed on a computer, cause the computer to perform the water resource security assessment method based on water quantity-water quality coupling for spatial matching of supply and demand mentioned in the above embodiments.

[0070] This application also provides a computer program product containing instructions that, when run on a computer, cause the computer to execute the water resource security assessment method based on water quantity-water quality coupling and supply-demand spatial matching mentioned in the above embodiments.

[0071] Although the invention has been specifically shown and described in conjunction with preferred embodiments, those skilled in the art should understand that various changes in form and detail may be made to the invention without departing from the spirit and scope of the invention as defined in the appended claims, all of which shall be within the scope of protection of the invention.

Claims

1. A water resource security assessment method based on water quantity-water quality coupling for spatial matching of supply and demand, characterized in that, include: Data preprocessing is performed on supply pixels and demand pixels respectively. The water quantity and water quality indicators of the supply pixels are dimensionless, and the social water demand and ecological water demand indicators of the demand pixels are dimensionless. Based on the dimensionless water quantity and water quality indicators, a water resource supply attraction function is constructed for each supply cell. The water resource supply attraction function is used to characterize the comprehensive supply capacity of the supply cell. A distance decay function is constructed to quantify the impact of the spatial distance between the demand pixel and the supply pixel on the intensity of supply-demand interaction. For each demand cell, the supply-demand interaction probability between the demand cell and each supply cell within the distance threshold is calculated based on the distance decay function and the water resource supply attractiveness. For each supply cell, the unit service capacity of the supply cell is calculated based on the supply-demand interaction probability and the demand intensity of the demand cell. For each demand cell, the water resource availability of the demand cell is calculated based on the supply-demand interaction probability and the unit service capacity of the supply cell. A water resource security index is constructed based on the ratio of water resource availability to the demand intensity of the demand pixels.

2. The method according to claim 1, characterized in that, The construction of the distance decay function includes: A Gaussian distance decay function is adopted, which uses the Euclidean distance between the demand pixel and the supply pixel as the independent variable and a preset distance decay parameter as the decay factor to calculate the decay weight within a preset distance threshold.

3. The method according to claim 1, characterized in that, The calculation of the supply-demand interaction probability between the demand pixel and each supply pixel within the distance threshold includes: Multiply the water resource supply attraction of the supplied pixel by the distance decay function to obtain the supply and demand interaction intensity value; Divide the supply-demand interaction intensity value by the sum of the supply-demand interaction intensity values ​​of all supply pixels within the distance threshold to obtain the supply-demand interaction probability of the demand pixel to the supply pixel.

4. The method according to claim 1, characterized in that, The calculation of the unit service capacity of the supplied pixel includes: The water resource supply attraction of the supplied pixels is taken as the molecule; The sum of the product of the demand intensity of the demand pixel and the probability of supply-demand interaction within the distance threshold is used as the denominator; The ratio of the numerator to the denominator is used as the unit service capacity of the supplied pixel.

5. The method according to claim 1, characterized in that, The demand intensity of the demand pixel is the sum of the dimensionless social water demand index and the dimensionless ecological water demand index.

6. The water resources security assessment method according to claim 1, characterized in that, The calculation of water resource availability for the required pixel includes: The supply contribution value is obtained by multiplying the unit service capacity of the supplied pixel by the supply-demand interaction probability. The water resource availability of the demand pixel is obtained by summing the supply contribution values ​​of all supply pixels within the distance threshold.

7. The water resources security assessment method according to claim 1, characterized in that, The construction of the water resource security index includes: The ratio of water resource availability to demand intensity of demand pixels is used as the water resource security index of demand pixels. When the demand intensity is zero, a very small positive number is added to the denominator to prevent division by zero error.

8. A water resource security assessment device based on water quantity-water quality coupling for spatial matching of supply and demand, characterized in that, include: The data preprocessing unit is used to preprocess the supply pixels and demand pixels respectively, and to make the water quantity index and water quality index of the supply pixels dimensionless, and to make the social water demand index and ecological water demand index of the demand pixels dimensionless. The supply attraction construction unit is used to construct a water resource supply attraction function for each supply cell based on the dimensionless water quantity index and water quality index. The water resource supply attraction function is used to characterize the comprehensive supply capacity of the supply cell. A distance attenuation construction unit is used to construct a distance attenuation function, which is used to quantify the impact of the spatial distance between the demand pixel and the supply pixel on the intensity of supply-demand interaction. An interaction probability calculation unit is used to calculate, for each demand pixel, the supply-demand interaction probability between the demand pixel and each supply pixel within a distance threshold, based on the distance decay function and the water resource supply attractiveness. The unit service capacity calculation unit is used to calculate the unit service capacity of each supply cell based on the supply-demand interaction probability and the demand intensity of the demand cell. A water availability calculation unit is used to calculate the water availability of each demand pixel based on the supply-demand interaction probability and the unit service capacity of the supply pixel. A safety index construction unit is used to construct a water resource safety index based on the ratio of the water resource availability to the demand intensity of the demand pixels.

9. An electronic device, characterized in that, The device includes: a processor and a memory; The memory is used to store instructions; The processor is configured to execute the instructions in the memory to perform the method according to any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, Including instructions that, when run on a computer, cause the computer to perform the method described in any one of claims 1-7 above.

Images