CN115619095A

CN115619095A - Method and system for evaluating homeland space suitability based on optimal balancing strategy

Info

Publication number: CN115619095A
Application number: CN202211266741.5A
Authority: CN
Inventors: 欧定华; 秦景; 龚上淇; 周瑞; 高雪松; 张琪; 邓欧平; 张艺辰; 罗李佳怡
Original assignee: Sichuan Agricultural University
Current assignee: Sichuan Agricultural University
Priority date: 2022-10-17
Filing date: 2022-10-17
Publication date: 2023-01-17

Abstract

The application discloses a homeland space suitability evaluation method based on an optimal balance strategy, which comprises the steps of evaluating the bearing capacity of a resource environment; calculating the initial suitability of the homeland space function; constructing a homeland space function game model, and calculating the homeland space function game model to obtain the selection probability of the homeland space function type under the optimal balance strategy; correcting the initial suitability of the homeland space function to obtain the final suitability of the homeland space function; and dividing the suitability level of the homeland space. The method can effectively solve the problem that multiple types of high-suitability soil space conflicts appear in evaluation results of a traditional soil space suitability evaluation method, neglects the comprehensive functional benefit maximization of the soil space pattern formed by constructing the high-suitability areas, and does not consider the technical defects that the influence of interaction relation between a planner on the configuration of the soil space and a soil space functional system on the suitability stability and the like.

Description

Method and system for evaluating homeland space suitability based on optimal balancing strategy

Technical Field

The invention relates to the technical field of territory space planning, in particular to a territory space suitability evaluation method and system based on an optimal balancing strategy.

Background

In the current regional development, along with the increase of the demand of urban and rural construction land, the agricultural and ecological land is greatly reduced, the contradiction between the urban land, the agricultural land and the ecological land is aggravated, and the national soil space development pattern is disordered. Meanwhile, the conflict between scarcity of the territorial space resources and the diversity of the functions causes the competition phenomenon of the territorial space functions in the resource allocation process. Unordered competition of various functions on the homeland space resources also causes unbalance of distribution of the homeland space resources, and causes a plurality of problems such as contradiction between urbanization and grain safety, contradiction between agricultural development and environmental safety and the like. Therefore, the formation of a high-quality collaborative development mode is the key point for solving the problems faced by the development of the homeland space, and the multifunctional coordination of the homeland space is the basic content and guarantee of the homeland space. The multifunctional coordination of the territorial space is one of the core contents of space planning, and is an important way for coping with the competition of the territorial space functions, realizing the reasonable utilization of the territorial space resources and guaranteeing the sustainable development of the society. The method develops the multifunctional coordination research of the homeland space, scientifically and effectively configures the functions of the homeland space, balances the distribution and the relation among the space functions, can fully excavate the potential value of the homeland space, improves the comprehensive benefit of utilization of the homeland space, and provides theoretical basis and working basis for constructing a new pattern for developing and protecting the homeland space with obvious main functions, complementary advantages and high-quality development.

Research aiming at the evaluation of the suitability of the homeland space function is not common, and the similar evaluation of the suitability of the homeland space development is the hot content of the current research and makes an important progress. In recent years, with the rapid development of information technology, evaluation data are richer and more diverse, the operability of the evaluation method and the accuracy of results are greatly improved, and students develop rich case research by using a multi-factor comprehensive evaluation method, a space interaction and trend simulation evaluation method, an ecological niche space supply and demand coupling evaluation method and other methods aiming at different targets, wherein the multi-factor comprehensive evaluation method is most widely applied. Although the method can scientifically complete regional territorial soil space function suitability evaluation and obtain space function suitability as the basis of territorial soil space reconstruction, the suitability has the problems of high-value region overlapping, poor stability, neglecting space utilization comprehensive function benefits and the like, great troubles are brought to territorial soil space planners, and territorial soil space planning and compiling are difficult to support well.

The problems of the prior art are as follows:

when space function suitability evaluation is adopted in the existing territorial space planning process, the problems that the territorial space utilization comprehensive function benefit and the functional system stability are ignored, and a plurality of high-suitability territorial space conflicts exist in an evaluation result exist, so that the high-suitability territorial space conflicts, the territorial space suitability stability is poor, and the territorial space utilization function comprehensive benefit is low in the territorial space planning.

Disclosure of Invention

The technical problems to be solved by the application are that when space function suitability evaluation is adopted in the existing territorial space planning process, comprehensive function benefits of territorial space utilization and the stability of a functional system are ignored, and a plurality of high-suitability territorial space conflicts exist in evaluation results.

The application is realized by the following technical scheme:

the application provides a territorial space suitability evaluation method based on an optimal balance strategy in a first aspect, which comprises the steps of

1. A homeland space suitability evaluation method based on an optimal balance strategy is characterized by comprising the following steps:

s1, evaluating the bearing capacity of the resource environment, and dividing the grade of the bearing capacity of the resource environment;

s2, calculating the initial suitability of the corresponding homeland space function according to the resource environment bearing capacity evaluation result;

s3, constructing a homeland space function game model, and solving and obtaining the selection probability of the homeland space function type under the optimal balance strategy through the homeland space function game model;

s4, taking the selection probability of the territorial space function type under the optimal balance strategy as a correction coefficient, and correcting the initial suitability of the corresponding territorial space function to obtain the final suitability of the territorial space function;

and S5, dividing the suitability level of the homeland space based on the final suitability of the homeland space function.

In the technical scheme, the resource environment bearing capacity is evaluated, the resource environment bearing capacity grade is divided, and the initial suitability of the corresponding homeland space function is calculated according to the resource environment bearing capacity evaluation result. Available territorial space can be obtained preliminarily through the initial suitability, the territorial space is generally planned according to the initial suitability in the prior art, and the initial suitability at the moment often has the problems of high-value region overlapping, poor stability and neglect of space utilization comprehensive functional benefits, so that high-suitability territorial space conflict, poor stability and low territorial space utilization comprehensive functional benefits are caused during territorial space planning.

According to the method and the device, a homeland space function game model is constructed based on the initial fitness of the homeland space function, and the selection probability of the homeland space function type under the optimal balancing strategy is obtained through calculation of the homeland space function game model. And taking the selection probability of the homeland space function type under the optimal balancing strategy as a correction coefficient to correct the initial suitability of the homeland space function. The method comprises the steps of constructing a homeland space function game model, considering space game comprehensive function benefits, integrating space utilization comprehensive function benefits through selection probability of homeland space function types under an optimal balance strategy obtained through calculation of the homeland space function model, and solving the problem of neglecting space utilization comprehensive function benefits when initial fitness evaluation is used; the high-value regions of the initial suitability degree are further distinguished through the selection probability of the territorial space function types, and the problem of high-value region overlapping is solved.

The suitability grade division of the homeland space is carried out through the final suitability of the homeland space function, compared with the initial suitability grade division of the homeland space only through the initial suitability, the suitability grade division of the homeland space is more efficient and stable, and the problems of high-suitability homeland space conflict, poor stability and low comprehensive benefit of the utilization function of the homeland space during the planning of the homeland space are solved.

The territorial space function game model fundamentally solves the problems of space conflict and competition in the traditional multi-factor comprehensive evaluation method, greatly improves the rationality and the practicability of the territorial space suitability evaluation result as territorial space reconstruction, and can better provide service for territorial space planning.

In an alternative embodiment, the territorial space function game model is as follows:

in the above formula, G represents a territorial space function game model, S _l The strategic space representing the participating decision maker,

a revenue function representing the participating decision maker; wherein, the first and the second end of the pipe are connected with each other,

the set of the decision-making persons participating in the game belongs to the group of l E { a, b }, wherein in the formula, a represents a first decision-making person, and b represents a second decision-making person;

the strategic space of each decider is

In the above formula, l represents the decision maker participating in the game, l = a, b; n represents the total number of the homeland space function types; k represents a homeland space function type, k =1,2, \8230;, n;

a kth specific strategy representing the l-th decision maker;

a strategy set representing the ith decision maker; if the decision maker l selects a specific strategy, vector

Representing a strategic combination;

the revenue function for each decision maker is

In the above formula, v _k Represents k-type native space function quantity, a _k Represents the initial suitability of k-type soil space function, a _k ∈[0,1]，p _ij Representing the probability of the transition of the homeland space i to the homeland space j.

In an alternative embodiment, the variable v in the revenue function of each decision maker _k 、a _k 、p _ij The determination method of (2) is as follows:

k-type homeland space functional quantity v _k By passing

Determining; wherein (r) _k ) _m The functional quantity of the Nth sub-function of the k-type territorial space function; (omega. Of _k ) _m The Nth sub-function weight value of the k-type territorial space function;

initial suitability a of k-type homeland space function _k Obtaining the initial adaptation of the homeland space function through calculationDetermining the suitability;

probability p of transformation from homeland space i to homeland space j _ij Calculating and obtaining a transfer probability matrix p for converting the territorial space function type i into the territorial space function type j in the period from t to t +1 through a Markov model _ij And (4) determining.

In an optional embodiment, the method for obtaining the selection probability of the homeland space function type under the optimal balancing strategy comprises the following steps:

s31, constructing a territorial space utilization function comprehensive benefit matrix of the first decision maker and the second decision maker;

s32, calculating the optimal pure strategy Nash equilibrium of the native space function game model based on the native space utilization function comprehensive benefit matrix

S33, if

Then will be

The selection probability of the corresponding homeland space function type is determined as 1, and the selection probabilities of other homeland space function types are determined as 0;

s34, if

Calculating the optimal mixed strategy Nash equilibrium solution of the national space function game model and utilizing the optimal mixed strategy Nash equilibrium solution (x) of the decision maker for the national space _i ) ^* As the selection probability of the homeland space function type.

In an optional embodiment, the method for calculating the optimal selection probability of the homeland space function type corresponding to the homeland space function selection of the first decision maker and/or the second decision maker is as follows:

s341, constructing an expected territorial space utilization function comprehensive benefit function of a first decision maker and a second decision maker;

s342, constructing a linear programming model by taking the comprehensive benefit function of the expected homeland space utilization function as a target function, and solving the optimal solution (x) of the model _i ) ^* And the selection probability of the function type of the territorial space under the optimal mixing strategy is used.

In an alternative embodiment, the expression of the integrated benefit function of the expected territorial space utilization functions of the first decision maker and the second decision maker is as follows:

in the above formula, the first and second carbon atoms are,

respectively representing the selection of a first decision maker a and a second decision maker b

The probability of (c).

In an alternative embodiment, the linear programming model of the expected homeland space utilization function synthetic benefit function is expressed as follows:

in the above formula, v _a A value that maximizes the function synthetic benefit for the minimum expected homeland space utilization of the first decision maker.

In an optional embodiment, based on the optimal objective function value of the linear programming model, the final suitability of the homeland space function is calculated as follows:

in the above formula: a is _k ,

The k-type native space function initial suitability and the corrected native space function final suitability of the evaluation unit are respectively obtained;

and (4) an initial fitness correction coefficient for the homeland space function represents the selection probability of the decision maker l participating in the game on the homeland space function type under the optimal balance strategy of the homeland space function type k.

The second aspect of the application provides a homeland space suitability evaluation system based on an optimal balance strategy, which comprises

The resource environment bearing capacity evaluation module is used for determining a resource environment bearing capacity single evaluation element and an evaluation unit, constructing an evaluation index system, developing resource environment bearing capacity single evaluation, performing dimensionless processing on an evaluation result, determining the weight of the single evaluation element, calculating a resource environment bearing capacity comprehensive index and dividing the resource environment bearing capacity grade;

the system comprises a territorial space function initial suitability calculation module, a resource environment bearing capacity analysis module and a resource environment function initial suitability calculation module, wherein the territorial space function initial suitability calculation module is used for calculating the territorial space function initial suitability on the basis of a resource environment bearing capacity analysis result; calculating the initial suitability of the function of the homeland space comprises dividing the function types of the homeland space and determining an analysis unit, establishing a homeland space function suitability analysis index system, calculating analysis index values, carrying out dimensionless processing, determining the weight of each evaluation index, calculating a comprehensive index of the homeland space function suitability, and determining the initial suitability of the homeland space function;

and the final suitability calculation module of the homeland space function is used for constructing a homeland space function game model, obtaining the selection probabilities of various homeland space function types under the optimal equilibrium strategy, and correcting the initial suitability of the corresponding homeland space function by applying the selection probabilities of the various homeland space function types to obtain the final suitability of the homeland space function for realizing the optimal strategic equilibrium.

The system comprises a homeland space function suitability grade dividing module, a classification module and a classification module, wherein the homeland space function suitability grade dividing module is used for selecting a proper classification threshold value according to the final suitability of the homeland space function and determining the homeland space suitability grade;

and the output module is used for outputting the evaluation result of the homeland space suitability.

A third aspect of the present application provides a computer storage medium having stored thereon a computer program that, when executed by a processor, implements the steps of the homeland space suitability evaluation method.

Compared with the prior art, the application has the following advantages and beneficial effects:

according to the method, the size, the suitability and the stability of the function of the territory space are comprehensively considered, the selection probability of the function type of the territory space under the optimal balance strategy is obtained by constructing a territory space function game model and is used as a correction function of the initial suitability, the final suitability of the territory space function is constructed by multiplying the selection probability of the function quantity of the territory space, the initial suitability and the function type of the territory space, the suitability grade of the territory space is divided by the final suitability of the territory space function, and the technical defects that the space function suitability of the territory space has high-value region overlapping, the stability is poor and the comprehensive functional benefit of the territory space utilization is neglected in the conventional territory space suitability evaluation method are overcome.

Drawings

In order to more clearly illustrate the technical solutions of the exemplary embodiments of the present invention, the drawings that are required to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present invention and therefore should not be considered as limiting the scope, and that for those skilled in the art, other related drawings can be obtained from these drawings without inventive effort. In the drawings:

fig. 1 is a flowchart of a method for evaluating the suitability of a territorial space based on an optimal balancing strategy according to an embodiment of the application;

FIG. 2 is a graph illustrating the load bearing capacity of the environment of the various kinds of resources in different markets according to an embodiment of the present disclosure;

FIG. 3 is a diagram illustrating the initial suitability of the functions of the local space in the city of the Xyle's map according to an embodiment of the present application;

FIG. 4 is a diagram illustrating the final suitability of the functions of the local space in the city of Xyle's map according to an embodiment of the present application;

fig. 5 is a comparison diagram of the result of partitioning the initial suitability and the final suitability of the town production function and the town life function provided in the embodiment of the present application.

Detailed Description

In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to examples and accompanying drawings, and the exemplary embodiments and descriptions thereof are only used for explaining the present invention and are not meant to limit the present invention.

Example 1

The application provides a homeland space suitability evaluation method based on an optimal balance strategy, and as shown in fig. 1, the homeland space suitability evaluation method comprises the following steps:

s1, evaluating the bearing capacity of the resource environment, and dividing the grade of the bearing capacity of the resource environment.

And S11, determining a single evaluation element and an evaluation unit of the resource environment bearing capacity, and constructing an evaluation index system.

And S12, carrying out single evaluation on the resource environment bearing capacity based on the evaluation index system, and carrying out dimensionless processing on the evaluation result.

Because each index factor of a single evaluation element has certain difference on a metering unit and is not comparable, and a uniform grade division standard needs to be established for the selected index factors, a threshold value method is used for weighting the indexes of the single element, wherein weights such as land resource evaluation, water resource evaluation, environment evaluation, ecological evaluation and the like are assigned, disaster evaluation is assigned by an AHP method, the influence degree of each index factor change on the corresponding single element is comprehensively obtained, and a group of values reflecting the attribute characteristics of the single element are obtained by weighting, stacking and summing, wherein the calculation method of the single evaluation element for resource environment bearing capacity evaluation is shown as the following table:

TABLE 1 resource environment bearing capacity evaluation list evaluation element calculation method

And S13, determining the weight of the single evaluation requirement on the basis of the step S12.

In order to objectively and scientifically divide the grade of the bearing capacity of the resource environment and find out the comprehensive supporting capacity of a research area on life and production activities such as urban construction, agricultural development and the like, the indexes are standardized by using a threshold method in the embodiment of the application, and the relative importance of each evaluation element is determined by using an entropy method, so that the weight value of each evaluation element relative to the evaluation result of the bearing capacity of the resource environment is obtained.

And normalizing the indexes by adopting a threshold value method.

In the formula, N _i ' is the normalized score, N, of the ith element _i Is the actual score, max (N), of the ith element _i ) Min (N) which is the maximum value of the actual value of the ith element _i ) Is the minimum score of the actual scores of the ith element.

In order to avoid the influence of subjective factors on the accuracy of an evaluation result by a subjective weighting method, an entropy weighting method is applied and information entropy is used for weighting. Mainly comprises the following steps:

calculating the information entropy E of the ith element _i The formula is as follows:

wherein k is a constant number, such that

Then 0 is not less than E _i ≤1；Y _i Is the specific gravity of the element of item i.

Wherein x is _i Is the actual score of the ith element.

Calculating the difference coefficient H of the i-th element _i The formula is as follows:

H _i ＝1-E _i

to index difference coefficient H _i Carrying out normalization processing to obtain the weight M of the ith element _i The formula is as follows:

and S14, calculating the comprehensive index of the bearing capacity of the resource environment and dividing the grade of the bearing capacity of the resource environment.

The resource environment bearing capacity is comprehensively evaluated by adopting a multi-factor comprehensive evaluation method, the comprehensive influence of evaluation elements in an evaluation unit on the soil space resource environment bearing capacity is converted into an evaluation unit total score by adopting a weighting index and a model according to a single evaluation result of the resource environment bearing capacity and a corresponding weight value, the total score is used as a regional resource environment bearing capacity index value, the higher the evaluation result score is, the better the suitability of the region is represented, and otherwise, the worse the suitability is. The calculation formula is as follows:

Q＝∑M _i ·N _i ′

in the formula, Q is the comprehensive resource environment supportA load index; m _i The index weight value of the ith element; n is a radical of _i ' is the normalized score of the ith element.

And (4) obtaining the comprehensive bearing capacity index of the resource environment by utilizing a Field Calculator tool in ArcGIS and combining the standardized value and the weight of the evaluation element.

Aiming at the resource environment bearing capacity evaluation result, the natural breakpoint method is adopted for automatic grading, and the resource environment bearing capacity evaluation result of the research area is automatically generated and output, so that the resource environment bearing capacity evaluation result is divided into four grades of a high-level area, a middle-level area and a low-level area.

And S2, calculating the initial suitability of the corresponding homeland space function according to the resource environment bearing capacity evaluation result.

On the basis of the resource environment bearing capacity evaluation result, calculating the initial suitability of the homeland space function according to the resource environment bearing capacity grade, wherein the calculating step specifically comprises the following steps:

and S21, dividing the function types of the homeland space, determining an evaluation unit, and establishing a homeland space function suitability evaluation index system.

The selection of the initial index has great subjective randomness, the determination of the evaluation index is greatly influenced by personal experience, and in order to enable the evaluation index to be more scientific and reasonable, the Pearson correlation evaluation of SPSS software is utilized in the embodiment of the application, the index with high correlation is removed according to the correlation of the evaluation factor, and the selected index set is representative and not repeated.

And S22, calculating an evaluation index value in the soil space function schematic index system, and performing dimensionless processing on the evaluation index value.

In order to ensure the accuracy of the corresponding weight, the evaluation indexes are subjected to standardized processing by using a threshold value method, and the relative importance of each layer index to the previous layer target is determined by using a method combining a layer evaluation method and an entropy weight method, so that the comprehensive weight value of each layer index relative to the total target is obtained, the influence caused by artificial subjective misjudgment in subjective weighting can be reduced, and the error caused by neglecting the importance degree of the index per se in objective weighting can be avoided. The calculation steps of the hierarchical evaluation method comprise:

A. establishing a hierarchical structure model: and drawing a hierarchical structure chart containing the space function type and the evaluation index according to the space function type and the corresponding index thereof.

B. Constructing a judgment matrix: and comparing every two factors of the single-layer next time, quantifying according to 9-level scales, and constructing a judgment matrix step by step.

C. Calculating a feature vector and a weight: and calculating the maximum characteristic root and the corresponding characteristic vector of the matrix to obtain a standard paired matrix, and calculating the arithmetic mean value of each row in the matrix to obtain the weight of each level factor and factor.

D, consistency check is carried out on the judgment matrix:

in the formula, CR is a random consistency ratio, CI is a judgment matrix consistency index, and RI is an average random consistency index. The results were more satisfactory when the random consistency ratio CR was < 0.10.

And step S23, determining each evaluation index weight on the basis of the dimension removing treatment.

And S24, calculating a comprehensive index of the suitability of the homeland space function according to the weight of each evaluation index, and determining the initial suitability of the homeland space function.

Wherein, a multi-factor comprehensive evaluation method is adopted to carry out comprehensive evaluation on the suitability of the national space function. And substituting the standard value of the evaluation index of the homeland space function suitability and the corresponding weight value into the weighting index and the model for calculation, and calculating the homeland space function suitability in each evaluation unit by weighting summation. The formula is as follows:

P＝∑W _i ·V _i ′

in the above formula, P is the function suitability of a certain type of homeland space; w is a group of _i Weighted value of i-th evaluation index；V _i ' is the normalized score of the i-th evaluation index.

Aiming at the evaluation result of the national soil space function suitability, the invention adopts a natural breakpoint method to automatically grade, and automatically generates and outputs seven national soil space function suitability evaluation result tables and evaluation grade maps, such as town production space, town living space, country production space, country living space, ecological supply service space, ecological regulation service space, ecological support service space and the like. The functional suitability of the soil space of seven countries is divided into high suitability, moderate suitability and low suitability. It can be seen from the results that suitable regions of each level with the same function in the local national space in different cities are obviously partitioned, the distribution difference of suitable regions of each level between different functions in the local space is obvious, and the highly suitable regions between different functions in the similar space have obvious overlapping phenomena, which indicates that the competition phenomenon between the functions in the local space is prominent.

And S3, constructing a homeland space function game model, and solving and obtaining the selection probability of the homeland space function type under the optimal equilibrium strategy through the homeland space function game model.

The territorial space function game model is as follows:

a revenue function representing the participating decision maker;

the set of decision-making users participating in the game is l e { a, b }, where a represents a first decision-making user and b represents a second decision-making user.

The strategic space of each decider is

In the above formula, l represents a decision maker participating in the game, and l = a, b; n represents the state, the earth and the airThe total number of function types; k represents a homeland space functional type, k =1,2, \ 8230;, n;

a kth specific strategy representing the l decision maker;

a strategy set representing the ith decision maker; if the decision maker l selects a specific strategy for each person, the vector

Representing a strategic combination;

in the territorial space function game model, the territorial space utilizes a function comprehensive benefit function to express a determined utility level which can be obtained by a decision maker under a specific territorial space function combination, and the expression is as follows:

in the above-mentioned formula, the compound has the following structure,

function combination expressed in homeland space

The territorial space utilization function comprehensive benefit function of the lower decision maker l; v. of _k The method is characterized by comprising the following steps of (1) representing k-type soil space functional quantity, wherein generally, the greater the functional quantity is, the more obvious the comprehensive benefit of the soil space utilization function is; a is _k For the land space function v _k The first correction coefficient (usually, the region with larger function quantity of the territory space is also the region with high suitability of the function of the territory space, but does not exclude that the region with large function quantity is not always high in suitability under the influence of certain restrictive factors, and the phenomenon of overlapping of the territory spaces with multiple types of high function quantities may occur, so that the function quantity of the territory space needs to be corrected with the suitability assisted by the first decision maker/the second decision maker to objectively reflect that the first decision maker/the second decision maker are in a certain stateActual profit under strategic combination) representing the fitness of k-type territorial space functions, the larger the value, the higher the fitness, and a _k ∈[0,1]；p _ij For the land space function v _k The second correction coefficient (generally, it is considered that if a certain type of homeland space function in the evaluation unit is large in quantity and high in suitability, the unit is more suitable for being developed into the type of homeland space function, but the unit is influenced by external social environment factors, such a space may be converted into a type of homeland space function with other function quantity and non-maximum suitability, so that the homeland space function is corrected by the homeland space function type transfer probability on the basis of the first correction coefficient correction result, and the probability of converting the homeland space function i into the homeland space function j is represented.

Wherein, the energy v of the homeland space _k Second correction coefficient p _ij The calculation method comprises the following steps: calculating and obtaining a transfer probability matrix p for converting the territorial space function type i into the territorial space function type j in the period from t to t +1 by using a Markov model _ij As a homeland space function v _k The second correction coefficient of (1).

The specific calculation method for constructing the homeland space function game model and obtaining the selection probability of the homeland space function type under the optimal balance strategy comprises the following steps:

and S31, constructing a first decision maker and a second decision maker territorial space utilization function comprehensive benefit matrix.

In the embodiment of the application, a territorial space planner is used as a first decision maker a, and a territorial space functional system is used as a second decision maker b.

The territorial space planning is that planners combine the development and change trend of regional society, economy and ecological environment and reasonably arrange and arrange territorial space function types in a certain period in the future based on the territorial space current situation. According to the space conservation law, each time a soil space decision maker determines the space type of an evaluation unit, the corresponding soil space function system loses the space type of the evaluation unit, and therefore the soil space function game model can be simplified into a zero-sum game model. And (4) respectively setting A and B as territorial space utilization function comprehensive benefit matrixes of decision makers a and B participating in the game, and obtaining A = -B according to the zero sum game principle. Wherein, the comprehensive benefit matrix of the utilization function of the territorial space of the decision maker a is as follows:

step S32, calculating optimal pure strategic Nash equilibrium of the native space function game model based on the native space utilization function comprehensive benefit matrix

Step S33, if

Then will be

step S34, if

Calculating the optimal mixed strategy Nash equilibrium solution of the native space function game model, and using the optimal mixed strategy Nash equilibrium solution (x) of the decision maker for the native space _i ) ^* As the selection probability of the homeland space function type.

Step S341, constructing an expected territorial space utilization function comprehensive benefit function of a first decision maker and a second decision maker;

in the above formula, the first and second carbon atoms are,

respectively representing the selection of the first decision maker a and the first decision maker b

The probability of (c).

Step S342, constructing a linear programming model by taking the comprehensive benefit function of the expected homeland space utilization function as an objective function, and solving the optimal solution (x) of the model _i ) ^* And the selection probability of the function type of the territorial space under the optimal mixing strategy is used.

In the non-optimal pure strategy nash equilibrium, when the first decision-maker a selects a certain mixed strategy, it is inevitable that the second decision-maker b will select a strategy in a targeted manner to minimize the expected yield, and thus the first decision-maker a is targeted to the probability X ^* When selecting a strategy, the minimum of the expected profit (maximum profit value v) can be made as large as possible no matter what strategy the second decision maker B adopts _a ) I.e. there is E (X) for each strategy of the second decision b ^* ,Y)≥v _a ，X ^* Called the non-optimal pure strategic nash equilibrium of the first decision maker a.

v _a >0，

The linear programming problem described above can be converted into:

solving the linear programming problem to obtain an optimal solution (x' _i ) ^* On the basis, the optimum mixing strategy Nash equilibrium solution of the soil space function game model can be obtained through the following conversion (because the soil space function suitability is determined from the perspective of a soil space utilization decision maker, only the optimum selection probability calculation model of the participant a is given here), and the specific conversion formula is as follows:

(x _i ) ^* ＝(x′ _i ) ^* ·(v _a ) ^*

in the formula (x) _i ) ^* Representing the optimal selection probability of a participant a (a territorial space utilization decision maker) for the i-type territorial space function type; (x' _i ) ^* Representing an optimal solution of the linear programming model; (v) of _a ) ^* For the optimal objective function value of the unconverted linear programming model, the saddle point of the zero sum game model is expressed by the formula

And S4, correcting the initial suitability of the corresponding homeland space function by taking the selection probability of the homeland space function type under the optimal balancing strategy as a correction coefficient to obtain the final suitability of the homeland space function.

And (2) correcting the initial suitability of the homeland space function corresponding to the evaluation unit by taking the solution of the homeland space function game model, namely the optimal mixed strategy (including the optimal pure strategy) Nash equilibrium solution as a correction coefficient to obtain the final suitability of the homeland space function, wherein the calculation formula is as follows:

in the above formula：a _k ,

The initial suitability of the k-type territorial space function of the evaluation unit and the final suitability of the corrected territorial space function are respectively obtained;

and (3) expressing the kth optimal mixing strategy of the ith participant decision maker, namely the optimal selection probability of the participant decision maker to the function type k of the homeland space, for the homeland space function initial fitness correction coefficient, and obtaining the result by solving a homeland space function game model.

And selecting a proper grading threshold value according to the final suitability degree of the homeland space function, and determining the suitability grade of the homeland space.

The selection probability of the type of the homeland space function under the optimal equilibrium strategy is obtained by utilizing a homeland space function game model, the initial suitability of the homeland space function is corrected to obtain the final suitability of the homeland space function, and a proper grading threshold value is selected to determine the suitability grade of the homeland space.

Example 2

In this embodiment, a method for evaluating the suitability of a homeland space based on an optimal equilibrium strategy is further described in combination with specific conditions of various departments in the market on the basis of embodiment 1.

And S12, carrying out single evaluation on the bearing capacity of the resource environment based on the evaluation index system, and carrying out dimensionless treatment on the evaluation result.

And comprehensively calculating to obtain the weights of evaluation elements such as the land availability, the water resource abundance, the environmental pollution capacity, the ecological background characteristic value and the disaster comprehensive index, wherein the weights are shown in the following table:

table 2 evaluation elements and weights of resource environment bearing capacity

Type of evaluation	Evaluation factors	Weight of
			Evaluation of land resources	Degree of land availability	0.45
Water resource evaluation	Abundance of water resources	0.27
			Environmental evaluation	Environmental pollutant-receiving capacity	0.08
Ecological evaluation	Characteristic value of ecological background	0.19
			Disaster evaluation	Hazard danger	0.01

The resource environment bearing capacity is comprehensively evaluated by adopting a multi-factor comprehensive evaluation method, the comprehensive influence of evaluation elements in an evaluation unit on the soil space resource environment bearing capacity is converted into an evaluation unit total score by adopting a weighting index and a model according to a single evaluation result of the resource environment bearing capacity and a corresponding weight value, the total score is used as a regional resource environment bearing capacity index value, the higher the evaluation result score is, the better the suitability of the region is represented, and otherwise, the poor suitability is represented. The calculation formula is as follows:

Q＝∑M _i ·N _i ′

in the formula, Q is a resource environment comprehensive bearing capacity index; m is a group of _i The index weight value of the ith element; n is a radical of hydrogen _i ' is the normalized score of the ith element.

By utilizing a Field Calculator tool in ArcGIS and combining the standardized values and weights of the evaluation elements, the resource environment comprehensive bearing capacity index is obtained through weighted calculation, and the formula is as follows:

resource and environment comprehensive bearing capacity index = land availability × 0.45+ water resource abundance × 0.27+ environmental pollution capacity × 0.08+ ecological background characteristic value × 0.19+ disaster risk × 0.01

Aiming at the resource environment bearing capacity evaluation result, the invention adopts a natural breakpoint method to automatically grade, and automatically generates and outputs a research area resource environment bearing capacity evaluation chart and a result table. The resource environment bearing capacity evaluation result is divided into 4 grades of a high-level area, a higher-level area, a middle-level area and a low-level area. Wherein, the resource environment bearing capacity evaluation result table of qiqiqizhu city is shown as the following table:

TABLE 3 QIONGLIE City resource environment bearing capacity evaluation result table

Resource environment bearer capability level	Area (km) ² )	Specific gravity of
			High level	401.60	28.94％
Higher level	373.05	26.88％
			Middle level	359.19	25.88％
Low level	254.07	18.31％

The evaluation of the suitability of the national space function is to judge and evaluate the suitability of the natural, economic and social conditions of the national space to seven types of functional types of town production, town life, country production, country life, ecological supply service, ecological regulation service and ecological support service by comprehensively considering the guidance and connotation of the national space function on the basis of the evaluation of the bearing capacity of the resource environment. For the resource environment bearing capacity high level area, the higher level area and the middle level area, seven kinds of soil space function suitability evaluations are carried out, and for the resource environment bearing capacity low level area, only three kinds of ecological function suitability evaluations are carried out.

The land resources have multiple suitability, so that the land development can be used for multiple targets to meet different requirements of human beings, and the land development shows that the same region space has multiple functions. The qiqi (qigan) local territory space function system comprises seven functions of a town production function, a town life function, a country production function, a country life function, an ecological supply service function, an ecological regulation service function and an ecological support service function. The combined wide city town, country and ecological situation, through theoretical evaluation, standard reference, literature evidences and consulting expert's opinion, consider that the suitability of seven major functions is influenced by multiple factors such as natural ecology of land and human social economy, select the index that can represent town production function, town life function, country production function, country life function, ecological supply service function, ecological regulation service function and ecological support service function as initial index, the evaluation of homeland space suitability selects the index as shown in the following table:

TABLE 4 preliminary selection index for evaluation of homeland space function suitability

The selection of the initial index has great subjective randomness, the determination of the evaluation index is greatly influenced by personal experience, and in order to enable the evaluation index to have scientificity and rationality, the Pearson correlation evaluation of SPSS software is utilized in the embodiment of the application, the index with high correlation is removed according to the correlation of the evaluation factor, and the selected index set is representative and not repeated.

And finally establishing a homeland space function suitability evaluation index system after the evaluation factor correlation is completed. The national soil space function suitability evaluation index system comprises 17 indexes of town production function, town life function, country production function, country life function, ecological supply service function, ecological regulation service function and ecological support service function, 7 large functions and plot concentration, economic comprehensive development level, traffic grid density, population gathering level, traffic dominance degree, public service function, soil nutrient, cultivation convenience degree, plot shape index, village occupation scale, vegetation net primary productivity, water resource conservation quantity, carbon fixation and oxygen release quantity, evapotranspiration quantity, environment purification capacity, soil retention quantity, habitat quality and the like. The evaluation index system of the homeland space functional suitability is shown in the following table:

TABLE 5 evaluation index for suitability of homeland space function

Wherein, the index values of the 17-item homeland space function suitability evaluation indexes mentioned in the step S21 are respectively calculated by applying a homeland space function suitability evaluation index calculation method, and the homeland space function suitability evaluation index calculation method is shown as the following table:

TABLE 6 calculation method of evaluation index of homeland space function suitability

In order to ensure the accuracy of the corresponding weight, the evaluation indexes are subjected to standardization processing by adopting a threshold value method, the relative importance of each layer index to the previous layer target is determined by adopting a method combining the layer evaluation method and the entropy weight method, so that the comprehensive weight value of each layer index relative to the total target is obtained, the influence caused by artificial subjective misjudgment in subjective weighting can be reduced, and the error caused by neglecting the importance degree of the index per se in objective weighting can be avoided.

And step S23, determining each evaluation index weight on the basis of dimension removal processing.

According to the evaluation indexes of the dimensionless processing in the step S22, the comprehensive weight of the evaluation indexes of suitability of seven types of territorial space functions, such as town production function, town life function, country production function, country life function, ecological supply service function, ecological regulation service function, ecological support service function and the like, is determined, as shown in the following table:

TABLE 7 evaluation index weight for homeland space function suitability

Aiming at the evaluation result of the national soil space function suitability, the invention adopts a natural breaking point method to automatically grade, and automatically generates and outputs seven national soil space function suitability evaluation result tables and evaluation grade charts, such as town production space, town living space, country production space, country living space, ecological supply service space, ecological regulation service space, ecological support service space and the like, as shown in figure 3. The functional suitability of the soil space of seven countries is divided into high suitability, moderate suitability and low suitability. It can be seen from the results that suitable regions of each level with the same function in the local national space in different cities are obviously partitioned, the distribution difference of suitable regions of each level between different functions in the local space is obvious, and the highly suitable regions between different functions in the similar space have obvious overlapping phenomena, which indicates that the competition phenomenon between the functions in the local space is prominent.

As shown in figure 3, the town space comprises a town production function and a town life function, and the overlapping area of the highly suitable areas of the two functions reaches 99.37km ² . The overlapping area accounts for 84.88 percent of the area of the highly suitable area for the urban production function and 24.36 percent of the area of the highly suitable area for the urban living function. According to the evaluation result chart of the suitability of the urban production function, the low-grade suitable area of the urban production function is more, and the area is 822.57km ² Occupying 59.27% of the area of the research region; the area of the moderate adequacy zone is 448.29km ² Accounting for 32.30% of the area of the research area; the height suitable area is distributed less and has an area of 117.06km ² And occupies 8.43 percent of the area of the research region. According to the evaluation result chart of the suitability of the urban life function, the scores of the areas of the three suitable grades of the urban life function are distributed equivalently, and the areas of the low-degree suitable area, the medium-degree suitable area and the high-degree suitable area are 466.93km respectively ² 、513.08km ² 、407.91km ² Respectively accounting for 33.64%, 36.97% and 29.39% of the area of the research region.

The rural space comprises the rural production function and the rural life function, and the overlapping area of the highly suitable regions of the two functions is 90.82km ² . The overlapping area accounts for 15.74 percent of the area of the highly suitable area for the rural production function and 84.35 percent of the area of the highly suitable area for the rural living function. According to the evaluation result chart of the suitability of the production function of the village, the low-degree suitable area and the high-degree suitable area of the production function of the village are distributed widely and the areas respectively reach 640.04km ² 、577.03km ² Respectively accounting for 46.12% and 41.58% of the area of the research region; the distribution of the moderate-fitness zone is more dispersed, the area is only 170.84km2, and the area of the moderate-fitness zone accounts for 12.31 percent of the area of the research zone. According to the evaluation result chart of the suitability of the rural life function, the area of the rural life function moderate suitable area is the largest, and is 740.37km ² 53.34% of the area of the research area; the area of the low-degree suitable region is 539.87km ² Occupying the area of the study area38.90 percent; the highly suitable region is distributed most dispersedly and has an area of only 107.67km ² And occupies 7.76% of the area of the research region.

The ecological space comprises an ecological supply service function, an ecological regulation service function and an ecological support service function, and the area of the overlapped part of the three areas with high suitable functions is only 30.11km ² This is mainly because the area of the highly functional area for the ecological regulation service is small. The overlapping area of the highly suitable areas of the ecological supply service function and the ecological support service function reaches up to 441.99km ² The overlapping area accounts for 86.27% of the area with the high suitability for the ecological supply service function and accounts for 68.35% of the area with the high suitability for the ecological support service function. According to the evaluation result chart of the ecological supply service function suitability, the areas with high ecological supply service function suitability and moderate ecological supply service function suitability are widely and intensively distributed and reach 512.31km2 and 764.22km respectively ² Respectively accounting for 36.91 percent and 55.06 percent of the area of the research region; the area of the low-degree suitable area is 111.38km ² And occupies 8.03 percent of the area of the research area. According to the evaluation result chart of the suitability of the ecological regulation service function, the areas of the moderate-degree suitable area and the low-degree suitable area of the ecological regulation service function are larger and reach 701.61km respectively ² 、648.26km ² Respectively accounting for 50.55% and 46.71% of the area of the research region; the highly suitable region has a small area, the area is only 38.04km2, and the highly suitable region occupies 2.74 percent of the area of the research region. According to the evaluation result chart of the suitability of the city ecological support service function of qiqiong, the area of the highly suitable area of the qiong' an city ecological support service function is the largest, which is 646.69km ² 46.59% of the area of the research area; the area of the medium suitable area is 584.64km ² 42.12% of the area of the study area; the area of the low-degree suitable area is the minimum, and is only 156.58km ² And occupies 11.28 percent of the area of the research region.

And S3, constructing a homeland space function game model, and solving and obtaining the selection probability of the homeland space function type under the optimal balance strategy through the homeland space function game model.

And programming program codes according to the established homeland space function game model to obtain the selection probability of the homeland space function type under the optimal balance strategy.

And S4, taking the selection probability of the territorial space function type under the optimal balance strategy as a correction coefficient, and correcting the initial suitability of the corresponding territorial space function to obtain the final suitability of the territorial space function.

The method comprises the steps of obtaining selection probability of a homeland space function type under an optimal balance strategy by using a homeland space function game model, correcting initial suitability of the homeland space function to obtain final suitability of the homeland space function, selecting a proper grading threshold value, and determining the suitability level of the homeland space, wherein the selection probability of the homeland space function type under the optimal balance strategy is obtained by using the homeland space function game model, and the grade is shown in figure 4. The optimized areas with highly suitable functions in the local space of the city can be distributed in obviously non-overlapping intersecting areas, namely non-functional competition areas. The result shows that the territorial space function game model fundamentally solves the problems of space conflict and competition in the traditional multi-factor comprehensive evaluation method, greatly increases the rationality and the practicability of the territorial space suitability evaluation result as territorial space reconstruction, and can better provide service for territorial space planning.

The high-suitability areas in the evaluation of the final suitability degree high-value distinguishing and distributing local space function suitability of the qiqiong city town production function are similar, the current situation industrial areas of a central urban area and an eastern Yanzhen district are mainly concentrated, and the areas are convenient to traffic and active in economic activities and are the centers of the economic development of the whole city; the final suitability high-value areas of the urban life functions are intensively distributed in central urban areas and various rural centers, the local areas are superior in geographic position, perfect in infrastructure and concentrated in population, and the areas can provide reliable guarantee for the life quality of urban citizens. The final suitability high-value areas of the country production function are intensively distributed in flat dam areas with flat and open eastern terrains, mountain valley areas and other areas, the terrains of the areas are flatter, the soil conditions are better, the areas are closer to roads and residential points, the country production function is suitable for agricultural production, and management and product transportation are facilitated; the final high-value areas with high suitability for the living functions of the villages are concentrated near the gathering places of the residents in the rural areas, and the agricultural population and villages in the areas are distributed intensively, so that the living functions of the villages are supported powerfully. The final suitability high-value areas of the ecological supply service functions are distributed scattered, no obvious concentrated distribution area exists, and probably in the game process, the ecological supply service functions are diluted by other functions to generate a distribution crushing phenomenon because the function quantity is small and the roll-out probability is large; the final suitability high-value areas of the ecological regulation service function are mainly distributed along the west portal mountains, in the middle deep-hill areas and in the east river water systems, vegetation along the west portal mountains and in the middle deep-hill areas is luxuriant, man-made interference is relatively less, water resources in the east river water systems are rich, and the areas have strong carbon fixing and oxygen releasing capacity, evapotranspiration capacity and environmental purification capacity, so that good conditions are provided for the ecological regulation service function; the final suitability high-value areas of the ecological regulation service function are mainly distributed in the regions of south treasure mountain and nine-top mountain in the northwest, the regions have large relief, mountain mountains and mountains overlap, and the gullies move vertically and horizontally, so that the ecological regulation service function can be well undertaken, as shown in fig. 5.

The territorial space function game model based on the game theory can effectively solve the problem of territorial space function competition. The multi-factor comprehensive evaluation method is most widely applied to suitability evaluation, and although the initial suitability of the homeland space obtained by the method is consistent with the distribution result of the current state of the homeland space functions, the method has the problem of overlapping highly suitable areas, so that the problems of conflict and competition of different homeland space functions in the same homeland space are caused. Compared with the initial suitability of the homeland space function obtained by the traditional evaluation method, the final suitability of the homeland space function obtained by combining the optimal solution optimization of the game model in the embodiment of the application has the advantages that the highly suitable distinguishing distribution difference among different functions is obvious, and the conflicting competition phenomenon does not exist. Taking the highly suitable area for the town production function and the highly suitable area for the town life function, for which the initial suitability competes most intensely, as shown in fig. 5, a represents the final suitability grade condition of the town production function, B represents the initial suitability grade condition of the town production function, C represents the final suitability grade condition of the town life function, D represents the initial suitability grade condition of the town life function, a represents the final suitability grade condition of the town production function in the central urban area, B represents the initial suitability grade condition of the town production function in the central urban area, C represents the final suitability grade condition of the town life function in the central urban area, and D represents the initial suitability grade condition of the town life function in the central urban area. As can be seen by comparing the layout of the highly suitable areas of A and C, and B and D, the highly suitable area of the final suitability of the production function in the town and the highly suitable area of the final suitability of the life function in the town have no intersecting area, namely no conflict area. And the overlapping area of the highly suitable area of the initial suitability of the urban production function and the highly suitable area of the initial suitability of the urban living function is more, so that the conflict area of the highly suitable areas of the initial suitability of the two functions is large, and the competition is obvious. Therefore, the final suitability of the homeland space function obtained by combining the optimal solution optimization of the game model better meets the requirement of practical application.

Example 3

The present embodiment provides a homeland space suitability evaluation system based on an optimal equilibrium strategy based on embodiments 1 and 2, and the system includes:

the resource environment bearing capacity evaluation module is used for determining a resource environment bearing capacity single evaluation element and an evaluation unit, constructing an evaluation index system, developing resource environment bearing capacity single evaluation, carrying out dimensionless processing on an evaluation result, determining the weight of the single evaluation element, calculating a resource environment bearing capacity comprehensive index and dividing the resource environment bearing capacity grade;

1. Is novel. In the mainstream resource environment bearing capacity and homeland space development suitability evaluation guide, the method for determining the homeland space suitability is the most common in a multi-factor comprehensive evaluation method, namely, after the dimensionless processing is carried out on each single-factor index value, the weighted arithmetic mean value of each single-factor is calculated and is used as the measurement index of the homeland space suitability. The evaluation method ignores the comprehensive functional benefit of the utilization of the territorial space and the stability of the territorial space functional system, and the evaluation result often has a plurality of territorial space overlapping phenomena with high suitability (namely, a space competition area is generated), so that a planner cannot determine the space type according to the obtained suitability during the territorial space planning and cannot provide technical support for the current territorial space planning compilation. The invention innovatively introduces a complete information static game theory, a soil space utilization decision maker and a soil space function system are used as a first decision maker and a second decision maker, the soil space function type is used as a strategic space, the soil space function size is used as a gain function, a soil space function game model is established, the selection probability of the soil space function type under the optimal balanced strategy is obtained, the soil space suitability (namely the initial suitability of the soil space) obtained by a conventional method is corrected by the selection probability to obtain the final suitability, the soil space suitability is determined accordingly, the technical problem which cannot be solved by the conventional evaluation method is solved, and the method has novelty.

2. Is inventive. The key technical innovation of the invention is that on the basis of conventional homeland space suitability evaluation (S2: calculating the initial suitability of the homeland space function), a homeland space function game model is constructed by applying a complete information static game theory, the selection probability of the homeland space function type under the optimal equilibrium strategy is obtained, and the selection probability is used as a correction coefficient to correct the initial suitability of the homeland space function corresponding to an evaluation unit to obtain the final suitability of the homeland space, so that the evaluation results with high comprehensive functional benefit and strong stability can be obtained, and various conflicts of the high-suitability homeland space function types can be avoided, and the evaluation results specifically comprise:

(1) Most of the existing evaluation methods for the suitability of the territory space are multi-factor comprehensive evaluation methods, and the evaluation technical methods have three defects: firstly, the stability of the suitability is ignored, and the influence of the interaction relationship between the configuration of the planners on the homeland space and the functional system of the homeland space on the stability of the suitability is not considered; secondly, the comprehensive benefit of the utilization function of the soil space is not considered in the obtained suitability, namely, the problem of low comprehensive functional benefit of the soil space pattern constructed according to the highly suitable region can exist; thirdly, the evaluation result may cause the problem of overlapping of a plurality of types of spaces highly suitable for the territory, and effective reference cannot be provided for determining the function type of the territory space.

(2) A homeland space utilization decision maker and a homeland space function system are used as a first decision maker and a second decision maker to construct a zero-sum game model, and an optimal mixed strategy (including an optimal pure strategy) Nash equilibrium solution is obtained, namely, the optimal selection probability of each homeland space function type is obtained when the benefits of both the homeland space utilization decision maker and the homeland space function system are maximized. According to the game theory, the selection probability is the best selection probability which is unwilling to break by both parties, has stronger stability, and can be used as a correction parameter of the initial suitability of the homeland space function, thereby solving the technical defect that the suitability stability is not considered in the current homeland space suitability evaluation method.

(3) The method comprehensively considers the size, the suitability and the stability of the function of the homeland space, takes the product of the functional quantity of the homeland space, the initial suitability and the transfer probability as a gain function, constructs a zero sum game model, obtains a Nash equilibrium solution (namely the optimal selection probability of the first decision maker a on the function type of the homeland space) which enables the gains of both game parties to reach the maximum, and takes the Nash equilibrium solution as a correction coefficient of the initial suitability of the function of the homeland space, thereby solving the technical defect that the comprehensive functional benefit of the utilization of the homeland space is neglected in the existing evaluation method of the suitability of the homeland space.

(4) Starting from the perspective of a territorial space utilization decision maker, converting a zero sum game model into a linear programming model, solving to obtain the optimal selection probability of the territorial space utilization decision maker on the function types of various territorial spaces, correcting the initial suitability of the corresponding territorial space function by using the selection probability to obtain the final suitability, determining the territorial space suitability grade, and solving the technical defect of high-suitability territorial space competition of evaluation results obtained by the existing territorial space suitability evaluation method.

3. And (4) applicability. The method solves the technical defects that the comprehensive function benefit of utilization of the national soil space and the stability of a functional system are neglected and a plurality of high-suitability national soil space conflicts exist in the evaluation result in the existing technical method for evaluating the suitability of the national soil space, the maximization of the comprehensive function benefit of utilization of the national soil space is realized by the obtained evaluation result, the stability of the suitability of the national soil space is ensured, and the conflict problem of the plurality of high-suitability national soil spaces is coordinated, so that the omnibearing technical support can be provided for various specific works of planning of the national soil spaces, such as optimization and reconstruction of the national soil space pattern, functional partitioning of the national soil space, marking of three regions (three types of space, agricultural space and ecological space, marking of a town development boundary, a permanent basic farmland protection red line and three control lines of the ecological protection red line). For example, highly suitable regions for ecological spaces such as supply services, regulation services, support services and the like can be used as the basis for defining an ecological red line, and highly suitable regions for rural production and rural living spaces can be used as the basis for optimizing the layout of agricultural spaces.

The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are merely exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims

s4, taking the selection probability of the type of the homeland space function under the optimal balancing strategy as a correction coefficient, and correcting the initial suitability of the corresponding homeland space function to obtain the final suitability of the homeland space function;

2. The territorial space suitability evaluation method based on the optimal balancing strategy of claim 1 is characterized in that the territorial space function game model is as follows:

a revenue function representing the participating decision maker; wherein the content of the first and second substances,

the set of decision-making persons participating in the game is l E { a, b }, wherein in the formula, a represents a first decision-making person, and b represents a second decision-making person;

the strategic space of each decider is

In the above formula, l represents the decision maker participating in the game, l = a, b; n represents the total number of the function types of the territorial space; k represents a homeland space functional type, k =1,2, \ 8230;, n;

a kth specific strategy representing the l decision maker;

Representing a strategic combination;

the revenue function for each decision maker is

In the above formula, v _k Represents k-type native space function quantity, a _k Represents the initial suitability of k-type soil space function, a _k ∈[0,1]，p _ij Representing the probability of a transition from the homeland space i to the homeland space j.

3. The method as claimed in claim 2, wherein the variable v in the revenue function of each decision maker is a variable _k 、a _k 、p _ij The determination method of (2) is as follows:

k-type homeland space functional quantity v _k By passing

Determining; wherein (r) _k ) _m The functional quantity of the Nth sub-function of the k-type homeland space function; (omega) _k ) _m The N-th sub-function weight value is the k-type homeland space function;

initial suitability a of k-type homeland space function _k Obtaining the initial suitability determination of the homeland space function through calculation;

probability p for converting homeland space function type i into homeland space function type j _ij Calculating and obtaining a transfer probability matrix p for converting the territorial space function type i into the territorial space function type j in the period from t to t +1 through a Markov model _ij And (4) determining.

4. The method for evaluating the suitability of the territorial space based on the optimal balancing strategy according to claim 2, wherein the method for obtaining the selection probability of the territorial space function type under the optimal balancing strategy is as follows:

S33, if

Then will be

s34, if

Calculating the optimal mixing strategy Nash equilibrium solution of the native space function game model, and utilizing the optimal mixing of the decision maker in the native spaceSynthetic strategy Nash equilibrium solution (x) _i ) Probability of selection as homeland space function type.

5. The method for evaluating the territorial space suitability based on the optimal equilibrium strategy according to claim 4, wherein the method for calculating the optimal selection probability of the territorial space function type corresponding to the territorial space function selection of the first decision maker and/or the second decision maker is as follows:

s342, constructing a linear programming model by taking the comprehensive benefit function of the expected homeland space utilization function as an objective function, and solving the optimal solution (x) of the model _i ) ^* And the selection probability of the function type of the territorial space under the optimal mixing strategy is used.

6. The territorial space suitability evaluation method based on the optimal equilibrium strategy of claim 5, wherein the expression of the expected territorial space utilization function comprehensive benefit function of the first decision maker and the second decision maker is as follows:

in the above formula, the first and second carbon atoms are,

The probability of (c).

7. The method for evaluating the suitability of the territorial space based on the optimal balancing strategy according to claim 6, wherein the linear programming model of the comprehensive benefit function of the expected territorial space utilization is expressed as follows:

in the above formula, v _a Showing when the second decision maker b takes strategy

And the first decision maker a obtains the maximum value of the minimum expected territorial space utilization function comprehensive benefit function.

8. The method for evaluating the suitability of the territorial space based on the optimal balancing strategy according to claim 1, wherein the method for calculating the final suitability of the territorial space function is as follows:

in the above formula: a is _k ,

9. A homeland space suitability evaluation system based on an optimal balance strategy is characterized by comprising:

and the final suitability calculation module of the homeland space function is used for constructing a homeland space function game model, obtaining the selection probabilities of various homeland space function types under the optimal balance strategy, and correcting the initial suitability of the corresponding homeland space function by applying the selection probabilities of various homeland space function types to obtain the final suitability of the homeland space function for realizing the optimal strategic balance.

The system comprises a homeland space function suitability grade dividing module, a classification module and a judging module, wherein the homeland space function suitability grade dividing module is used for selecting a proper grading threshold value according to the final suitability of the homeland space function and determining the homeland space suitability grade;

10. A computer storage medium, characterized in that the computer storage medium has stored thereon a computer program which, when executed by a processor, implements the steps of the homeland space suitability assessment method according to any one of claims 1 to 8.