Abstract

It may not be possible to collect adequate records of strong ground motions in a short period of time; hence microtremor survey is frequently conducted to reveal the stratum structure and earthquake characteristics at a specified construction site. This paper is therefore aimed at developing a neural network model, based on available microtremor measurement and on-site soil boring test data, for predicting peak ground acceleration at a site, in a science park of Taiwan. The four key parameters used as inputs for the model are soil values of the standard penetration test, the medium grain size, the safety factor against liquefaction, and the distance between soil depth and measuring station. The results show that a neural network model with four neurons in the hidden layer can achieve better performance than other models presently available. Also, a weight-based neural network model is developed to provide reliable prediction of peak ground acceleration at an unmeasured site based on data at three nearby measuring stations. The method employed in this paper provides a new way to treat this type of seismic-related problem, and it may be applicable to other areas of interest around the world.

1. Introduction

Earthquake problems are globally considered to be a research topic of importance since many countries are subject to this natural disaster. For instances, the recent big one with magnitude 9.0 on the Richter scale that occurred in Japan on 11 March 2011 and triggered a significant tsunami, caused approximately $35 billion in damage. On 12 January 2010, a devastating earthquake with magnitude 7.0 on the Richter scale struck Haiti in the Caribbean, and claimed more than 200 thousand lives have lost in the capital and surrounding areas. Other major earthquakes experienced by Chile, China, Indonesia, New Zealand, and Taiwan are listed in the archive of United States Geological Survey [1]. Without exception, all these earthquakes have caused tremendous casualties and property losses, requiring urgent attention to this calamitous problem.

Many seismic-related research issues have been investigated and published previously, with some focused on finding an early warning system, while others are based on records of historical seismic data (e.g., [28]). It is quite obvious that strong ground motions data cannot be collected in a short period of time, and also the records for a metropolitan area or a place with high population density are not easy to obtain. Whereas, microtremor surveys can be used to infer the stratum structure and earthquake characteristics at a specified construction site without destroying its ground surface. Thus, this fast and low cost measuring technique is often selected not only to provide useful information for an area which lacks seismic records, but also to effectively analyze potential liquefaction index for the construction site being investigated.

Further, it is worth mentioning that the microtremor measurements with appropriate transformation, such as Fourier transform or Nakamura technique, can estimate the key seismic parameter, that is, peak ground acceleration (PGA), which exhibits a tendency similar to the characteristics of strong ground motion [9, 10]. Note that the above literatures were focused on the development of neural network models based on actual seismic records, and microtremor measurements were used for the sake of comparison. Regarding the development of neural network model in accordance with microtremor measurement and soil profile was not examined in these studies. Also, some important factors such as dominant frequency, shear wave speed, and amplification can be explored by microtremor surveys [1118], and these can help to determine the distribution of soil layers, liquefaction hazard mapping, and earthquake site response. Although microtremor measurements can be easily carried out at a number of sites, the main limitation is the increase in cost with increased number of measuring stations. Hence, the development of a model for predicting microtremor information for other important but unmeasured sites is useful for economic reasons.

From the references mentioned above and other reports previously published in the field of earthquake engineering and soil dynamics, the microtremor measurements appear to be a function of the soil conditions at a specified site. However, prediction of microtremor information by using soil boring test result has rarely been reported up to now. Therefore, the purpose of this study is to develop a model for mapping soil boring test data to microtremor measurements by using a neural network approach. In particular, three key soil parameters; the standard penetration test value (STP-N), the medium grain size (D50), and the safety factor against liquefaction (FL), and one spatial factor, that is, the distance between soil layer and measuring station (DS), are used to evaluate PGA resulting from microtremor measurements. A weight-based neural network model is also developed to predict PGA at unmeasured sites by using values at three nearby measuring stations. The method developed in this study should provide a new approach for solving problems in the relevant engineering field.

2. Context and Rationale for the Research

Science parks are mainly occupied by many high-tech companies including some world class factories such as ACER, HTC, and TSMC, which play an important economic role in the island of Taiwan. In this study, the chosen Kaohsiung (Luchu) science park is one of the major parks located in the southern part of Taiwan. This park has a total area of 571 hectares, which started construction in the year of 2001 and was completed in 2010. Various high-tech industries, such as integrated circuits, precision machines, optoelectronic components, computer peripherals, communication and biotechnology products, were planned and developed in this park. These types of industries can be affected significantly by strong ground motions and are also sensitive to ambient vibrations. Therefore, it is necessary to consider antiearthquake design and to examine microtremor in the park from time to time to prevent different levels of damages.

Two crucial factors for evaluating the effect of ground motions in the science park are fault distribution and geological condition in Kaohsiung area, which are based on the information from Central Geological Survey shown in Figure 1 [19, 20]. It can be seen that there exists seven faults in this region, which are (1) Chishan fault, (2) Liukuei fault, (3) Tsaujou fault, (4) Hsiaokangshan fault, (5) Yuchang fault, (6) Jenwu fault, and (7) Fengshan fault. These faults may create strong ground motions and endanger the high-tech buildings and instruments. Also from this figure, it can be seen that alluvial soil occupies a large part of Kaohsiung area, particularly at the science park, which can have an influence on microtremor measurements.

The occurrence of strong ground motion is unpredictable, and it can cause serious structural damage within a very short period of time. Thus, a proper antiearthquake design is usually considered for high-tech factory constructions. In contrast, the existence of microtremor is easy to neglect as it is very small, but due to microtremor occurring very often on the earth surface, some precision instruments can be damaged during its operation process due to the constant continuous vibration frequency or peak ground acceleration. Consequently, microtremor can cause, for high-tech manufactures, an unexpected and significant financial loss, and thus environmental ambient vibration survey is a very important consideration for science parks.

Figure 2 shows the Kaohsiung science park, with four microtremor measuring stations MS1, MS2, MS3, and MS4. Also, there are twenty-seven soil boring test sites in the neighborhood of this park. As mentioned previously, this study is focused on the development of a model for predicting peak ground acceleration based on microtremor measurements and soil boring test data. Therefore, the records obtained from these measuring stations and boring test sites can provide useful information for developing the model by using neural network approach. In the next section, the processing of measured and test data is discussed, and then the results obtained from the developed model are presented.

3. Ambient Vibration Measurement and Soil Boring Test Data

Ambient vibration with very low amplitude (about 10−6 m) and acceleration (0.8–2.5 gal or cm2/s), which cannot in general be felt by humans [21, 22], is frequently found on the ground surface of the earth. This vibration, however, can be recorded by using a precise measuring instrument developed recently. In the present study, the results of ambient vibrations are measured and calculated by a set of ultrasensitive seismic accelerometer-Model 731A-made by Wilcoxon with other monitoring instrument and computer software [2325]. The original microtremor information is stored in the frequency and time domains, but only the data set in time domain is considered for analysis as the characteristics of ground motion is the primary concern in this study.

From microtremor data collected at the four measuring stations, it can be found that the accelerations for measuring station MS1 are basically in the range 0.43 gal to 51.71 gal in both the east-west (EW) and north-south (NS) directions, and the results in vertical (V) direction are all smaller than 0.58 gal. For measuring station MS2, the results are 0.33 gal–24.1 gal, 0.27 gal–47.49 gal, and 0.78 gal–80.19 gal, in EW, NS, and V directions, respectively. For measuring station MS3, the results are 0.33 gal–5.03 gal, 0.27 gal–4.29 gal, and 0.31 gal–6.13 gal, in EW, NS, and V directions, respectively. For measuring station MS4, the results are 0.58 gal–45.85 gal, 0.34 gal–15.46 gal and 1.05 gal–49.52 gal, in EW, NS and V directions, respectively.

The above numerical results show that the accelerations in the vertical direction are all relatively higher than those of the other directions, except at the measuring station MS1. As ambient vibrations may result from moving vehicles and construction work, some of the measured data may exhibit much higher values particularly between 8:00 AM and 5.00 PM, and these can affect the true microtremor response. Therefore, in this study, only ambient vibration data collected between 8:00 PM and 7:00 AM are taken for analysis so as to eliminate outside environmental factors as far as possible and increase the accuracy of natural microtremor response measurements.

For soil test data in the boreholes, the samplings are by auger boring method for soils above ground water level and by the method of wash boring for soils below ground water level. The laboratory tests conducted were for general physical, triaxial compression, shear strength, unconfined compressive strength, consolidation, compaction, California bearing ratio, resilient modulus, and groundwater quality. The in-situ tests were for lateral load and permeability. All of these test results can provide soil characteristics in detail for each of the drilling sites.

For a typical soil exploration and testing report in the research area, it can be seen that there are many items such as soil depth, soil profile, USCS (unified soil classification system) classification, standard penetration value, grain size analysis, water content, specific gravity, density, void ratio, liquid limit, plasticity index, and safety factor against liquefaction calculated from shear strength parameters [26]. Note that some of these soil test items may have their own physical meaning and also have a relationship with each other. Previous studies have found that the three important parameters relevant to the problem of actual earthquake response are STP-N, D50, and FL [2729]. Hence, these three parameters, with ambient vibration surveys and distance between measuring station and test layer within bore hole (DS), are considered for developing a PGA prediction model by the neural network approach. Typical values for STP-N, D50, FL, DS, and PGA in the four directions are shown in Table 1.

4. Neural Network Model and Analysis of Prediction Results

In the field of computational intelligence, neural network approach is widely applied in various engineering applications as it has some attractive features such as easiness to implement, strong pattern recognition capability, and good prediction performance [3032]. Basics of neural network modeling such as selecting a suitable architecture, learning algorithms, preprocessing of data, training, and testing of models have been comprehensively covered in many publications [3335]. Thus, further discussions of this method and the use of associated software tool are not included here, but only some of the key points for using this computational technique that are relevant to the present research problem are addressed below.

To develop a neural network model, it is essential to determine the number of neurons in the input layer, the hidden layer, and the output layer. In this study, five soil boring test data in the neighborhood of each microtremor measuring station are used in developing the model. As mentioned in the previous section, the soil input parameters included are STP-N, D50, and FL. Therefore, a total of 50 data sets are available for the five bore holes, as each bore hole has a 20 m depth, and the soil profile is divided into 10 layers. If the distance parameter DS defined previously is also included, then there are four neurons in the input layer. The PGA in each of the different directions (EW, NS, V, and H) obtained from microtremor measurements is used as the target, resulting in only one neuron in the output layer. The number of neurons in the hidden layer needs to be selected to provide a relatively better performing neural network model.

It is better to examine the input soil data sets in advance to find a suitable neural network prediction model. Initially, we consider a neural network model without the distance parameter as shown in Figure 3(a), and divide the normalized data sets into three groups, where 70% is for training, 20% for verification, and 10% for testing. These three calculation stages are performed in Matlab toolbox with the “train,” “adapt,” and “simulate” functions [36, 37]. The computational experiments, with the use of correlation coefficient ( ) as evaluation index, showed that the training result can achieve high values (from 0.656 to 0.900) with random data selection in the network calculation, but it has a poor performance in the verification cases (from 0.004 to 0.235) and the testing cases (from 0.001 to 0.361), as seen in Table 2. The poor performance of this model is due to random data selection with no rational basis for the association of the PGA values that are collected over time with the soil properties that are defined spatially within the bore holes. Thus, a rational basis for associating the data values for the input and output variables is required, and available domain knowledge is used as the basis for this association.

Because the STP-N value refers to soil hardness, it can play an important role in influencing the degree of liquefaction during an earthquake. Hence, by taking the STP-N as the primary factor and arranging its data set to increase from small to large values, with corresponding adjustments to the other input parameters, the target PGA data set is then arranged from large to small values. Again, without considering the distance parameter in the input layer, the performances of the neural network models with different number of neurons in the hidden layer were considered, and the model with four neurons in the hidden layer has a relatively better performance than the other models for all three calculation stages.

Note that only 50 data sets were used for developing neural network model in this study, so it is not suitable to choose too many layers or neurons in the hidden layer as it may cause ineffective learning during the training stage. The neural network model with four neurons in the hidden layer is found to be more reliable and will also be used with the distance parameter in the input layer.

Table 3 shows the performance of the preferred neural network model I4H4O1 (Figure 3(b) for the three microtremor measuring stations. It can be seen that the average values are quite high and up to 0.998 at the training stage. The average values range from 0.767 to 0.872 at the verification stage and from 0.714 to 0.812 at the testing stage, all exhibiting reasonably high coefficient of correlation between measurement and estimation. That is, the developed neural network model has a sufficient level of prediction capability and can be used for further investigation.

The above results from the neural network models demonstrate that microtremor measurements may have a relationship to the soil profile. It is crucial to check the capability and apply the developed model for predicting PGA at an unmeasured site. To perform this task, it can be assumed that microtremor measuring station MS4 is an unknown site, then the known PGA values at this station can be used for verifying the ability of neural network model. To estimate PGA in station MS4 from the three known stations MS1, MS2, and MS3, the straightforward method is by distributing the results of these three known stations with weighting factors based on the distances between stations and denoted here by “Model 1.”

Alternatively, a better way to estimate PGA at an unmeasured site is by taking a new set of soil data from five drilling holes nearby, and insert the data set solely to the neural network model developed for each of known measuring stations. Then by summing the results with weighting factors in accordance with the distances between the unknown site to the three known stations, the final estimation is obtained for the unknown site, and this method is denoted here as “Model 2”.

The comparison of prediction results for the two models and microtremor measurements in the different directions are shown in Figure 4. It can be seen that the neural network estimations are not too different for both models, but the results of “Model 2” seem to be slightly closer to the actual measurements. Because the total recorded ambient vibration surveys for MS4 are 70 data sets, the chosen 50 data sets cover a wider range of time interval compared to the other three stations. Therefore, some of the measurements exhibit higher values (PGA > 0.005 g), particularly for the east-west and vertical directions, and these may cause an error in the prediction. In general, the performances of both neural network methods have reasonable accuracy and are acceptable for the problem considered. This comparison of results provides confidence for using this method for prediction of PGA in an unmeasured but important site.

Installation of ambient vibration survey instrument is usually at a place where the density of high-tech buildings is not too high. Thus, this method is suitable for an unmeasured site (UMS) shown in Figure 2, which is closer to several important industrial buildings and is possibly sensitive to natural microtremor. The present approach of developing neural network model is thus very useful for predicting PGA at this site with the use of new soil data sets from nearby five bore holes. Figure 5 shows the prediction result for the unmeasured site. It can be seen that both models exhibit similar predictions, but the curve obtained from “Model 2” is not as smooth as for “Model 1.” The reason may be that the soil bore holes used in one of the known measuring stations, MS3, are too far away from the unmeasured site. Actually, this local instability problem is also found in the previous comparison shown in Figure 4, but it is believed that the present neural network model should still be sufficiently reliable. The sequential selection of the data sets for developing the neural network models results in the networks predicting the results for the test data sets by extrapolation. This accounts for the difference between the model predictions and the measurements on the right hand side of the curves in Figure 4. The predictions could be further improved by a random selection of data sets, for training, verification, and testing, after the reordering of the data sets as described previously.

If the measured results in the closest station MS2 are taken as reference, then the prediction result shows that PGA in vertical direction is larger than the other two directions, which is consistent with the measured results, and this is also true for predicted PGA in the horizontal direction. Since this science park is mostly alluvial soil, the prediction result is thus reasonable as there is no significant change of soil conditions in the unmeasured site. Overall, the neural network “Model 2” seems more preferable as the distance parameter and new soil data set are used in the calculation process of this model and, hence, may represent a more true response of the investigating site. In addition, it may be concluded that the model can be applied to predict PGA in any site of interest around the Kaohsiung science park.

5. Conclusion

Without a doubt, the seismic-related problems are very important research topic in the field of disaster prevention technology. This study presented a novel way of using neural network approach to develop a model for learning a relationship for linking two different types of parameters, that is, the ambient vibration measurement and the on-site soil boring test data. In addition, a weight-based neural network model is also developed for predicting peak ground acceleration at an unmeasured site, and is extendable for predicting natural microtremor on any other site of the science park investigated, as long as the soil conditions are suitably distributed.

Due to the limitation of soil boring test data for each microtremor measuring station, the present study picked up only fifty sets of peak ground acceleration from each of microtremor measuring station, to match with soil profiles from five drilling holes nearby for developing neural network model. More data sets might be required to develop a more accurate model for performing the prediction. Nonetheless, the results obtained in this study provide an insight into the seismic-related characteristics in the research area. One of the significant aspects of the present research is that even though the data collected at MS1, MS2, and MS3 are for different time periods, the three neural network models developed for these sites can be combined to predict the PGA distribution at a fourth site (MS4) for a different time period.

Further, it should be mentioned that the predicted results do prove the reliability of the developed model, but the choice of microtremor measurements with commonly found frequency or peak ground acceleration should be further investigated. This may occur repeatedly and hence can affect the accuracy of precision instrument and damage products such as semiconductors or biosensors during the manufacturing process within a high-tech company. Nevertheless, the method used in this study did provide a new way to treat this type of nonlinear problem, and may be applicable in other areas of interest around the world.

Acknowledgments

The support of National Science Council under Project NSC 101-2221-E-020-030 is greatly appreciated. The authors also wish to thank Dr. J. H. Lai of KaoYuan University for providing valuable inputs including ambient vibration surveys and soil boring test data from CECI engineering consultants. Dr. D. Gunaratnam is retired recently from Faculty of Architecture, Design and Planning, University of Sydney, Australia.