Measurement of Shear Strain in Map Liquefaction Area for Earthquake Mitigation in Bengkulu City

Bengkulu city is well-known as a high-risk area for an earthquake. One of the effects of the earthquake is liquefaction, the phenomenon in which the strength and stiffness of the soil are reduced by shock. The researchers measured the quantity of liquefaction by knowing the value of Vs30 values that can be obtained from the United States Geological Survey (USGS). In this study, however, we measured liquefaction based on the shear strain value (γ). Here, the shear strain values could be determined by measuring the seismic vulnerability index and the peak ground acceleration. The vulnerability index was obtained by measuring the resonance of the frequency and the amplification factor using a short period seismometer with 3 components (2 horizontal components and 1 vertical component), whereas the peak ground acceleration was obtained from the earthquake data that happened in Bengkulu City in 2000. The result showed that Vs30 and shear strain negatively correlated.


Introduction
Bengkulu Province often experiences natural disasters of an earthquake where causes liquefaction, a phenomenon whereby soil loses strength and behaves like a liquid due to cyclic force during an earthquake, so the pore water pressure exceeds the effective soil stress.As the water pressure increases due to the earthquake, the soil effective stress (σ') decreases.Liquefaction usually occurs in less dense and water-saturated soil (1), where the soil layer strength loss due to shocks caused by an earthquake or other vibration (2), so that if a shock occurs, soil layer behaves like a liquid.Liquefaction usually happens in water-saturated soil whose pores are full of water.When it gets vibration, the water gives pressure to the soil particles results in the decreasing soil density.Before the earthquake, the water pressure is low, but after that, the pressure increases so it can move the soil particles easily.Bengkulu City lies in the latitude of 3°45'-3°59' S and longitude of 102°14'-102°22' E, within the subduction area of Indo-Australia and Eurasia tectonic plates.Consequently, Bengkulu City is prone earthquake and liquefaction.
An earthquake of 7.9 on the Richter scale struck Bengkulu City in 2000 showed existing liquefaction in many locations.This liquefaction is related to the measurement of peak ground acceleration and shear wave velocity (3).In general, the condition of liquefaction represented by the peak ground acceleration should reach more than 300 gal (4).Our survey showed that the peak ground acceleration in all locations in Bengkulu City was higher than 300 gal.There is a high probability of liquefaction occurrence in Bengkulu if an earthquake occurs (5).One of liquefaction occurred during an earthquake in Bengkulu in 2000, where the impact is shown in Figure 1.
Geologically, almost all areas in Bengkulu have soft material such as sand, gravel, mud, and clay.Soft soil is not easily traveled by shear wave (Vs).In other words, soft soil has low Vs (6).This soft soil has high vulnerability index, meaning that this soil is highly prone to liquefaction.This statement is supported by the study of Huang and Tseng (2002) showing that in Yuan Lin, the most damaged area which experienced liquefaction was located in the zone of high seismic vulnerability index based on microtremor (7).The measurement of liquefaction area based on soft soil condition which has low shear wave velocity will be correlated with the area having a shear strain in the same location.The correlation has not been done, but it needs to be proven.The objective of this study was to measure the potential threat of liquefaction to areas in Bengkulu City using microtremor data.Data were acquired using short period seismometer with three components, namely horizontal component of north-south, the horizontal component of east-west and a vertical component.The equipment for seismometer was sampling frequency of 100 Hz, digitizer, solar cell panel, GPS, UPS, and laptop for data acquisition.The software DATAPRO was used for microtremor data acquisition, GEOPSY for HVSR analyses and ArcGIS for mapping.Acquisition of primary data was made by direct microtremor measurement in each zone.The microtremor measurement was done in 100 locations, and each location was measured for 30 minutes with a sampling frequency of 100 Hz.The technique for microtremor survey followed the standard of SESAME European Research Project 2004 as shown in Table 1.The seismometer for measurement is shown in Figure 2, while the example of microtremor data display from survey results is shown in Figure 3. -Determine the level gain (magnification) as high as possible without saturation of signal.
In-situ Coupling soil-censor -Adjust the sensor directly on the ground.
-Avoid placing seismograph censor on the soft soil surface (mud, shrub) or wet soil.
Artificial Coupling soil-censor -Avoid using soft plate material such as rubber or sponge.
-In steep slope where it is difficult to get a good level, place the sensor under the sand pile or in a container filled with sand.The presence of buildings/trees -Avoid measurement close to buildings, tall buildings, and tall trees if the wind speed is above ± 5 m/second, because it will affect the results of HVSR analyses as indicated by the presence of low frequency in the curve.-Avoid measurement in a parking place, water pipes, and culvert.
-Weather condition -Wind: Protect the sensor if the wind speed is above ± 5 m/second).
-Rain: Avoid measurement under heavy rain; light rain will not have much effect.
-Temperature: check the sensor temperature and follow the manufacturer's instruction.The analysis of tremor data resulted in seismic vulnerability index, which then was multiplied by the peak ground acceleration to get the shear strain value ().So far, the measurement of liquefaction vulnerability is usually done using an indicator of soil hardness represented by the shear wave velocity at a depth of 30 meters (Vs30) below the soil surface.To show that shear strain can be used as measurement indicator for the potential threat of liquefaction to an area, the shear strain indicator must be compared with Vs30.The ratio between shear strain and Vs30 was shown in a graph, so the measurement of the potential threat of liquefaction could be directly measured from the graph.

Research Method
Generally, liquefaction may occur due to certain factors, one of which is the increase in water pressure as a result of earthquake shocks, so that the effective stress (σ') decreases.This condition can be stated as follows: where: The soil effective stress is calculated using this equation: To measure the potential threat of liquefaction of the soil, two important variables must be calculated, namely (1) cyclic stress ratio (CSR), cyclic stress causing liquefaction and (2) cyclic resistance ratio (CRR), soil ability to resist liquefaction.The value of CSR is calculated using this equation:  9) as follows: 1.5 1.5 0.5 1.5 2 1.0 0.4113 0.04052 0.001753 1.0 0.4177 0.05729 0.006205 0.00121 The value of cyclic resistance ratio (CRR) is calculated using field data, namely standard penetration test (SPT), cone penetration test (CPT), and shear wave speed (Vs).Based on the values of CSR and CRR, the Factor of Safety of FS is calculated using this equation: Liquefaction happens if FS ≤ 1 and vice versa, it does not happen if FS > 1.In this study, whether or not liquefaction occurs was determined based on the shear strain (), which was calculated from the microtremor measurement in each study site.The microtremor data were analyzed with Geopsy using HVSR method.An empirical approach was used to calculate shear strain (10), and the result is formulated in the equation (7).
  In Equation (7)  is a shear strain, f 0 resonance frequency, A amplification factor and a max peak ground acceleration which was calculated using Kanai attenuation in (11), formulated in equation ( 8) (11).
In equation, T g is soil vibration period, M earthquake magnitude moment, and dan R is the distance from hypocenter to the station.Shear strain triggers deformation of soil surface layer.Liquefaction deformation occurs if the value of shear strain is higher than 1×10 -3 .Therefore, the measurement of liquefaction potential in this study was conducted using shear strain measurement (12).

Data Acquisition
Micrometer data acquisition was conducted in all zones in Bengkulu City using short period seismometer with three components, namely horizontal component of north-south, the This study also needed the earthquake data: the magnitude, the depth and the epicenter coordinates of the 2000 earthquake.The magnitude of that earthquake was 7.7 on Richter Scale, located at a latitude of-4.72 S and longitude of 102.09E. The earthquake data were needed to measure the peak ground acceleration using equation (8).In equation ( 8) the value of M=7.7, the value of T was calculated using this formula, T=1/f, and the value of R was determined using Pythagorean theorem as shown in Figure 5.  9) was used to calculate the distance between hypocenter and the station.
where R i is the distance between hypocenter to the i th station; h is the depth of hypocenter; (x,y) the coordinate of the epicenter, and (x i , y i ) the coordinates of the i th station.
After the values of R in each station were known, the value of a max in each station could be calculated.After the values of a max , A and f 0 were known, the value of  could be calculated.Then, the value of  in each station was integrated with the value of Vs30 in the same coordinates.The correlation between  and Vs30 was a mathematical correlation in the form of a function.If the mathematical correlation showed a certain function has a relatively high coefficient of regression, then the potential of liquefaction of an area can be determined using the value of  in that area.The flow of data from the beginning to the end is shown in Figure 6.

Results and Analysis
The potential of liquefaction in Bengkulu City determined from the distribution of shear strain  as shown in Table 2.The distribution of shear strain values indicating the potential for liquefaction is shown in Figure 7.The colors show the potential for liquefaction of an area, ranging from brown (very potential), red (potential), yellow (not potential or safe) and green (safe).The brown area in Figure 7 has shear strain value > 1×10 -2 , so it has very high potential for liquefaction.The red area has the potential for liquefaction because its shear strain value > 1×10 -3 as stated by ( 12) that liquefaction occurs if 1×10 -3 . Referring to Ishihara statement, the potential for liquefaction is distributed evenly in all areas of Bengkulu City, especially in the coastal area, and a certain area has very high potential for liquefaction as shown in Figure 5.The percentage of the area has the potential for liquefaction in Bengkulu City is very large as shown in Figure 8.Based on the potential for liquefaction, the areas in Bengkulu City can be categorized as very potential (1.83%), the potential (49.54%), safe (46.79%), and very safe (1.83%).The areas potentially and highly potentially to experience liquefaction in Bengkulu City are relatively large in term of land availability for settlement, offices, public facilities, and infrastructure.According to data published in "Bengkulu Dalam Angka" (13), the area for settlement gets larger over time, reaching 16% of the total city area in 2016.Most of the settlement area is in the coastal area.This is a bad situation because this area has high population density, many office buildings, and infrastructures.Liquefaction in this area will result in a large loss because most large physical structures will sink into the ground.
Liquefaction in this area will cause not only material loss but also social problems such as the presence homeless (refugee) and diseases.Geologically, the areas in Bengkulu City vulnerable to liquefaction have certain lithology, namely alluvial containing sand, clay, and mud.The study of (2) showed that the geology of Bengkulu City is dominated by Alluvium Terrace's rock which is soft and vulnerable to earthquake shocks.The study of (2) supports the theory that an area having high shear strain is vulnerable to liquefaction if an earthquake occurs (2).The result of Sugiyanto's study of the local geology of Bengkulu is shown in Figure 9.The results of this study are in accordance with the statement of ( 7) that the area in Yuan Lin which experienced liquefaction and suffered heavy damage was located in high shear strain zone based on micrometer data.The study of ( 14) which states that high shear strain value is located in the coastal area also supports our study (14).Almost all areas in Bengkulu have high shear strain value and are located in Alluvium Terraces.Nakamura et al. (10) conducted microseismic measurement in 400 locations in the area which suffered heavy damage during Kobe earthquake in 1965 (15).The measurement results showed a correlation between seismic vulnerability index and percentage of damage.The areas which suffered damage had a high seismic vulnerability, and vice versa.One of the damage is liquefaction.An area having high seismic vulnerability index has a high shear strain (16) automatically.
The areas which have high shear strain is located in the soft rock which has low shear wave velocity.Noutash et al. (17) state that measurement of liquefaction using wave shear velocity as an indicator may be conducted although it is not better than measurement using SPT (Standard Penetration Test).However, measurements of liquefaction potential using wave shear velocity as an indicator have been conducted by many researchers such as (3), ( 18) and (19).Theoretically, wave shear velocity is negatively correlated with shear strain.If the wave shear velocity is low, the shear strain is high.This correlation is shown in Figure 10.Where Y is a shear wave velocity and X is a shear strain.Having coefficient of regression r 2 =0.708, this function is reliable with strong confidence.The measurement of liquefaction potential of an area is important for land use planning of a city.Based on liquefaction potential map, the city administration is ready to make mitigation if an earthquake occurs.Since Bengkulu City lies in an earthquake-prone area, using liquefaction map as part land use planning shows that the city's policy has a high concern for the community's interest.

Conclusion
The measurement of liquefaction potential in an area can be done by measuring the shear strain in that area.This measurement may not be as good as that using the CSR and CRR methods, but the measurement of shear strain can be easily conducted because we just need to measure soil vibration and peak ground acceleration.The measurement of shear strain can produce not only the map of liquefaction potential, but also detection of the landslide, rockfall, beach erosion, and sediment thickness.The equipment to measure shear strain is portable and easy to operate.If it is correlated with the liquefaction using shear wave velocity data, the measurement of shear strain can be done directly using formula (10).

Figure 1 .
Figure 1.The building sink due to liquefaction

TELKOMNIKA
Vol. 16, No. 4, August 2018: 1597-1606 1600 ρ=bulk density of soil (kg/m 3 ) ρ'=bulk density of water (kg/m 3 ) h'=depth soil water surface (m) Measurement of Shear Strain in Map Liquefaction Area for... (Muhammad Farid) 1601 horizontal component of east-west and a vertical component.The data distribution in each zone is shown in Figure 4.

Figure 10 . 1605 Figure 10
Figure 10.The correlation between shear strain and shear wave velocity

Table 2 .
The Result of Micrometer Data Analysis Into Shear Strain