Inverse Evaluation of Monopile Pile–Soil Interaction Parameters Using Random Search

: To deal with the uncertainties in modeling offshore wind turbines, we propose a parameter inversion method for the pile–soil interaction model based on structural health monitoring results and the numerical model. The proposed parameter inversion method has a numerical model, an objective function selected using both the numerical and identiﬁed results


Introduction
Currently, most offshore wind farms have structural health monitoring systems that record extensive in situ operation data at a predefined time resolution.However, during monitoring, these operation data are often mixed with signal disturbances in the ocean environment.Furthermore, a large portion of monitoring data is usually archived as mean, median, or maximum values over a certain period to reduce data storage.Therefore, it is very challenging to analyze the monitoring data, neither identifying nor predicting possible faults in offshore wind turbines.
Nonetheless, plenty of exciting work has been conducted worldwide on structural health monitoring data analysis.Jahani [1] discussed the crucial concerns related to the structural dynamics of offshore wind turbines.The resonance frequencies and damping values of dominant mode shapes for a parked wind turbine are tracked using an automated state-of-the-art operational modal analysis technique [2].Modal parameters were used as a baseline for long-term observations [3].Estimates of natural frequencies, damping ratios, and mode shapes of offshore wind turbines were obtained by the poly-reference least squares complex frequency-domain estimator (PolyMAX) and the covariance-driven stochastic subspace identification (SSI-COV) method [4].The stochastic subspace identification method has been used to reveal the modal parameters of offshore wind turbines experiencing earthquakes [5].Substructure damping [6], aerodynamic damping [7], and damping of wind turbine towers [8] have all discussed, and Eigensystem Realization Algorithm (ERA), SSI-COV, the Enhanced Frequency Domain Decomposition (EFDD) [8], and time-frequency analysis [7] have all been used.Estimating offshore wind turbine damping using state-of-the-art operational modal analysis (OMA) techniques was summarized in [9].Besides, the relationship between modal parameters and environmental/operational factors was illustrated [10], and vibration-based damage detection was also discussed [11].
Meanwhile, analyzing and utilizing offshore wind turbine monitoring data is challenging, as the marine environment is full of interference signals.To reduce the adverse effects of interference, Dong [12] proposed a modified stochastic subspace identification (SSI) method considering known and time-invariant harmonic interference.A time-frequency method based on single-mode function (SMF) decomposition was presented to reduce mode mixing and data contamination effects in [13].Data containing high-energy components were identified by adopting iterative noise elimination by Liu [14], which was confirmed to be better than the stochastic subspace identification method.Noise cleansing and missing data imputation were found to increase the accuracy of fatigue assessment when mutually applied [15].Multi-dimensional information fusion was used based on Dempster-Shafer (D-S) evidence theory for the fault diagnosis of wind turbines in [16].In this paper, the monitoring data was confirmed to be mixed with trend terms in polynomial form.Therefore, the noise can be removed by a preprocessing algorithm for removing trend items; the details are not provided here, as they are common practice.
Offshore wind turbines are seldom equipped with environmental monitoring or load monitoring equipment.Therefore, the response of offshore wind turbines is always an output-only identification problem, with these related excitations being absent [17][18][19].In the meantime, as offshore wind turbines are always modeled using a linear beam and the soil-structure interaction is considered using American Petroleum Institute (API) based reaction springs attached to the monopile [20], the uncertainty in soil modeling is always overlooked [21][22][23][24].Andersen [21] studied the first natural frequency of a monopile foundation interacting with clay-type soil having a spatial variation of properties.They found that if the undrained shear strength of the soil has a log-normal distribution, the probability density functions (PDFs) for the pilehead stiffness components are also lognormal as well.However, this is different for the natural frequency.Damgaard [25] found that soil-structure interaction is critical in offshore wind turbine design.Therefore, it is possible to obtain good dynamic modeling of offshore wind turbines only after tuning the first natural frequency and damping.Arany [26] found that the cross-coupling spring term needs to be included in the offshore wind foundation analysis.Zania [27] found that soil properties mainly affected the modified soil-structure interaction eigenfrequency and damping of slender engineering structures, such as offshore wind turbines.The dynamic soil-structure interaction may drive even a conservative design to restrictive frequency ranges.Many more results show that the dynamic response of a wind turbine is highly sensitive to the properties of the soil where it has been built [22][23][24].To deal with the uncertainty in the modeling and excitation of offshore wind turbines, we propose an inverse evaluation method that can fit the numerical model to the in situ monitoring data using the parameters identified here, which can be considered the in situ parameters.
The motivation of this paper is to propose a method to identify the in situ parameters of offshore wind turbines and analyze the influence of uncertainty on the initial design parameters.In Section 2, the methods related to the inverse evaluation are provided.The results obtained from the inverse evaluation method are discussed in Section 3. Subsequently, Section 4 ends with conclusions.

Methods for Inverse Evaluation
In this section, we first briefly introduce the sensor deployment and the monitoring data.Then, the data are analyzed using different operational modal analysis methods.The results of the operational modal analysis are analyzed to check whether there are invariant (stable) modal frequencies during a relatively long time.After the invariant modal frequencies are obtained, they can be used to construct the inverse optimization objective function.Subsequently, the pile-soil interaction parameters of the monopile can be identified using the optimization objective function.

Monitoring Data and Sensors
The monitoring data was obtained from a monopile in an offshore wind farm in the eastern China Sea, where the water depth is about 10 m.The monopile foundation was installed with a 4.5 MW offshore wind turbine.The total length of the pile foundation is 54 m, the axial compression capacity is 27,887 kN, and the allowable frequency range for the wind turbine design is 0.24 to 0.32 Hz.The configuration of acceleration sensors on the wind turbine is shown in Figure 1.These accelerations are used in the following analysis.The monitoring data consisted of accelerations from four elevations on the wind turbine tower within 24 h.The data were stored in 24 files (sampling time of about 60 min each), with a sampling frequency F s = 50 Hz.During the monitoring period, the rotation speed of the wind turbine was 3.20 to 10.57 revolutions per minute (rpm).The monitoring points are listed in Table 1.

Location
Sensor ID Elevation Direction 1 12 m Y 1 X-direction: prevailing wind direction; Y-direction: vertical to prevailing wind direction.
From this wind farm, we obtained the prevailing wind direction according to the wind data observed over two years, using the definition that the direction with the longest accumulated occurrence time in a year or the direction with the longest petals in the annual wind rose-petal chart developed by the local meteorological department.Meanwhile, according to the official documentation of the offshore wind farm we investigated, no pile slope was recorded during pile driving, and the wind farm is not located in seismically active areas.Therefore, pile slope and soil liquefaction during pile driving is not considered.

Analysis of In Situ Modal Frequency Characteristics
This section uses operational modal analysis methods to analyze the monitoring data.However, the modal characteristics should be confirmed as stable before conducting the inverse evaluation of the monopile.Subsequently, inverse evaluation can be performed as follows: based on the stable modal characteristics, an objective function for parameter optimization is constructed, then an optimization algorithm is used for the inverse analysis of the in situ parameters.

Modal Identification Based on Subspace Identification
The stochastic subspace identification (SSI) method is a dynamic identification method that can deal with dynamic systems with uncertain inputs.This method constructs the Hankel matrix using the existing data, then forms the Toeplitz matrix to solve the coefficient matrices (SSI-COV, covariance-driven stochastic subspace identification) or forms the output projection matrix to solve the coefficient matrices (SSI-DATA, data-driven stochastic subspace identification).The SSI-COV algorithm is shown in the following steps.Assuming that there are l measuring points (possibly different sensor types), the following Hankel matrix is constructed by selecting r points as reference data [28]: where y ref 0 ∼ y ref i+s−2 and y i ∼ y 2i+s−2 are column vectors composed of r reference data and l measuring data, respectively, and s is the order of the model.Assuming that the dynamic system is ergodic, the block Toeplitz matrix T ref 1|i is composed of the covariance between the measurement data and the reference data, which satisfies The block Toeplitz matrix T ref 1|i can be decomposed into [28] where is the covariance of the state data and reference output data.
After the block Toeplitz matrix T ref 1|i is obtained from Equations ( 1) and (2), the observability matrix and controllability matrix, i.e., O i and C i in Equation (3), can be solved by singular value decomposition (SVD) of T ref 1|i .Subsequently, the state matrices (A and C) of the dynamic system can be obtained.The modal characteristics are the eigenvalues (eigenvectors) of the state matrices [28,29]: where Ψ A and λ A are the eigenvectors and eigenvalues, respectively, of the state space matrix and Ψ, λ, f , and ξ are the system's mode shape, complex eigenvalue, natural frequency, and damping ratio.Another modal identification method originating from control theory is the subspace identification method, which includes several variants, including N4SID (numerical algorithms for subspace state space system identification), MOESP (multivariable output error state space), and CVA (canonical variable analysis) [30].As the three subspace algorithms are special cases, with each of the algorithms a specific choice of the weighting matrices [30], only N4SID is selected as a concept demonstration example, shown below.
The N4SID algorithm begins with the discrete-form equation of a dynamic system [31,32] where The model order s should be larger than the system order.
s ≥ max{Rand(A), Rand(B), Rand(C)}. Let According to the above definition, the extended system equation is Assuming that the system states are ergodic and the system input is independent of the state and measurement noise, the matrices Γ s , H s , and G s can then be obtained by the convex optimization least squares method using Equation (11).Subsequently, the state matrices A, B, C, and D of the dynamic system in Equation ( 5) can be obtained using the block matrix projection relationship in Equations ( 6)- (10).The modal characteristics represent the eigenvalues (eigenvectors) of the state matrices.
To perform an inverse evaluation based on dynamic characteristics, it is necessary to check whether the modal characteristics are stable.The objective function of the inverse optimization can only be constructed based on these stable modal characteristics.Subsequently, inverse evaluation can be performed using an optimization algorithm according to the pre-mentioned inverse optimization function.

Methods for Inverse Evaluation
When offshore wind turbines are installed, their in situ working state can only be tracked by structural health monitoring systems.However, structural monitoring systems are rarely equipped with environmental monitoring devices, which are always extremely expensive.Therefore, it is not easy to numerically evaluate the in situ response of offshore wind turbines as the input (excitation load and environment data) and output (response) cannot be obtained simultaneously.To evaluate the in situ state of offshore wind turbines, we propose an inverse evaluation method based on modal characteristics.The implementation procedures are provided below.

Physical Model for Inverse Evaluation
A numerical model of the monopile should be established to perform the inverse evaluation procedure as a prerequisite of the optimization step.
For offshore wind turbines, as the main structural components such as the turbine foundation, tower, or blades are always long and slender, these structural components are always considered as beams or bars.In contrast, flanges and other local masses are modeled as concentrated mass points.
Monopiles with different marine soils are permanently installed far offshore, where wind, waves, and currents are present.According to the regular wave theory selection diagram in IEC 61400-3-1 [33] (Annex B), H/gT 2 is 0.01 and d/gT 2 is 0.019; therefore, the seventh-order stream function theory is used to model waves.The current is considered using the standard velocity profile in IEC 61400-3-1.However, the wind load is not calculated here; only the results provided by the wind turbine supplier are used.
As the physical properties of marine soil are highly dispersed [34][35][36][37], the soilstructure interaction of offshore wind turbines is one of the prominent sources of uncertainty in their dynamic response [22][23][24].These uncertainties could lead to disturbances in the condition evaluation of offshore wind turbines that cannot be ignored.Currently, soil reaction curves for piles under axial compression or lateral loads, i.e., t − z, Q − z and p − y in American Petroleum Institute (API) RP 2GEO [37] (or their recently-developed modified variants, such as the PISA model [38]) and classical soil constitutive models included in commercial finite element software are commonly used.As the soil reaction curves are easy to implement and have lower computational costs, they are used to show the proposed methodology.
To reduce the modeling uncertainty of offshore wind turbines, the inverse evaluation method is based on in situ optimization of the soil-structure interaction parameters, which are one of the primary sources of uncertainties in offshore wind turbine modeling and analysis.Recently, it has been argued that the soil reaction curves can be rather conservative with regard to the soil response [38].However, the proposed method may be enlightening, even though the parameters may suffer from a deficiency in model assumptions.
A schematic of the physical model proposed for offshore wind turbines is shown in Figure 2. In the numerical model, the structural geometry, materials, site conditions, and loads of the offshore wind turbines are first provided in several input files.A series of soil springs are formed to simulate the soil reaction response, and a required amount of wind, wave, current, and other external conditions are combined to obtain the environmental boundary conditions.Then, a collection of solution programs operate on these inputs, producing a solution containing joint displacements and internal element forces.Subsequently, another category of postprocessing programs evaluates the structure response using the solution shown above.

Objective Function for Inverse Evaluation
This paper combines the in situ modal identification results using structural monitoring data and numerical results obtained to form the objective function.In order to accomplish this, it should first be confirmed that the modal identification results are stable over a relatively long time, such as one day (twenty-four hours).Otherwise, the objective function based on the obtained results does not make sense, as it varies without parameter variations.
The objective function is selected as follows: where • denotes the 2-norm, f num and f are the first-order modal frequency obtained by the numerical model and that obtained by the monitored data, respectively, and f 1+α 1−α represents the first-order modal results that match f num , defined as the minimum f value satisfying where α is the variation ratio of the modal frequencies, which measures the span where the modal frequency varies, α = 0.05 ∼ 0.1 (DNV-ST-0126 [39]).
After the objective function is defined, to find the optimal value that best fits the in situ data, the inverse evaluation procedure should traverse as many values in the parameter space as possible, if not all.

Parameter Space for Inverse Evaluation
According to the geotechnical report of the wind farm, the investigated offshore wind turbine is installed on a site with eight strata of marine soil, with 42 total parameters comprising nine different types.The soil-structure interaction parameters are provided in Table 2.The statistical distribution of the soil parameters is highly spatially discrete.Therefore, the variation of the soil parameters can be assumed to follow a normal or log-normal distribution [34,35].It is assumed that the parameter X of the soil-structure interaction model in different soil layers obeys the normal distribution, with a mean value µ provided by the geological survey (the value X 0 used for foundation design), µ = X 0 , and coefficient of variation κ between 0 and 0.5, κ = 0 ∼ 0.5, that is, where σ X = κX 0 .The coefficient of variation is not the same for all parameters in reality.However, without loss of generality, an identical and constant coefficient of variation κ is used here to conduct a proof of concept for the inverse evaluation method.It is assumed that the parameters in different strata are independent and that only related parameters are considered in a single soil stratum, depending on the soil classification.In the sensitivity analysis, all parameters are kept constant in every other stratum.Meanwhile, relations between different parameters in a given stratum are employed to reduce the number of independent parameters per stratum.
After considering all of the above constraints, a screening process is first conducted to select parameters with higher single-factor sensitivity through sensitivity analysis to reduce the number of parameters for inverse optimization.Then, the chosen parameters are assumed to follow the distribution defined in Equation ( 14), and three to five representative values are randomly selected as the inputs for the optimization algorithm at each iteration.

Random Search Algorithm for Inverse Evaluation
Global optimization methods include grid search, random search, Bayesian optimization, gradient-based optimization, etc.As the efficiency of grid search is limited by the number of parameters and the values of parameters, Bayesian optimization relies on prior information to establish the probability model of the search parameters and objective function, and gradient-based optimization can quickly run into local optimization situations, a random optimization algorithm is used here for global optimization of the pile-soil interaction parameters.Moreover, the random optimization algorithm is a good choice, as there are many soil parameters for which the distribution is roughly unknown, and this algorithm can obtain globally optimized results in fewer iterations than grid search [40].
Considering the upper and lower bounds of the parameters, the random search algorithm takes a random combination of parameter values from the interval and carries out parameter updating through an iterative process.In addition, by designing a suitable parameter updating algorithm, random search can overcome the limitations caused by assumptions about the initial parameter distribution.
The random search algorithm can be implemented using the following steps: 1.
Define the objective function G(X), penalty function P( f , α), and coefficient α (if used).
The convergence criteria and the maximum number of iterations should be set as well.

2.
Determine the upper and lower bounds of the parameters [X] and estimate the distribution X ∼ f X (x).

3.
Set an initial parameter value X 0 and solve the initial objective function value G(X 0 ).

4.
Generate m random parameter values X1 , X2 , ... Xm according to parameter distribution f X (x) and solve G( Xi ), If there exists min G( Xi ) ≤ G(X 0 ), set X 0 = Xi .Repeat Steps 3 and 4 until the iteration reaches the convergence criteria or the maximum number of iterations.The final Xi is the optimized parameter value under the convergence criteria or the maximum number of iterations.
The flow chart of the random search algorithm is shown in Figure 3.

SHM System
Modal results

Figure 3. Flow chart of the random search algorithm
After the in situ modal frequencies of monopiles are identified as stable, an objective function can be selected using both the numerical and modal identification results.
Subsequently, inverse optimization can be performed using a random search algorithm in the assumed parameter space.The parameter results in the minimum optimization objective function are identified as the in situ parameter of the monopile, as this results in the numerical model best approximating the in situ structure.

Results and Discussion
To analyze the in situ state of offshore wind turbines, inverse evaluation is performed here to obtain their in situ parameters based on a random search algorithm using numerical modal results and those identified from monitored data.The parameters used in the numerical model are provided in Table 3.A scour depth of 6 m was used in the numerical model according to the site conditions and experience from offshore wind farms nearby.The soil parameters from the geological survey are listed in Table 4 (not all).This manuscript considers a total of eight strata, including forty-two parameters of nine different types.It is assumed that the parameters in different strata are independent, and only three to seven parameters are considered in a single soil stratum, depending on the soil classification.Therefore, in the sensitivity analysis, all parameters are kept constant in every other stratum.Meanwhile, relations between different parameters in a stratum are employed to reduce the number of independent parameters per stratum.After considering all the constraints above, a screening process is conducted to select higher-sensitivity parameters, after which the chosen parameters are assumed to follow the normal distribution, and three to five representative values are randomly chosen as the inputs for the optimization algorithm at each iteration.

Operational Modal Analysis
The operational modal results are first analyzed to check whether stable in situ modal frequencies exist, as this is a prerequisite for the inverse optimization step.

Verification of Stable Modes
Both SSI-COV, SSI-DATA, and N4SID are used in operational modal identification for mutual verification.In the modal identification procedures, the model order s is assumed to be 40.The modal results confirmed that the model order s is suitable (results not included here).The modal convergence criteria are set to an error in modal frequency and damping ratio no greater than one and 5%, respectively, and a mode shape consistency no less than 99%.The results obtained from different modal identification methods are shown in Figure 4.As shown, the results from N4SID are closely consistent with either SSI-COV or SSI-DATA.However, there are more modal results in SSI-COV than in SSI-DATA, which may result from differences in the solution method.
As shown in Figure 4, the Complex Mode Indicator Function (CMIF) is used for method verification, which is defined as [41] CMIF k (jω) = σ 2 k (jω) (15) where σ k (jω) is the k-th singular value of the frequency response function (FRF) matrix H( f ) at circular frequency ω = 2π f , i.e., where U( f ) and V( f ) are the complex unitary matrices and S( f ) is a rectangular diagonal matrix with non-negative real numbers on the diagonal.The diagonal entries of S( f ) are uniquely determined by H( f ) and are known as the singular values.The maximum of the CMIF corresponds to the characteristic frequency of the system.As shown in Figure 4, the CMIF maxima are the same in SSI-COV and SSI-DATA results, as the FRF results obtained from the structural monitoring data are identical for the same monopile.Meanwhile, the maxima of the CMIF results coincide with the N4SID results, which is evidence that the modal identification results are convincing.As the three modal identification methods provide satisfactory results here, the N4SID algorithm is used in the rest of the paper, considering its benefits in computation cost and postprocessing of the results.
To eliminate fake modes resulting from mathematical poles and to determine the mode order for modal matching in the inverse evaluation, the modal assurance criterion (MAC), defined as the cross-operation results between all identified mode shape vectors, is used here.The MAC is defined as [42] MAC where Ψ i and Ψ j are the mode shape, with order i and j.The first ten mode results from N4SID are shown in Table 5.The corresponding MAC results are shown in Figure 5. Multiple large non-diagonal values (greater than 0.5) appeared in the MAC results in the modal identification results, indicating the presence of several highly correlated modes.Furthermore, it can be seen from Table 5 that unstable modes with large damping ratios (greater than 10%) appear in the modal results, which are physically unrealizable and can be considered as results from highly correlated mode shapes.Figure 5b shows the results with the fake modes removed.The corresponding MAC results are satisfactory, indicating that the correlation between each mode meets the requirements for stable modal results.
The above discussion confirms that stable modal results can be identified and separated satisfactorily.However, as we need to use the modal results for inverse evaluation, it is necessary to check whether the modal frequencies are invariant or are only engaged with slow or minor changes, i.e., approximately constant, during a relatively long time (for example, one day).

Checking the Stationary Modal Frequencies
According to the modal identification and verification procedures proposed above, the results at different times in a day are shown in Figure 6 (only the results for the first ten orders are included).Time IDs are custom labels used to identify different times in a day, with no practical significance otherwise.In addition, Time IDs 1 to 4 correspond to the initial startup of the studied wind turbine, which is unstable in its working state, while Time IDs 5 to 24 correspond to power production.To eliminate the uncertainty in the working state of the wind turbine, the following discussion focuses on the results from Time IDs 5 to 24.
As shown in Figure 6, the first-order modal frequencies of the offshore wind turbine during power production are approximately constant, varying between 0.247 Hz and 0.291 Hz, with a median characteristic frequency of 0.274 Hz repeatedly occurring at a different time in one day.However, the corresponding frequencies are slightly changed at other times, with a median fluctuation range of -10% to 5.9% compared to the characteristic frequency.Meanwhile, the damping ratios show an irregular and drastic change between 0 and 10%.The above results quantitatively confirm that the first-order modal frequencies of offshore wind turbines are not a constant value during power production.On the other hand, it confirms that the first-order modal frequencies can be considered approximately constant here, and that a first-order modal frequency f 0 = 0.274 Hz can be used.
A further comment should be made here.Based on the above results, it is difficult to determine whether the first-order modal frequencies continue to exist in whole or in part for longer periods, such as 48 h or longer.However, it is sufficient for our purposes here to demonstrate the proposed inverse evaluation method.Further work is needed to reveal the modal characteristics of offshore wind turbines during their working state.

Selection of Sensitive Parameters
According to the geotechnical report of the wind farm, the investigated offshore wind turbine is installed on a site with eight marine soil strata.As shown in Table 4, each soil stratum is classified as silt, clay, or sand according to the grain size, with five, three, or seven parameters in each stratum, respectively.The eight soil strata are composed of forty-two parameters categorized into nine types.For the convenience of parameter selection, it is temporarily assumed that the coefficient of variation used to screen the sensitive parameters of the soil layer is κ = 0.5.
According to the single-factor sensitivity analysis results, the soil parameters can be classified into four categories, i.e., highly sensitive, sensitive, marginally related, and uncorrelated.Selected sensitivity analysis results for the first two categories are shown in Figure 7.The rest of the sensitivity results are omitted, as the idea behind the selection criterion is easy to understand.Meanwhile, as there are deficiencies in the geological measurement mechanisms of certain parameters that may lead to deviations in the geotechnical results, a total of thirty parameters for random search and inversion are considered here based on the sensitivity analysis.

Convergence Verification
Iterative optimization can be carried out by selecting all parameters (a total of fortytwo) or selected sensitive parameters (a total of thirty).Considering the computing capacity of our workstation (Intel i9 12950HX, NVIDIA RTX A3000, 128 Gigabyte RAM, 2 Terabyte Solid State Disk), the maximum number of iterations was set to n = 100 and the convergence criterion was set to a residual of the objective function not greater than 10 −4 .It was found that the computation time for one parameter case required approximately two to six hours.
During parameter iteration, we found that certain combinations of independent parameters led to cases in which the pile solution (either axial or lateral) failed to converge or in which the combined reduced structural stiffness matrix and pilehead stiffness matrix were non-positive definite.These parameter combinations can lead to physically meaningless or erroneous results.Therefore, these parameter cases were removed, with only the valid iterations plotted versus the objective function's residual.
The convergence curve of a typical calculation case is shown in Figure 8.As shown, the results confirm that this method can converge to a specific soil parameter combination, which is considered the actual soil parameter during the working state of the wind turbine.

Sensitivity Analysis of Initial Parameters
Optimization methods can be confusing when they are susceptible to initial values.To check the ability of the algorithm to always converge even under arbitrary parameter cases, five different parameter cases are proposed to show the convergence characteristics.The convergence curves are provided in Figure 9.

As shown in Figure all initial parameter cases (large objective function values corresponding to iteration numbers equaling zero) can converge to desirable objective function values (approximately zero)
. Meanwhile, although the convergence rate varies with different initial parameter cases, the ideal objective function value (less than 2.0‰) can be reached in four iterations at the least.All initial parameter cases shown can reach the convergence criterion in roughly ten iterations.Therefore, whatever the initial parameter values, the algorithm can always converge after several iterations corresponding to the in situ parameter values.
,WHUDWLRQQXPEHU 2EMHFWLYHIXQFWLRQ &DVH &DVH &DVH &DVH &DVH The initial parameters can be considered as parameters provided for the initial design of offshore wind turbines, which are always obtained by the geotechnical report of the wind farm.As a result of poor geotechnical investigation work, these parameters are sometimes found to deviate considerably from the actual values.
The results in this subsection confirm that whatever the provided parameters are, the proposed algorithm can always converge to a parameter combination after some number of iterations using the proposed inverse evaluation algorithm.The converged results correspond to the numerical results fitting the monitored results best, implying that the obtained parameters can be considered in situ.
Soil is not an ideally continuous medium, as its physical properties are discrete and spatially related.Therefore, the results of a geotechnical survey are reliable only at a certain probability level.However, the results of the geotechnical survey, i.e., the initial value of the inverse evaluation procedure, are vital, as they represent a significant input in foundation design for offshore wind turbines.Thus, it is necessary to evaluate whether all possible initial parameters obtained from a geotechnical survey can lead to identical or similar results.

Physical Characteristics of Inverse Evaluation Results
According to the inverse evaluation method, different initial parameter cases may converge to slightly different optimal parameters, i.e., multiple parameter combinations that all fit the monitored results well.In reality, only one combination of parameters should exist for the in situ parameters.Therefore, it is necessary to compare the results from different initial parameter cases to determine whether the structural response is numerically similar.
Based on the in situ monitoring data of an offshore wind farm in Jiangsu Province, China, five parameter combinations were obtained using the proposed inverse evaluation method regarding different initial parameter cases.Using these different parameter results from the inverse evaluation method, several design targets regarding monopiles, such as axial capacities and displacements, were analyzed, and the resulting deviations between the different results were compared.These results can be considered as the influence of initial geological parameters, which can shed light on the impact of uncertainty in geological surveys.

Deviation in Axial Capacities
As the axial bearings of monopile foundations are controlled mainly by compression, their tension capacity is not discussed here.Compression capacities solved by the physical model using different inverse evaluation methods are shown in Table 6, where Max.Error corresponds to the maximum relative error between compression capacities.The safety factor and UC values are shown in Table 6, and are respectively defined as follows: where R d and S d are the maximum axial force and design capacity.According to NB/T 10105-2018 (Code for Design of Wind Turbine Foundations for Offshore Wind Power Projects) [43] issued by the National Energy Administration of China, the ultimate axial capacity of the structure should satisfy α ≥ γ 0 (21) where the structural importance coefficient γ 0 is usually taken as 1.1.
Because the maximum axial force on the structure is far less than the design capacity, the safety factor is much larger than the importance coefficient γ 0 ; therefore, the difference in the design capacity resulting from different initial parameter cases does not cause any safety risk.However, an excessive safety factor may lead to redundancy and lack of economy in materials.
As shown in Table 6, the axial capacities are different under different parameter cases, with a maximum difference of about 28%.According to the design report of the offshore wind farm, the compression capacity is 27,887 kN; as the maximum compression force is 8356 kN, the compression capacities in the parameter cases all meet the design requirements.However, the design value of the compression capacity parameter in case 5 is less than the original design capacity.
These results suggest that the design axial capacities of offshore wind turbines are sensitive to geological parameters; thus, a sufficient allowance for bearing capacity should be reserved in the design.

Deviation in Pile Displacements
To ensure the normal operation and power generation of offshore wind turbines, the foundation design should meet the requirements for the safe operation of wind turbines.Displacement at several positions on the turbine foundation solved by the physical model using different inverse evaluation methods is shown in Table 7, where Max.Error corresponds to the maximum relative error between the results obtained by different parameter cases.Figure 10 shows the positions for which displacements are provided.In Table 7, there are differences in the displacement and rotation angles of the tower bottom, pile head, and pile bottom obtained from the inverse evaluation method under different initial parameter cases.The maximum difference in displacement value at the tower bottom under different cases can reach about 3.42 cm, while the maximum difference in rotation angle is about 0.28 degrees.Moreover, the maximum difference in displacement value at the pile head is approximately 1.80 cm and the maximum difference in rotation angle is about 0.28 degrees.The maximum results at the pile bottom are about 0.68 cm and 0.01 degrees, respectively.Although the absolute value differences at the pile bottom (C) are insignificant, these differences indicate the correlation between the pile bottom response and the soil parameters.The accumulated difference caused by the soil parameters results in a relative difference of over 400% under different initial parameter cases.According to [43], issued by the National Energy Administration of China, the displacement and rotation at the tower bottom are not specifically controlled, nor are the displacement at the pile head nor the rotation angle at the pile bottom.In accordance with the deformation control requirements for monopile foundations in the NB/T 10105-2018, the displacement at the pile head should be less than 10.8 cm, the rotation angle at the pile head should be less than 0.00436 radians (0.25 degrees), and the displacement at the pile bottom should be less than 10 mm.
According to Table 7, all of the results satisfy the requirements in [43] except for the X displacements at the pile head (position B in Figure 10).The displacements are further compared with the design results in Table 8.As shown, the maximum displacements and rotation angles are all larger in the inverse evaluation results than those in the design results.These deviations are possible because of soil stiffness distribution deviation resulting from the difference in soil parameters, as the design results are based on the parameters obtained from the geological survey, while the inverse evaluation results are obtained from parameters identified based on monitoring data.
The displacement and rotation angle difference in the two results show that it is possible to underestimate the pile displacement and rotation significantly due to uncertainties in the geological survey.

Deviation vs. Design Results
To evaluate the current assumptions used in the inverse evaluation procedure, the obtained optimal parameters are introduced to the physical model proposed in Section 2.3.1.The discrepancies in the response of pile foundations between the inverse evaluation results and those of the original design are discussed in this section.
According to the assumptions used during the inverse evaluation, these optimal parameters may correspond to the in situ parameters of offshore wind turbines.Therefore, the discrepancies in the inverse evaluation results and those of the original design can be considered to be issues raised by uncertainty in the design phase of offshore wind turbines.Meanwhile, as these discrepancies resulting from input uncertainty may bring about potentially significant risks, they should be carefully investigated by synthesizing other information on the wind farm.

Discrepancy in Lateral Deflections
The lateral deflections obtained from the inverse evaluation and those of the original design are shown in Figure 11.As shown in Figure 11, the lateral deflections of the pile foundation at different elevations are roughly different from the original design results; however, although the magnitude of difference may seem significant, the in situ status of the wind turbine is satisfactory, as the lateral deflections are small overall and mostly smaller than the original design, meeting the design requirements of the foundation.

Discrepancy in Axial Deflections
The axial deflections obtained from the inverse evaluation and those of the original design are shown in Figure 12.As shown in Figure 12, the axial deflections of pile foundations at different elevations are significantly different from the original design results.Though the trend is roughly similar, the magnitudes of maximum deflection are 0.131 cm and 0.998 cm, respectively, in the original design and in situ results.This discrepancy in axial deflection results from the difference in soil parameters.The variation in soil parameters can lead to different axial capacities.In this case, the compression capacity obtained from the inverse evaluation is 24,016.8kN, which is far below the design result (27,887 kN).As a result, the axial deflection is increased compared to the original design.According to [43], the in situ status of the wind turbine is compliant with the standard, as the total deflection at the pile head (pile elevation equal to 0) is less than 10.8 cm.
However, as significant discrepancies in axial deflections are revealed, attention should be paid to the in situ axial capacity and deflection of the offshore wind turbines in this wind farm.The difference indicates that the uncertainty in soil parameters can lead to significantly more considerable in situ axial deflections than the original design results.Therefore, although the current axial deflections are acceptable, the discrepancy between the two results, i.e., the in situ axial deflection and the original design result, implies that attention is required concerning the sensitivity of axial deflection to the soil parameters to avoid potential risks.

Discrepancy in Pile Rotation
The rotations obtained from the inverse evaluation and original design are shown in Figure 13.
As shown in Figure 13, the pile rotations at different elevations are different, though smaller than the original design results, with a maximum value of 0.00153 radians vs. 0.00251 radians.These results are all acceptable according to [43], as all the rotation angles at the pile head (pile elevation equals 0) are less than 0.00436 radians (0.25 degrees).However, it is notable that, with a 39.04% smaller in-situ pile rotation angle, the difference between the two results is quite significant.This means the original design can be optimized for the pile rotation should it be the design constraint.

Limitations of the Proposed Inverse Evaluation Method
Based on the measured results and the numerical model, we proposed a parameter inversion method for the pile-soil interaction model using soil reaction curves.This method is promising and straightforward to implement; however, there are several flaws that should be noted.
First, the obtained results are related to the selected numerical model.Therefore, the proposed method can only correct the relevant parameters of the numerical model, and cannot overcome the shortcomings of the numerical model itself.A more realistic numerical model would make for more realistic results.
Second, the parameter space is based on certain assumptions.However, the actual distribution of the model parameters is uncertain.A more realistic parameter distribution would mean a faster search for actual results, if any.
Finally, the calculation efficiency needs to be further improved.Surrogate models, such as the response surface models, could be considered in subsequent research to accelerate calculation under the premise of maintained reliability.

Conclusions
Using monitoring data from an offshore wind farm in the East China Sea and numerical models, we proposed a parameter inversion method for a pile-soil interaction model based on the measured results and the numerical model.The proposed parameter inversion method involves several steps.First, it is necessary to construct a numerical model.Then, after the in situ modal frequencies have been identified as stable, an objective function can be selected using both the numerical and identified results.Subsequently, inverse optimization is performed using a random search algorithm in the assumed parameter space.The parameter results in the minimum optimization objective function are identified as the in situ parameters of the monopile, which result in the numerical model that best approximates the in situ structure.Our convergence verification results confirmed that this method could converge to a specific soil parameter combination.Furthermore, our results showed that whatever the initial parameter values, the proposed method can always converge after some number of iterations corresponding to the in situ parameter values.However, it was found that different initial parameter cases may converge to slightly different optimal parameters, i.e., deviations in axial capacities or pile displacements, implying that the pile results are sensitive to geological parameters.Therefore, it is necessary to make a sufficient allowance for uncertainties in geological surveys.The results obtained by the proposed method were compared with the original design results to check for design redundancy or risks.Although the proposed method has several flaws, it can shed light on the influence of parameter uncertainties on offshore wind turbines.

Figure 1 .
Figure 1.Acceleration sensors on the tower: (a) front view and (b) top view of section A-A.

Figure 2 .
Figure 2. Physical model of offshore wind turbine.

Figure 4 .
Figure 4. Results from different methods: (a) SSI-COV and (b) SSI-DATA.The vertical dashed line shows N4SID results, while • and × respectively correspond to stable or unstable modal results (poles).

Figure 5 .
Figure 5. MAC results of first ten identified modal results: (a) original results and (b) with fake modes removed.

Figure 6 .
Figure 6.Modal frequencies at different times (Time IDs) in a day: (a) modal frequencies and (b) damping ratios.

Figure 7 .
Figure 7. Sensitivity analysis results of selected soil parameters (for illustration only): (a) internal friction angle and (b) effective weight.

Figure 8 .
Figure 8. Convergence curve of a typical calculation case.

Figure 9 .
Figure 9. Convergence curves of different initial parameter cases.Cases 1 and 2 correspond to different initial parameter combinations, while σ corresponds to the standard deviation of the parameters.

Figure 10 .
Figure 10.Schematic diagram of positions where displacements on offshore wind turbines were obtained.

Figure 11 .
Figure 11.Discrepancy in lateral deflections from different pile evaluations.

Figure 12 .
Figure 12.Discrepancy in axial deflections from different pile evaluations.

Figure 13 .
Figure 13.Discrepancy in pile rotation from different pile evaluations.

A, B, C, and D are
the coefficient matrices of the state space, u k , y k , and x k are the input vectors, output vectors, and state vectors, respectively, and w k and v k are the state noise and measurement noise vectors, respectively.Based on Equation (5), we can obtain

Table 3 .
Parameters used in the numerical model.

Table 4 .
Parameters for the offshore wind turbine installation site.

Table 5 .
Modal results obtained by N4SID (first ten orders).• and × correspond to stable or unstable modal results (poles).

Table 6 .
Comparison of compression capacity under different initial parameter cases (kN).

Table 7 .
Comparison of displacement (cm) and rotation angle (degree) under different initial parameter cases.

Table 8 .
Comparison of maximum displacement (cm) and rotation angle (degree) between different initial parameter cases and the design values.