Drainage basin resource allocation optimization method considering drainage basin water environment risk

Abstract

The invention relates to the technical field of water conservancy projects, and discloses a drainage basin resource allocation optimization method considering drainage basin water environment risks, comprising the following steps: obtaining drainage basin pollution discharge data; constructing a drainage basin resource allocation double-layer optimization model by taking minimization of the drainage basin water environment risk coefficient as an upper layer target and minimization of the drainage basin ecological compensation cost allocation Gini coefficient as a lower layer target; processing uncertain parameters of the basin water environment capacity by adopting an adjustable robust optimization method, and constructing a basin resource allocation double-layer robust optimization model; and solving the drainage basin resource allocation double-layer robust optimization model to obtain an optimal drainage basin ecological compensation cost and discharge capacity allocation result. Aiming at the problems of the traditional drainage basin pollution discharge distribution technology, the drainage basin ecological compensation method which comprehensively considers the sub-region water resource utilization efficiency, the water environment quality index and the water ecological service index is adopted to make up the deficiencies and defects of the original method, and the accuracy and feasibility of the distribution scheme are enhanced.

Description

Technical Field

[0001] The present invention relates to the technical field of water conservancy projects, and in particular to a method for optimizing watershed resource allocation taking into account water environment risks in the watershed. Background Art

[0002] Traditional methods for allocating pollution rights often rely on static indicators, ignoring the hydrological conditions of a river basin, the spatiotemporal heterogeneity of environmental capacity, and the migration and diffusion patterns of pollutants. This leads to a disconnect between total pollution discharge control and the region's actual carrying capacity. Furthermore, traditional methods employ relatively few decision-makers, have relatively consistent objectives, and offer relatively low decision-making complexity, which can easily lead to conflicts of interest among participating parties. Therefore, effectively implementing pollution discharge allocation, designing ecological compensation strategies, and controlling river basin water environmental risk indicators have become key issues in river water environmental governance. Summary of the Invention

[0003] In view of the above-mentioned deficiencies in the prior art, the present invention provides a method for optimizing watershed resource allocation taking into account water environment risks in the watershed.

[0004] In order to achieve the above-mentioned object of the invention, the technical solution adopted by the present invention is: A method for optimizing watershed resource allocation considering watershed water environment risks includes the following steps: Obtain basin pollution data; A two-level optimization model for watershed resource allocation is constructed with minimizing the water environment risk coefficient of the watershed as the upper level goal and minimizing the Gini coefficient of watershed ecological compensation cost allocation as the lower level goal. An adjustable robust optimization method is used to deal with the uncertain parameters of the water environment capacity of the basin, and a two-layer robust optimization model for basin resource allocation is constructed; The two-layer robust optimization model of watershed resource allocation is solved to obtain the optimal watershed ecological compensation cost and pollution discharge allocation results.

[0005] Furthermore, with minimizing the basin water environment risk coefficient as the upper level goal and minimizing the basin ecological compensation cost allocation Gini coefficient as the lower level goal, a two-level optimization model for basin resource allocation is constructed, including: The upper-level model is constructed with the minimization of the basin water environment risk coefficient as the upper-level goal and the upper limit constraint of water resource use, the safety margin constraint of water pollutants in the basin, the allocation constraint of the basin ecological compensation cost, and the water environment quality threshold constraint as the constraint conditions. Taking minimizing the Gini coefficient of watershed ecological compensation cost allocation as the lower-level objective, and taking the upper limit constraint of the Gini coefficient of watershed ecological compensation cost allocation and the total amount constraint of watershed ecological compensation cost as the constraint conditions, the lower-level model is constructed. A two-layer optimization model for watershed resource allocation is constructed based on the upper-layer model and the lower-layer model.

[0006] Furthermore, minimizing the risk factor of the water environment in the basin is the upper-level goal, specifically: ; Among them, min is the minimum value function; is the basin water environment risk coefficient; For the region water pollutants In the cycle emissions; Water pollutants for the basin In the cycle The upper limit value that can be accommodated; r is the total number of cycles; m is the amount of water pollutants; n is the number of regions.

[0007] Furthermore, the constraints are the upper limit of water resource use, the safety margin of water pollutants in the watershed, the allocation of ecological compensation costs in the watershed, and the water environment quality threshold. Specifically,

[0008]

[0009]

[0010]

[0011] in, Water pollutants In the cycle Concentration index; Water pollutants In the cycle Upper limit of grey water usage; Water pollutants In the cycle safety margin; Water pollutants In the cycle The water environment capacity of the basin; Water pollutants In the cycle The lower limit of the safety margin; Compensation costs for watershed ecology in the cycle the total amount; Water pollutants for the basin In the cycle Demand; For the region In the cycle The amount of ecological compensation cost allocated.

[0012] Furthermore, minimizing the Gini coefficient of watershed ecological compensation cost allocation is the lower-level goal, specifically:

[0013] in, assigning a Gini coefficient to the cost of ecological compensation for the watershed; For the region In the cycle Water features l index; For the region o In the cycle the amount of ecological compensation cost allocated; For the region o In the cycle Water features l index; For the region p In the cycle the amount of ecological compensation cost allocated; For the region p In the cycle Water features l index.

[0014] Furthermore, the upper limit constraint of the Gini coefficient optimization of the watershed ecological compensation cost distribution and the total amount constraint of the watershed ecological compensation cost are used as constraints, specifically:

[0015] in, is the Gini coefficient of watershed ecological compensation cost distribution before optimization.

[0016] Furthermore, the adjustable robust optimization method is used to deal with the uncertain parameters of the water environment capacity of the basin: The water environment capacity of the basin is converted into uncertainty parameters to obtain the water environment capacity interval of the basin; The relationship between pollutant emissions and water environment capacity of the basin is established using affine function.

[0017] Furthermore, the relationship between pollutant emissions and water environment capacity of the basin is established using the affine function:

[0018] in, It is a non-adjustable parameter; is the initial value of the water environment capacity of the basin; is the uncertainty parameter of the water environment capacity of the basin; λ For cycle Previous cycle.

[0019] Furthermore, a two-layer robust optimization model for watershed resource allocation is constructed as follows:

[0020] in, is a new auxiliary variable; Water pollutants In the cycle Uncertain parameters of water environment capacity of inner watershed; Water pollutants In the cycle Standard value of water environment capacity of inner watershed; is the fluctuation parameter of the water environment capacity of the basin; It is the water environment quality index.

[0021] Furthermore, the two-layer robust optimization model for watershed resource allocation is solved to obtain the optimal watershed ecological compensation cost and pollution discharge allocation results, including: Calculate the maximum and minimum values of the upper and lower layer objective functions and construct the satisfaction function of the upper and lower layer objective functions; Set the minimum value of the upper layer satisfaction function and the maximum value of the lower layer satisfaction function, and take maximizing the minimum value of the lower layer satisfaction function as the new goal, build a new single-layer single-objective model, and obtain the initial solution of the model; Set the minimum satisfaction of the upper objective function and the maximum satisfaction of the lower objective function, build a new single-level objective model, and solve the initial value while maximizing the minimum satisfaction of the lower objective function; By adjusting the fluctuation parameters of the water environment capacity of the basin, a series of optimal solutions of the optimization model under different uncertainty levels are obtained, and the optimal basin ecological compensation cost and pollution discharge allocation results are obtained.

[0022] The present invention has the following beneficial effects: The present invention addresses the problems existing in traditional watershed sewage allocation technology and adopts a watershed ecological compensation method that comprehensively considers the water resource utilization efficiency, water environment quality index and water ecological service index of the sub-region to make up for the shortcomings and defects of the original method. In addition, the uncertainty factors of key parameters that change over time in the decision-making process are processed using an adjustable robust optimization method, thereby enhancing the accuracy and feasibility of the allocation plan. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] Figure 1 This is a flow chart of a watershed resource allocation optimization method that takes into account water environment risks in the watershed; Figure 2 This is a schematic diagram of the framework of a watershed resource allocation optimization method that takes into account water environment risks in the watershed. DETAILED DESCRIPTION

[0024] The specific embodiments of the present invention are described below to facilitate understanding of the present invention by those skilled in the art. However, it should be clear that the present invention is not limited to the scope of the specific embodiments. For those skilled in the art, as long as various changes are within the spirit and scope of the present invention as defined and determined by the appended claims, these changes are obvious, and all inventions and creations utilizing the concepts of the present invention are protected.

[0025] like Figure 1 and Figure 2 As shown, an embodiment of the present invention provides a method for optimizing watershed resource allocation taking into account water environment risks in the watershed, comprising the following steps S1 to S4: S1. Obtain basin pollution data; In an optional embodiment of the present invention, step S1 obtains basin pollution discharge data as raw data to perform basin resource allocation optimization processing.

[0026] S2. Construct a two-level optimization model for watershed resource allocation, with minimizing the water environment risk coefficient of the watershed as the upper level goal and minimizing the Gini coefficient of watershed ecological compensation cost allocation as the lower level goal; In an optional embodiment of the present invention, step S2 takes minimizing the water environment risk coefficient of the watershed as the upper-level objective and minimizing the Gini coefficient of the watershed ecological compensation cost allocation as the lower-level objective, and constructs a two-level optimization model for watershed resource allocation, including: The upper-level model is constructed with the minimization of the basin water environment risk coefficient as the upper-level goal and the upper limit constraint of water resource use, the safety margin constraint of water pollutants in the basin, the allocation constraint of the basin ecological compensation cost, and the water environment quality threshold constraint as the constraint conditions. Taking minimizing the Gini coefficient of watershed ecological compensation cost allocation as the lower-level objective, and taking the upper limit constraint of the Gini coefficient of watershed ecological compensation cost allocation and the total amount constraint of watershed ecological compensation cost as the constraint conditions, the lower-level model is constructed. A two-layer optimization model for watershed resource allocation is constructed based on the upper-layer model and the lower-layer model.

[0027] In the context of ecological and environmental protection, step S2 focuses on the water environment risk issue in the process of water pollution discharge, with the goal of minimizing the water environment risk coefficient of the watershed. The water environment risk coefficient of the watershed is the ratio of the actual discharge of water pollutants in period t to the upper limit of water pollutants in the watershed. Therefore, minimizing the water environment risk coefficient of the watershed is the upper-level goal, specifically: ; Among them, min is the minimum value function; is the basin water environment risk coefficient; For the region water pollutants In the cycle emissions; Water pollutants for the basin In the cycle The upper limit value that can be accommodated; r is the total number of cycles; m is the amount of water pollutants; n is the number of regions.

[0028] Water resource usage upper limit constraint. Grey water refers to the amount of fresh water consumed to ensure that the water pollutant load in the basin meets the environmental water quality standards. This embodiment sets an upper limit on the use of grey water. The grey water volume is calculated by the ratio of the discharge of water pollutants to the water quality concentration standard, and each pollutant In the cycle All requirements must be met:

[0029] The safety margin of water pollutants in the river basin refers to the amount of water pollutants that are reserved in advance to prevent the adverse effects of internal and external factors on the water environment quality goals, so as to ensure the safety of the water environment quality in the river basin. The safety margin of water pollutants is equal to the safety margin of water pollutants in the previous period plus the water environment capacity in the current period, minus the discharge of water pollutants in the current period. In the cycle The safety margin must not be lower than the lower limit requirement:

[0030] Constraints on the allocation of watershed ecological compensation funds. In the process of pollutant discharge, it is necessary to consider the cycle The accounting problem after the completion is that the discharge of water pollutants in the basin will affect the allocation of ecological compensation costs. Therefore, the proportion of water pollutant emissions to the upper limit of water pollutant emissions in the basin should be less than or equal to the ratio of the ecological compensation costs received to the total amount:

[0031] Water environment quality bottom line constraints. The discharge of water pollutants in the entire river basin should not exceed the upper limit of the water pollutant capacity of the river basin, and at the same time should meet the basic development needs of the river basin:

[0032] in, Water pollutants In the cycle Concentration index; Water pollutants In the cycle Upper limit of grey water usage; Water pollutants In the cycle safety margin; Water pollutants In the cycle The water environment capacity of the basin; Water pollutants In the cycle The lower limit of the safety margin; Compensation costs for watershed ecology in the cycle the total amount; Water pollutants for the basin In the cycle demand; For the region In the cycle The amount of ecological compensation cost allocated.

[0033] In step S2, when allocating the cost of ecological compensation for water environment in the watershed, it is necessary to comprehensively consider the water resource utilization efficiency, water environment quality index and water ecological service index. This embodiment expands the Gini Coefficient and constructs the Gini Coefficient for the allocation of watershed ecological compensation funds to quantify and characterize the fairness of the allocation of watershed ecological compensation costs. Water resource utilization efficiency refers to the period The water consumption per 10,000 yuan of GDP in the region, the water environment quality index is the difference between the sub-region emission baseline value and the water pollutant emission (upper-level decision variable), that is, The water ecological service index refers to the percentage of the area of ecological land and protected areas in the sub-region to the watershed area of each city. The allocation function of fairness is:

[0034] in, assigning a Gini coefficient to the cost of ecological compensation for the watershed; For the region In the cycle Water features l index; For the region o In the cycle the amount of ecological compensation cost allocated; For the region o In the cycle Water features l index; For the region p In the cycle the amount of ecological compensation cost allocated; For the region p In the cycle Water features l index.

[0035] The fairness of the distribution of watershed ecological compensation funds should be improved after optimization; the total amount of watershed ecological compensation funds in period t is fixed and can be expressed as:

[0036] in, is the Gini coefficient of watershed ecological compensation cost distribution before optimization.

[0037] Finally, the two-layer robust optimization configuration model is expressed as a global model: .

[0038] S3. Adopt an adjustable robust optimization method to deal with the uncertain parameters of the water environment capacity of the basin and construct a two-layer robust optimization model for basin resource allocation; In an optional embodiment of the present invention, step S3 considers the problem of water environment capacity uncertainty in the decision-making process and uses an adjustable robust optimization method to solve it, including: S3-1. The water environment capacity parameters of the basin Replace with To characterize uncertainty; S3-2, assuming that in the cycle The discharge of water pollutants is calculated starting from the beginning of the period, and the exact value of the water environment capacity in the subsequent period is unknown. The range of the water environment capacity is known:

[0039] in, Uncertain parameters representing the water environment capacity of the basin; is the fluctuation parameter of water environment capacity; It is a pollutant In the cycle Nominal value of water environment capacity of the inner basin.

[0040] S3-3, in the cycle Decisions made at the beginning It is based on the historical data of water environment capacity in the early stage, i.e. and There is a correlation, where Represents the previous cycle, using affine function to build and The relationship between:

[0041] in, It is a non-adjustable parameter; is the initial value of the water environment capacity of the basin; is the uncertainty parameter of the water environment capacity of the basin; λ For cycle Previous cycle.

[0042] S3-4, through the above steps, the global model All of them are transformed equivalently, and finally a two-layer robust optimization model is constructed that considers the coordinated configuration of pollution discharge and ecological compensation costs:

[0043] in, It is a newly introduced auxiliary variable to facilitate the solution of the model; Water pollutants In the cycle Uncertain parameters of water environment capacity of inner watershed; Water pollutants In the cycle Standard value of water environment capacity of inner watershed; is the fluctuation parameter of the water environment capacity of the basin; It is the water environment quality index.

[0044] S4. Solve the two-layer robust optimization model of watershed resource allocation to obtain the optimal watershed ecological compensation cost and pollution discharge allocation results.

[0045] In an optional embodiment of the present invention, step S4 solves the two-layer robust optimization model for watershed resource allocation to obtain the optimal watershed ecological compensation cost and pollution discharge allocation results, including: S4-1. Construct the satisfaction function of the upper and lower layer objective functions by obtaining the maximum and minimum values of the upper and lower layer objective functions; S4-2. Set the minimum value of the upper-layer satisfaction function and the maximum value of the lower-layer satisfaction function, take maximizing the minimum value of the lower-layer satisfaction function as the new goal, construct a new single-layer single-objective model, and obtain the initial solution of the model; S4-3. Set the minimum satisfaction level for upper-level goals and maximum satisfaction of lower-level goals ; S4-4. Construct a new single-level objective model and solve the initial value while maximizing the minimum satisfaction of the lower-level decision makers; S4-5, Adjustment By taking the value of , we can obtain a series of optimal solutions of the optimization model under different uncertainty levels, and obtain the optimal watershed ecological compensation cost and pollution discharge allocation results.

[0046] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowcharts and / or block diagrams, as well as combinations of processes and / or blocks in the flowcharts and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowcharts and / or block diagrams. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0047] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0048] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0049] Specific embodiments are used in the present invention to illustrate the principles and implementation methods of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core ideas. At the same time, for those skilled in the art, according to the ideas of the present invention, there may be changes in the specific implementation methods and application scopes. In summary, the contents of this specification should not be understood as limiting the present invention.

[0050] Those skilled in the art will appreciate that the embodiments described herein are intended to help readers understand the principles of the present invention, and it should be understood that the scope of protection of the present invention is not limited to such specific descriptions and embodiments. Those skilled in the art can make various other specific variations and combinations based on the technical teachings disclosed in the present invention without departing from the essence of the present invention, and such variations and combinations are still within the scope of protection of the present invention.

Claims

1. A basin resource allocation optimization method considering basin water environment risks, characterized by: The following steps are involved: Obtain basin pollution data; A two-level optimization model for watershed resource allocation is constructed with minimizing the water environment risk coefficient of the watershed as the upper level goal and minimizing the Gini coefficient of watershed ecological compensation cost allocation as the lower level goal. An adjustable robust optimization method is used to deal with the uncertain parameters of the water environment capacity of the basin, and a two-layer robust optimization model for basin resource allocation is constructed; The two-layer robust optimization model of watershed resource allocation is solved to obtain the optimal watershed ecological compensation cost and pollution discharge allocation results.

2. A basin resource allocation optimization method considering basin water environment risks according to claim 1, characterized in that: With minimizing the water environment risk coefficient of the basin as the upper level goal and minimizing the Gini coefficient of the ecological compensation cost allocation of the basin as the lower level goal, a two-level optimization model for basin resource allocation is constructed, including: The upper-level model is constructed with the minimization of the basin water environment risk coefficient as the upper-level goal and the upper limit constraint of water resource use, the safety margin constraint of water pollutants in the basin, the allocation constraint of the basin ecological compensation cost, and the water environment quality threshold constraint as the constraint conditions. Taking minimizing the Gini coefficient of watershed ecological compensation cost allocation as the lower-level objective, and taking the upper limit constraint of the Gini coefficient of watershed ecological compensation cost allocation and the total amount constraint of watershed ecological compensation cost as the constraint conditions, the lower-level model is constructed. A two-layer optimization model for watershed resource allocation is constructed based on the upper-layer model and the lower-layer model.

3. A basin resource allocation optimization method considering basin water environment risks according to claim 2, characterized in that: The upper level goal is to minimize the risk factor of the water environment in the basin, specifically: ； Among them, min is the minimum value function; is the basin water environment risk coefficient; For the region water pollutants In the cycle emissions; Water pollutants for the basin In the cycle The upper limit value that can be accommodated; r is the total number of cycles; m is the amount of water pollutants; n is the number of regions.

4. A method for optimizing watershed resource allocation considering watershed water environment risks according to claim 3, characterized in that: The constraints are the upper limit of water resource use, the safety margin of water pollutants in the watershed, the allocation of ecological compensation costs in the watershed, and the water environment quality threshold. Specifically, in, Water pollutants In the cycle Concentration index; Water pollutants In the cycle Upper limit of grey water usage; Water pollutants In the cycle safety margin; Water pollutants In the cycle The water environment capacity of the basin; Water pollutants In the cycle The lower limit of the safety margin; Compensation costs for watershed ecology in the cycle the total amount; Water pollutants for the basin In the cycle Demand; For the region In the cycle The amount of ecological compensation cost allocated.

5. A basin resource allocation optimization method considering basin water environment risks according to claim 4, characterized in that: The lower-level goal is to minimize the Gini coefficient of watershed ecological compensation cost allocation, specifically: in, Allocate Gini coefficients for watershed ecological compensation costs; For the region In the cycle Water features l index; For the region o In the cycle the amount of ecological compensation cost allocated; For the region o In the cycle Water features l index; For the region p In the cycle the amount of ecological compensation cost allocated; For the region p In the cycle Water features l index.

6. A method for optimizing watershed resource allocation considering watershed water environment risks according to claim 5, characterized in that: The upper limit constraint of the Gini coefficient optimization of the watershed ecological compensation cost allocation and the total amount constraint of the watershed ecological compensation cost are used as the constraints, specifically: in, is the Gini coefficient of watershed ecological compensation cost distribution before optimization.

7. A method for optimizing watershed resource allocation considering watershed water environment risks according to claim 6, characterized in that: The specific parameters of the uncertain water environment capacity of the basin are handled by the adjustable robust optimization method: The water environment capacity of the basin is converted into uncertainty parameters to obtain the water environment capacity interval of the basin; The relationship between pollutant emissions and water environment capacity of the basin is established using affine function.

8. A method for optimizing watershed resource allocation considering watershed water environment risks according to claim 7, characterized in that: The relationship between pollutant emissions and water environment capacity of the basin is established using the affine function as follows: in, It is a non-adjustable parameter; is the initial value of the water environment capacity of the basin; is the uncertainty parameter of the water environment capacity of the basin; λ For cycle Previous cycle.

9. A method for optimizing watershed resource allocation considering watershed water environment risks according to claim 8, characterized in that: The specific steps of constructing a two-layer robust optimization model for watershed resource allocation are as follows: in, is a new auxiliary variable; Water pollutants In the cycle Uncertain parameters of water environment capacity of inner watershed; Water pollutants In the cycle Standard value of water environment capacity of inner watershed; is the fluctuation parameter of the water environment capacity of the basin; It is the water environment quality index.

10. A method for optimizing watershed resource allocation considering watershed water environment risks according to claim 9, characterized in that: The two-layer robust optimization model for watershed resource allocation is solved to obtain the optimal watershed ecological compensation cost and pollution discharge allocation results, including: Calculate the maximum and minimum values of the upper and lower layer objective functions and construct the satisfaction function of the upper and lower layer objective functions; Set the minimum value of the upper layer satisfaction function and the maximum value of the lower layer satisfaction function, and take maximizing the minimum value of the lower layer satisfaction function as the new goal, build a new single-layer single-objective model, and obtain the initial solution of the model; Set the minimum satisfaction of the upper objective function and the maximum satisfaction of the lower objective function, build a new single-level objective model, and solve the initial value while maximizing the minimum satisfaction of the lower objective function; By adjusting the fluctuation parameters of the water environment capacity of the basin, a series of optimal solutions of the optimization model under different uncertainty levels are obtained, and the optimal basin ecological compensation cost and pollution discharge allocation results are obtained.

Images