CN111652487A - Ecological influence evaluation system for road construction project in natural conservation area - Google Patents
Ecological influence evaluation system for road construction project in natural conservation area Download PDFInfo
- Publication number
- CN111652487A CN111652487A CN202010455025.6A CN202010455025A CN111652487A CN 111652487 A CN111652487 A CN 111652487A CN 202010455025 A CN202010455025 A CN 202010455025A CN 111652487 A CN111652487 A CN 111652487A
- Authority
- CN
- China
- Prior art keywords
- influence
- road
- grade
- evaluation index
- area
- 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.)
- Granted
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 350
- 238000010276 construction Methods 0.000 title claims abstract description 170
- 238000000034 method Methods 0.000 claims abstract description 40
- 241000894007 species Species 0.000 claims description 141
- 239000011159 matrix material Substances 0.000 claims description 136
- 238000004364 calculation method Methods 0.000 claims description 66
- 238000011835 investigation Methods 0.000 claims description 14
- 230000035945 sensitivity Effects 0.000 claims description 14
- 241001465754 Metazoa Species 0.000 claims description 13
- 230000002411 adverse Effects 0.000 claims description 13
- 239000002689 soil Substances 0.000 claims description 13
- 238000004162 soil erosion Methods 0.000 claims description 13
- 238000012360 testing method Methods 0.000 claims description 13
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 9
- 230000007613 environmental effect Effects 0.000 claims description 8
- 238000013467 fragmentation Methods 0.000 claims description 8
- 238000006062 fragmentation reaction Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 7
- 238000007726 management method Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 230000001681 protective effect Effects 0.000 claims description 4
- 239000000470 constituent Substances 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000005065 mining Methods 0.000 claims description 3
- 230000000644 propagated effect Effects 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 abstract description 3
- 238000007689 inspection Methods 0.000 abstract 1
- 241000196324 Embryophyta Species 0.000 description 20
- 241000722818 Aralia Species 0.000 description 4
- 235000014722 Aralia cordata Nutrition 0.000 description 4
- 235000004446 Aralia racemosa Nutrition 0.000 description 4
- 235000010658 Lavandula latifolia Nutrition 0.000 description 4
- 241000157862 Dicamptodontidae Species 0.000 description 3
- 238000013316 zoning Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000000903 blocking effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004321 preservation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 101100400452 Caenorhabditis elegans map-2 gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000000007 visual effect Effects 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 relates to the field of natural protection, in particular to a system for evaluating ecological influence of road construction projects in a natural protection area. The method comprises the steps of determining evaluation indexes, judging related to endangered species, determining influence levels of the influence evaluation indexes, determining membership functions of the influence evaluation indexes, calculating weight of the influence evaluation indexes, carrying out fuzzy comprehensive evaluation and judging ecological influence levels; the invention determines the membership degree of each element by using the 'semi-reduced trapezoidal distribution diagram', makes the influence of each index on the natural protection area more specific and definite, and carries out quick inspection by calculating the ecological influence index EI of the road construction project. The invention integrates the early warning index and the one-vote rejection index, can quickly evaluate the overall ecological influence of the road construction project, and can objectively and quantitatively reflect the ecological influence degree of the road construction project on the natural protection area.
Description
Technical Field
The invention relates to the field of natural protection, in particular to a system for evaluating ecological influence of road construction projects in a natural protection area.
Background
Road construction brings great benefits to the human society and also has direct or indirect influence on many ecological processes. A large amount of theoretical researches are developed aiming at the road ecology internationally, and the method is widely applied to the evaluation of the influence of the road ecology, and has important reference and reference values for developing the research of the influence of the road ecology in China. However, the study of the ecological influence of road construction on the natural protection area usually focuses on a specific aspect, such as the aspects of road influence domain, animal death, evaluation of influence of plants and soil in the road domain, biodiversity and the like, and the study of the size of the overall ecological influence of the natural protection area to judge whether the influence can be accepted is rarely reported.
In 2011, 4 months, the former environmental protection department issued ecological impact of environmental impact evaluation technical guide (HJ/19-2011). The guiding rule standardizes and systematizes the analysis content of the ecological influence of the engineering, analyzes and predicts the ecological influence possibly generated by the engineering construction through engineering analysis and ecological status investigation and evaluation, but most of the ecological influence analysis is qualitative descriptive prediction. In the same year of 10 months, the original environmental protection department issues 'Special ecological influence report formulation guide related to national level natural protection area construction projects', the guide can be applied to evaluation of ecological influence of road construction projects on natural protection areas, but the guide does not specially aim at road construction projects in the natural protection areas, the guide does not fully consider the particularity of the road construction projects, and evaluation results often cannot truly reflect the actual ecological influence of the road projects. The patent with the patent number of CN110245833A and the publication date of 2019, 9 and 17 discloses a method and a system for evaluating ecological influences of development and construction projects in a natural protected area, a remote sensing image of the natural protected area is obtained, remote sensing identification is carried out on one construction project based on the remote sensing image, each evaluation index is calculated respectively aiming at two different periods, specific evaluation indexes of each evaluation index in the two different periods are obtained, and a specific index change value is calculated according to the two specific indexes; substituting the specific index change value into an ecological influence evaluation model to obtain an influence index, and judging the influence degree of the influence index; the invention also discloses a system for evaluating the ecological influence of the development and construction project of the natural conservation area, which comprises a remote sensing identification subsystem, a project extraction subsystem, an ecological system conservation area influence evaluation subsystem and a wild animal conservation area influence evaluation subsystem. However, the distribution condition of endangered species in a natural protection area is ignored, indirect influences (water and soil loss, foreign invasive species, change of the habitat of a protected object and the like) caused by road construction cannot be fully considered by the evaluation indexes, and deviation of the evaluation result is possibly caused.
Therefore, it is an urgent need to solve the problem for those skilled in the art to design an evaluation system capable of objectively and quantitatively reflecting the ecological influence of road construction projects in natural reserve.
Disclosure of Invention
1. Problems to be solved
Aiming at the defect that the evaluation result of a road construction project ecological influence evaluation system in a natural preservation area is inaccurate in the prior art, the invention aims to provide the road construction project ecological influence evaluation system in the natural preservation area, which comprises the steps of determining evaluation indexes, judging related to endangered species, determining influence levels of the influence evaluation indexes, determining membership functions of the influence evaluation indexes, calculating weight of the influence evaluation indexes, carrying out fuzzy comprehensive evaluation and judging ecological influence levels; the invention utilizes the 'semi-trapezoid distribution map' to determine the membership degree of each element, so that the influence of each index on the natural protection area is more specific and definite, and the ecological influence degree of the road construction project on the natural protection area can be objectively and quantitatively reflected.
2. Technical scheme
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
the invention provides a natural reserve road construction project ecological impact assessment system, which comprises the following steps:
s100, determining an evaluation index; determining ecological influence evaluation indexes of road construction projects on natural protection areas; the ecological influence evaluation index comprises a road attribute influence evaluation index A1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4;
Road engineering influence evaluation index A1Including road lengths involving nature-preserving zonesDegree A11Area A of natural protection area invaded by road and accessory facilities12And road grade A13Wherein the road length A11Area A of natural protection area invaded by road and accessory facilities12Are all quantitative indexes, road grade A13Is a qualitative index; susceptibility influence evaluation index A2Functional area A including road construction project traversal21And to distribution of endangered species A22Wherein the road construction project passes through the functional area A21And to distribution of endangered species A22Are all qualitative indexes; direct influence evaluation index A3Including the direct loss of a primary protected object31Area of soil erosion32And disruption of habitat influence intensity A33Wherein the primary protected object directly loses A31Area of soil erosion32And disruption of habitat influence intensity A33Are all quantitative indicators; indirect influence evaluation index A4Including foreign invasive species status A41And indirect loss of protected objects A42Wherein the foreign invasive species status A41For qualitative index, indirect loss of protected object A42Is a quantitative index; when the quantitative indexes are determined, determining the numerical values of the quantitative indexes respectively;
wherein, the road grade A13Determining according to the traffic volume, the road use task and the property, wherein the higher the road grade is, the greater the influence degree is; concerning distribution of endangered species A22According to the endangered species diversity field investigation and determination, the distribution is from scattered to concentrated, and the influence degree is stronger; foreign invasive species condition A41Data consulting, accessing and field investigation and determination are adopted;
s200, judging related endangered species; according to actual investigation, related distribution condition A of endangered species22If IUCN is related along the road and the red famous and extremely dangerous species, the red famous and extremely dangerous species and the red famous and extremely dangerous species in China are intensively distributed, the road construction project is not acceptable, and the grade evaluation is finished; if IUCN is not involved along the road and the red directory of the Chinese endangered species is extremely dangerous, endangered and national grade I key protection species centralized distribution areaIf so, entering S300;
s300, determining the influence grade of the influence evaluation index; a is prepared by11、A12、A13、A21、A22、A31、A32、A33、A41And A42Dividing the ecological environment into three levels according to the influence degree on the ecology, wherein the three levels are level I, level II and level III respectively; wherein, level I shows that the ecological influence of the road project on the natural protection area is within an acceptable range, the influence level is green and the number of the influence level is 1, level II shows that the road project has obvious influence on the natural protection area, corresponding protection measures are needed to reduce the adverse influence, the influence level is yellow and the number of the influence level is 2, level III shows that the road project has serious adverse influence on the natural protection area, the influence level is red and the number of the influence level is 3;
s400, determining a membership function of the influence evaluation index; determining the membership degree of each ecological influence evaluation index, wherein,
membership matrix R1For road attribute influence evaluation index A1Is 3 × 3 matrix, and the membership matrixMembership matrix R1The a-th row and the j-th column of (A) are denoted as raj,rajFor road attribute influence evaluation index A1The membership degree of the middle alpha influence evaluation index to the grade with the number j; wherein a is a positive integer and a is 1, 2, 3, j is a positive integer and j is 1, 2, 3;
membership matrix R2Evaluation index A for susceptibility influence2Is a 2 × 3 matrix, a membership matrixMembership matrix R2The b-th row and j-th column element of (1) is denoted as rbj,rbjEvaluation index A for susceptibility influence2The membership degree of the middle b-th influence evaluation index to the grade with the number j; wherein b is a positive integer and b is 1,2. 3, j is a positive integer and j is 1, 2, 3;
membership matrix R3To directly influence the evaluation index A3Is 3 × 3 matrix, and the membership matrixMembership matrix R3The element of row c and column j is denoted as rcj,rcjTo directly influence the evaluation index A3The membership degree of the c-th influence evaluation index to the grade with the number j; wherein c is a positive integer and c is 1, 2, 3, j is a positive integer and j is 1, 2, 3;
membership matrix R4As an indirect influence evaluation index A4Is a 2 × 3 matrix, a membership matrixMembership matrix R4The d-th row and j-th column of (1) are denoted as rdj,rdjAs an indirect influence evaluation index A4The membership degree of the d-th influence evaluation index to the grade with the number j; wherein d is a positive integer and d is 1, 2, 3, j is a positive integer and j is 1, 2, 3;
when determining the membership degrees of the 6 quantitative indexes in the S100, determining a membership degree function by adopting trapezoidal distribution, and calculating the membership degree of the grade numbered j of each quantitative index by using the determined membership degree function;
for road grade A13Is determined according to the road grade A determined in S20013Is determined by the grade of R1e1 and R1e'0; wherein e is a positive integer and e is the road grade A determined in S20013The grade number of (2); e ' is a positive integer, e ' ═ 1, 2, 3, and e ' ≠ e;
functional area A for road construction project crossing21Is determined according to the functional area A traversed by the road construction project determined in S20021Is determined by the grade of R1f1 and R1f'0; wherein f is a positive integer and f is determined in S200Functional area A crossed by fixed road construction project21The grade number of (2); f ' is a positive integer, f ' is 1, 2, 3, and f ' ≠ f;
for distribution conditions A involving endangered species22Is determined according to the distribution condition A of the related endangered species determined in S20022Is determined by the grade of R1g1 and R1g'0; wherein g is a positive integer and g is the distribution A related to endangered species determined in S20022The grade number of (2); g 'is a positive integer, g' ≠ g, 1, 2, 3;
for foreign invasive species status A41Is determined according to the determined foreign invasive species condition A in S20041Is determined by the grade of R1h1 and R1h'0; wherein h is a positive integer and h is the alien invasive species status A determined in S20041The grade number of (2); h 'is a positive integer, h' ≠ h, 1, 2, 3;
s500, calculating the weight of the influence evaluation index; evaluation index A for influence on road attribute1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4The weight of the evaluation index is calculated, and the method for influencing the weight calculation of the evaluation index is an analytic hierarchy process; and the road attribute influence evaluation index A1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4The weights of (A) are respectively denoted as W1、W2、W3And W4To obtain a weight vector T, T ═ W1,W2,W3,W4);
Evaluation index A for influence on road attribute1Middle 3 ecological impact evaluation indexes A11、A12And A13Are calculated respectively, and A is calculated11、A12And A13The weights of (A) are respectively denoted as W11、W12And W13To obtain A1Weight vector T1,T1=(W11,W12,W13);
Evaluation index A for influence on sensitivity2Middle 2 ecological impact evaluation indexes A21And A22Are calculated respectively, and A is calculated21And A22The weights of (A) are respectively denoted as W21And W22To obtain A2Weight vector T2,T2=(W21,W22);
For direct influence evaluation index A3Middle 3 ecological impact evaluation indexes A31、A32And A33Are calculated respectively, and A is calculated31、A32And A33The weights of (A) are respectively denoted as W31、W32And W33To obtain A3Weight vector T3,T3=(W31,W32,W33);
Evaluation index A for indirect influence4Middle 2 ecological impact evaluation indexes A41And A42Are calculated respectively, and A is calculated41And A42The weights of (A) are respectively denoted as W42And W42To obtain A4Weight vector T4,T4=(W42,W42);
S600, fuzzy comprehensive evaluation; according to the membership degree matrix R determined in S4001、R2、R3And R4And the weight vector T determined in S5001、T2、T3And T4Calculating road attribute influence evaluation index A1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4Membership degree vector of (1), road attribute influence evaluation index A1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4Are respectively B1、B2、B3And B4(ii) a Wherein the road attribute influences the evaluation index A1The calculation formula of the membership degree vector is B1=T1R1(ii) a Susceptibility influence evaluation index A2The calculation formula of the membership degree vector is B2=T2R2(ii) a Direct influence evaluation index A3The calculation formula of the membership degree vector is B3=T3R3(ii) a Indirect influence evaluation index A4The calculation formula of the membership degree vector is B4=T4R4(ii) a According to B1、B2、B3And B4Obtaining an influence evaluation membership degree matrix R,calculating a fuzzy evaluation matrix B according to the formula B ═ TR, (B ═ TR)1,b2,b3);
S700, judging the ecological influence level; evaluating the influence degree of the road construction project on the natural protection area by using the fuzzy evaluation matrix B; b1、b2、b3Respectively corresponding to the evaluated specific gravities of the grade I, the grade II and the grade III; level I indicates that the ecological influence of the road project on the natural conservation area is within an acceptable range; level II shows that the road project has obvious influence on the natural protected area, and protective measures need to be taken to relieve the adverse influence; grade iii indicates that the road project has a serious adverse effect on the natural reserve, which is unacceptable.
Preferably, the calculation of the influence evaluation index weight in S500 is an influence evaluation index weight comprehensive calculation, and the steps are as follows,
s501, subjective weight calculation of the impact evaluation indexes; evaluation index A for influence on road attribute1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4The subjective weights of (a) are calculated to obtain a subjective weight vector C, C ═ μ1,μ2,μ3) (ii) a The constructed judgment matrix is a first judgment matrix D1First decision matrix D1Is a 4-order matrix;
evaluation index A for influence on road attribute1Middle 3 ecological impact evaluation indexes A11、A12And A13Respectively calculating the subjective weights of A and B11、A12And A13The subjective weights of (A) are respectively denoted as mu11、μ12And mu13To obtain A1Subjective weight vector C1,C1=(μ11,μ12,μ13) (ii) a The constructed judgment matrix is a second judgment matrix D2Second decision matrix D2Is a 3-order matrix;
evaluation index A for influence on sensitivity2Middle 2 ecological impact evaluation indexes A21And A22Respectively calculating the subjective weights of A and B21And A22The subjective weights of (A) are respectively denoted as mu22And mu22To obtain A2Subjective weight vector C2,C2=(μ22,μ22) (ii) a The constructed judgment matrix is a third judgment matrix D3The third judgment matrix D3Is a 2-order matrix;
for direct influence evaluation index A3Middle 3 ecological impact evaluation indexes A31、A32And A33Respectively calculating the subjective weights of A and B31、A32And A33The subjective weights of (A) are respectively denoted as mu31、μ32And mu33To obtain A3Subjective weight vector C3,C3=(μ31,μ32,μ33) (ii) a The constructed judgment matrix is a fourth judgment matrix D4Fourth decision matrix D4Is a 3-order matrix;
evaluation index A for indirect influence4Middle 2 ecological impact evaluation indexes A41And A42Respectively calculating the subjective weights of A and B41And A42The subjective weights of (A) are respectively denoted as mu42And mu42To obtain A4Subjective weight vector C4,C4=(μ42,μ42) (ii) a The constructed judgment matrix is a fifth judgment matrix D5Fifth decision matrix D5Is a 2-order matrix;
s502, calculating objective weight of the influence evaluation index; using the first determination matrix D1Calculating to obtain road attribute influence evaluation index by adopting entropy weight methodMark A1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4Objective weight λ of1、λ2、λ3And λ4;
Using the second decision matrix D2Evaluation index A for influence of entropy weight method on road attribute1Middle 3 ecological impact evaluation indexes A11、A12And A13Respectively calculating the objective weights of A and B11、A12And A13The objective weights of (A) are respectively denoted as λ11、λ12And λ13;
Using the third determination matrix D3Evaluation index A for sensitivity influence by entropy weight method2Middle 2 ecological impact evaluation indexes A21And A22Respectively calculating the objective weights of A and B21And A22The objective weights of (A) are respectively denoted as λ22And λ22;
Using the fourth determination matrix D4Evaluating index A by direct influence of entropy weight method3Middle 3 ecological impact evaluation indexes A31、A32And A33Respectively calculating the objective weights of A and B31、A32And A33The objective weights of (A) are respectively denoted as λ31、λ32And λ33;
Using the fifth decision matrix D5Evaluation index A for indirect influence by entropy weight method4Middle 2 ecological impact evaluation indexes A41And A42Respectively calculating the objective weights of A and B41And A42The objective weights of (A) are respectively denoted as λ42And λ42;
S503, calculating the combined weight of the influence evaluation indexes; calculation of A1、A2、A3And A4Combining weight W of1、W2、W3And W4To obtain a weight vector T, T ═ W1,W2,W3,W4);
The calculation formula is as follows,
Wm=α1λm+(1-α1)μm
wherein,
Wmis AmM is a positive integer and m is more than or equal to 1 and less than or equal to 4;
μmis AmM is a positive integer and is more than or equal to 1 and less than or equal to 4;
λmis AmM is a positive integer and is more than or equal to 1 and less than or equal to 4;
Calculating road attribute influence evaluation index A1Combining weight W of11、W12And W13To obtain A1Weight vector T1,T1=(W11,W12,W13);
The calculation formula is as follows,
W1m=α2λ1m+(1-α2)μ1m
wherein,
W1mis A1mM is a positive integer and m is more than or equal to 1 and less than or equal to 3;
μ1mis A1mM is a positive integer and is more than or equal to 1 and less than or equal to 3;
λ1mis A1mM is a positive integer and is more than or equal to 1 and less than or equal to 3;
Calculation of sensitivity influence evaluation index A2Combining weight W of21And W22To obtain A2Weight vector T2,T2=(W21,W22);
The calculation formula is as follows,
W2m=α3λ2m+(1-α3)μ2m
wherein,
W2mis A2mM is a positive integer and m is more than or equal to 1 and less than or equal to 2;
μ2mis A2mM is a positive integer and is more than or equal to 1 and less than or equal to 2;
λ2mis A2mM is a positive integer and is more than or equal to 1 and less than or equal to 2;
α3is A2And α3=W21+W22-1.5;
Calculating direct influence evaluation index A3Combining weight W of31、W32And W33To obtain A3Weight vector T3,T3=(W31,W32,W33);
The calculation formula is as follows,
W3m=α2λ3m+(1-α3)μ3m
wherein,
W3mis A3mM is a positive integer and m is more than or equal to 1 and less than or equal to 3;
μ3mis A3mM is a positive integer and is more than or equal to 1 and less than or equal to 3;
λ3mis A3mM is a positive integer and is more than or equal to 1 and less than or equal to 3;
Calculating an indirect influence evaluation index A4Combining weight W of42And W42To obtain A4Weight vector T4,T4=(W42,W42);
The calculation formula is as follows,
W4m=α2λ4m+(1-α3)μ4m
wherein,
W4mis A4mM is a positive integer and m is more than or equal to 1 and less than or equal to 2;
μ4mis A4mM is a positive integer and is more than or equal to 1 and less than or equal to 2;
λ4mis A4mM is a positive integer and is more than or equal to 1 and less than or equal to 2;
α5is A4And α5=W41+W42-1.5。
Preferably, a validity check is also included in S600; the validity test is to test the fuzzy comprehensive evaluation, the calculation formula is,
wherein alpha is the maximum membership principle validity index;
β is b1,b2,b3Maximum value of (1);
gamma is b1,b2,b3The second largest value of;
when alpha is less than 0.5, evaluating the influence degree of the road construction project on the natural protection area by adopting a weighted average principle; otherwise, the influence degree of the road construction project on the natural protection area is evaluated by adopting the maximum membership principle.
Preferably, the road length A related to the nature reserve11Is divided into standard functions ofd1The formula for calculating (a) is as follows,
wherein d is1The influence strength of the road length in the natural protection area is obtained;
l is the length of the road, km, related to the natural reserve;
p is the perimeter of the natural protected area, km;
area A of natural protection area invaded by auxiliary facilities12Is divided into standard functions ofd2The formula for calculating (a) is as follows,
wherein d is2The influence strength of the area of the road and the accessory facilities thereof;
s is the area of the road and its subsidiary facilities invading the natural protection area hm2;
A is the total area of the natural reserve, hm2;
Road grade A13The division standard of (1) is that when the road construction project is a highway and a first-level road, the influence level is level III; the road construction projects are second-level and third-level roads, and the influence level is level II; the road construction project is a four-level road and other roads, and the influence level is level I;
functional area A traversed by road construction project21The division standard of (1) is that when the road construction project relates to a core area and a buffer area, the influence level is level III; when the road construction project relates to a test area, the influence level is level II; when the road construction project relates to other areas, the influence level is level I;
concerning distribution of endangered species A22The division standard is that when the road construction project relates to IUCN, the red directory of Chinese endangered species is extremely dangerous, the national class I species is scattered distribution area or the national class II centralized distribution area, the influence level is level III; when the road construction project relates to IUCN and red name records of endangered species in China, and the national level II is scattered, the influence grade is level II; when the road construction project does not relate to the distribution area of the species, the influence level is level I;
direct loss of primary protected object A31Is divided into standard functions ofWherein d is3The formula for calculating (a) is as follows,
wherein d is3Directly impact strength for the primary protected object;
when the main protected object is a plant, M is the direct loss amount of the plant individual number, the plant cluster number or the plant distribution area caused by the road construction project; when the main protected object is an animal, M is the habitat area directly invading the main protected object species;
when the main protected object is a plant, T is the plant individual number, the plant clump number or the total plant distribution area of the main protected object in the natural protected area; when the main protected object is an animal, T is the total habitat area of the main protected object in the natural protected area;
area of soil erosion32Is divided into standard functions ofWherein d is4Setting the percentage of the water and soil loss area caused by road construction in the range of 30m at both sides of the road in the natural protection area to the total area of 30m buffer areas at both sides of the road;
intensity of influence of disruption of habitat A33Is divided into standard functions ofd5The formula for calculating (a) is as follows,
wherein d is5For the environment fragmentation change rate, t is after road construction, and t-1 is before road construction;
FI is an environmental disruption index;
Skis the area of the kth habitat patch, where k is a positive integer, hm2;
A is the total area of the natural reserve, hm2;
Foreign invasive species condition A41The division standard of (1) is that when external invasive species are found along the road construction project and are propagated in a large area, the influence grade is grade III; the road construction project is sporadically distributed with external invasive species along the line, but the influence level is level II when large-area diffusion is not achieved at present; when no foreign invasive species are found along the road construction project, the influence level is level I;
indirect loss of protected objects a42The division criterion function of (a) is,d6the formula for calculating (a) is as follows,
wherein d is6Indirectly influencing the intensity for the main protected object;
r is the rate of change of the importance (or frequency of occurrence) of the main protected species in the community.
Preferably, road class A13The standard for the classification of the road grade in (1) is,
the highway is a highway with particularly important political and economic significance, is provided with four or more lanes, is provided with a central separation zone, is completely intersected and has perfect traffic safety facilities, management facilities and service facilities, and is a special highway specially used for high-speed driving of automobiles, wherein the ADT is more than or equal to 25000;
the first-level highway is a highway connecting the important political and economic culture center and part of overpasses, and ADT is more than or equal to 10000 and less than 25000;
the secondary highway is a main highway connected with politics and economic centers or a great industrial and mining area or a suburban highway with busy transportation, and ADT is more than or equal to 2000 and less than 10000;
the third level road is a branch road for communicating with a city in or above a county, and ADT is more than or equal to 200 and less than 2000;
the four-level road is a branch road for communicating counties, towns and villages, and ADT < 200.
Preferably, in S400, the road attribute influences the evaluation index a1The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
area A of natural protection area invaded by auxiliary facilities12The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
direct loss of primary protected object A31The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
area of soil erosion32Formula for calculating membership of grade I, grade II and grade IIIIn order to realize the purpose,
intensity of influence of disruption of habitat A33The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
indirect loss of protected objects a42The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
preferably, the protection object indirectly loses A42The protected object in (1) refers to a specific protected object species in a natural protected area, a rare or endangered species or a main constituent species in a representative ecosystem.
Preferably, the method also comprises an ecological influence index test, wherein the calculation formula is,
the EI is an ecological influence index of a road construction project;
Diis AiI is an integer and i is greater than or equal to 1 and less than or equal to 4;
wiis AiI is an integer and is more than or equal to 1 and less than or equal to 4;
when the EI is less than or equal to 50, the ecological influence level of the road construction project is green; when EI is more than 50 and less than 80, the ecological influence level of the road construction project is yellow, and when EI is more than or equal to 80, the ecological influence level of the road construction project is red.
Preferably, DiIncluding a road attribute index D1Sensitivity index D2Direct influence index D3And an indirect influence index D4;
Road attribute index D1The formula for calculating (a) is as follows,
wherein D is1Is a road attribute index;
d1the influence strength of the road length in the natural protection area is obtained;
d2the influence strength of the area of the road and the accessory facilities thereof;
x is the road grade influence strength; grade A of the current road13When the classified influence grade is grade III, x is 100; grade A of the current road13When the classified influence grade is grade II, x is 61.8; grade A of the current road13When the classified influence grade is grade I, x is 38.2;
sensitivity index D2The formula for calculating (a) is as follows,
wherein D is2Is a sensitivity index;
y is the influence strength of the related functional regions; functional area A crossed by road construction project21When the classification level of (1) is level III, y is 100; functional area A crossed by road construction project21When the grade of (1) is grade II, y is 38.2; functional area A crossed by road construction project21When the classification level of (1) is level I, y is 0;
z is the impact intensity related to endangered species; when concerning distribution of endangered species A22When the classification level of (1) is level III, z is 100; when concerning distribution of endangered species A22When the grade is II, z is 61.8; when concerning distribution of endangered species A22When the grade is I, z is 38.2;
direct influence index D3The formula for calculating (a) is as follows,
wherein D is3Is a direct impact index;
d3directly impact strength for the primary protected object;
d4setting the percentage of the water and soil loss area caused by road construction in the range of 30m at both sides of the road in the natural protection area to the total area of 30m buffer areas at both sides of the road;
d5to an environmental disruption index;
index of indirect influence D4The formula for calculating (a) is as follows,
wherein D is4Is an indirect influence index;
n is the foreign invasive species impact strength; when foreign invaded species status A41When the grade of (1) is grade III, n is 100; when foreign invaded species status A41When the grade of (1) is grade II, n is 61.8; when foreign invaded species status A41When the classification level of (1) is level I, n is 0;
d6the intensity is indirectly influenced for the main protected object.
Preferably, wiDerived from analytic hierarchy process and expert consulting method, w1=0.2,w2=0.3,w3=0.25,w4=0.25。
3. Advantageous effects
Compared with the prior art, the invention has the beneficial effects that:
(1) the invention discloses an ecological influence evaluation system for road construction projects in a natural reserve area, which comprises the steps of determining evaluation indexes, judging related to endangered species, determining influence levels of the influence evaluation indexes, determining membership functions of the influence evaluation indexes, calculating the weight of the influence evaluation indexes, carrying out fuzzy comprehensive evaluation and judging the ecological influence levels; the invention utilizes the 'semi-trapezoid distribution map' to determine the membership of each element, so that the influence of each index on the natural protection area is more specific and definite, and the ecological influence degree of the road construction project on the natural protection area can be objectively and quantitatively reflected.
(2) According to the ecological influence evaluation system for the road construction project in the natural conservation area, the difference coefficient between subjective and objective weights is introduced, and the combined weight method of an analytic hierarchy process and an entropy weight method is adopted to weight each influence evaluation index, so that the advantages of subjective weighting and objective weighting can be fully exerted; the method is combined with a fuzzy comprehensive evaluation method to evaluate the ecological influence of road construction projects in the natural conservation area, the evaluation time is short, the evaluation result is objective and accurate, and the application prospect is wide.
(3) According to the system for evaluating the ecological influence of the road construction project in the natural protection area, disclosed by the invention, the ecological influence index EI of the road construction project is calculated, the ecological influence of the road construction project on the natural protection area can be quantitatively reflected, the evaluation result is quickly verified, the objectivity of the evaluation result is improved, and the system is favorable for a natural protection area manager to make a decision according to the judged grade.
Drawings
FIG. 1 is a flow chart of the evaluation of the present invention;
FIG. 2 is a diagram illustrating a location relationship between a new road in a natural reserve and a functional area;
FIG. 3 is a diagram illustrating water and soil loss caused by new roads along the natural conservation area;
FIG. 4 is a graph showing a comparison of the fragmentation of a natural reserve environment before and after road construction.
Detailed Description
The invention is further described with reference to specific examples.
The invention discloses a natural reserve road construction project ecological impact assessment system, as shown in figure 1, comprising the following steps:
s100, determining an evaluation index; determining ecological influence evaluation indexes of road construction projects on natural protection areas; the ecological influence evaluation index comprises a road attribute influence evaluation index A1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4;
Road engineering influence evaluation index A1Including road length A relating to natural reserve11Area A of natural protection area invaded by road and accessory facilities12And road grade A13Wherein the road length A11Area A of natural protection area invaded by road and accessory facilities12Are all quantitative indexes, road grade A13Is a qualitative index; susceptibility influence evaluation index A2Functional area A including road construction project traversal21And to distribution of endangered species A22Wherein the road construction project passes through the functional area A21And to distribution of endangered species A22Are all qualitative indexes; direct influence evaluation index A3Including the direct loss of a primary protected object31Area of soil erosion32And disruption of habitat influence intensity A33Wherein the primary protected object directly loses A31Area of soil erosion32And disruption of habitat influence intensity A33Are all quantitative indicators; indirect influence evaluation index A4Including foreign invasive species status A41And indirect loss of protected objects A42Wherein the foreign invasive species status A41For qualitative index, indirect loss of protected object A42Is a quantitative index; when the quantitative indexes are determined, determining the numerical values of the quantitative indexes respectively;
road grade A13Determining according to the traffic volume, the road use task and the property, wherein the higher the road grade is, the greater the influence degree is; the detailed division criteria are: the highway is a highway with particularly important political and economic significance, has four or more than four lanes, is provided with a central separation zone, is completely intersected and has perfect traffic safety facilities, management facilities and service facilities, and is a special highway specially used for high-speed driving of automobiles, wherein the ADT (average daily traffic volume in the year) is more than or equal to 25000; the first-level highway is a highway connecting the important political and economic culture center and part of overpasses, and ADT is more than or equal to 10000 and less than 25000; the secondary highway is a main highway connected with politics and economic centers or a great industrial and mining area or a suburban highway with busy transportation, and ADT is more than or equal to 2000 and less than 10000; the third level road is a branch road for communicating with a city in or above a county, and ADT is more than or equal to 200 and less than 2000; the four-level road is a branch road for communicating counties, towns and villages, and is ADT<200;
Functional area A traversed by road construction project21The method comprises the steps of determining according to technical rules of functional zoning of a natural protection zone, and dividing the natural protection zone into a core zone, a buffer zone and an experimental zone according to technical rules of functional zoning of the natural protection zone (GB/T35822-; wherein, the core area is a centralized distribution area of a natural ecosystem, rare or endangered wild animals and plants and natural vestige which are well preserved in a natural protection area; the buffer area is an area which is defined at the periphery of the core area and used for slowing down the interference of the outside world to the core area; the experimental area is an effective combination area of natural protection and sustainable resource utilization, and can be used for developing traditional production, scientific experiments, propaganda and education, ecological tourism, management service and natural recovery;
concerning distribution of endangered species A22According to the field investigation and determination of the variety of endangered species, the endangered species refers to IUCN and the extremely dangerous or endangered species in the red directory of the Chinese endangered species, and the wild animal and plant species are mainly protected in the first grade and the second grade of the country, the distribution of the endangered species is from scattered to concentrated, and the influence degree is stronger;
foreign invasive species condition A41According to the famous records of foreign invasive species in China (first batch and second batch), the data is consultedDetermining the state of the foreign invasive species caused by road construction in modes of access, field investigation and the like;
s200, judging related endangered species; according to actual investigation, related distribution condition A of endangered species22If IUCN and Chinese endangered species red directory are extremely dangerous, endangered and national level I key protection species centralized distribution areas along the road, the road construction project is unacceptable, and the level evaluation is finished; if IUCN is not involved along the road, and the red directory of the Chinese endangered species is extremely dangerous, endangered and the concentrated distribution area of the national class I key protection species, the method enters S300;
s300, determining the influence grade of the influence evaluation index; a is prepared by11、A12、A13、A21、A22、A31、A32、A33、A41And A42Dividing the ecological environment into three levels according to the influence degree on the ecology, wherein the three levels are level I, level II and level III respectively; wherein, level I shows that the ecological influence of the road project on the natural protection area is within an acceptable range, the influence level is green and the number of the influence level is 1, level II shows that the road project has obvious influence on the natural protection area, corresponding protection measures are needed to reduce the adverse influence, the influence level is yellow and the number of the influence level is 2, level III shows that the road project has serious adverse influence on the natural protection area, the influence level is red and the number of the influence level is 3;
road length A relating to natural conservation zones11Is divided into standard functions ofd1The formula for calculating (a) is as follows,
wherein d is1The influence strength of the road length in the natural protection area is obtained;
l is the length of the road, km, related to the natural reserve;
p is the perimeter of the natural protected area, km;
area A of natural protection area invaded by auxiliary facilities12Is divided into standard functions ofd2The formula for calculating (a) is as follows,
wherein d is2The influence strength of the area of the road and the accessory facilities thereof;
s is the area of the road and its subsidiary facilities invading the natural protection area hm2;
A is the total area of the natural reserve, hm2;
Road grade A13The division standard of (1) is that when the road construction project is a highway and a first-level road, the influence level is level III; the road construction projects are second-level and third-level roads, and the influence level is level II; the road construction project is a four-level road and other roads, and the influence level is level I;
functional area A traversed by road construction project21The division standard of (1) is that when the road construction project relates to a core area and a buffer area, the influence level is level III; when the road construction project relates to a test area, the influence level is level II; when the road construction project relates to other areas, the influence level is level I;
concerning distribution of endangered species A22The division standard is that when the road construction project relates to IUCN, the red directory of Chinese endangered species is extremely dangerous, the national class I species is scattered distribution area or the national class II centralized distribution area, the influence level is level III; when the road construction project relates to IUCN and red name records of endangered species in China, and the national level II is scattered, the influence grade is level II; when the road construction project does not relate to the distribution area of the species, the influence level is level I;
direct loss of primary protected object A31Is divided into standard functions ofWherein d is3The formula for calculating (a) is as follows,
wherein d is3Directly impact strength for the primary protected object;
when the main protected object is a plant, M is the direct loss amount of the plant individual number, the plant cluster number or the plant distribution area caused by the road construction project; when the main protected object is an animal, M is the habitat area directly invading the main protected object species;
when the main protected object is a plant, T is the plant individual number, the plant clump number or the total plant distribution area of the main protected object in the natural protected area; when the main protected object is an animal, T is the total habitat area of the main protected object in the natural protected area;
area of soil erosion32Is divided into standard functions ofWherein d is4Setting the percentage of the water and soil loss area caused by road construction in the range of 30m at both sides of the road in the natural protection area to the total area of 30m buffer areas at both sides of the road;
intensity of influence of disruption of habitat A33Is divided into standard functions ofd5The formula for calculating (a) is as follows,
wherein d is5For the environment breaking change rate, t is the road after the road construction, and t-1 is the roadBefore road construction;
Skis the area of the kth habitat patch, where k is a positive integer, hm2;
A is the total area of the natural reserve, hm2;
Foreign invasive species condition A41The division standard of (1) is that when external invasive species are found along the road construction project and are propagated in a large area, the influence grade is grade III; the road construction project is sporadically distributed with external invasive species along the line, but the influence level is level II when large-area diffusion is not achieved at present; when no foreign invasive species are found along the road construction project, the influence level is level I;
indirect loss of protected objects a42Is divided into standard functions ofd6The formula for calculating (a) is as follows,
wherein d is6Indirectly influencing the intensity for the main protected object;
r is the rate of change of the importance (or frequency of occurrence) of the main protected species in the community.
S400, determining a membership function of the influence evaluation index; determining the membership degree of each ecological influence evaluation index, wherein,
membership matrix R1For road attribute influence evaluation index A1Is 3 × 3 matrix, and the membership matrixMembership matrix R1The a-th row and the j-th column of (A) are denoted as raj,rajFor road attribute influence evaluation index A1The membership degree of the middle alpha influence evaluation index to the grade with the number j; wherein a is a positive integer and a is 1, 2, 3, j is a positive integer and j is 1, 2, 3;
degree of membershipMatrix R2Evaluation index A for susceptibility influence2Is a 2 × 3 matrix, a membership matrixMembership matrix R2The b-th row and j-th column element of (1) is denoted as rbj,rbjEvaluation index A for susceptibility influence2The membership degree of the middle b-th influence evaluation index to the grade with the number j; wherein b is a positive integer and b is 1, 2, 3, j is a positive integer and j is 1, 2, 3;
membership matrix R3To directly influence the evaluation index A3Is 3 × 3 matrix, and the membership matrixMembership matrix R3The element of row c and column j is denoted as rcj,rcjTo directly influence the evaluation index A3The membership degree of the c-th influence evaluation index to the grade with the number j; wherein c is a positive integer and c is 1, 2, 3, j is a positive integer and j is 1, 2, 3;
membership matrix R4As an indirect influence evaluation index A4Is a 2 × 3 matrix, a membership matrixMembership matrix R4The d-th row and j-th column of (1) are denoted as rdj,rdjAs an indirect influence evaluation index A4The membership degree of the d-th influence evaluation index to the grade with the number j; wherein d is a positive integer and d is 1, 2, 3, j is a positive integer and j is 1, 2, 3;
when determining the membership degrees of the 6 quantitative indexes in the S100, determining a membership degree function by adopting trapezoidal distribution, and calculating the membership degree of the grade numbered j of each quantitative index by using the determined membership degree function;
for road grade A13Is determined according to the road grade A determined in S20013Is determined by the grade of R1e1 and R1e'0; wherein e is a positive integer and e is the road grade A determined in S20013The grade number of (2); e ' is a positive integer, e ' ═ 1, 2, 3, and e ' ≠ e;
functional area A for road construction project crossing21Is determined according to the functional area A traversed by the road construction project determined in S20021Is determined by the grade of R1f1 and R1f'0; wherein f is a positive integer and f is the functional area A traversed by the road construction project determined in S20021The grade number of (2); f ' is a positive integer, f ' is 1, 2, 3, and f ' ≠ f;
for distribution conditions A involving endangered species22Is determined according to the distribution condition A of the related endangered species determined in S20022Is determined by the grade of R1g1 and R1g'0; wherein g is a positive integer and g is the distribution A related to endangered species determined in S20022The grade number of (2); g 'is a positive integer, g' ≠ g, 1, 2, 3;
for foreign invasive species status A41Is determined according to the determined foreign invasive species condition A in S20041Is determined by the grade of R1h1 and R1h'0; wherein h is a positive integer and h is the alien invasive species status A determined in S20041The grade number of (2); h 'is a positive integer, h' ≠ h, 1, 2, 3;
wherein the road attribute influences the evaluation index A1The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
area A of natural protection area invaded by auxiliary facilities12The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
direct loss of primary protected object A31The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
area of soil erosion32The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
intensity of influence of disruption of habitat A33The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
indirect loss of protected objects a42The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
s500, calculating the weight of the influence evaluation index; evaluation index A for influence of analytic hierarchy process on road attribute1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4The weight of the road is calculated, and the road attribute influence evaluation index A is calculated1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4The weights of (A) are respectively denoted as W1、W2、W3And W4To obtain a weight vector T, T ═ W1,W2,W3,W4);
Evaluation index A for influence on road attribute1Middle 3 ecological impact evaluation indexes A11、A12And A13Are calculated respectively, and A is calculated11、A12And A13The weights of (A) are respectively denoted as W11、W12And W13To obtain A1Weight vector T1,T1=(W11,W12,W13);
Evaluation index A for influence on sensitivity2Middle 2 ecological impact evaluation indexes A21And A22Are calculated respectively, and A is calculated21And A22The weights of (A) are respectively denoted as W21And W22To obtain A2Weight vector T2,T2=(W21,W22);
For direct influence evaluation index A3Middle 3 ecological impact evaluation indexes A31、A32And A33Are calculated respectively, and A is calculated31、A32And A33The weights of (A) are respectively denoted as W31、W32And W33To obtain A3Weight vector T3,T3=(W31,W32,W33);
Evaluation index A for indirect influence4Middle 2 ecological impact evaluation indexes A41And A42Are calculated respectively, and A is calculated41And A42The weights of (A) are respectively denoted as W42And W42To obtain A4Weight vector T4,T4=(W42,W42);
S600, fuzzy comprehensive evaluation; according to the membership degree matrix R determined in S4001、R2、R3And R4And the weight vector T determined in S5001、T2、T3And T4Calculating road attribute influence evaluation index A1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4Membership degree vector of (1), road attribute influence evaluation index A1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4Are respectively B1、B2、B3And B4(ii) a Wherein the road attribute influences the evaluation index A1The calculation formula of the membership degree vector is B1=T1R1(ii) a Susceptibility influence evaluation index A2The calculation formula of the membership degree vector is B2=T2R2(ii) a Direct influence evaluation index A3The calculation formula of the membership degree vector is B3=T3R3(ii) a Indirect influence evaluation index A4The calculation formula of the membership degree vector is B4=T4R4(ii) a According to B1、B2、B3And B4Obtaining an influence evaluation membership degree matrix R,calculating a fuzzy evaluation matrix B according to the formula B ═ TR, (B ═ TR)1,b2,b3);
S600, validity check is also included; the validity test is to test the fuzzy comprehensive evaluation, the calculation formula is,
wherein alpha is the maximum membership principle validity index;
β is b1,b2,b3Maximum value of (1);
gamma is b1,b2,b3The second largest value of;
when alpha is less than 0.5, evaluating the influence degree of the road construction project on the natural protection area by adopting a weighted average principle; otherwise, the influence degree of the road construction project on the natural protection area is evaluated by adopting the maximum membership principle.
S700, judging the ecological influence level; evaluating the influence degree of the road construction project on the natural protection area by using the fuzzy evaluation matrix B; b1、b2、b3Respectively corresponding to the evaluated specific gravities of the grade I, the grade II and the grade III; level I indicates that the ecological influence of the road project on the natural conservation area is within an acceptable range; level II shows that the road project has obvious influence on the natural protected area, and protective measures need to be taken to relieve the adverse influence; grade iii indicates that the road project has a serious adverse effect on the natural reserve, which is unacceptable.
It is worth mentioning that the indirect loss of protected objects A42The protected object in (1) refers to a specific protected object species in a natural protected area, a rare or endangered species or a main constituent species in a representative ecosystem. Except for a protection area with specific protection object species, other protection object protection objects in the protection area generally refer to rare endangered species and a representative ecosystem, and particularly take species habitat as a core, for example, the fir in Qinling mountains does not belong to the endangered species, but the ecosystem of the fir in Qinling mountains belongs to the representative ecosystem, and the protection object is measured according to the change of important values of the fir in Qinling mountains during specific measurement and calculation.
The invention also comprises an ecological impact index test for verifying the evaluation result, the calculation formula is,
the EI is an ecological influence index of a road construction project;
Diis AiI is an integer and i is greater than or equal to 1 and less than or equal to 4;
wiis AiI is an integer and i is greater than or equal to 1 and less than or equal to 4.
Further description of wiDerived from analytic hierarchy process and expert consulting method, and w1=0.2,w2=0.3,w3=0.25,w40.25. The evaluation criteria for EI are: when the EI is less than or equal to 50, the ecological influence level of the road construction project is green; when EI is more than 50 and less than 80, the ecological influence level of the road construction project is yellow, and when EI is more than or equal to 80, the ecological influence level of the road construction project is red.
DiIncluding a road attribute index D1Sensitivity index D2Direct influence index D3And an indirect influence index D4;
Road attribute index D1The formula for calculating (a) is as follows,
wherein D is1Is a road attribute index;
d1the influence strength of the road length in the natural protection area is obtained;
d2the influence strength of the area of the road and the accessory facilities thereof;
x is the road grade influence strength; grade A of the current road13When the classified influence grade is grade III, x is 100; grade A of the current road13When the classified influence grade is grade II, x is 61.8; grade A of the current road13When the classified influence grade is grade I, x is 38.2;
sensitivity index D2The formula for calculating (a) is as follows,
wherein D is2Is a sensitivity index;
y is the influence strength of the related functional regions; functional area A crossed by road construction project21When the classification level of (1) is level III, y is 100; functional area A crossed by road construction project21When the grade of (1) is grade II, y is 38.2; functional area A crossed by road construction project21When the classification level of (1) is level I, y is 0;
z is the impact intensity related to endangered species; when concerning distribution of endangered species A22When the classification level of (1) is level III, z is 100; when concerning distribution of endangered species A22When the grade is II, z is 61.8; when concerning distribution of endangered species A22When the grade is I, z is 38.2;
direct influence index D3The formula for calculating (a) is as follows,
wherein D is3Is a direct impact index;
d3directly impact strength for the primary protected object;
d4setting the percentage of the water and soil loss area caused by road construction in the range of 30m at both sides of the road in the natural protection area to the total area of 30m buffer areas at both sides of the road;
d5to an environmental disruption index;
index of indirect influence D4The formula for calculating (a) is as follows,
wherein D is4Is an indirect influence index;
n is the foreign invasive species impact strength; when foreign invaded species status A41When the classification of (1) is class III, n100; when foreign invaded species status A41When the grade of (1) is grade II, n is 61.8; when foreign invaded species status A41When the classification level of (1) is level I, n is 0;
d6the intensity is indirectly influenced for the main protected object.
In this embodiment, taking an example of a country-level natural protection area of a blimp ray lawn, the evaluation system provided by the present invention is used for evaluating the ecological impact of a newly added road construction project on the natural protection area after carrying out the ecological impact, and the method for evaluating the ecological impact of the road construction project on the natural protection area includes the following steps:
s100, determining an evaluation index;
s200, judging related endangered species; the road construction project does not relate to IUCN and the concentrated distribution area of the red famous and extremely dangerous, endangered and national I-level key protection species in China, and the road construction project enters S300;
s300, determining the influence grades of the influence evaluation indexes as a grade I, a grade II and a grade III;
calculating the road length based on ArcGIS according to the extracted new road vector data, wherein the road length is 15.38km, and the road length A related to the natural protection area is combined11Will relate to the road length A of the natural reserve11The influence grade of (2) is classified into grade II;
through actual investigation, the length of the road section with the width of 4m of the newly added road is 12.26km, and the length of the road section with the width of 6m is 3.12km, so that the area of the newly added road is 12.26km + 4m +3.12km + 6m which is 6.68hm2The area of the road auxiliary facilities (road maintenance station) is 0.06hm2Combined with subsidiary facilities encroaching on the area A of the natural protected area12Dividing standard function of (1) to encroach subsidiary facilities into natural protected area A12The influence grade of (2) is classified into grade II;
according to the traffic volume, the road use task and the property, the newly added road in the natural protection area is a branch road which is communicated with villages and towns in the packet hometown range, has small traffic flow and belongs to the category of four-level roads; combined with road grade a13The division standard of (1) road class A13The influence grade of (A) is classified into grade I;
extracting road construction information of the natural protection area by comparing the remote sensing images in 2014 and 2017; the extracted road vector data are superposed to a functional zoning map of the natural protection zone, namely, the road construction project passes through a buffer zone of the natural protection zone by combining the map 2; functional area A for crossing in combination with road construction project21The division standard of (A) is a functional area A which is crossed by the road construction project21The influence grade of (1) is classified into grade III;
the newly added road passes through a scattered distribution area of the second-level key protection wild plant species Dabie mountain spikenard of the main protection object country near the threshold ridge of the natural protection area; binding to distribution status related to endangered species A22Will relate to distribution of endangered species A22The influence grade of (2) is classified into grade II;
road construction projects penetrate scattered distribution areas of main protection objects of the macrolepis gigantea in natural protection areas, and investigation proves that the individual macrolepis gigantea is not directly invaded in the project construction process, and direct loss of the main protection object of the macrolepis gigantea is not caused; direct loss of bound primary protected object A31Will directly lose a primary protected object by a division criterion of31The influence grade of (A) is classified into grade I;
the construction of the road can cause water and soil loss on two sides of the roadside, and particularly, the water and soil loss is more serious when the slopes of two sides of the road are larger; by visual interpretation of the remote sensing image, as shown in FIG. 3, and combined with field investigation, the loss area of soil and water within 30m of buffer on both sides of the newly added road is calculated to be 2.21hm2The area of the 30m area of the buffer on the two sides is 15.38 km/0.06 km/92.28 hm2(ii) a Combined soil and water loss area A32Dividing the soil erosion area A into standard functions32The influence grade of (2) is classified into grade II;
making a natural habitat distribution map of a natural protected area before and after road construction, as shown in fig. 4, before road project construction, the number of natural habitat patches is 4, and patch areas are respectively: 5194.26hm2、3713.85hm2、1251.46hm2、1065.7hm2According to the fragmentation algorithm, the fragmentation index Fi is calculated to be 0.72. Method for constructing roadsAnd then, the number of the original natural habitat patches is increased to 5, and the areas are respectively as follows: 3297.04hm2、1897.17hm2、3713.85hm2、1251.46hm2、1065.7hm2The fragmentation index FI is 0.81; after a road is newly built, the fragmentation of the natural habitat of the natural protection area is increased by 12.5 percent; effect of Combined habitat disruption A33By breaking the habitat into influence intensities A33The influence grade of (2) is classified into grade II;
determining the state of the foreign invasive species through field investigation and by combining interview data; no foreign invasive species are found along the road construction project; binding to foreign invasive species Condition A41The division criteria of (2) is to classify the foreign invasive species status A41The influence grade of (A) is classified into grade I;
only 1 Dabie mountain spikenard is distributed around the newly-repaired road near the threshold ridge, survey samples are arranged at a distance of 20m and 100m from the road, and the Dabie mountain spikenard of the protected object is not recorded, so that the change of the important value of the Dabie mountain spikenard is not caused; however, considering that the giant salamander which is a national second-level key protection animal is also a main protection object in a natural protection area, according to interview data of local residents, after the road is built, because the human activities along the road are increased, the frequency of finding the giant salamander in a nearby river channel is reduced compared with that before the road is built, and the indirect loss A of the protection object is combined42Dividing the standard function into a number of blocks, and indirectly losing the protected object by A42The influence grade of (A) is classified into grade I;
s400, determining a membership function of the influence evaluation index;
s500, calculating the weight of the influence evaluation index; evaluation index A for influence of analytic hierarchy process on road attribute1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4Calculating the weight of (a); and a11、A12、A13、A21、A22、A31、A32、A33、A41And A42The weights of (a) are calculated respectively;
s600, fuzzy comprehensive evaluation; in the validity test, the maximum membership principle validity index alpha is more than or equal to 0.5, and a fuzzy evaluation matrix B is obtained;
s700, judging the ecological influence level; evaluating the influence degree of the road construction project on the natural protection area by using the fuzzy evaluation matrix B; the evaluation result shows that the ecological influence level of the road construction project is yellow; the evaluation result is checked in combination with an ecological impact index, wherein the road attribute index D166.82, sensitivity index D283.12, direct impact index D366.62, indirect influence index D4When 43.69, the ecological impact level of the road construction project was "yellow" as a result of the evaluation. The results of the two are the same, which indicates that the newly added road construction project in the protected area has obvious influence on the natural protected area, and corresponding protective measures need to be taken to alleviate the adverse influence. Such as: the new construction of the road causes the disruption of the natural protected area to be intensified, the blocking effect is obvious, and animal channels are established to slow down the blocking effect according to the distribution characteristics of animals along the road; the habitat fragmentation aggravates, has also increased the edge effect simultaneously, and the possibility that the resident of peripheral community got into the protected area through this road and carried out the interference increases, and the probability of discovering the giant salamander like near river course reduces, consequently should strengthen daily patrol and protect, improves the management and maintenance dynamics along the road line.
Example 2
The basic contents of this embodiment are different from those of embodiment 1 in that the calculation of the influence evaluation index weight in S500 is a comprehensive calculation of the influence evaluation index weight, and the steps are as follows,
s501, subjective weight calculation of the impact evaluation indexes; evaluation index A for influence on road attribute1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4The subjective weights of (a) are calculated to obtain a subjective weight vector C, C ═ μ1,μ2,μ3) (ii) a The constructed judgment matrix is a first judgment matrix D1First decision matrix D1Is a 4-order matrix;
evaluation index A for influence on road attribute1Middle 3 ecological impact evaluation indexes A11、A12And A13Respectively calculating the subjective weights of A and B11、A12And A13The subjective weights of (A) are respectively denoted as mu11、μ12And mu13To obtain A1Subjective weight vector C1,C1=(μ11,μ12,μ13) (ii) a The constructed judgment matrix is a second judgment matrix D2Second decision matrix D2Is a 3-order matrix;
evaluation index A for influence on sensitivity2Middle 2 ecological impact evaluation indexes A21And A22Respectively calculating the subjective weights of A and B21And A22The subjective weights of (A) are respectively denoted as mu22And mu22To obtain A2Subjective weight vector C2,C2=(μ22,μ22) (ii) a The constructed judgment matrix is a third judgment matrix D3The third judgment matrix D3Is a 2-order matrix;
for direct influence evaluation index A3Middle 3 ecological impact evaluation indexes A31、A32And A33Respectively calculating the subjective weights of A and B31、A32And A33The subjective weights of (A) are respectively denoted as mu31、μ32And mu33To obtain A3Subjective weight vector C3,C3=(μ31,μ32,μ33) (ii) a The constructed judgment matrix is a fourth judgment matrix D4Fourth decision matrix D4Is a 3-order matrix;
evaluation index A for indirect influence4Middle 2 ecological impact evaluation indexes A41And A42Respectively calculating the subjective weights of A and B41And A42The subjective weights of (A) are respectively denoted as mu42And mu42To obtain A4Subjective weight vector C4,C4=(μ42,μ42) (ii) a The constructed judgment matrix is a fifth judgment matrix D5Fifth decision matrix D5Is a 2-order matrix;
s502, influencing objective weight of evaluation indexCalculating; using the first determination matrix D1Calculating to obtain road attribute influence evaluation index A by adopting an entropy weight method1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4Objective weight λ of1、λ2、λ3And λ4;
Using the second decision matrix D2Evaluation index A for influence of entropy weight method on road attribute1Middle 3 ecological impact evaluation indexes A11、A12And A13Respectively calculating the objective weights of A and B11、A12And A13The objective weights of (A) are respectively denoted as λ11、λ12And λ13;
Using the third determination matrix D3Evaluation index A for sensitivity influence by entropy weight method2Middle 2 ecological impact evaluation indexes A21And A22Respectively calculating the objective weights of A and B21And A22The objective weights of (A) are respectively denoted as λ22And λ22;
Using the fourth determination matrix D4Evaluating index A by direct influence of entropy weight method3Middle 3 ecological impact evaluation indexes A31、A32And A33Respectively calculating the objective weights of A and B31、A32And A33The objective weights of (A) are respectively denoted as λ31、λ32And λ33;
Using the fifth decision matrix D5Evaluation index A for indirect influence by entropy weight method4Middle 2 ecological impact evaluation indexes A41And A42Respectively calculating the objective weights of A and B41And A42The objective weights of (A) are respectively denoted as λ42And λ42;
S503, calculating the combined weight of the influence evaluation indexes; calculation of A1、A2、A3And A4Combining weight W of1、W2、W3And W4To obtain a weight vector T, T ═ W1,W2,W3,W4);
The calculation formula is as follows,
Wm=α1λm+(1-α1)μm
wherein,
Wmis AmM is a positive integer and m is more than or equal to 1 and less than or equal to 4;
μmis AmM is a positive integer and is more than or equal to 1 and less than or equal to 4;
λmis AmM is a positive integer and is more than or equal to 1 and less than or equal to 4;
Calculating road attribute influence evaluation index A1Combining weight W of11、W12And W13To obtain A1Weight vector T1,T1=(W11,W12,W13);
The calculation formula is as follows,
W1m=α2λ1m+(1-α2)μ1m
wherein,
W1mis A1mM is a positive integer and m is more than or equal to 1 and less than or equal to 3;
μ1mis A1mM is a positive integer and is more than or equal to 1 and less than or equal to 3;
λ1mis A1mM is a positive integer and is more than or equal to 1 and less than or equal to 3;
Calculation of sensitivity influence evaluation index A2Combining weight W of21And W22To obtain A2Weight vector T2,T2=(W21,W22);
The calculation formula is as follows,
W2m=α3λ2m+(1-α3)μ2m
wherein,
W2mis A2mM is a positive integer and m is more than or equal to 1 and less than or equal to 2;
μ2mis A2mM is a positive integer and is more than or equal to 1 and less than or equal to 2;
λ2mis A2mM is a positive integer and is more than or equal to 1 and less than or equal to 2;
α3is A2And α3=W21+W22-1.5;
Calculating direct influence evaluation index A3Combining weight W of31、W32And W33To obtain A3Weight vector T3,T3=(W31,W32,W33);
The calculation formula is as follows,
W3m=α2λ3m+(1-α3)μ3m
wherein,
W3mis A3mM is a positive integer and m is more than or equal to 1 and less than or equal to 3;
μ3mis A3mM is a positive integer and is more than or equal to 1 and less than or equal to 3;
λ3mis A3mM is a positive integer and is more than or equal to 1 and less than or equal to 3;
Calculating an indirect influence evaluation index A4Combining weight W of42And W42To obtain A4Weight vector T4,T4=(W42,W42);
The calculation formula is as follows,
W4m=α2λ4m+(1-α3)μ4m
wherein,
W4mis A4mM is a positive integer and m is more than or equal to 1 and less than or equal to 2;
μ4mis A4mM is a positive integer and is more than or equal to 1 and less than or equal to 2;
λ4mis A4mM is a positive integer and is more than or equal to 1 and less than or equal to 2;
α5is A4And α5=W41+W42-1.5。
In this embodiment, taking an Anhui ray lawn country-level natural protection area as an example, the evaluation system provided in this embodiment is used for evaluating ecological influence on a newly-added built road project in the natural protection area, the ecological influence of the evaluation road construction project on the natural protection area is large or small, the evaluation result shows that the ecological influence level of the road construction project is yellow, the evaluation result is tested by combining with an ecological influence index, and the results are the same, which indicates that the method can quickly evaluate the overall ecological influence of the road construction project, and is short in evaluation time, objective and accurate in evaluation result and wide in application prospect.
Finally, it should be noted that: it should be understood that the above examples are only for clearly illustrating the present application and are not intended to limit the embodiments. Other variations and modifications will be apparent to persons skilled in the art in light of the above description. And are neither required nor exhaustive of all embodiments. And obvious variations or modifications of this type are intended to be covered by the present invention.
Claims (10)
1. A natural conservation area road construction project ecological influence assessment system is characterized by comprising the following steps:
s100, determining an evaluation index; for road constructionDetermining ecological influence evaluation indexes of the natural protection area by the project; the ecological influence evaluation index comprises a road attribute influence evaluation index A1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4;
The road engineering influence evaluation index A1Including road length A relating to natural reserve11Area A of natural protection area invaded by road and accessory facilities12And road grade A13Wherein the road length A11Area A of natural protection area invaded by road and accessory facilities12Are all quantitative indexes, road grade A13Is a qualitative index; the sensitivity influence evaluation index A2Functional area A including road construction project traversal21And to distribution of endangered species A22Wherein the road construction project passes through the functional area A21And to distribution of endangered species A22Are all qualitative indexes; the direct influence evaluation index A3Including the direct loss of a primary protected object31Area of soil erosion32And disruption of habitat influence intensity A33Wherein the primary protected object directly loses A31Area of soil erosion32And disruption of habitat influence intensity A33Are all quantitative indicators; the indirect influence evaluation index A4Including foreign invasive species status A41And indirect loss of protected objects A42Wherein the foreign invasive species status A41For qualitative index, indirect loss of protected object A42Is a quantitative index; when the quantitative indexes are determined, determining the numerical values of the quantitative indexes respectively;
wherein the road class A13Determining according to the traffic volume, the road use task and the property, wherein the higher the road grade is, the greater the influence degree is; the distribution condition A of the species involved in endangerment22According to the endangered species diversity field investigation and determination, the distribution is from scattered to concentrated, and the influence degree is stronger; said alien invasive species condition A41Data consulting, accessing and field investigation and determination are adopted;
s200, judging related endangered species; according to actual investigation, related distribution condition A of endangered species22If IUCN is related along the road and the red famous and extremely dangerous species, the red famous and extremely dangerous species and the red famous and extremely dangerous species in China are intensively distributed, the road construction project is not acceptable, and the grade evaluation is finished; if IUCN is not involved along the road, and the red directory of the Chinese endangered species is extremely dangerous, endangered and the concentrated distribution area of the national class I key protection species, the method enters S300;
s300, determining the influence grade of the influence evaluation index; a is prepared by11、A12、A13、A21、A22、A31、A32、A33、A41And A42Dividing the ecological environment into three levels according to the influence degree on the ecology, wherein the three levels are level I, level II and level III respectively; wherein, level I shows that the ecological influence of the road project on the natural protection area is within an acceptable range, the influence level is green and the number of the influence level is 1, level II shows that the road project has obvious influence on the natural protection area, corresponding protection measures are needed to reduce the adverse influence, the influence level is yellow and the number of the influence level is 2, level III shows that the road project has serious adverse influence on the natural protection area, the influence level is red and the number of the influence level is 3;
s400, determining a membership function of the influence evaluation index; determining the membership degree of each ecological influence evaluation index, wherein,
membership matrix R1For road attribute influence evaluation index A1Is a 3 × 3 matrix, the membership matrixMembership matrix R1The a-th row and the j-th column of (A) are denoted as raj,rajFor road attribute influence evaluation index A1The membership degree of the middle alpha influence evaluation index to the grade with the number j; wherein a is a positive integer and a is 1, 2, 3, j is a positive integer and j is 1, 2, 3;
membership matrix R2To be sensitiveEvaluation index of influence A2Is a 2 × 3 matrix, the membership matrixMembership matrix R2The b-th row and j-th column element of (1) is denoted as rbj,rbjEvaluation index A for susceptibility influence2The membership degree of the middle b-th influence evaluation index to the grade with the number j; wherein b is a positive integer and b is 1, 2, 3, j is a positive integer and j is 1, 2, 3;
membership matrix R3To directly influence the evaluation index A3Is a 3 × 3 matrix, the membership matrixMembership matrix R3The element of row c and column j is denoted as rcj,rcjTo directly influence the evaluation index A3The membership degree of the c-th influence evaluation index to the grade with the number j; wherein c is a positive integer and c is 1, 2, 3, j is a positive integer and j is 1, 2, 3;
membership matrix R4As an indirect influence evaluation index A4Is a 2 × 3 matrix, the membership matrixMembership matrix R4The d-th row and j-th column of (1) are denoted as rdj,rdjAs an indirect influence evaluation index A4The membership degree of the d-th influence evaluation index to the grade with the number j; wherein d is a positive integer and d is 1, 2, 3, j is a positive integer and j is 1, 2, 3;
when determining the membership degrees of the 6 quantitative indexes in the S100, determining a membership degree function by adopting trapezoidal distribution, and calculating the membership degree of each quantitative index with the grade of j by using the determined membership degree function;
for road grade A13Is determined according to the road grade A determined in S20013Grade (1) ofIs determined, R1e1 and R1e'0; wherein e is a positive integer and e is the road grade A determined in S20013The grade number of (2); e ' is a positive integer, e ' ═ 1, 2, 3, and e ' ≠ e;
functional area A for road construction project crossing21Is determined according to the functional area A traversed by the road construction project determined in S20021Is determined by the grade of R1f1 and R1f'0; wherein f is a positive integer and f is the functional area A traversed by the road construction project determined in S20021The grade number of (2); f ' is a positive integer, f ' is 1, 2, 3, and f ' ≠ f;
for distribution conditions A involving endangered species22Is determined according to the distribution condition A of the related endangered species determined in S20022Is determined by the grade of R1g1 and R1g'0; wherein g is a positive integer and g is the distribution A related to endangered species determined in S20022The grade number of (2); g 'is a positive integer, g' ≠ g, 1, 2, 3;
for foreign invasive species status A41Is determined according to the determined foreign invasive species condition A in S20041Is determined by the grade of R1h1 and R1h'0; wherein h is a positive integer and h is the alien invasive species status A determined in S20041The grade number of (2); h 'is a positive integer, h' ≠ h, 1, 2, 3;
s500, calculating the weight of the influence evaluation index; evaluation index A for influence of the road attribute1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4The weight of the influence evaluation index is calculated, and the method for calculating the weight of the influence evaluation index is an analytic hierarchy process; and the road attribute influence evaluation index A1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4The weights of (A) are respectively denoted as W1、W2、W3And W4To obtain a weight vector T, T ═ W1,W2,W3,W4);
Evaluation index A for influence on road attribute1Middle 3 ecological impact evaluation indexes A11、A12And A13Are calculated respectively, and A is calculated11、A12And A13The weights of (A) are respectively denoted as W11、W12And W13To obtain A1Weight vector T1,T1=(W11,W12,W13);
Evaluation index A for influence on sensitivity2Middle 2 ecological impact evaluation indexes A21And A22Are calculated respectively, and A is calculated21And A22The weights of (A) are respectively denoted as W21And W22To obtain A2Weight vector T2,T2=(W21,W22);
For direct influence evaluation index A3Middle 3 ecological impact evaluation indexes A31、A32And A33Are calculated respectively, and A is calculated31、A32And A33The weights of (A) are respectively denoted as W31、W32And W33To obtain A3Weight vector T3,T3=(W31,W32,W33);
Evaluation index A for indirect influence4Middle 2 ecological impact evaluation indexes A41And A42Are calculated respectively, and A is calculated41And A42The weights of (A) are respectively denoted as W42And W42To obtain A4Weight vector T4,T4=(W42,W42);
S600, fuzzy comprehensive evaluation; according to the membership degree matrix R determined in S4001、R2、R3And R4And the weight vector T determined in S5001、T2、T3And T4Calculating road attribute influence evaluation index A1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4The road attribute influences the evaluation index A1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4Are respectively B1、B2、B3And B4(ii) a Wherein the road attribute influences the evaluation index A1The calculation formula of the membership degree vector is B1=T1R1(ii) a Susceptibility influence evaluation index A2The calculation formula of the membership degree vector is B2=T2R2(ii) a Direct influence evaluation index A3The calculation formula of the membership degree vector is B3=T3R3(ii) a Indirect influence evaluation index A4The calculation formula of the membership degree vector is B4=T4R4(ii) a According to B1、B2、B3And B4Obtaining an influence evaluation membership degree matrix R,calculating a fuzzy evaluation matrix B according to the formula B ═ TR, (B ═ TR)1,b2,b3);
S700, judging the ecological influence level; evaluating the influence degree of the road construction project on the natural protection area by using the fuzzy evaluation matrix B; b1、b2、b3Respectively corresponding to the evaluated specific gravities of the grade I, the grade II and the grade III; level I indicates that the ecological influence of the road project on the natural conservation area is within an acceptable range; level II shows that the road project has obvious influence on the natural protected area, and protective measures need to be taken to relieve the adverse influence; grade iii indicates that the road project has a serious adverse effect on the natural reserve, which is unacceptable.
2. The ecological impact assessment system for road construction projects in natural reserve according to claim 1, characterized in that: the calculation of the influence evaluation index weight in S500 is an influence evaluation index weight comprehensive calculation, and includes the steps of,
s501, influence evaluationCalculating subjective weight of price index; evaluation index A for influence on road attribute1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4The subjective weights of (a) are calculated to obtain a subjective weight vector C, C ═ μ1,μ2,μ3) (ii) a The constructed judgment matrix is a first judgment matrix D1The first judgment matrix D1Is a 4-order matrix;
evaluation index A for influence on road attribute1Middle 3 ecological impact evaluation indexes A11、A12And A13Respectively calculating the subjective weights of A and B11、A12And A13The subjective weights of (A) are respectively denoted as mu11、μ12And mu13To obtain A1Subjective weight vector C1,C1=(μ11,μ12,μ13) (ii) a The constructed judgment matrix is a second judgment matrix D2The second determination matrix D2Is a 3-order matrix;
evaluation index A for influence on sensitivity2Middle 2 ecological impact evaluation indexes A21And A22Respectively calculating the subjective weights of A and B21And A22The subjective weights of (A) are respectively denoted as mu22And mu22To obtain A2Subjective weight vector C2,C2=(μ22,μ22) (ii) a The constructed judgment matrix is a third judgment matrix D3The third determination matrix D3Is a 2-order matrix;
for direct influence evaluation index A3Middle 3 ecological impact evaluation indexes A31、A32And A33Respectively calculating the subjective weights of A and B31、A32And A33The subjective weights of (A) are respectively denoted as mu31、μ32And mu33To obtain A3Subjective weight vector C3,C3=(μ31,μ32,μ33) (ii) a The constructed judgment matrix is a fourth judgment matrix D4The fourth determination matrix D4Is a 3-order matrix;
evaluation index A for indirect influence4Middle 2 ecological impact evaluation indexes A41And A42Respectively calculating the subjective weights of A and B41And A42The subjective weights of (A) are respectively denoted as mu42And mu42To obtain A4Subjective weight vector C4,C4=(μ42,μ42) (ii) a The constructed judgment matrix is a fifth judgment matrix D5The fifth judgment matrix D5Is a 2-order matrix;
s502, calculating objective weight of the influence evaluation index; using the first determination matrix D1Calculating to obtain road attribute influence evaluation index A by adopting an entropy weight method1Sensitivity influence evaluation index A2Evaluation index A of direct influence3And indirect influence evaluation index A4Objective weight λ of1、λ2、λ3And λ4;
Using the second decision matrix D2Evaluation index A for influence of entropy weight method on road attribute1Middle 3 ecological impact evaluation indexes A11、A12And A13Respectively calculating the objective weights of A and B11、A12And A13The objective weights of (A) are respectively denoted as λ11、λ12And λ13;
Using the third determination matrix D3Evaluation index A for sensitivity influence by entropy weight method2Middle 2 ecological impact evaluation indexes A21And A22Respectively calculating the objective weights of A and B21And A22The objective weights of (A) are respectively denoted as λ22And λ22;
Using the fourth determination matrix D4Evaluating index A by direct influence of entropy weight method3Middle 3 ecological impact evaluation indexes A31、A32And A33Respectively calculating the objective weights of A and B31、A32And A33The objective weights of (A) are respectively denoted as λ31、λ32And λ33;
Using the fifth decision matrix D5Evaluation index A for indirect influence by entropy weight method4Middle 2 ecological impact evaluation indexes A41And A42Respectively calculating the objective weights of A and B41And A42The objective weights of (A) are respectively denoted as λ42And λ42;
S503, calculating the combined weight of the influence evaluation indexes; calculation of A1、A2、A3And A4Combining weight W of1、W2、W3And W4To obtain a weight vector T, T ═ W1,W2,W3,W4);
The calculation formula is as follows,
Wm=α1λm+(1-α1)μm
wherein,
Wmis AmM is a positive integer and m is more than or equal to 1 and less than or equal to 4;
μmis AmM is a positive integer and is more than or equal to 1 and less than or equal to 4;
λmis AmM is a positive integer and is more than or equal to 1 and less than or equal to 4;
Calculating road attribute influence evaluation index A1Combining weight W of11、W12And W13To obtain A1Weight vector T1,T1=(W11,W12,W13);
The calculation formula is as follows,
W1m=α2λ1m+(1-α2)μ1m
wherein,
W1mis A1mM is a positive integer and m is more than or equal to 1 and less than or equal to 3;
μ1mis A1mM is a positive integer and is more than or equal to 1 and less than or equal to 3;
λ1mis A1mM is a positive integer and is more than or equal to 1 and less than or equal to 3;
Calculation of sensitivity influence evaluation index A2Combining weight W of21And W22To obtain A2Weight vector T2,T2=(W21,W22);
The calculation formula is as follows,
W2m=α3λ2m+(1-α3)μ2m
wherein,
W2mis A2mM is a positive integer and m is more than or equal to 1 and less than or equal to 2;
μ2mis A2mM is a positive integer and is more than or equal to 1 and less than or equal to 2;
λ2mis A2mM is a positive integer and is more than or equal to 1 and less than or equal to 2;
α3is A2And α3=W21+W22-1.5;
Calculating direct influence evaluation index A3Combining weight W of31、W32And W33To obtain A3Weight vector T3,T3=(W31,W32,W33);
The calculation formula is as follows,
W3m=α2λ3m+(1-α3)μ3m
wherein,
W3mis A3mM is a positive integer and m is more than or equal to 1 and less than or equal to 3;
μ3mis A3mM is a positive integer and is more than or equal to 1 and less than or equal to 3;
λ3mis A3mM is a positive integer and is more than or equal to 1 and less than or equal to 3;
Calculating an indirect influence evaluation index A4Combining weight W of42And W42To obtain A4Weight vector T4,T4=(W42,W42);
The calculation formula is as follows,
W4m=α2λ4m+(1-α3)μ4m
wherein,
W4mis A4mM is a positive integer and m is more than or equal to 1 and less than or equal to 2;
μ4mis A4mM is a positive integer and is more than or equal to 1 and less than or equal to 2;
λ4mis A4mM is a positive integer and is more than or equal to 1 and less than or equal to 2;
α5is A4And α5=W41+W42-1.5。
3. The ecological impact assessment system for road construction projects in natural reserve according to claim 1, characterized in that: the S600 also comprises validity check; the validity test is to test the comprehensive evaluation of the fuzzy, the calculation formula is,
wherein alpha is the maximum membership principle validity index;
β is b1,b2,b3Maximum value of (1);
gamma is b1,b2,b3The second largest number inA value;
when alpha is less than 0.5, evaluating the influence degree of the road construction project on the natural protection area by adopting a weighted average principle; otherwise, the influence degree of the road construction project on the natural protection area is evaluated by adopting the maximum membership principle.
4. The ecological impact assessment system for road construction projects in natural reserve according to claim 1, characterized in that: the road length A related to the natural reserve11Is divided into standard functions ofd1The formula for calculating (a) is as follows,
wherein d is1The influence strength of the road length in the natural protection area is obtained;
l is the length of the road, km, related to the natural reserve;
p is the perimeter of the natural protected area, km;
the area A of the natural protection area invaded by the auxiliary facilities12Is divided into standard functions ofd2The formula for calculating (a) is as follows,
wherein d is2The influence strength of the area of the road and the accessory facilities thereof;
s is the area of the road and its subsidiary facilities invading the natural protection area hm2;
A is the total area of the natural reserve, hm2;
The road grade A13The division standard is that when the road construction project is expressway, first class highwayOn the road, the influence grade is grade III; the road construction projects are second-level and third-level roads, and the influence level is level II; the road construction project is a four-level road and other roads, and the influence level is level I;
functional area A traversed by road construction project21The division standard of (1) is that when the road construction project relates to a core area and a buffer area, the influence level is level III; when the road construction project relates to a test area, the influence level is level II; when the road construction project relates to other areas, the influence level is level I;
the distribution condition A of the species involved in endangerment22The division standard is that when the road construction project relates to IUCN, the red directory of Chinese endangered species is extremely dangerous, the national class I species is scattered distribution area or the national class II centralized distribution area, the influence level is level III; when the road construction project relates to IUCN and red name records of endangered species in China, and the national level II is scattered, the influence grade is level II; when the road construction project does not relate to the distribution area of the species, the influence level is level I;
the primary protected object directly loses A31Is divided into standard functions ofWherein d is3The formula for calculating (a) is as follows,
wherein d is3Directly impact strength for the primary protected object;
when the main protected object is a plant, M is the direct loss amount of the plant individual number, the plant cluster number or the plant distribution area caused by the road construction project; when the main protected object is an animal, M is the habitat area directly invading the main protected object species;
when the main protected object is a plant, T is the plant individual number, the plant clump number or the total plant distribution area of the main protected object in the natural protected area; when the main protected object is an animal, T is the total habitat area of the main protected object in the natural protected area;
the area of soil erosion32Is divided into standard functions ofWherein d is4Setting the percentage of the water and soil loss area caused by road construction in the range of 30m at both sides of the road in the natural protection area to the total area of 30m buffer areas at both sides of the road;
the disruption of the habitat affects intensity A33Is divided into standard functions ofd5The formula for calculating (a) is as follows,
wherein d is5For the environment fragmentation change rate, t is after road construction, and t-1 is before road construction;
FI is an environmental disruption index;
Skis the area of the kth habitat patch, where k is a positive integer, hm2;
A is the total area of the natural reserve, hm2;
Said alien invasive species condition A41The division standard of (1) is that when external invasive species are found along the road construction project and are propagated in a large area, the influence grade is grade III; the road construction project is sporadically distributed with external invasive species along the line, but the influence level is level II when large-area diffusion is not achieved at present; when no foreign invasive species are found along the road construction project, the influence level is level I;
the indirect loss of protected object A42The division criterion function of (a) is,d6the formula for calculating (a) is as follows,
wherein d is6Indirectly influencing the intensity for the main protected object;
r is the rate of change of the importance (or frequency of occurrence) of the main protected species in the community.
5. The ecological impact assessment system for road construction projects in natural reserve according to claim 1, characterized in that: the road grade A13The standard for the classification of the road grade in (1) is,
the highway is a highway with particularly important political and economic significance, is provided with four or more lanes, is provided with a central separation zone, is completely intersected and has perfect traffic safety facilities, management facilities and service facilities, and is a special highway specially used for high-speed driving of automobiles, wherein the ADT is more than or equal to 25000;
the first-level highway is a highway connecting the important political and economic culture center and part of overpasses, and ADT is more than or equal to 10000 and less than 25000;
the secondary highway is a main highway connected with politics and economic centers or a great industrial and mining area or a suburban highway with busy transportation, and ADT is more than or equal to 2000 and less than 10000;
the third level road is a branch road for communicating with a city in or above a county, and ADT is more than or equal to 200 and less than 2000;
the four-level road is a branch road for communicating counties, towns and villages, and ADT < 200.
6. The ecological impact assessment system for road construction projects in natural reserve according to claim 4, characterized in that: in S400, the road attribute influence evaluation index a1The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
the area A of the natural protection area invaded by the auxiliary facilities12The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
the primary protected object directly loses A31The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
the area of soil erosion32The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
the disruption of the habitat affects intensity A33The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
the indirect loss of protected object A42The calculation formula of the membership degrees of the grade I, the grade II and the grade III is as follows,
7. the ecological impact assessment system for road construction projects in natural reserve according to claim 4, characterized in that: the indirect loss of protected object A42The protected object in (1) refers to a specific protected object species in a natural protected area, a rare or endangered species or a main constituent species in a representative ecosystem.
8. The ecological impact assessment system for road construction projects in natural reserve according to claims 1 to 7, characterized in that: also comprises an ecological influence index test, the calculation formula is as follows,
the EI is an ecological influence index of a road construction project;
Diis AiI is an integer and i is greater than or equal to 1 and less than or equal to 4;
wiis AiI is an integer and is more than or equal to 1 and less than or equal to 4;
when the EI is less than or equal to 50, the ecological influence level of the road construction project is green; when EI is more than 50 and less than 80, the ecological influence level of the road construction project is yellow, and when EI is more than or equal to 80, the ecological influence level of the road construction project is red.
9. The ecological impact assessment system for road construction projects in natural reserve according to claim 8, characterized in that: said DiIncluding a road attribute index D1Sensitivity index D2Direct influence index D3And an indirect influence index D4;
The road attribute index D1The formula for calculating (a) is as follows,
wherein D is1Is a road attribute index;
d1the influence strength of the road length in the natural protection area is obtained;
d2the influence strength of the area of the road and the accessory facilities thereof;
x is the road grade influence strength; grade A of the current road13When the classified influence grade is grade III, x is 100; grade A of the current road13When the classified influence grade is grade II, x is 61.8; grade A of the current road13When the classified influence grade is grade I, x is 38.2;
the sensitivity index D2The formula for calculating (a) is as follows,
wherein D is2Is a sensitivity index;
y is the influence strength of the related functional regions; functional area A crossed by road construction project21When the classification level of (1) is level III, y is 100; functional area A crossed by road construction project21When the grade of (1) is grade II, y is 38.2; functional area A crossed by road construction project21When the classification level of (1) is level I, y is 0;
z is the impact intensity related to endangered species; when concerning distribution of endangered species A22When the classification level of (1) is level III, z is 100; when concerning distribution of endangered species A22When the grade is II, z is 61.8; when concerning distribution of endangered species A22When the grade is I, z is 38.2;
the direct influence index D3The formula for calculating (a) is as follows,
wherein D is3Is a direct impact index;
d3directly impact strength for the primary protected object;
d4setting the percentage of the water and soil loss area caused by road construction in the range of 30m at both sides of the road in the natural protection area to the total area of 30m buffer areas at both sides of the road;
d5to an environmental disruption index;
the indirect influence index D4The formula for calculating (a) is as follows,
wherein D is4Is an indirect influence index;
n is the foreign invasive species impact strength; when foreign invaded species status A41When the grade of (1) is grade III, n is 100; when foreign invaded species status A41When the grade of (1) is grade II, n is 61.8; when foreign invaded species status A41When the classification level of (1) is level I, n is 0;
d6the intensity is indirectly influenced for the main protected object.
10. The ecological impact assessment system for road construction projects in natural reserve according to claim 8, characterized in that: said wiAnalysis according to hierarchyLaw and expert consulting law to derive, w1=0.2,w2=0.3,w3=0.25,w4=0.25。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010455025.6A CN111652487B (en) | 2020-05-26 | 2020-05-26 | Ecological influence evaluation system for road construction project in natural conservation area |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010455025.6A CN111652487B (en) | 2020-05-26 | 2020-05-26 | Ecological influence evaluation system for road construction project in natural conservation area |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111652487A true CN111652487A (en) | 2020-09-11 |
CN111652487B CN111652487B (en) | 2021-07-13 |
Family
ID=72348605
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010455025.6A Active CN111652487B (en) | 2020-05-26 | 2020-05-26 | Ecological influence evaluation system for road construction project in natural conservation area |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111652487B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2362186C1 (en) * | 2007-12-27 | 2009-07-20 | Государственное унитарное предприятие города Москвы-объединенный эколого-технологический и научно-исследовательский центр по обезвреживанию РАО и охране окружающей среды (ГУП МосНПО "Радон") | Method for radio-environmental monitoring of industrialised region |
CN102136103A (en) * | 2011-03-14 | 2011-07-27 | 天津农学院 | Natural preservation area management effectiveness quantitative evaluation method and system |
CN102622651A (en) * | 2012-02-21 | 2012-08-01 | 天津市市政工程设计研究院 | Method for building and improving green ecological road |
CN105469196A (en) * | 2015-11-18 | 2016-04-06 | 山东科技大学 | Comprehensive evaluation method and comprehensive evaluation system for evaluating mine construction project process |
CN106446586A (en) * | 2016-10-21 | 2017-02-22 | 重庆大学 | River health evaluation method based on natural and social influence |
CN110245833A (en) * | 2019-05-15 | 2019-09-17 | 生态环境部卫星环境应用中心 | A kind of Development in Nature Reserve construction project Ecological Evaluation method and its system |
-
2020
- 2020-05-26 CN CN202010455025.6A patent/CN111652487B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2362186C1 (en) * | 2007-12-27 | 2009-07-20 | Государственное унитарное предприятие города Москвы-объединенный эколого-технологический и научно-исследовательский центр по обезвреживанию РАО и охране окружающей среды (ГУП МосНПО "Радон") | Method for radio-environmental monitoring of industrialised region |
CN102136103A (en) * | 2011-03-14 | 2011-07-27 | 天津农学院 | Natural preservation area management effectiveness quantitative evaluation method and system |
CN102622651A (en) * | 2012-02-21 | 2012-08-01 | 天津市市政工程设计研究院 | Method for building and improving green ecological road |
CN105469196A (en) * | 2015-11-18 | 2016-04-06 | 山东科技大学 | Comprehensive evaluation method and comprehensive evaluation system for evaluating mine construction project process |
CN106446586A (en) * | 2016-10-21 | 2017-02-22 | 重庆大学 | River health evaluation method based on natural and social influence |
CN110245833A (en) * | 2019-05-15 | 2019-09-17 | 生态环境部卫星环境应用中心 | A kind of Development in Nature Reserve construction project Ecological Evaluation method and its system |
Non-Patent Citations (4)
Title |
---|
YING LI等: "Study on Spatial Variability of Soil Infiltration Coefficient in the Plain of the Yili River Valley in Xinjiang", 《METEOROLOGICAL AND ENVIRONMENTAL RESEARCH》 * |
宗跃光 等: "道路生态学研究进展", 《生态学报》 * |
易浪等: "不同等级道路网络对常宁市景观格局生态效应的影响", 《湖南林业科技》 * |
李月辉 等: ""道路生态研究进展"", 《应用生态学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111652487B (en) | 2021-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Aburas et al. | Land suitability analysis of urban growth in Seremban Malaysia, using GIS based analytical hierarchy process | |
Li et al. | Early warning of illegal development for protected areas by integrating cellular automata with neural networks | |
Anelli et al. | Urban resilience against natural disasters: Mapping the risk with an innovative indicators-based assessment approach | |
CN113487123B (en) | Hyperspectral monitoring and GIS (geographic information System) coupled torrential flood disaster dynamic risk assessment method | |
CN104899437A (en) | Early-warning method for heavy-rainfall type landslide hazard | |
CN102289734A (en) | Dynamic risk evaluation and management method for foundation pit engineering | |
Zhang et al. | How road network transformation may be associated with reduced carbon emissions: An exploratory analysis of 19 major Chinese cities | |
Azari et al. | A GIS-based approach for accident hotspots mapping in mountain roads using seasonal and geometric indicators | |
CN114550445A (en) | Urban area traffic safety state evaluation method and device | |
Pusuluri et al. | Road crash zone identification and remedial measures using GIS | |
CN112508336B (en) | Space and environmental efficiency correlation measurement method based on structural equation model | |
CN117522151A (en) | Coal mine enterprise security risk assessment method and system | |
CN111652487B (en) | Ecological influence evaluation system for road construction project in natural conservation area | |
Rossouw | A review of methods and generic criteria for determining impact significance | |
Aggarwal et al. | Development of pavement management system for Indian national highway network | |
CN115438547A (en) | Overall evaluation method and system based on pavement service state | |
Dunne et al. | A large scale method for extracting geographical features on bus routes from OpenStreetMap and assessment of their impact on bus speed and reliability | |
Suraji et al. | Correlation of Expansive Soil and Road Pavement Conditions Using Data Mining from GIS Portal | |
Li et al. | Study on road damage assessment based on RS and GIS | |
Haghparast et al. | Comprehensive Environmental Monitoring based on Stations of Environmental Pollutants (Air, Water and Soil) in Tehran | |
Aquil et al. | Comparative Study on Spatial Clustering Methods for Identifying Traffic Accident Hotspots | |
Tsapakis et al. | Statewide Implementation of Innovative Safety Analysis Tools in Identifying Highway Safety Improvement Projects: Technical Report | |
Anabestani et al. | Analysis of Key Proponents Affecting Land Use Change of Pri-urban Settlements in Mashhad metropolitan with approach future study | |
Yu et al. | Simulation-based safety analysis of diversion and confluence area in urban underground loop | |
Ali et al. | Identification of accident black-spots in 18 Michigan freeways using GIS |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |