Regional Security Assessment of Integrated Energy Systems with Renewables in China: A Grid-Connected Perspective

: Stable and reliable integrated energy systems are one of the major issues related to sustainable regional and national energy development. Because most existing studies are conducted on whole countries, few address the e ﬀ ects of regional interaction and renewable energy. Therefore, a natural disaster risk assessment model (NDRAM) combined with spatial models is used as a general systematic tool to assess and resolve regional energy security, based on a framework of resources, generation, transmission, marketing and consumption, with 17 metrics. In particular, energy systems were treated as organic connected-units and their security status was regarded as a combined result of potential hazard and system vulnerability. The proposed method was applied to evaluate and classify the security situation of 31 Chinese provinces in 2016. The results showed that transmission had the most signiﬁcant impact among ﬁve major risk sources. The closer grid connections have a stronger ability to deal with risks among regions, where renewables consumption could be better stimulated cross-regionally. In terms of a regional perspective, there is still a gap among di ﬀ erent regions, and eastern China presented higher energy risk status. The most energy-hazard provinces are mainly in the east provinces with well-developed levels in Beijing, Tianjin and Shanghai. The least energy-vulnerable provinces are mainly in the abundant natural resources regions such as Inner Mongolia, Sichuan and Xinjiang. The NDRAM-based general model provides a systematic tool for quantitative assessment of regional energy security with a full accounting of regional interaction and renewable energy issues, which may help to develop clean energy, optimize system infrastructure and improve scientiﬁc management.


Introduction
Energy security has always been a high priority in stable and sustainable energy management for national and regional development. Especially in the era of global climate change and environmental concerns, security issues with renewable energy have been increasingly outlined from both a scientific and political perspective. In the last decade, renewable energy has rapidly increased and provided an effective clean guarantee for energy security. A total of 181 GW of renewable power was added, a consistent pace compared to 2017, and the number of countries integrating high shares of variable renewable energy keeps rising. (REN21, 2019). On the other hand, their rapid expansion poses new challenges for energy systems, such as cross-region scheduling, low reliability, occurring overproduction, etc., which remind us that sufficient security attention needs to be paid to energy systems with renewable energy. To further confirm the energy security situation of different countries, regional integrated power systems have been seen as vital units to provide operational flexibility on national security management. In particular, because the availability of renewable energy generally reflects large differences in resources and demand based on geography, the role of cross-regional grids in energy security should not be ignored. Therefore, systematic studies of the security situations of integrated energy systems with large renewable penetration and grid connections should be conducted, in order to form a greater synergic advantage in the process of energy transition and clean development.
Despite its importance, limited attention has been paid to the perspective of different regions within the country based on cross-regional grids. Fortunately, the studies conducted on energy security generally include definitions and assessments on global, national and regional scales which can provide helpful references.
The definitions of energy system security are varied as provided by different researchers in academic papers, governmental and industrial publications [1][2][3]. According to the objectives and emphasis of specific studies, energy systems can be illustrated as objects exposed to security threats, or referred to as security of demand and supply or a subject generating or enhancing insecurity caused by economic, political, technology and environmental risk factors [4,5]. The indicators generally concern two inter-related dimensions: a physical dimension, referring to available, reliable and/or assessable energy supply [6,7], and an economic dimension that incorporates aspects such as price volatility and affordability [8,9]. Concerning the security studies focused on renewable energy participation, the availability and affordability of energy security are still most impactful, and the promotion of renewable energy is compelling for national energy security strategies, which could bring environmental benefits and sustainable security [10]. Although the specific descriptions vary, the main consensus is that any jurisdiction that has an integrated indispensable energy system is responsible for meeting its end-use energy demands, safely and reliably, regardless of its level of development or size [11].
In terms of security assessment and evaluation, they have been quantitatively measured by a number of scientific assessment methods from various angles based on an adequate energy security framework [12]. Indicator approaches commonly focus on empirical calculations of indicators [13] or simulation tools to analyze interdependence among various energy sources [14]. Moreover, the indexes widely used in the disciplines of ecology, cybernetics and genetics have been cross-applied in energy security, such as the Shannon-Wiener index in ecology [15] and the Herfindahl-Hirschman index in competition economics [16]. Although various indicators and methods have been developed, certain issues still need to be addressed. Specifically, since energy systems with the characteristics of geographical dependency change dynamically, there are some studies conducted from spatial and/or temporal perspectives [17], focusing on specific areas at the national [18], regional [19], or interstate [20] levels. With the prevailing scientific attitude on the global economic and geopolitical scene, a spatial geo-economic perspective has been proposed gradually in energy security research [21]. The effectiveness of existing research is somewhat limited since interactions among grid-connected regional energy systems were not systematically considered, which may consequently be considered an important factor in energy security evaluation [3].
China is facing pressure on energy security [22] and is home to nearly 30% of the world's renewable power capacity, totaling approximately 647 GW and achieving 26.4% of the renewable share of total energy generation of China in 2017 [23]. Meanwhile, China has a great spatial disparity of energy endowment, economic development and industrial structure, and even of social and cultural customs at the provincial level [24,25]. Therefore, limited research on regional energy security issues in China has been carried out in recent years. The results show that all provinces confronted threats related to energy availability and diversity, and also identified the most energy-secure regions. From a grid-connected perspective, China has a complete and efficient regional transmission and distribution network and plans to build 16 new ultra-high voltage direct current lines, in order to connect wind and solar power regions to demand centers. Therefore, China is a most suitable sample for a case study of regional energy system security with rapid renewable energy expansion and developed power grid dispatch.
Moreover, grid interconnection among provinces from a new spatial perspective for determining the power security performance of Chinese provincial energy systems would be novel and meaningful [26].
In summary, existing studies on energy security have the limitation of paying little grid-systematic attention to regions' interaction within the country accompanied by a high proportion of renewable penetration. Therefore, this paper aims to develop a comprehensive energy security assessment criteria framework for the 31 Chinese provinces that consists of 17 indicators from 5 dimensions, thereby highlighting spatial disparities and interaction issues and providing a new insight to consider the energy systems as interconnected and recoverable systems in security issues. In an effort to achieve this goal, a natural disaster risk assessment model is used in this paper as a fundamental system-oriented framework, in order to establish a comprehensive framework matched with the definition of resilient energy systems. Moreover, a systematical assessment framework with spatial grid transmission dimension and high renewable energy penetration factor is likely to be more reliable, and would enrich the research interests and understanding of regional security issues in China.
The remainder of the paper is organized as follows: The next section introduces the assessment framework and methodology that have been used to measure energy security. This includes the methods selected and how the main methods are aggregated to form system comprehensive risk index (SCRI) sub-indexes. Section 3 presents the security assessment framework and processing procedure of the integrated energy system. Section 4 discusses the main results concerning the Chinese empirical research of the proposed framework. Finally, conclusions and suggested policy and research are presented in Section 5.

Methodology
The security assessment framework and methodology are presented briefly in this section, as shown in Figure 1. This study constructed the assessment framework and empirical study in four steps. First, a 17-metric framework was established based on the literature review from 5 risk sources: resources, generation, marketing, consumption and transmission. The transmission indexes were modeled and described by a spatial weight matrix and Moran's I test models, which are expounded in Section 2.1. Next, the paper is designed to provide a system comprehensive risk index (SCRI) and two sub-indexes, a hazard index (HI) and vulnerability index (VI), based on a natural disaster risk assessment method for traceable analysis, which would better break down regional risk situations as illustrated in Section 2.2. Then, the selected indicators are normalized. The normalized indicators are objectively weighted through an entropy method. Finally, empirical studies and classification results for the 31 Chinese provinces were carried out.
Sustainability 2020, 12, x 3 of 17 from a new spatial perspective for determining the power security performance of Chinese provincial energy systems would be novel and meaningful [26].
In summary, existing studies on energy security have the limitation of paying little gridsystematic attention to regions' interaction within the country accompanied by a high proportion of renewable penetration. Therefore, this paper aims to develop a comprehensive energy security assessment criteria framework for the 31 Chinese provinces that consists of 17 indicators from 5 dimensions, thereby highlighting spatial disparities and interaction issues and providing a new insight to consider the energy systems as interconnected and recoverable systems in security issues.
In an effort to achieve this goal, a natural disaster risk assessment model is used in this paper as a fundamental system-oriented framework, in order to establish a comprehensive framework matched with the definition of resilient energy systems. Moreover, a systematical assessment framework with spatial grid transmission dimension and high renewable energy penetration factor is likely to be more reliable, and would enrich the research interests and understanding of regional security issues in China.
The remainder of the paper is organized as follows: The next section introduces the assessment framework and methodology that have been used to measure energy security. This includes the methods selected and how the main methods are aggregated to form system comprehensive risk index (SCRI) sub-indexes. Section 3 presents the security assessment framework and processing procedure of the integrated energy system. Section 4 discusses the main results concerning the Chinese empirical research of the proposed framework. Finally, conclusions and suggested policy and research are presented in Section 5.

Methodology
The security assessment framework and methodology are presented briefly in this section, as shown in Figure 1. This study constructed the assessment framework and empirical study in four steps. First, a 17-metric framework was established based on the literature review from 5 risk sources: resources, generation, marketing, consumption and transmission. The transmission indexes were modeled and described by a spatial weight matrix and Moran's I test models, which are expounded in Section 2.1. Next, the paper is designed to provide a system comprehensive risk index (SCRI) and two sub-indexes, a hazard index (HI) and vulnerability index (VI), based on a natural disaster risk assessment method for traceable analysis, which would better break down regional risk situations as illustrated in Section 2.2. Then, the selected indicators are normalized. The normalized indicators are objectively weighted through an entropy method. Finally, empirical studies and classification results for the 31 Chinese provinces were carried out.

Moran's I Test Models
Moran's I test models are commonly used to study spatial autocorrelation, and are used to describe and test the grid connection among regions in this paper. Global Moran's I is the earliest method applied to global clustering tests [27], which can judge whether the distribution features of the adjacent regions' attributes are random, clustered or dispersed. Global Moran's I is expressed as where, n is total number of study areas; x i and x j represent the variable at location i and j; x = 1 n n i=1 x i is the mean value of the variable with the sample number; 2 is the variance of the variable; and w ij is the spatial weight, which can be represented based on inverse distance, fixed distance, K nearest neighbors, contiguity edges only/corners and customize. Local Moran's I index is used to test whether there are similar or different observed values in local areas. The local Moran's I index of the region i measures the relevance degree of a region with its connected regions. It is defined as A positive local Moran's I value refers to a study region that has high or low values similar to its neighbors, which means spatial cluster. On the other hand, a negative local Moran's I value indicates a potential spatial outlier that is different from the values of its surrounding locations [28].
A custom grid interconnected spatial weight matrix was established based on the situation of Chinese grid subdivision and UNV transmission lines, in order to accurately reflect the system connectivity in various provinces. The result of the local Moran's I index was selected as an important index to measure the tightness of transmission of each province after it was pre-checked by the global Moran's I index.

Natural Disaster Risk Assessment Method
Risk assessment of natural disasters is defined as an assessment of both the occurrence probability and damage degree of natural disasters. With the continuous expansion of research, this method has been widely used in global climate change, human health, energy systems and other research fields [28][29][30][31][32]. Combining different research objects and perspectives, the concept of risk gradually evolved into multiple conceptual sets of concepts such as inclusion, sensitivity, adaptability and resilience. The common concept was proposed by the United Nations Development Programme [33] as Function (4). Hazard is the dangerousness of disaster-causing factors. Vulnerability can be understood as the vulnerability of the hazard-affected body or the vulnerability of the disaster system according to the user's requirements.
Considering the applicability of the natural disaster method to an energy system security field, the hazard of energy systems refers to the possibility that the system cannot guarantee operation with the existing conditions. The hazard is various among regions with different basic conditions. As for the same degree of hazard in different regions, the degree of damage caused by risk is also different; namely, the vulnerability varies. Vulnerability is used to represent a region's susceptibility to risks and the ability to respond. A greater hazard index and lower vulnerability index would reflect a poorer energy security situation. System comprehensive risk index (SCRI) is expressed as where R s is the system comprehensive risk index (SCRI); H is the hazard index (HI); V is the vulnerability index (VI); R s , H , V ∈ [0, 1]; x h ,y n are the standardized values of the system hazard and vulnerability index; w h , w v are the weight of the two indexes. The metrics from the 5 risk dimensions would be divided into the hazard index (HI) and vulnerability index (VI) to carry out the assessment and empirical analyses of the integrated energy systems.

Integrated Energy System Security Assessment Framework
To provide a comprehensive and systematic security assessment of an integrated energy system at a provincial level, this study synthesized various dimensions and metrics from the literature and official publications. The paper presents the risk sources of interprovincial integrated energy systems as resources, generation, transmission, marketing and consumption, which take the classification example of the study of traditional energy system security risk appraisal [34]. The framework is a synthesis of the existing indexes, taking full account of the electricity grid and renewable energy. As Table 1 presents, these 5 risk sources can be effectively broken down into 17 components, and thereby measured with 17 metrics. The selection of components for each dimension is based on an adequate literature review, with detailed discussion of each dimension in Sections 3.1.1-3.1.5.

Resources
Energy availability is the top consideration in energy security definitions over time. The self-sufficiency and developed degree of exploitability can reflect the capacity of available resources in a certain area. Since there is a growing need of awareness to utilize energy resources in an environmentally-friendly way, environmental and sustainability indicators have increasingly become part of energy security considerations. A diversity factor is the key issue to determine the availability of energy resources [47]. The source of diversity is significantly dependent on the degree of the technical exploitation of the renewable energy.

Generation
Generation risks can be caused by the capacity or the pollution in a power generation process. The generation capacity can reflect the annual output and the utilization characteristics of generation equipment. Renewable sources can significantly enhance the sustainability situation, and their generation amount represents the cleanliness level of regional generation. In particular, utilization and cleanliness metrics should be of reasonable scope; though there has much room to improve in Chinese actual operating situations, the metrics were set as greater preferred indexes.

Transmission
The infrastructure of transmission is integral in providing stable and uninterrupted energy supply. Transmission risk is more common because the grid covers a wide and complex area with multiple nodes and electricity losses. To reflect the trans-regional, multi-channel and large-scale transmission and deployment of power, validated spatial correlation analysis [40] was applied in the transmission dimension. In depicting the transmission risks, this paper considered adequate investments in infrastructure to be proof and guarantee of energy security, while the connected situation and the convergence degree among adjacent areas would reflect resistance [48].

Marketing
Economic development is dependent on the energy system; in turn, improvement of the energy system is dependent on economic and electricity marketing factors. Basically, energy price determines the affordability of energy supplies [1]. The more electricity required per GDP means a stronger capacity to absorb, and an integrated system has less risk when there is strong energy demand. The elasticity coefficient of energy production and consumption expresses the advanced or backward conditions between the power industry and regional economic development. Generally, the development of electric power should be faster than that of the regional economy. If this is backward, the risk to the energy system will be greatly increased.

Consumption
Consumption risk can result from a rapid expansion of the amount of consumption and the pollution created in this process. Consumption expansion can be measured by annual electricity consumption and the difference between supply and demand, while the consumption situation for renewable energy can reflect the potential environmental pollution caused by consumption.

Data Collection and Processing
Considering the incomplete statistics of Hong Kong, Macao and Taiwan, the basic data of the 31 Chinese provinces in 2016 were collected from the 2016 China Statistics Yearbook [49]. Table 2 presents the characteristics of the provincial energy security original data, the standardized mean and the weight of each metric. The weights of the indicators can be assigned based on expert opinions or other subjective procedures. The study selected the entropy evaluation method, a more objective approach, to give weight to the metrics, and aggregated the weighted indicators into a composite index. Since the selected metrics are in different units of measurement, three steps are needed to aggregate them to form a composite index. They are (a) normalizing the indicators, (b) weighting the normalized indicators, and (c) aggregating the normalized indicators. The common practice in normalization is using one of the following methods: min-max, distance to reference and standardization. The metrics data are normalized on a min-max approach by the following equations: Greater preferred metrics Smaller preferred metrics where x ib is the normalized value of the province i with respect to the metric b, x ib represents the original value and maxx b , minx b are the maximum and minimum value of the 31 provinces with respect to the metric b.
Based on the normalized original data, the weight of the metrics was estimated with Equations (9) and (10) where w b is the weight of metric b, w b ∈ [0, 1]; H b is the information entropy; and f b is a metric's weight.
The comprehensive index of a given province is calculated in Equation (6).

Convergence Characteristics of Interprovincial Grid Connection
In order to analyze the components of connectivity (A31) and convergence (A32), the situation of the grid network in the 31 Chinese provinces is collected and scientifically simplified into regional power grids and cross-regional ultra-high voltage (UHV) transmission lines.
As for the regional power grids, they are composed of adjacent provincial power grids with a dispatch and management center. As a consequence, the provinces in the interprovincial grid are strongly connected and can act as reliable security providers reciprocally. Table 3 lists six Chinese regional power grids and their comprised provinces. Table 3. Chinese six regional power grids and their comprised provinces.

Regional Power Grid Provinces Number
Northeast In terms of cross-regional UHV lines, China had 11 operated UHV lines until 2016 as shown in Table 4. These UHV lines delivered 172.4 billion kilowatt-hours of renewable energy accounting for 74% of total transmission capacity. Among them, lines 1-4 are AC-UHV lines with substations. They can make up the national UHV backbone network according to the actual needs of power distribution, load distribution, transmission power and power exchange. The lines 5-11 are DC-UHV transmissions without intermediate point, namely, the power can be transferred directly to the load center by point-to-point, high power and long distance. Based on the situations of Chinese grid subdivisions and UHV transmission lines as drawn in Figure 2, a custom grid interconnected spatial weight matrix was established. In the matrix, the connected relations of every province assume that one province is closely associated with other provinces in the regional grid and through UHV transmission lines. Using the weight matrix to carry out the Moran tests, the pretest results of global Moran's I model for various attributes are shown in Table 5. provinces in the regional grid and through UHV transmission lines. Using the weight matrix to carry out the Moran tests, the pretest results of global Moran's I model for various attributes are shown in Table 5.  Significant weight was given to the P value obtained by the statistical significance test method; p < 0.05 means statistic difference, and p < 0.01 means statistically significant. According to the spatial autocorrelation report, renewable electricity consumption amount and non-hydro renewable electricity consumption amount are clustered and distributed at a significant level of 0.05; the ratio of non-hydro renewable energy consumption is clustered at a significant level of 0.01, as shown in Figure 3, with the most significant level.  Significant weight was given to the P value obtained by the statistical significance test method; p < 0.05 means statistic difference, and p < 0.01 means statistically significant. According to the spatial autocorrelation report, renewable electricity consumption amount and non-hydro renewable electricity consumption amount are clustered and distributed at a significant level of 0.05; the ratio of non-hydro renewable energy consumption is clustered at a significant level of 0.01, as shown in Figure 3, with the most significant level. The test result is consistent with the current situation of renewable energy, especially illustrating that wind and photovoltaic power can be exploited and utilized by the electricity grid. Further, the links among provinces through transmission infrastructure are better reflected by a non-hydro renewable electricity consumption ratio. Therefore, the local Moran's I model is performed as a resulting code of the metric B32-Local spatial autocorrelation, as shown in Figure 4. According to the results of the local Moran's I index of the 31 provinces, the northeastern region is a typical high-value cluster, which is closely related through grid connection and generally at a high ratio of non-hydro renewable energy consumption. Central China and East China are typical low-value clusters. The intra-and inter-regional grid connections are closely related, and there is generally a low ratio of non-hydro renewable energy consumption. The level of the Beijing-Tianjin-Hebei in North China is quite different from that in Shandong and Shaanxi. In the southern regions with low general prevalence, Hainan and Yunnan are regional outliers with high value based on their consumption. The northwest region is at a higher level except in Shaanxi. In conclusion, the UHV grid has a significant effect on regional linkages, and the regional outliers have weaker connections  The test result is consistent with the current situation of renewable energy, especially illustrating that wind and photovoltaic power can be exploited and utilized by the electricity grid. Further, the links among provinces through transmission infrastructure are better reflected by a non-hydro renewable electricity consumption ratio. Therefore, the local Moran's I model is performed as a resulting code of the metric B32-Local spatial autocorrelation, as shown in Figure 4. The test result is consistent with the current situation of renewable energy, especially illustrating that wind and photovoltaic power can be exploited and utilized by the electricity grid. Further, the links among provinces through transmission infrastructure are better reflected by a non-hydro renewable electricity consumption ratio. Therefore, the local Moran's I model is performed as a resulting code of the metric B32-Local spatial autocorrelation, as shown in Figure 4. According to the results of the local Moran's I index of the 31 provinces, the northeastern region is a typical high-value cluster, which is closely related through grid connection and generally at a high ratio of non-hydro renewable energy consumption. Central China and East China are typical low-value clusters. The intra-and inter-regional grid connections are closely related, and there is generally a low ratio of non-hydro renewable energy consumption. The level of the Beijing-Tianjin-Hebei in North China is quite different from that in Shandong and Shaanxi. In the southern regions with low general prevalence, Hainan and Yunnan are regional outliers with high value based on their consumption. The northwest region is at a higher level except in Shaanxi. In conclusion, the UHV grid has a significant effect on regional linkages, and the regional outliers have weaker connections and capabilities in responding to regional energy system risks.  According to the results of the local Moran's I index of the 31 provinces, the northeastern region is a typical high-value cluster, which is closely related through grid connection and generally at a high ratio of non-hydro renewable energy consumption. Central China and East China are typical low-value clusters. The intra-and inter-regional grid connections are closely related, and there is generally a low ratio of non-hydro renewable energy consumption. The level of the Beijing-Tianjin-Hebei in North China is quite different from that in Shandong and Shaanxi. In the southern regions with low general prevalence, Hainan and Yunnan are regional outliers with high value based on their consumption. The northwest region is at a higher level except in Shaanxi. In conclusion, the UHV grid has a significant effect on regional linkages, and the regional outliers have weaker connections and capabilities in responding to regional energy system risks.

Weight of Grid-Connected Integrated Energy System Security Metrics
Considering the objectivity and reliability of the metrics, the weight of the 17 metrics was determined by the standardized processing and entropy method. The weights of the 17 metrics are shown in Figure 5, and the weights of the 5 dimensions (resources, generation, transmission, marketing and consumption) were presented as follows:

Weight of Grid-Connected Integrated Energy System Security Metrics
Considering the objectivity and reliability of the metrics, the weight of the 17 metrics was determined by the standardized processing and entropy method. The weights of the 17 metrics are shown in Figure 5, and the weights of the 5 dimensions (resources, generation, transmission, marketing and consumption) were presented as follows: According to the data and weights analysis results, the metrics were sorted at small intervals. Among nine hazard indexes, B42-Electricity consumption per unit GDP (0.157) and B52-Consumption cleanliness (0.153) had a greater influence on the risks of the system, with electric power consumption and production given as the risk source, while the impact result from B11-Selfsufficiency (0.082) and B21-Utilization hours of installed capacity (0.043) were on a small scale. This can be interpreted as showing that the supply and development of energy resources are relatively safe in China. Among eight vulnerability indexes, B32-Local spatial autocorrelation (0.190) and B34-Electricity transmission losses (0.174) have the greatest impact on system vulnerability, showing that the enhancement of power grid connection and transmission efficiency can improve the resistance ability of energy systems to face random risk. B43-Generation maturity is the least influential factor, which indicates that the current generation level can better withstand risk.
From the perspective of risk source analysis, transmission (0.615) had the most significant impact. After this, the marketing dimension (0.416) and consumption dimension (0.404) also had larger impacts on grid-connected integrated energy security performance. It can be seen that marketing has a great influence on the risk of the integrated energy system in terms of price affordability, consumption capacity, consumption structure, demand difference and economic and social development. In the future, rational consumption structure adjustment and electric energy substitution can be used to reduce the risk to the system.

Energy Security Performance of Integrated Energy System
According to the results of the data and weight processing, the hazard index (HI), vulnerability index (VI) and system comprehensive risk index (SCRI) of the 31 Chinese provinces were evaluated. According to the data and weights analysis results, the metrics were sorted at small intervals. Among nine hazard indexes, B42-Electricity consumption per unit GDP (0.157) and B52-Consumption cleanliness (0.153) had a greater influence on the risks of the system, with electric power consumption and production given as the risk source, while the impact result from B11-Self-sufficiency (0.082) and B21-Utilization hours of installed capacity (0.043) were on a small scale. This can be interpreted as showing that the supply and development of energy resources are relatively safe in China. Among eight vulnerability indexes, B32-Local spatial autocorrelation (0.190) and B34-Electricity transmission losses (0.174) have the greatest impact on system vulnerability, showing that the enhancement of power grid connection and transmission efficiency can improve the resistance ability of energy systems to face random risk. B43-Generation maturity is the least influential factor, which indicates that the current generation level can better withstand risk.
From the perspective of risk source analysis, transmission (0.615) had the most significant impact. After this, the marketing dimension (0.416) and consumption dimension (0.404) also had larger impacts on grid-connected integrated energy security performance. It can be seen that marketing has a great influence on the risk of the integrated energy system in terms of price affordability, consumption capacity, consumption structure, demand difference and economic and social development. In the future, rational consumption structure adjustment and electric energy substitution can be used to reduce the risk to the system.

Energy Security Performance of Integrated Energy System
According to the results of the data and weight processing, the hazard index (HI), vulnerability index (VI) and system comprehensive risk index (SCRI) of the 31 Chinese provinces were evaluated. The 31 provinces are divided into potential risk areas, low security areas, moderate security areas and high security areas by using the natural breakpoint method for the SCRI. The HI/VI index score of the 31 provinces in each subdivision are presented in Figure 6. The 31 provinces are divided into potential risk areas, low security areas, moderate security areas and high security areas by using the natural breakpoint method for the SCRI. The HI/VI index score of the 31 provinces in each subdivision are presented in Figure 6. From the origin of the risks, the results show that the security areas of the 31 Chinese provinces are consistent with the regional economic and energy construction development levels. For example, Beijing is classified into potential risk areas with a lower energy self-sufficiency rate and proportion of renewable energy generation and consumption, and the development of its electricity market consumption lags behind the level of its economic and social development, resulting in higher system risk. Meanwhile, its potential development of renewable energy is lower, and the amounts of installed capacity and energy generation are small. Namely, the ability to resist risks is low and the integrated energy system is relatively fragile. Under the effects of systemic risk and resistance, the integrated energy system is relatively fragile. In the high security area, Inner Mongolia is fully selfsufficient and the technical amount of renewable energy ranks first among provinces. Further, its production of renewable energy is well-consumed, so the system risk is low. According to the hazard and vulnerability situation of the integrated energy system, the risk situation and causes of each province could be traced.
From the perspective of geographical distribution, the study combined ArcGIS and the provincial energy system risk results to obtain a corresponding spatial distribution profile as shown in Figure 7. In terms of hazard, the eastern region is more dangerous than the central and western regions, and the overall decline is from the east to the west. In terms of vulnerability, the central region is more vulnerable because of the close grid-connection. Figure 7c visually presents the results of the risk assessment of the 31 Chinese provinces under the combined effects of hazard and vulnerability. From the origin of the risks, the results show that the security areas of the 31 Chinese provinces are consistent with the regional economic and energy construction development levels. For example, Beijing is classified into potential risk areas with a lower energy self-sufficiency rate and proportion of renewable energy generation and consumption, and the development of its electricity market consumption lags behind the level of its economic and social development, resulting in higher system risk. Meanwhile, its potential development of renewable energy is lower, and the amounts of installed capacity and energy generation are small. Namely, the ability to resist risks is low and the integrated energy system is relatively fragile. Under the effects of systemic risk and resistance, the integrated energy system is relatively fragile. In the high security area, Inner Mongolia is fully self-sufficient and the technical amount of renewable energy ranks first among provinces. Further, its production of renewable energy is well-consumed, so the system risk is low. According to the hazard and vulnerability situation of the integrated energy system, the risk situation and causes of each province could be traced.
From the perspective of geographical distribution, the study combined ArcGIS and the provincial energy system risk results to obtain a corresponding spatial distribution profile as shown in Figure 7. In terms of hazard, the eastern region is more dangerous than the central and western regions, and the overall decline is from the east to the west. In terms of vulnerability, the central region is more vulnerable because of the close grid-connection. Figure 7c visually presents the results of the risk assessment of the 31 Chinese provinces under the combined effects of hazard and vulnerability.

Conclusions
The security of regional integrated energy systems has great significance for improving energy utilization efficiency and promoting large-scale development of renewable energy. This paper comprehensively addressed energy security issues at the regional level, which is vital for systematic and large-scale energy development and was not sufficiently addressed previously. Because of the lack of a systematic approach for regional energy security, natural disaster risk assessment was introduced to establish a comprehensive assessment framework of regional energy systems. The empirical study of Chinese provinces quantitatively measured the hazards and vulnerability with 5 dimensions and 17 metrics. The regional inter-relationship and interactions were highlighted and fully considered through a spatial method. The proposed method can solve the problem of energy security with renewable energy and grid connection, and help evaluate energy security for different regions or countries.
First, a spatial weight matrix enabled us to portray the connection aggregation relationship among regions. The results of the Chinese provinces showed that renewable energy consumption shows a significant aggregation, that is, the grid connection among each province has a positive effect on renewable energy development. From the perspective of each region, the Northeast, Central China and East China regions are typical clustering regions. There intra-and inter-regional grids are closely connected, and the grid connection has a strong ability to cope with risks. Therefore, the spatial weight matrix is worthy of further discussion on grid connectivity to realize integrated development.
Second, the security assessment framework takes full account of the renewables and grid connection, and regards the regional energy systems as organic units. The natural disaster risk assessment method demonstrated its applicability to break down energy system risks into hazard and vulnerability aspects, which helped to analyze and explain systemic risks in a traceable way. In the objective weighting process, transmission had the most significant impact, which should be paid more attention in further energy research.
Third, a subdivisions study was conducted with a hazard index and vulnerability index, based on which suggestions could develop and a more targeted analysis of risk in various areas can be carried out. Further, the causes of risks can be circumvented and weakened, and the ability of the system to resist risks can be improved. The 31 Chinese provinces were divided into four subdivisions based on the results of risk measurement. In terms of hazards, the eastern region exhibits more than the central and western regions, and the overall decline is from the east to the west. The analysis results are in good agreement with practical experience, and would help with systematic evaluation and make up for relative research from the perspective of resources and transmission.
Overall, the assessment framework proposed in this paper provides a general tool for the systematic evaluation of regional energy security issues. In this paper, the model was specifically applied to 31 Chinese provinces. It should be noted that China has some typical features, such as its energy security has relatively stronger support, and the grid connection and operations are tightly associated among provinces. Nevertheless, the research framework could also be applied to other regions or nations around the world. Based on the different features in other regions, different results may be produced and compared. The assessment method could be regarded as a start for using a spatial method in systematic energy security evaluation and a comprehensive framework for various energy resources, in order to support optimal and intelligent national or regional energy management.

Conflicts of Interest:
The authors declare no conflict of interest.