CN115271500A - Comprehensive assessment method for remediation effect of heavy metal contaminated farmland soil - Google Patents
Comprehensive assessment method for remediation effect of heavy metal contaminated farmland soil Download PDFInfo
- Publication number
- CN115271500A CN115271500A CN202210944130.5A CN202210944130A CN115271500A CN 115271500 A CN115271500 A CN 115271500A CN 202210944130 A CN202210944130 A CN 202210944130A CN 115271500 A CN115271500 A CN 115271500A
- Authority
- CN
- China
- Prior art keywords
- index
- soil
- heavy metal
- indexes
- evaluation
- 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
- 239000002689 soil Substances 0.000 title claims abstract description 103
- 230000000694 effects Effects 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000005067 remediation Methods 0.000 title claims abstract description 63
- 229910001385 heavy metal Inorganic materials 0.000 title claims abstract description 52
- 238000011156 evaluation Methods 0.000 claims abstract description 75
- 239000011159 matrix material Substances 0.000 claims abstract description 58
- 230000008569 process Effects 0.000 claims abstract description 16
- 238000012163 sequencing technique Methods 0.000 claims abstract description 12
- 230000008439 repair process Effects 0.000 claims description 21
- 230000008901 benefit Effects 0.000 claims description 12
- 238000005516 engineering process Methods 0.000 claims description 11
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 7
- 238000010606 normalization Methods 0.000 claims description 6
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 4
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 claims description 4
- 238000005904 alkaline hydrolysis reaction Methods 0.000 claims description 4
- 244000005700 microbiome Species 0.000 claims description 4
- 229910052757 nitrogen Inorganic materials 0.000 claims description 4
- 229910052698 phosphorus Inorganic materials 0.000 claims description 4
- 239000011574 phosphorus Substances 0.000 claims description 4
- 239000011591 potassium Substances 0.000 claims description 4
- 229910052700 potassium Inorganic materials 0.000 claims description 4
- 238000005341 cation exchange Methods 0.000 claims description 3
- 239000004016 soil organic matter Substances 0.000 claims description 3
- 229910052793 cadmium Inorganic materials 0.000 description 11
- BDOSMKKIYDKNTQ-UHFFFAOYSA-N cadmium atom Chemical compound [Cd] BDOSMKKIYDKNTQ-UHFFFAOYSA-N 0.000 description 11
- 241000209140 Triticum Species 0.000 description 7
- 235000021307 Triticum Nutrition 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- QJZYHAIUNVAGQP-UHFFFAOYSA-N 3-nitrobicyclo[2.2.1]hept-5-ene-2,3-dicarboxylic acid Chemical compound C1C2C=CC1C(C(=O)O)C2(C(O)=O)[N+]([O-])=O QJZYHAIUNVAGQP-UHFFFAOYSA-N 0.000 description 2
- 235000001674 Agaricus brunnescens Nutrition 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 239000003124 biologic agent Substances 0.000 description 2
- 235000013339 cereals Nutrition 0.000 description 2
- 239000003337 fertilizer Substances 0.000 description 2
- 239000004021 humic acid Substances 0.000 description 2
- 239000005416 organic matter Substances 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- VTYYLEPIZMXCLO-UHFFFAOYSA-L Calcium carbonate Chemical compound [Ca+2].[O-]C([O-])=O VTYYLEPIZMXCLO-UHFFFAOYSA-L 0.000 description 1
- 238000003723 Smelting Methods 0.000 description 1
- ZOYQVVVFGFZXLB-UHFFFAOYSA-N [K].[Mg].[Ca].[Si] Chemical compound [K].[Mg].[Ca].[Si] ZOYQVVVFGFZXLB-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 230000002262 irrigation Effects 0.000 description 1
- 238000003973 irrigation Methods 0.000 description 1
- 229910052752 metalloid Inorganic materials 0.000 description 1
- 150000002738 metalloids Chemical class 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000000575 pesticide Substances 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- 238000003900 soil pollution Methods 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 239000000126 substance Substances 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression 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/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- General Business, Economics & Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Probability & Statistics with Applications (AREA)
- General Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Animal Husbandry (AREA)
- Primary Health Care (AREA)
- Mining & Mineral Resources (AREA)
- Agronomy & Crop Science (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
Abstract
The invention discloses a comprehensive evaluation method for the remediation effect of heavy metal contaminated farmland soil, which comprises the following steps: determining an evaluation index of the remediation effect of the heavy metal contaminated farmland soil; obtaining scales among all evaluation indexes by using an analytic hierarchy process, constructing a judgment matrix, and calculating to obtain the weight of each evaluation index according to the judgment matrix; obtaining sample data repaired by adopting a plurality of different repairing modes in the target area according to the evaluation index; and comprehensively evaluating and sequencing the repairing effects of various different repairing modes in the target area by using an approximate ideal solution sequencing method based on the sample data and the weight corresponding to each evaluation index. The method adopts an analytic hierarchy process combined with an approximate ideal solution ordering method to evaluate and order different remediation modes adopted by the heavy metal polluted farmland soil, thereby realizing comprehensive evaluation of the remediation effect of the heavy metal polluted farmland soil.
Description
Technical Field
The invention belongs to the technical field of soil remediation evaluation, and particularly relates to a comprehensive evaluation method for the remediation effect of heavy metal contaminated farmland soil.
Background
In recent years, farmland soil is affected by mining and smelting, industrial production, atmospheric sedimentation, sewage irrigation, excessive use of fertilizers and pesticides and the like for a long time, and serious soil heavy metal and metalloid pollution problems are formed. Among them, the remediation of heavy metal contaminated farmland soil gradually becomes a research hotspot in the field of soil remediation, and scholars at home and abroad make a great deal of experimental research on the heavy metal contaminated farmland remediation technology, and certain effect is achieved.
The existing heavy metal polluted farmland soil remediation technology can effectively reduce the heavy metal content of edible parts of crops to a certain extent, and achieves a certain remediation effect. The key of restoring heavy metal contaminated farmland soil is the selection of the restoration mode, and the main factors influencing the restoration effect of the heavy metal contaminated farmland soil comprise: repairing mode, soil heavy metal form and bioavailability, soil physical and chemical properties, soil type, soil texture, crop species and the like. Therefore, how to evaluate the repairing effect of different repairing modes on the heavy metal polluted farmland soil so as to screen out a suitable repairing method is an urgent problem to be solved in the field of farmland soil heavy metal pollution repairing at present.
Most of the existing evaluation methods are used for evaluating single indexes of soil pollution conditions, and mainly comprise a single-factor index method, an inner-Meiluo index method, a soil accumulation index method and the like. The single-factor index method is only suitable for evaluating a single-factor polluted area; the internal Metro index method compares the soil environment as a whole with a historical background value, and is difficult to reflect the quality change characteristics of pollution; the cumulative index method focuses on single heavy metals, and does not consider bioavailability and different pollution contribution ratios of factors.
In fact, in the evaluation of the farmland polluted soil remediation effect, various factors such as soil quality, crop growth safety, remediation economy and the like can be involved, and the existing evaluation method is difficult to realize the comprehensive evaluation of the heavy metal polluted farmland soil remediation effect.
Disclosure of Invention
In order to solve the defects of the prior art, the invention provides a comprehensive assessment method for the remediation effect of heavy metal contaminated farmland soil, which adopts an Analytic Hierarchy Process (AHP) combined approach to an ideal solution ordering process (TOPSIS) to evaluate and order different remediation modes adopted by the heavy metal contaminated farmland soil, realizes the comprehensive assessment of the remediation effect of the heavy metal contaminated farmland soil, and screens to obtain an optimal remediation mode.
In order to achieve the purpose, the disclosure provides a comprehensive assessment method for the remediation effect of heavy metal contaminated farmland soil, which comprises the following steps:
determining an evaluation index of the remediation effect of the heavy metal contaminated farmland soil;
obtaining scales among all evaluation indexes by using an analytic hierarchy process, constructing a judgment matrix, and calculating to obtain the weight of each evaluation index according to the judgment matrix;
obtaining sample data repaired by adopting a plurality of different repairing modes in the target area based on the evaluation index;
and comprehensively evaluating and sequencing the repairing effects of various different repairing modes in the target area by using an approximate ideal solution sequencing method based on the sample data and the weight corresponding to each evaluation index.
According to the further technical scheme, the evaluation indexes comprise criterion layer indexes and index layer indexes, and the criterion layer indexes comprise restoration technical effect indexes, soil quality indexes and economic benefit indexes.
According to a further technical scheme, index layer indexes of the technical effect repairing indexes comprise the content of heavy metal in an effective state in soil and the content of heavy metal in edible parts of crops;
the indexes of the index layer of the soil quality index comprise a soil pH value, a soil organic matter content, a soil alkaline hydrolysis nitrogen content, a soil available phosphorus content, a soil quick-acting potassium content, a soil cation exchange capacity and a soil microorganism diversity index;
the index layer indexes of the economic benefit indexes comprise crop yield and repair cost.
According to a further technical scheme, the method for obtaining the scale among the evaluation indexes by using the analytic hierarchy process specifically comprises the following steps: and comparing each index in each level index pairwise according to expert experience, and determining the scale of each index in the level through comparison by analyzing the importance degree of each index factor in the level relative to the index factor of the previous level.
According to a further technical scheme, the weight of each evaluation index is calculated according to the judgment matrix, and the method specifically comprises the following steps:
carrying out normalization processing on the judgment matrix, wherein the formula is as follows:
in the above formula, a ji ＝1/a ij ，
Is a feature vector; i, j =1, \ 8230, n, n is constant, representing the number of all indexes;
and calculating the weight of each evaluation index according to the judgment matrix, wherein the formula is as follows:
in the above formula, w i Is the eigenvalue of the given matrix, i.e. represents the weight of the ith index in all indices.
According to the further technical scheme, after the weight of each evaluation index is obtained through calculation, consistency check is carried out on the obtained weight, and the method specifically comprises the following steps: and calculating a consistency ratio CR based on the judgment matrix, and if the consistency ratio CR is less than 0.1, determining that the judgment matrix passes the consistency test and the weight of each evaluation index is reasonably calculated.
In a further technical solution, the calculation formula of the consistency ratio CR is:
in the above formula, w i Represents the weight of the ith index in all indexes, n is a constant and represents the number of all indexes, and lambda max In order to judge the maximum eigenvalue of the matrix, CI is a consistency index, and RI is a random one-time index.
According to the further technical scheme, based on sample data and weight corresponding to each evaluation index, comprehensive evaluation and sequencing are carried out on the repairing effects of various different repairing modes in the target area by utilizing an approximate ideal solution sequencing method, and the method specifically comprises the following steps:
carrying out normalization weighting processing on the obtained sample data to obtain a weighting normalized decision matrix;
determining an ideal solution and a negative ideal solution of the matrix according to the weighted normalized decision matrix;
calculating Euclidean distance values from each repairing mode to an ideal solution and a negative ideal solution based on sample data corresponding to each evaluation index;
calculating the relative closeness of each repairing mode based on the Euclidean distance value;
and performing comprehensive evaluation sorting according to the relative closeness.
According to the further technical scheme, normalization weighting processing is carried out on the obtained sample data to obtain a weighting normalized decision matrix, and the method specifically comprises the following steps:
constructing a multi-parameter decision matrix based on sample data, and standardizing the multi-parameter decision matrix to obtain a standardized decision matrix;
and multiplying the normalized decision matrix by a weight vector formed by the weights of the evaluation indexes to obtain a weighted normalized decision matrix.
According to a further technical scheme, based on the Euclidean distance value, the relative closeness of each repairing mode is calculated, and specifically:
in the above formula,. DELTA. i In order to be a relative degree of proximity,
for the euclidean distance values between each repair mode and the ideal solution,
for the Euclidean distance values between each repair method and the negative ideal solution, m represents the number of repair methods used.
The technical scheme has the following beneficial effects:
1. the comprehensive evaluation system for the remediation effect of the heavy metal contaminated farmland soil is provided based on an Analytic Hierarchy Process (AHP) and an approximate ideal solution ordering method (TOPSIS), the remediation effect comprehensive evaluation system for the heavy metal contaminated farmland soil is provided, the remediation technical effect, the soil quality and the economic benefit are comprehensively selected as criteria layer indexes, then the optimal technology suitable for the remediation of the local cadmium contaminated farmland is selected, and scientific guidance is provided for the formulation of the remediation decision of the heavy metal contaminated farmland.
2. The comprehensive assessment method adopts an Analytic Hierarchy Process (AHP) combined approach to an ideal solution ordering process (TOPSIS), carries out assessment ordering on different remediation modes adopted by heavy metal polluted farmland soil, realizes comprehensive assessment on the remediation effect of the heavy metal polluted farmland soil, screens to obtain an optimal remediation mode, and provides a theoretical basis for later-stage remediation treatment and ecological environment protection of the heavy metal polluted farmland soil.
3. After the cadmium-polluted farmland soil remediation effect is scientifically evaluated, the cadmium-polluted farmland soil remediation effect can be efficiently solved by formulating reasonable remediation means or researching and developing new remediation technologies.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, are included to provide a further understanding of the invention, and are incorporated in and constitute a part of this specification, illustrate exemplary embodiments of the invention and together with the description serve to explain the invention and not to limit the invention.
FIG. 1 is a flow chart of a comprehensive evaluation method according to an embodiment of the present invention;
FIG. 2 is a schematic structural diagram of an evaluation index system framework in an embodiment of the present invention.
Detailed Description
It is to be understood that the following detailed description is exemplary and is intended to provide further explanation of the invention as claimed. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the invention. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof, unless the context clearly indicates otherwise.
Aiming at the problem that the existing assessment method is difficult to comprehensively consider and assess the remediation effect of heavy metal contaminated farmland soil in multiple aspects and multiple levels, the invention provides a comprehensive assessment method for the remediation effect of heavy metal contaminated farmland soil.
The embodiment provides a comprehensive evaluation method for the remediation effect of heavy metal contaminated farmland soil, as shown in fig. 1, comprising the following steps:
determining an evaluation index of the remediation effect of the heavy metal contaminated farmland soil;
obtaining scales among all evaluation indexes by using an analytic hierarchy process, constructing a judgment matrix, and calculating to obtain the weight of each evaluation index according to the judgment matrix;
obtaining sample data repaired by adopting a plurality of different repairing modes in the target area according to the evaluation index;
and comprehensively evaluating and sequencing the repairing effects of various different repairing modes in the target area by using an approximate ideal solution sequencing method based on the sample data and the weight corresponding to each evaluation index.
In the embodiment, a certain cadmium-polluted farmland is selected as a test site, the pH value of the soil of the farmland is 5.69, the CEC is 20.63cmol/kg, the organic matter content is 13.29g/kg, the total cadmium content is 0.44mg/kg, the effective cadmium content is 0.136mg/kg, and the soil cadmium pollution risk of the farmland exists according to the standard of cadmium pollution risk screening values (0.3 mg/kg, pH is less than or equal to 6.5) of the soil of the farmland. Therefore, the test site is selected for repairing the cadmium-polluted farmland soil.
Firstly, determining the evaluation index of the remediation effect of the cadmium-polluted farmland soil. As shown in fig. 2, in this embodiment, a questionnaire is developed for domestic ecological restoration, soil restoration experts, and the like, and literature research is performed to determine criteria layer indexes and index layer indexes, where the criteria layer indexes include restoration technical effect indexes, soil quality indexes, and economic benefit indexes, and an evaluation index system frame is constructed by using these three indexes.
Further, on the basis of determining the index of the criterion layer, the index of the index layer is determined: selecting the cadmium content of the soil in an effective state and the cadmium content of the wheat grains as index layer indexes of the effect indexes of the repairing technology; selecting the pH value of soil, the content OM of soil organic matter, the content AN of soil alkaline hydrolysis nitrogen, the content AP of soil available phosphorus, the content AK of soil quick-acting potassium, the soil cation exchange quantity (namely soil CEC) and the soil microorganism diversity index as indexes of AN index layer of a soil quality index; the crop yield and the repair cost are selected as index layer indexes of economic benefit indexes, wherein the crops planted in the farmland are wheat, and therefore the wheat yield and the repair cost are selected as index layer indexes of the economic benefit indexes.
Next, the weights of the evaluation indexes are calculated. Specifically, the scale between each evaluation index is obtained by using an Analytic Hierarchy Process (AHP), namely, each index in each level index is compared pairwise according to expert experience, and the scale of each index in the level is determined by comparing and analyzing the importance degree of each index factor in the level relative to the previous level.
Taking a criterion layer index as an example, the criterion layer index comprises a restoration technical effect index, a soil quality index and an economic benefit index, firstly, the restoration technical effect index and the soil quality index are compared, as shown in the following table 1, the importance degree of the restoration technical effect factor and the soil quality factor relative to the previous level (namely, the target layer and the restoration technology screening) in the restoration effect evaluation process is compared and judged according to the expert experience, and if the importance degrees of the restoration technical effect factor and the soil quality factor are the same, the scale a of the restoration technical effect index relative to the soil quality index is recorded ij Is 1; if the effect factors of the restoration technology are more important than the soil quality factors, dividing the importance degrees into slightly important, obviously important, strongly important and extremely important, and determining the scale a of the effect indexes of the restoration technology relative to the soil quality indexes according to different importance degrees ij Is 3, 5, 7, 9; similarly, if the importance of the remediation technique effect factor compared to the soil quality factor lies between the importance of the divisions, the scale a of the remediation technique effect index relative to the soil quality index can be recorded ij Is 2, 4, 6 and 8.
TABLE 1 evaluation index Scale determination
And then, comparing and determining the scale of the restoration technical effect index relative to the economic benefit index and the scale of the soil quality index relative to the economic benefit index by the comparison scheme.
By the scheme, the scale among the evaluation indexes in each level is obtained by using an Analytic Hierarchy Process (AHP), so that a standard judgment matrix of the criterion layer indexes and the index layer indexes is constructed:
A＝{a ij },i,j＝1,…n
in the above formula, a ij >0，a ji ＝1/a ij I and j are respectively expressed as any two indexes in a certain level index, and n is a constant and expresses the number of all evaluation indexes in the level.
Normalizing the judgment matrix, wherein the formula is as follows:
in the above formula, a ji ＝1/a ij ，
Is a feature vector.
And calculating the weight of each evaluation index according to the judgment matrix:
in the above formula, w i Is the eigenvalue of a given matrix, in fact, w i Representing the weight of the ith index in all indexes.
Each evaluation is obtained through calculationWeight of index w = (w) 1 ,w 2 ,…,w n ) T And then, carrying out consistency check on the obtained weight to ensure that the obtained weight is scientific and credible. Specifically, the consistency ratio CR is calculated based on the judgment matrix, and if the consistency ratio CR is less than 0.1, the judgment matrix is considered to pass the consistency test, and the obtained evaluation indexes have reasonable weights.
The process of calculating the consistency ratio CR is: first, the maximum eigenvalue λ of the decision matrix is calculated max ：
In the above formula, (Aw) i The i-th index of the vector Aw is represented, a represents the decision matrix, and w represents the eigenvalue.
Then, calculating to obtain a consistency index based on the maximum characteristic value, and further calculating to obtain a consistency ratio:
in the above formula, CI is a consistency index, and RI is a random one-time index.
In this embodiment, the criterion layer weight results are sequentially denoted as B1, B2, and B3, the index layer weight results are sequentially denoted as C1, C2, C3, C4, C5, C6, C7, C8, C9, C10, and C11, the comprehensive criterion layer and index layer weight results are sequentially denoted as W1, W2, W3, W4, W5, W6, W7, W8, W9, W10, and W11, and the comprehensive evaluation index weight determination results of different repair modes (i.e., different repair technologies) determined by the above scheme are shown in table 2 below.
TABLE 2 weight determination results for different comprehensive evaluation indexes
And then, acquiring sample data repaired by adopting a plurality of different repairing modes in the target area according to the evaluation indexes.
In this embodiment, a certain farmland at a test site is selected as a target area, the target area is divided into 8 small areas, and eight different restoration methods are respectively adopted to perform restoration in the same time period. Wherein, the materials of eight repairing modes are selected as follows: the basic conditions of the repair materials are obtained simultaneously by using a wheat variety T1 with low cadmium accumulation, humic acid T2, a silicon-calcium-potassium-magnesium fertilizer T3, a biological agent T4, mushroom dregs T5, oyster shell powder T6, wormcast T7 and a leaf surface resistance control agent T8, as shown in the following table 3.
TABLE 3 repair materials base case
At the same time, based on the evaluation indexes, collecting sample data, and respectively obtaining actual data after each repairing mode in a target area is repaired, namely performing index analysis on a crop sample (namely a wheat sample) and a soil sample to obtain the cadmium content of wheat grains, the wheat yield, the effective cadmium content of soil, the pH value of the soil, the organic matter content, the alkaline hydrolysis nitrogen content, the effective phosphorus content, the quick-acting potassium content, the soil microorganism diversity index and the repairing cost of the repairing mode.
And comprehensively evaluating and sequencing the repairing effects of a plurality of different repairing modes in the target area by using an approximate ideal solution sequencing method TOPSIS (technique for order preference) based on the sample data corresponding to each evaluation index and the weight of each evaluation index.
Firstly, carrying out normalization weighting processing on the acquired sample data to obtain a weighting normalized decision matrix. Setting a multi-parameter decision matrix asY＝{y ij Normalizing the normalized decision matrix to obtain a normalized decision matrix Z = { Z = ij And that is:
multiplying the normalized decision matrix Z by a weight vector W composed of evaluation index weights to obtain a weighted normalized matrix X = { X = ij And i.e.:
x ij ＝w j ·z ij ,i＝1,…,m；j＝1,…,n
in the above formula, m represents the number of repair methods used, and n represents the number of evaluation indexes.
In this embodiment, the weighted normalized decision matrix shown in table 4 below is obtained through the above scheme
TABLE 4 weight normalized decision matrix for comprehensive evaluation system
Secondly, after a weighted normalized decision matrix is obtained through calculation, an ideal solution x of the normalized matrix is determined + And negative ideal solution x - . Specifically, let the ideal solution
The jth parameter of
The jth parameter of the negative ideal solution is
In the above formula, I + Is a benefit type parameter; I.C. A - Are cost-based parameters.
Then, according to the sample data of each repairing mode, calculating the Euclidean distance value from each repairing mode to the ideal solution and the negative ideal solution, wherein the formula is as follows:
in the process of coming into the market,
is the Euclidean distance value between the ideal solution and the target solution;
is the euclidean distance value from the negative ideal solution.
And finally, calculating the relative closeness of each repairing mode based on the Euclidean distance value, wherein the formula is as follows:
in the above formula,. DELTA. i Relative proximity, Δ i The order of the repair modes, delta, is arranged from large to small i The larger the size, the better the comprehensive repair effect of the repair mode.
And obtaining comprehensive evaluation ranking of the remediation effect of the heavy metal polluted farmland region according to the relative closeness of each remediation mode calculated by using the TOPSIS (approximate ideal solution ordering method) so as to screen out the optimal technology suitable for the remediation of the local cadmium polluted farmland and provide scientific guidance for the decision making of the remediation of the heavy metal polluted farmland.
In this example, the relative proximity ranking of eight different repair methods shown in table 5 below was obtained by the above scheme, and it can be seen from table 5 that the relative proximity ranking of the mushroom dregs T5 treatment is 0.752 at the first, which is the most ideal repair method under the present test conditions, and the relative proximity ranking is 0.679 and 0.668 for the biological agent T4 and the humic acid T2, respectively. The analysis result can provide a theoretical basis for the soil heavy metal remediation treatment and ecological environment protection of the test farmland area in the later period.
TABLE 5 Final evaluation ranking results for different repair modes
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Although the embodiments of the present invention have been described with reference to the accompanying drawings, it is not intended to limit the scope of the invention, and it should be understood by those skilled in the art that various modifications and variations can be made without inventive changes in the technical solutions of the present invention.
Claims (10)
1. A comprehensive assessment method for heavy metal contaminated farmland soil remediation effect is characterized by comprising the following steps:
determining an evaluation index of the remediation effect of the heavy metal contaminated farmland soil;
obtaining scales among all evaluation indexes by using an analytic hierarchy process, constructing a judgment matrix, and calculating to obtain the weight of each evaluation index according to the judgment matrix;
obtaining sample data repaired by adopting a plurality of different repairing modes in the target area based on the evaluation index;
and comprehensively evaluating and sequencing the repairing effects of various different repairing modes in the target area by using an approximate ideal solution sequencing method based on the sample data and the weight corresponding to each evaluation index.
2. The method for comprehensively evaluating the soil remediation effect of heavy metal contaminated farmland as claimed in claim 1, wherein the evaluation indexes comprise criterion layer indexes and index layer indexes, and the criterion layer indexes comprise a remediation technology effect index, a soil quality index and an economic benefit index.
3. The comprehensive assessment method for the soil remediation effect of heavy metal contaminated farmland as claimed in claim 2, wherein the index layer indexes of the remediation technical effect index comprise the soil available heavy metal content and the crop edible part heavy metal content;
the indexes of the index layer of the soil quality index comprise a soil pH value, a soil organic matter content, a soil alkaline hydrolysis nitrogen content, a soil available phosphorus content, a soil quick-acting potassium content, a soil cation exchange amount and a soil microorganism diversity index;
the index layer indexes of the economic benefit indexes comprise crop yield and repair cost.
4. The method for comprehensively evaluating the remediation effect of heavy metal contaminated farmland soil according to claim 1, wherein the obtaining of the scale between the evaluation indexes by using an analytic hierarchy process specifically comprises: and comparing every two indexes in each level index according to expert experience, analyzing the importance degree of each index factor in the level relative to the index factor of the previous level, and comparing to determine the scale of each index in the level.
5. The comprehensive assessment method for the soil remediation effect of heavy metal contaminated farmland as claimed in claim 1, wherein the weights of the evaluation indexes are calculated according to the judgment matrix, and the method specifically comprises the following steps:
carrying out normalization processing on the judgment matrix, wherein the formula is as follows:
in the above formula, a ji ＝1/a ij ，
Is a feature vector; i, j =1, \ 8230, m, n are constants, representing the number of all indexes;
and calculating the weight of each evaluation index according to the judgment matrix, wherein the formula is as follows:
in the above formula, w i Is the eigenvalue of the given matrix, i.e. represents the weight of the ith index among all the indices.
6. The comprehensive assessment method for the soil remediation effect of heavy metal contaminated farmland as claimed in claim 1, wherein after the weights of the evaluation indexes are calculated, consistency check is performed on the obtained weights, and the method specifically comprises the following steps: and calculating a consistency ratio CR based on the judgment matrix, and if the consistency ratio CR is less than 0.1, determining that the judgment matrix passes the consistency test and the weight of each evaluation index is reasonably calculated.
7. The method for comprehensively evaluating the soil remediation effect of heavy metal contaminated farmland as claimed in claim 6, wherein the consistency ratio CR is calculated by the following formula:
in the above formula, w i Represents the weight of the ith index in all indexes, n is a constant and represents the number of all indexes, and lambda max In order to judge the maximum eigenvalue of the matrix, CI is a consistency index, and RI is a random one-time index.
8. The comprehensive assessment method for the remediation effect of the heavy metal contaminated farmland soil according to claim 1, which is characterized in that based on sample data and weights corresponding to each evaluation index, an approximate ideal solution ranking method is utilized to comprehensively assess and rank the remediation effects of a plurality of different remediation modes in a target area, and specifically comprises the following steps:
carrying out normalization weighting processing on the obtained sample data to obtain a weighting normalized decision matrix;
determining an ideal solution and a negative ideal solution of the matrix according to the weighted normalized decision matrix;
calculating Euclidean distance values from each repairing mode to an ideal solution and a negative ideal solution based on sample data corresponding to each evaluation index;
calculating the relative closeness of each repairing mode based on the Euclidean distance value;
and performing comprehensive evaluation ranking according to the relative closeness.
9. The comprehensive assessment method for the soil remediation effect of heavy metal contaminated farmland according to claim 8, characterized by performing normalized weighting processing on the acquired sample data to obtain a weighted normalized decision matrix, specifically comprising:
constructing a multi-parameter decision matrix based on sample data, and standardizing the multi-parameter decision matrix to obtain a standardized decision matrix;
and multiplying the normalized decision matrix by a weight vector formed by the weights of the evaluation indexes to obtain a weighted normalized decision matrix.
10. The comprehensive assessment method for the soil remediation effect of heavy metal contaminated farmland as claimed in claim 8, wherein the relative closeness of each remediation mode is calculated based on the Euclidean distance value, specifically:
in the above formula,. DELTA. i For the purpose of the relative degree of proximity,
for the euclidean distance values between each repair mode and the ideal solution,
for the euclidean distance values between each repair method and the negative ideal solution, m represents the number of repair methods employed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210944130.5A CN115271500A (en) | 2022-08-05 | 2022-08-05 | Comprehensive assessment method for remediation effect of heavy metal contaminated farmland soil |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210944130.5A CN115271500A (en) | 2022-08-05 | 2022-08-05 | Comprehensive assessment method for remediation effect of heavy metal contaminated farmland soil |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115271500A true CN115271500A (en) | 2022-11-01 |
Family
ID=83749675
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210944130.5A Pending CN115271500A (en) | 2022-08-05 | 2022-08-05 | Comprehensive assessment method for remediation effect of heavy metal contaminated farmland soil |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115271500A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115825392A (en) * | 2022-12-13 | 2023-03-21 | 云南大学 | Mining area heavy metal contaminated soil ecological restoration technology evaluation method and system |
| CN116011745A (en) * | 2022-12-20 | 2023-04-25 | 速度时空信息科技股份有限公司 | Ecological restoration scheme optimization method and system based on cloud service |
-
2022
- 2022-08-05 CN CN202210944130.5A patent/CN115271500A/en active Pending
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115825392A (en) * | 2022-12-13 | 2023-03-21 | 云南大学 | Mining area heavy metal contaminated soil ecological restoration technology evaluation method and system |
| CN115825392B (en) * | 2022-12-13 | 2024-03-26 | 云南大学 | Evaluation method and system for ecological restoration technology of mining area heavy metal contaminated soil |
| CN116011745A (en) * | 2022-12-20 | 2023-04-25 | 速度时空信息科技股份有限公司 | Ecological restoration scheme optimization method and system based on cloud service |
| CN116011745B (en) * | 2022-12-20 | 2024-02-13 | 速度科技股份有限公司 | Ecological restoration scheme optimization method and system based on cloud service |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115271500A (en) | Comprehensive assessment method for remediation effect of heavy metal contaminated farmland soil | |
| Dengiz et al. | Multi-criteria assessment and geostatistical approach for determination of rice growing suitability sites in Gokirmak catchment | |
| Mousavi-Avval et al. | Combined application of Life Cycle Assessment and Adaptive Neuro-Fuzzy Inference System for modeling energy and environmental emissions of oilseed production | |
| CN109284921B (en) | Agricultural irrigation water quality dynamic evaluation method based on hidden Markov model | |
| CN112893427B (en) | Intelligent decision-making method for heavy metal polluted farmland restoration treatment | |
| Montazar et al. | An analytical hierarchy model for assessing global water productivity of irrigation networks in Iran | |
| CN108520345A (en) | Evaluation for cultivated-land method and system based on GA-BP neural network models | |
| CN113785676A (en) | Saline-alkali soil improvement system and construction method thereof | |
| CN109064039B (en) | Farmland soil health evaluation method | |
| CN109858180B (en) | Regional crop water consumption space pattern optimization design method | |
| CN110991921A (en) | Three-dimensional magic cube-based farmland ecological quality comprehensive evaluation method | |
| CN110942242A (en) | Method for evaluating operation performance of incubator of scientific and technological enterprise | |
| Song et al. | Spatial variability of soil properties in red soil and its implications for site-specific fertilizer management | |
| CN109604316B (en) | Screening method for radioactive contaminated soil remediation technology | |
| CN110782131A (en) | Performance evaluation method for water cooperation organization of farmers | |
| CN113822770A (en) | Farmland construction area intelligent demarcation method and system | |
| CN111428985B (en) | Assessment method for regional water ecological bearing capacity | |
| CN113344409A (en) | Evaluation method and system for facility continuous cropping soil quality | |
| Chang et al. | Soil health assessment methods: Progress, applications and comparison | |
| CN111462834A (en) | Method and system for predicting probability value of excessive cadmium in plant | |
| CN104778368A (en) | Pareto set individual ranking method aiming at high-dimensional multi-objective optimization problem | |
| CN115775042B (en) | Crop heavy metal enrichment risk prediction method and system based on Bayesian theory | |
| Mandal et al. | Soil health measurement techniques | |
| Dehaghia et al. | Modified TODIM method for water reuse application: a case study from Iran | |
| Zangeneh et al. | Assessment of agricultural mechanization status of potato production by means of artificial neural network model |
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 |