Risk Assessment of Distribution Lines in Typhoon Weather Considering Socio-economic Factors

: In recent years, the frequent occurrence of typhoon weather has posed a signiﬁcant threat to the stable operation of the distribution network in the southeastern coastal areas of China. Ensuring the safety of distribution lines is crucial for the normal functioning of the distribution network. Therefore, this paper proposes a risk assessment method for distribution lines in typhoon weather. Firstly, the risk assessment system for distribution lines is constructed by considering three perspectives: line structure, line state, and social economic factors. Secondly, the weight of each evaluation index is calculated using the analytic hierarchy process and CRITIC weight method. The cooperative game method is then employed to combine the calculation results, and the results are further optimized using variable weight theory. Finally, a cloud model-based risk assessment model for distribution lines is established. The analysis and calculation of distribution network data in a speciﬁc area indicate that the risk assessment level, which takes into account social and economic factors, is more accurate compared to other methods discussed in this paper. It is observed that the multi-model approach yields higher accuracy than the single-model approach. Therefore, the proposed method holds signiﬁcant reference value for evaluating the risk level of distribution lines.


Introduction
Typhoons pose a significant security threat to coastal power grids, resulting in severe social and economic losses [1,2].To mitigate the impact of typhoon disasters, conducting advanced risk assessments and implementing appropriate measures on coastal distribution networks is crucial.Power outages during typhoon events are primarily caused by cascading failures in distribution lines [3].Hence, conducting risk assessments on distribution lines under typhoon conditions plays a vital role in preventing distribution network disasters.
Several scholars have conducted research on the evaluation index system of distribution line vulnerability [4].Ref. [5] summarizes the application of the maximum flow theory method and the improved betweenness method in identifying vulnerable lines based on complex network theory.In Ref. [6], the line's apparent power is used as the line flow, and the electrical in-degree and out-degree centrality are proposed to identify vulnerable lines.However, the influence of the nodes at both ends of the line on the line vulnerability assessment is not considered.Ref. [7] establishes a comprehensive line vulnerability index based on grid topology information, taking into account the node information at both ends of the line.In Refs.[8,9], an improved transmission betweenness method is proposed to quantify the vulnerability of transmission lines.This method combines edge betweenness and power transmission distribution factor.Ref. [10] extends the traditional structural vulnerability index by considering factors such as the overall load rate of the system and

Line Degree
The degree is a measure of the importance of the connection between a node and other nodes in the topology.In the case of distribution lines, the degree of importance of the line is higher when the nodes at both ends of the line are more important.
Based on the traditional node degree, this paper considers the node input power of the distribution line [12], so the line degree is defined as:  The degree is a measure of the importance of the connection between a node and other nodes in the topology.In the case of distribution lines, the degree of importance of the line is higher when the nodes at both ends of the line are more important.
Based on the traditional node degree, this paper considers the node input power of the distribution line [12], so the line degree is defined as: Energies 2023, 16, 6664 3 of 15 D i represents the number of edges connected to the node i; V represents the set of all nodes adjacent to the node i; P i indicates the actual power flowing into the node i.

Line Betweenness
To assess the power flow propagation in the line, we consider the law of power flow propagation between nodes.In this study, we utilize the following expressions for calculation: In the equation, ω ij (k) is the current value transmitted by line k when a unit current source is applied between the bus node i and the load node j, assuming that the active power transmitted from the node m to the node n is P mn (i, j).
The traditional betweenness calculation, which relies on the shortest path, is not applicable to distribution networks due to their radial open-loop operation.Therefore, it is more practical to use ω ij (k) instead of the shortest path.

Line Failure Rate
The failure rate of distribution lines under typhoon weather is strongly influenced by the typhoon wind field.In this study, we utilize the enhanced Batts wind field model [18] to estimate the wind speed along the distribution line.The wind speed at a specific point is determined by the geographical location of the typhoon center and the point itself.The relationship between the probability of a single tower failure in the distribution line and the typhoon wind speed can be represented by an exponential curve function [19]: (3) Among them, V z is the wind speed value of the typhoon at a certain moment, and V min is the design wind speed of the tower.According to the design wind speed standard of the 35 kV distribution line, the value of this paper is 30 m/s, and the value range of the parameter K is (0, 0.4).
Several towers form a distribution line.According to the reliability mathematical model of the series system, the fault probability of a distribution line k is [18]:

Line Loss Value
In the identification process of vulnerable distribution lines, the electricity sales loss caused by the fault and its consequences are multiplied to evaluate the vulnerable lines, that is: In the equation E k is the loss value of the line k caused by the typhoon disaster.U k represents the unit loss cost when the fault occurs in the line k, and Q k represents the load loss of the fault load node caused by the fault of the line k.According to the importance of load nodes, the economic losses of average unit power can be divided into two categories, namely 0.79 CNY per kilowatt hour and 0.57 CNY per kilowatt hour.

Population Size
The failure of distribution lines during typhoon weather can result in power outages within a specific time range in the region.This can significantly impact the lives of the population in the area and may even lead to casualties.This indicator highlights the significance of routes in the region based on population density.The importance of a line passing through an area is directly proportional to the population density, and a higher weight value is assigned to such lines.

Industrial Output Value
In order to mitigate economic losses, enterprises in the region may need to halt production during typhoons, which can have an impact on the region's economic development to some extent.The industrial output value serves as an indicator of the significance of the lines passing through the region in terms of industrial production.A higher industrial output value indicates greater importance of the line passing through the region and, therefore, assigns a higher weight value to it.

Gross Production
In addition to industrial production, the production of agriculture and the tertiary industry may also be affected by certain procedures during typhoon weather.The impact of typhoons varies depending on the industrial structure of different regions.The gross domestic product (GDP) is a reflection of the significance of the economic route passing through a region.The higher the economic aggregate, the more important the route becomes, and thus, it is assigned a greater weight value.

Integrated Assessment Model
In this paper, we combine the subjective analytic hierarchy process and CRITIC weight method to comprehensively determine the weight coefficients of the seven indicators.We utilize the cloud model theory to construct the evaluation model and compare the standard cloud model with the comprehensive cloud model to determine the final risk level.The framework of the risk assessment model is depicted in Figure 2. the line passing through the region and, therefore, assigns a higher weight value to it.

Gross Production
In addition to industrial production, the production of agriculture and the tertiary industry may also be affected by certain procedures during typhoon weather.The impact of typhoons varies depending on the industrial structure of different regions.The gross domestic product (GDP) is a reflection of the significance of the economic route passing through a region.The higher the economic aggregate, the more important the route becomes, and thus, it is assigned a greater weight value.

Integrated Assessment Model
In this paper, we combine the subjective analytic hierarchy process and CRITIC weight method to comprehensively determine the weight coefficients of the seven indicators.We utilize the cloud model theory to construct the evaluation model and compare the standard cloud model with the comprehensive cloud model to determine the final risk level.The framework of the risk assessment model is depicted in Figure 2.

AHP Method
The Analytic Hierarchy Process (AHP) [20] is a multi-objective decision analysis method that combines qualitative and quantitative analysis.It utilizes pairwise comparison to assess the judgment of experts and determine the weight of each element.The calculation steps of the AHP are as follows [21]: (1) According to the expert experience, the nine-scale method is used to construct the judgment matrix:

AHP Method
The Analytic Hierarchy Process (AHP) [20] is a multi-objective decision analysis method that combines qualitative and quantitative analysis.It utilizes pairwise comparison to assess the judgment of experts and determine the weight of each element.The calculation steps of the AHP are as follows [21]: Energies 2023, 16, 6664 5 of 15 (1) According to the expert experience, the nine-scale method is used to construct the judgment matrix: Among them, b ij represents the importance of indicator i to indicator j.
(2) Hierarchical single sorting.After passing the consistency test of the judgment matrix, the normalized feature vector is used as the weight vector of this level: n is the index number of this layer.λ max is the maximum eigenvalue of the judgment matrix.
In the equation, CR is the random consistency ratio of the comparison matrix; CI is the general consistency index; RI is the random average consistency index; when CR < 0.1, the judgment matrix is considered to meet the consistency.If CR ≥ 0.1, the matrix needs to be modified until the condition is satisfied.
(3) Hierarchical total sorting is used to calculate the weight vector of the scheme layer to the target layer and test its consistency.The decision is made based on the total ranking weight vector, taking into account the weight of each index of the analytic hierarchy process:

CRITIC Method
The CRITIC weight method [22] is used to comprehensively measure the weight of evaluation indicators based on their volatility and conflict and represents the amount of information using the product of the two.The specific steps of the CRITIC weight method are as follows [23]: (1) Using data standardization processing: (2) Calculate the variability quantitative index of the index j and other indicators.The standard deviation is used to represent the contrast of the four evaluation indexes of line failure rate, line loss value, line betweenness, and line degree.The calculation equation is shown in Equation (11).The study uses 32 lines as the evaluation object, so the value range of parameter i is [1,32].
where x j is the mean value of the index j, and S j is the standard deviation of the index j.
(3) Seven evaluation indicators have been established above to further calculate conflicting quantitative values between them.Calculate the conflict quantitative value between the index j and other evaluation indexes.The calculation equation is: In the equation, r ij is the correlation coefficient between the index i and the evaluation index j.cov(Z i , Z j ) denotes the covariance between columns i and j of the standard matrix Z; σ j is the standard deviation of each index.
(4) Calculate the amount of information contained in the indicator j.
(5) Calculate the weight ω j of the index j.

Cooperation Game Model
By minimizing the deviation between the combination weight and the basic weight [24], the combination weighting is more reasonable.Two groups of basic weights are obtained by using analytic hierarchy process and CRITIC weight method.According to the method of cooperative game [25], the combination coefficients λ 1 * and λ 2 * are solved.Finally, the coefficients are substituted into the combination weight W j * .
The objective function is as shown in Equation ( 15): According to the differential principle, the solution of the objective function is transformed into solving linear equations: The obtained linear combination coefficients λ 1 and λ 2 are standardized: Further, the combined weight coefficient is:

Variable Weight Theory
The comprehensive weight obtained from the above calculation is a constant weight, which remains unchanged once determined for each index.However, the power flow distribution of distribution lines changes with the variation in typhoon wind speed.This change becomes particularly significant when individual indicators reach critical values.In such cases, the constant weight fails to accurately assess the risk of the region, resulting in unreasonable assessment results.To address this limitation, this paper adopts the variable weight theory to reflect the balance of the state of each factor in the comprehensive evaluation [26].The calculation equation is as follows: In the equation, ω i v is the variable weight of the index i; W j * is the index j combination weight; x i is the index value i; ∂ is an equilibrium function, and its value depends on the relative importance of each comprehensive state quantity.The value range is [0, 1].The smaller the value, the higher the attention to the corresponding state quantity.[27].For the line structure and line state index, it can reflect the operation state of the distribution line and the empirical value ∂= 0.3; for the social factor index, it can reflect the economic loss and the number of people affected by the typhoon when the distribution line is affected by the typhoon, which needs to be paid attention to, and the empirical value ∂= 0.2.

Cloud Model 4.1. Definition of the Cloud Model
The cloud model can realize the mutual transformation between qualitative concept and quantitative representation, and its properties are represented by three digital features: expectation Ex, entropy En, and hyper entropy He [28].The expectation En is the center of the cloud, and the entropy En represents the degree of dispersion of the qualitative concept, which further reflects the relationship between the randomness and fuzziness of the qualitative concept [29].He is the uncertainty measure of entropy.Taking the cloud model with digital feature (1, 0.09, 0.007) as an example, a one-dimensional cloud model of 5000 cloud droplets is constructed, as shown in Figure 3.
the balance of the state of each factor in the comprehensive evaluation [26].The calculation equation is as follows: In the equation, ω v i is the variable weight of the index i; * Wj is the index j combination weight; xi is the index value i; ∂ is an equilibrium function, and its value depends on the relative importance of each comprehensive state quantity.The value range is [0,1].The smaller the value, the higher the attention to the corresponding state quantity.[27].For the line structure and line state index, it can reflect the operation state of the distribution line and the empirical value =0.3 ∂ ; for the social factor index, it can reflect the economic loss and the number of people affected by the typhoon when the distribution line is affected by the typhoon, which needs to be paid attention to, and the empirical value =0.2 ∂ .

Definition of the Cloud Model
The cloud model can realize the mutual transformation between qualitative concept and quantitative representation, and its properties are represented by three digital features: expectation Ex , entropy En , and hyper entropy He [28].The expectation En is the center of the cloud, and the entropy En represents the degree of dispersion of the qualitative concept, which further reflects the relationship between the randomness and fuzziness of the qualitative concept [29].He is the uncertainty measure of entropy.Taking the cloud model with digital feature (1,0.09,0.007)as an example, a one-dimensional cloud model of 5000 cloud droplets is constructed, as shown in Figure 3.  Suppose that U is used to represent a universe of discourse described by quantitative values; C(Ex, En, He) is a qualitative concept described by qualitative language in U space.If U has a quantitative value x, x is a random realization of C, then the membership degree of x to C is expressed as [30]: The forward cloud generator is a mapping from qualitative concepts to quantitative representations, while the reverse cloud generator is a transformation from quantitative concepts to qualitative concepts.Figure 4 illustrates the cloud generator.
The forward cloud generator is a mapping from qualitative concepts to quantitative representations, while the reverse cloud generator is a transformation from quantitative concepts to qualitative concepts.Figure 4

Establish a Standard Cloud Model
The risk level of distribution lines is divided into 'extremely low risk', 'low risk', 'medium risk', 'higher risk', and 'highest risk'.According to the golden ratio segmentation method, the comments are graded [30], and the numerical characteristics of the evaluation standard cloud are calculated, as shown in Table 1.The standard cloud model is shown in Figure 5.

Establish a Standard Cloud Model
The risk level of distribution lines is divided into 'extremely low risk', 'low risk', 'medium risk', 'higher risk', and 'highest risk'.According to the golden ratio segmentation method, the comments are graded [30], and the numerical characteristics of the evaluation standard cloud are calculated, as shown in Table 1.The standard cloud model is shown in Figure 5.To address the challenge of distinguishing the relative position between the comprehensive evaluation cloud and the standard cloud, this study proposes a method to calculate the similarity and obtain the final risk evaluation level.The vector of the standard cloud is expressed as (Ex , En , He ) , and the vector of the integrated cloud is expressed as (Ex , En , He ) . The similarity is calculated according to the following equation [31], and sim is the similarity function.According to the principle of maximum similarity, the risk level is finally evaluated.sim( , ) In order to simplify the calculation results of the similarity function, Di is To address the challenge of distinguishing the relative position between the comprehensive evaluation cloud and the standard cloud, this study proposes a method to calculate the similarity and obtain the final risk evaluation level.The vector of the standard cloud is expressed as → F i = (Ex i , En i , He i ), and the vector of the integrated cloud is expressed as → F j = (Ex j , En j , He j ).The similarity is calculated according to the following equation [31], and sim is the similarity function.According to the principle of maximum similarity, the risk level is finally evaluated.
Energies 2023, 16, 6664 In order to simplify the calculation results of the similarity function, Di is used in this paper.D 1 ~D5 is used to represent the similarity between the risk value of the comprehensive cloud and the standard cloud from extremely low to extremely high.

Experimental Model and Data
The landing path and time of Typhoon Muifa [32], as well as the layout of the distribution line system [33], are shown in Figure 6, using 'Severe Typhoon Muifa' as an example.The maximum wind speed, central pressure, moving speed, and landing path data of the typhoon center were obtained from the typhoon network of the Central Meteorological Observatory.The population, total area, GDP, and total industrial output value of the study area were obtained from the statistical yearbook of the city (2021) and the statistical yearbook of each county [34].The data were classified into ten grades using the natural discontinuity point classification method, with an assignment interval of [0, 1].The weight value increases as the color gets darker.Figure 7 shows the population density of the region, Figure 8 shows the industrial output value, and Figure 9 shows the GDP.The population, total area, GDP, and total industrial output value of the study area were obtained from the statistical yearbook of the city (2021) and the statistical yearbook of each county [34].The data were classified into ten grades using the natural discontinuity point classification method, with an assignment interval of [0, 1].The weight value increases as the color gets darker.Figure 7 shows the population density of the region, Figure 8 shows the industrial output value, and Figure 9 shows the GDP.The population, total area, GDP, and total industrial output value of the study area were obtained from the statistical yearbook of the city (2021) and the statistical yearbook of each county [34].The data were classified into ten grades using the natural discontinuity point classification method, with an assignment interval of [0, 1].The weight value increases as the color gets darker.Figure 7 shows the population density of the region, Figure 8 shows the industrial output value, and Figure 9 shows the GDP.

Weight Calculation
(1) The calculation results of subjective and objective weighting methods.
According to the equation, the results of the two weighting methods are calculated as shown in Table 2.

Weight Calculation
(1) The calculation results of subjective and objective weighting methods.According to the equation, the results of the two weighting methods are calculated as shown in Table 2.

Weight Calculation
(1) The calculation results of subjective and objective weighting methods.
According to the equation, the results of the two weighting methods are calculated as shown in Table 2. (2) Calculate the combined weight.
According to the Equations ( 15)-( 17), the combined weight coefficient is calculated as: According to Equations ( 18) and ( 19), the combined weight coefficients of the four evaluation indexes are obtained as shown in Table 3.The weight results of the factor layer and the criterion layer are shown in Figure 10.
According to the Equations ( 15)-( 17), the combined weight coefficient is calculated as: According to Equations ( 18) and ( 19), the combined weight coefficients of the four evaluation indexes are obtained as shown in Table 3.

The Cloud Model of Each Index Is Constructed
The calculation results are shown in Table 4.

The Cloud Model of Each Index Is Constructed
The calculation results are shown in Table 4. Figure 12 illustrates a comparison diagram between the cloud model of the criterion layer and the standard cloud model based on the numerical characteristics of the criterion layer shown in Table 4.

The Cloud Model of Each Index Is Constructed
The calculation results are shown in Table 4.

Normal Level Index
Digital Features of Criterion Layer line structure B1 (0.6848, 0.1240, 0.0295) line status B2 (0.6529, 0.1113, 0.0045) socio-economic factor B3 (0.3917, 0.0760, 0.0560) Figure 12 illustrates a comparison diagram between the cloud model of the criterion layer and the standard cloud model based on the numerical characteristics of the criterion layer shown in Table 4.By applying the criterion layer and the variable weight coefficient, the digital eigenvalues of the comprehensive cloud model can be determined.The cloud model representing the target layer is illustrated in Figure 13.The characteristic value of the comprehensive cloud model is (Ex, En, He ) = (0.057, 0.1052, 0.0354).In order to further accurately judge the risk level of the distribution line, the similarity between the cloud models is calculated by Equation ( 21).The similarity calculation results of the comprehensive cloud model and each standard cloud In order to further accurately judge the risk level of the distribution line, the similarity between the cloud models is calculated by Equation (21).The similarity calculation results of the comprehensive cloud model and each standard cloud model are as follows: D 1 ~D5 , the calculation results are D 1 = 0.1867, D 2 = 0.9780, D 3 = 0.9935, D 4 = 0.9949 and D 5 = 0.9909.According to the calculation results, D 4 > D 3 > D 5 > D 2 > D 1 .The evaluation results indicate a higher risk level, as seen in the comprehensive cloud model diagram.Therefore, early warning measures should be implemented for the distribution line.To further validate the effectiveness of the method proposed in this paper, the AHP method and CRITIC method are combined with the cloud model separately, resulting in a medium risk assessment level.The cloud model method, which considers social and economic factors and variable weight theory, aligns better with the actual situation.The comparison results of the model evaluation are presented in Table 5.

Discussion
This paper first establishes seven distribution line vulnerability evaluation indexes from the three aspects of line structure, line status, and socio-economic factors and uses the method of cooperative game-variable power theory to optimize the assignment results of the subjective assignment method and objective assignment method and, finally, uses the cloud model to realize the assessment of the risk level.The experimental results show that when a typhoon occurs, the risk level of the region is 'higher risk', which is consistent with the actual situation and can truly and objectively reflect the risk level of the distribution lines in the region.
First, compared with existing studies, this paper takes into account the impact of socio-economic factors on risk assessment when establishing evaluation indicators.When a typhoon occurs, if a line is located in a densely populated area with high economic output value, it will bring serious casualties and economic losses if it is disturbed by the typhoon; therefore, even if the structure and operation of the line are at a medium level, it should be given some attention; if the structure and operation of the line are in a dangerous state, it is more important to combine the socio-economic factors to make an accurate risk level assessment.Therefore, it is one-sided to analyze only the structure and operation status of distribution lines.
Secondly, this paper uses cooperative game combined with the variable weighting theory method to optimize the basic assignment result.This paper adopts the cooperative game method to consider the imbalance between the data, avoiding the subjectivity of the traditional linear weighting method; on this basis, this paper adopts the variable weighting theory to optimize the result, taking into account the problem of sudden changes in the distribution of tidal currents of the distribution line due to the typhoon, and gives different values of equilibrium function to the evaluation indexes to further optimize the experimental results.
Finally, through the comparative analysis of the single evaluation model and the comprehensive model proposed in this paper, it can be concluded that both the method of this paper and the traditional method consider that the distribution line in this area is risky, but the results of the model in this paper are more in line with the actual situation, and it is considered that the distribution line has a higher risk level, as shown in Table 5.
The research in this paper also has certain limitations, not fully considering the impact of a small area of building shading on the risk assessment of distribution lines, and it is hoped that fieldwork can be carried out so that more accurate data can be obtained and a more accurate risk assessment of the study area can be realized.

Conclusions
This paper presents a distribution line risk assessment model based on variable weight and cloud models to address the shortcomings of the current evaluation index and model.The proposed model takes into account not only the physical structure of the distribution line but also the social and economic factors of the area where the distribution line is located as evaluation indices.By incorporating both power and social-economic systems, the evaluation model can effectively assess the risk level of the distribution line.The cloud model enables the conversion of qualitative and quantitative concepts, facilitating the visualization of risk levels in distribution lines.The results of the example show that the proposed method can realize the judgment of the risk level of the distribution line and can be used to guide the typhoon disaster prevention of the distribution network in the coastal area.

Figure 1 .
Figure 1.Risk assessment index system of distribution network in typhoon weather.

Figure 1 .
Figure 1.Risk assessment index system of distribution network in typhoon weather.

Figure 5 .
Figure 5. Risk assessment standard cloud model diagram.

Figure 5 .
Figure 5. Risk assessment standard cloud model diagram.

Figure 7 .Figure 6 .
Figure 7. Population density map of a certain area.

Figure 7 .
Figure 7. Population density map of a certain area.

Figure 8 .
Figure 8. Map of the distribution of total industrial output value in the region.

Figure 9 .
Figure 9. Map of the distribution of GDP in the region.

Figure 8 .
Figure 8. Map of the distribution of total industrial output value in the region.

Figure 9 .
Figure 9. Map of the distribution of GDP in the region.

Figure 9 .
Figure 9. Map of the distribution of GDP in the region.
of the factor layer and the criterion layer are shown in Figure 10.

Figure 10 .
Figure 10.Evaluation index layer weight comparison result diagram.The criterion layer weight comparison result diagram are shown in Figure 11.

Figure 10 . 18 Figure 11 .
Figure 10.Evaluation index layer weight comparison result diagram.The criterion layer weight comparison result diagram are shown in Figure 11.Energies 2023, 16, x FOR PEER REVIEW 14 of 18

Energies 2023 ,
16, x FOR PEER REVIEW 15 of 18 By applying the criterion layer and the variable weight coefficient, the digital eigenvalues of the comprehensive cloud model can be determined.The cloud model representing the target layer is illustrated in Figure 13.The characteristic value of the comprehensive cloud model is (Ex,En,He) 0.057 0.1052 0.0354 = ( ， ， ) .

Table 1 .
Eigenvalues of standard cloud model.

Table 1 .
Eigenvalues of standard cloud model.

Table 2 .
Subjective and objective weight results.Map of the distribution of total industrial output value in the region.

Table 2 .
Subjective and objective weight results.

Table 2 .
Subjective and objective weight results.

Table 3 .
Constant weight and variable weight coefficient.

Table 3 .
Constant weight and variable weight coefficient.

Table 4 .
Numerical characteristics of the criterion layer.

Table 4 .
Numerical characteristics of the criterion layer.

Table 4 .
Numerical characteristics of the criterion layer.

Table 5 .
Model comparison results.
Author Contributions: Methodology, A.H. and X.F.; formal analysis, X.F.; resources, S.S.; data curation, A.H.; writing-original draft preparation, X.F.; writing-review and editing, A.H. and X.F.; visualization, S.S.; supervision, D.H.; funding acquisition, F.Z. and D.H.All authors have read and agreed to the published version of the manuscript.This research was funded by the Shanghai Municipal Science and Technology Commission Local Colleges and Universities Capacity Building Project, grant number 20020500700; this research was funded by the Huaneng Group Headquarters Science and Technology Project "Research on Site Selection and Support Technology of Offshore Wind Power Based on Deep Sea", grant number HNKJ20-H66.