Extension of REBMIX algorithm to von Mises parametric family for modeling joint distribution of wind speed and direction
Introduction
For the long-span bridges, the environmental wind load will cause relatively large vibration and deformation with the increase of flexibility and span. As stress cycles would occur under frequent wind load induced by wind speed, critical components of long-span steel bridges may suffer an enormous amount of stress cycles, which will lead to fatigue damage accumulation and finally structural failure may appear without extreme climate condition. In addition to wind speed, wind direction is another important factor when considering the effect of wind load because wind in the along and transverse directions of a bridge will cause a completely different effect [23]. Due to the wind-induced fatigue damage of long-span steel bridges depending on wind speed and wind direction, the statistical model of field monitoring of wind data near the bridge site will lay a solid foundation for the accurate fatigue evaluation. A suitable model of joint distribution and an efficient parameter estimation method are major challenges for constructing the statistical model of wind speed and direction.
Nowadays, structural health monitoring systems have been installed on structures to measure environmental load and actual structural response of the existing structures. Field monitoring of wind data includes wind speed variable and wind direction variable. For the wind speed variable, the two-parameter Weibull distribution has been commonly used and accepted to model wind speed distribution [2]. Compared with the literature on wind speed, the literature on wind direction modeling is more limited because the distribution of the wind direction variable is a circular distribution [8]. Masseran et al. [12] used the circular distribution based on nonnegative trigonometric sums as the probability density function of wind direction. Carta et al. [5] proposed a finite mixture model comprised of a finite mixture of von Mises distribution to represent the distribution of directional wind speed. Based on the distribution functions of wind speed and direction, there are some approaches to construct the joint distribution model. Feng et al. [9] constructed the joint distribution of wind speed and direction based on the parameters of the sector-wise Weibull distribution and interpolations between direction sectors. Additionally, the angular-linear distribution models are used to construct joint distributions with specified marginal distributions [6]. While the models described above are all comprised of one-dimensional marginal distributions that will lead to the deviation between the predictive model and the distribution of observations, very little has been written in the literature about the finite mixture distribution models that are multidimensional.
Many optimization algorithms have been proposed to optimize the joint distribution models by estimating the parameters in models and selecting the optimal joint distribution model, such as the least square method [4], [5], the expectation maximization (EM) algorithm[21], [3] and the genetic algorithm [10]. Carta et al. [5] used the least squares method to estimate the parameters in a mixture of von Mises distributions. Masseran [13] separated wind direction data into eight groups based on geographical directions and used Markov chain models to model the distribution. Nasr et al. [20] applied three maximum likelihood technics based on Newton Raphson algorithm, Levenberg Marquardt algorithm and Trust Region reflective algorithm to the parameter estimation of the Weibull distribution. Calderara et al. [3] proposed a new specific EM algorithm for estimating the parameters in a mixture of von Mises distributions. In addition, the EM algorithm has been used widely to estimate parameters for the Weibull distribution [22], [15]. However, the EM algorithm is prone to converge to the local maximum of the likelihood function and has a limitation for the parameter estimation of multidimensional models.
The basic purpose of this paper is to develop an effective and efficient parameter estimation algorithm for the multivariate and multimodal circular distribution. The finite mixture of von Mises distributions is selected as the predictive model of the circular distribution. The unknown parameters in the predictive models are estimated by the extended REBMIX algorithm using an iterative process to assign the component densities one after another to the empirical mixture density. The optimal model will be selected from the models of different component numbers by judging the value of Akaike’s information criterion (AIC). This paper is organized as follows: the basic finite mixture model is introduced firstly in Section 2. Moreover, the mixture models of wind direction and the joint mixture models of wind speed and direction are formulated by the mixture of von Mises distributions and the mixture of Weibull-von Mises distributions, respectively. Then, Section 3 presents an extended REBMIX algorithm for estimating the unknown parameters in the predictive models and a detailed description of the flow of the algorithm. In Section 4, the extended algorithm is verified by being compared to the simulated data. Meanwhile, its performance of modeling with the field monitoring data, as compared to the performance calculated by the EM algorithm, is presented in Section 5. Finally, the conclusions are listed in Section 6.
Section snippets
Finite mixture distribution model
The finite mixture distribution model, in statistics, is a probability distribution of random variables whose populations contain two or more subpopulations. The mixture model represents subpopulations within the global population with no need to identify the subpopulation that each individual observation belongs to. The probability density function of the mixture model can be regarded as a convex combination of other different functions with non-negative weights, which sum to one. The
Parameter estimation framework
The REBMIX algorithm for conditionally independent normal, lognormal, Weibull, gamma, binomial, Poisson and Dirac component densities is explained in detail in Nagode [17], Nagode and Fajdiga [18] and Nagode and Fajdiga [19]. It is depicted in Fig. 1 and is implemented in a freely available R software package rebmix [16].
The paper is focused on the mixtures of von Mises component densities. In the REBMIX algorithm, parameter estimation is largely independent of parametric families. This means
Numerical examples
To verify the applicability and effectiveness of the proposed REBMIX algorithm in parameters estimation of finite mixture circular distribution models, two examples of numerical simulations have been tested here. Herein, the proposed algorithm is applied to estimate the parameters of mixture models without being compared to the EM algorithm. In example 1, the datasets of size n = 10,000 generated from the one, two and three components mixture of von Mises distributions are applied for
Field monitoring application
In order to demonstrate the effectiveness of the REBMIX algorithm in univariate and bivariate circular distribution modeling, one-year wind field monitoring data, including a total of 52,560 mean wind speed and wind direction data during ten minutes, measured by the ultrasonic anemometer in structural health monitoring system are employed for the construction of the distribution of wind direction and the joint distribution of wind speed and direction by the proposed algorithm and comparison
Conclusions
In this paper, an improved parameter estimation algorithm for estimating unknown parameters in multivariate and multimodal circular distributions is presented in detail. The predictive models select the finite mixture of von Mises distributions to characterize the distributions of circular variables. The unknown parameters for the mixture components are estimated through the proposed REBMIX algorithm. The estimated model with the lowest AIC value is selected as the optimal mixture model from
Acknowledgments
The work described in this paper was jointly supported by the National Science Foundation of China (Grant Nos. 51822810, 51778574), the Zhejiang Provincial Natural Science Foundation of China (Grant No. LR19E080002), and the Slovenian program Nr. P2-0182 entitled development evaluation financed by the Slovenian Ministry of Education, Science and Sport.
References (23)
- et al.
An investigation on wind energy potential and small scale wind turbine performance at İncek region - Ankara, Turkey
Energy Convers Manage
(2015) - et al.
Statistical modelling of directional wind speeds using mixtures of von Mises distributions: case study
Energy Convers Manage
(2008) - et al.
A joint probability density function of wind speed and direction for wind energy analysis
Energy Convers Manage
(2008) - et al.
ARMA based approaches for forecasting the tuple of wind speed and direction
Appl Energy
(2011) Markov chain model for the stochastic behaviors of wind-direction data
Energy Convers Manage
(2015)- et al.
Handbook of mathematical functions with formulas
(1972) - et al.
Mixtures of von Mises distributions for people trajectory shape analysis
Trans Circuits Syst Video Technol
(2011) - et al.
Measurement of the distribution of residual stresses in layered thick-walled GFRP pipes
Exp Mech
(2014) - et al.
Pattern classification and scene analysis
(1973) - et al.
Modelling wind for wind farm layout optimization using joint distribution of wind speed and wind direction
Energies
(2015)
Adaptation in natural and artificial system
Cited by (16)
Fault detection and classification with the rebmix R package
2023, Computers and Industrial EngineeringGumbel mixture modelling for multiple failure data
2023, Reliability Engineering and System SafetyCitation Excerpt :However, this assumes that origination failure mechanisms and the optimal separation rule are known in advance, which is often not the case [4,7,12]. Another popular solution approach is to formulate composite probability distributions, also known as mixture models [13–18]. Mixture models are constructed from several simple parametric distributions called “mixture components”.
Analysis of the wind field characteristics induced by the 2019 Typhoon Bailu for the high-speed railway bridge crossing China's southeast bay
2021, Journal of Wind Engineering and Industrial AerodynamicsA Swarm intelligence approach for statistical modeling of wind speed and direction: A case study of new york bight
2024, Computing in Civil Engineering 2023: Data, Sensing, and Analytics - Selected Papers from the ASCE International Conference on Computing in Civil Engineering 2023Joint Model of Wind Speed and Corresponding Direction Based on Wind Rose for Wind Energy Exploitation
2022, Journal of Ocean University of China