A calculation method for comprehensive index of regional resource and environmental carrying capacity

Abstract

The present invention discloses a method for calculating a comprehensive index of regional resource and environmental carrying capacity, comprising the following steps: obtaining regional resource and environmental carrying capacity measurement data composed of social and economic development data and natural resource and environmental data; constructing a windmill quadrilateral model for the resource carrying capacity index RCI, environmental carrying capacity index ECI, ecological adaptability index EAI, and social and economic development capacity index SDI of the measured area, obtaining a comprehensive index REBCI of resource and environmental carrying capacity of the measured area, and then determining the resource and environmental carrying capacity state of the measured area by comparing with the standard equilibrium index. The present invention calculates REBCI by modeling and fusing RCI, ECI, EAI, and SDI, and then obtains the restrictive state of resource and environmental carrying capacity. The present invention has important scientific and practical guiding significance for the establishment and improvement of the early warning and monitoring mechanism of regional resource and environmental carrying capacity, and the formulation and implementation of national land space planning.

Description

Technical Field

[0001] The present invention relates to the technical field of statistical data, geographic information data and urban information data processing, and in particular to a method for calculating a comprehensive index of regional resource and environmental carrying capacity. Background Art

[0002] The resource and environmental carrying capacity is one of the most basic endogenous forces of the natural ecological environment. The basis of human social production activities almost all comes from the natural environment. However, the world's resources are not infinite and have a relatively clear carrying capacity bottom line. The population and related social activities that can be accommodated and carried have a unique and clear threshold. Once the boundaries of the natural resource environment are broken, the resource environment will be broken and destroyed. Therefore, it is necessary to calculate and measure the resource and environmental carrying capacity of a region. The materials that a specific region can provide for people's production and life are limited. How to promote stable and sustainable economic development under limited conditions is worthy of in-depth discussion and research.

[0003] In the existing technology, there are many calculations on the carrying capacity of a single resource type, such as water resources, land resources, etc. The carrying capacity of a single resource is used as a representation of the regional resource and environmental carrying capacity to provide a reference for resource utilization measures. However, it should be taken into account that the natural resource environment, as a complex system, includes various types of resources such as water resources and land resources. The various types interact and influence each other, and their combined effect is the systematic representation of the regional resource and environmental carrying capacity. Therefore, it is inaccurate to only use the carrying capacity of a single resource type as a measure of the regional resource and environmental carrying capacity, and the resulting resource and environmental carrying capacity index cannot provide reasonable and reliable guidance for the formulation of regional resource utilization policies and the implementation of measures.

[0004] With the continuous development of social economy, technological revolutions continue to emerge, new technologies such as data analysis and data processing continue to emerge, and more and more fields are beginning to get involved in problem analysis based on objective data, which makes the comprehensive analysis of data modeling an important means to solve quantitative analysis. Resource and environmental carrying capacity involves many aspects such as nature, resources, ecology, society, and economy, and is a systematic and comprehensive work. System modeling based on statistical data and geographic information data, and the use of computers to process massive data can achieve efficient calculation of regional resource and environmental carrying capacity; in complex and complex data information, corresponding operation models are established for different attributes, which also makes the evaluation of regional resource and environmental carrying capacity more comprehensive and accurate. Summary of the invention

[0005] In view of this, the purpose of the present invention is to provide a method for calculating the comprehensive index of regional resource and environmental carrying capacity, by multi-dimensionally integrating various data on natural resources and environment and social and economic development, establishing a highly readable data analysis model, and improving computer processing efficiency.

[0006] In order to achieve the above-mentioned invention object, the technical solution adopted is as follows:

[0007] A method for calculating a comprehensive index of regional resource and environmental carrying capacity, comprising:

[0008] The resource pressure index PI and resource elasticity index RI evaluation index system of the tested area are constructed from four aspects: water resources, land resources, biological resources and energy resources, and the two indexes are standardized; the standardized resource pressure index PI and resource elasticity index RI are calculated using the entropy method to obtain the corresponding index; at the same time, a triangle model is constructed with the calculated resource pressure index PI and resource elasticity index RI as the two coordinate axes of the windmill model, and the resource carrying capacity index RCI of the tested area is obtained according to the area ratio of the triangle model OA1A1' and the standard triangle OAA';

[0009] The hydrological environment destructiveness index HEDI, atmospheric environment destructiveness index AEDI, soil environment destructiveness index SEDI, hydrological environment stability index HESI, atmospheric environment stability index AESI, soil environment stability index SESI were normalized respectively by using the log function conversion method; and the full arrangement polygonal graphic index evaluation method was used to construct a plane regular hexagonal model to calculate the environmental destructiveness index DI and environmental stability index SI respectively; at the same time, the calculated environmental destructiveness index DI and environmental stability index SI were used as the two coordinate axes of the windmill model to construct a triangle model, and the environmental carrying capacity index ECI of the area to be tested was obtained according to the area ratio of the triangle model OB1B1' and the standard triangle OBB';

[0010] The corresponding calculation method of the sub-index is used to calculate the ecological consumption index CNP, land stress index LSI, vegetation destruction index VDI, biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP; and the ecological stress index CI and ecological adaptability index AI are calculated respectively using different weights; at the same time, a triangle model is constructed with the calculated ecological stress index CI and ecological adaptability index AI as the two coordinate axes of the windmill model, and the ecological adaptability index EAI of the tested area is obtained according to the area ratio of the triangle model OC1C1' and the standard triangle OCC';

[0011] Perform basic calculations on the social unemployment index URI, population mortality index PMI, Engel index EGI, economic innovation index EII, green economy index GEI, and sustainable economy index SEI of the area to be tested; and use the social unemployment index URI, population mortality index PMI, and Engel index EGI to construct a three-dimensional pyramid model, and calculate the economic and social development recession index WI according to the volume ratio of the three-dimensional pyramid to the volume ratio of the standard three-dimensional pyramid; use the economic innovation index EII, green economy index GEI, and sustainable economy index SEI to construct a three-dimensional cylinder-cone model, and calculate the economic and social development promotion index EI according to the volume ratio of the three-dimensional cylinder-cone to the volume ratio of the standard cylinder; at the same time, use the calculated economic and social development recession index WI and economic and social development promotion index EI as the two coordinate axes of the windmill model to construct a triangle model, and obtain the socio-economic development index SDI of the area to be tested according to the area ratio of the triangle model OD1D1' and the standard triangle ODD';

[0012] A quadrilateral model is constructed based on the calculated resource carrying capacity index RCI, environmental carrying capacity index ECI, ecological adaptability index EAI, and socio-economic development index SDI, and a comprehensive resource and environmental carrying capacity index REBCI of the area to be measured is obtained based on the area ratio of the quadrilateral model W1X1Y1Z1 to the standard quadrilateral WXYZ; the comprehensive resource and environmental carrying capacity index REBCI of the measured region is used to determine the regional resource and environmental carrying status, and to construct an early warning monitoring status of load, balance, and surplus;

[0013] The resource pressure index PI and resource resilience index RI evaluation index system of the tested area are constructed for the four aspects of water resources, land resources, biological resources and energy resources. The regional county is used as the calculation unit, and water resources, land resources, biological resources and energy resources are used as the four primary indicators. The resource pressure index PI contains 8 secondary indicators, and the resource resilience index RI contains 8 secondary indicators. The data is standardized, mainly using the z-score standardization method, and the 16 indicators of the resource pressure index PI and the resource resilience index RI are input into the first formula to obtain the standardized value.

[0014] The first formula is:

[0015]

[0016] In the formula, x j is the average value of the jth index; S j is the standard deviation of the jth item;

[0017] Normalizing the hydrological environment destructiveness index HEDI, the atmospheric environment destructiveness index AEDI, the soil environment destructiveness index SEDI, and the hydrological environment stability index HESI, the atmospheric environment stability index AESI, and the soil environment stability index SESI respectively, including: using a log function conversion method to input the hydrological environment destructiveness index HEDI, the atmospheric environment destructiveness index AEDI, the soil environment destructiveness index SEDI, and the hydrological environment stability index HESI, the atmospheric environment stability index AESI, and the soil environment stability index SESI of the measured area into a second formula to obtain a normalized value;

[0018] The second formula is:

[0019]

[0020] Where, X i is the standardized data, x i is the original data, x max The maximum value of the original data in this column;

[0021] The calculation method of the resource and environmental carrying capacity comprehensive index REBCI of the tested area includes the following steps:

[0022] Construct a quadrilateral model W1X1Y1Z1: take the origin of the plane rectangular coordinate system as the center point of the model, and the distances from the origin of the plane rectangular coordinate system to the four directions are the resource carrying capacity index RCI, the environmental carrying capacity index ECI, the ecological adaptability index EAI, and the socio-economic development index SDI of the area to be measured to construct a quadrilateral model, and obtain the area of ​​the quadrilateral model;

[0023] The resource and environmental carrying capacity comprehensive index REBCI of the measured area is the ratio between the area of ​​the quadrilateral W1X1Y1Z1 calculated by the quadrilateral model and the area of ​​the standard quadrilateral WXYZ with a distance of 1 from the origin to the four directions.

[0024] As a further improvement of the present invention, when the regional resource and environmental carrying capacity comprehensive index REBCI is between 0.6 and 1, the resource and environmental carrying capacity warning state of the tested region is a surplus state;

[0025] When the regional resource and environmental carrying capacity comprehensive index REBCI is between 0 and 0.4, the resource and environmental carrying capacity warning state of the tested region is a load state;

[0026] When the regional resource and environmental carrying capacity comprehensive index REBCI is between 0.4 and 0.6, the resource and environmental carrying capacity warning state of the tested area is a balanced state.

[0027] As a further improvement of the present invention, the entropy method is used to calculate the resource pressure index PI and the resource elasticity index RI for the standardized index data, and the windmill triangle model is constructed with the resource pressure index PI and the resource elasticity index RI, as follows:

[0028] Construct an entropy value calculation model, and use the standardized data to calculate the resource pressure index PI and resource elasticity index RI according to the following formula;

[0029] ① Determine the indicator weights: The entropy weight method objectively assigns weights and calculates the indicator weights:

[0030]

[0031] Where W i 、E i , P ij Represent the weight of the indicator, information entropy and feature ratio respectively; define ln0=0;

[0032] ② Calculate the resource pressure index PI and resource elasticity index RI:

[0033] PI=W PI ×M PI RI = W RI ×M RI

[0034] Where W PI , W RI are the weights of the resource pressure index PI and resource elasticity index RI, M PI 、M RI They are the standardized indicator data of resource pressure index PI and resource elasticity index RI respectively;

[0035] Construct a windmill triangle model, take point O as the origin of the windmill single-leaf triangle model, take 45° as the angle, take the resource pressure index PI and the resource elasticity index RI as the two sides of the angle, form the first windmill single-leaf triangle model, and calculate the area of ​​the first windmill single-leaf triangle model The calculation formula is as follows:

[0036]

[0037] The resource carrying capacity index RCI of the area to be measured is obtained by the ratio of the area of ​​the first windmill single-blade triangle model to the area of ​​the first standard triangle model with an angle side of 1.

[0038] As a further improvement of the present invention, a full-array polygonal graphical index evaluation method is used to construct a plane regular hexagonal model to calculate the environmental destructive force index DI and the environmental stability index SI respectively; at the same time, a triangular model is constructed with the calculated environmental destructive force index DI and environmental stability index SI as the two coordinate axes of the windmill model, as follows:

[0039] Characterize the hydrological environment destructiveness index HEDI, atmospheric environment destructiveness index AEDI, soil environment destructiveness index SEDI, as well as the hydrological environment stability index HESI, atmospheric environment stability index AESI, soil environment stability index SESI;

[0040] The total wastewater discharge X1 and chemical oxygen demand X2 are used to represent the hydrological environment destructiveness index HEDI, the sulfur dioxide emission X3 and carbon emission X4 are used to represent the atmospheric environment destructiveness index AEDI, and the industrial hazardous solid waste generation X5 and industrial general solid waste generation X6 are used to represent the soil environment destructiveness index SEDI; the sewage treatment rate Y1 and wastewater treatment volume Y2 are used to represent the hydrological environment stability index HESI, the industrial sulfur dioxide removal volume Y3 and annual carbon emission reduction volume Y4 are used to represent the atmospheric environment stability index AESI, and the industrial general solid waste disposal volume Y5 and industrial hazardous solid waste disposal volume Y6 are used to represent the soil environment stability index SESI;

[0041] A regular hexagonal model of the environmental destructive force index DI and the environmental stability index SI was constructed. The normalized total wastewater discharge X1, chemical oxygen demand X2, sulfur dioxide emissions X3, carbon emissions X4, industrial hazardous solid waste generation X5, and industrial general solid waste generation X6 were used as the six axes radiating from the origin of the regular hexagonal plane. The environmental destructive force index DI was calculated using the following formula:

[0042]

[0043] Where, N represents the number of sub-items of the environmental destructive index DI; k i , k j is the index value of item i and item j in the index system of environmental destructiveness index DI; M represents the number of sub-indexes of environmental stability index SI; h m 、h n It is the index value of item m and item n in the environmental stability index SI index system;

[0044] Construct a windmill triangle model, take point O as the origin of the windmill single-leaf triangle model, take 45° as the angle, take the environmental destructive force index DI and the environmental stability force index SI as the two sides of the angle, form a second windmill single-leaf triangle model, and calculate the area of ​​the second windmill single-leaf triangle model The calculation formula is as follows:

[0045]

[0046] The environmental capacity index ECI of the area to be measured is obtained by the ratio of the area of ​​the second windmill single-blade triangle model to the area of ​​the second standard triangle model with an angle side of 1.

[0047] As a further improvement of the present invention, the corresponding calculation method of the sub-index is used to calculate the ecological consumption index CNP, the land stress index LSI, the vegetation destruction index VDI, the biodiversity index BI, the ecological carrying index EBI, and the ecological supply index SNP; and the ecological stress index CI and the ecological adaptability index AI are calculated respectively using different weights; at the same time, the calculated ecological stress index CI and the ecological adaptability index AI are used as the two coordinate axes of the windmill model to construct a triangular model, as follows:

[0048] Calculate the ecological consumption index CNP, land stress index LSI, vegetation destruction index VDI, biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP. The steps include:

[0049] The ecological consumption index CNP is represented by the normalized ecological consumption, and the land stress index LSI is calculated according to the following formula:

[0050]

[0051] In the formula, Aero is the normalized coefficient of land stress index, SE is the area of ​​slight erosion, ME is the area of ​​moderate erosion, LE is the area of ​​severe erosion, and AS is the area of ​​calculation area;

[0052] Calculate the vegetation destruction index VDI, the calculation formula is as follows:

[0053] VDI=Aveg(0.5FD+0.3GD+0.2FA) / AS

[0054] In the formula, Aveg is the normalized coefficient of vegetation destruction index, FD is the forest degradation area, GD is the grassland degradation area, FA is the farmland area, and AS is the calculation area;

[0055] Calculate the biodiversity index BI, the calculation formula is as follows:

[0056]

[0057] Where N is the total number of individuals; n i is the number of individuals in the ith population;

[0058] The normalized ecological supply quantity is used to represent the ecological supply index SNP; the ecological carrying index EBI is calculated according to the following formula;

[0059]

[0060] Where Pa is the total population of the measurement area, CNPP ave is the per capita ecological consumption, SNP is the ecological supply index;

[0061] The ecological stress index CI and ecological adaptability index AI are calculated using different weights. The steps include:

[0062] Different weights are assigned according to the internationally accepted urbanization level classification, as follows:

[0063] When the urbanization rate is between 0% and 30%, the weights of the ecological consumption index CNP, land stress index LSI, and vegetation destruction index VDI are 0.60, 0.25, and 0.15 respectively, and the weights of the biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP are 0.5, 0.3, and 0.2 respectively;

[0064] When the urbanization rate is between 30% and 70%, the weights of the ecological consumption index CNP, land stress index LSI, and vegetation destruction index VDI are 0.5, 0.3, and 0.2 respectively, and the weights of the biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP are 0.4, 0.3, and 0.3 respectively;

[0065] When the urbanization rate is between 70% and 100%, the weights of the ecological consumption index CNP, land stress index LSI, and vegetation destruction index VDI are 0.3, 0.6, and 0.1 respectively, and the weights of the biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP are 0.2, 0.55, and 0.25 respectively;

[0066] According to the corresponding weights, the ecological stress index CI and ecological adaptability index AI are calculated;

[0067] CI=W CNP ·CNP+W LSI ·LSI+W VDI ·VDI

[0068] AI=W BI BI+W EBI ·EBI+W SNP ·SNP

[0069] Where W CNP , W LSI , W VDI are the weights of ecological consumption index CNP, land stress index LSI, vegetation destruction index VDI, and W BI , W EBI, W SNP They are the weights of the biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP;

[0070] Construct a windmill triangle model, take point O as the origin of the windmill single-leaf triangle model, take 45° as the angle, take the ecological stress index CI and the ecological adaptability index AI as the two sides of the angle, form a third windmill single-leaf triangle model, and calculate the area S of the third windmill single-leaf triangle model △ oC1C1', the calculation formula is as follows:

[0071]

[0072] The ecological adaptability index EAI of the area to be measured is obtained by calculating the ratio of the area of ​​the third windmill single-blade triangle model to the area of ​​the third standard triangle model with an angle of 1.

[0073] As a further improvement of the present invention, basic calculations are performed on the social unemployment index URI, the population mortality index PMI, the Engel index EGI, the economic innovation index EII, the green economic index GEI, and the sustainable economic index SEI of the area to be measured; and a three-dimensional pyramid model is constructed using the social unemployment index URI, the population mortality index PMI, and the Engel index EGI, and the economic and social development recession index WI is calculated according to the volume ratio of the three-dimensional pyramid to the volume ratio of the standard three-dimensional pyramid; a three-dimensional cylinder-cone model is constructed using the economic innovation index EII, the green economic index GEI, and the sustainable economic index SEI, and the economic and social development promotion index EI is calculated according to the volume ratio of the three-dimensional cylinder-cone to the volume ratio of the standard three-dimensional cylinder; and a triangular model is constructed using the calculated economic and social development recession index WI and the economic and social development promotion index EI as the two coordinate axes of the windmill model, as follows:

[0074] The steps of representing and calculating the social unemployment index URI, the population mortality index PMI, the Engel index EGI, the economic innovation index EII, the green economy index GEI, and the sustainable economy index SEI include:

[0075] The unemployment rate represents the social unemployment index URI, the mortality rate represents the population mortality index PMI, the Engel coefficient represents the Engel index EGI, and the economic innovation index EII is calculated according to the following formula;

[0076] EII=0.5NPG+0.5EPD R&D

[0077] Where NPG is the normalized regional patent authorization data, EPD is R&D Normalize R&D expenditure data for the region;

[0078] Calculate the green economy index GEI using the following formula:

[0079] GEI = 0.65EPA + 0.35EPFE

[0080] In the formula, EPA is the normalized number of regional environmental protection agencies, and EPFE is the normalized regional environmental protection fiscal expenditure data;

[0081] Calculate the Sustainable Economy Index SEI using the following formula:

[0082] SEI=(OV eti +E i +PC GDP ) / 3,E i =(PSE+SSE+USE) / 3

[0083] Where, OV eti is the normalized high-tech industry output value, E i For the education index, PC GDP is the normalized GDP per capita, PSE, SSE, and USE are primary school, secondary school, and university enrollment rates, respectively;

[0084] A three-dimensional pyramid model is constructed, wherein the origin of the three-dimensional space is a vertex of the bottom surface of the three-dimensional pyramid model, the distances from the origin of the three-dimensional space to the bottom surface in two-dimensional directions are respectively the social unemployment index URI and the population death index PMI of the area to be measured, and the Engel index EGI is the height of the pyramid. The pyramid model is constructed, and the volume of the pyramid model is obtained; the economic and social development recession force index WI of the area to be measured is the ratio of the volume of the pyramid model to the volume of a standard pyramid whose bottom side length and height are both 1;

[0085] Construct a three-dimensional cylinder-cone model, select any point as the center point of the bottom surface of the three-dimensional cylinder-cone, use the economic innovation index EII as the radius of the lower bottom surface of the cylinder-cone, use the green economic index GEI as the radius of the upper bottom surface of the cylinder, and use the sustainable economic index SEI as the height of the cylinder-cone to construct a three-dimensional cylinder-cone model, and obtain the volume of the three-dimensional cylinder-cone model; the economic and social development promotion index EI of the area to be measured is the ratio of the volume of the three-dimensional cylinder-cone model to the volume of a standard three-dimensional cylinder model with upper and lower bottom radii and a height of 1;

[0086] Construct a windmill triangle model, take point O as the origin of the windmill single-leaf triangle model, take 45° as the included angle, take the economic and social development recession force index WI and the economic and social development promotion force index EI as the two sides of the included angle, form a fourth windmill single-leaf triangle model, and calculate the area of ​​the fourth windmill single-leaf triangle model The calculation formula is:

[0087]

[0088] The socio-economic development index SDI of the area to be measured is obtained by comparing the area of ​​the fourth windmill single-blade triangle model to the area of ​​the standard triangle model with an angle of 1.

[0089] As a further improvement of the present invention, a quadrilateral model is constructed based on the calculated resource carrying capacity index RCI, environmental carrying capacity index ECI, ecological adaptability index EAI, and socio-economic development index SDI, and the resource and environmental carrying capacity comprehensive index REBCI of the measured area is obtained based on the area ratio of the quadrilateral model W1X1Y1Z1 to the standard quadrilateral WXYZ, as follows:

[0090] Constructing a quadrilateral model W1X1Y1Z1, including: taking the origin of the plane rectangular coordinate system as the center point of the model, the distances from the origin of the plane rectangular coordinate system to the four directions are the resource carrying capacity index RCI, the environmental carrying capacity index ECI, the ecological adaptability index EAI, and the socio-economic development index SDI of the area to be measured to construct a quadrilateral model, and calculating the area of ​​the quadrilateral model The calculation formula is:

[0091]

[0092] The resource and environmental carrying capacity comprehensive index REBCI of the measured area is the ratio of the area of ​​the quadrilateral model to the area of ​​a standard quadrilateral with a distance from the origin to four directions of 1.

[0093] The beneficial effects of the present invention are as follows: the present invention provides a method for calculating data modeling and data processing of a comprehensive index of resource and environmental carrying capacity based on statistical data and geographic information data, and establishes a highly readable data analysis model by multi-dimensionally integrating various data on natural resources, environment, and social and economic development, thereby improving computer processing efficiency. In addition, the patent of the present invention calculates the comprehensive index REBCI of regional resource and environmental carrying capacity by modeling and integrating RCI, ECI, EAI, and SDI. The index obtains the restrictive state of resource and environmental carrying capacity through the numerical distance from the standard equilibrium value of 1, and has important scientific and practical guiding significance for the establishment and improvement of the early warning and monitoring mechanism of regional resource and environmental carrying capacity, and the formulation and implementation of national land space planning. BRIEF DESCRIPTION OF THE DRAWINGS

[0094] The drawings constituting a part of this application are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings:

[0095] Figure 1It is a schematic flow chart of a method for calculating a comprehensive index of regional resource and environmental carrying capacity of the present invention;

[0096] Figure 2 It is a schematic diagram of the evaluation of the full-array polygonal graphic index according to the present invention, wherein (a) is the environmental destructive force index DI, and (b) is the environmental stability index SI;

[0097] Figure 3 is a schematic diagram of the relationship between the three-dimensional pyramid according to the present invention and the standard three-dimensional pyramid;

[0098] Figure 4 It is a schematic diagram of the relationship between the three-dimensional cylinder-truncated cone body and the standard three-dimensional cylinder according to the present invention;

[0099] Figure 5 is a schematic diagram of a windmill model of a sub-index of regional resource and environmental carrying capacity according to the present invention;

[0100] Figure 6 It is a schematic diagram of the relationship between the quadrilateral model constructed according to the present invention and the standard quadrilateral. DETAILED DESCRIPTION

[0101] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present invention will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0102] In order to enable those skilled in the art to better understand the solution of the present application, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present application.

[0103] like Figure 1-6 As shown, a method for calculating the comprehensive index of regional resource and environmental carrying capacity includes:

[0104] (a). Construct the resource pressure index PI and resource elasticity index RI evaluation index system of the tested area from the four aspects of water resources, land resources, biological resources and energy resources, and standardize the two indexes; use the entropy method to calculate the standardized resource pressure index PI and resource elasticity index RI to obtain the corresponding index; at the same time, use the calculated resource pressure index PI and resource elasticity index RI as the two coordinate axes of the windmill model to construct a triangle model, and obtain the resource carrying capacity index RCI of the tested area according to the area ratio of the triangle model OA1A1' and the standard triangle OAA';

[0105] The resource pressure index PI and resource elasticity index RI evaluation index system of the tested area is constructed for the four aspects of water resources, land resources, biological resources and energy resources. The regional county is used as the calculation unit, and water resources, land resources, biological resources and energy resources are used as the four primary indicators. The resource pressure index PI contains 8 secondary indicators, and the resource elasticity index RI contains 8 secondary indicators, as shown in Table 1;

[0106] Table 1 Evaluation index system of resource pressure index PI and resource elasticity index RI

[0107]

[0108] The data were standardized, mainly using the z-score standardization method, and the 16 indicators of the resource pressure index PI and the resource elasticity index RI were input into the first formula to obtain the standardized values;

[0109] The first formula is:

[0110]

[0111] In the formula, x j is the average value of the jth index; S j is the standard deviation of the jth item;

[0112] The entropy method is used to calculate the resource pressure index PI and the resource elasticity index RI for the standardized indicator data, and the windmill triangle model is constructed with the resource pressure index PI and the resource elasticity index RI, as follows:

[0113] Construct an entropy value calculation model, and use the standardized data to calculate the resource pressure index PI and resource elasticity index RI according to the following formula;

[0114] ① Determine the indicator weights: The entropy weight method objectively assigns weights and calculates the indicator weights:

[0115]

[0116] Where Wi, Ei, Pi j Represent the weight of the indicator, information entropy and feature ratio respectively; define ln0=0;

[0117] ② Calculate the resource pressure index PI and resource elasticity index RI:

[0118] PI=W PI ×M PI RI = W RI ×M RI

[0119] Where W PI , WRI are the weights of the resource pressure index PI and resource elasticity index RI, M PI 、M RI They are the standardized indicator data of resource pressure index PI and resource elasticity index RI respectively;

[0120] Construct a windmill triangle model, take point O as the origin of the windmill single-leaf triangle model, take 45° as the angle, take the resource pressure index PI and the resource elasticity index RI as the two sides of the angle, form the first windmill single-leaf triangle model, and calculate the area S of the first windmill single-leaf triangle model △ oA1A1', the calculation formula is as follows:

[0121]

[0122] The resource carrying capacity index RCI of the area to be measured is obtained by the ratio of the area of ​​the first windmill single-blade triangle model to the area of ​​the first standard triangle model with an angle side of 1.

[0123] (b). Use the log function conversion method to normalize the hydrological environment destructiveness index HEDI, the atmospheric environment destructiveness index AEDI, the soil environment destructiveness index SEDI, the hydrological environment stability index HESI, the atmospheric environment stability index AESI, and the soil environment stability index SESI respectively; and use the full arrangement polygonal graphic index evaluation method to construct a plane regular hexagonal model to calculate the environmental destructiveness index DI and the environmental stability index SI respectively; at the same time, use the calculated environmental destructiveness index DI and environmental stability index SI as the two coordinate axes of the windmill model to construct a triangle model, and obtain the environmental carrying capacity index ECI of the area to be tested according to the area ratio of the triangle model OB1B1' and the standard triangle OBB';

[0124] The measurement indicators are constructed for the hydrological environment destructiveness index HEDI, the atmospheric environment destructiveness index AEDI, the soil environment destructiveness index SEDI, the hydrological environment stability index HESI, the atmospheric environment stability index AESI, and the soil environment stability index SESI, as shown in Table 2:

[0125] Table 2 Evaluation index system of environmental destructive index DI and environmental stability index SI

[0126]

[0127]

[0128] The total wastewater discharge X1 and chemical oxygen demand X2 are used to represent the hydrological environment destructiveness index HEDI, the sulfur dioxide emission X3 and carbon emission X4 are used to represent the atmospheric environment destructiveness index AEDI, and the industrial hazardous solid waste generation X5 and industrial general solid waste generation X6 are used to represent the soil environment destructiveness index SEDI; the sewage treatment rate Y1 and wastewater treatment volume Y2 are used to represent the hydrological environment stability index HESI, the industrial sulfur dioxide removal volume Y3 and annual carbon emission reduction volume Y4 are used to represent the atmospheric environment stability index AESI, and the industrial general solid waste disposal volume Y5 and industrial hazardous solid waste disposal volume Y6 are used to represent the soil environment stability index SESI;

[0129] Normalizing the hydrological environment destructiveness index HEDI, the atmospheric environment destructiveness index AEDI, the soil environment destructiveness index SEDI, the hydrological environment stability index HESI, the atmospheric environment stability index AESI, and the soil environment stability index SESI respectively, including: using a log function conversion method to input the hydrological environment destructiveness index HEDI, the atmospheric environment destructiveness index AEDI, the soil environment destructiveness index SEDI, and the hydrological environment stability index HESI, the atmospheric environment stability index AESI, and the soil environment stability index SESI of the measured area into a second formula to obtain a normalized value;

[0130] The second formula is:

[0131]

[0132] Where, X i is the standardized data, x i is the original data, x max The maximum value of the original data in this column.

[0133] A regular hexagonal model of the environmental destructive force index DI and the environmental stability index SI was constructed. The normalized total wastewater discharge X1, chemical oxygen demand X2, sulfur dioxide emissions X3, carbon emissions X4, industrial hazardous solid waste generation X5, and industrial general solid waste generation X6 were used as the six axes radiating from the origin of the regular hexagonal plane. The environmental destructive force index DI was calculated using the following formula:

[0134]

[0135] Where, N represents the number of sub-items of the environmental destructive index DI; k i , k j is the index value of item i and item j in the index system of environmental destructiveness index DI; M represents the number of sub-indexes of environmental stability index SI; h m 、h nIt is the index value of item m and item n in the environmental stability index SI index system;

[0136] Construct a windmill triangle model, take point O as the origin of the windmill single-leaf triangle model, take 45° as the angle, take the environmental destructive force index DI and the environmental stability force index SI as the two sides of the angle, form a second windmill single-leaf triangle model, and calculate the area of ​​the second windmill single-leaf triangle model The calculation formula is as follows:

[0137]

[0138] The environmental capacity index ECI of the area to be measured is obtained by the ratio of the area of ​​the second windmill single-blade triangle model to the area of ​​the second standard triangle model with an angle side of 1.

[0139] (c). Using the corresponding calculation method of the sub-index, calculate the ecological consumption index CNP, land stress index LSI, vegetation destruction index VDI, biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP; and use different weights to calculate the ecological stress index CI and ecological adaptability index AI respectively; at the same time, use the calculated ecological stress index CI and ecological adaptability index AI as the two coordinate axes of the windmill model to construct a triangle model, and obtain the ecological adaptability index EAI of the area to be tested according to the area ratio of the triangle model OC1C1' and the standard triangle OCC'; the steps are as follows:

[0140] The corresponding calculation method of the sub-index is used to calculate the ecological consumption index CNP, land stress index LSI, vegetation destruction index VDI, biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP; and the ecological stress index CI and ecological adaptability index AI are calculated respectively using different weights; at the same time, the calculated ecological stress index CI and ecological adaptability index AI are used as the two coordinate axes of the windmill model to construct a triangular model, as follows:

[0141] Calculate the ecological consumption index CNP, land stress index LSI, vegetation destruction index VDI, biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP. The steps include:

[0142] The ecological consumption index CNP is represented by the normalized ecological consumption, and the land stress index LSI is calculated according to the following formula:

[0143]

[0144] In the formula, Aero is the normalized coefficient of land stress index, SE is the area of ​​slight erosion, ME is the area of ​​moderate erosion, LE is the area of ​​severe erosion, and AS is the area of ​​calculation area;

[0145] Calculate the vegetation destruction index VDI, the calculation formula is as follows:

[0146] VDI=Aveg(0.5FD+0.3GD+0.2FA) / AS

[0147] In the formula, Aveg is the normalized coefficient of vegetation destruction index, FD is the forest degradation area, GD is the grassland degradation area, FA is the farmland area, and AS is the calculation area;

[0148] Calculate the biodiversity index BI, the calculation formula is as follows:

[0149]

[0150] Where N is the total number of individuals; n i is the number of individuals in the ith population;

[0151] The normalized ecological supply quantity is used to represent the ecological supply index SNP; the ecological carrying index EBI is calculated according to the following formula;

[0152]

[0153] Where Pa is the total population of the measurement area, CNPP ave is the per capita ecological consumption, SNP is the ecological supply index;

[0154] The ecological stress index CI and ecological adaptability index AI are calculated using different weights. The steps include:

[0155] Different weights are assigned according to the internationally accepted urbanization level classification, as follows:

[0156] When the urbanization rate is between 0% and 30%, the weights of the ecological consumption index CNP, land stress index LSI, and vegetation destruction index VDI are 0.60, 0.25, and 0.15 respectively, and the weights of the biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP are 0.5, 0.3, and 0.2 respectively;

[0157] When the urbanization rate is between 30% and 70%, the weights of the ecological consumption index CNP, land stress index LSI, and vegetation destruction index VDI are 0.5, 0.3, and 0.2 respectively, and the weights of the biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP are 0.4, 0.3, and 0.3 respectively;

[0158] When the urbanization rate is between 70% and 100%, the weights of the ecological consumption index CNP, land stress index LSI, and vegetation destruction index VDI are 0.3, 0.6, and 0.1 respectively, and the weights of the biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP are 0.2, 0.55, and 0.25 respectively;

[0159] According to the corresponding weights, the ecological stress index CI and ecological adaptability index AI are calculated;

[0160] CI=W CNP ·CNP+W LSI ·LSI+W VDI ·VDI

[0161] AI=W BI BI+W EBI ·EBI+W SNP ·SNP

[0162] Where W CNP , W LSI , W VDI are the weights of ecological consumption index CNP, land stress index LSI, vegetation destruction index VDI, and W BI , W EBI , W SNP They are the weights of the biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP;

[0163] Construct a windmill triangle model, take point O as the origin of the windmill single-leaf triangle model, take 45° as the angle, take the ecological stress index CI and the ecological adaptability index AI as the two sides of the angle, form a third windmill single-leaf triangle model, and calculate the area S of the third windmill single-leaf triangle model △ oC1C1', the calculation formula is as follows:

[0164]

[0165] The ecological adaptability index EAI of the area to be measured is obtained by calculating the ratio of the area of ​​the third windmill single-blade triangle model to the area of ​​the third standard triangle model with an angle of 1.

[0166] (d) Perform basic calculations on the social unemployment index URI, population mortality index PMI, Engel index EGI, economic innovation index EII, green economy index GEI, and sustainable economy index SEI of the area to be tested; and use the social unemployment index URI, population mortality index PMI, and Engel index EGI to construct a three-dimensional pyramid model, and calculate the economic and social development recession index WI according to the volume ratio of the three-dimensional pyramid to the standard three-dimensional pyramid; use the economic innovation index EII, green economy index GEI, and sustainable economy index SEI to construct a three-dimensional cylinder-cone model, and calculate the economic and social development promotion index EI according to the volume ratio of the three-dimensional cylinder-cone to the volume ratio of the standard cylinder; at the same time, use the calculated economic and social development recession index WI and economic and social development promotion index EI as the two coordinate axes of the windmill model to construct a triangle model, and obtain the socio-economic development index SDI of the area to be tested according to the area ratio of the triangle model OD1D1' and the standard triangle ODD';

[0167] The steps of representing and calculating the social unemployment index URI, the population mortality index PMI, the Engel index EGI, the economic innovation index EII, the green economy index GEI, and the sustainable economy index SEI include:

[0168] The unemployment rate represents the social unemployment index URI, the mortality rate represents the population mortality index PMI, the Engel coefficient represents the Engel index EGI, and the economic innovation index EII is calculated according to the following formula;

[0169] EII=0.5NPG+0.5EOD R&D

[0170] Where NPG is the normalized regional patent authorization data, EPD is R&D Normalize R&D expenditure data for the region;

[0171] Calculate the green economy index GEI using the following formula:

[0172] GEI = 0.65EPA + 0.35EPFE

[0173] In the formula, EPA is the normalized number of regional environmental protection agencies, and EPFE is the normalized regional environmental protection fiscal expenditure data;

[0174] Calculate the Sustainable Economy Index SEI using the following formula:

[0175] SEI=(OV eti +E i +PC GDP ) / 3,E i =(PSE+SSE+USE) / 3

[0176] Where, OV eti is the normalized high-tech industry output value, E i For the education index, PC GDP is the normalized GDP per capita, PSE, SSE, and USE are primary school, secondary school, and university enrollment rates, respectively;

[0177] A three-dimensional pyramid model is constructed, wherein the origin of the three-dimensional space is a vertex of the bottom surface of the three-dimensional pyramid model, the distances from the origin of the three-dimensional space to the bottom surface in two-dimensional directions are respectively the social unemployment index URI and the population death index PMI of the area to be measured, and the Engel index EGI is the height of the pyramid. The pyramid model is constructed, and the volume of the pyramid model is obtained; the economic and social development recession force index WI of the area to be measured is the ratio of the volume of the pyramid model to the volume of a standard pyramid whose bottom side length and height are both 1;

[0178] Construct a three-dimensional cylinder-cone model, select any point as the center point of the bottom surface of the three-dimensional cylinder-cone, use the economic innovation index EII as the radius of the lower bottom surface of the cylinder-cone, use the green economic index GEI as the radius of the upper bottom surface of the cylinder, and use the sustainable economic index SEI as the height of the cylinder-cone to construct a three-dimensional cylinder-cone model, and obtain the volume of the three-dimensional cylinder-cone model; the economic and social development promotion index EI of the area to be measured is the ratio of the volume of the three-dimensional cylinder-cone model to the volume of a standard three-dimensional cylinder model with upper and lower bottom radii and a height of 1;

[0179] Construct a windmill triangle model, take point O as the origin of the windmill single-leaf triangle model, take 45° as the included angle, take the economic and social development recession force index WI and the economic and social development promotion force index EI as the two sides of the included angle, form a fourth windmill single-leaf triangle model, and calculate the area of ​​the fourth windmill single-leaf triangle model The calculation formula is:

[0180]

[0181] The socio-economic development index SDI of the area to be measured is obtained by comparing the area of ​​the fourth windmill single-blade triangle model to the area of ​​the standard triangle model with an angle of 1.

[0182] (e) Construct a quadrilateral model based on the calculated resource carrying capacity index RCI, environmental carrying capacity index ECI, ecological adaptability index EAI, and socio-economic development index SDI, and obtain the resource and environmental carrying capacity comprehensive index REBCI of the measured area based on the area ratio of the quadrilateral model W1X1Y1Z1 and the standard quadrilateral WXYZ; the measured regional resource and environmental carrying capacity comprehensive index REBCI is used to determine the regional resource and environmental carrying status and construct the early warning monitoring status of load, balance, and surplus; the steps are as follows:

[0183] Constructing a quadrilateral model W1X1Y1Z1, including: taking the origin of the plane rectangular coordinate system as the center point of the model, the distances from the origin of the plane rectangular coordinate system to the four directions are the resource carrying capacity index RCI, the environmental carrying capacity index ECI, the ecological adaptability index EAI, and the socio-economic development index SDI of the area to be measured to construct a quadrilateral model, and calculating the area of ​​the quadrilateral model The calculation formula is:

[0184]

[0185] The resource and environmental carrying capacity comprehensive index REBCI of the measured area is the ratio of the area of ​​the quadrilateral model to the area of ​​a standard quadrilateral with a distance from the origin to four directions of 1.

[0186] Build early warning monitoring status of load, balance, and surplus, including:

[0187] When the regional resource and environmental carrying capacity comprehensive index REBCI is between 0.6 and 1, the resource and environmental carrying capacity warning state of the tested region is a surplus state, which is a green warning;

[0188] When the regional resource and environmental carrying capacity comprehensive index REBCI is between 0 and 0.4, the resource and environmental carrying capacity warning state of the tested area is a load state, which is a red warning;

[0189] When the regional resource and environmental carrying capacity comprehensive index REBCI is between 0.4 and 0.6, the resource and environmental carrying capacity warning state of the tested area is a balanced state, which is a blue warning.

[0190] The present invention provides a method for calculating data modeling and data processing of a comprehensive index of resource and environmental carrying capacity based on statistical data and geographic information data. By multi-integrating various data on natural resources and environment and social and economic development, a highly readable data analysis model is established to improve computer processing efficiency. In addition, the patent of the present invention calculates the comprehensive index REBCI of regional resource and environmental carrying capacity through modeling fusion of RCI, ECI, EAI and SDI. The index obtains the restrictive state of resource and environmental carrying capacity through the numerical distance from the standard equilibrium value of 1, which has important scientific and practical guiding significance for the establishment and improvement of the early warning and monitoring mechanism of regional resource and environmental carrying capacity, and the formulation and implementation of national land space planning.

[0191] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, component splitting or combination, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for calculating the comprehensive index of regional resource and environmental carrying capacity, characterized in that: include: The resource pressure index PI and resource elasticity index RI evaluation index system of the tested area are constructed from four aspects: water resources, land resources, biological resources and energy resources, and the two indexes are standardized; the standardized resource pressure index PI and resource elasticity index RI are calculated using the entropy method to obtain the corresponding index; at the same time, a triangle model is constructed with the calculated resource pressure index PI and resource elasticity index RI as the two coordinate axes of the windmill model, and the resource carrying capacity index RCI of the tested area is obtained according to the area ratio of the triangle model OA1A1' and the standard triangle OAA'; The hydrological environment destructiveness index HEDI, atmospheric environment destructiveness index AEDI, soil environment destructiveness index SEDI, hydrological environment stability index HESI, atmospheric environment stability index AESI, soil environment stability index SESI were normalized respectively by using the log function conversion method; and the full arrangement polygonal graphic index evaluation method was used to construct a plane regular hexagonal model to calculate the environmental destructiveness index DI and environmental stability index SI respectively; at the same time, the calculated environmental destructiveness index DI and environmental stability index SI were used as the two coordinate axes of the windmill model to construct a triangle model, and the environmental carrying capacity index ECI of the area to be tested was obtained according to the area ratio of the triangle model OB1B1' and the standard triangle OBB'; The corresponding calculation method of the sub-index is used to calculate the ecological consumption index CNP, land stress index LSI, vegetation destruction index VDI, biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP; and the ecological stress index CI and ecological adaptability index AI are calculated respectively using different weights; at the same time, a triangle model is constructed with the calculated ecological stress index CI and ecological adaptability index AI as the two coordinate axes of the windmill model, and the ecological adaptability index EAI of the tested area is obtained according to the area ratio of the triangle model OC1C1' and the standard triangle OCC'; Perform basic calculations on the social unemployment index URI, population mortality index PMI, Engel index EGI, economic innovation index EII, green economy index GEI, and sustainable economy index SEI of the area to be tested; and use the social unemployment index URI, population mortality index PMI, and Engel index EGI to construct a three-dimensional pyramid model, and calculate the economic and social development recession index WI according to the volume ratio of the three-dimensional pyramid to the volume ratio of the standard three-dimensional pyramid; use the economic innovation index EII, green economy index GEI, and sustainable economy index SEI to construct a three-dimensional cylinder-cone model, and calculate the economic and social development promotion index EI according to the volume ratio of the three-dimensional cylinder-cone to the volume ratio of the standard cylinder; at the same time, use the calculated economic and social development recession index WI and economic and social development promotion index EI as the two coordinate axes of the windmill model to construct a triangle model, and obtain the socio-economic development index SDI of the area to be tested according to the area ratio of the triangle model OD1D1' and the standard triangle ODD'; A quadrilateral model is constructed based on the calculated resource carrying capacity index RCI, environmental carrying capacity index ECI, ecological adaptability index EAI, and socio-economic development index SDI, and a comprehensive resource and environmental carrying capacity index REBCI of the area to be measured is obtained based on the area ratio of the quadrilateral model W1X1Y1Z1 to the standard quadrilateral WXYZ; the comprehensive resource and environmental carrying capacity index REBCI of the measured region is used to determine the regional resource and environmental carrying status, and to construct an early warning monitoring status of load, balance, and surplus; The resource pressure index PI and resource resilience index RI evaluation index system of the tested area are constructed for the four aspects of water resources, land resources, biological resources and energy resources. The regional county is used as the calculation unit, and water resources, land resources, biological resources and energy resources are used as the four primary indicators. The resource pressure index PI contains 8 secondary indicators, and the resource resilience index RI contains 8 secondary indicators. The data is standardized, mainly using the z-score standardization method, and the 16 indicators of the resource pressure index PI and the resource resilience index RI are input into the first formula to obtain the standardized value. The first formula is: In the formula, x j is the average value of the jth index; S j is the standard deviation of the jth item; Normalizing the hydrological environment destructiveness index HEDI, the atmospheric environment destructiveness index AEDI, the soil environment destructiveness index SEDI, and the hydrological environment stability index HESI, the atmospheric environment stability index AESI, and the soil environment stability index SESI respectively, including: using a log function conversion method to input the hydrological environment destructiveness index HEDI, the atmospheric environment destructiveness index AEDI, the soil environment destructiveness index SEDI, and the hydrological environment stability index HESI, the atmospheric environment stability index AESI, and the soil environment stability index SESI of the measured area into a second formula to obtain a normalized value; The second formula is: In the formula, X i is the standardized data, x i is the original data, x max The maximum value of the original data in this column; The calculation method of the resource and environmental carrying capacity comprehensive index REBCI of the tested area includes the following steps: Construct a quadrilateral model W1X1Y1Z1: take the origin of the plane rectangular coordinate system as the center point of the model, and the distances from the origin of the plane rectangular coordinate system to the four directions are the resource carrying capacity index RCI, the environmental carrying capacity index ECI, the ecological adaptability index EAI, and the socio-economic development index SDI of the area to be measured to construct a quadrilateral model, and obtain the area of ​​the quadrilateral model; The resource and environmental carrying capacity comprehensive index REBCI of the measured area is the ratio between the area of ​​the quadrilateral W1X1Y1Z1 calculated by the quadrilateral model and the area of ​​the standard quadrilateral WXYZ with a distance of 1 from the origin to the four directions.

2. The method for calculating the comprehensive index of regional resource and environmental carrying capacity according to claim 1, characterized in that: When the regional resource and environmental carrying capacity comprehensive index REBCI is between 0.6 and 1, the resource and environmental carrying capacity warning state of the tested region is a surplus state; When the regional resource and environmental carrying capacity comprehensive index REBCI is between 0 and 0.4, the resource and environmental carrying capacity warning state of the tested region is a load state; When the regional resource and environmental carrying capacity comprehensive index REBCI is between 0.4 and 0.6, the resource and environmental carrying capacity warning state of the tested area is a balanced state.

3. The method for calculating the comprehensive index of regional resource and environmental carrying capacity according to claim 1, characterized in that: The entropy method is used to calculate the resource pressure index PI and the resource elasticity index RI for the standardized indicator data, and the windmill triangle model is constructed with the resource pressure index PI and the resource elasticity index RI, as follows: Construct an entropy value calculation model, and use the standardized data to calculate the resource pressure index PI and resource elasticity index RI according to the following formula; ① Determine the indicator weights: The entropy weight method objectively assigns weights and calculates the indicator weights: Where W i 、E i , P ij They represent the weight, information entropy and feature ratio of the indicators respectively; Define ln0 = 0; ② Calculate the resource pressure index PI and resource elasticity index RI: PI=W PI ×M PI ;RI=W RI ×M RI Where W PI , W RI are the weights of the resource pressure index PI and resource elasticity index RI, M PI 、M RI They are the standardized indicator data of resource pressure index PI and resource elasticity index RI respectively; Construct a windmill triangle model, take point O as the origin of the windmill single-leaf triangle model, take 45° as the angle, take the resource pressure index PI and the resource elasticity index RI as the two sides of the angle, form the first windmill single-leaf triangle model, and calculate the area of ​​the first windmill single-leaf triangle model The calculation formula is as follows: The resource carrying capacity index RCI of the area to be measured is obtained by the ratio of the area of ​​the first windmill single-blade triangle model to the area of ​​the first standard triangle model with an angle side of 1.

4. The method for calculating the comprehensive index of regional resource and environmental carrying capacity according to claim 1, characterized in that: The full arrangement polygonal graphic index evaluation method is used to construct a plane regular hexagonal model to calculate the environmental destructive index DI and the environmental stability index SI respectively; at the same time, the calculated environmental destructive index DI and environmental stability index SI are used as the two coordinate axes of the windmill model to construct a triangular model, as follows: Characterize the hydrological environment destructiveness index HEDI, atmospheric environment destructiveness index AEDI, soil environment destructiveness index SEDI, as well as the hydrological environment stability index HESI, atmospheric environment stability index AESI, soil environment stability index SESI; The total wastewater discharge X1 and chemical oxygen demand X2 are used to represent the hydrological environment destructiveness index HEDI, the sulfur dioxide emission X3 and carbon emission X4 are used to represent the atmospheric environment destructiveness index AEDI, and the industrial hazardous solid waste generation X5 and industrial general solid waste generation X6 are used to represent the soil environment destructiveness index SEDI; the sewage treatment rate Y1 and wastewater treatment volume Y2 are used to represent the hydrological environment stability index HESI, the industrial sulfur dioxide removal volume Y3 and annual carbon emission reduction volume Y4 are used to represent the atmospheric environment stability index AESI, and the industrial general solid waste disposal volume Y5 and industrial hazardous solid waste disposal volume Y6 are used to represent the soil environment stability index SESI; A regular hexagonal model of the environmental destructive force index DI and the environmental stability index SI was constructed. The normalized total wastewater discharge X1, chemical oxygen demand X2, sulfur dioxide emissions X3, carbon emissions X4, industrial hazardous solid waste generation X5, and industrial general solid waste generation X6 were used as the six axes radiating from the origin of the regular hexagonal plane. The environmental destructive force index DI was calculated using the following formula: Where, N represents the number of sub-items of the environmental destructive index DI; k i , k j is the index value of item i and item j in the index system of environmental destructiveness index DI; M represents the number of sub-indexes of environmental stability index SI; h m 、h n It is the index value of item m and item n in the environmental stability index SI index system; Construct a windmill triangle model, take point O as the origin of the windmill single-leaf triangle model, take 45° as the angle, take the environmental destructive force index DI and the environmental stability force index SI as the two sides of the angle, form a second windmill single-leaf triangle model, and calculate the area of ​​the second windmill single-leaf triangle model The calculation formula is as follows: The environmental capacity index ECI of the area to be measured is obtained by the ratio of the area of ​​the second windmill single-blade triangle model to the area of ​​the second standard triangle model with an angle side of 1.

5. The method for calculating the comprehensive index of regional resource and environmental carrying capacity according to claim 1, characterized in that: The corresponding calculation method of the sub-index is used to calculate the ecological consumption index CNP, land stress index LSI, vegetation destruction index VDI, biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP; and the ecological stress index CI and ecological adaptability index AI are calculated respectively using different weights; at the same time, the calculated ecological stress index CI and ecological adaptability index AI are used as the two coordinate axes of the windmill model to construct a triangular model, as follows: Calculate the ecological consumption index CNP, land stress index LSI, vegetation destruction index VDI, biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP. The steps include: The ecological consumption index CNP is represented by the normalized ecological consumption, and the land stress index LSI is calculated according to the following formula: In the formula, Aero is the normalized coefficient of land stress index, SE is the area of ​​slight erosion, ME is the area of ​​moderate erosion, LE is the area of ​​severe erosion, and AS is the area of ​​calculation area; Calculate the vegetation destruction index VDI, the calculation formula is as follows: VDI=Aveg(0.5FD+0.3GD+0.2FA) / AS In the formula, Aveg is the normalized coefficient of vegetation destruction index, FD is the forest degradation area, GD is the grassland degradation area, FA is the farmland area, and AS is the calculation area; Calculate the biodiversity index BI, the calculation formula is as follows: Where N is the total number of individuals; n i is the number of individuals in the ith population; The normalized ecological supply quantity is used to represent the ecological supply index SNP; the ecological carrying index EBI is calculated according to the following formula; Where Pa is the total population of the measurement area, CNPP ave is the per capita ecological consumption, SNP is the ecological supply index; The ecological stress index CI and ecological adaptability index AI are calculated using different weights. The steps include: Different weights are assigned according to the internationally accepted urbanization level classification, as follows: When the urbanization rate is between 0% and 30%, the weights of the ecological consumption index CNP, land stress index LSI, and vegetation destruction index VDI are 0.60, 0.25, and 0.15 respectively, and the weights of the biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP are 0.5, 0.3, and 0.2 respectively; When the urbanization rate is between 30% and 70%, the weights of the ecological consumption index CNP, land stress index LSI, and vegetation destruction index VDI are 0.5, 0.3, and 0.2 respectively, and the weights of the biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP are 0.4, 0.3, and 0.3 respectively; When the urbanization rate is between 70% and 100%, the weights of the ecological consumption index CNP, land stress index LSI, and vegetation destruction index VDI are 0.3, 0.6, and 0.1 respectively, and the weights of the biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP are 0.2, 0.55, and 0.25 respectively; According to the corresponding weights, the ecological stress index CI and ecological adaptability index AI are calculated; CI=W CNP ·CNP+W LSI ·LSI+W VDI ·VDI AI=W BI ·BI+W EBI ·EBI+W SNP ·SNP Where W CNP , W LSI , W VDI are the weights of ecological consumption index CNP, land stress index LSI, vegetation destruction index VDI, and W BI , W EBI , W SNP They are the weights of the biodiversity index BI, ecological carrying index EBI, and ecological supply index SNP; Construct a windmill triangle model, take point O as the origin of the windmill single-leaf triangle model, take 45° as the angle, take the ecological stress index CI and the ecological adaptability index AI as the two sides of the angle, form a third windmill single-leaf triangle model, and calculate the area S of the third windmill single-leaf triangle model △ oC1C1', the calculation formula is as follows: The ecological adaptability index EAI of the area to be measured is obtained by calculating the ratio of the area of ​​the third windmill single-blade triangle model to the area of ​​the third standard triangle model with an angle of 1.

6. The method for calculating the comprehensive index of regional resource and environmental carrying capacity according to claim 1, characterized in that: The social unemployment index URI, the population mortality index PMI, the Engel index EGI, the economic innovation index EII, the green economy index GEI, and the sustainable economy index SEI of the tested area are calculated based on the basic calculations; and a three-dimensional pyramid model is constructed using the social unemployment index URI, the population mortality index PMI, and the Engel index EGI, and the economic and social development recession index WI is calculated based on the volume ratio of the three-dimensional pyramid to the standard three-dimensional pyramid; a three-dimensional cylinder-cone model is constructed using the economic innovation index EII, the green economy index GEI, and the sustainable economy index SEI, and the economic and social development promotion index EI is calculated based on the volume ratio of the three-dimensional cylinder-cone to the standard three-dimensional cylinder; at the same time, a triangular model is constructed using the calculated economic and social development recession index WI and the economic and social development promotion index EI as the two coordinate axes of the windmill model, as follows: The steps of representing and calculating the social unemployment index URI, the population mortality index PMI, the Engel index EGI, the economic innovation index EII, the green economy index GEI, and the sustainable economy index SEI include: The unemployment rate represents the social unemployment index URI, the mortality rate represents the population mortality index PMI, the Engel coefficient represents the Engel index EGI, and the economic innovation index EII is calculated according to the following formula; EII=0.5NPG+0.5EPD R&D Where NPG is the normalized regional patent authorization data, EPD is R&D Normalize R&D expenditure data for the region; Calculate the green economy index GEI using the following formula: GEI = 0.65EPA + 0.35EPFE In the formula, EPA is the normalized number of regional environmental protection agencies, and EPFE is the normalized regional environmental protection fiscal expenditure data; Calculate the Sustainable Economy Index SEI using the following formula: SIX=(OV eti +E i +PC GDP ) / 3,And t =(PSE+SSE+USE) / 3 In the formula, OV eti is the normalized high-tech industry output value, E i For the education index, PC GDP is the normalized GDP per capita, PSE, SSE, and USE are primary school, secondary school, and university enrollment rates, respectively; A three-dimensional pyramid model is constructed, wherein the origin of the three-dimensional space is a vertex of the bottom surface of the three-dimensional pyramid model, the distances from the origin of the three-dimensional space to the bottom surface in two-dimensional directions are respectively the social unemployment index URI and the population death index PMI of the area to be measured, and the Engel index EGI is the height of the pyramid. The pyramid model is constructed, and the volume of the pyramid model is obtained; the economic and social development recession force index WI of the area to be measured is the ratio of the volume of the pyramid model to the volume of a standard pyramid whose bottom side length and height are both 1; Construct a three-dimensional cylinder-cone model, select any point as the center point of the bottom surface of the three-dimensional cylinder-cone, use the economic innovation index EII as the radius of the lower bottom surface of the cylinder-cone, use the green economic index GEI as the radius of the upper bottom surface of the cylinder, and use the sustainable economic index SEI as the height of the cylinder-cone to construct a three-dimensional cylinder-cone model, and obtain the volume of the three-dimensional cylinder-cone model; the economic and social development promotion index EI of the area to be measured is the ratio of the volume of the three-dimensional cylinder-cone model to the volume of a standard three-dimensional cylinder model with upper and lower bottom radii and a height of 1; Construct a windmill triangle model, take point O as the origin of the windmill single-leaf triangle model, take 45° as the included angle, take the economic and social development recession force index WI and the economic and social development promotion force index EI as the two sides of the included angle, form a fourth windmill single-leaf triangle model, and calculate the area of ​​the fourth windmill single-leaf triangle model The calculation formula is: The socio-economic development index SDI of the area to be measured is obtained by comparing the area of ​​the fourth windmill single-blade triangle model to the area of ​​the standard triangle model with an angle of 1.

7. The method for calculating the comprehensive index of regional resource and environmental carrying capacity according to claim 1, characterized in that: A quadrilateral model is constructed based on the calculated resource carrying capacity index RCI, environmental carrying capacity index ECI, ecological adaptability index EAI, and socio-economic development index SDI. The resource and environmental carrying capacity comprehensive index REBCI of the tested area is obtained based on the area ratio of the quadrilateral model W1X1Y1Z1 to the standard quadrilateral WXYZ, as follows: Constructing a quadrilateral model W1X1Y1Z1, including: taking the origin of the plane rectangular coordinate system as the center point of the model, the distances from the origin of the plane rectangular coordinate system to the four directions are the resource carrying capacity index RCI, the environmental carrying capacity index ECI, the ecological adaptability index EAI, and the socio-economic development index SDI of the area to be measured to construct a quadrilateral model, and calculating the area of ​​the quadrilateral model The calculation formula is: The resource and environmental carrying capacity comprehensive index REBCI of the measured area is the ratio of the area of ​​the quadrilateral model to the area of ​​a standard quadrilateral with a distance from the origin to four directions of 1.

Images