Urban territorial space expansion integrated regulation and control method and system

Abstract

The invention discloses a town territorial space expansion integrated regulation and control method and system, and the method comprises the steps: carrying out the calculation of the comprehensive value of farmland resources of a region to be regulated and controlled based on the public information of a government, and carrying out the correction based on resource scarcity; the agricultural land resource comprehensive value comprises an agricultural land resource economic output value, an agricultural land resource ecological service value and an agricultural land resource social security value; based on an interpolation method and a statistical test method, establishing a town space land price database; calculating marginal income and green marginal income of urban territorial space expansion based on carbon source space negative extrinsic property; regulating and delimiting a town space expansion boundary based on the corrected farmland resource comprehensive value, the town space low-price database, the marginal income of town territorial space expansion and the green marginal income; the invention provides the urban territorial space expansion integrated regulation and control method and system capable of considering land price balance and low-carbon constraint.

Description

Technical Field

[0001] The present invention relates to an integrated control method and system, and in particular to an integrated control method and system for urban land space expansion, belonging to the technical field of land space planning and urban planning. Background Art

[0002] The rapid expansion of urban land space has met the needs of economic construction, but it also threatens food security and ecological security. Although the government has adopted a variety of macro-control strategies to limit the disorderly expansion of urban space, the situation of urban space expansion remains severe. At present, the methods of identifying the current boundaries of urban space through multi-spectral remote sensing data and night light remote sensing data are relatively mature, but the demarcation of urban development boundaries is mostly based on 1.3 times the existing urban space scale, which is difficult to meet the actual development needs of different regions. The land price equilibrium model is an important method for regulating the expansion of urban land space, that is, a method to find the spatial equilibrium position of urban space land prices and non-urban space land prices, so as to determine the boundary that meets the optimal comprehensive benefits of urban space expansion.

[0003] However, in the past, the delineation of urban development boundaries did not take into account the negative ecological externalities brought about by urban space as a carbon source space, resulting in an imbalance in the carbon balance of the national land space. Therefore, the ideas for regulating the expansion of urban land space and the methods for delineating urban development boundaries that take into account both food security and low-carbon development goals are not yet mature. It is urgent to propose an integrated regulatory model and method for urban land space expansion that takes into account both land price equilibrium and low-carbon constraints. Summary of the invention

[0004] Purpose of the invention: The purpose of the present invention is to provide an integrated control method and system for urban land space expansion that can take into account land price equilibrium and low-carbon constraints.

[0005] Technical solution: The integrated control method for urban land space expansion described in the present invention comprises the following steps:

[0006] Step 1: Calculate the comprehensive value of farmland resources in the area to be regulated based on government public information and make corrections based on resource scarcity; the comprehensive value of farmland resources includes the economic output value of farmland resources, the ecological service value of farmland resources and the social security value of farmland resources;

[0007] Step 2: Based on the interpolation method and statistical test method, establish the urban spatial land price database;

[0008] Step 3: Based on the negative externalities of carbon source space, calculate the marginal benefits and green marginal benefits of urban land space expansion;

[0009] Step 4: Based on the revised comprehensive value of agricultural land resources, the urban space low-price database, the marginal benefits of urban land space expansion and the green marginal benefits, regulate and define the boundaries of urban space expansion.

[0010] Furthermore, the step 2 is specifically as follows: using the benchmark land price map of the urban research year as the working base map for digitization, and establishing an urban spatial graphic database; collecting information on the plots of land sold in the primary market of the town, locating and establishing an attribute database, and based on the interpolation method and statistical test method, considering the representativeness of the sample point type, spatially interpolating the urban spatial land value, and forming an urban spatial land price database.

[0011] Furthermore, the step 3 comprises:

[0012] Step 3.1: Collect the annual average carbon trading price of the national carbon market in the study year P C ;

[0013] Step 3.2: Calculate the marginal benefit MP of urban spatial expansion GDP :

[0014] MP GDP =d∑GDP / ds

[0015] Among them, s represents the urban spatial scale; ΣGDP is the sum of economic output within the urban spatial scope;

[0016] Step 3.3: Calculate the marginal benefit MP′ considering urban spatial expansion GDP ; Based on the negative externalities of urban space, deducting the loss of ecological value caused by increased carbon emissions due to urban expansion:

[0017] MP′ GDP =d{∑GDP-∑CO 2 ×P C} / ds

[0018] Among them, ΣCO 2 It is the cumulative carbon emissions within the urban space.

[0019] Furthermore, the step 4 comprises:

[0020] Step 4.1: Based on the urban spatial land price database, the Ordinary Kring method is used in the ArcGIS Geostatistical Analysis module to complete the spatial interpolation of urban spatial land price sample points;

[0021] Step 4.2: The preset multiples of the urban spatial scale based on the planning starting year are used as the ideal boundary;

[0022] Step 4.3: Constructing an ideal efficiency model for urban spatial expansion1 , Moderate Efficiency Model JP 2 , Efficiency Limit Model JP 3 , Carbon Reduction Ideal Model JP 4 , Carbon Reduction Moderation Model JP 5 and the Carbon Reduction Limit Model JP 6 as follows:

[0023] JP 1 =MP GDP =(V′ e +V′ c +V′ s )

[0024] JP 2 =MP GDP =(V′ e +V′ c )

[0025] JP 3 =MP GDP =V′ e

[0026]

[0027] Among them, MP GDP is the marginal benefit of town spatial expansion, MP′ GDP is the green marginal benefit of town space expansion, V e ′ represents the corrected economic value of agricultural land resources, V c ′ represents the corrected ecological value of farmland resources, V s ′ represents the corrected social value of agricultural land resources, s is the urban space scale, ∑CO 2 is the cumulative carbon emissions within the urban space, P C is the annual average price of carbon trading in the national carbon market;

[0028] Step 4.4: According to the calculation results of the ideal boundary model, select the construction land value as JP in the spatial database of urban construction land value after spatial interpolation. 1 , JP 2 , JP 3 , JP 4 , JP 5 and JP 6 The six closed curves are used as the ideal efficiency improvement boundary, moderate efficiency improvement boundary, extreme efficiency improvement boundary, ideal carbon reduction boundary, moderate carbon reduction boundary and extreme carbon reduction boundary to complete the regulation of the urban space expansion boundary.

[0029] Furthermore, the step 4.2 is specifically as follows:

[0030] Step 4.2.1: Based on the multi-objective optimization model, on the basis of determining that the urban space scale is n times the current urban space scale, and satisfying the constraints that the agricultural land area is not less than the permanent basic farmland red line area set by the government, and the ecological land area is not less than the permanent ecological protection red line area set by the government, the optimization goal is to maximize the economic and ecological benefits of land use, and obtain the quantitative structure of various types of urban space, where n is the preset multiple;

[0031] Step 4.2.2: Using the PLUS model, based on the land use change data of the two periods, the LEAS module is used to capture the historical land use change strategy, and then the CARS module is used to simulate the urban spatial layout after the urban space has expanded n times;

[0032] Step 4.2.3: Extract the urban spatial boundary in the simulation results as the ideal boundary for urban land space expansion.

[0033] Based on the same inventive concept, the present invention also provides an integrated control system for urban land space expansion, comprising:

[0034] The farmland calculation module is used to calculate the comprehensive value of farmland resources in the area to be regulated based on government public information and make corrections based on resource scarcity; the comprehensive value of farmland resources includes the economic output value of farmland resources, the ecological service value of farmland resources and the social security value of farmland resources;

[0035] Database module: used to establish urban spatial land price database based on interpolation method and statistical test method;

[0036] The benefit calculation module is used to calculate the marginal benefits and green marginal benefits of urban land space expansion based on the negative externalities of carbon source space;

[0037] The control module is used to control and delineate the boundaries of urban space expansion based on the revised comprehensive value of agricultural land resources, the urban space low-price database, the marginal benefits of urban land space expansion, and the green marginal benefits.

[0038] Furthermore, the specific implementation method of the database module is to digitize the benchmark land price map of the urban research year as the working base map to establish an urban spatial graphic database; collect information on the land parcels sold in the primary market of the town, locate and establish an attribute database, and based on the interpolation method and statistical test method, consider the representativeness of the sample point type, perform spatial interpolation on the urban spatial land value to form an urban spatial land price database.

[0039] Furthermore, the revenue toilet module is specifically implemented as follows:

[0040] Collect the annual average carbon trading price of the national carbon market in the research year P C ;

[0041] Calculating the marginal benefits of urban spatial expansion GDP :

[0042] MP GDP =d∑GDP / ds

[0043] Among them, s represents the urban spatial scale; ∑CDP is the sum of economic output within the urban spatial scope;

[0044] Calculate the marginal benefit MP′ considering urban spatial expansion GDP ; Based on the negative externalities of urban space, deducting the loss of ecological value caused by increased carbon emissions due to urban expansion:

[0045] MP′ GDP =d{∑GDP-∑CO 2 ×P C} / ds

[0046] Among them, ∑CO 2 It is the cumulative carbon emissions within the urban space.

[0047] Furthermore, the control module is specifically implemented as follows:

[0048] Based on the urban spatial land price database, the Ordinary Kring method is used in the ArcGIS Geostatistical Analysis module to complete the spatial interpolation of urban spatial land price sample points;

[0049] The preset multiples of the urban spatial scale based on the planning starting year are the ideal boundaries;

[0050] Constructing an ideal model for increasing the efficiency of urban spatial expansion JP 1 , Moderate Efficiency Model JP 2 , Efficiency Limit Model JP 3 , Carbon Reduction Ideal Model JP 4 , Carbon Reduction Moderation Model JP 5 and the Carbon Reduction Limit Model JP 6 as follows:

[0051] JP 1 =MP GDP =(V′ e +V′ c +V′ s )

[0052] JP 2 =MP GDP =(V′ e +V′ c )

[0053] JP3 =MP GDP =V′ e

[0054]

[0055] Among them, MP GDP is the marginal benefit of town spatial expansion, MP′ GDP is the green marginal benefit of town spatial expansion, V′ e represents the corrected economic value of agricultural land resources, V c ′ represents the corrected ecological value of farmland resources, V s ′ represents the corrected social value of agricultural land resources, s is the urban space scale, ∑CO 2 is the cumulative carbon emissions within the urban space, P C The annual average price of carbon trading in the national carbon market;

[0056] According to the calculation results of the ideal boundary model, the construction land value of JP is selected in the spatial database of urban construction land value after spatial interpolation. 1 , JP 2 , JP 3 , JP 4 , JP 5 and JP 6 The six closed curves are used as the ideal efficiency improvement boundary, moderate efficiency improvement boundary, extreme efficiency improvement boundary, ideal carbon reduction boundary, moderate carbon reduction boundary and extreme carbon reduction boundary to complete the regulation of the urban space expansion boundary.

[0057] Furthermore, the preset multiples of the urban space scale based on the planning starting year are the ideal boundaries, specifically including:

[0058] Based on the multi-objective optimization model, on the basis of determining that the urban space scale is n times the current urban space scale, and satisfying the constraints that the agricultural land area is not less than the permanent basic farmland red line area set by the government, and the ecological land area is not less than the permanent ecological protection red line area set by the government, the optimization goal is to maximize the economic and ecological benefits of land use, and the quantitative structure of various types of urban space is obtained, where n is a preset multiple;

[0059] Using the PLUS model, based on the land use change data of two periods, the LEAS module is used to capture the historical land use change strategy, and then the CARS module is used to simulate the urban spatial layout after the urban space has expanded n times;

[0060] The urban spatial boundaries in the simulation results are extracted as the ideal boundaries for urban territorial spatial expansion.

[0061] Beneficial effects: Compared with the existing technology, the integrated regulation and control model and method of urban land space expansion provided by the present invention comprehensively balances multiple factors such as land space and economic efficiency, social security and dual carbon goals, and focuses on the negative externalities of urban space as a carbon source space. The newly proposed 6 boundaries theoretically define the spatial indifference curves of land resources as agricultural land and urban construction land, and in practice find a spatial reference system for multiple integrated regulation and control schemes of urban land space expansion for economic efficiency and dual carbon goals; the present invention also integrates and constructs an integrated three-dimensional model of "strategy-mechanism-governance-function-background" from a multidisciplinary perspective, and proposes an integrated regulation and control model and method for urban land space expansion from multiple schemes (planning, efficiency improvement and carbon reduction), which can provide regulation and control schemes for regional land space governance and high-quality development. BRIEF DESCRIPTION OF THE DRAWINGS

[0062] Figure 1 The figure is a flow chart of a method according to an embodiment of the present invention. DETAILED DESCRIPTION

[0063] The present invention will be further described in detail below in conjunction with the accompanying drawings and examples of the present invention.

[0064] As attached Figure 1 As shown, this embodiment provides an integrated control method for urban land space expansion, which is specifically implemented according to the following steps:

[0065] Step 1: Calculate the comprehensive value of farmland resources in the area to be regulated based on government public information and make corrections based on resource scarcity; the comprehensive value of farmland resources includes the economic output value of farmland resources, the ecological service value of farmland resources and the social security value of farmland resources;

[0066] Step 2: Based on the interpolation method and statistical test method, establish the urban spatial land price database;

[0067] Step 3: Based on the negative externalities of carbon source space, calculate the marginal benefits and green marginal benefits of urban land space expansion;

[0068] Step 4: Based on the revised comprehensive value of agricultural land resources, the urban space low-price database, the marginal benefits of urban land space expansion and the green marginal benefits, regulate and define the boundaries of urban space expansion.

[0069] Specifically, step 2 is to digitize the benchmark land price map of the urban research year as the working base map to establish an urban spatial graphic database; collect information on the land parcels sold in the urban primary market, locate and establish an attribute database, and based on the interpolation method and statistical test method, consider the representativeness of the sample point type, perform spatial interpolation on the urban spatial land value to form an urban spatial land price database.

[0070] The step 3 comprises:

[0071] Step 3.1: Collect the annual average carbon trading price of the national carbon market in the study year P C ;

[0072] Step 3.2: Calculate the marginal benefit MP of urban spatial expansion GDP :

[0073] MP GDP =d∑GDP / ds

[0074] Among them, s represents the urban space scale; ∑GDP is the sum of economic output within the urban space; the urban space scale is determined by using ArcGIS to calculate the urban space scale of the study area based on the Chinese multi-period land use / land cover remote sensing detection dataset of the Resource and Environmental Science Data Center of the Chinese Academy of Sciences. The land types included in the urban space are urban and industrial and transportation construction land.

[0075] Step 3.3: Calculate the marginal benefit MP′ considering urban spatial expansion GDP ; Based on the negative externalities of urban space, deducting the loss of ecological value caused by increased carbon emissions due to urban expansion:

[0076] MP′ GDP =d{∑GDP-∑CO 2 ×P C} / ds

[0077] Among them, ∑CO 2 is the cumulative carbon emissions within the scope of urban space. The calculation method of cumulative carbon emissions within the scope of urban space is as follows: based on the carbon emission data of the province where the study area is located (derived from the China Carbon Accounting Database CAEDs) and the night light remote sensing data, a relationship model between the carbon emission data of the province in the study area and the night light brightness value is established. Generally, a linear model can be used. On the basis of good model simulation effect, the carbon emission value of each grid within the study area is calculated through the night light brightness value of the study area and the relationship model between the obtained carbon emission data and the night light brightness value, and then the cumulative carbon emissions of the urban space are obtained through ArcGIS spatial statistics.

[0078] The step 4 comprises:

[0079] Step 4.1: Based on the urban spatial land price database, the Ordinary Kring method is used in the ArcGIS Geostatistical Analysis module to complete the spatial interpolation of urban spatial land price sample points;

[0080] Step 4.2: The preset multiples of the urban spatial scale based on the planning starting year are used as the ideal boundary;

[0081] Step 4.2.1: Based on the multi-objective optimization model, on the basis of determining that the urban space scale is n times the current urban space scale, and satisfying the constraints that the agricultural land area is not less than the permanent basic farmland red line area set by the government, and the ecological land area is not less than the permanent ecological protection red line area set by the government, the optimization goal is to maximize the economic and ecological benefits of land use, and obtain the quantitative structure of various types of urban space, where n is the preset multiple;

[0082] Step 4.2.2: Using the PLUS model, based on the land use change data of the two periods, the LEAS module is used to capture the historical land use change strategy, and then the CARS module is used to simulate the urban spatial layout after the urban space has expanded n times;

[0083] Step 4.2.3: Extract the urban spatial boundary in the simulation results as the ideal boundary for urban land space expansion.

[0084] Step 4.3: Constructing an ideal efficiency model for urban spatial expansion 1 , Moderate Efficiency Model JP 2 , Efficiency Limit Model JP 3 , Carbon Reduction Ideal Model JP 4 , Carbon Reduction Moderation Model JP 5 and the Carbon Reduction Limit Model JP 6 as follows:

[0085] JP 1 =MP GDP =(V′ e +V′ c +V′ s )

[0086] JP 2 =MP GDP =(V′ c +V′ c )

[0087] JP 3 =MP GDP =V′ e

[0088]

[0089] Among them, MP GDP is the marginal benefit of town spatial expansion, MP′ GDP is the green marginal benefit of town space expansion, V e ′ represents the corrected economic value of agricultural land resources, V c ′ represents the corrected ecological value of farmland resources, V s′ represents the corrected social value of agricultural land resources, s is the urban space scale, ∑CO 2 is the cumulative carbon emissions within the urban space, P C is the annual average price of carbon trading in the national carbon market;

[0090] Land prices can reflect the utility of space. When the utility of urban space is equal to the utility of agricultural land, the urban land price should be equal to the agricultural land price. The urban development boundary found based on this land price equilibrium is conducive to the coordination of urban expansion and farmland protection, and realizes the optimal comprehensive benefits of national space.

[0091] However, urban land prices do not take into account the negative externalities of urban space as a carbon source space. Therefore, the high land prices reflect the utility of urban space and ignore the negative externalities of urban space on the ecological environment.

[0092] Among them, r is the reduction interest rate, which refers to the bank loan interest rate. Because the annual carbon emissions brought by urban space as a carbon source space are continuous, the negative impact of carbon emissions on the ecological environment is continuous.

[0093] Step 4.4: According to the calculation results of the ideal boundary model, select the construction land value as JP in the spatial database of urban construction land value after spatial interpolation. 1 , JP 2 , JP 3 , JP 4 , JP 5 and JP 6 The six closed curves are used as the ideal efficiency improvement boundary, moderate efficiency improvement boundary, extreme efficiency improvement boundary, ideal carbon reduction boundary, moderate carbon reduction boundary and extreme carbon reduction boundary to complete the regulation of the urban space expansion boundary.

[0094] Based on the same inventive concept, this embodiment also provides an integrated control system for urban land space expansion, including:

[0095] The farmland calculation module is used to calculate the comprehensive value of farmland resources in the area to be regulated based on government public information and make corrections based on resource scarcity; the comprehensive value of farmland resources includes the economic output value of farmland resources, the ecological service value of farmland resources and the social security value of farmland resources;

[0096] Database module: used to establish urban spatial land price database based on interpolation method and statistical test method;

[0097] The benefit calculation module is used to calculate the marginal benefits and green marginal benefits of urban land space expansion based on the negative externalities of carbon source space;

[0098] The control module is used to control and delineate the boundaries of urban space expansion based on the revised comprehensive value of agricultural land resources, the urban space low-price database, the marginal benefits of urban land space expansion, and the green marginal benefits.

[0099] Furthermore, the specific implementation method of the database module is to digitize the benchmark land price map of the urban research year as the working base map to establish an urban spatial graphic database; collect information on the land parcels sold in the primary market of the town, locate and establish an attribute database, and based on the interpolation method and statistical test method, consider the representativeness of the sample point type, perform spatial interpolation on the urban spatial land value to form an urban spatial land price database.

[0100] Furthermore, the revenue toilet module is specifically implemented as follows:

[0101] Collect the annual average carbon trading price of the national carbon market in the research year P C ;

[0102] Calculating the marginal benefits of urban spatial expansion GDP :

[0103] MP GDP =d∑GDP / ds

[0104] Among them, s represents the urban spatial scale; ∑GDP is the sum of economic output within the urban spatial scope;

[0105] Calculate the marginal benefit MP′ considering urban spatial expansion GDP ; Based on the negative externalities of urban space, deducting the loss of ecological value caused by increased carbon emissions due to urban expansion:

[0106] MP GDP =d{ΣGDP-ΣCO 2 ×P C} / ds

[0107] Among them, ∑CO 2 It is the cumulative carbon emissions within the urban space.

[0108] Furthermore, the control module is specifically implemented as follows:

[0109] Based on the urban spatial land price database, the Ordinary Kring method is used in the ArcGIS Geostatistical Analysis module to complete the spatial interpolation of urban spatial land price sample points;

[0110] The preset multiples of the urban spatial scale based on the planning starting year are the ideal boundaries;

[0111] Constructing an ideal model for increasing the efficiency of urban spatial expansion JP 1, Moderate Efficiency Model JP 2 , Efficiency Limit Model JP 3 , Carbon Reduction Ideal Model JP 4 , Carbon Reduction Moderation Model JP 5 and the Carbon Reduction Limit Model JP 6 as follows:

[0112] JP 1 =MP GDP =(V′ e +V′ c +V′ s )

[0113] JP 2 =MP GDP =(V′ e +V′ c )

[0114] JP 3 =MP GDP =V′ e

[0115]

[0116] Among them, MP GDP is the marginal benefit of town spatial expansion, MP′ GDP is the green marginal benefit of town spatial expansion, V′ e represents the corrected economic value of agricultural land resources, V c ′ represents the corrected ecological value of farmland resources, V s ′ represents the corrected social value of agricultural land resources, s is the urban space scale, ∑CO 2 is the cumulative carbon emissions within the urban space, P C The annual average price of carbon trading in the national carbon market;

[0117] According to the calculation results of the ideal boundary model, the construction land value of JP is selected in the spatial database of urban construction land value after spatial interpolation. 1 , JP 2 , JP 3 , JP 4 , JP 5 and JP 6 The six closed curves are used as the ideal efficiency improvement boundary, moderate efficiency improvement boundary, extreme efficiency improvement boundary, ideal carbon reduction boundary, moderate carbon reduction boundary and extreme carbon reduction boundary to complete the regulation of the urban space expansion boundary.

[0118] Furthermore, the preset multiples of the urban space scale based on the planning starting year are the ideal boundaries, specifically including:

[0119] Based on the multi-objective optimization model, on the basis of determining that the urban space scale is n times the current urban space scale, and satisfying the constraints that the agricultural land area is not less than the permanent basic farmland red line area set by the government, and the ecological land area is not less than the permanent ecological protection red line area set by the government, the optimization goal is to maximize the economic and ecological benefits of land use, and the quantitative structure of various types of urban space is obtained, where n is a preset multiple;

[0120] Using the PLUS model, based on the land use change data of two periods, the LEAS module is used to capture the historical land use change strategy, and then the CARS module is used to simulate the urban spatial layout after the urban space has expanded n times;

[0121] The urban spatial boundaries in the simulation results are extracted as the ideal boundaries for urban territorial spatial expansion.

Claims

1. A method for integrated control of urban land space expansion, characterized in that: The following steps are involved: Step 1: Calculate the comprehensive value of farmland resources in the area to be regulated based on government public information and make corrections based on resource scarcity; the comprehensive value of farmland resources includes the economic output value of farmland resources, the ecological service value of farmland resources and the social security value of farmland resources; Step 2: Based on the interpolation method and statistical test method, establish the urban spatial land price database; Step 3: Based on the negative externalities of carbon source space, calculate the marginal benefits and green marginal benefits of urban land space expansion; Step 4: Based on the revised comprehensive value of agricultural land resources, the urban space low-price database, the marginal benefits of urban land space expansion and the green marginal benefits, regulate and define the boundaries of urban space expansion.

2. The integrated control method for urban land space expansion according to claim 1 is characterized in that: Specifically, the step 2 is to digitize the benchmark land price map of the urban research year as the working base map to establish an urban spatial graphic database; collect information on the land parcels sold in the primary market of the urban area, locate and establish an attribute database, and based on the interpolation method and statistical test method, consider the representativeness of the sample point type, perform spatial interpolation on the urban spatial land value to form an urban spatial land price database.

3. The integrated control method for urban land space expansion according to claim 2 is characterized in that: The step 3 comprises: Step 3.1: Collect the annual average carbon trading price of the national carbon market in the study year P C ; Step 3.2: Calculate the marginal benefit MP of urban spatial expansion GDP : MP GDP =d∑GDP / ds Among them, s represents the urban spatial scale; ∑GDP is the sum of economic output within the urban spatial scope; Step 3.3: Calculate the marginal benefit MP′ considering urban spatial expansion GDP ; Based on the negative externalities of urban space, deducting the loss of ecological value caused by increased carbon emissions due to urban expansion: MP′ GDP =d(∑GDP-∑CO2×P c } / ds Among them, ∑CO2 is the cumulative carbon emissions within the urban space.

4. The integrated control method for urban land space expansion according to claim 3 is characterized in that: The step 4 comprises: Step 4.1: Based on the urban spatial land price database, the Ordinary Kring method is used in the ArcGIS Geostatistical Analysis module to complete the spatial interpolation of urban spatial land price sample points; Step 4.2: The preset multiples of the urban spatial scale based on the planning starting year are used as the ideal boundary; Step 4.3: Construct the ideal efficiency model JP1, the moderate efficiency model JP2, the extreme efficiency model JP3, the ideal carbon reduction model JP4, the moderate carbon reduction model JP5 and the extreme carbon reduction model JP6 of urban spatial expansion as follows: JP1=MP GDP =(V′ c +V′ c +V′ s ) JP2=MP GDP =(V′ c +V′ c ) JP3=MP GDP =V′ c Among them, MP GDP is the marginal benefit of town spatial expansion, MP′ GDP is the green marginal benefit of town spatial expansion, V′ e represents the corrected economic value of agricultural land resources, V′ c represents the corrected ecological value of farmland resources, V′ s represents the corrected social value of agricultural land resources, s is the urban space scale, ∑CO2 is the cumulative carbon emissions within the urban space, and P C The annual average price of carbon trading in the national carbon market; Step 4.4: According to the calculation results of the ideal boundary model, six closed curves with construction land values ​​of JP1, JP2, JP3, JP4, JP5 and JP6 are selected from the spatial database of urban construction land value after spatial interpolation as the ideal efficiency-enhancing boundary, moderate efficiency-enhancing boundary, extreme efficiency-enhancing boundary, ideal carbon reduction boundary, moderate carbon reduction boundary and extreme carbon reduction boundary to complete the regulation of the urban spatial expansion boundary.

5. The integrated control method for urban land space expansion according to claim 4 is characterized in that: The step 4.2 is specifically as follows: Step 4.2.1: Based on the multi-objective optimization model, on the basis of determining that the urban space scale is n times the current urban space scale, and satisfying the constraints that the agricultural land area is not less than the permanent basic farmland red line area set by the government, and the ecological land area is not less than the permanent ecological protection red line area set by the government, the optimization goal is to maximize the economic and ecological benefits of land use, and obtain the quantitative structure of various types of urban space, where n is the preset multiple; Step 4.2.2: Using the PLUS model, based on the land use change data of the two periods, the LEAS module is used to capture the historical land use change strategy, and then the CARS module is used to simulate the urban spatial layout after the urban space has expanded n times; Step 4.2.3: Extract the urban spatial boundary in the simulation results as the ideal boundary for urban land space expansion.

6. An integrated control system for urban land space expansion, characterized in that: include: The farmland calculation module is used to calculate the comprehensive value of farmland resources in the area to be regulated based on government public information and make corrections based on resource scarcity; the comprehensive value of farmland resources includes the economic output value of farmland resources, the ecological service value of farmland resources and the social security value of farmland resources; Database module: used to establish urban spatial land price database based on interpolation method and statistical test method; The benefit calculation module is used to calculate the marginal benefits and green marginal benefits of urban land space expansion based on the negative externalities of carbon source space; The control module is used to control and delineate the boundaries of urban space expansion based on the revised comprehensive value of agricultural land resources, the urban space low-price database, the marginal benefits of urban land space expansion, and the green marginal benefits.

7. The integrated control system for urban land space expansion according to claim 6 is characterized in that: The specific implementation method of the database module is to digitize the benchmark land price map of the urban research year as the working base map to establish an urban spatial graphic database; collect information on the land parcels sold in the primary market of the urban area, locate and establish an attribute database, and based on the interpolation method and statistical test method, consider the representativeness of the sample point type, perform spatial interpolation on the urban spatial land value to form an urban spatial land price database.

8. The integrated control system for urban land space expansion according to claim 7 is characterized in that: The specific implementation method of the revenue toilet module is as follows: Collect the annual average carbon trading price of the national carbon market in the research year P C ; Calculating the marginal benefits of urban spatial expansion GDP : MP GDP =d∑GDP / ds Among them, s represents the urban spatial scale; ∑GDP is the sum of economic output within the urban spatial scope; Calculate the marginal benefit MP′ considering urban spatial expansion GDP ; Based on the negative externalities of urban space, deducting the loss of ecological value caused by increased carbon emissions due to urban expansion: MP′ GDP =d{∑GDP-∑CO2×P C } / ds Among them, ∑CO2 is the cumulative carbon emissions within the urban space.

9. The integrated control system for urban land space expansion according to claim 8 is characterized in that: The control module is specifically implemented as follows: Based on the urban spatial land price database, the Ordinarykring method is used in the ArcGIS Geostatistical Analysis module to complete the spatial interpolation of urban spatial land price sample points; The preset multiples of the urban spatial scale based on the planning starting year are the ideal boundaries; The ideal efficiency model JP1, moderate efficiency model JP2, extreme efficiency model JP3, ideal carbon reduction model JP4, moderate carbon reduction model JP5 and extreme carbon reduction model JP6 of urban spatial expansion are constructed as follows: JP1=MP GDP =(V′ e +V′ c +V′ s ) JP2=MP GDP =(V′ e +V′ c ) JP3=MP GDP =V′ e Among them, MP′ GDP is the marginal benefit of town spatial expansion, MP′ GDP is the green marginal benefit of town spatial expansion, V′ e represents the corrected economic value of agricultural land resources, V′ c represents the corrected ecological value of farmland resources, V′ s represents the corrected social value of agricultural land resources, s is the urban space scale, ∑CO2 is the cumulative carbon emissions within the urban space, and P C is the annual average price of carbon trading in the national carbon market; According to the calculation results of the ideal boundary model, six closed curves with construction land values ​​of JP1, JP2, JP3, JP4, JP5 and JP6 were selected in the spatial database of urban construction land value after spatial interpolation as the ideal efficiency-enhancing boundary, moderate efficiency-enhancing boundary, extreme efficiency-enhancing boundary, ideal carbon reduction boundary, moderate carbon reduction boundary and extreme carbon reduction boundary to complete the regulation of the urban spatial expansion boundary.

10. The integrated control system for urban land space expansion according to claim 9 is characterized in that: The preset multiples of the urban space scale based on the planning starting year are the ideal boundaries, specifically including: Based on the multi-objective optimization model, on the basis of determining that the urban space scale is n times the current urban space scale, and satisfying the constraints that the agricultural land area is not less than the permanent basic farmland red line area set by the government, and the ecological land area is not less than the permanent ecological protection red line area set by the government, the optimization goal is to maximize the economic and ecological benefits of land use, and the quantitative structure of various types of urban space is obtained, where n is a preset multiple; Using the PLUS model, based on the land use change data of two periods, the LEAS module is used to capture the historical land use change strategy, and then the CARS module is used to simulate the urban spatial layout after the urban space has expanded n times; The urban spatial boundaries in the simulation results are extracted as the ideal boundaries for urban territorial spatial expansion.

Images