Skip to main content
Records of precipitation in Macau are collected from three local monitoring stations from 1993 to 2003. Chemical compositions of precipitation are analyzed to reveal the causes of wet deposition inMacau. The review of volume weighted... more
Records of precipitation in Macau are collected from three local monitoring stations from 1993 to 2003. Chemical compositions of precipitation are analyzed to reveal the causes of wet deposition inMacau. The review of volume weighted annual pH values over this period shows that the entire area of Macau 2 is subjected to acidic precipitation. The concentrations of three ionic species H+, SO42- and NO3- are analyzed statistically and their distributions are found to have positive skewness and kurtosis. Lognormal cumulative distribution function fits each sample cumulative distribution function for each ionic species well implying that their temporal behaviours are governed according to the lognormal distribution. Nonetheless, the distribution of each ionic species exhibits spatial variability across the stations. Therefore, it is suspected that the sources of wet acid deposition indifferent parts of Macau are different and further correlation analysis on the values of pH, sulfate, and...
Abstract Extended Kalman filter (EKF) is widely adopted for state estimation and parametric identification of dynamical systems. In this algorithm, it is required to specify the covariance matrices of the process noise and measurement... more
Abstract Extended Kalman filter (EKF) is widely adopted for state estimation and parametric identification of dynamical systems. In this algorithm, it is required to specify the covariance matrices of the process noise and measurement noise based on prior knowledge. However, improper assignment of these noise covariance matrices leads to unreliable estimation and misleading uncertainty estimation on the system state and model parameters. Furthermore, it may induce diverging estimation. To resolve these problems, we propose a Bayesian probabilistic algorithm for online estimation of the noise parameters which are used to characterize the noise covariance matrices. There are three major appealing features of the proposed approach. First, it resolves the divergence problem in the conventional usage of EKF due to improper choice of the noise covariance matrices. Second, the proposed approach ensures the reliability of the uncertainty quantification. Finally, since the noise parameters are allowed to be time-varying, nonstationary process noise and/or measurement noise are explicitly taken into account. Examples using stationary/nonstationary response of linear/nonlinear time-varying dynamical systems are presented to demonstrate the efficacy of the proposed approach. Furthermore, comparison with the conventional usage of EKF will be provided to reveal the necessity of the proposed approach for reliable model updating and uncertainty quantification.
The soil–water characteristic curve (SWCC) of fine-grained soils is usually determined experimentally. In the design of mine waste covers and landfill liners, the unsaturated hydraulic conductivity function, k(h), is often derived... more
The soil–water characteristic curve (SWCC) of fine-grained soils is usually determined experimentally. In the design of mine waste covers and landfill liners, the unsaturated hydraulic conductivity function, k(h), is often derived theoretically from the measured SWCC. Implicit in these derivations is the transformation of the SWCC to a pore-size distribution (PSD), typically assumed to be constant and monomodal. However, PSD measurements of a clayey till compacted at various water contents after compaction, after flexible-wall permeability testing and before and after SWCC tests show that the PSD of the same material varies significantly under the stated physical conditions. Predictions of the SWCCs using PSDs measured both before and after the SWCC tests significantly underpredicted the values measured. By applying a simple transformation to the PSD to account for the scaling effect from the porosimetry samples (approximately 1 g dry weight) to the SWCC test samples (approximately 200 g dry weight), the ...
AbstractThe pullout resistance of soil nail is a key parameter in soil nailing design. Pullout resistance is affected by many factors, such as overburden pressure, grouting pressure, soil dilation, and degree of saturation of soil.... more
AbstractThe pullout resistance of soil nail is a key parameter in soil nailing design. Pullout resistance is affected by many factors, such as overburden pressure, grouting pressure, soil dilation, and degree of saturation of soil. Because of the complexity of the pullout mechanism, some factors have not been well incorporated in the current soil nail design methods. In this study, Bayesian analysis is performed to investigate the relative importance of several key factors and to build a new design formula to estimate maximum pullout shear stress of grouted soil nails. By using a series of laboratory soil nail pullout test data, Bayesian analysis is performed to select a predictive formula with suitable complexity and to identify its parameters. It is found that the most important factors are the degree of saturation and the product of grouting pressure and overburden pressure. It is shown that the proposed optimal model exhibits significantly stronger correlation with measurements than the existing effec...
Evaluation of the cyclic shear modulus of soils is a crucial but challenging task for many geotechnical earthquake engineering and soil dynamic issues. Improper determination of this property unnecessarily drives up design and maintenance... more
Evaluation of the cyclic shear modulus of soils is a crucial but challenging task for many geotechnical earthquake engineering and soil dynamic issues. Improper determination of this property unnecessarily drives up design and maintenance costs or even leads to the construction of unsafe structures. Due to the complexities involved in the direct measurement, empirical curves for estimating the cyclic shear modulus have been commonly adopted in practice for simplicity and economical considerations. However, a systematic and robust approach for formulating a reliable model and empirical curve for cyclic shear modulus prediction for clayey soils is still lacking. In this study, the Bayesian model class selection approach is utilized to identify the most significant soil parameters affecting the normalized cyclic shear modulus and a reliable predictive model for normally to moderately over-consolidated clays is proposed. Results show that the predictability and reliability of the proposed model out performs the well-known empirical models. Finally, a new design chart is established for practical usage.
The present study introduces the novel Bayesian approach to solve a difficult problem that modellers face when doing linear regression for air quality prediction; i.e. the non-uniqueness of the input variables and the functional forms to... more
The present study introduces the novel Bayesian approach to solve a difficult problem that modellers face when doing linear regression for air quality prediction; i.e. the non-uniqueness of the input variables and the functional forms to be selected. Using the historical information of ozone, nitrogen dioxide, respirable particulates and eight meteorological elements recorded during the high ozone seasons (May-October) in Macau between 2006 and 2007 as the training data, the Bayesian approach was applied for model selection from 16383 candidates. The model with the best efficiency-robustness balance selected is not the most complicated one. It was then examined against the most complicated model with the data between 2008 and 2009. Results show that the selected model yields better performance with the root-mean-squared error (RMSE) and the coefficient of determination (r2) equal to 19.18 µg/m3 and 0.81, respectively. In addition, it is capable to capture 87% of the ozone episodes w...
The aim of this paper is to present a novel attempt for parametric estimation in the hydrostatic-season-time (HST) model. The empirical HST-model has been widely used for the analysis of different measurement data types on dams. The... more
The aim of this paper is to present a novel attempt for parametric estimation in the hydrostatic-season-time (HST) model. The empirical HST-model has been widely used for the analysis of different measurement data types on dams. The significance of individual parameters or their sub-groups for modelling the influence of the water level, air and water temperature, and irreversible deformations due to the ageing of the dam, depends on the structure itself. The process of finding an accurate HST-model for a given data set, which remains robust to outliers, cannot only be demanding but also time consuming. The Bayesian model class selection approach imposes a penalisation against overly complex model candidates and admits a selection of the most plausible HST-model according to the maximum value of model evidence provided by the data or relative plausibility within a set of model class candidates. The potential of Bayes interference and its efficiency in an HST-model are presented on geodetic time series as a result of a permanent monitoring system on a rock-fill embankment dam. The method offers high potential for engineers in the decision making process, whilst the HST-model can be promptly adapted to new information given by new measurements and can enhance the safety and reliability of dams.
This study presents a method of predicting the soil water retention curve (SWRC) of a soil using a set of measured SWRC data from a soil with the same texture but different initial void ratio. The relationships of the volumetric water... more
This study presents a method of predicting the soil water retention curve (SWRC) of a soil using a set of measured SWRC data from a soil with the same texture but different initial void ratio. The relationships of the volumetric water contents and the matric suctions between two samples with different initial void ratios are established. An adjustment parameter (β) is introduced to express the relationships between the matric suctions of two soil samples. The parameter β is a function of the initial void ratio, matric suction or volumetric water content. The function can take different forms, resulting in different predictive models. The optimal predictive models of β are determined for coarse-grained and fine-grained soils using the Bayesian method. The optimal models of β are validated by comparing the estimated matric suction and measured data. The comparisons show that the proposed method produces more accurate SWRCs than do other models for both coarse-grained and fine-grained soils. Furthermore, the influence of the model parameters of β on the predicted matric suction and SWRC is evaluated using Latin Hypercube sampling. An uncertainty analysis shows that the reliability of the predicted SWRC decreases with decreasing water content in fine-grained soils, and the initial void ratio has no apparent influence on the reliability of the predicted SWRCs in coarse-grained and fine-grained soils.
A probabilistic approach for damage detection is presented using noisy incomplete input and response measurements that is an extension of a Bayesian system identification approach developed by the authors. This situation may be... more
A probabilistic approach for damage detection is presented using noisy incomplete input and response measurements that is an extension of a Bayesian system identification approach developed by the authors. This situation may be encountered, for example, during low-level ambient vibrations when a structure is instrumented with accelerometers that measure the input ground motion and structural response at a few locations but the wind excitation is not measured. A substructuring approach is used for the parameterization of the mass and stiffness distributions. Damage is defined to be a reduction of the substructure stiffness parameters compared with those of the undamaged structure. By using the proposed probabilistic methodology, the probability of various damage levels in each substructure can be calculated based on the available data. A four-story benchmark building subjected to wind and ground shaking is considered in order to demonstrate the proposed approach.
Outlier detection is an important problem in statistics. In this paper, we introduce a novel concept of outlier probability for outlier detection and robust linear regression. First, the Mahalanobis distance is utilized to identify the... more
Outlier detection is an important problem in statistics. In this paper, we introduce a novel concept of outlier probability for outlier detection and robust linear regression. First, the Mahalanobis distance is utilized to identify the leverage points. By excluding the leverage points, the maximum trimmed likelihood estimation will be used to obtain an initial set of regular data points while the remaining and the leverage points are included in the initial suspicious data set. Then, each suspicious data point and the combinations are evaluated by a novel outlier probability that depends not only on the residuals but also the size of the data set. Incorporating the data size is important as it controls the probability of the existence of a data point (or a batch or data points) exceeding a given value of the normalized residual. This outlier probability is robust because it incorporates also the posterior uncertainty quantified using the Bayesian approach. Then, the data points with...
The Pearl River Delta (PRD) region is located on the southeast coast of mainland China and it is an important economic hub. The high levels of particulate matter (PM) in the atmosphere, however, and poor visibility have become a complex... more
The Pearl River Delta (PRD) region is located on the southeast coast of mainland China and it is an important economic hub. The high levels of particulate matter (PM) in the atmosphere, however, and poor visibility have become a complex environmental problem for the region. Air quality modeling systems are useful to understand the temporal and spatial distribution of air pollution, making use of atmospheric emission data as inputs. Over the years, several atmospheric emission inventories have been developed for the Asia region. The main purpose of this work is to evaluate the performance of the air quality modeling system for simulating PM concentrations over the PRD using three atmospheric emission inventories (i.e., EDGAR, REAS and MIX) during a winter and a summer period. In general, there is a tendency to underestimate PM levels, but results based on the EDGAR emission inventory show slightly better accuracy. However, improvements in the spatial and temporal disaggregation of em...
Abstract This study was devoted to investigating stochastic model updating in a Bayesian inference framework based on a frequency response function (FRF) vector without any post-processing such as smoothing and windowing. The statistics... more
Abstract This study was devoted to investigating stochastic model updating in a Bayesian inference framework based on a frequency response function (FRF) vector without any post-processing such as smoothing and windowing. The statistics of raw FRFs were inferred with a multivariate complex-valued Gaussian ratio distribution. The likelihood function was formulated by embedding the theoretical FRFs that contained the model parameters to be updated in the class of the probability model of the raw FRFs. The Transitional Markov chain Monte Carlo (TMCMC) used to sample the posterior probability density function implies considerable computational toll because of the large batch of repetitive analyses of the forward model and the increasing expense of the likelihood function calculations with large-scale loop operations. The vectorized formula was derived analytically to avoid time-consuming loop operations involved in the likelihood function evaluation. Furthermore, a distributed parallel computing scheme was developed to allow the TMCMC stochastic simulation to run across multiple CPU cores on multiple computers in a network. The case studies demonstrated that the fast-computational scheme could exploit the availability of high-performance computing facilities to drastically reduce the time-to-solution. Finally, parametric analysis was utilized to illustrate the uncertainty propagation properties of the model parameters with the variations of the noise level, sampling time, and frequency bandwidth.
Abstract Reliable verification and evaluation of the mechanical properties of a layered composite ensemble are critical for industrially relevant applications, however it still remains an open engineering challenge. In this study, a fast... more
Abstract Reliable verification and evaluation of the mechanical properties of a layered composite ensemble are critical for industrially relevant applications, however it still remains an open engineering challenge. In this study, a fast Bayesian inference scheme based on multi-frequency single shot measurements of wave propagation characteristics is developed to overcome the limitations of ill-conditioning and non-uniqueness associated with the conventional approaches. A Transitional Markov chain Monte Carlo (TMCMC) algorithm is employed for the sampling process. A Wave and Finite Element (WFE)-assisted metamodeling scheme in lieu of expensive-to-evaluate explicit FE analysis is proposed to cope with the high computational cost involved in TMCMC sampling. For this, the Kriging predictor providing a surrogate mapping between the probability spaces of the model predictions for the wave characteristics and the mechanical properties in the likelihood evaluations is established based on the training outputs computed using a WFE forward solver, coupling periodic structure theory to conventional FE. The valuable uncertainty information of the prediction variance introduced by the use of a surrogate model is also properly taken into account when estimating the parameters’ posterior probability distribution by TMCMC. A numerical study as well as an experimental study are conducted to verify the computational efficiency and accuracy of the proposed methodology. Results show that the TMCMC algorithm in conjunction with the WFE forward solver-aided metamodeling can sample the posterior Probability Density Function (PDF) of the updated parameters at a very reasonable cost. This approach is capable of quantifying the uncertainties of recovered independent characteristics for each layer of the composite structure under investigation through fast and inexpensive experimental measurements on localized portions of the structure.
AbstractEstimation on the uniaxial compressive strength (UCS) of rock is an important issue in geotechnical engineering. Empirical relation establishment for UCS estimation is particularly favorabl...
AbstractNatural fibers are environment-friendly and efficient for soil reinforcement. Many studies have reported the influences of fiber percentage on the shear strength of fiber reinforced soil. H...
AbstractStructural health indicators, such as modal frequencies, have been commonly utilized to interpret the health condition of monitored structures. This study modeled the relationship between s...
Abstract Of late, various constitutive models have been proposed in the literature for the purpose of capturing the various complex physical mechanisms governing the creep behavior of soft soil. However, the more complex the model, the... more
Abstract Of late, various constitutive models have been proposed in the literature for the purpose of capturing the various complex physical mechanisms governing the creep behavior of soft soil. However, the more complex the model, the greater the number of associated uncertain parameters it has, and the less robustness it is. In this study, the Bayesian model class selection approach is applied to select the most plausible/suitable model describing the creep behavior of soft soil using laboratory measurements. In total, one elastic plastic (EP) model and eight elastic viscoplastic (EVP) models are investigated. To assess the performance of the different models in the prediction of creep behavior of soft soils, Bayesian model class selection is respectively performed using the oedometer test data from the intact samples of Vanttila clay and reconstituted samples of Hong Kong Marine Clay collected from the literature. All unknown model parameters are identified simultaneously by adopting the transitional Markov Chain Monte Carlo (TMCMC) method, and their uncertainty is quantified through the obtained posterior probability density functions (PDFs). The result shows that the proposed method is an excellent candidate for identifying the most plausible model and its associated parameters for different kinds of soft soils. The approach also provides uncertainty evaluation of the model prediction based on the given data.
In ground motion prediction, the key is to develop a suitable and reliable GMPE (ground motion prediction equation) characterizing the ground motion pattern of the target seismic region. There are two critical goals encountered in GMPE... more
In ground motion prediction, the key is to develop a suitable and reliable GMPE (ground motion prediction equation) characterizing the ground motion pattern of the target seismic region. There are two critical goals encountered in GMPE development. Proposing a suitable predictive formula applicable to target seismic region has attracted much of the attention in previous studies. On the other hand, dependence between prediction–error variance and ground motion data has been observed and the study on this kind of heterogeneous relation becomes an important task yet to be explored. In this article, a novel HEteRogeneous BAyesian Learning (HERBAL) approach is proposed for achieving these two goals simultaneously. The homogeneity assumption on error in the traditional learning approach is relaxed, so the proposed approach is applicable for more general heterogeneous cases. With the generalization made on the traditional Bayesian learning by embedding the derived closed form expression for error variance parameter optimization component into the hyperparameter optimization of ARD (automatic relevance determination) prior, the proposed learning approach is capable of performing continuous model training on a prescribed predictive formula with unknown error pattern. A database of strong ground motion records in the Tangshan region of China is obtained for the analysis. It is shown that the trained optimal model class by the proposed approach is promising as that, the trained optimal model class retains model simplicity of the predictive formula with capability on both robustness enhancement ground motion prediction and precise determination of the error pattern.
Abstract The techniques of bias correction are commonly used for improving the performance of deterministic air quality forecasting systems. One issue not addressed in previous studies is how to select systematically and objectively the... more
Abstract The techniques of bias correction are commonly used for improving the performance of deterministic air quality forecasting systems. One issue not addressed in previous studies is how to select systematically and objectively the best correction model from a pool of candidates. In this study, a method that could evaluate the probabilities of all model candidates based on a set of training data is proposed to select the most accurate and robust model by finding the one with the maximum probability. The Bayesian method was applied to select bias correction models at 12 monitoring stations for improving the forecasts of daily averaged PM 10 concentrations given by the deterministic air quality forecasting system WRF-EURAD in Porto, Portugal. At each station, 4095 (=2 12 –1) correction model candidates were systematically formed by adopting different linear combinations of 12 input variables. Selection of the best model was processed based on one year of monitoring and WRF-EURAD data. Based on the 2012 data, the selected model at each station was found to have significantly higher probability than the other candidates, and it is also much simpler than the full model. These selected models were then used to correct the raw forecasts by WRF-EURAD in the following year. The corrected forecasts show significant improvement on the performance indicators ( RMSE by 35.8%, R by 58.5%, MFB by 68%, EDR by 38.3%, FAR by 51.8%, CSI by 30.8%) over the raw outputs of WRF-EURAD, confirming the success of the proposed technique.
Modal frequency is an important indicator for structural health assessment. Previous studies have shown that this indicator is substantially affected by the fluctuation of ambient conditions, such as temperature and humidity. Therefore,... more
Modal frequency is an important indicator for structural health assessment. Previous studies have shown that this indicator is substantially affected by the fluctuation of ambient conditions, such as temperature and humidity. Therefore, recognizing the pattern between modal frequency and ambient conditions is necessary for reliable long-term structural health assessment. In this article, a novel machine-learning algorithm is proposed to automatically select relevance features in modal frequency-ambient condition pattern recognition based on structural dynamic response and ambient condition measurement. In contrast to the traditional feature selection approaches by examining a large number of combinations of extracted features, the proposed algorithm conducts continuous relevance feature selection by introducing a sophisticated hyperparameterization on the weight parameter vector controlling the relevancy of different features in the prediction model. The proposed algorithm is then u...
The undrained shear strength (su) of cohesive soils is a crucial parameter for many geotechnical engineering applications. Due to the complexities and uncertainties associated with laboratory and in situ tests, it is a challenging task to... more
The undrained shear strength (su) of cohesive soils is a crucial parameter for many geotechnical engineering applications. Due to the complexities and uncertainties associated with laboratory and in situ tests, it is a challenging task to obtain the undrained shear strength in a reliable and economical manner. In this study, a probabilistic model for the su of moderately overconsolidated clays is developed using the Bayesian model class selection approach. The model is based on a comprehensive geotechnical database compiled for this study with field measurements of field vane strength (su), plastic limit (PL), natural water content (Wn), liquid limit (LL), vertical effective overburden stress ($$\sigma_{\nu }^{\prime }$$σν′), preconsolidation pressure ($$\sigma_{\text{p}}^{\prime }$$σp′) and overconsolidated ratio (OCR). Comparison study shows that the proposed model is superior to some well-known empirical relationships for OC clays. The proposed probabilistic model not only provides reliable and economical estimation of su but also facilitates reliability-based analysis and design for performance-based engineering applications.
This study compared two types of (offline vs. online) bias correction models applied to correct the mismatch between the predicted daily averaged PM10 concentration by the deterministic air quality forecasting system WRF-EURAD and the... more
This study compared two types of (offline vs. online) bias correction models applied to correct the mismatch between the predicted daily averaged PM10 concentration by the deterministic air quality forecasting system WRF-EURAD and the measured concentration at the air quality monitoring station in Porto, Portugal. The WRF-EURAD is a Eulerian system consisting of a Weather Research Forecasting (WRF) model and a European Air Pollution Dispersion (EURAD) model. Both bias correction models were linear statistical models developed with the same set of input variables. The major difference between the online or the offline models is the adaptiveness. While the coefficients of the offline bias correction model are fixed after training with the ordinary least squares, those in the online bias correction model are updated adaptively with the Kalman filter. Comparison of these bias correction models was made at an urban traffic station Senhora da Hora in Porto, Portugal within 2013. The fract...
A daily PM2.5 forecasting model based on multiple linear regression (MLR) and backward trajectory clustering of HYSPLIT was designed for its application to small cities where PM2.5 level is easily affected by regional transport. The... more
A daily PM2.5 forecasting model based on multiple linear regression (MLR) and backward trajectory clustering of HYSPLIT was designed for its application to small cities where PM2.5 level is easily affected by regional transport. The objective of this study is to investigate the regions that affect the fine particulate concentration of Macau and to develop an effective forecasting system to enhance the capture of PM2.5 episodes. By clustering the HYSPLIT 24-hr backward trajectories originated at Macau from 2015 to 2017, five potential transportation paths of PM2.5 were found. A cluster based statistical model was developed and trained with air quality and meteorological data of2015 and 2016. Then, the trained model was evaluated with data of 2017. Comparing to an ordinary model without backward trajectory clustering, the cluster based PM2.5 forecasting model yielded similar general forecast performance in 2017. However, the critical success index of the cluster based model was 11% hi...
This paper presents a comprehensive literature review on point set registration. The state-of-the-art modeling methods and algorithms for point set registration are discussed and summarized. Special attention is paid to methods for... more
This paper presents a comprehensive literature review on point set registration. The state-of-the-art modeling methods and algorithms for point set registration are discussed and summarized. Special attention is paid to methods for pairwise registration and groupwise registration. Some of the most prominent representative methods are selected to conduct qualitative and quantitative experiments. From the experiments we have conducted on 2D and 3D data, CPD-GL pairwise registration algorithm and JRMPC groupwise registration algorithm seem to outperform their rivals both in accuracy and computational complexity. Furthermore, future research directions and avenues in the area are identified.
In this article, a novel propagative broad learning approach is proposed for nonparametric modeling of the ambient effects on structural health indicators. Structural health indicators interpret the structural health condition of the... more
In this article, a novel propagative broad learning approach is proposed for nonparametric modeling of the ambient effects on structural health indicators. Structural health indicators interpret the structural health condition of the underlying dynamical system. Long-term structural health monitoring on in-service civil engineering infrastructures has demonstrated that commonly used structural health indicators, such as modal frequencies, depend on the ambient conditions. Therefore, it is crucial to detrend the ambient effects on the structural health indicators for reliable judgment on the variation of structural integrity. However, two major challenging problems are encountered. First, it is not trivial to formulate an appropriate parametric expression for the complicated relationship between the operating conditions and the structural health indicators. Second, since continuous data stream is generated during long-term structural health monitoring, it is required to handle the gr...
It is a normal practice to consider the reduction of modal frequencies as an indicator for structural damage but some long-term monitoring studies revealed that the structural modal frequencies exhibited substantial seasonal variation... more
It is a normal practice to consider the reduction of modal frequencies as an indicator for structural damage but some long-term monitoring studies revealed that the structural modal frequencies exhibited substantial seasonal variation without evidence of structural damage. Therefore, in order to perform reliable and consistent assessment on the health condition of civil engineering infrastructures, it is an indispensable task to quantify the modal frequency fluctuation induced by change of ambient conditions, such as temperature and humidity. In other words, it is necessary to distinguish the inevitable ambient interference from the abnormal changes due to structural damage and deterioration. In this paper, one-year daily measurement of a 22-storey reinforced concrete building is utilized to trace the variation of its modal frequencies, which are identified using the Bayesian spectral density approach with the ambient vibration data. Then, the Timoshenko beam model is considered to construct the mathematical model for the modal frequency-ambient condition relationship. Since the Timoshenko beam model incorporates both the shear and flexural deformation, it is appropriate for the behavior of buildings with different aspect ratios and configurations. Finally, Bayesian analysis is conducted to quantify the uncertain parameters in the modal frequency-ambient condition model. Results show that explicit consideration of the ambient temperature and relative humidity is essential for long-term structural health monitoring.
In this article, a novel Bayesian framework is proposed for real‐time system identification with calibratable model classes. This self‐calibrating scheme adaptively reconfigures the model classes to achieve reliable real‐time estimation... more
In this article, a novel Bayesian framework is proposed for real‐time system identification with calibratable model classes. This self‐calibrating scheme adaptively reconfigures the model classes to achieve reliable real‐time estimation for the system state and model parameters. At each time step, the plausibilities of the model classes are computed and they serve as the cue for calibration. Once calibration is triggered, all model classes will be reconfigured. Thereafter, identification will continue to propagate with the calibrated model classes until the next recalibration. Consequently, the model classes will evolve and their deficiencies can be corrected adaptively. This remarkable feature of the proposed framework stimulates the accessibility of reliable real‐time system identification. Examples are presented to demonstrate the efficacy of the proposed approach using noisy response measurement of linear and nonlinear time‐varying dynamical systems under stationary condition.
... 989 Market Street, San Francisco, CA 94103-1741, USA Wiley-VCH Verlag GmbH, Boschstrasse 12, D-69469 Weinheim, Germany John Wiley & Sons Australia Ltd, 42 McDougall Street, Milton ... To my family To my wife,... more
... 989 Market Street, San Francisco, CA 94103-1741, USA Wiley-VCH Verlag GmbH, Boschstrasse 12, D-69469 Weinheim, Germany John Wiley & Sons Australia Ltd, 42 McDougall Street, Milton ... To my family To my wife, Vanessa, my daughter, Hoi-Man and my son, Kin-Hang ...
... I would also like to thank my committee members, Prof. Joel W. Burdick, Prof. John F. Hall, Prof. ... coefficients and participation factors, of a linear multiple-degree-of-freedom (MDOF) system (Gersch and Foutch 1974 Beck 1978... more
... I would also like to thank my committee members, Prof. Joel W. Burdick, Prof. John F. Hall, Prof. ... coefficients and participation factors, of a linear multiple-degree-of-freedom (MDOF) system (Gersch and Foutch 1974 Beck 1978 Hoshiya and Saito 1984 Quek et al. 1999 Shi et al. ...
ABSTRACT
In the present study, the Bayesian probabilistic approach is proposed to estimate the process noise and measurement noise parameters in the Kalman filter for the case when the input is a zero-mean Gaussian white noise process and limited... more
In the present study, the Bayesian probabilistic approach is proposed to estimate the process noise and measurement noise parameters in the Kalman filter for the case when the input is a zero-mean Gaussian white noise process and limited output measurements are available. The methodology presented in this study is based on the idea that the process noise and measurement noise parameters directly affect the state estimate and the covariance matrix of the Kalman filter. Therefore, the noise parameters can be estimated by expressing the likelihood factor of the measurements as a function of the state estimate and the covariance matrix for a given set of process noise and measurement noise parameters. By evaluating the likelihood factor for all the noise parameters within the search domain, the optimal estimates of the noise parameters are chosen by the maximum likelihood criterion. Through the demonstration of an example, the optimal estimates of the noise parameters are close to the a...
Outlier detection is an important problem in statistics. In this paper, we introduce a novel concept of outlier probability for outlier detection and robust linear regression. First, the Mahalanobis distance is utilized to identify the... more
Outlier detection is an important problem in statistics. In this paper, we introduce a novel concept of outlier probability for outlier detection and robust linear regression. First, the Mahalanobis distance is utilized to identify the leverage points. By excluding the leverage points, the maximum trimmed likelihood estimation will be used to obtain an initial set of regular data points while the remaining and the leverage points are included in the initial suspicious data set. Then, each suspicious data point and the combinations are evaluated by a novel outlier probability that depends not only on the residuals but also the size of the data set. Incorporating the data size is important as it controls the probability of the existence of a data point (or a batch or data points) exceeding a given value of the normalized residual. This outlier probability is robust because it incorporates also the posterior uncertainty quantified using the Bayesian approach. Then, the data points with...
Spatial and temporal variation of fine particulate pollution in a complex small coastal city, Macau, was investigated by using a mobile air quality monitoring platform. Real-time street level measurements of aerosol size distribution were... more
Spatial and temporal variation of fine particulate pollution in a complex small coastal city, Macau, was investigated by using a mobile air quality monitoring platform. Real-time street level measurements of aerosol size distribution were carried out three times per day by using an on-board wide-range particle spectrometer along a designed route covering the main area of the city for one week. The measured distributions were then transformed into PM2.5 mass concentrations for assessment based on an empirical relation estimated from a co-located measurement calibration process. Within this week, the averaged PM2.5 mass concentrations of the morning, afternoon, and evening sessions are 51.76 μg m-3, 28.87 μg m-3 and 30.77 μg m-3, respectively. By comparing the spatial statistics of the estimated PM2.5 concentrations with the US 24-hour primary NAAQS as a reference, a significant probability of exceedance (> 86%) was observed in 1/3 of the monitoring sessions, particularly in the mo...
The present study introduces the novel Bayesian approach to solve a difficult problem that modellers face when doing linear regression for air quality prediction; i.e. the non-uniqueness of the input variables and the functional forms to... more
The present study introduces the novel Bayesian approach to solve a difficult problem that modellers face when doing linear regression for air quality prediction; i.e. the non-uniqueness of the input variables and the functional forms to be selected. Using the historical information of ozone, nitrogen dioxide, respirable particulates and eight meteorological elements recorded during the high ozone seasons (May-October) in Macau between 2006 and 2007 as the training data, the Bayesian approach was applied for model selection from 16383 candidates. The model with the best efficiency-robustness balance selected is not the most complicated one. It was then examined against the most complicated model with the data between 2008 and 2009. Results show that the selected model yields better performance with the root-mean-squared error (RMSE) and the coefficient of determination (r2) equal to 19.18 μg/m3 and 0.81, respectively. In addition, it is capable to capture 87% of the ozone episodes w...
The goal of this study is to investigate the structural performance of reinforced concrete building under the influence of severe typhoon. For this purpose, full-scale monitoring of a 22-story reinforced concrete building was conducted... more
The goal of this study is to investigate the structural performance of reinforced concrete building under the influence of severe typhoon. For this purpose, full-scale monitoring of a 22-story reinforced concrete building was conducted during the entire passage process of a severe typhoon "Vicente." Vicente was the eighth tropical storm developed in the Western North Pacific Ocean and the South China Sea in 2012. Moreover, it was the strongest and most devastating typhoon that struck Macao since 1999. The overall duration of the typhoon affected period that lasted more than 70 hours and the typhoon eye region covered Macao for around one hour. The wind and structural response measurements were acquired throughout the entire typhoon affected period. The wind characteristics were analyzed using the measured wind data including the wind speed and wind direction time histories. Besides, the structural response measurements of the monitored building were utilized for modal iden...
Estimation of peak ground acceleration is one of the main issues in civil and earthquake engineering practice. The Boore-Joyner-Fumal empirical formula [1] is well known for this purpose. In this paper we propose to use the Bayesian... more
Estimation of peak ground acceleration is one of the main issues in civil and earthquake engineering practice. The Boore-Joyner-Fumal empirical formula [1] is well known for this purpose. In this paper we propose to use the Bayesian probabilistic model class selection approach to obtain the most suitable prediction model class for the seismic attenuation formula. The optimal model class is
ABSTRACT A Bayesian probabilistic method is proposed for online estimation of the process noise and measurement noise parameters for Kalman filter. Kalman filter is a well-known recursive algorithm for state estimation of dynamical... more
ABSTRACT A Bayesian probabilistic method is proposed for online estimation of the process noise and measurement noise parameters for Kalman filter. Kalman filter is a well-known recursive algorithm for state estimation of dynamical systems. In this algorithm, it is required to prescribe the covariance matrices of the process noise and measurement noise. However, inappropriate choice of these covariance matrices substantially deteriorates the performance of the Kalman filter. In this paper, a probabilistic method is proposed for online estimation of the noise parameters which govern the noise covariance matrices. The proposed Bayesian method not only estimates the optimal noise parameters but also quantifies the associated estimation uncertainty in an online manner. By utilizing the estimated noise parameters, reliable state estimation can be accomplished. Moreover, the proposed method does not assume any stationarity condition of the process noise and/or measurement noise. By removing the stationarity constraint, the proposed method enhances the applicability of the state estimation algorithm for nonstationary circumstances generally encountered in practice. To illustrate the efficacy and efficiency of the proposed method, examples using a fifty-story building with different stationarity scenarios of the process noise and measurement noise are presented.
Research Interests:
Research Interests:
ABSTRACT
A probabilistic approach for damage detection is presented using noisy incomplete input and response measurements that is an extension of a Bayesian system identification approach developed by the authors. This situation may be... more
A probabilistic approach for damage detection is presented using noisy incomplete input and response measurements that is an extension of a Bayesian system identification approach developed by the authors. This situation may be encountered, for example, during low-level ambi- ent vibrations when a structure is instrumented with accelerometers that measure the input ground motion and structural response at a few
... to Macau. In fact, previous studies have shown that the change of wind direction can affect the pollutant concentrations of the Pearl River Delta region (Lee & Hills, 2003; Mok & Hoi, 2003; Chang et al., 2007).... more
... to Macau. In fact, previous studies have shown that the change of wind direction can affect the pollutant concentrations of the Pearl River Delta region (Lee & Hills, 2003; Mok & Hoi, 2003; Chang et al., 2007). In addition, this ...
ABSTRACT
... Permissions & Reprints. An efficient simulation method for reliability analysis of linear dynamical systems using simple additive rules of probability. Ka-Veng Yuen a and Lambros S. Katafygiotis b , Corresponding Author... more
... Permissions & Reprints. An efficient simulation method for reliability analysis of linear dynamical systems using simple additive rules of probability. Ka-Veng Yuen a and Lambros S. Katafygiotis b , Corresponding Author Contact Information , E-mail The Corresponding Author. ...
... Single-Degree-of-Freedom Systems. Bayesian System Identification Formulation. ... Consider a system with degrees of freedom (DOFs) whose displacement response is modeled using the equation of motion: where (known) mass matrix,... more
... Single-Degree-of-Freedom Systems. Bayesian System Identification Formulation. ... Consider a system with degrees of freedom (DOFs) whose displacement response is modeled using the equation of motion: where (known) mass matrix, nonlinear restoring force ...
ABSTRACT Increase in pile capacity after initial driving has been well observed in clays and sands over decades. The phenomenon is referred to as pile set-up by geotechnical engineers. More economical pile design may benefit from this... more
ABSTRACT Increase in pile capacity after initial driving has been well observed in clays and sands over decades. The phenomenon is referred to as pile set-up by geotechnical engineers. More economical pile design may benefit from this time-dependent increase subject to a reliable prediction. Simple empirical relations of the current capacity with the initial capacity and elapse time after driving are available in the literature with different model parameters being suggested for clays and sands, respectively. Nevertheless, appropriateness of the relations and confidence interval of the model parameters are rarely investigated and this hinders the application of these formulae. In this study, a revised single-parameter empirical relation is proposed based on the existing formulae. A comprehensive database from pile field tests data in clayey and sandy ground in the literature is compiled and Bayesian analysis is conducted on both these empirical formulae independently for clays and sands. Bayesian inference allows not only the estimation of the uncertain parameter but also the quantification of the associated uncertainty in the form of probability distribution. This study sheds lights on the confidence interval of the model parameter and it provides designers more reliable prediction of the additional capacity due to pile set-up.
ABSTRACT
... In an ambient vibration survey, the naturally occurring vibrations of the structure (due to wind, tra c, micro-tremors, etc.) are measured and then a system ... Many time-domain methods have been developed to tackle this problem. One... more
... In an ambient vibration survey, the naturally occurring vibrations of the structure (due to wind, tra c, micro-tremors, etc.) are measured and then a system ... Many time-domain methods have been developed to tackle this problem. One example is the random decrement ...