CN113469563A - Method for determining background area of surface water environment of drainage basin - Google Patents
Method for determining background area of surface water environment of drainage basin Download PDFInfo
- Publication number
- CN113469563A CN113469563A CN202110822523.4A CN202110822523A CN113469563A CN 113469563 A CN113469563 A CN 113469563A CN 202110822523 A CN202110822523 A CN 202110822523A CN 113469563 A CN113469563 A CN 113469563A
- Authority
- CN
- China
- Prior art keywords
- index
- area
- background
- background area
- surface water
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000002352 surface water Substances 0.000 title claims abstract description 51
- 238000000034 method Methods 0.000 title claims abstract description 48
- 239000002689 soil Substances 0.000 claims abstract description 29
- 239000003864 humus Substances 0.000 claims abstract description 24
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 24
- 244000025254 Cannabis sativa Species 0.000 claims abstract description 20
- 238000009439 industrial construction Methods 0.000 claims abstract description 20
- 239000010865 sewage Substances 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims description 73
- 230000007704 transition Effects 0.000 claims description 5
- 238000012544 monitoring process Methods 0.000 abstract description 3
- 238000011156 evaluation Methods 0.000 abstract description 2
- 238000011835 investigation Methods 0.000 abstract description 2
- 239000011159 matrix material Substances 0.000 description 10
- 241000282414 Homo sapiens Species 0.000 description 8
- 230000000694 effects Effects 0.000 description 8
- 239000002131 composite material Substances 0.000 description 7
- 238000007621 cluster analysis Methods 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000003911 water pollution Methods 0.000 description 2
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- 239000011575 calcium Substances 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 239000003309 forest litter Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000002957 persistent organic pollutant Substances 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 238000012847 principal component analysis method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000013077 scoring method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Educational Administration (AREA)
- Economics (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Primary Health Care (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for determining a background area of a watershed surface water environment, which comprises the steps of dividing a watershed into sub-watershed catchment area units of different levels according to the catchment relation of a river system, using the sub-watershed catchment area units as computing units, determining index information in each computing unit, carrying out iterative discrimination on each computing unit according to the index information until the difference between the area of the background area discriminated by the current generation and the area of the background area discriminated by the previous generation is less than 5%, finishing discrimination, and using the background area discriminated by the current generation as a background area of the watershed surface water environment, wherein the index information is specifically a sewage outlet, a district-level government site, a large reservoir, a cultivated land area ratio, an industrial construction land area ratio, a population density, a high-humus soil area ratio and a forest grass ground area ratio, and scientifically determines the watershed surface water environment background area, provides technical support for researching water environment background value investigation and monitoring, influence evaluation and the like.
Description
Technical Field
The invention belongs to the technical field of water environment protection, and particularly relates to a method for determining a background area of a watershed surface water environment.
Background
The concept of the surface water environment background value is derived from geochemical background concentration, and the meaning of the concept is the normal concentration level of elements in a surface water environment medium under natural conditions, so that the surface water environment background value is defined as the chemical components and the content of the surface water body in an area which is not or little affected by artificial activities, namely the original content of the water environment elements and the normal value of the water environment quality distribution of the water environment elements under the condition that the water environment elements are not affected by pollution, and the surface water environment background value is generated by natural environment conditions in a certain space range and area, therefore, the concept of 'the surface water environment background area' is provided and defined, and the concept refers to the certain space range and area which are not or little affected by human activities.
However, in the prior art, the definition and division of the surface water environment background area are only expressed as an area which is not or less affected by human beings in principle, the range division of the surface water environment background area lacks standardization and operation rigor, the selection of the surface water environment background area is artificially and subjectively large, all areas subjectively judge the range and distribution of the background area according to self understanding, the influence degrees of the surface water environment background in the areas are respectively one word, the artificial pollution degree of the surface water environment in the areas affected by the background value is difficult to objectively evaluate and examine, and interference is caused to reasonably formulating water pollution prevention and control measures.
For example, in northeast regions, due to organic erosion of forest land sources, river water humus is far higher than that of a common river, organic pollutant projects are complex, background concentrations of surface water such as chemical oxygen demand, permanganate index and ammonia nitrogen index are higher all the year round, the problem of the background value of the surface water environment is extremely outstanding and typical in China, the influence range is wide and the influence degree is large, the high background value of the surface water in the northeast regions is originated from river source regions and is formed by loss of forest litter, soil humus and the like and is transported to the downstream, sampling point distribution in background value monitoring research is based on subjective judgment, and scientific definition of the background region of the surface water environment in the northeast regions is still blank.
Therefore, how to determine the background area of the surface water environment of the drainage basin, thereby determining the water environment management control target, making water pollution prevention and treatment measures and perfecting the water environment quality standard technical system in China is a technical problem to be solved by technical personnel in the field.
Disclosure of Invention
The invention aims to determine a background area of a surface water environment of a drainage basin, and provides a method for determining the background area of the surface water environment of the drainage basin.
The technical scheme of the invention is as follows: a method for determining background areas of surface water environments of drainage basins comprises the following steps:
s1, dividing the watershed into sub watershed area units of different levels according to the river water system catchment relation, and taking the sub watershed area units as calculation units;
s2, determining index information in each computing unit, wherein the index information is artificial activity information and natural condition information in the computing unit;
s3, performing iterative discrimination on each calculation unit according to the index information until the difference between the area of the background region discriminated by the current generation and the area of the background region discriminated by the previous generation is less than 5 percent, and finishing discrimination;
and S4, taking the background area determined by the current generation as the watershed surface water environment background area.
Further, the index information is a sewage outlet, a district-county-level government location, a large reservoir, a cultivated land area ratio, an industrial construction land area ratio, population density, a high-humus soil area ratio and a forest grass land area ratio.
Further, the sewage draining outlet, the district-county-level government site and the large reservoir are used as limiting indexes, and the high-humus soil area ratio, the forest and grass land area ratio, the population density, the cultivated land area ratio and the industrial construction land area ratio are used as constraint indexes.
Further, the step S3 specifically includes the following sub-steps:
s31, establishing a comprehensive index according to the constraint index;
s32, judging each calculation unit through the index information and the comprehensive index to obtain 9 judgment results of the corresponding calculation unit;
and S33, judging whether the corresponding computing unit is a background area according to the 9 judging results, if so, finishing the judgment of the corresponding computing unit, if not, dividing the corresponding computing unit to obtain a plurality of next-stage computing units, judging all the next-stage computing units and the computing units judged as the background areas through the index information and the comprehensive index, finishing the judgment until the difference between the area of the currently-judged background area of the surface water and the area of the previously-judged background area is less than 5%, and taking the last-generation judged background area as the background area of the watershed surface water.
Further, the 9 discrimination results include 3 restriction index discrimination results, 5 restriction index discrimination results, and a composite index discrimination result, where the restriction index discrimination result is specifically whether there is a corresponding restriction index, the restriction index discrimination result is specifically a category to which the corresponding restriction index belongs, the categories include a first category, a second category, and a third category, the first category is a background region, the second category is a non-background region, the third category is a transition region, and the composite index discrimination result is specifically a composite index value of the corresponding restriction index.
Further, the step S33 of determining whether the corresponding computing unit is a background area according to the 9 determination results specifically includes the following sub-steps:
s331, judging whether one or more than one of the 3 limiting index judging results have corresponding limiting indexes, if so, judging the corresponding calculating unit as a non-background area, and if not, executing a step S332;
s332, judging whether one or more corresponding constraint indexes in the 5 constraint index judgment results are of the second type, judging the corresponding calculation units as non-background areas if the corresponding constraint indexes in the results are of the second type, and executing the step S333 if the corresponding constraint indexes in the results are not of the second type;
and S333, judging whether the comprehensive index value in the comprehensive index judgment result is smaller than a preset threshold value, and judging the corresponding calculation unit as a non-background area if the comprehensive index value in the comprehensive index judgment result is smaller than the preset threshold value.
Further, the integrated index value is determined by the following formula:
wherein I is the above-mentioned combination index, wiIs the index weight, y, of the i-th calculation unitiThe normalized index value of the ith calculation unit;
wherein the normalized index values include a normalized index value of a forward index and a normalized index value of a reverse index, the forward index including a high humus soil area fraction and a forest grass area fraction, the reverse index including the population density, the plough area and the industrial construction area fraction, the normalized index values of the forward index being determined by the following formula:
in the formula, yiNormalized index value, x, of the forward index of the i-th calculation unitiIs the index value of the x index in the ith calculation unit, xminThe minimum index value of x indexes in all calculation units of the current stage, xmaxFor the current stage allCalculating the maximum index value of the x index in the unit;
the normalized index value of the reverse index is 1-yi。
Compared with the prior art, the invention has the following beneficial effects:
(1) the method comprises dividing the river basin into sub-basin catchment area units of different levels according to the river water catchment relation, determining index information in each calculation unit by using the sub-basin catchment area units as calculation units, iteratively judging each computing unit according to the index information until the difference between the area of the background area judged by the current generation and the area of the background area judged by the previous generation is less than 5 percent, finishing the judgment, taking the background area judged by the current generation as the background area of the surface water environment of the drainage basin, the index information specifically comprises a sewage discharge outlet, a district-county-level government location, a large reservoir, a cultivated land area ratio, an industrial construction land area ratio, population density, a high-humus soil area ratio and a forest grass land area ratio, the determination of the surface water environment background area of the drainage basin is scientifically realized, and technical support is provided for carrying out water environment background value investigation monitoring, influence evaluation research and the like.
(2) The non-background area is determined by the limiting index, the cluster analysis method and the comprehensive index, and the background area of the surface water environment of the watershed can be determined more accurately.
Drawings
Fig. 1 is a schematic flow chart of a method for determining a background area of a surface water environment of a watershed according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
The application provides a method for determining a background area of a watershed surface water environment, which is used for determining the background area, and fig. 1 shows a schematic flow chart of the method for determining the background area of the watershed surface water environment, which is provided by the embodiment of the application, and the method comprises the following steps:
and step S1, dividing the river basin into sub-river basin catchment area units of different levels according to the river water system catchment relation, and taking the sub-river basin catchment area units as calculation units.
Specifically, the sub-basin catchment area units of different levels are divided according to the river water catchment relation, and the number of samples of the calculation units is more than 20 as the calculation units for judgment.
And step S2, determining index information in each calculation unit.
In the embodiment of the application, the index information is a sewage outlet, a district-level government site, a large reservoir, a cultivated land area ratio, an industrial construction land area ratio, population density, a high-humus soil area ratio and a forest grass land area ratio.
Specifically, the technical scheme of the application is constructed from two aspects of human activity intensity and water environment background value source intensity, including 8 indexes of 2 classes (natural conditions and human activities) and 2 levels (restrictive and restrictive), wherein the 8 indexes are respectively high-humus soil area ratio, forest and grass land area ratio, population density, cultivated land area ratio, industrial construction land area ratio, sewage discharge outlet, district and county level government site and large reservoir, and the indexes are obtained based on a forming mechanism of a surface water environment background value of a set area given by a river basin environment. The human activity intensity is represented by 6 indexes, namely 6 indexes of population density, cultivated land area ratio, industrial construction land area ratio, sewage discharge outlet, district-county-level government location and large reservoir, and the indexes belong to reverse index types. 2 indexes of the high-humus soil area ratio and the forest grass land area ratio are forward index types, the size of the humus loss source in the natural environment is represented, and the influence degree of the natural condition on the water environment background value is disclosed. Wherein, the index of the sewage outlet, the district and county level government location and the large reservoir 3 is a restrictive index, and a ticket denial system is adopted when the background area is identified; the area ratio of the other 5 high-humus soil, the area ratio of forest and grass, the area ratio of cultivated land, the area ratio of industrial construction land and the population density are restrictive indexes, and whether the area is a background area is judged according to respective threshold values.
The calculation and determination of the 8 index values are all relative values by taking a certain space range as a statistic and calculation unit, usually dividing the calculation unit according to a drainage basin water system, wherein the numerator is a natural environment element value of a specified year or a specified period and can reflect natural condition states at different times; the denominator is the area of a calculation unit obtained by dividing according to the catchment area of the river system, and whether the calculation unit is a background area is quantitatively judged by adopting an objective method.
The high-humus soil area ratio is the ratio of the sum of the soil areas of dark brown soil, meadow soil, brown coniferous forest soil, black soil, white clay, black calcium soil and the like to the catchment area of the water functional area, and the calculation formula is the area of each type of high-humus soil/the area of a calculation unit.
The area occupation ratio of the forest and grass land is the proportion of the sum of the areas of the forest land, the grassland and the marshland land to the catchment area of the water functional area, and the calculation formula is (the area of the forest land, the grassland and the marshland land)/the area of the unit.
The proportion of the cultivated land area is the proportion of the cultivated land area to the catchment area of the water functional area, and the calculation formula is the cultivated land area/the calculation unit area.
The area ratio of the industrial construction land is the ratio of the industrial and mining storage land to the catchment area of the water functional area, and the calculation formula is the industrial and mining storage land area/the calculation unit area.
The population density is the number of population living on the land in unit area, and the calculation formula is the population/calculation unit area in the area.
In the embodiment of the application, the sewage draining outlet, the district and county level government and the large reservoir are used as limiting indexes, and the high-humus soil area ratio, the forest and grass land area ratio, the population density, the cultivated land area ratio and the industrial construction land area ratio are used as constraint indexes.
And step S3, performing iterative discrimination on each calculation unit according to the index information until the difference between the area of the background region discriminated by the current generation and the area of the background region discriminated by the previous generation is less than 5%, and finishing discrimination.
In this embodiment, the step S3 specifically includes the following sub-steps:
the step S3 specifically includes the following sub-steps:
s31, establishing a comprehensive index according to the constraint index;
s32, judging each calculation unit through the index information and the comprehensive index to obtain 9 judgment results of the corresponding calculation unit;
and S33, judging whether the corresponding computing unit is a background area according to the 9 judging results, if so, finishing the judgment of the corresponding computing unit, if not, dividing the corresponding computing unit to obtain a plurality of next-stage computing units, judging all the next-stage computing units and the computing units judged as the background areas through the index information and the comprehensive index, finishing the judgment until the difference between the area of the currently-judged background area of the surface water and the area of the previously-judged background area is less than 5%, and taking the last-generation judged background area as the background area of the watershed surface water.
Specifically, the judgment in the application is an iterative process, the non-background area is further subdivided after the current stage judges the non-background area, then the subdivided next stage calculation unit and the calculation unit judged as the background area are continuously judged through 8 indexes and comprehensive indexes until the difference between the area of the current generation judged background area and the area of the previous generation judged background area is less than 5%, the judgment is finished, and the last generation judged background area is used as the watershed surface water environment background area.
And controlling the precision of the judgment result, subdividing the calculation units according to a smaller flow domain scale for the calculation units of the non-background region, combining the calculation units which are judged as the background regions, forming a new calculation unit set, repeating and circulating the identification and judgment processes until the relative error of the difference between the areas of the background regions judged twice is less than 5%, and taking the background region judged last time as the final division result of the ground surface water environment background region, so that the technical scheme of the application can be more accurate in judging the background region.
In this embodiment of the application, the 9 discrimination results include 3 restriction index discrimination results, 5 restriction index discrimination results, and a composite index discrimination result, where the restriction index discrimination result is specifically whether there is a corresponding restriction index, the restriction index discrimination result is specifically a category to which the corresponding restriction index belongs, the categories include a first category, a second category, and a third category, the first category is a background region, the second category is a non-background region, the third category is a transition region, and the composite index discrimination result is specifically a composite index value of the corresponding restriction index.
In this embodiment of the present application, the step S33 of determining whether the corresponding computing unit is a background area according to the 9 determination results specifically includes the following sub-steps:
s331, judging whether one or more than one of the 3 limiting index judging results have corresponding limiting indexes, if so, judging the corresponding calculating unit as a non-background area, and if not, executing a step S332;
s332, judging whether one or more corresponding constraint indexes in the 5 constraint index judgment results are of the second type, judging the corresponding calculation units as non-background areas if the corresponding constraint indexes in the results are of the second type, and executing the step S333 if the corresponding constraint indexes in the results are not of the second type;
and S333, judging whether the comprehensive index value in the comprehensive index judgment result is smaller than a preset threshold value, and judging the corresponding calculation unit as a non-background area if the comprehensive index value in the comprehensive index judgment result is smaller than the preset threshold value.
In the embodiment of the present application, the comprehensive index value is determined by the following formula:
wherein I is the above-mentioned combination index, wiIs the index weight, y, of the i-th calculation unitiThe normalized index value of the ith calculation unit;
wherein the normalized index values include a normalized index value of a forward index and a normalized index value of a reverse index, the forward index including a high humus soil area fraction and a forest grass area fraction, the reverse index including the population density, the plough area and the industrial construction area fraction, the normalized index values of the forward index being determined by the following formula:
in the formula, yiNormalized index value, x, of the forward index of the i-th calculation unitiIs the index value of the x index in the ith calculation unit, xminThe minimum index value of x indexes in all calculation units of the current stage, xmaxThe maximum index value of the x index in all the calculation units of the current stage is obtained, wherein the index value is an actual value of the corresponding constraint index;
the normalized index value of the reverse index is 1-yi。
Specifically, 5 restrictive indexes, namely high-humus soil area ratio, forest grass land area ratio, population density, cultivated land area ratio and industrial construction land area ratio, are classified and grouped respectively by applying a statistical cluster analysis method, and the cluster is divided into three types: the first category is a region where human activity is low and background values are sourced large (background region), the second category is a region where human activity is high and background values are sourced small (non-background region), and the third category is a transition region between the first and second categories. And the second type human motion is judged as a non-background area by the calculating unit with high intensity and low intensity of the background value source.
Specifically, a comprehensive index is established for 5 restrictive indexes (high humus soil area ratio, forest and grass land area ratio, population density, cultivated land area ratio and industrial construction land area ratio), and a non-background area is identified according to the size of the comprehensive index value.
The calculation formula of the comprehensive index I is as follows:
in the formula, wiIs the index weight, y, of the i-th calculation unitiThe numerical range of the comprehensive index is 0-1, and the smaller the value of I is, the more the value tends to be in a non-background area.
And respectively calculating the comprehensive indexes of each calculating unit at the current stage, arranging the comprehensive indexes from large to small, drawing a comprehensive index curve graph of each calculating unit, taking the inflection point of the curve as a preset threshold, and judging the calculating unit smaller than the threshold as a non-background area.
Further, the 5 restrictive indexes are subjected to dimensionless normalization processing. Wherein two indexes of high humus soil area ratio and forest grass land area ratio are positive indexes; three indexes of population density, arable land area ratio and industrial construction land area ratio are reverse indexes. The normalization processing mode of the forward index is shown in the following formula, and the normalization processing mode of the reverse index is 1-yi。
The weights of the 5 constraint indexes are determined by an analytic hierarchy process. Firstly, a judgment matrix is constructed, the judgment matrix is the comparison of the relative importance of each index in the layer relative to one index in the previous layer, and the element a of the judgment matrixijGiven the 1-9 scale of Santy, the numerical size is used to measure the degree of importance of each index relative to the other in the index system, where the number 1 indicates that two indices are equally important and the number 9 indicates that the former is extremely important relative to the latter. Secondly, a square root method is adopted to calculate the maximum characteristic root lambda of the judgment matrix, the characteristic vector corresponding to the maximum characteristic root is normalized and then serves as a weight vector W, and each element of the weight vector W is the weight of each index. Finally, calculating the consistency index CI of the judgment matrix, and obtaining the consistencyAfter passing the inspection, the result shows that the calculated weight result is usable, and the weight value of each index information in the application is specifically shown in the following table 1:
TABLE 1
And (4) solving a union set of the calculation units which are judged by the three methods and are used as non-background areas to obtain a set of the non-background areas, and judging the rest calculation unit sets to be used as surface water environment background areas, namely background areas.
In addition, in a specific application scenario, the selection and composition of the indexes can be added or changed according to actual conditions, for example, the indexes such as land coverage, hydrogeology, litter, microorganisms, economic structures, water and soil loss, urbanization level, road density, pollutant discharge amount and the like are added, and the method for determining the index weight can be divided into three types, namely subjective weighting method, objective weighting method and combined weighting method, which are commonly used and include expert survey method, analytic hierarchy process, binomial coefficient method, cyclic ratio scoring method, least squares method, principal component analysis method, entropy value method, dispersion and mean square deviation method, multi-objective planning method, multiplication integration method, addition integration method and the like, and the protection range of the application is not affected.
Taking the northeast region as an example, firstly dividing the watershed into sub-watershed area units of different levels according to the river water catchment relation in the northeast region, then constructing an index system based on a watershed environment endowment and a surface water environment background value forming mechanism, and determining 3 limiting indexes: drain, district and county level government place and large reservoir, and 5 restraint indexes: the method comprises the following steps of judging each calculation unit according to a high-humus soil area ratio, a forest and grass land area ratio, population density, a cultivated land area ratio, an industrial construction land area ratio and a comprehensive index established according to 5 constraint indexes by 8 indexes and the comprehensive index, and obtaining 9 judgment results, wherein the judgment method comprises the following steps:
1. the limiting index is a non-background area by adopting a vote rejection method, and the corresponding computing unit is judged as the non-background area as long as the corresponding computing unit contains one or more limiting indexes;
2. the 5 constraint indexes are respectively judged through a cluster analysis method, and the judgment results are as the first class, the second class and the third class, wherein the first class is a background area, the second class is a non-background area, and the third class is a transition area;
3. establishing and judging comprehensive indexes:
the composite index value is determined by the formula:
wherein I is the above-mentioned combination index, wiIs the index weight, y, of the i-th calculation unitiThe normalized index value of the ith calculation unit;
wherein the normalized index values include a normalized index value of a forward index and a normalized index value of a reverse index, the forward index including a high humus soil area fraction and a forest grass area fraction, the reverse index including the population density, the plough area and the industrial construction area fraction, the normalized index values of the forward index being determined by the following formula:
in the formula, yiNormalized index value, x, of the forward index of the i-th calculation unitiIs the index value of the x index in the ith calculation unit, xminThe minimum index value of x index in all calculation units of the current generation, xmaxThe maximum index value of the x index in all the calculation units of the current generation is obtained;
the normalized index value of the reverse index is 1-yi。
The weights w of the 5 constraint indexes are determined by an analytic hierarchy process. Firstly, a judgment matrix is constructed, the judgment matrix is a comparison of the relative importance of each index in the layer relative to one index in the previous layer, the element aij of the judgment matrix is given by adopting a 1-9 scale method of Santy, the number is used for measuring the importance degree of each index relative to other indexes in an index system, wherein the number 1 represents that the two indexes have equal importance, and the number 9 represents that the former is extremely important relative to the latter. Secondly, a square root method is adopted to calculate the maximum characteristic root lambda of the judgment matrix, the characteristic vector corresponding to the maximum characteristic root is normalized and then serves as a weight vector W, and each element of the weight vector W is the weight of each index. And finally, calculating a consistency index CI of the judgment matrix, and after passing consistency check, indicating that the calculated weight result is available. The weight values of the indexes obtained by final determination and calculation are shown in table 1.
Then, after a generation of discrimination is performed on all the computing units, a non-background region and a background region which are discriminated by a current generation can be determined, it should be noted that the discrimination process of the present application is iterative discrimination, the non-background region which is discriminated by the current generation needs to be further subdivided to be used as a next-stage computing unit of a smaller unit, then all the next-stage computing units and the computing unit which is discriminated as the background region by the current generation are discriminated again until the difference between the area of the background region which is discriminated by the current generation or the area of the non-background region and the corresponding area which is discriminated by the previous generation is less than 5%, the discrimination process is ended, and then the background region which is discriminated by the last discrimination process is used as a final environment background region of the northeast region.
It will be appreciated by those of ordinary skill in the art that the embodiments described herein are intended to assist the reader in understanding the principles of the invention and are to be construed as being without limitation to such specifically recited embodiments and examples. Those skilled in the art can make various other specific changes and combinations based on the teachings of the present invention without departing from the spirit of the invention, and these changes and combinations are within the scope of the invention.
Claims (7)
1. A method for determining background areas of surface water environments of drainage basins is characterized by comprising the following steps:
s1, dividing the watershed into sub watershed area units of different levels according to the river water system catchment relation, and taking the sub watershed area units as calculation units;
s2, determining index information in each calculation unit;
s3, performing iterative discrimination on all the calculation units according to the index information until the difference between the area of the background region discriminated by the current generation and the area of the background region discriminated by the previous generation is less than 5%, and finishing discrimination;
and S4, taking the background area determined by the current generation as the watershed surface water environment background area.
2. The method for determining the background area of the watershed surface water environment according to claim 1, wherein the index information is a sewage outlet, a district-level government site, a large reservoir, a cultivated land area ratio, an industrial construction land area ratio, population density, a high-humus soil area ratio and a forest grass land area ratio.
3. The method for determining the background area of the watershed surface water environment according to claim 2, wherein the sewage discharge outlet, the district-county-level government site and the large reservoir are used as limiting indexes, and the high-humus soil area ratio, the forest and grass land area ratio, the population density, the cultivated land area ratio and the industrial construction land area ratio are used as limiting indexes.
4. The method for determining the background area of the watershed surface water environment as claimed in claim 3, wherein the step S3 comprises the following sub-steps:
s31, establishing a comprehensive index according to the constraint index;
s32, judging each calculation unit through the index information and the comprehensive index to obtain 9 judgment results of the corresponding calculation unit;
and S33, judging whether the corresponding computing unit is a background area according to the 9 judging results, if so, finishing the judgment of the corresponding computing unit, if not, dividing the corresponding computing unit to obtain a plurality of next-stage computing units, judging all the next-stage computing units and the computing units judged as the background areas through the index information and the comprehensive index, finishing the judgment until the difference between the area of the currently-judged background area of the surface water and the area of the previously-judged background area is less than 5%, and taking the last-generation judged background area as the background area of the watershed surface water.
5. The method for determining the background area of the watershed surface water environment according to claim 4, wherein the 9 discrimination results include 3 restriction index discrimination results, 5 restriction index discrimination results and a comprehensive index discrimination result, wherein the restriction index discrimination result is specifically whether a corresponding restriction index exists, the restriction index discrimination result is specifically a category to which the corresponding restriction index belongs, the category includes a first category, a second category and a third category, the first category is the background area, the second category is a non-background area, the third category is a transition area, and the comprehensive index discrimination result is specifically a comprehensive index value of the corresponding restriction index.
6. The method for determining the background area of the watershed surface water environment as claimed in claim 5, wherein the step S33 of determining whether the corresponding computing unit is the background area according to the 9 determination results specifically includes the following sub-steps:
s331, judging whether one or more than one of the 3 limiting index judging results have corresponding limiting indexes, if so, judging the corresponding calculating unit as a non-background area, and if not, executing a step S332;
s332, judging whether one or more corresponding constraint indexes in the 5 constraint index judgment results are of the second type, judging the corresponding calculation units as non-background areas if the corresponding constraint indexes in the results are of the second type, and executing the step S333 if the corresponding constraint indexes in the results are not of the second type;
and S333, judging whether the comprehensive index value in the comprehensive index judgment result is smaller than a preset threshold value, and judging the corresponding calculation unit as a non-background area if the comprehensive index value in the comprehensive index judgment result is smaller than the preset threshold value.
7. The method for determining the background area of the watershed surface water environment as claimed in claim 5, wherein the comprehensive index value is determined by the following formula:
wherein I is the above-mentioned combination index, wiIs the index weight, y, of the i-th calculation unitiThe normalized index value of the ith calculation unit;
wherein the normalized index values include a normalized index value of a forward index and a normalized index value of a reverse index, the forward index including a high humus soil area fraction and a forest grass area fraction, the reverse index including the population density, the plough area and the industrial construction area fraction, the normalized index values of the forward index being determined by the following formula:
in the formula, yiNormalized index value, x, of the forward index of the i-th calculation unitiIs the index value of the x index in the ith calculation unit, xminThe minimum index value of x index in all calculation units of the current generation, xmaxThe maximum index value of the x index in all the calculation units of the current generation is obtained;
the normalized index value of the reverse index is 1-yi。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110822523.4A CN113469563A (en) | 2021-07-21 | 2021-07-21 | Method for determining background area of surface water environment of drainage basin |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110822523.4A CN113469563A (en) | 2021-07-21 | 2021-07-21 | Method for determining background area of surface water environment of drainage basin |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113469563A true CN113469563A (en) | 2021-10-01 |
Family
ID=77881767
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110822523.4A Pending CN113469563A (en) | 2021-07-21 | 2021-07-21 | Method for determining background area of surface water environment of drainage basin |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113469563A (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108681974A (en) * | 2018-05-29 | 2018-10-19 | 中国水利水电科学研究院 | A kind of water resource multidimensional equilibrium allocation method towards Arid Region |
| CN110717257A (en) * | 2019-09-23 | 2020-01-21 | 中国水利水电科学研究院 | Optimization method, device, equipment and medium of distributed hydrological model |
-
2021
- 2021-07-21 CN CN202110822523.4A patent/CN113469563A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108681974A (en) * | 2018-05-29 | 2018-10-19 | 中国水利水电科学研究院 | A kind of water resource multidimensional equilibrium allocation method towards Arid Region |
| CN110717257A (en) * | 2019-09-23 | 2020-01-21 | 中国水利水电科学研究院 | Optimization method, device, equipment and medium of distributed hydrological model |
Non-Patent Citations (2)
| Title |
|---|
| MAOQING DUAN等: "A Revised Method of SurfaceWater Quality Evaluation Based on Background Values and Its Application to Samples Collected in Heilongjiang Province, China", 《WATER》 * |
| 段茂庆等: "特殊区域地表水环境背景值表征技术方法研究与应用", 《中国环境科学》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wu et al. | Evaluating uncertainty estimates in distributed hydrological modeling for the Wenjing River watershed in China by GLUE, SUFI-2, and ParaSol methods | |
| CN106339536B (en) | Comprehensive Evaluation of Water Quality based on water pollution index's method and cloud model | |
| Wu et al. | Diagnosis of wetland ecosystem health in the Zoige Wetland, Sichuan of China | |
| Can et al. | Assessing impacts of different land use scenarios on water budget of Fuhe River, China using SWAT model | |
| Al-Ansari et al. | Long term effect of climate change on rainfall in northwest Iraq | |
| Costa et al. | Clustering and forecasting of dissolved oxygen concentration on a river basin | |
| CN113887972A (en) | Comprehensive drought monitoring and evaluating method based on hydrological process | |
| Pardo et al. | An invertebrate predictive model (NORTI) for streams and rivers: sensitivity of the model in detecting stress gradients | |
| Emamgholizadeh et al. | Comparison of artificial neural networks, geographically weighted regression and Cokriging methods for predicting the spatial distribution of soil macronutrients (N, P, and K) | |
| Voyslavov et al. | Hasse diagram technique as a tool for water quality assessment | |
| Giovanardi et al. | Development of ecological classification criteria for the Biological Quality Element phytoplankton for Adriatic and Tyrrhenian coastal waters by means of chlorophyll a (2000/60/EC WFD) | |
| Wang et al. | Real-time water quality monitoring and estimation in AIoT for freshwater biodiversity conservation | |
| CN113033997A (en) | Urban water quality grade determination method, device and medium based on improved set pair analysis | |
| CN111428985B (en) | Assessment method for regional water ecological bearing capacity | |
| CN111090831A (en) | Lake area change key driving factor identification method | |
| Zhang et al. | A time-varying parameter estimation approach using split-sample calibration based on dynamic programming | |
| Tran et al. | Quantification of predictive uncertainty with a metamodel: Toward more efficient hydrologic simulations | |
| Sotomayor et al. | Implications of macroinvertebrate taxonomic resolution for freshwater assessments using functional traits: The Paute River Basin (Ecuador) case | |
| CN115409374A (en) | Method for early warning influence of urbanization difference on watershed water ecological service function | |
| Alemneh et al. | Modeling the impact of highland settlements on ecological disturbance of streams in Choke Mountain Catchment: Macroinvertebrate assemblages and water quality | |
| Liu et al. | Copula-based framework for integrated evaluation of water quality and quantity: A case study of Yihe River, China | |
| Donyaii | Evaluation of climate change impacts on the optimal operation of multipurpose reservoir systems using cuckoo search algorithm | |
| Zhang | A spatially explicit model for estimating annual average loads of nonpoint source nutrient at the watershed scale | |
| Yan et al. | Water quality evaluation based on entropy coefficient and blind number theory measure model | |
| Eslamian et al. | Monitoring and prediction of drought by Markov chain model based on SPI and new index in Isfahan |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211001 |
|
| RJ01 | Rejection of invention patent application after publication |