CN107886216A - A kind of land carrying capacity analysis method based on Evaluation of Land Use Suitability - Google Patents
A kind of land carrying capacity analysis method based on Evaluation of Land Use Suitability Download PDFInfo
- Publication number
- CN107886216A CN107886216A CN201710947590.2A CN201710947590A CN107886216A CN 107886216 A CN107886216 A CN 107886216A CN 201710947590 A CN201710947590 A CN 201710947590A CN 107886216 A CN107886216 A CN 107886216A
- Authority
- CN
- China
- Prior art keywords
- land
- evaluation
- area
- soil
- suitability
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 135
- 238000004458 analytical method Methods 0.000 title claims abstract description 27
- 239000002689 soil Substances 0.000 claims abstract description 41
- 238000010276 construction Methods 0.000 claims abstract description 33
- 238000011161 development Methods 0.000 claims abstract description 23
- 238000011160 research Methods 0.000 claims description 38
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 36
- 238000000034 method Methods 0.000 claims description 24
- 238000013439 planning Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000035945 sensitivity Effects 0.000 claims description 12
- 230000003628 erosive effect Effects 0.000 claims description 8
- 238000001556 precipitation Methods 0.000 claims description 8
- 239000004576 sand Substances 0.000 claims description 8
- 238000013316 zoning Methods 0.000 claims description 7
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 5
- 229910052799 carbon Inorganic materials 0.000 claims description 5
- 239000004927 clay Substances 0.000 claims description 5
- 238000004162 soil erosion Methods 0.000 claims description 5
- 230000002265 prevention Effects 0.000 claims description 4
- 238000001704 evaporation Methods 0.000 claims description 3
- 230000008020 evaporation Effects 0.000 claims description 3
- 239000000843 powder Substances 0.000 claims description 3
- 238000003860 storage Methods 0.000 claims description 3
- 230000003746 surface roughness Effects 0.000 claims description 3
- 230000002194 synthesizing effect Effects 0.000 claims description 2
- 238000009826 distribution Methods 0.000 abstract description 5
- 238000013480 data collection Methods 0.000 abstract description 2
- 238000011835 investigation Methods 0.000 abstract 1
- 238000004519 manufacturing process Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 244000144972 livestock Species 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000012732 spatial analysis Methods 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
- Y02A40/22—Improving land use; Improving water use or availability; Controlling erosion
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of land carrying capacity analysis method based on Evaluation of Land Use Suitability：Data collection and investigation are carried out to survey region；Eco-land designing evaluation points and appraisement system are determined, determines the weight of each evaluation points；Divide evaluation unit；Eco-land designing single pollution assessment；The superposition of factor figure layer and evaluation result analysis；Land supply amount calculates；Each evaluation unit land resources endowment calculates；Each evaluation unit land resource supply pressure calculates；The analysis of survey region Carrying Capacity of Land Resources.The present invention is from land supply angle, by eco-land designing subregion, it is determined that the soil amount of can be supplied to of suitable development ＆ construction, on this basis, the sustainable population of land resource and economic scale under different land use scene are determined, Carrying Capacity of Land Resources and its Study on Spatial Distribution are measured with this.
Description
Technical Field
The invention belongs to the field of ecological resource evaluation, and particularly relates to a land bearing capacity analysis method based on land utilization suitability evaluation.
Background
The research on land bearing capacity mostly focuses on calculating population bearing capacity of future production capacity of land, and the calculation methods include two types, namely a direct calculation method and a production potential deduction method. The direct calculation method is to use a system dynamics method to comprehensively consider various factors influencing the land production potential, regard the bearing capacity of the region as a whole and dynamically and quantitatively calculate the population number. The production potential deduction algorithm is used for calculating the number of land resource carrying population according to certain production conditions and living standard, and the core of the production potential deduction algorithm is the calculation of the land production potential. The traditional measuring and calculating method is mainly based on the land population bearing capacity method supplied by basic human farm and livestock husbandry products, and the land utilization suitability zoning method is selected to measure and calculate the land resource supply capacity.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provides a land bearing capacity analysis method based on land utilization suitability evaluation.
The purpose of the invention is realized by the following technical scheme.
A land bearing capacity analysis method based on land use suitability evaluation comprises the following steps:
step one, collecting and investigating data of a research area;
determining evaluation factors and evaluation systems of land ecological suitability, and determining the weight of each evaluation factor;
dividing evaluation units;
step four, land ecological suitability single factor evaluation: respectively determining a scheme according to the evaluation system and the weight of each evaluation factor, and analyzing and evaluating each common evaluation factor;
step five, superposition of the factor layer and analysis of evaluation results:
under the support of a geographic information system, synthesizing the single-factor evaluation results and superposing layers by adopting a factor weighted superposition method to obtain the total ecological suitability value in the research area; ecological suitability is divided into four levels according to the calculation result of suitability scores for each grid: suitable, more suitable, less suitable and less suitable;
step six, calculating the land supply amount: calculating the land area available for each evaluation unit, namely the area of the land available for construction by everyone, by combining the result of suitability zoning and GIS and comprehensively analyzing and comparing;
step seven, the evaluation units calculate the soil resources through innate endowments:
r = area of land for which any evaluation unit person can construct/area of land for which research area person can construct
When the area of the construction land available for any evaluation unit person is larger than or equal to the area of the construction land available for people in the research region, R =1; when the area of the all-person-constructable land of any evaluation unit is smaller than that of the all-person-constructable land of the research area, R is more than 0 and less than 1;
step eight, calculating land resource supply pressure of each evaluation unit:
p = area of construction used by any evaluation unit/area of area available for construction by evaluation unit
P is the land resource supply pressure of each evaluation unit, 0<P is less than or equal to 1, and when the used construction land area of each evaluation unit is larger than the available construction land area of each evaluation unit, P =1 is taken;
step nine, analyzing the bearing capacity of the land resources in the research area:
analyzing and researching the bearing capacity of regional land resources according to the pressure response relation between the endowment of the land resources and the supply pressure of the land resources; the land resource bearing capacity of each evaluation unit is defined as a function of the endowment of land resources and the pressure of the land resources, and the formula is as follows:
Q＝R×(1-P)
wherein Q is the bearing capacity of the land resource, and the value is between 0 and 1.
The data collected in the first step comprise the land utilization history and current data of the research area, the land utilization planning condition of the area and the land utilization planning principle.
In the second step, the evaluation factors of the land ecological suitability are divided into limiting factors and common evaluation factors,
the limiting factors comprise a basic farmland area, a water area, a buffer area of 200m around the water area, a natural protection area, an ecological sensitive area, an important mountain territory, an earthquake fracture zone, a protection forest zone, a water source area and a flood regulation and storage functional area;
the common evaluation factors comprise importance of soil protection function, importance of water source conservation function, importance of biodiversity protection function, sensitivity to water and soil loss, importance of wind prevention and sand fixation function, distance from a main road, distance from a river water area, an existing land utilization mode and gradient.
Importance of the soil protection function S pro Calculated according to the following formula:
S pro ＝NPP mean ×(1-K)×(1-F slo )
wherein, NPP mean The average value of net primary productivity of a perennial ecosystem in a research area; k is a soil erodability factor; f slo The method comprises the steps of normalizing to an evaluation area gradient grid graph between 0 and 1 according to a maximum and minimum value method;
the importance WR of the water source conservation function is calculated according to the following formula:
WR＝NPP mean ×F sic ×F pre ×(1-F slo )
wherein, F sic Is a soil seepage capability factor, F pre Interpolating and normalizing the average annual precipitation data of a plurality of years (more than 30 years) to be between 0 and 1;
protective functional importance of biodiversity S bio Calculated according to the following formula:
S bio ＝NPP mean ×F pre ×F tem ×(1-F alt )
wherein, F tem The temperature parameter is obtained by interpolation of annual average precipitation data of a plurality of years (10-30 years), and the obtained result is normalized to be between 0 and 1; f alt The altitude parameter is obtained by normalizing the altitude of the research area;
the soil erosion sensitivity is calculated according to the following formula:
wherein SS i Is the i space unit water and soil loss sensitivity index, R i For erosive power of rainfall, K i For soil erodibility, LS i Is the long and steep slope of the slope, C i Covering ground vegetation;
wind-proof sand-fixation function importance S ws Calculated according to the following formula:
S ws ＝NPP mean ×K×F q ×D
ETP i ＝0.19(20+T i ) 2 ×(1-r i )
K＝{0.2+0.3exp[-0.0256SAN(1-SIL/100)]}
×[SIL/(CLA+SIL)] 0.3 ×{1.0-0.25C/[C+exp(3.72-2.59C)]}
×{1.0-0.7SIL/[(1-SAN/100)+exp(-5.51+22.9(1-SAN/100))]}
D＝1/cos(θ)
wherein, NPP mean For researching the average value of net primary productivity of a perennial ecosystem in an area, K is a soil erodability factor, SAN, SIL and CLA are the contents (%) of sand grains, powder grains and clay grains of soil, C is the content (%) of organic carbon in soil, and F is the content (%) of organic carbon in soil q Average weather erosion force for years, u is the average monthly wind speed at 2m height, ETP i Is the monthly latent evaporation (mm), P i The monthly rainfall (mm), d is the number of days of the month, T i Mean monthly air temperature, r i The monthly average relative humidity (%), D the surface roughness factor, theta the slope, and radians.
The process of dividing the evaluation unit in the third step: according to the land utilization planning, the region development direction and the positioning characteristics of the research region, a point-by-point comprehensive evaluation method is adopted under the support of a GIS, namely, the research region is divided into grids, and each grid is used as a basic evaluation unit.
Compared with the prior art, the technical scheme of the invention has the following beneficial effects:
(1) The invention provides a novel method for measuring and calculating land bearing capacity, and the land bearing capacity is positioned in the environment protection angle and can be used for determining the threshold value of a development and construction land. And a land utilization suitability zoning method is selected to measure and calculate regional land resource supply capacity.
(2) For regional land resource bearing capacity evaluation, the traditional land resource bearing capacity evaluation relates to the number of land resources, the land resource production potential and the per capita consumption level, generally refers to the number of population which can be continuously fed by land in a certain region, namely the land resource population bearing capacity, and is essentially used for researching the balance relationship between population consumption and food production, human demand and resource supply. At present, the land resource emphasizes the land supply function for industrial and urban construction, the production mode and the value structure of land are obviously changed, and the land bearing capacity research based on the human-grain relationship can not meet the requirements of continuous regional development and space pattern optimization.
(3) Based on the actual demand of the urbanization development and the requirement of research and development of the bearing capacity, the development of comprehensive evaluation of the bearing capacity of the land resources and the research of the spatial pattern diversity rule have multiple meanings, so that a land resource bearing capacity evaluation system in a regional collaborative development mode is established, the spatial diversity rule of the bearing condition of the land resources is determined, and a beneficial reference is provided for improving the regional coordinated development level and optimizing the spatial pattern of the land.
Drawings
FIG. 1 is a flow chart of a land bearing capacity analysis method based on land use suitability evaluation of the present invention;
Detailed Description
The invention is further described below with reference to the accompanying drawings.
The invention adopts a land utilization suitability zoning method and combines the spatial analysis technology of a geographic information system to obtain a land ecological suitability zoning map. And then, calculating the land area suitable for development according to the result of the zoning of the ecological suitability of the land, thereby calculating the bearing capacity of land resources.
The land bearing capacity analysis method based on land utilization suitability evaluation is realized by a series of processes of screening and dividing evaluation factors, determining weight, evaluating a single factor, comprehensively evaluating multiple factors and the like. The method is characterized in that a land utilization current situation diagram is combined, the influence of natural attributes of land on land utilization capacity and land utilization suitability is taken as a main evaluation scale, the influence of social and economic factors is considered, the suitability of the land on a certain current situation application or a preset application is researched, the land suitability in the whole urban area range is analyzed and evaluated, the utilizability grade and the range of the land are determined, the available capacity of the existing land resources is estimated, and the land resource bearing capacity is measured and calculated. Mainly comprising the following, as shown in fig. 1.
(one) data collection
And collecting and surveying data of the research area, wherein the data comprises the land utilization history and current data of the research area, the land utilization planning condition of the area, the land utilization planning principle and the like.
(II) determining evaluation factors and evaluation systems of land ecological suitability, and determining the weight of each evaluation factor
The selection of the evaluation factors of the land ecological suitability follows the principles of dominance, stability, pertinence, measurability and actuality, and factors which have the most obvious and stable influence on the land ecological suitability and can be obtained from the existing data, data and planning are selected as much as possible. Some factors, such as ecological sensitive points, basic farmland protection areas, natural protection areas and the like, are protected or cannot be occupied in the development of areas, can cause the ecological suitability degree of the land to be zero, and are limiting factors in the evaluation of the ecological suitability of the land. When a certain evaluation unit satisfies the condition of the limiting factor, the evaluation unit is considered as a protected area and is considered as an unsuitable development area regardless of the suitability of the remaining evaluation factors. And the other factors influence evaluation unit land ecological suitability evaluation grading, and are common evaluation factors.
In combination with regional land utilization planning, the construction land is the land type mainly developed in the future, so the land ecological suitability is evaluated in a subarea mode aiming at the construction land. The selected limiting factors comprise a basic farmland area, a water area and a buffer area of 200m around the water area, a natural protection area, an ecological sensitive area, an important mountain land, an earthquake fracture zone, a shelter zone, a water source land and a flood regulation and storage functional area. The selected common evaluation factors comprise the importance of soil maintenance function, the importance of water source conservation function, the importance of biodiversity protection function, water and soil loss sensitivity, the importance of wind prevention and sand fixation function, the distance from a main road, the distance from a river water area, the existing land utilization mode and the gradient. The restriction factors are not suitable for construction and development, and the common evaluation factors establish grading standards according to relevant criteria for evaluation.
TABLE 1 evaluation factor division for land ecological suitability
Obtaining a common evaluation factor and dividing weights:
importance of the soil protection function S pro Calculated according to the following formula:
S pro ＝NPP mean ×(1-K)×(1-F slo ) (1)
wherein, NPP mean The average value of net primary productivity of a perennial ecosystem in a research area; k is a soil erodability factor; f slo The evaluation area gradient grid map (calculated by DEM by using geographic information system software) is normalized to be between 0 and 1 by a maximum-minimum method.
The importance of the water source conservation function: by referring to a method for evaluating the importance of the water source conservation function at home and abroad and combining the actual situation of a research area, 6 factors of vegetation coverage, precipitation, gradient, slope length, soil type and coverage type are selected as main factors influencing the evaluation of the importance of the water source conservation function. The influence of each factor on the water source conservation function is analyzed, the weight of the factor is determined, and the importance of the water source conservation function is comprehensively evaluated by using a weighted summation method.
WR＝NPP mean ×F sic ×F pre ×(1-F slo ) (2)
Wherein WR is a water source conservation service capability index of the ecological system; f sic For the soil seepage capability factor, according to the soil texture classification of the United States Department of Agriculture (USDA), the 13 soil texture types are respectively and equally assigned with values between 0 and 1, such as clay (heavy) of 1/13, silty clay of 2/13, … … and sand of 1.F pre Interpolating and normalizing the average annual precipitation data of a plurality of years (more than 30 years) to be between 0 and 1; NPP mean And F slo The meaning and algorithm of (a) are as above.
Importance of protection function of biodiversity S bio And (3) calculating: when species data is ambiguous or inadequate, the following methods based on habitat diversity can be employed.
S bio ＝NPP mean ×F pre ×F tem ×(1-F alt ) (3)
Wherein, F tem The temperature parameter is obtained by interpolation of annual average precipitation data of a plurality of years (10-30 years), and the obtained result is normalized to be between 0 and 1; f alt The altitude parameter is obtained by normalizing the altitude of the research area; NPP mean And F slo The meaning and algorithm of (a) are as above.
Susceptibility to soil erosion:
and evaluating the water and soil loss sensitivity by selecting evaluation indexes such as precipitation erosion force, soil erodibility, slope length, surface vegetation coverage and the like, and correspondingly adjusting the grading evaluation standard according to the actual condition of the research area. And carrying out product operation and weighted superposition on the single-factor distribution map reflecting the sensitivity of each factor to water loss and soil erosion by using a geographic information system technology to obtain a water loss and soil erosion sensitivity area distribution map of the research area.
Wherein SS i Is the i space unit water and soil loss sensitivity index, R i For erosive power of rainfall, K i For soil erodibility, LS i Is a long slope, C i Covering the vegetation on the ground.
Wind-prevention sand-fixation function importance S ws Calculated according to the following formula:
S ws ＝NPP mean ×K×F q ×D (5)
ETP i ＝0.19(20+T i ) 2 ×(1-r i ) (7)
K＝{0.2+0.3exp[-0.0256SAN(1-SIL/100)]}
×[SIL/(CLA+SIL)] 0.3 ×{1.0-0.25C/[C+exp(3.72-2.59C)]}
×{1.0-0.7SIL/[(1-SAN/100)+exp(-5.51+22.9(1-SAN/100))]} (8)
D＝1/cos(θ) (9)
wherein, NPP mean The average value of net primary productivity of a perennial ecosystem in a research area; k is a soil erodability factor; SAN, SIL and CLA are the contents (%) of sand grains, powder grains and clay grains of the soil; c is the organic carbon content (%) of the soil; f q The weather erosion force is the average weather erosion force for many years; u is the monthly mean wind speed at 2m altitude, ETP i Monthly latent evaporation (mm); p i Monthly rainfall (mm); d is the number of days in the month; t is i The average monthly air temperature; r is i Monthly average relative humidity (%); d is a surface roughness factor; theta is the slope and is in radians. K. F q Normalizing the parameters D to be between 0 and 1 and substituting the normalized parameters D into the parameters S ws And calculating the wind prevention and sand fixation service capability index.
Distance from main road: the spatial distance of land to major roads is the key to determining land development costs and location dominance. In general, the closer the land is to the road, the easier the development of the land is to form the space accumulation effect.
Distance from water area: the construction and development of the land have great influence on the ecological environment of nearby rivers and lakes. In principle, the exploitation site should be as far away from rivers and lakes as possible so as not to cause damage to the river ecosystem and pollution to the water body.
The existing land utilization mode is as follows: the existing land utilization mode also determines the land ecological suitability of the evaluation unit to a certain extent.
Gradient: the terrain is one of limiting primers influencing land construction and development, and an excessively steep slope is generally not suitable for a construction land.
Example (b): on the basis of the determination of the land use suitability evaluation factor index system and the analysis method, the land use general planning is combined, and the land use suitability is evaluated in a subarea mode according to the requirement of the construction land based on the target of the construction land which is the land type mainly developed in the future. The land suitability of the restriction factor is 0. The general evaluation factors are classified according to the degree of development and utilization, and the classification is shown in table 2. TABLE 2 suitability rating and weight assignment for each evaluation factor
(III) partition evaluation unit
The evaluation unit is a land block with relatively consistent land characters, and is the most basic unit for reading information of evaluation factors and performing statistical analysis in land suitability evaluation. According to the characteristics of land utilization planning, area development direction, positioning and the like of the research area, a point-by-point comprehensive evaluation method is adopted under the support of a GIS, namely, the research area is subjected to grid division, each grid is used as a basic evaluation unit, and comprehensive evaluation is carried out on the grids one by one. The superposition operation of the weight and the score of each factor (layer) relatively reduces the subjectivity of people.
(IV) Single factor evaluation of land ecological suitability
And analyzing and evaluating the common evaluation factors according to the evaluation system and the weight determination scheme of the evaluation factors. The result is used as the basis of the superposition analysis of the land suitability. The method comprises a gradient factor suitability analysis diagram, a water and soil loss sensitivity space distribution diagram, a land utilization mode factor suitability analysis diagram, a distance from road factor suitability analysis diagram, a water source conservation function importance space distribution diagram and the like.
Superposition of factor layer and analysis of evaluation result
Ecological suitability analysis for the land in the research area can be carried out by adopting a factor weighted superposition method under the support of a geographic information system. Aiming at various land forms set by development planning, according to the existing ecological single-factor evaluation results and corresponding weight values, the single-factor grading evaluation results are integrated and the layers are superposed to obtain the total ecological suitability value in the area, and the calculation formula of the integrated value is as follows:
E＝∑C·B (10)
where C represents the score and B represents the weight.
According to the calculation result of the suitability score of each grid, deducting the limiting factors, and dividing the ecological suitability into four levels: suitable, more suitable, less suitable and less suitable.
And selecting a reasonable development scheme and a land supply scheme through comprehensive analysis and comparison. And finally, obtaining a land ecological suitability subarea graph by using a geographic information system overlay analysis method.
(VI) calculation of land supply quantity
And analyzing and evaluating each evaluation factor respectively by combining the suitability partition result with GIS, comprehensively analyzing and comparing, calculating the land area available for each evaluation unit, namely the area of the land which can be built by everyone, and selecting a reasonable development scheme and a land supply scheme.
(VII) land resource endowment calculation of each evaluation unit
And the endowment of the land resources can be obtained through the calculation result of the land supply amount. The formula is as follows:
r = area of construction land available for any evaluation unit person/area of construction land available for study region person (11), wherein R =1 when the area of construction land available for any evaluation unit person is greater than or equal to the area of construction land available for study region person; when the area of the per-capita constructable land of any evaluation unit is smaller than that of the per-capita constructable land of the research area, R is more than 0 and less than 1.
(VIII) land resource supply pressure calculation for each evaluation unit
P = area of construction land used by any evaluation unit/area of construction land available for evaluation unit (12), wherein P is land resource supply pressure of each evaluation unit, 0<P is less than or equal to 1, and P =1 is taken when the area of construction land used by the evaluation unit is larger than the area of construction land available for evaluation unit.
(nine) analysis of land resource bearing capacity of a research area:
and analyzing and researching the bearing capacity of the regional land resources according to the pressure response relation between the endowment of the land resources and the supply pressure of the land resources. Here, the land resource bearing capacity of each evaluation unit is defined as a function of the endowment of land resources and the pressure of land resources, and the formula is as follows:
Q＝R×(1-P) (13)
wherein Q is the land resource bearing capacity, the numerical value is between 0 and 1, and the larger the numerical value is, the larger the development space is.
While the present invention has been described in terms of its functions and operations with reference to the accompanying drawings, it is to be understood that the invention is not limited to the precise functions and operations described above, and that the above-described embodiments are illustrative rather than restrictive, and that various changes and modifications may be effected therein by one skilled in the art without departing from the scope or spirit of the invention as defined by the appended claims.
Claims (5)
1. A land bearing capacity analysis method based on land use suitability evaluation is characterized by comprising the following steps:
step one, collecting and investigating data of a research area;
determining land ecological suitability evaluation factors and an evaluation system, and determining the weight of each evaluation factor;
dividing evaluation units;
step four, land ecological suitability single factor evaluation: respectively determining schemes according to the evaluation system and the weight of each evaluation factor, and analyzing and evaluating each common evaluation factor;
step five, superposition of the factor map layer and analysis of evaluation results:
under the support of a geographic information system, synthesizing the single-factor evaluation results and superposing layers by adopting a factor weighted superposition method to obtain the total ecological suitability value in the research area; according to the calculation result of the suitability score for each grid, the ecological suitability is divided into four levels: suitable, more suitable, less suitable and less suitable;
step six, calculating the land supply amount: calculating the land area available for each evaluation unit, namely the area of the land available for construction by everyone, by combining the result of suitability zoning and GIS and comprehensively analyzing and comparing;
step seven, the evaluation units calculate the soil resources through innate endowments:
r = area of land for which any evaluation unit person can construct/area of land for which research area person can construct
When the area of the construction land available for any evaluation unit person is larger than or equal to the area of the construction land available for people in the research region, R =1; when the area of the all-person-constructable land of any evaluation unit is smaller than that of the all-person-constructable land of the research area, R is more than 0 and less than 1;
step eight, calculating land resource supply pressure of each evaluation unit:
p = area of construction used by any evaluation unit/area of area available for construction by evaluation unit
P is the land resource supply pressure of each evaluation unit, 0<P is less than or equal to 1, and when the used construction land area of each evaluation unit is larger than the available construction land area of each evaluation unit, P =1 is taken;
step nine, analyzing the bearing capacity of the land resources in the research area:
analyzing and researching the bearing capacity of regional land resources according to the pressure response relation between the endowment of the land resources and the supply pressure of the land resources; the land resource bearing capacity of each evaluation unit is defined as a function of the endowment of land resources and the pressure of the land resources, and the formula is as follows:
Q＝R×(1-P)
wherein Q is the bearing capacity of the land resource, and the value is between 0 and 1.
2. The land bearing capacity analysis method based on land use suitability evaluation as claimed in claim 1, wherein the data collected in the first step comprises the data of land use history and current situation of the research area, the land use planning condition of the area and the land use planning principle.
3. The land bearing capacity analysis method based on land use suitability evaluation according to claim 1, wherein the land ecological suitability evaluation factors in the second step are divided into restriction factors and general evaluation factors,
the limiting factors comprise a basic farmland area, a water area, a buffer area of 200m around the water area, a natural protection area, an ecological sensitive area, an important mountain territory, an earthquake fracture zone, a protection forest zone, a water source area and a flood regulation and storage functional area;
the common evaluation factors comprise importance of soil protection function, importance of water source conservation function, importance of biodiversity protection function, sensitivity to water and soil loss, importance of wind prevention and sand fixation function, distance from a main road, distance from a river water area, an existing land utilization mode and gradient.
4. Land bearing capacity analysis method based on land use suitability evaluation according to claim 3, characterized by the importance of the soil protection function S pro Calculated according to the following formula:
S pro ＝NPP mean ×(1-K)×(1-F slo )
wherein, NPP mean The average value of net primary productivity of the perennial ecosystems in the research area; k is a soil erodability factor; f slo The method comprises the steps of normalizing to an evaluation area gradient grid graph between 0 and 1 according to a maximum and minimum value method;
the importance WR of the water source conservation function is calculated according to the following formula:
WR＝NPP mean ×F sic ×F pre ×(1-F slo )
wherein, F sic Is a soil seepage capability factor, F pre Interpolating and normalizing to be between 0 and 1 by average annual precipitation data of multiple years (more than 30 years);
importance of protection function of biodiversity S bio Calculated according to the following formula:
S bio ＝NPP mean ×F pre ×F tem ×(1-F alt )
wherein, F tem The temperature parameter is obtained by interpolation of annual average precipitation data of a plurality of years (10-30 years), and the obtained result is normalized to be between 0 and 1; f alt The altitude parameter is obtained by normalizing the altitude of the research area;
the soil erosion sensitivity is calculated according to the following formula:
wherein SS i Is the i space unit water and soil loss sensitivity index, R i For erosive power of rainfall, K i For soil erodibility, LS i Is a long slope, C i Covering ground vegetation;
wind-proof sand-fixation function importance S ws Calculated according to the following formula:
S ws ＝NPP mean ×K×F q ×D
ETP i ＝0.19(20+T i ) 2 ×(1-r i )
K＝{0.2+0.3exp[-0.0256SAN(1-SIL/100)]}
×[SIL/(CLA+SIL)] 0.3 ×{1.0-0.25C/[C+exp(3.72-2.59C)]}
×{1.0-0.7SIL/[(1-SAN/100)+exp(-5.51+22.9(1-SAN/100))]}
D＝1/cos(θ)
wherein, NPP mean For researching the average value of net primary productivity of a perennial ecosystem in an area, K is a soil erodability factor, SAN, SIL and CLA are the contents (%) of sand grains, powder grains and clay grains of soil, C is the content (%) of organic carbon in soil, and F is the content (%) of organic carbon in soil q Average weather erosive power over years, u is the average monthly wind speed at 2m height, ETP i Is the monthly latent evaporation (mm), P i Is the monthly rainfall (mm), d is the number of days in the month, T i Mean monthly air temperature, r i The monthly average relative humidity (%), D the surface roughness factor, theta the slope, and radians.
5. A land bearing capacity analysis method based on land use suitability evaluation according to claim 1, characterized by the process of dividing evaluation units in the third step: according to the land utilization planning, the region development direction and the positioning characteristics of the research region, a point-by-point comprehensive evaluation method is adopted under the support of a GIS, namely, the research region is divided into grids, and each grid is used as a basic evaluation unit.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710947590.2A CN107886216A (en) | 2017-10-12 | 2017-10-12 | A kind of land carrying capacity analysis method based on Evaluation of Land Use Suitability |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710947590.2A CN107886216A (en) | 2017-10-12 | 2017-10-12 | A kind of land carrying capacity analysis method based on Evaluation of Land Use Suitability |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107886216A true CN107886216A (en) | 2018-04-06 |
Family
ID=61781401
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710947590.2A Pending CN107886216A (en) | 2017-10-12 | 2017-10-12 | A kind of land carrying capacity analysis method based on Evaluation of Land Use Suitability |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107886216A (en) |
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109543950A (en) * | 2018-10-24 | 2019-03-29 | 上海市政工程设计研究总院（集团）有限公司 | A kind of sponge urban planning area ecological suitability evaluation analysis method based on GIS |
| CN111047221A (en) * | 2019-12-30 | 2020-04-21 | 武汉理工大学 | Modeling method for regional water resource ecological safety evaluation |
| CN111210142A (en) * | 2020-01-03 | 2020-05-29 | 北京大学 | Landscape integrated management unit dividing method |
| CN111310674A (en) * | 2020-02-20 | 2020-06-19 | 生态环境部卫星环境应用中心 | Remote sensing monitoring method and device for river ecological buffer zone based on grid division technology |
| CN111582686A (en) * | 2020-04-27 | 2020-08-25 | 国家海洋信息中心 | Ocean space development suitability evaluation method |
| CN111639806A (en) * | 2020-05-29 | 2020-09-08 | 河南大学 | GIS-based territorial space planning optimization method and system |
| CN111798032A (en) * | 2020-06-05 | 2020-10-20 | 武汉大学 | Refined grid evaluation method for supporting dual evaluation of homeland space planning |
| CN111881333A (en) * | 2020-07-21 | 2020-11-03 | 中国水利水电科学研究院 | Method for determining gravity center of regional soil erosion space |
| CN111950942A (en) * | 2020-10-19 | 2020-11-17 | 平安国际智慧城市科技股份有限公司 | Model-based water pollution risk assessment method and device and computer equipment |
| CN112036721A (en) * | 2020-08-21 | 2020-12-04 | 洛阳众智软件科技股份有限公司 | GDAL-based land resource development double-evaluation method, device and equipment |
| CN112598206A (en) * | 2020-03-20 | 2021-04-02 | 江西农业大学 | Ecological niche theory-based method for evaluating development and utilization suitability of selenium-rich soil resources |
| CN112836912A (en) * | 2019-11-25 | 2021-05-25 | 天津大学 | GIS-based soil pollution control partition defining method |
| CN113158591A (en) * | 2021-03-17 | 2021-07-23 | 南京信息职业技术学院 | Method for determining utilization bearing capacity of river basin land development |
| CN113340286A (en) * | 2021-06-28 | 2021-09-03 | 武汉财源通网络科技有限公司 | Method and equipment for analyzing surveying and mapping information data of land planning surveying and mapping project and computer storage medium |
| CN113361899A (en) * | 2021-06-01 | 2021-09-07 | 中煤航测遥感集团有限公司 | Farmland planning processing method and device, electronic equipment and readable storage medium |
| CN113379155A (en) * | 2021-06-29 | 2021-09-10 | 哈尔滨工业大学 | Method for estimating development suitability of biomass energy based on village and town population prediction |
| CN113537793A (en) * | 2021-07-22 | 2021-10-22 | 西南科技大学 | Method for ecological hydrological zoning of drainage basin |
| CN114549223A (en) * | 2022-02-24 | 2022-05-27 | 安徽耘联大数据有限公司 | Grid cultivated land protection method and system based on big data |
| CN115204547A (en) * | 2022-01-17 | 2022-10-18 | 农业农村部南京农业机械化研究所 | Method for determining suitability of farmland in hilly and mountainous areas based on high-resolution remote sensing image |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104794350A (en) * | 2015-04-23 | 2015-07-22 | 中国科学院地理科学与资源研究所 | System and method for evaluating comprehensive carrying capacity of region |
| CN106251020A (en) * | 2016-08-05 | 2016-12-21 | 辜寄蓉 | Land_use change layout method for optimizing based on resosurces environment loading capacity |
| CN106934488A (en) * | 2017-02-15 | 2017-07-07 | 中国农业大学 | A kind of land resources utilization calculation of pressure method and system |
| CN107220768A (en) * | 2017-05-26 | 2017-09-29 | 北京师范大学 | A kind of construction land Space Expanding method and system |
-
2017
- 2017-10-12 CN CN201710947590.2A patent/CN107886216A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104794350A (en) * | 2015-04-23 | 2015-07-22 | 中国科学院地理科学与资源研究所 | System and method for evaluating comprehensive carrying capacity of region |
| CN106251020A (en) * | 2016-08-05 | 2016-12-21 | 辜寄蓉 | Land_use change layout method for optimizing based on resosurces environment loading capacity |
| CN106934488A (en) * | 2017-02-15 | 2017-07-07 | 中国农业大学 | A kind of land resources utilization calculation of pressure method and system |
| CN107220768A (en) * | 2017-05-26 | 2017-09-29 | 北京师范大学 | A kind of construction land Space Expanding method and system |
Non-Patent Citations (5)
| Title |
|---|
| A336661132: "国家生态保护红线—生态功能红线划定技术指南(试行)", 《HTTPS://WWW.DOC88.COM/P-5147621387649.HTML?S=REL&ID=1》 * |
| 叶映 等: "基于GIS的铜鼓岭旅游用地生态适宜性评价", 《中国科技信息》 * |
| 李晓东 等: "《生态、生产、生活空间发展战略规划环境影响评价——以长沙大河西先导区为例》", 31 December 2014, 湖南大学出版社 * |
| 王成新 等: "基于GIS的土地建设承载水平研究--以福州市为例", 《中国环境科学学会环境规划专业委员会2013 年学术年会会议论文集》 * |
| 蔡东燕: "基于适宜性评价的土地承载力研究——以宁波市为例", 《中国优秀硕士学位论文全文数据库经济与管理科技辑》 * |
Cited By (28)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109543950A (en) * | 2018-10-24 | 2019-03-29 | 上海市政工程设计研究总院（集团）有限公司 | A kind of sponge urban planning area ecological suitability evaluation analysis method based on GIS |
| CN112836912A (en) * | 2019-11-25 | 2021-05-25 | 天津大学 | GIS-based soil pollution control partition defining method |
| CN111047221B (en) * | 2019-12-30 | 2023-09-12 | 武汉理工大学 | Modeling method for regional water resource ecological safety evaluation |
| CN111047221A (en) * | 2019-12-30 | 2020-04-21 | 武汉理工大学 | Modeling method for regional water resource ecological safety evaluation |
| CN111210142A (en) * | 2020-01-03 | 2020-05-29 | 北京大学 | Landscape integrated management unit dividing method |
| CN111210142B (en) * | 2020-01-03 | 2023-12-26 | 北京大学 | Landscape integrated management unit division method |
| CN111310674B (en) * | 2020-02-20 | 2020-11-10 | 生态环境部卫星环境应用中心 | Remote sensing monitoring method and device for river ecological buffer zone based on grid division technology |
| CN111310674A (en) * | 2020-02-20 | 2020-06-19 | 生态环境部卫星环境应用中心 | Remote sensing monitoring method and device for river ecological buffer zone based on grid division technology |
| CN112598206A (en) * | 2020-03-20 | 2021-04-02 | 江西农业大学 | Ecological niche theory-based method for evaluating development and utilization suitability of selenium-rich soil resources |
| CN112598206B (en) * | 2020-03-20 | 2023-10-03 | 江西农业大学 | Ecological niche theory-based selenium-rich soil resource development and utilization suitability evaluation method |
| CN111582686A (en) * | 2020-04-27 | 2020-08-25 | 国家海洋信息中心 | Ocean space development suitability evaluation method |
| CN111639806B (en) * | 2020-05-29 | 2022-08-26 | 河南大学 | GIS-based territorial space planning optimization method and system |
| CN111639806A (en) * | 2020-05-29 | 2020-09-08 | 河南大学 | GIS-based territorial space planning optimization method and system |
| CN111798032B (en) * | 2020-06-05 | 2024-03-01 | 武汉大学 | Fine grid evaluation method for supporting dual evaluation of homeland space planning |
| CN111798032A (en) * | 2020-06-05 | 2020-10-20 | 武汉大学 | Refined grid evaluation method for supporting dual evaluation of homeland space planning |
| CN111881333A (en) * | 2020-07-21 | 2020-11-03 | 中国水利水电科学研究院 | Method for determining gravity center of regional soil erosion space |
| CN112036721A (en) * | 2020-08-21 | 2020-12-04 | 洛阳众智软件科技股份有限公司 | GDAL-based land resource development double-evaluation method, device and equipment |
| CN111950942A (en) * | 2020-10-19 | 2020-11-17 | 平安国际智慧城市科技股份有限公司 | Model-based water pollution risk assessment method and device and computer equipment |
| CN113158591B (en) * | 2021-03-17 | 2024-01-23 | 南京信息职业技术学院 | Method for determining bearing capacity of basin land development and utilization |
| CN113158591A (en) * | 2021-03-17 | 2021-07-23 | 南京信息职业技术学院 | Method for determining utilization bearing capacity of river basin land development |
| CN113361899A (en) * | 2021-06-01 | 2021-09-07 | 中煤航测遥感集团有限公司 | Farmland planning processing method and device, electronic equipment and readable storage medium |
| CN113340286A (en) * | 2021-06-28 | 2021-09-03 | 武汉财源通网络科技有限公司 | Method and equipment for analyzing surveying and mapping information data of land planning surveying and mapping project and computer storage medium |
| CN113340286B (en) * | 2021-06-28 | 2022-03-08 | 广西国土资源规划设计集团有限公司 | Method and equipment for analyzing surveying and mapping information data of land planning surveying and mapping project and computer storage medium |
| CN113379155A (en) * | 2021-06-29 | 2021-09-10 | 哈尔滨工业大学 | Method for estimating development suitability of biomass energy based on village and town population prediction |
| CN113537793A (en) * | 2021-07-22 | 2021-10-22 | 西南科技大学 | Method for ecological hydrological zoning of drainage basin |
| CN115204547A (en) * | 2022-01-17 | 2022-10-18 | 农业农村部南京农业机械化研究所 | Method for determining suitability of farmland in hilly and mountainous areas based on high-resolution remote sensing image |
| CN115204547B (en) * | 2022-01-17 | 2023-09-26 | 农业农村部南京农业机械化研究所 | Hilly mountain area cultivated land mechanized determination method based on high-resolution remote sensing image |
| CN114549223A (en) * | 2022-02-24 | 2022-05-27 | 安徽耘联大数据有限公司 | Grid cultivated land protection method and system based on big data |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107886216A (en) | A kind of land carrying capacity analysis method based on Evaluation of Land Use Suitability | |
| Ye et al. | Construction of an ecological resistance surface model and its application in urban expansion simulations | |
| Güntner et al. | Modeling spatial patterns of saturated areas: An evaluation of different terrain indices | |
| Bagherzadeh | Estimation of soil losses by USLE model using GIS at Mashhad plain, Northeast of Iran | |
| Matomela et al. | Siting of rainwater harvesting potential sites in arid or semi-arid watersheds using GIS-based techniques | |
| CN113077153B (en) | Ecological space control partition method based on ecological system service supply and demand relationship | |
| Shi et al. | Evaluation system of coastal wetland ecological vulnerability under the synergetic influence of land and sea: A case study in the Yellow River Delta, China | |
| Wu et al. | Ecological risk assessment and difference analysis of pit ponds under different ecological service functions-A case study of Jianghuai ecological Economic Zone | |
| Boori et al. | Natural and eco-environmental vulnerability assessment through multi-temporal satellite data sets in Apodi valley region, Northeast Brazil | |
| CN111882244A (en) | Construction method of multi-source homeland development risk assessment system based on hierarchical framework | |
| Liu et al. | Dynamic evaluation of ecological vulnerability in a lake watershed based on RS and GIS technology | |
| CN111445116A (en) | Auxiliary compiling system for territorial space planning | |
| CN113240257A (en) | Territorial space partitioning method and device based on minimum cumulative resistance model | |
| Ashenafi | Modeling hydrological responses to changes in land cover and climate in Geba River Basin, Northern Ethiopia | |
| Qi et al. | Comprehensive evaluation and indicator system of land desertification in the Heihe River Basin | |
| Fang et al. | Research on the influence of land use change to habitat of cranes in Shengjin Lake wetland | |
| Powell et al. | Predicting loss and fragmentation of habitat of the vulnerable subtropical rainforest tree Macadamia integrifolia with models developed from compiled ecological data | |
| Wang et al. | Understanding relationships between landscape multifunctionality and land-use change across spatiotemporal characteristics: Implications for supporting landscape management decisions | |
| CN114925974A (en) | Classification gradual double-evaluation method for city and county territory space planning | |
| Chen et al. | Distribution and land use characteristics of alluvial fans in the Lhasa River Basin, Tibet | |
| Chen et al. | Water quality assessment and spatial–temporal variation analysis in Erhai lake, southwest China | |
| CN114493182A (en) | Ecological health evaluation method and device for urban river channels in plain river network area | |
| Zhu et al. | Evolution of landscape pattern and response of ecosystem service value in international wetland cities: A case study of Nanchang City | |
| Zhong et al. | Topographic effects on soil organic carbon in Louisiana watersheds | |
| Yang et al. | Multi-scale evaluation of ecological restoration effects in the riparian zone using Landsat series images from 1980 to 2019 |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180406 |