#### Abstract

For long-span cable-stayed bridges, cables serve as one of the most important components to guarantee structural integrity. Forces of stay cables indicate not only the performance of cables themselves but also the overall condition of bridges. In order to help stakeholders to make maintenance decisions, an extreme cable force estimation method was proposed based on cable force measurements and traffic data from the weighing system. First, raw monitoring data were preprocessed based on a median filtering to obtain usable cable force signals. The multiresolution wavelet method was used to extract traffic-induced force component from mixed signals. Then, a Monte Carlo-based random vehicle model was developed using traffic data from the weighing system. Based on field temperature measurements and simulation of traffic-induced effects, extreme cable forces with respect to vehicle loads and temperature effects were predicted by extreme value theory. The Generalized Pareto Distribution (GPD) was adopted to establish the probability distribution models of the daily maximum cable force. Then, the extreme value within a return period of 100 years was determined and compared with the design loading demand. Finally, the effectiveness of the proposed method was validated through a cable-stayed bridge in China. As a result, the low-frequency varying component of cable force response had positive correlation with environmental temperatures, and the extreme value of the predicted cable force under prospective traffic volumes was within limit interval value according to the design code. The conclusions can be utilized by bridge owners to make maintenance decisions.

#### 1. Introduction

With rapid growth of infrastructure investment in recent years, more and more long-span cable-stayed bridges have been constructed over rivers and valleys [1]. Meanwhile, traffic volumes of bridges have increased significantly compared to those when bridges were initially opened to public traffics, especially in some emerging countries [2–4]. For instance, China’s annual growth rate of freight traffic volume approached 6% in 2018 [5], which was a critical parameter for evaluation of structural performances. Structural responses might exceed the designed ones in the following several decades owing to the pressure of rapid increase in traffic volumes [6]. This phenomenon will inevitably impact operational and structural safety of long-span cable-stayed bridges. In this regard, measures should be taken to investigate the actual responses of bridge components (stay cables, towers, girder, etc.). Moreover, further prediction is necessary based on existing information and extreme value estimation methods. Nowadays, structural health monitoring (SHM) systems have been installed in plenty of long-span cable-stayed bridges [7–9]. With the assistance of SHM systems that consist of various sensors, real-time responses of cable-stayed bridges under different actions can be recorded, so a large amount of monitoring data has been accumulated for further discussion. It is possible to take advantages of these monitoring data to investigate the actual loading effect and to evaluate and predict the performance of structural components [10–12].

For cable-stayed bridges, the stay cable is regarded as one of the most important components to ensure the allowable displacement and distribution of bending moments along the bridge deck with prestressing force [13]. Sensors installed in cables are used to collect real-time signals in a SHM system of cable-stayed bridges. Then, the cable forces can be obtained directly or can be calculated from indirect measurements such as frequency or acceleration [14]. It is primary to analyze the characteristics of recorded cable force signals under actual loading conditions. Signals induced by temperature actions and vehicle loadings are two main ingredients in the mixed signal, namely, slow-varying ingredient and dynamic ingredient, respectively [15–17]. To extract different components, blind separation methods are usually adopted [18]. Xu et al. [19] used a multiresolution wavelet-based method to separate thermal effects from bridge responses based on the distinguished frequency bandwidths. A practical multivariate linear-based model was also presented to simulate and separate thermal effects from the cable-stayed bridge response [20]. Kromanis and Kripakaran [21] suggested novel data-driven methods referred to as Regression-Based Thermal Response Prediction and Traffic-Induced Response Prediction to predict thermal and vehicular response. The proposed methods were evaluated by experimental data collected from a laboratory truss. Zhou and Sun [22] established a sound understanding of thermal response in cable-stayed bridges and found that the thermal response of cables was mainly determined by environmental temperature. In addition, thermal response of cables subject to fire was also studied [23, 24].

To predict the performance of cables, the random vehicle model should be established to simulate traffic-induced cable force responses under prospective traffic volumes. Previously, the passing vehicles were only tracked and identified by monitoring camera system on bridges [25]. Now, the weigh-in-motion (WIM) system is widely used in long-span bridges for ascertaining vehicle parameters, including total weight, axle distance, speed, and driving lane [26]. In this field, Ma et al. [27] obtained basic vehicle information based on traffic data analysis, and the statistical distribution functions of gross vehicle weight under different loading conditions and wheelbase for each type of vehicle were also analyzed to develop the vehicle model. Wang et al. [28] proposed a vehicle model for long-span bridges in the framework of Eurocode Load Model 1, in which multiple lane factors and impact factors were taken into account. Zhou et al. [29] used micro- and macrosimulation with different load model varieties of traffic loading in different bridge regions to propose an efficient multiscale traffic modelling approach. Guo et al. [30] developed a vehicle load model based on weigh-in-motion data, which considered the uncertainties associated with the number of axles, axle weights, axle spacing, and transversal position of vehicles. Then, the model was used to evaluate the fatigue reliability of steel bridge details. It is also a significant work to forecast future traffic volumes based on current traffic volumes as both seasonality and growth trend can be included [31]. However, there are rare studies published in evaluating prospective performance of cables combining monitoring data with extreme value calculated from random vehicle model [32].

In this paper, the cable forces of cable-stayed bridges were discussed and the extreme value was predicted based on monitoring data and the developed random vehicle model. First, a median-filtering-based preprocessing was conducted to obtain usable cable force measurements. The multiresolution wavelet method was used to extract traffic-induced force. Meanwhile, the Monte Carlo based random vehicle model was developed using traffic data from bridge weighing system. Then, the extreme cable forces with respect to vehicle loads and temperature effects for the return period of 100 years were predicted by applying the Generalized Pareto Distribution (GPD) model. Finally, the reliability of the proposed method was validated by a case study of a long-span cable-stayed bridge, the third Nanjing Yangtze River Bridge in China.

#### 2. Methodology

At present, monitoring systems are widely installed in long-span cable-stayed bridges. The monitoring data, including cable force measurements from SHM system and traffic data from weighing system, can achieve continuous collection. However, further investigation is still needed to obtain valuable information from mass data. Considering the probable performance of stay cables within the return period (which can be regarded as an approximate frequency of occurrence [33]), an analysis and estimation method was developed based on current monitoring data and the average daily traffic, as shown in Figure 1. Firstly, the raw data was preprocessed. Then, effectiveness of the developed random vehicle model was verified by the separated vehicle-induced cable forces. Finally, the GPD model was adopted to predict extreme values of cable forces in the future.

##### 2.1. Processing for Monitoring Data

At present, the recorded cable force data by SHM systems includes direct and indirect measurements. Compared with the common vibration-based measurements, the direct force measurements help to know the operational condition conveniently. The influence of conversion error was also eliminated to ensure data accuracy. However, direct measuring methods skip the transforming procedure such as fast Fourier transform (FFT), which can reduce signal noise simultaneously. Therefore, it is necessary to perform preprocessing to eliminate the influence of data noise before analyzing the direct cable force measurements. Figure 2(a) shows the typical cable force measurements collected by sensors. It indicates that the signals are considerably affected by salt-and-pepper noise, which is caused by the uncertainty of sensor measurements itself and transmission channels as well as the decoding process. In this paper, the median filtering was used to deal with the raw signals, and the output of filter is denoted as follows:where deontes the output signals, is the filtering function, is the input signals, and *N* is the selection number to ensure the window width of 2*N* + 1. The basic theory of median filtering is to replace the value of a certain point in the data series with the median value of the adjacent several points. For instance, Figure 2(b) shows the results of median filtering with a window width of 5. Thus, the signals were filtered without the effects of noise.

**(a)**

**(b)**

To simplify, the monitoring cable forces in bridge service period are considered as the linear superposition of temperature effect and dead load-induced force and vehicle load-induced force, corresponding to slow-varying ingredient and dynamic ingredient in signal, respectively. Compared with traffic load, the wind effect to cable force is ignored in this paper for it does not significantly influence the results, which was also validated by existing research [34, 35]. The signal is denoted as follows:where is the monitoring cable force signal, is the total cable force induced by environmental temperature effect and dead load, is the cable force induced by vehicle load, denotes the environmental temperature measurements, and *k* and *c* are constant parameters.

To separate signals of aforementioned two ingredients, the multiresolution wavelet method based on distinguished frequency bandwidths is adopted [19]. The core idea of wavelet method is to use a finite length or fast decay wave to reconstruct signals. A time-domain signal can be described as the convolution of a wavelet coefficient and a wavelet function . After mathematical transformations, the wavelet function can be expressed as follows:where *a* and *b* are stretch parameter and translation parameters. The wavelet coefficient explains the similarity between signals and wavelet function , which is

The multiresolution wavelet method transforms the signals to approximation and detail coefficients that depict the frequency distribution in the time and frequency domains, as shown in Figure 3. Ai denotes the *i*^{th} term of approximation, while Di denotes the *i*^{th} term of detail. Thus, the original signals can be separated. For instance, Figure 4 shows the one-day monitoring cable force data for a typical stay cable with a sampling frequency of 1 Hz. The wavelet separation was conducted with Db12 wavelet to obtain the slow-varying ingredient with low frequency, namely, the total cable force induced by environmental temperature and dead load. The rest was regarded as the high-frequency ingredient, that is, the total cable force induced by vehicle load.

##### 2.2. Establishment of the Random Vehicle Model

Similar to other codes, the Chinese codes mention that the vehicle load should be modified based on the actual situation during the assessment [36]. The mentioned modified factors that significantly influence the extreme value of vehicle-induced cable forces include the annual average daily traffic (AADT) and the proportion of heavy vehicles. In the early operational period of bridges, the traffic volume gradually increased but usually would not reach saturation. Thus, the accumulated data in early period should not be directly used to predict the extreme cable forces in the future. Therefore, the traffic growth is often considered a critical factor of the future traffic effect [37]. In this paper, a random vehicle model was developed to predict the extreme value of vehicle loads with prospective traffic volume. To establish the random vehicle model, the major vehicles in operation were divided into different types according to axle parameters based on the weighing data and monitoring video. The maximum likelihood estimation and expectation maximization (EM) algorithm were used to build the probability models of major vehicle parameters, such as vehicle type, total weight, axle weight, arrival time, wheelbase, driving lane, and speed.

Then, the establishment procedure for dynamic random vehicle model was developed in several stages, as shown in Figure 5. Firstly, traffic volume and simulation duration were determined. Vehicle arrival time and number were also generated via Poisson processing. Then, the probability models of vehicle parameters were utilized to develop a Monte Carlo based vehicle database, which was regarded as inputs of vehicle-following program, hence the dynamic traffic flow. Finally, the random vehicle model was acted on structural influence lines to calculate the time history of vehicle-induced responses.

##### 2.3. Extreme Value Analysis and Estimation

To realize extreme cable force estimation in the field of statistics for the return period, it is required to estimate the extreme cable forces from current records [38]. As a supplement to the normal distribution, the extreme value method was proposed to predict confidence limits on extreme cable forces [39]. In view of the complicacy of the signals, two preprocessed ingredients were analyzed, respectively. Assume that these two ingredients data series (slow-varying ingredient is analyzed seasonally) were stationary, independent, and in accord with GPD model. The common methods to obtain the extreme value of stationary data series include interval maximum value and peaks over threshold (POT). Due to the insufficiency of the accumulated data, the assumptions of interval maximum model were difficult to meet. Therefore, the POT model described by the GPD was eventually utilized to estimate the return level of cable force [40, 41]. For a stationary series , …, , is denoted aswhere is the selected threshold, and is a subsample indicating that exceeds . Then, at least subsequent *r* consecutive data are lower than , where *n* is the total amount of data. Based on the GPD model, the POT extreme value analysis and estimation for stationary data series mainly involve the four following:(i)Step 1: determination of subsample Based on abundant monitoring data, subsample corresponding to the certain threshold and number *r* should be determined. To select the suitable threshold , the average excess function *e*(*u*) of the excess quantity *u* is computed. The threshold can be determined by focusing on the slope change characteristics of *e*(*u*). The turning point of *e*(*u*) can be used as the threshold.(ii)Step 2: maximum excess quantity calculation Then, the maximum excess quantity for each subsample is obtained through calculating deviations of series data and the selected threshold .(iii)Step 3: GPD fitting and parameter estimation Assuming that it is independent, each maximum excess quantity can be fitted by the GPD model with 2 parameters, which can be expressed as follows: where is the shape parameter and is the scale parameter. The variable can be denoted as follows: The GPD fitting is conducted with data to estimate the shape parameter and the scale parameter .(iv)Step 4: extreme value estimation With estimated parameters, the POT model described by the GPD is established to estimate maximum excess quantity. Therefore, the GPD-based extreme value of data series can be calculated as where *p* is return level, is the number of subsamples that exceed threshold , and *n* is the amount of data. *p* can be determineded considering the local code. Therefore, the extreme value of temperature-related and vehicle-induced cable force can be obtained, respectively.

#### 3. Case Study

The monitoring data in this paper were collected from an SHM system and a weighing system on a highway cable-stayed bridge crossing Yangtze River in Nanjing, China. Several typical cables were analyzed to select the critical position. Usually, cables designed with maximum force or the longest length were emphatically studied. The traffic data from weighing system was used as a probabilistic database for random simulation. In addition, all the sensors on the bridge are brand new in early years so that their measurements can be utilized for analysis.

##### 3.1. Descriptions of Bridge and Its Monitoring System

The third Nanjing Yangtze River Bridge was built in 2005, which connects Nanjing City and its Pukou District. It is a cable-stayed bridge with a main span of 648 m and a tower height of 215 m and has bidirectional 6 lanes. The main girder is supported by 168 stay cables, each consisting of 109 to 241 wires with a diameter of 7 mm. The appearance of the bridge is shown in Figure 6. A sophisticated monitoring system including an SHM system and weighing system is then devised to continually track the operational situation and structural health status. The advanced anchor load cell (ALC) sensors, which have higher precision compared to other cable force measurement approaches, are installed in the anchors of stay cables [42]. All the 168 stay cables are equipped with ALC sensors with an acquisition frequency of 10 Hz and a relative error of ±1%. Different from common acceleration sensors, the ALC received cable force measurements directly. Meanwhile, environmental temperature data can be collected by temperature sensors (TEM) around towers with R. M. YOUNG-41372 thermometers. The range of this sensor is [−50°C, 50°C] with the precision of ±0.3%. Temperature data are also collected with an acquisition frequency of 10 Hz.

The bridge weighing system is set at the bridge toll stations (both south and north banks). For each vehicle passing the bridge, the system will record its pass time, vehicle class, total weight, axle weight, axle distance, and driving lane. The recorded traffic data are used to establish random vehicle models combined with other parameters, such as speed and acceleration.

##### 3.2. Estimation of Temperature Effect

For a practical cable-stayed bridge, the longest cables near the mid-span are usually designed to bear the maximum force. Figure 7 shows the designed forces and length of the third Nanjing Yangtze River Bridge for all stay cables. It can be seen that cable NJ21 is the longest cable (with a length of 354.53 m) of this bridge with almost the maximum force (designed as 3708 kN). Different from cable NA21, which is also a long cable but only in the side span, cable NJ21 directly supports the girder in the main span, so more attention has been paid to cable NJ21 by the bridge management department [43]. Hence, cable NJ21 in the third Nanjing Yangtze River Bridge was selected for a case study to validate the effectiveness of the proposed method. The structure of cable NJ21 is shown in Figure 8; it constituted 241 wires, each with a diameter of 7 mm, so the area of total wires can be calculated as 9275 mm^{2}. The wires are galvanized and have an ultimate strength of 1670 MPa.

In Figure 9, the total of 530 daily pieces of cable force monitoring data of cable from January 2007 to October 2009 were used and then separated into 15 layers with db12 wavelet basis functions. Total cable force induced by environmental temperature effect and dead load was obtained, as shown in Figure 9(b). The abnormality in early 2008 was caused by a sudden snow disaster, which resulted in the significant increase of cable force monitoring data. The monitoring cable forces vary around the initial value of 3708 kN. To ensure the safety of cables and the total structure, a designed safety coefficient of more than 3 was adopted for this bridge, which accounted for the measurements being far below the limit material strength.

**(a)**

**(b)**

**(c)**

The obvious linear relationship can be observed between and environmental temperature, as shown in Figure 10. The monitoring temperature data in the corresponding period were also extracted to obtain the constant parameters *k* and *c* in equation (3). Finally, can be expressed as

Extreme value analysis was then conducted for temperature data with the GPD model. The threshold was selected as 32°C, while *r* was set to 3. Since the extreme high temperature must appear in summer, the monitoring temperature data of summer (June, July, and August) from 2006 to 2010 was utilized for extreme value analysis. The MATLAB platform was adopted to calculate the shape parameter = −0.221 and scale parameter = 3.076. Figure 11 shows the probability density function (PDF) results from the GPD fitting.

The two parameters of the GPD processing were then substituted into equation (8) to calculate the extreme temperature with a return period of 100 years, and the result was 45.095°C. According to equation (9), the temperature-dead load-induced extreme value cable force with a return period of 100 years was 3800.40 kN.

##### 3.3. Estimation of Vehicle Effect

The weighing system data of the third Nanjing Yangtze River Bridge from January 2006 to April 2009 were extracted as a sample to analyze the vehicle characteristics. According to different vehicle axle distribution, the data were divided into 12 types, as shown in Table 1. *P* denotes the proportion of each vehicle type within the total vehicle number. As for the vehicle axle, 2 denotes one axle and two wheels, 4 denotes one axle and four wheels, 8 denotes two axles and eight wheels, and 12 denotes three axles and twelve wheels. Thus, the parameters of random vehicle model can be determined.

Based on the associated information of vehicle types, the mixture normal distribution (MND), mixture lognormal distribution (MLD), and EM algorithm were applied to establish the statistical models of vehicle weight of different vehicle types, and the t-distribution was used to establish the statistical models of vehicle speed. The vehicle arrival time intervals were determined according to Poisson process theory. Figure 12 provides some examples of the PDF of vehicle weight (mixture lognormal distribution), vehicle speed (t-distribution), and arrival time interval (exponential distribution), respectively. According to the process of determining the parameters of random vehicle model, it is obvious that all the parameters are obtained through the analysis of vehicle data collected from the weighing system in real bridge. Then, the random vehicle model was developed with the Monte Carlo method, which can be used to simulate dynamic traffic flow for the third Nanjing Yangtze River Bridge.

**(a)**

**(b)**

**(c)**

In order to verify the reliability of the random vehicle model developed in this study, the simulated vehicle flow was acted on structural influence lines to calculate the time history of vehicle-induced responses, which was then compared with the monitoring cable forces. Similarly, the cable NJ21 was selected to be investigated. The cable force influence line for cable NJ21 and the layout of lanes are shown in Figure 13.

**(a)**

**(b)**

To evaluate the accuracy of simulation results, the actual average hourly traffic volumes of March 2009 were used as an input of the established random vehicle model. Meanwhile, cable force monitoring data for the same period were extracted and separated. A comparison for the calculated result and the monitoring data of cable NJ21 was conducted during the period from 11:00 to 12:00 on March 6, 2009, as shown in Figure 14. The time history of simulated vehicle-induced cable forces matched well with the monitoring data in normal operation condition. The specific values for comparison are shown in Table 2. In terms of maximum, minimum, and average values, the ratios of calculated result to monitoring data were 0.964, 0.881, and 1.036, respectively. The simulation results were proved to be in line with the actual situation. The slight deviation may be caused by the neglected wind effect, which will be considered in further research.

Moreover, the Kullback–Leibler (KL) divergence was introduced to quantitatively evaluate the simulation. The calculated KL divergence should be a positive value, and the more its value is close to zero, the more accurate the simulation method used is. Assuming that *P* and *Q* are two random variable series, the KL divergence can be expressed aswhere *p*(*x*) and *q*(*x*) are probability density functions of *P* and *Q*. The KL divergence between the calculated cable forces with the monitoring cable forces was calculated based on equation (10), which resulted to be 0.1219. The analysis of simulated and actual vehicle-induced cable forces shows that the developed random vehicle model was reliable for simulating response in a specific operational condition when the available records were limited.

The average daily traffic of the third Nanjing Yangtze River Bridge in the analysis period was about 27,000. To predict the extreme value of vehicle-induced responses, prospective traffic volume should be analyzed. In this respect, influencing factors including population, industry, and business development need to be considered. In this study, the traffic volume of another long-span bridge (Jiangyin Bridge) located at Yangtze River near the case bridge was investigated to analyze prospective traffic, for the operational situation of this bridge was close to the saturation level. The average daily traffic of Jiangyin Bridge was about 70,000 vehicles, while its ADTT (average daily truck traffic) was about 11,000. The traffic jam occurs in Jiangyin Bridge when the daily traffic exceeded 90,000. When the daily traffic volume reached the peak value of 130,000 vehicles on a holiday, a severe traffic jam happened. Mean while, the corresponding ADTT was only 5,000 in the same period. This traffic condition was almost the operational limitation of Jiangyin Bridge. Since ADTT is an important parameter that influences the vehicle-induced cable force, the passenger cars dominating traffic volume are not suitable to calculate response extreme value. As these two bridges are located in adjacent regions of the middle-lower Yangtze plains, the traffic data of Jiangyin Bridge was used as a sample to determine the prospective average daily traffic volume of the third Nanjing Yangtze River Bridge. The result was estimated to be 80,000 with the ADTT of 15,000.

Then, the received prospective traffic volume was taken as an input parameter of random vehicle model to calculate prospective vehicle-induced responses of cable NJ21. In order to consider various adverse conditions that may appear in service periods, vehicle-induced daily maximum cable force of cable NJ21 was simulated for a long period of five years; see Figure 15. With a determined threshold = 400 kN and *r* = 3, the GPD extreme value analysis was conducted via MATLAB platform. Thus, the shape parameter and the scale parameter can be obtained as = 0.0467 and = 24.1651, respectively. The PDF results from the GPD fitting are shown in Figure 16.

The values of the two parameters of the GPD were substituted into equation (8) to calculate the extreme value of vehicle-induced cable force with the return period of 100 years, and the result was obtained as 616.86 kN. According to the current Chinese design codes (JTG D60-2015), the design value of the vehicle-induced cable force was 626.10 kN in Highway-I vehicle load. It indicates that, in operational condition of prospective traffic volume, the extreme value of vehicle-induced cable force with the return period of 100 years was still lower than the design value.

##### 3.4. Estimation Results

To illustrate the predicted extreme value of cable forces in prospective operational conditions, Figure 17 plots the variation of the predicted extreme values of temperature, vehicle induced responses, and total responses with respect to the return period obtained using the proposed methods.

**(a)**

**(b)**

The extreme value of total cable force in prospective operational conditions with the return period of 100 years can be calculated by equation (2); namely, = 4417.26 kN. As for two ingredients of cable force, only a small percentage of responses were induced by vehicle (616.86 kN and 13.96%). The corresponding extreme stress of wires was 476.25 MPa, which was less than the limitation of 668 MPa, calculated from the Chinese design code. The estimation result shows that the longest cable NJ21 of the third Nanjing Yangtze River Bridge will be in a safe status without considering the material degradation. However, it should be mentioned that when more monitoring data becomes available in the future, it needs to be included to update the prediction in order to reduce the uncertainties of parameters. The method can also provide reliable estimates for the cable forces in specific return periods.

#### 4. Conclusions

This paper proposed a practical approach to analyze and predict the performance of stay cables of long-span cable-stayed bridges based on cable force monitoring data and traffic data. After raw data preprocessing, a random vehicle model was developed and applied to predict the prospective cable force with extreme value theory. To validate the reliability of the proposed method, measurements from the SHM system and the weighing system of the third Nanjing Yangtze River Bridge were involved for a case study. According to the results of the study, the following conclusions can be drawn:(1)The cable forces for a cable-stayed bridge are considered as the combination of different action effects. The wavelet transform method is suitable to separate the monitoring cable force induced by different effects.(2)The multiparameter vehicle modelling method was developed to generate Monte Carlo random traffic flow based on the collected bridge weighing system data and then calculated vehicle-induced response. The vehicles can be divided into specific several types according to the numbers of axles and wheels.(3)As a supplement to the normal distribution, the GPD-based extreme value method was utilized to predict extreme cable force in a specific return period. In the application, the estimated extreme cable force in prospective operational conditions with the return period of 100 years was 4239.04 kN.(4)The proposed analysis and estimation methods were applied to the third Nanjing Yangtze River Bridge. Based on the actual cable force monitoring data and the bridge weighing system data, the extreme cable force in prospective traffic operational condition for cable NJ21 was obtained. The method provided reliable estimates for the cable forces in specific return periods with limited data.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

The research was supported by the Natural Science Foundation of Jiangsu Province (no. BK20181278), Transportation Science Research Project in Jiangsu (no. 2019Z02), the Postgraduate Research & Practice Innovation Program of Jiangsu Province (no. KYCX19_0099), and the Academician Special Science Research Project of CCCC (no. YSZX-03-2020-01-B).