CN116562698A - Method, system, equipment and storage medium for evaluating coordinated development of flood storage area - Google Patents

Method, system, equipment and storage medium for evaluating coordinated development of flood storage area Download PDF

Info

Publication number
CN116562698A
CN116562698A CN202310532850.5A CN202310532850A CN116562698A CN 116562698 A CN116562698 A CN 116562698A CN 202310532850 A CN202310532850 A CN 202310532850A CN 116562698 A CN116562698 A CN 116562698A
Authority
CN
China
Prior art keywords
evaluation
water
flood
index
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202310532850.5A
Other languages
Chinese (zh)
Inventor
王宇晖
姬金雪
付政
张俊
崔申申
韩戈
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Donghua University
Original Assignee
Donghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Donghua University filed Critical Donghua University
Priority to CN202310532850.5A priority Critical patent/CN116562698A/en
Publication of CN116562698A publication Critical patent/CN116562698A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Economics (AREA)
  • Tourism & Hospitality (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Theoretical Computer Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a coordinated development evaluation method for a diapause area, which is characterized in that an evaluation index set of the diapause area, which relates to flood, economy and ecology, is determined based on a fuzzy comprehensive evaluation method, and an evaluation value of the evaluation index set is obtained through calculation according to historical data of the evaluation index set and a pre-constructed value accounting model; on the basis of the evaluation index set, an AHP algorithm is adopted to construct a coordinated development evaluation system of the flood storage area, an AHP-entropy weight method is adopted to calculate and obtain weights of all elements in the coordinated development evaluation system, a final decision set is obtained according to the weights, finally, based on a maximum membership principle, a comment corresponding to the maximum value in the decision set is selected as an evaluation result, and an objective and more accurate analysis result is provided for sustainable development of the flood storage area.

Description

Method, system, equipment and storage medium for evaluating coordinated development of flood storage area
Technical Field
The invention relates to the field of data analysis, in particular to a method, a system, equipment and a storage medium for evaluating coordinated development of a flood storage area.
Background
The flood storage area refers to a low-lying area and a lake which are used for temporarily storing flood or separating flood peak outside the back surface of a river levee including a flood diversion port, and comprises four components of a flood delivery area, a flood diversion area, a flood storage area and a flood stagnation area. Generally speaking, the flood storage areas in China are mostly lakes and flood-flooded wetlands, the population is relatively sparse, and the flood storage areas play an important role in regulating flood, purifying water quality and protecting ecology.
For a long time, the economic development level in the flood storage area has been limited due to the pressure of the flood storage water; meanwhile, along with population growth and economic development in the area, human activities have serious influence on the ecological environment in the flood storage area, and how to evaluate the development of the flood storage area according to investigation data, coordinate the relations among flood control, economy and ecology according to evaluation results, and exert various functions and benefits of the flood storage area to the maximum extent, so that the flood storage area is on a sustainable development road, and the method is a problem facing to the present.
Disclosure of Invention
In view of the above-mentioned drawbacks of the prior art, the present invention aims to solve the technical problems of how to evaluate the development of the flood storage area according to the investigation data, coordinate the relationship among flood control, economy and ecology according to the evaluation result, and exert various functions and benefits of the flood storage area to the maximum extent, so that the flood storage area can go on the road of sustainable development, which is a problem facing to the present.
In order to achieve the above object, the present invention provides a method for evaluating coordinated development of a flood storage area, comprising:
determining an evaluation index set and a comment set of the flood storage area based on a fuzzy comprehensive evaluation method; wherein the evaluation index set comprises flood index, economic index and ecological index;
Determining a membership function between the evaluation index set and the evaluation set;
according to the historical data of the evaluation index set, calculating to obtain an evaluation value of the evaluation index set;
establishing an evaluation matrix according to the evaluation value and the membership function;
constructing a coordinated development evaluation system of a flood storage area according to an AHP algorithm, wherein the coordinated development evaluation system comprises three layers, a first layer is a target to be realized, a second layer and a third layer comprise at least one element, and the elements of the third layer are the evaluation index set;
according to the evaluation value, calculating to obtain the weight of each element by adopting an AHP-entropy weight method;
according to the evaluation matrix and the weight, calculating to obtain a synthetic matrix;
and selecting a comment corresponding to the maximum value in the synthesis matrix as an evaluation result.
In a preferred embodiment of the present invention, the calculating, according to the historical data of the evaluation index set, an evaluation value of the evaluation index set includes:
according to the historical data of the flood index, taking the numerical value meeting a first preset condition in a preset time as an evaluation value of the flood index;
according to the historical data of the economic index, taking the numerical value meeting a second preset condition in the preset time as an evaluation value of the economic index;
And selecting a numerical value meeting a third preset condition in a preset time according to the historical data of the ecological index, and inputting the selected numerical value into a pre-constructed value accounting model to obtain an evaluation value of the ecological index.
In another preferred embodiment of the present invention, the ecological index includes water conservation, climate control, air purification and water purification;
the value accounting model comprises a water source conservation sub-model, a climate adjustment sub-model, an air purification sub-model and a water quality purification sub-model;
the water conservation sub-model includes:
Q l =Q b +Q x +Q f +Q j -Q t -Q s -Q z
V w =Q l ×P s ×C s +Q l ×P y ×C y +Q l ×P g ×C g
P s +P y +P g ≤100%
wherein Q is l For accumulating water resource quantity, Q of flood area b Is the surface water quantity, Q x Make-up quantity, Q for groundwater penetration f Is of flood quantity, Q j Is precipitation amount, Q t For flood removal, Q s To maintain the amount of water, Q, consumed by the ecosystem z Vapor deposition for the area; v (V) w Evaluation value for conservation of water source, P s Percentage of water needed to maintain the ecosystem and the amount of water available, C s Market price, P, of water resource trade for maintaining water consumed by ecosystem y The water consumption is the percentage of the applicable water consumption in the cultivation water consumption occupying area, C y Trade market price, P, for water resources of aquaculture water g For the percentage of available water resources and C in the irrigation water occupying area g Trading market prices for water resources of irrigation water;
the climate conditioner sub-model comprises:
Q z =(E f1 +E f2 )×A
E f1 =E s +E w +E p
E f2 =S×T
V q =E tt ×P e
E tt =E z ×q×10 3 /3600+E z ×y
wherein Q is z For total evaporation of stagnant zone, E f1 For the evapotranspiration of the area after flood withdrawal, E f2 The evaporation capacity of the area during flood diversion, and the area A is the area of the area; e (E) s Soil vapor emission amount E for region w Vapor emission amount of water surface of region E p The method is characterized in that the method comprises the steps of transpiration of plants in a region, S is the evaporation rate of water surface in the region when flood storage areas are used for flood diversion, and T is the evaporation time; v (V) q For climate control evaluation value, E tt Is an ecological systemEnergy, P, consumed by evaporation e For local electricity price, E z Is the evaporation capacity, q is the latent heat of volatilization, y is 1m 3 Power consumption of water to steam;
the air purification submodel comprises:
wherein Q is w Is the total emission quantity, Q of the atmospheric pollutants i Is the discharge amount of the i-type atmospheric pollutants and Q y Air purifying capacity, Q of regional ecological system ij Purification amount of the ith treatment unit for the jth class of atmospheric pollutants, A i Is the area of the ith processing unit; v (V) i Value, Q for purifying class i atmospheric environment i Purifying amount of class i atmospheric pollutants, C i The treatment cost and V of the i-th type atmospheric pollutants k An evaluation value for air purification;
the water quality purifying sub-model comprises:
Wherein Q is i To purify the i-th pollutant, P i For the treatment capacity of each treatment unit of the flood storage area to the ith pollutant, S i For the area of the processing unit, T is the processing time, Q c Is the total purification amount of pollutants, V i For purifying the water quality of the i-th pollutant and C i Unit treatment cost for the ith pollutant, V z An evaluation value for water quality purification.
In another preferred embodiment of the present invention, the ecological index further comprises soil conservation, carbon fixation and species conservation; the value accounting model also comprises a soil conservation sub-model, a carbon fixation sub-model and a species conservation sub-model;
the soil conservation sub-model comprises:
V s =V f +V w
V f =(Q d -Q l )×P 1
V w =(Q d -Q l )×P 2
wherein V is s Maintaining an evaluation value, V, for soil f Is the nutritive value of sediment, V w Can be converted into the value of soil, Q d Amount of deposit, Q l The soil quantity for loss, P1 is the price of fertilizer which needs to be manually input, P 2 The price of the soil which needs to be manually input is set;
the carbon fixation sub-model comprises:
V c =Q c ×C c
wherein V is c For the evaluation value of carbon sequestration, Q c Is the total carbon content, C c For carbon trade price, S i Ability to fix carbon for individual land types, A i Utilizing the area of land for each land type;
The species conservation sub-model comprises:
V s =V b +V e
V e =Q e ×P e
wherein V is s Conservation evaluation value, V for total species b Is of value for species conservation, V e For the cost of the artificially added species, S i For the number of various animals and plants, P i Cost, Q of artificial conservation for each species of conservation e P, an amount of additional valuable species e Is the value of the species requiring manual input.
In another preferred embodiment of the present invention, the calculating, according to the evaluation value, the weight of each element by using an AHP-entropy weight method includes:
calculating subjective weight of each element by adopting an AHP algorithm;
calculating objective weights of the elements based on the evaluation values by adopting an entropy weight method;
and determining a combined weight according to the subjective weight and the objective weight, and taking the calculated combined weight as the weight of each element.
In another preferred embodiment of the present invention, the ecological assessment system includes a target layer, a criterion layer, and an index layer;
the subjective weight of each element is calculated by AHP, which comprises the following steps:
respectively constructing a first judgment matrix of a target layer alignment rule layer and a second judgment matrix of a criterion layer pair index layer according to the relative importance degree among the elements;
Respectively carrying out hierarchical single sequencing and consistency check on the first judgment matrix, and respectively carrying out hierarchical single sequencing and consistency check on the second judgment matrix;
performing hierarchical total sequencing and consistency verification on the weights of the target layers passing the consistency verification and the target layers passing the consistency verification, and performing hierarchical total sequencing and consistency verification on the weights of the index layers by the criterion layers passing the consistency verification;
and taking the weight passing the consistency check as the subjective weight.
In another preferred embodiment of the present invention, the calculating the objective weight of each element based on the evaluation value using the entropy weight method includes:
constructing a third judgment matrix of n evaluation indexes of m samples based on the evaluation values;
sequentially carrying out standardization processing and normalization processing on the third judgment matrix to obtain a normalized third judgment matrix;
calculating entropy of each evaluation index according to the normalized third judgment matrix;
and taking the calculated entropy weight of each element as the objective weight according to the entropy of each evaluation index.
In a second aspect, the present application also provides a coordinated development assessment system for a diapause area. The system comprises:
The index determining module is used for determining an evaluation index set and a comment set of the flood storage area based on a fuzzy comprehensive evaluation method; wherein the evaluation index set comprises flood index, economic index and ecological index;
the computing module is used for determining a membership function between the evaluation index set and the evaluation set; according to the historical data of the evaluation index set, calculating to obtain an evaluation value of the evaluation index set; establishing an evaluation matrix according to the evaluation value and the membership function;
the evaluation system construction module is used for constructing a coordinated development evaluation system of the flood storage area according to an AHP algorithm, wherein the coordinated development evaluation system comprises three layers, a first layer is a target to be realized, a second layer and a third layer comprise at least one element, and the elements of the third layer are the evaluation index set;
the calculation module is further used for calculating the weight of each element by adopting an AHP-entropy weight method according to the evaluation value;
the evaluation module is used for calculating to obtain a synthesis matrix according to the evaluation matrix and the weight; and selecting a comment corresponding to the maximum value in the synthesis matrix as an evaluation result.
In a third aspect, the present application also provides a computer device. The computer device comprises a memory storing a computer program and a processor implementing the method steps of any of the first aspects when the processor executes the computer program.
In a fourth aspect, the present application also provides a computer-readable storage medium. The computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements the method steps of any of the first aspects.
The method, the system, the equipment and the storage medium for evaluating the coordinated development of the flood storage area have at least the following advantages:
1. the evaluation indexes comprise flood indexes, economic indexes and ecological indexes, and fully consider dynamic ecological characteristics of flood advance and retreat in a flood storage area and economic characteristics caused by the dynamic ecological characteristics, so that the evaluation indexes are comprehensive and balanced;
2. the weight of each element in the coordinated development evaluation system is calculated by adopting an AHP-entropy weight method, so that the calculated weight distribution is more reasonable under the conditions of multiple indexes and multiple factors;
3. the fuzzy comprehensive evaluation method can divide the development change interval of the problem to be evaluated into a plurality of stages and states, then analyze various aspects of the problem and the current state in deep detail, analyze the relative degree of the overall state belonging to the stage class of each stage, and enable the description of things to be more objective and the analysis result to be more accurate.
The conception, specific structure, and technical effects of the present invention will be further described with reference to the accompanying drawings to fully understand the objects, features, and effects of the present invention.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention.
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for the description of the embodiments will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of the method for evaluating the coordinated development of a flood storage area according to the invention;
FIG. 2 is a flow chart of the invention for calculating the weights of elements;
FIG. 3 is a schematic diagram of the coordinated development assessment system for a flood storage area of the present invention;
FIG. 4 is an internal block diagram of a computer device in one embodiment.
Detailed Description
Other advantages and effects of the present invention will become apparent to those skilled in the art from the following disclosure, which describes the embodiments of the present invention with reference to specific examples. The invention may be practiced or carried out in other embodiments that depart from the specific details, and the details of the present description may be modified or varied from the spirit and scope of the present invention. It should be noted that the following embodiments and features in the embodiments may be combined with each other without conflict.
It should be noted that the illustrations provided in the following embodiments merely illustrate the basic concept of the present invention by way of illustration, and only the components related to the present invention are shown in the illustrations, not according to the number, shape and size of the components in actual implementation, and the form, number and proportion of each component in actual implementation may be arbitrarily changed, and the layout of the components may be more complex.
Some exemplary embodiments of the invention have been described for illustrative purposes, it being understood that the invention may be practiced otherwise than as specifically shown in the accompanying drawings.
Referring to fig. 1, fig. 1 is a flow chart of a method for estimating the development of a flood storage area according to the present embodiment, which specifically includes the following steps:
step S102, determining an evaluation index set and a comment set of the flood storage area based on a fuzzy comprehensive evaluation method; wherein the evaluation index set comprises flood index, economic index and ecological index.
Step S104, determining a membership function of the evaluation index set.
Step S106, according to the historical data of the evaluation index set, an evaluation value of the evaluation index set is obtained through calculation.
And S108, establishing an evaluation matrix according to the evaluation value and the membership function.
Step S110, constructing a coordinated development evaluation system of the flood storage area according to an AHP algorithm, wherein the coordinated development evaluation system comprises three layers, the first layer is a target to be realized, the second layer and the third layer comprise at least one element, and the third layer element is an evaluation index set.
Step S112, calculating to obtain the weight of each element by adopting an AHP-entropy weight method according to the evaluation value;
step S114, calculating to obtain a synthetic matrix according to the evaluation matrix and the weight; selecting the comment corresponding to the maximum value in the synthesis matrix as an evaluation result.
Specifically, the fuzzy comprehensive evaluation method is a comprehensive evaluation method based on fuzzy mathematics. According to membership theory of fuzzy mathematics, the comprehensive evaluation method converts qualitative evaluation into quantitative evaluation, namely, the fuzzy mathematics are used for carrying out overall evaluation on things or objects limited by various factors. The method has the characteristics of clear result and strong systematicness, can better solve the problems of ambiguity and difficult quantization, and is suitable for solving various nondeterminacy problems.
In particular, in this embodiment, in order to realize the requirement of sustainable development in the flood storage area, the present application establishes an evaluation index set of an object to be evaluated from three aspects of flood control, economy and ecology, and the selected evaluation index set needs to satisfy the following requirements: firstly, objectively describing and representing the current situation of the development of the flood storage area; secondly, the change rate of each aspect of the development of the flood storage area can be dynamically represented, and the change trend can be reflected through the change rate; thirdly, the coordinated scheduling of various aspects of the development of the flood storage area can be embodied.
Optionally, the evaluation index set U in the present embodiment includes a flood index U 1 Economic index u 2 And an ecological index u 3 I.e. u= { U 1 ,u 2 ,u 3 }。
Wherein the flood control aspect considers the influence degree of the application of the flood storage area on the area and the disaster resistance of the area, thereby determining the flood control subsystem u 1 The corresponding flood indicator includes an activation frequency u 11 Maximum submerged depth u 12 Flood protection area u 13 Flood storage amount u 14 Emergency transfer setting population u 15 Qualification rate u of dike 16 Quantity u of sluice 17 Any one or more of the following.
The economic aspect mainly considers the current development degree of the flood storage area, describes the development level of the economy in the area by input and output, mainly relates to factors influencing the economy, such as population, resources and society, and the population factors consider the influence of population quality and quantity on the development of the flood storage area, the resource factors consider the factors, such as the existing resources, the average resource possession and development potential, in the flood storage area, and the social factors consider the factors, such as the social stability degree and the normal living level of people. The economical subsystem u thus determined 2 The corresponding economic index comprises natural population growth rate u 21 Rate u of illiterate population 22 Area u of arable land 23 Water resource u available for everyone 24 Average person GDP u 25 GDP growth Rate u 26 The average resident person can control income u 27 Area u of average residence 28 Urban registration loss rate u 29 Any one or more of the following.
The ecological aspect mainly considers the value brought by the functions of ecological products, ecological services and the like, and the ecological subsystem u is determined by the ecological aspects 3 The corresponding ecological index comprises water source conservation u 31 Climate control u 32 U for purifying air 33 U for purifying water 34 Water and soil conservation u 35 Carbon fixation u 36 Species conservation u 37 Any one or more of the following.
Further, a comment set v= { V1, V2, …, vm } is established, which is a set of comments made by the object to be evaluated, for representing the degree of merit of the evaluation factor. If the flood storage area can be continuously developed, five grades of V= { good, better, general, worse and bad } are adopted.
Further, the membership is a mathematical scale that characterizes the magnitude of uncertainty in membership of the fuzzy set. In fuzzy comprehensive evaluation, selecting or constructing a proper fuzzy membership function is a key step for obtaining a reasonable evaluation result. Because the meaning and the measurement unit of each index are different, an appropriate membership function must be constructed to perform standardization processing on index data. The membership function of each index is established by soliciting expert opinion and the existing research results and combining the existing average level of part of the indexes.
In this embodiment, the membership degree of the evaluation index set to the comment set characterizes the degree that each index in the evaluation index set belongs to the comment set. Generally, determining the membership function mainly includes three calculation methods: the fuzzy statistics method, the assignment method and the borrowing of the existing objective scale can select a proper calculation method according to the needs in practical application China, and the embodiment does not limit the determination mode of the membership function.
By the above way, the corresponding membership function of each index is obtained, such as the economic subsystem u 2 Natural population growth rate u in (a) 21 And the illiterate rate u of population 22 The established membership function:
substituting the evaluation value of each index in the evaluation index set into the corresponding membership function to obtain an evaluation matrix R for representing the index relation in different functional subsystems, wherein the expression is as follows:
wherein m is the number of flood storage areas; n is the number of selected indexes. For example, in a specific m embodiments, flood control subsystem u 1 The corresponding evaluation matrix R1 is a matrix of 7×m, and the economic subsystem u 2 The corresponding evaluation matrix R2 is a matrix of 9×m, and the ecological subsystem u 3 The corresponding evaluation matrix R3 is a 7×m matrix.
Optionally, calculating, according to the historical data of the evaluation index set, an evaluation value of the evaluation index set includes:
And taking the numerical value meeting the first preset condition in the preset time as an evaluation value of the flood index according to the historical data of the flood index.
And taking the numerical value meeting the second preset condition in the preset time as an evaluation value of the economic index according to the historical data of the economic index.
According to the historical data of the ecological indexes, selecting the numerical value meeting a third preset condition in the preset time, and inputting the selected numerical value into a pre-constructed value accounting model to obtain the evaluation value of each ecological index.
The first preset condition, the second preset condition and the third preset condition are determined according to each index parameter of the flood index, the economic index and the ecological index. In addition, the predicted values of flood indicators and economic indicators can be directly or indirectly obtained through historical dataObtaining, for example, the activation frequency u in the flood indicator 11 And maximum submerged depth u 12 The starting frequency can be obtained by calling the starting times in the preset time; and obtaining the maximum submerged water depth by calling and comparing the submerged water depth data within the preset time.
The ecological index is different from the ecological value accounting of the general area to a certain extent because the flood storage area relates to the dynamic process of flood advance and retreat, the ecological characteristics of the flood storage area in the normal state are considered, and the state of the whole ecological system in the flood storage state and the influence change of the flood advance and retreat process on the area are considered. The principle of accounting should follow science, practicality, systematicness and openness, and more should follow the principle of specificity and dynamics, and accounting the ecological service functional value can bring ecological benefits into the economic and social development evaluation system of the flood storage area, and provide basis for sustainable development of the whole area.
Further, when the ecological index includes water conservation, climate adjustment, air purification, and water purification, the value accounting model includes a water conservation sub-model, a climate adjustment sub-model, an air purification sub-model, and a water purification sub-model.
The water source conservation service is the function of an ecological system for intercepting stagnant flood, enhancing soil infiltration and accumulation, conserving soil moisture, regulating storm runoff and supplementing groundwater, and increasing available water resources. The area with large water conservation amount not only meets the water source requirements of production and life in the accounting area, but also continuously provides water resources outside the area. The conservation value of the water source is mainly represented by the economic value of water storage and water retention. The water conservation value of the flood storage area is mainly reflected in the utilization of regional water resources, and the regional water resources are mainly used for maintaining the self-demand, cultivation and irrigation of the ecological system due to the special development characteristics and requirements of the flood storage area.
The conception of the water source conservation sub-model is as follows: the water conservation evaluation value is equal to the product of the water conservation amount in the accounting area and the water resource trading market price. The water conservation amount of the ecological system is calculated by physical quantity, and the water conservation amount is calculated by adopting a water balance equation preferentially. The water balance equation means that the water keeps the quality conservation in the ecological system in a certain time and space, namely, the water conservation amount of the ecological system is the difference value between the input and flood separation amount of precipitation and the flood removal amount and the water consumption amount of the ecological system.
The water source conservation sub-model constructed by the method comprises the following steps:
Q l =Q b +Q x +Q f +Q j -Q t -Q s -Q z
V w =Q l ×P s ×C s +Q l ×P y ×C y +Q l ×P g ×C g
P s +P y +P g ≤100%
wherein Q is l For accumulating water resource quantity, Q of flood area b Is the surface water quantity, Q x Make-up quantity, Q for groundwater penetration f Is of flood quantity, Q j Is precipitation amount, Q t For flood removal, Q s To maintain the amount of water, Q, consumed by the ecosystem z Vapor deposition for the area; v (V) w Evaluation value for conservation of water source, P s Percentage of water needed to maintain the ecosystem and the amount of water available, C s Market price, P, of water resource trade for maintaining water consumed by ecosystem y The water consumption is the percentage of the applicable water consumption in the cultivation water consumption occupying area, C y Trade market price, P, for water resources of aquaculture water g For the percentage of available water resources and C in the irrigation water occupying area g Market prices are traded for water resources of irrigation water.
The climate regulation service of the ecological system refers to the ecological function of the ecological system that the ecological system absorbs solar energy through vegetation transpiration and water surface evaporation, reduces air temperature, increases air humidity and improves the comfort level of the living environment. The energy consumed in the evaporation process of the ecological system is selected as an evaluation index of the climate regulating service of the ecological system.
The evaporation capacity of the flood storage area is mainly divided into two states, wherein the first state is a general state when no flood is stored, and the second state is a state when flood is separated. The evaporation in the normal state mainly comprises water surface evaporation, soil evaporation and plant emission of various water bodies (rivers, lakes, reservoirs and the like). Land evaporation is related to the structure of the soil, the moisture content and the situation of plant coverage. Wherein the water surface evaporation is related to the phase state of water (water or ice), temperature, water surface size and the condition of impurities contained in water; soil evaporation is mainly affected by factors such as soil water content, groundwater burial depth, soil structure, soil surface characteristics and topography, and is mainly divided into three evaporation stages: a capillary operation stage, a film operation stage, and a diffusion operation stage; the transpiration of plants is not only influenced by external environmental conditions, but also by the regulation and control of the plants themselves, mainly related to the type of plants in the area and the vegetation cover. The evaporation capacity of the flood storage area in flood diversion is mainly related to the evaporation rate and evaporation time.
The climate control sub-model is conceived in that: the evaluation value of the climate regulation of the ecological system is the product of the total energy consumed by the temperature or humidity regulation of the ecological system and the local electricity price, and the cost replacement method, namely the power consumption required by the manual temperature and humidity regulation is mainly applied to account the value of the transpiration regulation temperature or humidity of the ecological system and the value of the evaporation regulation temperature or humidity of the water surface. The energy consumption of the ecosystem for regulating the temperature or the humidity is obtained by accounting the real object quantity. The climate regulating service accounting can be performed by actually measuring the temperature difference between the inside and outside of the ecological system, the solar energy consumed by the ecological system and the total evapotranspiration of the ecological system, wherein the actual measuring method is preferentially selected, and then the solar energy consumed by the ecological system or the total evapotranspiration of the ecological system is selected according to the data availability to perform the accounting. Electricity prices are obtained from related documents or power departments released by the development and reform committee in the accounting place, generally referenced to industrial electricity prices.
The climate conditioner sub-model thus constructed comprises:
Q z =(E f1 +E f2 )×A
E f1 =E s +E w +E p
E f2 =S×T
V q =E tt ×P e
E tt =E z ×q×10 3 /3600+E z ×y
wherein Q is z For total evaporation of stagnant zone, E f1 For the evapotranspiration of the area after flood withdrawal, E f2 The evaporation capacity of the area during flood diversion, and the area A is the area of the area; e (E) s Soil vapor emission amount E for region w Vapor emission amount of water surface of region E p The method is characterized in that the method comprises the steps of transpiration of plants in a region, S is the evaporation rate of water surface in the region when flood storage areas are used for flood diversion, and T is the evaporation time; v (V) q For climate control evaluation value, E tt Energy, P, consumed for evaporation of ecosystems e For local electricity price, E z Is the evaporation capacity, q is the latent heat of volatilization, y is 1m 3 The water is converted into electricity consumption of steam.
The air purifying function refers to the function of an ecological system for absorbing, filtering, blocking and decomposing atmospheric pollutants (such as sulfur dioxide, nitrogen oxides, particulate matters and the like), purifying the air pollutants and improving the atmospheric environment. The air purifying function is mainly embodied in the aspects of purifying pollutants and blocking particulate matters.
The ecological system air purification evaluation value refers to the ecological effect generated by the ecological system for absorbing, filtering, blocking and decomposing and reducing the atmospheric pollutants (such as sulfur dioxide, nitrogen oxides, particulate matters and the like) so as to improve the atmospheric environment. The air purification value of the ecological system is calculated by adopting an alternative cost method, namely, the cost of atmospheric pollutants is treated by industry.
The method for calculating the purification value of sulfur dioxide, nitrogen oxides and particulate matters comprises the following steps: the three pollutant air purification physical quantities of sulfur dioxide, nitrogen oxide and particulate matter are respectively multiplied by the unit cost of sulfur dioxide, nitrogen oxide and particulate matter treatment to calculate the air purification value. The value of the ecological system for purifying the atmospheric environment is the sum of the purifying amount of various atmospheric pollutants and the treatment cost of the corresponding atmospheric pollutants.
Air purification function accounting different methods are chosen depending on whether the contaminant concentration exceeds the ambient air function area quality criteria. If the pollutant concentration does not exceed the quality standard of the environmental air functional area, calculating pollutant purifying amount by adopting pollutant emission amount, namely, a Q_w calculation method; if the pollutant concentration exceeds the quality standard of the environmental air functional area, the self-cleaning capacity of the ecological system is adopted to calculate the pollutant purifying amount, namely a Q_y calculating method. The unit treatment cost adopts a pollution discharge fee collection standard which is printed in the accounting place, has no local standard, and can refer to the tax standard in the pollution discharge fee collection standard and calculation method published by the national development and reform committee or the tax standard in the environmental protection tax law of the people's republic of China.
The air purification submodel constructed by the method comprises the following steps:
wherein Q is w Is the total emission quantity, Q of the atmospheric pollutants i Is the discharge amount of the i-type atmospheric pollutants and Q y Air purifying capacity, Q of regional ecological system ij Purification amount of the ith treatment unit for the jth class of atmospheric pollutants, A i Is the area of the ith processing unit; v (V) i Value, Q for purifying class i atmospheric environment i Purifying amount of class i atmospheric pollution, C i The treatment cost and V of the i-th type atmospheric pollutants k An evaluation value for air purification.
The water quality purifying function refers to the functions of adsorbing, degrading and converting water pollutants by a wetland ecosystem in water areas such as lakes, rivers, marshes and the like and purifying water environment. The water quality purifying service evaluation value accounting mainly uses monitoring data, and selects proper indexes to quantitatively account according to pollutant composition and concentration change in a water ecological system. Common indicators include ammonia nitrogen, COD, total nitrogen, total phosphorus, partial heavy metals and the like.
Due to the influence of factors such as soil structure, vegetation coverage, terrain, flood residence time and the like in the flood storage area, the flood storage area can be divided into various processing units, and the pollutant removal capacity of each processing unit is different. The water quality purifying function calculates and selects different methods according to whether the pollutant concentration exceeds the water area environment function and the protection target, and if the pollutant concentration does not exceed the water area environment function standard limit value, the pollutant emission quantity is adopted to calculate the physical quantity; if the pollutant emission concentration exceeds the water area environment function standard limit value, the physical quantity is calculated according to the self-cleaning capability of the ecological system.
The water quality purification evaluation value is calculated by adopting a substitution cost method, and the water quality purification value of the ecological system is calculated by industrial treatment of water pollutant cost. Such as chemical oxygen demand and ammonia nitrogen purification value calculation method: the water purification value is calculated by multiplying the water purification physical quantity of two pollutants, namely the chemical oxygen demand and the ammonia nitrogen, by the unit chemical oxygen demand and the ammonia nitrogen treatment cost respectively. The value of the water quality purification of the ecological system is the sum of the purification amount of each treatment unit to various water pollutants and the unit treatment cost of the corresponding water pollutants.
The water quality purifying sub-model constructed by the method comprises the following steps:
wherein Q is i To purify the i-th pollutant, P i For the treatment capacity of each treatment unit of the flood storage area to the ith pollutant, S i For the area of the processing unit, T is the processing time, Q c Is the total purification amount of pollutants, V i For purifying the water quality of the i-th pollutant and C i Unit treatment cost for the ith pollutant, V z An evaluation value for water quality purification.
Optionally, when the ecological index further includes soil conservation, carbon fixation, and species conservation, the value accounting model further includes a soil conservation sub-model, a carbon fixation sub-model, a species conservation sub-model.
The soil conservation evaluation value is mainly expressed in the conservation of soil fertility, and the cost required by manually increasing the soil fertility, namely the value amount brought by flood, is required by applying the replacement cost method. For example, in the downstream region of yellow river, there is much saline-alkali soil along the coast, which is unfavorable for the growth of crops, while the sediment of yellow river is rich in many minerals, especially the sediment in flood season has fertilizer efficiency, and can be changed into a fertile farmland by covering the surface layer of saline-alkali soil by silting. The flood is a accident, and meanwhile, thick fertilizer is also brought, and the flood is generated in some places, and even the fertilizer is not needed for years. Meanwhile, the ecological system can reduce the non-point source pollution by keeping the soil and reducing the entry of nitrogen, phosphorus and other soil nutrients into the water body and the air. According to the soil holding quantity and the content of nitrogen and phosphorus in the soil, an alternative cost method is applied, namely cost accounting of pollutant treatment reduces the value of non-point source pollution.
The soil conservation sub-model thus constructed comprises:
V s =V f +V w
V f =(Q d -Q l )×P 1
V w =(Q d -Q l )×P 2
wherein V is s Maintaining an evaluation value, V, for soil f Is the nutritive value of sediment, V w Can be converted into the value of soil, Q d Amount of deposit, Q l For the amount of lost soil, P 1 To the price and P of the fertilizer which needs to be manually input 2 The cost of the soil which needs to be manually input is high.
The carbon fixing function of the ecological system refers to the function of the natural ecological system for absorbing carbon dioxide in the atmosphere to synthesize organic matters and fixing carbon in plants or soil. The function is beneficial to reducing the concentration of carbon dioxide in the atmosphere and slowing down the greenhouse effect. The carbon fixation function of the ecological system has important significance for reducing the emission reduction pressure. And selecting the fixed carbon dioxide as an evaluation index of the carbon fixation function of the ecological system.
The biological system fixed carbon evaluation value can adopt an alternative cost method (forestation cost method and industrial emission reduction cost) and a market value method (carbon trade price) to calculate the value of the ecological system fixed carbon. The carbon sequestration value of the ecosystem is the product of the total amount of carbon sequestration and the carbon price of the ecosystem. The carbon fixation capacity of ecosystems, such as forests, grasslands, wetlands, farmlands, and soil, varies from one land type to another. The carbon fixation amount is obtained by accounting the physical amount. The carbon sequestration cost for unit forestation, the emission reduction cost for industrial carbon and the market price for carbon trade refer to the relevant literature and the market price for carbon market is recommended.
The carbon stator model thus constructed includes:
V c =Q c ×C c
wherein V is c For the evaluation value of carbon sequestration, Q c Is the total carbon content, C c For carbon trade price, S i Ability to fix carbon for individual land types, A i The area of land is utilized for each land type.
Species diversity is the most prominent structural and functional unit of biodiversity, can provide necessary species and genetic resources for biological system succession and biological evolution, and is the basis for human survival and development. The conservation service refers to the role and value of the ecological system for providing living and reproduction places for rare animal and plant species, thereby playing a role in conservation. And calculating an evaluation value of the species conservation service by using the conservation cost per unit area, and multiplying the physical quantity by the species conservation cost to obtain the evaluation value. In addition, some fish and other species of interest in the area can be used as baits for some animals in the area, and the added value of the baits which are manually input, namely the added value which can be brought by wetland type flood storage areas, is usually required.
The thus constructed species conservation sub-model comprises:
V s =V b +V e
V e =Q e ×P e
wherein V is s Conservation evaluation value, V for total species b Is of value for species conservation, V e For the cost of the artificially added species, S i For the number of various animals and plants, P i Cost, Q of artificial conservation for each species of conservation e P, an amount of additional valuable species e Is the value of the species requiring manual input.
In the above embodiment, the evaluation values of the flood index and the economic index may be obtained through the history data; meanwhile, an evaluation value of an ecological index which can embody the dynamic characteristics of the diapause area and has higher accuracy can be obtained through a value accounting model which is built in advance, and a reliable basis is provided for the evaluation of the sustainable development of the diapause area.
Further, according to the demand of coordinated development in the flood storage area, constructing a coordinated development evaluation system of the flood storage area on the basis of an evaluation index set according to an AHP algorithm, wherein the coordinated development evaluation system comprises a target layer, a criterion layer and an index layer, wherein the target layer is the sustainable development of the flood storage area and characterizes the problem to be solved; elements in the criteria layer include flood control subsystem, economic subsystem and ecological subsystem, characterizing the strategies employed to solve the problem; the elements in the index layer are evaluation index sets capable of reflecting the conditions of all subsystems, and represent specific indexes or specific influencing factors selected for solving the problems.
It should be appreciated that the analytic hierarchy process (analytic hierarchy process, AHP), which primarily addresses the problem of evaluating classes. Such as which scheme is best chosen, which employee performs best, etc. The method is a relatively subjective evaluation method, and when the weight vector is obtained by weighting, the subjective factor is large in duty ratio. Therefore, on the basis of AHP, the method also adopts an entropy weight method to obtain objective weight, synthesizes the weight vector obtained by the subjective method with the weight vector obtained by the objective method to obtain a comprehensive weight vector, and performs subsequent operation.
Referring to fig. 2, optionally, according to the evaluation value, an AHP-entropy weighting method is adopted to calculate the weight of each element, including:
step S202, calculating subjective weight of each element by adopting AHP;
step S204, calculating objective weights of all elements based on the evaluation values by adopting an entropy weight method;
step S206, determining a combined weight according to the subjective weight and the objective weight, and taking the calculated combined weight as the weight of the evaluation index set.
Optionally, the subjective weight of each element is calculated by using the AHP, including:
respectively constructing a first judgment matrix of a target layer alignment rule layer and a second judgment matrix of a criterion layer pair index layer according to the relative importance degree among the elements;
respectively carrying out hierarchical single sequencing and consistency check on the first judgment matrix, and respectively carrying out hierarchical single sequencing and consistency check on the second judgment matrix;
performing hierarchical total sequencing and consistency verification on the weights of the target layers passing the consistency verification and the target layers passing the consistency verification, and performing hierarchical total sequencing and consistency verification on the weights of the index layers by the criterion layers passing the consistency verification;
the weight passing the consistency check is taken as the subjective weight.
Specifically, the method compares the importance of each element of the same layer with respect to a certain criterion in the previous layer in pairs to obtain a first judgment matrix and a second judgment matrix. Illustratively, the relative importance of the evaluation indexes of each level is qualitatively described by using a 1-9 scale method, and the relative importance is represented by accurate numbers, so that a first judgment matrix and a second judgment matrix are obtained. It should be understood that the first judgment matrix and the second judgment matrix include at least one judgment matrix, such as a group of judgment matrices of the target layer alignment rule layer and three groups of judgment matrices of the criterion layer pair index layer. The definition of the scale is referred to in Table I. The expression of the judgment matrix A-B is shown in the second table, wherein A is as follows 1 B is a certain criterion in the last hierarchy 1 -B n Is the elements of the next hierarchy.
List one
Watch II
A k B 1 B 2 B i B n
B 1 1 b 12 b 1i b 1n
B 2 b 21 1 b 2i b 2n
B i b i1 b i2 1 b in
B n b n1 b n2 b ni 1
Further, the purpose of the hierarchical single ranking is to determine the weight value of each evaluation index importance order that the present hierarchy has a relationship with an element in the upper hierarchy. Maximum characteristic root lambda of judgment matrix max Is normalized and is marked as W.
I.e. for the decision matrix B, the calculation satisfies bw=λ max Feature root and feature vector of W.
In the above formula lambda max To determine the largest eigenvalue of matrix B, W is the corresponding lambda max The component of W is the weight value corresponding to the rank of the evaluation index sheet.
Further, the index defining the consistency is:
if ci=0, there is complete agreement; CI is close to 0, and satisfactory consistency is achieved; the larger the CI, the more serious the inconsistency. To measure the size of CI, a random uniformity index RI is introduced:
the random consistency index RI is related to the order of the judgment matrix, and in general, the larger the matrix order is, the greater the probability of occurrence of consistency random deviation is, and the corresponding relationship is shown in table three.
Watch III
Considering that the deviation of the consistency may be caused by random reasons, when checking whether the judging matrix has satisfactory consistency, the CI and the random consistency index RI are also required to be compared to obtain a checking coefficient CR, and the formula is as follows:
If CR <0.1, the judgment matrix is considered to pass the consistency test, otherwise, the consistency is not satisfied, and the first judgment matrix and the second judgment matrix need to be corrected.
Further, on the basis of the above processing, the weight of all factors of a certain hierarchy for the relative importance of the highest hierarchy (total target) is continuously calculated, which is called the hierarchy total ranking. This process is performed sequentially from the highest level to the lowest level.
Element a, assuming that the total ordering of the layers of the previous layer a has been completed 1 、A 2 、…、A m The obtained weight values are respectively a 1 、a 2 、…、a m The method comprises the steps of carrying out a first treatment on the surface of the And A is a j Corresponding element B of the hierarchy 1 、B 2 、…、B n The hierarchical single ordering result of (2) isAnd when B i And A is a j No contact->The overall ranking result for level B at this time is shown in table four.
Table four
Thus, subjective weight of the evaluation index set can be obtained.
Optionally, an entropy weighting method is adopted, and objective weights of the elements are calculated based on the evaluation values, including:
constructing a third judgment matrix of n evaluation indexes of m samples based on the evaluation values;
sequentially carrying out standardization processing and normalization processing on the third judgment matrix to obtain a normalized third judgment matrix;
calculating entropy of each evaluation index according to the normalized third judgment matrix;
And taking the entropy weight of each element obtained by calculation as objective weight according to the entropy of each evaluation index.
Specifically, the entropy weighting method is an objective weighting method, and is based on the following principle: the smaller the variation degree of the index, the smaller the reflected information quantity, and the lower the corresponding weight value, and the entropy weight method is used for distributing the weight according to the variation degree of one index.
Optionally, before constructing the third judgment matrix, the method further includes: and (3) correcting each evaluation value in the evaluation index set by adopting a Heidecke method, and processing positive and negative indexes to normalize data, so that the calculation and evaluation of the weight sustainable development capability are facilitated, and the data are in normal distribution.
The expression of the forward index processing is as follows:
the expression of the negative index processing is:
above x ij For each index data of evaluation value, x imax For maximum value, x in each index data imin V, which is the minimum value in each item of data ij Is normalized data.
The standardized data are arranged, and a third judgment matrix X of n evaluation indexes of m samples is constructed:
carrying out normalization processing on the third judgment matrix X, and constructing a normalization matrix Y of m rows and n columns by utilizing maximum data and minimum data of each row and each column:
Calculating entropy H of jth evaluation index j
Wherein y is ij For normalizing the elements in matrix Y. In the above formula, if f ij =0, take lnf ij =0。
Thus, according to the entropy of the j-th evaluation index, the entropy weight of the j-th evaluation index is calculated,
wherein the j-th evaluation index is any one of the n evaluation indexes:
when j=1, 2,..n, and 0.ltoreq.ω j ≤1,∑ω j =1。
Combining the subjective weight obtained by AHP and the objective weight obtained by entropy weight method to obtain comprehensive weight a j
Wherein a is j Is the comprehensive weight of the jth index, h j To calculate subjective weight, ω, of the jth index by AHP j In order to calculate objective weight of the jth index through the entropy weight method, n is the index number.
And synthesize weight a j Satisfy the following requirements
In the above embodiment, the weight of the evaluation index set is corrected by adopting the method of combining the AHP and the entropy weight method, so that the weight distribution of the evaluation index set is more reasonable.
Optionally, according to the evaluation matrix R and the weight, a synthesis matrix is calculated, including:
and multiplying the evaluation matrix R and the weight to obtain a composite matrix B=A×R. Further, in order to comprehensively consider the influence of each factor on the sustainable development level, a weighted average type compound operation is carried out on the synthesis matrix B, a final decision set is obtained, and a comment corresponding to the maximum value in the decision set is selected as an evaluation result based on the maximum membership principle.
According to the method for evaluating the coordinated development of the flood storage area, the evaluation index set related to flood, economy and ecology of the flood storage area is determined based on the fuzzy comprehensive evaluation method, and the evaluation value of the evaluation index set is obtained through calculation according to the historical data of the evaluation index set and a pre-constructed value accounting model; on the basis of the evaluation index set, an AHP algorithm is adopted to construct a coordinated development evaluation system of the flood storage area, an AHP-entropy weight method is adopted to calculate and obtain weights of all elements in the coordinated development evaluation system, a final decision set is obtained according to the weights, finally, based on a maximum membership principle, a comment corresponding to the maximum value in the decision set is selected as an evaluation result, and an objective and more accurate analysis result is provided for sustainable development of the flood storage area.
It should be understood that, although the steps in the flowcharts related to the embodiments described above are sequentially shown as indicated by arrows, these steps are not necessarily sequentially performed in the order indicated by the arrows. The steps are not strictly limited to the order of execution unless explicitly recited herein, and the steps may be executed in other orders. Moreover, at least some of the steps in the flowcharts described in the above embodiments may include a plurality of steps or a plurality of stages, which are not necessarily performed at the same time, but may be performed at different times, and the order of the steps or stages is not necessarily performed sequentially, but may be performed alternately or alternately with at least some of the other steps or stages.
Based on the same inventive concept, the embodiment of the application also provides a diapause area coordinated development evaluation system for realizing the diapause area coordinated development evaluation method. The implementation of the solution provided by the system is similar to the implementation described in the above method, so the specific limitation in the embodiments of the system for evaluating the coordinated development of the flood storage area provided below may be referred to the limitation of the method for evaluating the coordinated development of the flood storage area hereinabove, and will not be repeated here.
Referring to fig. 3, in one possible embodiment, an embodiment of the present application provides a coordinated development evaluation system for a flood storage area, including: the system comprises an index determining module, a calculating module, an evaluation system constructing module, an evaluation value calculating module and an evaluation module, wherein:
the index determining module is used for determining an evaluation index set and a comment set of the flood storage area based on a fuzzy comprehensive evaluation method; wherein the evaluation index set includes a flood index, an economic index, and an ecological index.
The computing module is used for determining a membership function of the evaluation index set; according to the historical data of the evaluation index set, calculating to obtain an evaluation value of the evaluation index set; and establishing an evaluation matrix according to the evaluation value and the membership function.
The evaluation system construction module is used for constructing a coordinated development evaluation system of the flood storage area according to the AHP algorithm, wherein the coordinated development evaluation system comprises three layers, a first layer is a target to be realized, a second layer and a third layer comprise at least one element, and the third layer element is an evaluation index set.
And the calculation module is also used for calculating the weight of each element by adopting an AHP-entropy weight method according to the evaluation value.
The evaluation module is used for calculating to obtain a synthetic matrix according to the evaluation matrix and the weight; selecting the comment corresponding to the maximum value in the synthesis matrix as an evaluation result.
Specifically, the evaluation index set U determined by the index determination module includes the flood index U 1 Economic index u 2 And an ecological index u 3 I.e. u= { U 1 ,u 2 ,u 3 }. Wherein, flood control subsystem u 1 The corresponding flood indicator includes an activation frequency u 11 Maximum submerged depth u 12 Flood protection area u 13 Flood storage amount u 14 Emergency transfer setting population u 15 Qualification rate u of dike 16 Quantity u of sluice 17 Any one or more of the following; economic subsystem u 2 The corresponding economic index comprises natural population growth rate u 21 Rate u of illiterate population 22 Area u of arable land 23 Water resource u available for everyone 24 GDP of people's average u 25 GDP growth Rate u 26 The average resident person can control income u 27 Area u of average residence 28 Urban registration loss rate u 29 Any one or more of the following; ecological subsystem u 3 The corresponding ecological index comprises water source conservation u 31 Climate control u 32 U for purifying air 33 U for purifying water 34 Water and soil conservation u 35 Carbon fixation u 36 Species conservation u 37 Any one or more of the following.
Further, the index determination module establishes a comment set v= { V1, V2, …, vm }, which is a set of comments made by the object to be evaluated, for representing the degree of merit of the evaluation factor. If the flood storage area can be continuously developed, five grades of V= { good, better, general, worse and bad } are adopted.
Optionally, the calculating module calculates an evaluation value of the evaluation index set according to the historical data of the evaluation index set, including:
and taking the numerical value meeting the first preset condition in the preset time as an evaluation value of the flood index according to the historical data of the flood index.
And taking the numerical value meeting the second preset condition in the preset time as an evaluation value of the economic index according to the historical data of the economic index.
According to the historical data of the ecological indexes, selecting the numerical value meeting a third preset condition in the preset time, and inputting the selected numerical value into a pre-constructed value accounting model to obtain the evaluation value of each ecological index.
Further, when the ecological index includes water conservation, climate adjustment, air purification, and water purification, the value accounting model includes a water conservation sub-model, a climate adjustment sub-model, an air purification sub-model, and a water purification sub-model.
The constructed water source conservation sub-model comprises:
Q l =Q b +Q x +Q f +Q j -Q t -Q s -Q z
V w =Q l ×P s ×C s +Q l ×P y ×C y +Q l ×P g ×C g
P s +P y +P g ≤100%
wherein Q is l For accumulating water resource quantity, Q of flood area b Is the surface water quantity, Q x Make-up quantity, Q for groundwater penetration f Is of flood quantity, Q j Is precipitation amount, Q t For flood removal, qs for maintaining the water quantity consumed by the ecological system, Q z Vapor deposition for the area; v (V) w Evaluation value for conservation of water source, P s Percentage of water needed to maintain the ecosystem and the amount of water available, C s Market price, P, of water resource trade for maintaining water consumed by ecosystem y The water consumption is the percentage of the applicable water consumption in the cultivation water consumption occupying area, C y Trade market price, P, for water resources of aquaculture water g For the percentage of available water resources and C in the irrigation water occupying area g Market prices are traded for water resources of irrigation water.
The constructed climate conditioner sub-model comprises:
Q z =(E f1 +E f2 )×A
E f1 =E s +E w +E p
E f2 =S×T
V q =E tt ×P e
E tt =E z ×q×10 3 /3600+E z ×y
wherein Q is z For total evaporation of stagnant zone, E f1 For the evapotranspiration of the area after flood withdrawal, E f2 The evaporation capacity of the area during flood diversion, and the area A is the area of the area; e (E) s Soil vapor emission amount E for region w Vapor emission amount of water surface of region E p The method is characterized in that the method comprises the steps of transpiration of plants in a region, S is the evaporation rate of water surface in the region when flood storage areas are used for flood diversion, and T is the evaporation time; v (V) q For climate control evaluation value, E tt Energy, P, consumed for evaporation of ecosystems e To the local electricity price、E z Is the evaporation capacity, q is the latent heat of volatilization, y is 1m 3 The water is converted into electricity consumption of steam.
The constructed air purification submodel comprises:
wherein Q is w Is the total emission quantity, Q of the atmospheric pollutants i Is the discharge amount of the i-type atmospheric pollutants and Q y Air purifying capacity, Q of regional ecological system ij Purification amount of the ith treatment unit for the jth class of atmospheric pollutants, A i Is the area of the ith processing unit; v (V) i Value, Q for purifying class i atmospheric environment i Purifying amount of class i atmospheric pollutants, C i The treatment cost and V of the i-th type atmospheric pollutants k An evaluation value for air purification.
The constructed water quality purification sub-model comprises the following components:
wherein Q is i To purify the i-th pollutant, P i For the treatment capacity of each treatment unit of the flood storage area to the ith pollutant, S i For the area of the processing unit, T is the processing time, Q c Is the total purification amount of pollutants, V i For purifying the water quality of the i-th pollutant and C i Unit treatment cost for the ith pollutant, V z An evaluation value for water quality purification.
Optionally, when the ecological index further includes soil conservation, carbon fixation, and species conservation, the value accounting model further includes a soil conservation sub-model, a carbon fixation sub-model, a species conservation sub-model.
The constructed soil conservation sub-model comprises the following steps:
V s =V f +V w
V f =(Q d -Q l )×P 1
V w =(Q d -Q l )×P 2
wherein V is s Maintaining an evaluation value, V, for soil f Is the nutritive value of sediment, V w Can be converted into the value of soil, Q d Amount of deposit, Q l For the amount of lost soil, P 1 To the price and P of the fertilizer which needs to be manually input 2 The cost of the soil which needs to be manually input is high.
The constructed carbon stator model comprises:
V c =Q c ×C c
wherein V is c For the evaluation value of carbon sequestration, Q c Is the total carbon content, C c For carbon trade price, S i Ability to fix carbon for individual land types, A i The area of land is utilized for each land type.
The constructed species conservation sub-model comprises:
V s =V b +V e
V e =Q e ×P e
wherein V is s Conservation evaluation value, V for total species b Is of value for species conservation, V e For the cost of the artificially added species, S i For the number of various animals and plants, P i Cost, Q of artificial conservation for each species of conservation e P, an amount of additional valuable species e Is the value of the species requiring manual input.
Optionally, the evaluation system construction module constructs a coordinated development evaluation system of the flood storage area according to an AHP algorithm, including:
the evaluation system construction module constructs a coordinated development evaluation system of the flood storage area on the basis of an evaluation index set according to an AHP algorithm, wherein the coordinated development evaluation system comprises a target layer, a criterion layer and an index layer, the target layer is the sustainable development of the flood storage area, and the problem to be solved is represented; elements in the criteria layer include flood control subsystem, economic subsystem and ecological subsystem, characterizing the strategies employed to solve the problem; the elements in the index layer are evaluation index sets capable of reflecting the conditions of all subsystems, and represent specific indexes or specific influencing factors selected for solving the problems.
Further, the calculating module calculates the weight of each element by adopting an AHP-entropy weight method according to the evaluation value, and the calculating module comprises the following steps:
step one, calculating subjective weights of all elements by adopting AHP;
step two, calculating objective weights of all elements based on the evaluation values by adopting an entropy weight method;
and thirdly, carrying out weighted summation on the subjective weight and the objective weight, and taking the comprehensive weight obtained through calculation as the weight of the evaluation index set.
Optionally, the subjective weight of each element is calculated by using the AHP, including:
respectively constructing a first judgment matrix of a target layer alignment rule layer and a second judgment matrix of a criterion layer pair index layer according to the relative importance degree among the elements;
respectively carrying out hierarchical single sequencing and consistency check on the first judgment matrix, and respectively carrying out hierarchical single sequencing and consistency check on the second judgment matrix;
performing hierarchical total sequencing and consistency verification on the weights of the target layers passing the consistency verification and the target layers passing the consistency verification, and performing hierarchical total sequencing and consistency verification on the weights of the index layers by the criterion layers passing the consistency verification;
the weight passing the consistency check is taken as the subjective weight.
Optionally, an entropy weighting method is adopted, and objective weights of the elements are calculated based on the evaluation values, including:
constructing a third judgment matrix of n evaluation indexes of m samples based on the evaluation values;
sequentially carrying out standardization processing and normalization processing on the third judgment matrix to obtain a normalized third judgment matrix;
calculating entropy of each evaluation index according to the normalized third judgment matrix;
and taking the entropy weight of each element obtained by calculation as objective weight according to the entropy of each evaluation index.
Optionally, the evaluation module calculates to obtain a composite matrix according to the evaluation matrix and the weight; selecting the comment corresponding to the maximum value in the synthesis matrix as an evaluation result, wherein the evaluation result comprises:
and multiplying the evaluation matrix R and the weight to obtain a composite matrix B=A×R. Further, in order to comprehensively consider the influence of each factor on the sustainable development level, a weighted average type compound operation is carried out on the synthesis matrix B, a final decision set is obtained, and a comment corresponding to the maximum value in the decision set is selected as an evaluation result based on the maximum membership principle.
The coordinated development evaluation system for the flood storage areas determines an evaluation index set related to flood, economy and ecology of the flood storage areas based on a fuzzy comprehensive evaluation method, and calculates an evaluation value of the evaluation index set according to historical data of the evaluation index set and a pre-constructed value accounting model; on the basis of the evaluation index set, an AHP algorithm is adopted to construct a coordinated development evaluation system of the flood storage area, an AHP-entropy weight method is adopted to calculate and obtain weights of all elements in the coordinated development evaluation system, a final decision set is obtained according to the weights, finally, based on a maximum membership principle, a comment corresponding to the maximum value in the decision set is selected as an evaluation result, and an objective and more accurate analysis result is provided for sustainable development of the flood storage area.
The modules in the flood storage area coordinated development assessment system can be fully or partially implemented by software, hardware and combinations thereof. The above modules may be embedded in hardware or may be independent of a processor in the computer device, or may be stored in software in a memory in the computer device, so that the processor may call and execute operations corresponding to the above modules.
In one possible embodiment, a computer device is provided, which may be a server, the internal structure of which may be as shown in fig. 4. The computer device includes a processor, a memory, an Input/Output interface (I/O) and a communication interface. The processor, the memory and the input/output interface are connected through a system bus, and the communication interface is connected to the system bus through the input/output interface. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The database of the computer device is used to store XX data. The input/output interface of the computer device is used to exchange information between the processor and the external device. The communication interface of the computer device is used for communicating with an external terminal through a network connection. The computer program when executed by the processor implements a method for coordinated development assessment of a flood storage area.
In one possible embodiment, a computer device is provided, comprising a memory and a processor, the memory having stored therein a computer program, the processor, when executing the computer program, implementing the method steps in the above-described method for assessing coordinated development of a flood area.
In one possible embodiment, a computer readable storage medium is provided, on which a computer program is stored which, when executed by a processor, implements the method steps in the above-described method for assessing the coordinated development of a flood area.
The above embodiments are merely illustrative of the principles of the present invention and its effectiveness, and are not intended to limit the invention. Modifications and variations may be made to the above-described embodiments by those skilled in the art without departing from the spirit and scope of the invention. Accordingly, it is intended that all equivalent modifications and variations of the invention be covered by the claims, which are within the ordinary skill of the art, be within the spirit and scope of the present disclosure.

Claims (10)

1. A method for assessing coordinated development of a flood storage area, comprising:
determining an evaluation index set of the flood storage area based on a fuzzy comprehensive evaluation method; wherein the evaluation index set comprises flood index, economic index and ecological index;
Determining a membership function of the evaluation index set;
according to the historical data of the evaluation index set, calculating to obtain an evaluation value of the evaluation index set;
establishing an evaluation matrix according to the evaluation value and the membership function;
constructing a coordinated development evaluation system of a flood storage area according to an AHP algorithm, wherein the coordinated development evaluation system comprises three layers, a first layer is a target to be realized, a second layer and a third layer comprise at least one element, and the elements of the third layer are the evaluation index set;
according to the evaluation value, calculating to obtain the weight of each element by adopting an AHP-entropy weight method;
according to the evaluation matrix and the weight, calculating to obtain a synthetic matrix;
and selecting a comment corresponding to the maximum value in the synthesis matrix as an evaluation result.
2. The method of claim 1, wherein calculating an evaluation value of the evaluation index set based on historical data of the evaluation index set comprises:
according to the historical data of the flood index, taking the numerical value meeting a first preset condition in a preset time as an evaluation value of the flood index;
according to the historical data of the economic index, taking the numerical value meeting a second preset condition in the preset time as an evaluation value of the economic index;
And selecting a numerical value meeting a third preset condition in a preset time according to the historical data of the ecological index, and inputting the selected numerical value into a pre-constructed value accounting model to obtain an evaluation value of the ecological index.
3. The method of claim 2, wherein the ecological indicators include water conservation, climate control, air purification, and water purification;
the value accounting model comprises a water source conservation sub-model, a climate adjustment sub-model, an air purification sub-model and a water quality purification sub-model;
the water conservation sub-model includes:
Q l =Q b +Q x +Q f +Q j -Q t -Q s -Q z
V w =Q l ×P s ×C s +Q l ×P y ×C y +Q l ×P g ×C g
P s +P y +P g ≤100%
wherein Q is l For accumulating water resource quantity, Q of flood area b Is the surface water quantity, Q x Make-up quantity, Q for groundwater penetration f Is of flood quantity, Q j Is precipitation amount, Q t For flood removal, Q s To maintain the amount of water, Q, consumed by the ecosystem z Vapor deposition for the area; v (V) w Evaluation value for conservation of water source, P s Percentage of water needed to maintain the ecosystem and the amount of water available, C s Market price, P, of water resource trade for maintaining water consumed by ecosystem y The water consumption is the percentage of the applicable water consumption in the cultivation water consumption occupying area, C y Trade market price, P, for water resources of aquaculture water g For the percentage of available water resources and C in the irrigation water occupying area g Trading market prices for water resources of irrigation water;
the climate conditioner sub-model comprises:
Q z =(E f1 +E f2 )×A
E f1 =E s +E w +E p
E f2 =S×T
V q =W tt ×P e
E tt =E z ×q×10 3 /3600+E z ×y
wherein Q is z For total evaporation of stagnant zone, E f1 For the evapotranspiration of the area after flood withdrawal, E f2 The evaporation capacity of the area during flood diversion, and the area A is the area of the area; e (E) s Soil vapor emission amount E for region w Vapor emission amount of water surface of region E p The method is characterized in that the method comprises the steps of transpiration of plants in a region, S is the evaporation rate of water surface in the region when flood storage areas are used for flood diversion, and T is the evaporation time; v (V) q For climate control evaluation value, E tt Energy, P, consumed for evaporation of ecosystems e For local electricity price, E z Is the evaporation capacity, q is the latent heat of volatilization, y is 1m 3 Power consumption of water to steam;
the air purification submodel comprises:
wherein Q is w Is the total emission quantity, Q of the atmospheric pollutants i Is the discharge amount of the i-type atmospheric pollutants and Q y Air purifying capacity, Q of regional ecological system ij Purification amount of the ith treatment unit for the jth class of atmospheric pollutants, A i Is the area of the ith processing unit; v (V) i Value, Q for purifying class i atmospheric environment i Purifying amount of class i atmospheric pollutants, C i The treatment cost and V of the i-th type atmospheric pollutants k An evaluation value for air purification;
the water quality purifying sub-model comprises:
Wherein Q is i To purify the i-th pollutant, P i For the treatment capacity of each treatment unit of the flood storage area to the ith pollutant, S i For the area of the processing unit, T is the processing time, Q c Is the total purification amount of pollutants, V i For purifying the water quality of the i-th pollutant and C i Unit treatment cost for the ith pollutant, V z An evaluation value for water quality purification.
4. The method of claim 3, wherein the ecological index further comprises soil conservation, carbon fixation, and species conservation; the value accounting model also comprises a soil conservation sub-model, a carbon fixation sub-model and a species conservation sub-model;
the soil conservation sub-model comprises:
V s =V f +V w
V f =(Q d -Q l )×P 1
V w =(Q d -Q l )×P 2
wherein V is s Maintaining an evaluation value, V, for soil f Is the nutritive value of sediment, V w Can be converted into the value of soil, Q d Amount of deposit, Q l For the amount of lost soil, P 1 To the price and P of the fertilizer which needs to be manually input 2 The price of the soil which needs to be manually input is set;
the carbon fixation sub-model comprises:
V c =Q c ×C c
wherein V is c For the evaluation value of carbon sequestration, Q c Is the total carbon content, C c For carbon trade price, S i Ability to fix carbon for individual land types, A i Utilizing the area of land for each land type;
The species conservation sub-model comprises:
V s =V b +V e
V e =Q e ×P e
wherein V is s Conservation evaluation value, V for total species b Is of value for species conservation, V e For the cost of the artificially added species, S i For the number of various animals and plants, P i Cost, Q of artificial conservation for each species of conservation e P, an amount of additional valuable species e Is the value of the species requiring manual input.
5. The method of claim 1, wherein said calculating weights for each of said elements using an AHP-entropy weight method based on said evaluation values comprises:
calculating subjective weight of each element by adopting an AHP algorithm;
calculating objective weights of the elements based on the evaluation values by adopting an entropy weight method;
and determining a combined weight according to the subjective weight and the objective weight, and taking the calculated combined weight as the weight of each element.
6. The method of claim 5, wherein the ecological assessment system comprises a goal layer, a criterion layer, and an indicator layer;
the subjective weight of each element is calculated by AHP, which comprises the following steps:
respectively constructing a first judgment matrix of a target layer alignment rule layer and a second judgment matrix of a criterion layer pair index layer according to the relative importance degree among the elements;
Respectively carrying out hierarchical single sequencing and consistency check on the first judgment matrix, and respectively carrying out hierarchical single sequencing and consistency check on the second judgment matrix;
performing hierarchical total sequencing and consistency verification on the weights of the target layers passing the consistency verification and the target layers passing the consistency verification, and performing hierarchical total sequencing and consistency verification on the weights of the index layers by the criterion layers passing the consistency verification;
and taking the weight passing the consistency check as the subjective weight.
7. The method of claim 5, wherein said calculating objective weights for each of said elements based on said evaluation values using an entropy weight method comprises:
constructing a third judgment matrix of n evaluation indexes of m samples based on the evaluation values;
sequentially carrying out standardization processing and normalization processing on the third judgment matrix to obtain a normalized third judgment matrix;
calculating entropy of each evaluation index according to the normalized third judgment matrix;
and taking the calculated entropy weight of each element as the objective weight according to the entropy of each evaluation index.
8. A system for coordinated development assessment of a flood storage region, the system comprising:
The index determining module is used for determining an evaluation index set and a comment set of the flood storage area based on a fuzzy comprehensive evaluation method; wherein the evaluation index set comprises flood index, economic index and ecological index;
the computing module is used for determining a membership function between the evaluation index set and the evaluation set; according to the historical data of the evaluation index set, calculating to obtain an evaluation value of the evaluation index set; establishing an evaluation matrix according to the evaluation value and the membership function;
the evaluation system construction module is used for constructing a coordinated development evaluation system of the flood storage area according to an AHP algorithm, wherein the coordinated development evaluation system comprises three layers, a first layer is a target to be realized, a second layer and a third layer comprise at least one element, and the elements of the third layer are the evaluation index set;
the calculation module is further used for calculating the weight of each element by adopting an AHP-entropy weight method according to the evaluation value;
the evaluation module is used for calculating to obtain a synthesis matrix according to the evaluation matrix and the weight; and selecting a comment corresponding to the maximum value in the synthesis matrix as an evaluation result.
9. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor implements the steps of the method of any of claims 1 to 7 when the computer program is executed.
10. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of the method of any of claims 1 to 7.
CN202310532850.5A 2023-05-11 2023-05-11 Method, system, equipment and storage medium for evaluating coordinated development of flood storage area Pending CN116562698A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310532850.5A CN116562698A (en) 2023-05-11 2023-05-11 Method, system, equipment and storage medium for evaluating coordinated development of flood storage area

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310532850.5A CN116562698A (en) 2023-05-11 2023-05-11 Method, system, equipment and storage medium for evaluating coordinated development of flood storage area

Publications (1)

Publication Number Publication Date
CN116562698A true CN116562698A (en) 2023-08-08

Family

ID=87501397

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310532850.5A Pending CN116562698A (en) 2023-05-11 2023-05-11 Method, system, equipment and storage medium for evaluating coordinated development of flood storage area

Country Status (1)

Country Link
CN (1) CN116562698A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117391316A (en) * 2023-12-13 2024-01-12 长江水资源保护科学研究所 Pre-evaluation method for water purification capacity of flood storage area

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117391316A (en) * 2023-12-13 2024-01-12 长江水资源保护科学研究所 Pre-evaluation method for water purification capacity of flood storage area
CN117391316B (en) * 2023-12-13 2024-03-19 长江水资源保护科学研究所 Pre-evaluation method for water purification capacity of flood storage area

Similar Documents

Publication Publication Date Title
Liu et al. Measurement and assessment of carrying capacity of the environment in Ningbo, China
Longyang Assessing the effects of climate change on water quality of plateau deep-water lake-A study case of Hongfeng Lake
CN110969346B (en) River basin water ecological function partition treatment demand evaluation method based on index screening
Feng et al. A risk assessment model of water shortage based on information diffusion technology and its application in analyzing carrying capacity of water resources
O'Bannon et al. Globalization of agricultural pollution due to international trade
CN102831297A (en) Integration method for cause diagnosis of lake pollution
CN105389742B (en) Water environment sensitive area remote sensing identification method and system and bank side band construction method
Tan et al. A non-probabilistic programming approach enabling risk-aversion analysis for supporting sustainable watershed development
Li et al. Measurement methods and applications for beneficial and detrimental effects of ecological services
Hu et al. Explore the relationship between fish community and environmental factors by machine learning techniques
CN115034425B (en) Method and device for dynamically simulating and scheduling aggregate community of coupling water quantity, water quality and aquatic organisms
CN111428985B (en) Assessment method for regional water ecological bearing capacity
CN116562698A (en) Method, system, equipment and storage medium for evaluating coordinated development of flood storage area
CN115953061B (en) Water resource shortage degree estimation method, device and computer readable storage medium
CN115630802A (en) Ecological restoration space planning method combining with ecological system service supply and demand
Jácome et al. Assessment of water quality monitoring for the optimal sensor placement in lake Yahuarcocha using pattern recognition techniques and geographical information systems
Kim et al. Self-organizing map and species abundance distribution of stream benthic macroinvertebrates in revealing community patterns in different seasons
Das et al. Polluted waters of the reclaimed islands of Indian Sundarban promote more greenhouse gas emissions from mangrove ecosystem
Tursunova et al. Water resources of Kazakhstan in conditions of uncertainty
Yang et al. Optimal allocation of best management practices based on receiving water capacity constraints
CN108122077A (en) A kind of water environment safety evaluation method and device
Donyaii Evaluation of climate change impacts on the optimal operation of multipurpose reservoir systems using cuckoo search algorithm
Colloff et al. The use of historical environmental monitoring data to test predictions on cross-scale ecological responses to alterations in river flows
Small et al. A direct distribution model for regional aquatic acidification
Srđević et al. Group model for evaluating the importance of Ramsar sites in Vojvodina Province of Serbia

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