Land planning optimization method based on water environment quality constraint, medium and equipment

Abstract

The invention relates to the technical field of urban planning, and provides a water environment quality constraint-based land planning optimization method, medium and equipment, and the method comprises the steps: dividing a research region into n drainage subregions; estimating a runoff non-point source pollution load total amount and a point source pollution load total amount; estimating the total pollution load; performing feasibility evaluation on the land planning scheme, and judging whether the land planning scheme meets a rigid constraint condition based on a water environment target or not; and optimizing a land utilization planning scheme, firstly determining a target function of a receiving water body of the research area by taking the pollutant concentration of the receiving water body as a target, judging whether the target function meets a water quality control target of the receiving water body or not, and if not, constructing an optimal distribution model by adopting a linear programming function to optimize the land utilization planning scheme. According to the method, the corresponding high-risk drainage partition can be precisely traced, the water environment problem caused by urban development is avoided from the planning level, and urban land utilization planning is more reasonable.

Description

Technical Field

[0001] The present invention relates to the technical field of urban planning, and in particular to a land planning optimization method, medium and equipment based on water environment quality constraints. Background Art

[0002] Urban planning is the premise of a city's development and construction. Urban planning is a comprehensive discipline that focuses on promoting urban development. It creates places suitable for people to live, study, work, etc. through reasonable regional construction and reasonable allocation. Urban planning is essentially a standardized division of urban functional areas. According to the spatial positioning of urban development, the needs of urban development, and the density of population residence, the corresponding functional areas are standardized to further promote the expansion of urban scale and urban development. Urban planning is not only normative, but also scientific, forward-looking, strategic and holistic. These characteristics and its sustainable development are the guarantee of good urban development.

[0003] Due to the expansion of urban hardened underlying surface area, the historical debt of underground rainwater and sewage pipeline construction and serious lack of investment, surface rainwater runoff and combined system overflow and other non-point source pollution problems are prominent, the intensity of damage to the watershed water environment is increasing, various water problems are emerging, and the conflict between people, water and cities is becoming increasingly acute. The excessive development and irrational use of land caused by the traditional urban planning model exceeds the carrying capacity of water resources and water environment, which is one of the main reasons for the crisis of water resources, water environment and water ecology. Summary of the invention

[0004] The purpose of the present invention is to provide a land planning optimization method based on water environment quality constraints, which can avoid water environment pollution from the planning level, improve the scientific nature of urban land use, and reduce the investment in subsequent water environment improvement. The specific technical scheme is as follows:

[0005] The present invention provides a land planning optimization method based on water environment quality constraints, comprising the following steps:

[0006] S1: Divide the study area into n drainage zones and number each drainage zone, denoted as i, where i is a natural number from 1 to n; and obtain the total planned construction land area and land use planning scheme of each drainage zone;

[0007] S2: Estimation of the total amount of runoff non-point source pollution load, specifically: based on the total planned construction land area of ​​each drainage zone and the average concentration of runoff pollutants, the total amount of runoff non-point source pollution load in the study area is obtained;

[0008] The total amount of point source pollution load is estimated as follows: the total amount of point source pollution load in the study area is calculated based on the average flow of point source pollution sewage and the point source pollution emission concentration in each drainage zone;

[0009] S3: Estimation of the total amount of pollution load, specifically: Calculate the total amount of pollution load w0 in the study area based on the drainage system, point source discharge outlet location, total amount of runoff surface source pollution load and total amount of point source pollution load in the study area;

[0010] S4: Feasibility assessment of land use planning scheme, specifically: use the zero-dimensional model to calculate the allowable discharge level M of the receiving water body in the study area. If w0>M, it means that the land use planning scheme in the study area does not meet the rigid constraint conditions based on the water environment quality as the goal, and enter S5; if w0≤M, it means that the land use planning scheme is feasible;

[0011] S5: Optimization of land use planning scheme, specifically: Taking the pollutant concentration of receiving water as the target, first determine the objective function C of the receiving water in the study area i , if C i ≤S S This means that the receiving water body of the i-th drainage zone and its upstream land use planning scheme are feasible, and the i-th drainage zone still has sufficient pollution absorption capacity; if C i >S S The linear programming function is used to construct an optimal allocation model, and the receiving water body of the i-th drainage zone and the high-risk drainage zone land use planning scheme in its upstream are adjusted through the optimal allocation model, and S2 is returned, where: S S It is the water quality control target of receiving water bodies.

[0012] Optionally, the method for obtaining the total amount of runoff surface source pollution load in S2 includes:

[0013] According to the relevant classification standards of land planning and the emission characteristics of major pollutants in surface runoff, the planned construction land in each drainage zone is reclassified into four types: I, II, III, and IV;

[0014] Determine the comprehensive runoff coefficient φ of the drainage zone according to the planned town type and relevant standards, and obtain the average annual rainfall of the drainage zone;

[0015] The average concentration of runoff pollutants is determined by calculating the average concentration of runoff pollutants EMC in the drainage zone based on area weighting. The calculation formula is as follows:

[0016]

[0017] Among them: F I The land area of ​​type I in the drainage zone, F II The land area of ​​Category II in the drainage zone, F III The land area of ​​Category III in the drainage zone, F IV Class IV land area for drainage zones; EMCI The average concentration of pollutants in runoff from Class I land in the drainage zone, EMC II The average concentration of pollutants in runoff from Class II land in the drainage zone, EMC III The average concentration of pollutants in runoff from Class III land use in the drainage zone, EMC IV is the average concentration of pollutants in runoff from Class IV land in the drainage zone; F is the total area of ​​planned construction land in the drainage zone;

[0018] The volumetric method is used to calculate the annual runoff volume of each drainage zone. The calculation formula is as follows:

[0019] V Yi =10H·φ·F;

[0020] Where: V Yi is the annual runoff volume of the ith drainage zone; H is the annual rainfall in a wet year or a dry year;

[0021] The total amount of runoff non-point source pollution load in the study area is calculated based on the annual runoff volume of each drainage zone and the average runoff pollutant concentration EMC, specifically:

[0022] Calculate the annual runoff pollution load of each drainage zone using the following formula:

[0023] w Yi =10 -3 ·V Yi EMC;

[0024] Where: w Yi is the annual runoff pollution load of the ith drainage zone;

[0025] The total runoff pollution load w in the study area is calculated based on the annual runoff pollution load of each drainage zone. Y0 , the calculation formula is as follows:

[0026]

[0027] Optionally, the S2 further includes:

[0028] The total annual runoff volume of the study area is calculated based on the annual runoff volume of each drainage zone, specifically:

[0029] The annual average runoff flow of each drainage zone is calculated based on the annual runoff volume of each drainage zone. The calculation formula is as follows:

[0030]

[0031] Where: q Yi is the average annual runoff of the ith drainage zone, T is the dry year or wet year;

[0032] The total annual runoff flow of the study area is calculated based on the annual average runoff flow of each drainage zone. The calculation formula is as follows:

[0033]

[0034] Where: q Y0 is the total annual runoff volume in the study area.

[0035] Optionally, the total point source pollution load w in the study area in S2 w0 The methods for obtaining include:

[0036] The average flow rate of point source pollution sewage in the i-th drainage area is q wi The calculation formula is as follows:

[0037]

[0038] Where: γ i is the water consumption index of different land use types in the i-th drainage zone; a i is the area of ​​different land use types in the i-th drainage zone; β is the sewage discharge coefficient;

[0039] The annual sewage flow rate q in the study area w0 The calculation formula is as follows:

[0040]

[0041] The annual point source pollution load w of the i-th drainage area wi The calculation formula is as follows:

[0042] w wi =10 -3 ·q wi ·c wi T;

[0043] Where: c wi is the point source pollution emission concentration in the i-th drainage zone;

[0044] The total amount of point source pollution load w in the study area w0 The calculation formula is as follows:

[0045]

[0046] Optionally, in S3, the drainage system includes a separate system and a combined system;

[0047] If the drainage system in the study area is a diversion system, and there is no single point source flowing into the receiving water body in the study area, and the influent pollution load only considers the runoff surface source pollution load, then the calculation formula for the total flow of rainwater and sewage discharged into the receiving water body, the pollutant concentration, and the total pollution load in the study area is as follows:

[0048] q0=q Y0 ;

[0049] c0=EMC;

[0050] w0=w Y0 ;

[0051] Where: q0 is the total flow of rainwater and sewage discharged into the receiving water body in the study area, q Y0 is the total annual runoff discharged into the receiving water body in the study area; c0 is the pollutant concentration discharged into the receiving water body in the study area;

[0052] If the drainage system of the study area is a diversion system, and there is a separate point source in the study area that discharges into the receiving water body, the influent pollution load considers the runoff surface source pollution load and the point source pollution load. The calculation formula for the total flow of rainwater and sewage discharged into the receiving water body, the pollutant concentration, and the total pollution load of the study area are as follows:

[0053] q0=q Y0 +q w0 ;

[0054]

[0055] w0=w Y0 +w w0 ;

[0056] If the drainage system of the study area is a combined system, and there is no separate point source in the study area that discharges into the receiving water body, the runoff non-point source pollution load and point source pollution are discharged into the receiving water body in the form of combined system overflow pollution. The calculation formula for the total flow of rainwater and sewage discharged into the receiving water body, the pollutant concentration, and the total pollution load of the study area are as follows:

[0057] q0=q Y0 -mq w0 ;

[0058] c0=EMC;

[0059] w0=(q Y0 -mq w0 )EMC;

[0060] Where: m is the interception multiple;

[0061] If the drainage system in the study area is a combined system, and there is a separate point source in the study area that discharges into the receiving water body, the influent pollution load takes into account the overflow pollution of the combined system and the pollution load of the point source. The calculation formula for the total flow of rainwater and sewage discharged into the receiving water body, the pollutant concentration, and the total pollution load of the study area are as follows:

[0062] q0=q Y0 +qw0 ;

[0063]

[0064] w0=(q Y0 -mq w0 )·EMC+(m+1)·q w0 ·C w .

[0065] Optionally, in S4, when the receiving water body in the study area is a river, the calculation formula for the allowable discharge level of the receiving water body is as follows:

[0066] M=31.54(Q0+q0)*S s -31.54Q0S0-L;

[0067] Where: is the pollutant concentration of the upstream water of the receiving water body; Q0 ​​is the upstream water flow of the receiving water body; q0 is the total flow of rainwater and sewage discharged into the receiving water body in the study area; L is the safety boundary value; K is the first-order attenuation rate constant of the receiving water body; V is the average volume of water in the receiving water body;

[0068] When the receiving water body in the study area is a lake, the calculation formula for the allowable discharge level of the receiving water body is as follows:

[0069] M=31.54(Q0+q0)*S s +KS s V-3.65×10 -4 Q0S0-L.

[0070] Optionally, in S5, a zero-dimensional model is used to generalize the impact of pollution load inflow on the receiving water body. When the receiving water body of the drainage partition is a river, the calculation formula of the objective function is as follows:

[0071]

[0072] Where: N0 is the water quality concentration of the upstream water of the receiving water body;

[0073] When the receiving water body of the drainage zone is a lake, the calculation formula of the objective function is as follows:

[0074]

[0075] Optionally, when the linear programming function in S5 cannot obtain a feasible solution for the planned construction land development plan, the concept of sponge city construction is adopted in the drainage zoning development process to control the non-point source pollution in the study area;

[0076] The expression of the optimization allocation model in S5 is as follows:

[0077] When the receiving water body of the drainage zone is a river:

[0078]

[0079]

[0080] Where: C E is the pollutant concentration of the upstream water of the river; Q E is the water flow from the upstream of the river; EMC i is the average concentration of runoff pollutants in the ith drainage area; F Ii is the area of ​​Class I land in the i-th drainage zone, F IIi is the area of ​​Class II land in the ith drainage zone, F IIIi is the land area of ​​Class III in the i-th drainage zone, F IVi is the land area of ​​Class IV in the ith drainage zone; a i is the area of ​​various types of land in the i-th drainage zone; θ max is the upper limit of the total area of ​​each type of land use in the i-th drainage zone, θ min The lower limit of the area of ​​each type of land use in the i-th drainage zone;

[0081] When the receiving water body of the drainage zone is a lake:

[0082]

[0083] Among them: F max is the upper limit of the area of ​​each type of land in the i-th drainage zone, F min It is the lower limit of the area of ​​each type of land use in the i-th drainage zone.

[0084] The present invention also provides a computer-readable storage medium having computer program instructions stored thereon, and when the computer program instructions are executed by a processor, the land planning optimization method based on water environment quality constraints as described above is implemented.

[0085] The present invention also provides a computer electronic device, comprising: at least one processor, at least one memory and computer program instructions stored in the memory, when the computer program instructions are executed by the processor as described in the land planning optimization method based on water environment quality constraints.

[0086] The present invention constructs a functional relationship between planned construction land and water environment quality, which can not only quickly determine whether the land planning scheme meets the rigid constraints based on water environment goals, but also accurately trace the corresponding high-risk drainage zones, which is of guiding significance for evaluating the feasibility of planned land use schemes; moreover, based on the linear programming function, a feasible solution of the planned construction land development scheme is obtained, which provides data support for the optimization of land planning schemes, avoids water environment problems caused by urban development from the planning level, reduces the investment in subsequent water environment remediation projects, makes urban land use planning more reasonable, and further implements the urban development strategy of determining cities, land, people and production based on water.

[0087] In addition to the above-described purposes, features and advantages, the present invention has other purposes, features and advantages. The present invention will be further described in detail with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0088] The drawings constituting a part of this application are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings:

[0089] Figure 1 is a flow chart of a land planning optimization method based on water environment quality constraints in an embodiment of the present invention;

[0090] Figure 2 It is a schematic diagram of the relationship between drainage partitions and receiving water bodies in an embodiment of the present invention. DETAILED DESCRIPTION

[0091] The embodiments of the present invention are described in detail below with reference to the accompanying drawings, but the present invention can be implemented in many different ways as defined and covered by the claims.

[0092] In one embodiment, if Figure 1 and Figure 2 As shown, a land planning optimization method based on water environment quality constraints includes the following steps:

[0093] S1: Divide the study area into n drainage zones and number each drainage zone, denoted as ; and obtain the total planned construction land area and land use planning scheme of each drainage zone.

[0094] When dividing drainage zones, the results of existing drainage special planning and drainage zones can be directly used; or drainage zones can be automatically generated using GIS software based on high-precision DEM elevation data.

[0095] S2: Estimation of the total amount of runoff non-point source pollution load, specifically: based on the total planned construction land area of ​​each drainage zone and the average concentration of runoff pollutants, the total amount of runoff non-point source pollution load in the study area is obtained;

[0096] The methods for obtaining the total amount of runoff non-point source pollution load in S2 include:

[0097] (I) According to the relevant classification standards of land planning and the emission characteristics of major pollutants in surface runoff, the planned construction land in each drainage zone is reclassified into four types: I, II, III, and IV. Taking into account the similarity of the pollution degree of runoff rainwater on the underlying surface of the same land use nature in cities and towns, and referring to the relevant research on runoff water quality, according to the code category of GB50137 of "Urban Land Classification and Planning and Construction Land Standards", the planned construction land types with relatively similar emission characteristics of major pollutants in surface runoff are merged into the same category, and reclassified into four types: I, II, III, and IV, as shown in Table 1.

[0098] Table 1 Classification of planned construction land based on surface runoff pollutant emission characteristics

[0099]

[0100] (II) Determine the comprehensive runoff coefficient φ of the drainage zone according to the planned town type and relevant standards, and obtain the average annual rainfall of the drainage zone; according to the planned town type, refer to the "Outdoor Drainage Design Standard" GB50014 to determine the value of the comprehensive runoff coefficient of the drainage zone, among which the urban building-dense area is 0.6-0.7, the urban building-dense area is 0.45-0.6, and the urban building-sparse area is 0.2-0.45. Comprehensively consider the characteristics of non-point source pollution in wet years and dry years to determine its impact on the water environment, and the designed annual average rainfall H (mm) can be determined according to the annual rainfall in wet years or dry years.

[0101] (III) Determination of average concentration of runoff pollutants: Specifically, the average concentration of runoff pollutants EMC of drainage zones is calculated based on area weighting. The calculation formula is as follows:

[0102]

[0103] Among them: F I The land area of ​​type I in the drainage zone, F II The land area of ​​Category II in the drainage zone, F III The land area of ​​Category III in the drainage zone, F IV Class IV land area for drainage zones; EMC I The average concentration of pollutants in runoff from Class I land in the drainage zone, EMC II The average concentration of pollutants in runoff from Class II land in the drainage zone, EMCIII The average concentration of pollutants in runoff from Class III land use in the drainage zone, EMC IV is the average concentration of runoff pollutants from Class IV land in the drainage zone, in mg / L; F is the total planned construction land area in the drainage zone, in hm2 2 ;

[0104] (IV) The annual runoff volume of each drainage zone is calculated using the volumetric method. The calculation formula is as follows:

[0105] V Yi =10H·φ·F;

[0106] Where: V Yi is the annual runoff volume of the ith drainage zone, in m 3 ; H is the annual rainfall in a wet year or a dry year, in mm;

[0107] (V) Calculate the total runoff non-point source pollution load in the study area based on the annual runoff volume of each drainage zone and the average runoff pollutant concentration EMC, specifically:

[0108] 1. Calculate the annual runoff pollution load of each drainage zone. The calculation formula is as follows:

[0109] w Yi =10 -3 ·V Yi EMC;

[0110] Where: w Yi is the annual runoff pollution load of the ith drainage zone, in kg / a;

[0111] 2. Calculate the total runoff pollution load w in the study area based on the annual runoff pollution load of each drainage zone Y0 , the calculation formula is as follows:

[0112]

[0113] (VI) Calculate the total annual runoff volume of the study area based on the annual runoff volume of each drainage zone, that is, the total annual runoff volume discharged into the receiving water body in the study area, specifically:

[0114] 1. Calculate the annual average runoff flow of each drainage zone according to the annual runoff volume of each drainage zone. The calculation formula is as follows:

[0115]

[0116] Where: q Yi is the annual average runoff flow from the ith drainage area to the receiving water body, and T is the dry year or wet year;

[0117] 2. Calculate the total annual runoff flow of the study area based on the annual average runoff flow of each drainage zone. The calculation formula is as follows:

[0118]

[0119] Where: q Y0 is the total annual runoff flow of the study area, that is, the total annual runoff flow discharged into the receiving water body in the study area, in m 3 / s.

[0120] S2: Estimation of point source pollution load, specifically: according to the average flow of point source pollution sewage and the point source pollution emission concentration of each drainage zone, the total annual point source pollution load discharged into the receiving water body in the study area is calculated;

[0121] The total amount of point source pollution load w in the study area in S2 w0 The methods for obtaining include:

[0122] (I) Average flow rate of point source pollution sewage in the i-th drainage area q wi The calculation formula is as follows:

[0123]

[0124] Where: q wi is the average sewage flow rate of the ith drainage area, in m 3 / s; γ i is the water consumption index of different land use types in the i-th drainage zone, in m 3 / hm 2 ·d, water consumption indicators for different types of land use refer to Table 2; a i is the area of ​​different land types in the ith drainage zone, according to the statistical table 2 of the code category of "Urban Land Classification and Planning and Construction Land Standard" GB50137, the area of ​​different types of land in the drainage zone, the unit is hm 2 ; β is the sewage discharge coefficient;

[0125] Table 2 Water consumption indicators for different types of land use (m 3 / hm 2 d)

[0126]

[0127]

[0128] (II) Annual sewage flow q in the study area w0 The calculation formula is as follows:

[0129]

[0130] (III) Annual point source pollution load w for the i-th drainage area wi The calculation formula is as follows:

[0131] w wi =10 -3 ·q wi ·c wi T;

[0132] Where: c wi is the point source pollution emission concentration of the i-th drainage zone, which is determined according to the nature of the point source pollution and the current national standards or local standards, in mg / L;

[0133] (IV) Total point source pollution load w in the study area w0 The calculation formula is as follows:

[0134]

[0135] The unit is kg / a.

[0136] S3: Estimation of the total amount of pollution load, specifically: Calculate the total amount of pollution load w0 in the study area based on the drainage system, point source discharge outlet location, total amount of runoff surface source pollution load and total amount of point source pollution load in the study area;

[0137] The drainage system includes separate system and combined system;

[0138] If the drainage system in the study area is a diversion system, and there is no single point source flowing into the receiving water body in the study area, and the influent pollution load only considers the runoff surface source pollution load, then the calculation formula for the total flow of rainwater and sewage discharged into the receiving water body, the pollutant concentration, and the total pollution load in the study area is as follows:

[0139] q0=q Y0 ;

[0140] c0=EMC;

[0141] w0=w Y0 ;

[0142] Where: q0 is the total flow of rainwater and sewage discharged into the receiving water body in the study area, q Y0 is the total annual runoff discharged into the receiving water body in the study area; c0 is the pollutant concentration discharged into the receiving water body in the study area; w0 is the total pollution load discharged into the receiving water body in the study area;

[0143] If the drainage system of the study area is a diversion system, and there is a separate point source discharge from the sewage treatment plant outlet into the receiving water body in the study area, and the influent pollution load considers the runoff surface source pollution load and the point source pollution load, then the calculation formula for the total flow of rainwater and sewage discharged into the receiving water body, the pollutant concentration and the total pollution load of the study area are as follows:

[0144] q0=q Y0 +q w0 ;

[0145]

[0146] w0=w Y0 +w w0 ;

[0147] If the drainage system in the study area is a combined system, and there is no separate point source discharge from the sewage treatment plant in the study area into the receiving water body, then the runoff non-point source pollution load and point source pollution are discharged into the receiving water body in the form of combined system overflow pollution. The calculation formula for the total flow of rainwater and sewage discharged into the receiving water body, the pollutant concentration, and the total pollution load in the study area are as follows:

[0148] q0=q Y0 -mq w0 ;

[0149] c0=EMC;

[0150] w0=(q Y0 -mq w0 )EMC;

[0151] Where: m is the interception multiple;

[0152] If the drainage system of the study area is a combined system, and there is a separate point source discharge from the sewage treatment plant outlet into the receiving water body in the study area, and the influent pollution load considers the overflow pollution of the combined system and the point source pollution load, then the calculation formula for the total flow of rainwater and sewage discharged into the receiving water body, the pollutant concentration and the total pollution load of the study area are as follows:

[0153] q0=q Y0 +q w0 ;

[0154]

[0155] w0=(q Y0 -mq w0 )·EMC+(m+1)·q w0 ·C w .

[0156] S4: Feasibility assessment of land use planning scheme, specifically: use the zero-dimensional model to calculate the allowable discharge level M of the receiving water body in the study area. If w0>M, it means that the land use planning scheme in the study area does not meet the rigid constraint conditions based on the water environment quality as the goal, and enter S5; if w0≤M, it means that the land use planning scheme is feasible;

[0157] In S4, when the receiving water body in the study area is a river, the calculation formula for the allowable discharge level of the receiving water body is as follows:

[0158] M=31.54(Q0+q0)*S s -31.54Q0S0-L;

[0159] Where: S s is the water quality control target of the receiving water body, in mg / L; S0 is the pollutant concentration of the upstream water of the receiving water body, in mg / L; Q0 is the upstream water flow of the receiving water body, in m 3 / s; q0 is the total flow of rainwater and sewage discharged into the receiving water body in the study area, unit is m 3 / s; L is the safety boundary value, unit is t / a; K is the first-order attenuation rate constant of the receiving water body, unit is 1 / d; V is the average volume of water in the receiving water body, unit is m 3 ; The unit of M is t / a;

[0160] When the receiving water body in the study area is a lake, the calculation formula for the allowable discharge level of the receiving water body is as follows:

[0161] M=31.54(Q0+q0)*S s +KS s V-3.65×10 -4 Q0S0-L.

[0162] S5: Optimization of land use planning scheme, specifically: Taking the pollutant concentration of receiving water as the target, first determine the objective function C of the receiving water in the study area i , if C i ≤S S This means that the receiving water body of the i-th drainage zone and its upstream land use planning scheme are feasible, and the i-th drainage zone still has sufficient pollution absorption capacity; if C i >S S The linear programming function is used to construct an optimal allocation model, and the land use planning scheme of the receiving water body of the i-th drainage zone and the high-risk drainage zone in its upstream is adjusted through the optimal allocation model, and S2 is returned.

[0163] In S5, a zero-dimensional model is used to generalize the impact of pollution load inflow on the receiving water body. When the receiving water body of the drainage partition is a river, the calculation formula of the objective function is as follows:

[0164]

[0165] Where: N0 is the water quality concentration of the upstream water of the receiving water body;

[0166] When the receiving water body of the drainage zone is a lake, the calculation formula of the objective function is as follows:

[0167]

[0168] When the linear programming function in S5 cannot obtain a feasible solution for the planned construction land development plan, the sponge city construction concept is used to control the non-point source pollution in the study area during the development of the drainage zone; the initial rainwater interception facilities are used to reduce mid-course or the end-of-pipe treatment facilities such as storage ponds and artificial wetlands are used to reduce runoff pollutants to meet the water environment goals.

[0169] The expression of the optimized allocation model in S5 is as follows:

[0170] When the receiving water body of the drainage zone is a river:

[0171]

[0172]

[0173] Where: C E is the pollutant concentration of the upstream water of the river; Q E is the water flow from the upstream of the river; EMC i is the average concentration of runoff pollutants in the ith drainage area; F Ii is the area of ​​Class I land in the i-th drainage zone, F IIi is the area of ​​Class II land in the ith drainage zone, F IIIi is the land area of ​​Class III in the i-th drainage zone, F IVi is the area of ​​Class IV land in the i-th drainage zone; θ max is the upper limit of the total area of ​​each type of land use in the i-th drainage zone, θ min The lower limit of the area of ​​each type of land use in the i-th drainage zone;

[0174] When the receiving water body of the drainage zone is a lake:

[0175]

[0176] Among them: F max is the upper limit of the area of ​​each type of land in the i-th drainage zone, F min It is the lower limit of the area of ​​each type of land use in the i-th drainage zone.

[0177] According to the "Urban Land Classification and Planning and Construction Land Standards" GB 50137, the proportion of residential land R, public administration and public service facilities land A, industrial land M, road and transportation facilities land S, and green land (G1+G2) in urban construction land should be 25%~40%, 5%~8%, 15%~30%, 10%~25%, and 10%~15% respectively.

[0178] This embodiment, by constructing a functional relationship between planned construction land and water environment quality, can not only quickly determine whether the land planning scheme meets the rigid constraints based on water environment goals, but also accurately trace the corresponding high-risk drainage zones, which is of guiding significance for evaluating the feasibility of planned land use schemes; moreover, based on the construction of linear programming functions, a feasible solution for the development plan of planned construction land is obtained, which provides data support for the optimization of land planning schemes, avoids water environment problems caused by urban development from the planning level, reduces the investment in subsequent water environment remediation projects, makes urban land use planning more reasonable, and further implements the urban development strategy of determining cities, land, people, and production based on water.

[0179] This embodiment also includes a computer-readable storage medium on which computer program instructions are stored. When the computer program instructions are executed by a processor, the above land planning optimization method based on water environment quality constraints is implemented.

[0180] It should be noted that the device embodiments described above are merely schematic, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. In addition, in the accompanying drawings of the device embodiments provided by the present invention, the connection relationship between the modules indicates that there is a communication connection between them, which may be specifically implemented as one or more communication buses or signal lines. A person of ordinary skill in the art may understand and implement it without paying any creative effort.

[0181] This embodiment also includes an electronic device, including: at least one processor, at least one memory, and computer program instructions stored in the memory. When the computer program instructions are executed by the processor, the mining slope stability analysis method as described above is performed.

[0182] Exemplarily, the computer program may be divided into one or more modules / units, which are stored in the memory and executed by the processor to implement the present invention. The one or more modules / units may be a series of computer program instruction segments capable of implementing specific functions, which are used to describe the execution process of the computer program in the electronic device.

[0183] The electronic device may be a computing device such as a mobile phone, a desktop computer, a notebook, a PDA, a cloud server, etc. The electronic device may include, but is not limited to, a processor and a memory. For example, the electronic device may also include an input / output device, a network access device, a bus, etc.

[0184] The processor may be a central processing unit (CPU), other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc. The processor is the control center of the electronic device, and uses various interfaces and lines to connect various parts of the entire electronic device.

[0185] The memory can be used to store the computer program and / or module, and the processor implements the computer program by running or executing the computer program and / or module stored in the memory, and calling the data stored in the memory. The memory can mainly include a program storage area and a data storage area, wherein the program storage area can store an operating system, an application required for at least one function (such as a sound playback function, an image playback function, etc.), etc.; the data storage area can store data created according to the use of the mobile phone (such as audio data, a phone book, etc.), etc. In addition, the memory can include a high-speed random access memory, and can also include a non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (SecureDigital, SD) card, a flash card (Flash Card), at least one disk storage device, a flash memory device, or other volatile solid-state storage devices.

[0186] Wherein, if the module / unit integrated in the electronic device is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on such an understanding, the present invention implements all or part of the processes in the above-mentioned embodiment method, and can also be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a readable storage medium, and the computer program can implement the steps of the above-mentioned various method embodiments when executed by the processor. Wherein, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, U disk, mobile hard disk, disk, optical disk, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electric carrier signal, telecommunication signal and software distribution medium, etc. It should be noted that the content contained in the computer-readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, computer-readable media do not include electric carrier signals and telecommunication signals.

[0187] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A land planning optimization method based on water environment quality constraints, characterized in that: The steps include: S1: Divide the study area into n drainage zones and number each drainage zone, denoted as i, where i is a natural number from 1 to n; and obtain the total planned construction land area and land use planning scheme of each drainage zone; S2: Estimation of the total amount of runoff non-point source pollution load, specifically: based on the total planned construction land area of ​​each drainage zone and the average concentration of runoff pollutants, the total amount of runoff non-point source pollution load in the study area is obtained; The total amount of point source pollution load is estimated as follows: the total amount of point source pollution load in the study area is calculated based on the average flow of point source pollution sewage and the point source pollution emission concentration in each drainage zone; S3: Estimation of the total amount of pollution load, specifically: Calculate the total amount of pollution load w0 in the study area based on the drainage system, point source discharge outlet location, total amount of runoff surface source pollution load and total amount of point source pollution load in the study area; S4: Feasibility assessment of land use planning scheme, specifically: use the zero-dimensional model to calculate the allowable discharge level M of the receiving water body in the study area. If w0>M, it means that the land use planning scheme in the study area does not meet the rigid constraint conditions based on the water environment quality as the goal, and enter S5; if w0≤M, it means that the land use planning scheme is feasible; S5: Optimization of land use planning scheme, specifically: Taking the pollutant concentration of receiving water as the target, first determine the objective function C of the receiving water in the study area i , if C i ≤S S This means that the receiving water body of the i-th drainage zone and its upstream land use planning scheme are feasible, and the i-th drainage zone still has sufficient pollution absorption capacity; if C i >S S The linear programming function is used to construct an optimal allocation model, and the receiving water body of the i-th drainage zone and the high-risk drainage zone land use planning scheme in its upstream are adjusted through the optimal allocation model, and S2 is returned, where: S S It is the water quality control target of receiving water bodies.

2. The land planning optimization method based on water environment quality constraints according to claim 1 is characterized in that: The method for obtaining the total amount of runoff non-point source pollution load in S2 includes: According to the relevant classification standards of land planning and the emission characteristics of major pollutants in surface runoff, the planned construction land in each drainage zone is reclassified into four types: I, II, III, and IV; Determine the comprehensive runoff coefficient φ of the drainage zone according to the planned town type and relevant standards, and obtain the average annual rainfall of the drainage zone; The average concentration of runoff pollutants is determined by calculating the average concentration of runoff pollutants EMC in the drainage zone based on area weighting. The calculation formula is as follows: Among them: F I The land area of ​​type I in the drainage zone, F II The land area of ​​Category II in the drainage zone, F III The land area of ​​Category III in the drainage zone, F IV Class IV land area for drainage zones; EMC I The average concentration of pollutants in runoff from Class I land in the drainage zone, EMC II The average concentration of pollutants in runoff from Class II land in the drainage zone, EMC III The average concentration of pollutants in runoff from Class III land use in the drainage zone, EMC IV is the average concentration of pollutants in runoff from Class IV land in the drainage zone; F is the total area of ​​planned construction land in the drainage zone; The volumetric method is used to calculate the annual runoff volume of each drainage zone. The calculation formula is as follows: 5 Yi =10H·φ·F; Where: V Yi is the annual runoff volume of the ith drainage zone; H is the annual rainfall in a wet year or a dry year; The total amount of runoff non-point source pollution load in the study area is calculated based on the annual runoff volume of each drainage zone and the average runoff pollutant concentration EMC, specifically: Calculate the annual runoff pollution load of each drainage zone using the following formula: w Yi =10 -3 ·V Yi ·EMC; Where: w Yi is the annual runoff pollution load of the ith drainage zone; The total runoff pollution load w in the study area is calculated based on the annual runoff pollution load of each drainage zone. Y0 , the calculation formula is as follows:

3. The land planning optimization method based on water environment quality constraints according to claim 2 is characterized in that: The S2 further includes: The total annual runoff volume of the study area is calculated based on the annual runoff volume of each drainage zone, specifically: The annual average runoff flow of each drainage zone is calculated based on the annual runoff volume of each drainage zone. The calculation formula is as follows: Where: q Yi is the average annual runoff of the ith drainage zone, T is the dry year or wet year; The total annual runoff flow of the study area is calculated based on the annual average runoff flow of each drainage zone. The calculation formula is as follows: Where: q Y0 is the total annual runoff volume in the study area.

4. The land planning optimization method based on water environment quality constraints according to claim 3 is characterized in that: The total amount of point source pollution load w in the study area in S2 w0 The methods for obtaining include: The average flow rate of point source pollution sewage in the i-th drainage area is q wi The calculation formula is as follows: Where: γ i is the water consumption index of different land use types in the i-th drainage zone; a i is the area of ​​different land use types in the i-th drainage zone; β is the sewage discharge coefficient; The annual sewage flow rate q in the study area w0 The calculation formula is as follows: The annual point source pollution load w of the i-th drainage area wi The calculation formula is as follows: w wi =10 -3 ·q wi ·c wi ·T; Where: c wi is the point source pollution emission concentration in the i-th drainage zone; The total amount of point source pollution load w in the study area w0 The calculation formula is as follows:

5. The land planning optimization method based on water environment quality constraints according to claim 4 is characterized in that: In S3, the drainage system includes a separate system and a combined system; If the drainage system in the study area is a diversion system, and there is no single point source flowing into the receiving water body in the study area, and the influent pollution load only considers the runoff surface source pollution load, then the calculation formula for the total flow of rainwater and sewage discharged into the receiving water body, the pollutant concentration, and the total pollution load in the study area is as follows: q0=q Y0 ; c0=EMC; w0=w Y0 ; Where: q0 is the total flow of rainwater and sewage discharged into the receiving water body in the study area, q Y0 is the total annual runoff discharged into the receiving water body in the study area; c0 is the pollutant concentration discharged into the receiving water body in the study area; If the drainage system of the study area is a diversion system, and there is a separate point source in the study area that discharges into the receiving water body, the influent pollution load considers the runoff surface source pollution load and the point source pollution load. The calculation formula for the total flow of rainwater and sewage discharged into the receiving water body, the pollutant concentration, and the total pollution load of the study area are as follows: q0=q Y0 +q w0 ; w0=w Y0 +in w0 ; If the drainage system of the study area is a combined system, and there is no separate point source in the study area that discharges into the receiving water body, the runoff non-point source pollution load and point source pollution are discharged into the receiving water body in the form of combined system overflow pollution. The calculation formula for the total flow of rainwater and sewage discharged into the receiving water body, the pollutant concentration, and the total pollution load of the study area are as follows: q0=q Y0 -mq w0 ; c0=EMC; w0=(q Y0 -mq w0 )·EMC; Where: m is the interception multiple; If the drainage system in the study area is a combined system, and there is a separate point source in the study area that discharges into the receiving water body, the influent pollution load takes into account the overflow pollution of the combined system and the pollution load of the point source. The calculation formula for the total flow of rainwater and sewage discharged into the receiving water body, the pollutant concentration, and the total pollution load of the study area are as follows: q0=q Y0 +q w0 ; w0=(q Y0 -mq w0 )·EMC+(m+1)·q w0 ·C w 。 6. The land planning optimization method based on water environment quality constraints according to any one of claims 1 to 5, characterized in that: In S4, when the receiving water body in the study area is a river, the calculation formula for the allowable discharge level of the receiving water body is as follows: M=31.54(Q0+q0)*S s -31.54Q0S0-L; Among them: S0 is the pollutant concentration of the upstream water of the receiving water body; Q0 ​​is the upstream water flow of the receiving water body; q0 is the total flow of rainwater and sewage discharged into the receiving water body in the study area; L is the safety boundary value; K is the first-order attenuation rate constant of the receiving water body; V is the average volume of water in the receiving water body; When the receiving water body in the study area is a lake, the calculation formula for the allowable discharge level of the receiving water body is as follows: M=31.54(Q0+q0)*S s +KS s V-3.65×10 -4 Q0S0-L。 7. The land planning optimization method based on water environment quality constraints according to any one of claims 1 to 5 is characterized in that: In S5, a zero-dimensional model is used to generalize the impact of pollution load inflow on the receiving water body. When the receiving water body of the drainage partition is a river, the calculation formula of the objective function is as follows: Where: N0 is the water quality concentration of the upstream water of the receiving water body; When the receiving water body of the drainage zone is a lake, the calculation formula of the objective function is as follows:

8. The land planning optimization method based on water environment quality constraints according to any one of claims 1 to 5 is characterized in that: When the linear programming function in S5 cannot obtain a feasible solution for the planned construction land development plan, the concept of sponge city construction is adopted in the drainage zoning development process to control the non-point source pollution in the study area; The expression of the optimization allocation model in S5 is as follows: When the receiving water body of the drainage zone is a river: Where: C E is the pollutant concentration of the upstream water of the river; Q E is the water flow from the upstream of the river; EMC i is the average concentration of runoff pollutants in the ith drainage area; F Ii is the area of ​​Class I land in the i-th drainage zone, F IIi is the area of ​​Class II land in the ith drainage zone, F IIIi is the land area of ​​Class III in the i-th drainage zone, F IVi is the land area of ​​Class IV in the ith drainage zone; a i is the area of ​​various types of land in the i-th drainage zone; θ max is the upper limit of the total area of ​​each type of land use in the i-th drainage zone, θ min The lower limit of the area of ​​each type of land use in the i-th drainage zone; When the receiving water body of the drainage zone is a lake: Among them: F max is the upper limit of the area of ​​each type of land in the ith drainage zone, F min It is the lower limit of the area of ​​each type of land use in the i-th drainage zone.

9. A computer-readable storage medium, characterized in that: Computer program instructions are stored thereon, and when the computer program instructions are executed by a processor, a land planning optimization method based on water environment quality constraints as described in any one of claims 1 to 8 is implemented.

10. A computer electronic device, characterized in that: include: At least one processor, at least one memory and computer program instructions stored in the memory, when the computer program instructions are executed by the processor, it is a land planning optimization method based on water environment quality constraints as described in any one of claims 1 to 8.

Images