CN106485076A - A kind of method of basin chemical integrity assessment - Google Patents

Abstract

The invention discloses a kind of method of basin chemical integrity assessment, belong to a kind of chemical integrity appraisal procedure, including three steps, respectively step, perch to be located on river and riverbank and be correspondingly provided with some detection zones, arranging some test points in each detection zone carries out ecosystem integrity evaluation by this appraisal procedure to river, step 2, development of judgment matrix, step 3, carry out fuzzy evaluation, the present invention is the most important condition for carrying out ecosystem restoration, the effect of ecological recovery is carried out according to its evaluation result assessment Damaged River System, can provide fundamental basis for the selection of concrete control measures, be conducive to the management of river ecosystem, its evaluation result is not only the sustainable development offer decision support that river ecological builds social economy, also following to the river basins control and plan has important directive significance.

Description

A kind of method of basin chemical integrity assessment

Technical field

The present invention relates to a kind of integrity assessment method, especially a kind of method of basin chemical integrity assessment.

Background technology

River is played an important role to the development of human society, and the river in basin provides substantial amounts of light for the mankind Water resource and food, also have generating, shipping, the recharge of ground water, receive dirt and degradation of contaminant, flood discharge, landscape, amusement and recreation, Multiple Ecosystem Service such as sports.It can be seen that, river is had an important significance to the survival and development of the mankind.In recent years, Due to the high speed development of social economy and the continuous upgrading of human development activity, the river of countries in the world all receives different journeys The pollution of degree and destruction, wherein, the phenomenon destroyed by urban river ecosystem is projected the most.The ecosystem in river is not because Disconnected interference and destruction by mankind's activity, river ecological environment occur in that Cut off, water degradation, hydrologic regime change with And river morphology condition degenerate etc. a series of problems, such as, cause the degeneration of river ecosystem, the health for having a strong impact on river is sent out Exhibition.This problem, has become the biggest crisis of human society sustainable development, has caused whole world water conservancy working person and Sheng The extensive concern of state researcher.

River ecosystem is one kind of river aquatic ecosystem, is that biocoene is interacted with river environment in river Entity；River ecosystem is also the ecosystem of a flowing, together with it is closely connected with the land of surrounding, formation The ecosystem of one opposing open, it are the ties that land is contacted with ocean, rise emphatically in the material circulation in biosphere The effect that wants.River ecosystem includes land riparian ecosystem, aquatic ecosystem, related wetland and marsh etc. interior A series of subsystems, with multiple ecosystems such as filtration, channeling, shielding action, habitat function and source remittance functions System service function.

Content of the invention

The technical assignment of the present invention is for above the deficiencies in the prior art, and provides a kind of assessment of basin chemical integrity Method.

The technical solution adopted for the present invention to solve the technical problems is：A kind of method of basin chemical integrity assessment, Wherein it is divided into three below step, specific as follows:

Step one, perches to be located on river and riverbank and is correspondingly provided with some detection zones, arranges some inspections in each detection zone Measuring point；

Step 2, determines each index weights according to proportion quotiety reference table；

(1) according to comparative result, river water environmental state and riverbank habitat state are compared, it is believed that River water environmental state is slightly more important than riverbank habitat state, draws single finger thus according to classification proportion quotiety reference table Target relative importance value, and development of judgment matrix A-B；

A	B1	B2
			B1	1	3
B2	1/3	1

According to above-mentioned development of judgment matrix A-B, obtained by formula：

b₂=0.577

Can be calculated by formula：

w₂=0.250

Therefore the characteristic vector of matrix B is：W=(0.750,0.250) T；Then matrix B=[0.75 0.25]

(2) weight in river and each test point institute Testing index in the habitat of riverbank is determined according to proportion quotiety reference table, Thus weight parameter matrix A is built, and wherein the heavy weight parameter matrix for river is A1, for the heavy weight parameter of riverbank habitat Matrix is A2；

(3) initially with power multiplication or the uniformity of geometric average method or area method judgment matrix, then again to matrix one Cause property is verified；

Step 3, fuzzy evaluation；

(1) one-level fuzzy evaluation, according to chemical integrity Assessment for classification standard, evaluation index is divided into five ranks, profit Use Cauchy's distribution functionCalculate same index be located at different stage when parameter, so as to Development of judgment matrix R, the wherein judgment matrix for river are R1, and the judgment matrix for riverbank habitat is R2；

Wherein when index is in the IIth class of evaluation criterion, III class, IV class：

When index is in I class of evaluation criterion：

When index is in V class of evaluation criterion：

In formula：x_uFor the upper boundary values of corresponding level relevant parameter, x_iFor the lower border value of corresponding level relevant parameter, x is Actually measured parameter value；

N1 is the one-level fuzzy evaluation collection in river, wherein N1=A1*R1；

N2 is the one-level fuzzy evaluation collection of riverbank habitat, wherein N2=A2*R2

(2) N1 is combined with the collection of N2 and is constituted new judgment matrix N by two-level appraisement first, then allows matrix B and matrix N Multiplication is B*N, and obtained set is the result of two-level appraisement, and the rank that last maximum number is located just represents the river Degree of membership.

Optimize further：In step 3, same index be located at different stage when criteria for classification table as follows；

Criteria for classification table

Optimize further：In step 2, the uniformity of matrix is determined using geometric average method；

1. product bi of each row element of judgment matrix is calculated,

2. calculateN th Root b_i,

3. to vectorial b_i=(b₁b₂… b_n)^TMake normalized,

Then w_iFor required weight coefficient,

4. the Maximum characteristic root λ of judgment matrix is calculated_max,

5. Mode of Level Simple Sequence and consistency check,

If the element of judgment matrix A meets following formula:

Matrix A is then claimed to have uniformity；

Whether matrix A is verified with uniformity：Using below equation

In formula：The random Consistency Ratio of CR- judgment matrix；

N is the exponent number of judgment matrix；CI is the general coincident indicator of judgment matrix；RI is the mean random of judgment matrix Coincident indicator；Occurrence is shown in Table；

Rank matrix RI value

6. synthetic weight of each layer element to aims of systems is calculated, and is ranked up, as random Consistency Ratio CR=CI/ During RI ＜ 0.10, it is believed that Mode of Level Simple Sequence result has satisfied uniformity, otherwise, being considered as the preliminary judgment matrix that sets up is Unsafty, need to readjust element scale value, till there is satisfied uniformity.

Optimize further：18 test points are provided with same detection zone, wherein have 9 points to be located in river, respectively PH, electrical conductivity, dissolved oxygen, COD, BOD, total nitrogen, total phosphorus, heavy metal in water, Volatile Phenols, in addition 9 test points are positioned at riverbank In habitat, respectively soil organic matter content, total nitrogen content of soil, Soil total nitrogen, total potassium content of soil, soil are effective Nitrogen content, soil available phosphorus content, available potassium in soils content, BHC, DDT.

Optimize further：The matrix R1 is the matrix that 9 row 5 is arranged, and detects on where each row and river and riverbank habitat The index detected in area is corresponding, and each row represent the indicator layer of different stage.

Optimize further：In step 2, according to proportion quotiety reference table, this 18 point index weights tables are drawn, as follows；

Index weights table

Obtain matrix：

A1=[0.028 0.028 0.05 0.0234 0.096 0.0234 0.0112 0.082 0.136]

A2=[0.098 0.046 0.043 0.043 0.098 0.098 0.098 0.0158 0.158].

Advantages of the present invention：Ecosystem integrity evaluation is carried out to river by this appraisal procedure, is to carry out ecosystem The most important condition that system recovers, carries out the effect of ecological recovery, can administer for concrete according to its evaluation result assessment Damaged River System The selection of measure is provided fundamental basis, and is conducive to the management of river ecosystem, and its evaluation result is not only river ecological construction The sustainable development of social economy provides decision support, and also following to the river basins control and plan has important guidance meaning Justice.

Description of the drawings

Fig. 1 is the block diagram of assessment result of the present invention.

Specific embodiment

A kind of method of basin chemical integrity assessment, is wherein divided into three below step, specific as follows:

Step one, perches to be located on river and riverbank and is correspondingly provided with some detection zones, arranges some inspections in each detection zone Measuring point；In whole piece river, 17 detection zones are wherein provided with altogether, 18 test points are provided with each detection zone, wherein have 9 Individual point is located in river, respectively PH, electrical conductivity, dissolved oxygen, COD, BOD, total nitrogen, total phosphorus, heavy metal in water, Volatile Phenols, Its detection the data obtained such as following table：

Tables of data 1

Other 9 test points are located in the habitat of riverbank, respectively soil organic matter content, total nitrogen content of soil, soil Content of tatal phosphorus, total potassium content of soil, soil available nitrogen content, soil available phosphorus content, available potassium in soils content, BHC, drop DDT, its detection the data obtained is following table:

Tables of data 2

Step 2, determines each index weights according to proportion quotiety reference table；Wherein

Proportion quotiety reference table

(1) according to comparative result, river water environmental state and riverbank habitat state are compared, it is believed that River water environmental state is slightly more important than riverbank habitat state, draws single finger thus according to classification proportion quotiety reference table Target relative importance value, and development of judgment matrix A-B；

A	B1	B2
			B1	1	3
B2	1/3	1

According to above-mentioned development of judgment matrix A-B, obtained by formula：

b₂=0.577

W can be calculated by formula_i

w₂=0.250

Therefore the characteristic vector of matrix B is：W=(0.750,0.250) T；Then matrix B=[0.75 0.25]

(2) determine the power in river and each test point institute Testing index in the habitat of riverbank in the same manner according to proportion quotiety reference table Weight such as following table,

Index weights table

Thus weight parameter matrix A is built, and wherein the heavy weight parameter matrix for river is A1, for riverbank habitat Weight weight parameter matrix is A2, and：

A1=[0.028 0.028 0.05 0.0234 0.096 0.0234 0.0112 0.082 0.136]

A2=[0.098 0.046 0.043 0.043 0.098 0.098 0.098 0.0158 0.158]

(3) because second-order matrix always has uniformity, this judgment matrix need not be tested, it is therefore desirable to right The uniformity of weight parameter matrix is judged again, can adopt power multiplication, geometric average method, the uniformity of area method judgment matrix, Here is judged using geometric average method：

1. product bi of each row element of judgment matrix is calculated,

2. calculateN th Root b_i,

3. to vectorial b_i=(b₁b₂… b_n)^TMake normalized,

Then w_iFor required weight coefficient,

4. the Maximum characteristic root λ of judgment matrix is calculated_max,

5. Mode of Level Simple Sequence and consistency check,

If the element of judgment matrix A meets following formula:

Matrix A is then claimed to have uniformity；

Then, whether matrix A is verified with uniformity：Using below equation

In formula：The random Consistency Ratio of CR- judgment matrix；

N is the exponent number of judgment matrix；CI is the general coincident indicator of judgment matrix；RI is the mean random of judgment matrix Coincident indicator；Occurrence is shown in Table.

Rank matrix RI value

6. synthetic weight of each layer element to aims of systems is calculated, and is ranked up, as random Consistency Ratio CR=CI/ During RI ＜ 0.10, it is believed that Mode of Level Simple Sequence result has satisfied uniformity, otherwise, being considered as the preliminary judgment matrix that sets up is Unsafty, need to readjust element scale value, till there is satisfied uniformity.

Step 3, fuzzy evaluation；

(1) one-level fuzzy evaluation, according to chemical integrity Assessment for classification standard scale, evaluation index is divided into five ranks, Using Cauchy's distribution functionCalculate same index be located at different stage when parameter, from And development of judgment matrix R, the wherein judgment matrix for river are R1, the judgment matrix for riverbank habitat is R2；

When index is in the IIth class of evaluation criterion, III class, IV class：

When index is in I class of evaluation criterion：

When index is in V class of evaluation criterion：

In formula：x_uFor the upper boundary values of corresponding level relevant parameter, x_iFor the lower border value of corresponding level relevant parameter, x is Actually measured parameter value；

Assessment for classification standard scale

By taking the detected value at Ford shop as an example, through calculating

N1 is the one-level fuzzy evaluation collection in river, wherein N1=A1*R1；

N1=[0.2957 0.1610 0.5335 0.2735 0.4115]

N2 is the one-level fuzzy evaluation collection of riverbank habitat, wherein N2=A2*R2

N2=[0.1318 0.1422 0.3130 0.2130 0.1467]

(2) N1 is combined with the collection of N2 and is constituted new judgment matrix N by two-level appraisement first, then allows matrix B and matrix N Multiplication is M=B*N,

Now

M=[0.2547 0.1563 0.4785 0.2584 0.3453]

It is maximum positioned at tertial number 0.4785 that matrix from obtained by M is can be seen that in the matrix, therefore should The row that number is located represent the degree of membership in river at the detection zone, and therefore the evaluation result at Ford shop is general.

Because representing same index parameter per a line, therefore from the beginning of the first row in R1, last column in R2, per Row successively represent PH, electrical conductivity, dissolved oxygen ... BHC, DDT, for example：The 0.743 of the first row first row in R1 matrix Calculating process.

1. it is the PH index of river water environmental state as the position is corresponding, Ford shop PH number is understood according to tables of data 1 It is worth for 8.117, from Assessment for classification standard scale, when 8.117 in the first estate, chooses formula such as (2-3)；Now a= x_u=9 (maximums)；That is α=4/ (9-6) 2=0.444；X value is for actually measuring the number of acquisition in tables of data 1 According to herein referring to pH value, i.e., 8.117 in the river water environmental of Ford shop；

2. a and α is substituted into formula (2-1) to obtain Approximate 0.743.

Embodiments of the invention are the foregoing is only, the scope of the claims of the present invention is not thereby limited, every using this Equivalent structure or equivalent flow conversion that bright description is made, or directly or indirectly it is used in other related technology necks Domain, is included within the scope of the present invention.

Claims (6)

1. a kind of method that basin chemical integrity is assessed, it is characterised in that：It is divided into three below step, specific as follows:

Step one, perches to be located on river and riverbank and is correspondingly provided with some detection zones, arranges some test points in each detection zone；

Step 2, determines each index weights according to proportion quotiety reference table；

(1) according to comparative result, river water environmental state and riverbank habitat state are compared, it is believed that river Water environment state is slightly more important than riverbank habitat state, draws single index thus according to classification proportion quotiety reference table Relative importance value, and development of judgment matrix A-B；

A B1 B2 B1 1 3 B2 1/3 1

According to above-mentioned development of judgment matrix A-B, obtained by formula：

$b_1=\underset{j=1}{\overset{n}{\Π}}{a}_{1j}={(1\×3)}^{\frac{1}{2}}=1.732$

b₂=0.577

Can be calculated by formula：

$w_1=\frac{b_1}{{\mathrm{\Σ}}_{i=1}^2{b}_i}=\frac{1.732}{(1.732+0.577)}=0.750$

w₂=0.250

Therefore the characteristic vector of matrix B is：W=(0.750,0.250) T；Then matrix B=[0.75 0.25]

(2) weight in river and each test point institute Testing index in the habitat of riverbank is determined according to proportion quotiety reference table, thus Weight parameter matrix A is built, wherein the heavy weight parameter matrix for river is A1, for the heavy weight parameter matrix of riverbank habitat For A2；

(3) initially with power multiplication or the uniformity of geometric average method or area method judgment matrix, then uniformity again to matrix Verified；

Step 3, fuzzy evaluation；

(1) one-level fuzzy evaluation, according to chemical integrity Assessment for classification standard, evaluation index is divided into five ranks, using Ke Western distribution functionCalculate same index be located at different stage when parameter, so as to build Judgment matrix R, the wherein judgment matrix for river are R1, and the judgment matrix for riverbank habitat is R2；

Wherein when index is in the IIth class of evaluation criterion, III class, IV class：

$\begin{array}{cc}a=\frac{x_u+x_i}{2}& \mathrm{\α}=\frac{4}{{(x_u-x_i)}^2}\end{array}---(2-2)$

When index is in I class of evaluation criterion：

$\begin{array}{cc}a=x_u& \mathrm{\α}=\frac{4}{{(x_u-x_i)}^2}\end{array}---(2-3)$

When index is in V class of evaluation criterion：

$\begin{array}{cc}a=x_i& \mathrm{\α}=\frac{4}{{(x_u-x_i)}^2}\end{array}---(2-4)$

In formula：x_uFor the upper boundary values of corresponding level relevant parameter, x_iFor the lower border value of corresponding level relevant parameter, x is actual The parameter value for measuring；

N1 is the one-level fuzzy evaluation collection in river, wherein N1=A1*R1；

N2 is the one-level fuzzy evaluation collection of riverbank habitat, wherein N2=A2*R2

(2) N1 is combined with the collection of N2 and is constituted new judgment matrix N, then allows matrix B be multiplied with matrix N by two-level appraisement first That is B*N, obtained set are the result of two-level appraisement, and the rank that last maximum number is located just represents the person in servitude in the river Category degree.

2. the method that a kind of basin chemical integrity according to claim 1 is assessed, it is characterised in that：In step 3, Same index be located at different stage when criteria for classification table as follows；

Classification of assessment standard scale

3. the method that a kind of basin chemical integrity according to claim 1 is assessed, it is characterised in that：In step 2, The uniformity of matrix is determined using geometric average method；

1. product bi of each row element of judgment matrix is calculated,

$\stackrel{\‾}{b_i}=\underset{j=1}{\overset{n}{\mathrm{\Π}}}{a}_{ij},(i,j=1,2,3,...,n)---(3-1)$

2. calculateN th Root b_i,

$b_i=\sqrt[n]{\stackrel{\‾}{b_i}}$

3. to vectorial b_i=(b₁b₂… b_n)^TMake normalized,

$w_i=b_i/\underset{i=1}{\overset{n}{\Σ}}{b}_i---(3-3)$

Then w_iFor required weight coefficient,

4. the Maximum characteristic root λ of judgment matrix is calculated_max,

${\mathrm{\λ}}_{\max }=\frac{1}{n}\underset{i=1}{\overset{n}{\Σ}}\frac{\underset{j=1}{\overset{n}{\Σ}}{a}_{ij}{w}_j}{w_i}---(3-4)$

5. Mode of Level Simple Sequence and consistency check,

If the element of judgment matrix A meets following formula:

$a_{ij}=\frac{a_{ik}}{a_{jk}},i,j,k=1,2,...,n---(3-5)$

Matrix A is then claimed to have uniformity；

Whether matrix A is verified with uniformity：Using below equation

$CI=\frac{1}{(n-1)}({\mathrm{\λ}}_{max}-n)---(3-6)$

$CR=\frac{CI}{RI}---(3-7)$

In formula：The random Consistency Ratio of CR- judgment matrix；

N is the exponent number of judgment matrix；CI is the general coincident indicator of judgment matrix；RI is consistent for the mean random of judgment matrix Property index；Occurrence is shown in Table.

Rank matrix RI value

6. synthetic weight of each layer element to aims of systems is calculated, and is ranked up, as random Consistency Ratio CR=CI/RI ＜ When 0.10, it is believed that Mode of Level Simple Sequence result has satisfied uniformity, otherwise, the judgment matrix for being considered as preliminary foundation is can not Satisfactorily, need to readjust element scale value, till there is satisfied uniformity.

4. the method that a kind of basin chemical integrity according to claim 1-3 wherein any one is assessed, its feature exist In：Be provided with 18 test points in same detection zone, wherein have 9 points be located at river in, respectively PH, electrical conductivity, dissolved oxygen, COD, BOD, total nitrogen, total phosphorus, heavy metal in water, Volatile Phenols, in addition 9 test points are in the habitat of riverbank, respectively native The earth content of organic matter, total nitrogen content of soil, Soil total nitrogen, total potassium content of soil, soil available nitrogen content, soil available phosphorus Content, available potassium in soils content, BHC, DDT.

5. the method that a kind of basin chemical integrity according to claim 4 is assessed, it is characterised in that：The matrix R1 is The matrix that 9 row 5 is arranged, where each row are corresponding with the index detected in detection zone on river and riverbank habitat, each row generation The indicator layer of table different stage.

6. the method that a kind of basin chemical integrity according to claim 4 is assessed, it is characterised in that：In step 2, According to proportion quotiety reference table, this 18 point index weights tables are drawn, as follows

Index weights table

Obtain matrix：

A1=[0.028 0.028 0.05 0.0234 0.096 0.0234 0.0112 0.082 0.136]

A2=[0.098 0.046 0.043 0.043 0.098 0.098 0.098 0.0158 0.158].